GROMACS使用教程.doc

GROMACS教程一Gromacs基本模拟流程 31下载pdb文件 32用pdb2gmx处理pdb文件 33建立盒子 35设置能量最小化 46用grompp程序进行文件处理 67使用genion和tpr文件添加离子 68 69在后台运行能量最小化(在命令后加&) 7二设置位置限制性动力学模拟 7三设置非限制性动力学模拟 91如何重启一个计算 112如何延长一个计算 113如何设置并行计算 11五模拟结果分析 121如何将特定帧的轨迹保存成*.pdb文件 122用ngmx观察轨迹文件(也能够用VMD观察轨迹文件) 123比较常用的分析工具 14g_covar计算斜方差 16g_energy能量数据作图,如压力、体积、密度等 16g_gyrate测量回旋半径 17g_rms与g_rmsdist计算结构的RMSD值 17g_rmsf计算原子位置的根均方波动(rmsf) 18do_dssp计算模型的二级结构 20g_hbond计算模拟过程中分子间的氢键的数目、距离或角度 21g_saltbr分析模拟中残基间的盐桥 21GROMACS是一个使用经典分子动力学理论研究蛋白质动力学的高端的高效的工具。GROMACS是遵守GNU许可的免费软件,能够从以下站点下载:,并且能够在linux和Windows上使用。在本教程中,将研究一个从漏斗形蜘蛛的毒液中分离的***。我们将使用显性溶剂动力学的方法来进行研究。首先比较真空中和溶解的模型。我们将把***肽溶在水盒子里,紧接着用牛顿运动定律加以平衡。我们还将比较偿离子在显性溶剂动力学中的影响。更全面的用法指导请参考官网的GROMACS用户手册注意:在本教程中,将要生成的gromacs(*.gro)结构文件,能够用VMD(下)查看。一Gromacs基本模拟流程1下载pdb文件()2用pdb2gmx处理pdb文件pdb2gmx–ignh–ffG43a1–f–o–p–waterspcepdb2gmx此命令将pdb文件转换成gromacs文件并产生拓扑文件。-ignh因为本pdb文件是由NMR产生的,含有氢原子,因此用-ignh选项忽略文件中的氢原子。-ff指定力场(G43a1是Gromos96力场,一个通用原子力场)。-f读入pdb文件,-o指定一个新产生的pdb文件(也能够是其它多种类型文件)的文件名。-p指定新产生的拓扑文件名。拓扑文件包含了所有力场参数(基于一开始选择的力场),因此非常重要。-water来指定水模型研究表明SPC/E水模型在水盒子模拟中表现最好。用SPC/E水模型研究长程静电相互作用较好。#注:对于下面将要用到的任何命令,都能够使用“-h”查看该命令的使用方法,比如,对于命令pdb2gmx能够使用:pdb2gmx–h3建立盒子editconf-btcubic–f–o–d用上面的命令建立了一个简单的立方体盒子.-d决定了盒子的尺寸,即盒子边缘距离分子边缘(9Å)。理论上在绝大多数系统中,-d。注:editconf也能够用来进行gromacs文件(*.gro)和pdb文件(*.pdb)的相互转化。例如:editconf–f–o转换为现在就能够用产生的文件进行真空模拟了。真空模拟就是先能量最小化,然后进行动态模拟。4在盒子中放入溶剂genbox–cp–cs–o–pgenbox命令在editconf产生的盒子基础上生成水盒子。上面的命令行指定了SPC水盒子。genbox命令能够在给定尺寸的盒子中加入正确数目的水分子。5设置能量最小化:Gromacs用*.mdp文件指定所有计算的参数。它用最速下降法消除原子位置碰撞。编辑文件,将nsteps变成400。如果最小化不能收敛,就用nsteps=500再做一次。(最小化在400步内一般是能收敛的,但不同的平台可能结果会不一样。)要重做的话,必须重新运行grompp(注意:预处理器的位置在你的机器上可能不同,用which命令来定位,即whichcpp):title–标题随便取(最长64个字,简单点好)cpp–指定预处理器的位置define–传递给预处理器的一些定义。–DFLEXIBLE告诉grompp将spc水模型而非刚性SPC包含进拓扑文件,以便用最陡下降法进一步最小化能量。constraints–设置模型约束integrator–steep,告诉gompp使用最速下降法进行能量最小化。cg则代表使用共轭梯度法。dt–能量最小化用不用。只在动力学模拟中用(如md)。nsteps–在能量最小化中,指定最大运行步数。nstlist–更新邻居列表的频率。nstlist=10表示每10步更新一次。rlist–短程邻居列表的阈值。coulombtype–告诉gromacs如何计算静电。PME为particlemeshewald法(参见Gromacs用户手册)rcou