Gromac命令锦集.docVIP

目录 TOC \o 1-3 \h \z \u HYPERLINK \l _Toc2703673931、make_ndx PAGEREF _Toc270367393 \h 2 HYPERLINK \l _Toc2703673942、g_traj PAGEREF _Toc270367394 \h 3 HYPERLINK \l _Toc2703673953、g_energy另外一种用法 PAGEREF _Toc270367395 \h 3 HYPERLINK \l _Toc2703673964、pdb2gmx PAGEREF _Toc270367396 \h 3 HYPERLINK \l _Toc2703673975、genion PAGEREF _Toc270367397 \h 5 HYPERLINK \l _Toc2703673986、editconf PAGEREF _Toc270367398 \h 6 HYPERLINK \l _Toc2703673997、g_rama PAGEREF _Toc270367399 \h 6 HYPERLINK \l _Toc2703674008、g_cluster PAGEREF _Toc270367400 \h 7 HYPERLINK \l _Toc2703674019、genbox PAGEREF _Toc270367401 \h 7 1、make_ndx Gromacs的索引文件，即index文件，由make_ndx程序生成，文件后缀为.ndx。 索引文件是gromacs最重要的概念文件，使用它可以在模拟过程中为所欲为。举一个简单的例子，比如想详细了解HIV整合酶切割DNA的反应机理，使用量子力学模拟反应位点的反应过程，而分子其他部位使用一般分子动力学模拟。于是我们就面临一个对模拟系统进行分割定义的问题，在gromacs中，就要用到索引文件。 基本的思路是这样的，在索引文件中，定义一个独立的组，这个组包括反应位点处所有原子。在模拟的.mdp文件中，对这个组定义量子力学模拟，事情就是这么简单。对蛋白进行量子力学模拟时，一般使用洋葱模型。所谓洋葱模型，就是对反应位点使用量子机制，在反应位点一定的半径内，使用半量子力学机制，然后分子部分使用分子机制。那么索引文件就定义一个使用量子力学的组，把需要引进量子机制的原子都放到这个组中；再定义一个半量子机制的组，同时放进需要半量子力学机制模拟的原子，再在.mdp文件中独自定义即可。 再举一个例子，比如说在进行SMD（Steered MolecularDynamics，这个我一直没有想到或者找到切恰的中文翻译方法，或许可以叫做牵引分子动力学？？别扭！！）中，要对蛋白莫一个原子或者残基作用力，那么可以建立一个索引文件，在该文件中定义一个组，把要施力的残基或者原子放到该组中。然后在.ppa文件中使用该组就行了。 如果我还没有说明白，那么看看gromacs的参考文件吧。如果还是不明白，可以来找我，我免费培训。^_^ 索引文件使用make_ndx命令产生，make_ndx -h可以看到全部的参数。运行make_ndx后，可以使用 r 命令选择残基， a 命令选择原子， name 命令多组进行改名。可以使用 | 表示或运算， 表示与运算。下面是几个简单的例子： 1.选择56号残基 r 56 2.选择3至45号残基 r 3-45 3.选择3至15，23至67号残基 r 3-15 | r 23-67 4.选择3至15号残基的主干链原子 r 3-15 4 #在索引文件中，4号组为默认的主干链。 组合是灵活的，使用的时候好好发挥聪明的大脑啦。 个人觉得gromacs把索引文件概念做得非常好，并独立成一类文件是一个不小的创举啊。在这个概念很值钱的年代，基本上使gromacs多一个大大的卖点。 ? 2、g_traj 几乎gromacs的所有分析数据都可以输出为xmgrace的数据文件，g_traj可以产生gromacs轨迹的各个组的坐标，速度，受力和边界等等。使用” -com “参数可以求出轨迹中各个组的质心的坐标，速度，受力等；使用” -mol “则可以求取系统中各个分子的信息；” -ot “则可以求出系统中各个组的温度。还有几个其他参数，比如” -cv “可以求平均速度，” -cf “可以求平均受力等。其他参数同其他命令无异，参加说明文件。 3、g_energy另外一种用法 Gromacs的各个工具都很有个性，如果互相结合，可以做很多事情。 g_energy求系统轨迹各个能量的，一般跑完MD之后，使用g_energy处理ener.edr只能得到系统的各个能量项。但是如果想求系统中两个不同部分在模拟过程中的相互作用能量，那就要使用一些小窍