| 查看: 3374 | 回复: 8 | |||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||
御剑江湖荣誉版主 (著名写手)
天池冶海
|
[交流]
【转帖】GROMACS分析程序讨论，在gromacs中添加Amber力场已有7人参与
|
||||||
|
gromacs的分析程序很多，不知到各位平时都用到哪些呢？现在遇到的困惑是学会了做MD，但不会深入分析结果。我先说说我自己的,比较常规，抛砖引玉： g_energy：能量提取和分析 g_dist：研究两个基团之间的距离波动 g_rms：研究某些基团的位置偏离 g_rmsf：计算温度因子，查看哪些残基波动较强 g_rmsf可取平均结构，用-ox参数输出 g_gyrate：分析蛋白质的回旋半径变化 蛋白质的回旋半径反应了蛋白质分子的体积和形状。同一体系的回旋半径越大，说明体系发生了膨胀。 g_gyrate –f md.trr –s md.tpr –o gyrate.xvg g_sas：分析溶剂的可及表面积（SASA） 它是描述蛋白质疏水性的重要手段，氨基酸残基的疏水性是影响蛋白质折叠的重要物理作用。 g_sas -f md.trr -s md.tpr -o area.xvg -or resarea.xvg -oa atomarea.xvg do_dssp：分析蛋白质的二级结构变化 group务必选择protein，否则…… do_dssp –s md.tpr –f md.xtc –o ss.xpm g_dist和g_bond：提取两个原子间距离随时间变化 提取两个原子距离随时间变化，用g_dist和g_bond都可以。 g_dist的index写法是，在里面加入两个group，内容分别是这两个原子序号，运行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后会提示选择group，分别选择那两个group就行了，结果输出在dist.xvg，第二列是距离，后面三列是x/y/z的距离。 如果用g_bond，则index里新加入一个group，把两个原子都写进那一个group下。用的时候选这个group就行了，结果和g_dist的第二列的内容一样。 gromacs的分析程序很多，不知到各位平时都用到哪些呢？现在遇到的困惑是学会了做MD，但不会深入分析结果。我先说说我自己的,比较常规，抛砖引玉： g_energy：能量提取和分析g_dist：研究两个基团之间的距离波动g_rms：研究某些基团的位置偏离g_rmsf：计算温度因子，查看哪些残基波动较强 g_rmsf可取平均结构，用-ox参数输出g_gyrate：分析蛋白质的回旋半径变化蛋白质的回旋半径反应了蛋白质分子的体积和形状。同一体系的回旋半径越大，说明体系发生了膨胀。 g_gyrate –f md.trr –s md.tpr –o gyrate.xvg g_sas：分析溶剂的可及表面积（SASA）它是描述蛋白质疏水性的重要手段，氨基酸残基的疏水性是影响蛋白质折叠的重要物理作用。 g_sas -f md.trr -s md.tpr -o area.xvg -or resarea.xvg -oa atomarea.xvgdo_dssp：分析蛋白质的二级结构变化 group务必选择protein，否则…… do_dssp –s md.tpr –f md.xtc –o ss.xpmg_dist和g_bond：提取两个原子间距离随时间变化 提取两个原子距离随时间变化，用g_dist和g_bond都可以。g_dist的index写法是，在里面加入两个group，内容分别是这两个原子序号，运行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后会提示选择group，分别选择那两个group就行了，结果输出在dist.xvg，第二列是距离，后面三列是x/y/z的距离。如果用g_bond，则index里新加入一个group，把两个原子都写进那一个group下。用的时候选这个group就行了，结果和g_dist的第二列的内容一样。 gromacs的分析程序很多，不知到各位平时都用到哪些呢？现在遇到的困惑是学会了做MD，但不会深入分析结果。我先说说我自己的,比较常规，抛砖引玉：g_energy：能量提取和分析g_dist：研究两个基团之间的距离波动g_rms：研究某些基团的位置偏离g_rmsf：计算温度因子，查看哪些残基波动较强g_rmsf可取平均结构，用-ox参数输出g_gyrate：分析蛋白质的回旋半径变化蛋白质的回旋半径反应了蛋白质分子的体积和形状。同一体系的回旋半径越大，说明体系发生了膨胀。g_gyrate –f md.trr –s md.tpr –o gyrate.xvgg_sas：分析溶剂的可及表面积（SASA）它是描述蛋白质疏水性的重要手段，氨基酸残基的疏水性是影响蛋白质折叠的重要物理作用。g_sas -f md.trr -s md.tpr -o area.xvg -or resarea.xvg -oa atomarea.xvgdo_dssp：分析蛋白质的二级结构变化group务必选择protein，否则……do_dssp –s md.tpr –f md.xtc –o ss.xpmg_dist和g_bond：提取两个原子间距离随时间变化提取两个原子距离随时间变化，用g_dist和g_bond都可以。g_dist的index写法是，在里面加入两个group，内容分别是这两个原子序号，运行比如g_dist -f a-sol-md.xtc -s a-sol-md.tpr -n index.ndx，之后会提示选择group，分别选择那两个group就行了，结果输出在dist.xvg，第二列是距离，后面三列是x/y/z的距离。如果用g_bond，则index里新加入一个group，把两个原子都写进那一个group下。用的时候选这个group就行了，结果和g_dist的第二列的内容一样。 在gromacs中添加Amber力场 （GROMACS4.5自带Amber力场） Install FF: 所有信息来自http://folding.stanford.edu/ffamber/ 1.下载力场 ： http://folding.stanford.edu/ffamber/ffamber_v3.3.1.tar.gz (490 kB) http://folding.stanford.edu/ffamber/ffamber_v3.3.1-doc.tar.gz (6.54 MB)带文献例子 2.解压缩 3.把aminoacides.dat、vdwradii.dat拷到top文件夹下覆盖，实际上一股脑考过去就完了。 4.子文件夹下的东西都拷到外面，也就是top里面。 5.修改FF.dat,修改第一行的数字，加上新添的力场数。 6.为每个新的力场增加新行，如： ffamber94 AMBER94 Cornell protein/nucleic forcefield ffamber99 AMBER99 Wang protein/nucleic acid forcefield ffamber99p AMBER99p protein/nucleic forcefield ffamber03 AMBER03 Duan protein/nucleic forcefield [ Last edited by 御剑江湖 on 2011-4-9 at 09:48 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 精华 | GROMACS | gromacs | 力场 |
» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
求Gromacs中Glycam力场数据包。。非常感谢
已经有5人回复- Gromacs中加入离子
已经有5人回复
- 【讨论】gromacs电荷不为整数
已经有6人回复
- 【讨论】gromacs 并行 mpirun
已经有4人回复
- 【求助】如何使用gromacs中的constraint
已经有8人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【求助】GROMACS分析二面角随时间的变化
已经有4人回复
- 【原创】在gromacs的gro中插入小分子
已经有12人回复
- 【求助】Gromacs和NAMD比Amber快很多吗？如何选择？
已经有7人回复
7楼2012-02-16 10:32:45
changle
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1800.6
- 红花: 1
- 帖子: 95
- 在线: 97.5小时
- 虫号: 1203702
- 注册: 2011-02-15
- 性别: GG
- 专业: 生物大分子结构与功能
2楼2011-04-08 17:21:39
奋斗1s
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1258.1
- 散金: 109
- 红花: 1
- 帖子: 415
- 在线: 106.7小时
- 虫号: 980571
- 注册: 2010-03-24
- 专业: 有机合成
3楼2011-05-11 12:07:11
changle
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1800.6
- 红花: 1
- 帖子: 95
- 在线: 97.5小时
- 虫号: 1203702
- 注册: 2011-02-15
- 性别: GG
- 专业: 生物大分子结构与功能
4楼2011-05-12 09:45:10