Gromacs的genion命令
在给蛋白质添加水环境之后，一般要在水环境中添加金属离子，使模拟系统更加接近真实系统。如果系统中蛋白质本身已经带了静电量，那么就更要给系统加几个带相反电量的金属离子，使系统处于电中性。
gromacs中添加金属离子的命令是genion，使用" genion -h "可以得到其使用的参数，其中有几个比较常用：
---------------------------
" -s: "指定系统tpr文件。
" -p: "指定系统拓扑文件，在往系统中添加金属离子时，genion会往拓扑文件最后的分子类型中写入添加的离子数，并修改拓扑文件中系统原子数。
" -o: "指定输出文件，genion的输出是pdb文件或者gro等结构文件。也就是说你产生这个文件之后，还要再用这个文件产生tpr文件。
" -np/-nn: "带正/负电金属粒子的数目。这个数目有一点讲究，一般需要看个人的应用。假如想要得到" 0.1mol/L " 的离子浓度到底要加多少，可以自己算一下（很简单，方法很多，比如看课本）。也可以直接使用" -conc " 参数直接指定离子浓度，在使用" -conc "参数时，建议使用" -neutral "参数配合，即使系统的最后处于电中性。嗯，gromacs开发组想得不要太周到哦（南京话）。
" -pn/-nn "指定正负金属离子的名字，比如" NA+ "或者" CL- "。可以看看gromacs安装途径" share/gromacs/top/ " 下面你用的力场文件中离子到底用什么名字，也可以使用新的离子，但是要在力场中定义，或者把新离子的itp文件使用" include "添加到系统拓扑文件中。
" -random "随机位置添加离子。这个比较有说法，如果不用该参数，那么离子就会添加在势能最低处，即靠近蛋白质相反电量的部位。有的人说这样可以节省时间，但是这样一跑MD，离子就会抱跟着蛋白，死死不放，不是很好。
" -seed "有随机，就是随机种子，如果发现使用" -random "添加离子自由，离子离蛋白太近（比如说小于0.1nm），那么可以指定新的seed。
---------------------------
嗯，给一个例子：
------
" genion -s topol.tpr -o system_ion.pdb -p system.top -np 100 -pname Na -nn 100 -nname Cl -random "
------
说明一下：加了100个Na和100个Cl，随机加，输出文system_ion.pdb文件。
再使用grompp生成一个新的tpr，就可以计算了。
搜索更多相关主题的帖子:
Gromacs genion 命令
