| 查看: 3698 | 回复: 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 | 力场 |
» 猜你喜欢
288求调剂
已经有6人回复
-
复试调剂
已经有12人回复
-
计算机11408，286分求调剂
已经有3人回复
-
材料专硕322分
已经有8人回复
-
一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业083200
已经有5人回复
-
081200-11408-276学硕求调剂
已经有5人回复
-
材料调剂
已经有10人回复
-
313求调剂
已经有7人回复
-
304求调剂（085602，过四级，一志愿985）
已经有14人回复
-
材料工程302分求调剂
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 求Gromacs中Glycam力场数据包。。非常感谢
已经有5人回复
- Gromacs中加入离子
已经有5人回复
- 【讨论】gromacs电荷不为整数
已经有6人回复
- 【讨论】gromacs 并行 mpirun
已经有4人回复
- 【求助】如何使用gromacs中的constraint
已经有8人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【求助】GROMACS分析二面角随时间的变化
已经有4人回复
- 【原创】在gromacs的gro中插入小分子
已经有12人回复
- 【求助】Gromacs和NAMD比Amber快很多吗？如何选择？
已经有7人回复
changle
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1800.6
- 红花: 1
- 帖子: 95
- 在线: 97.5小时
- 虫号: 1203702
- 注册: 2011-02-15
- 性别: GG
- 专业: 生物大分子结构与功能
2楼2011-04-08 17:21:39
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
5楼2011-12-21 09:39:10
6楼2011-12-21 11:23:33
7楼2012-02-16 10:32:45
fanc232
新虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: 1238.4
- 散金: 206
- 红花: 3
- 帖子: 253
- 在线: 108.7小时
- 虫号: 1046781
- 注册: 2010-06-24
- 性别: MM
- 专业: 细胞生物学研究中的新方法
8楼2015-03-19 16:23:32
白玉浴血
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2513.6
- 散金: 102
- 红花: 8
- 帖子: 838
- 在线: 243.2小时
- 虫号: 2684358
- 注册: 2013-09-27
- 性别: MM
- 专业: 高分子组装与超分子结构
9楼2015-06-03 17:43:52
