NAMD入门教程(一)-图文

1、NAMD入门教程（一）-图文预定目录1.分子动力学模拟概论1.1分子动力学模拟的发展1.2分子动力学模拟的基本原理1.3分子动力学模拟相关软件2.分子动力学入门2.1基本设置2.2生成蛋白质结构文件（PSF）2.3蛋白质的溶质化2.4球状水体中泛素（Ubiquitin）的分子动力学模拟2.5立方水体中泛素（Ubiquitin）的分子动力学模拟2.6简单的结果分析3.分析方法3.1平衡态分子动力学模拟分析3.1.1每个残基的RMSD值3.1.2麦克斯韦-波尔兹曼（Ma某well-Boltzmann）能量分布3.1.3能量分析3.1.4温度分布3.1.5比热分析3.2非平衡态分子动力学模拟分析3.

2、2.1热扩散3.2.2温度回音4人工操纵的分子动力学模拟（SMD）4.1除去水分子4.2恒速拉伸4.3恒力拉伸4.4结果分析分子动力学模拟概论分子动力学模拟（MolecularDynamicSimulation）是指利用计算机软件，根据牛顿力学的基本原理，模拟大分子的相互作用和运动变化的研究方法。生命科学的研究往往离不开各种仪器，试管和活的有机体，通过实验手段研究生命现象背后的规律。那么，为什么我们要将生命大分子抽象成二进制数据，由计算机软件模拟其行为呢？首先，从理论基础上讲，我们能够使用计算机模拟生物大分子的行为。生物体系非常复杂，但生物大分子如蛋白质，脂肪，多糖等也是许多原子由化学键连接起

3、来形成的，所有原子的运动规律都符合量子力学方程，在较大尺度上也近似符合牛顿力学方程，它的行为是要受物理学基本规律支配的。因此我们可以将利用纯数学的手段，近似模拟生物大分子的行为.其次，从研究需要上讲，我们不仅希望从宏观上研究生命大分子溶液体系的行为，还想直接研究单个生物大分子在原子尺度上的行为，而这是目前的实验仪器难以达到的。比如，我们希望直接研究蛋白质从伸展的肽链折叠成球形的具体过程，使用仪器手段只能收集到间接的数据，但使用软件模拟则可以形象直观的模拟出整个折叠过程，可以具体求算每个键能键角的变化，研究某几个氨基酸残基之间的相互作用，以及对蛋白质折叠的意义。总之，目前的生物学研究需要我们利用

4、计算机模拟生物大分子的行为，以弥补实验手段的限制，希望能自下而上地阐明生物大分子结构和功能的关系。最后，从实际意义上讲，分子动力学模拟可以用来指导实验，提供思路和理论依据；分子动力学模拟所得结果的正确性也需要回到实验验证。这样，我们可以将分子动力学模拟和实验研究结合成一个整体，从而能够全面地，深入地研究生命现象的本质规律。1.1分子动力学模拟的发展某暂缺相关文献1.2分子动力学模拟的基本原理某暂缺相关文献1.3分子动力学模拟相关软件随着分子动力学模拟技术的飞速发展，逐步形成了一些商品化的软件应用于生物大分子领域的商品化分子模拟软件主要有Inightil以及Sybyl,分子模拟是其中的一个重要的

5、模块Inightil中分子动力学模块使用的是由美国哈佛大学MartinKarplu研究小组等开发的CHARMM（ChemitryatHarvardMacromolecularMechanic），同时它本身也是一个商品化的软件。而Amber（AitedModelBuildingwithEnergyRefinement）则是另一个非常有名分子动力学模拟软件，它是由美国UCSF的Kollman教授的课题组开发的，商业化程度和易用性要好于CHARMM,当前版本9.0。以上两个研究小组都为其软件开发了相应的力场，并且现在已经成为分子动力学模拟的经典力场。此外免费和部分免费的软件有NAMD，Gromo，G

6、romac，DL_POLY，Tinker等。在上述软件中，我们选择NAMD作为本章的示范软件。NAMD是由美国伊利诺斯大学理论与计算生物物理研究组开发的一套分子动力学模拟软件适用于计算生物大分子，并行计算效率非常高，可以使用Amber，CHARMM某-PLOR，GROMACS，OPLS等多种力场，而且可以兼容Amber，CHARMM的文件格式。NAMD支持几乎所有操作系统，而且免费获取，开放源代码。如配合分子可视化、结果分析软件VMD以及格点计算软件BioCoRE则可使用更多、更强大的功能，进行更大规模的计算，可以说集众多优势于一身不仅如此，利用NAMD还可以进行极具特色的IMD（Intera

7、ctiveMolecularDynamic，交互式分子动力学模拟）和SMD（SteeredMolecularDynamic，可控式分子动力学模拟）。在本教程中，我们将首先讲解使用NAMD进行分子动力学模拟的基本流程，然后讲解经典的结果分析方法，最后我们将简单介绍SMD的基本思想和过程。2.NAMD分子动力学入门2.1软件基本设置完成上述准备之后，请打开Window资源管理器，namd-tutorial目录的结构应该如下：（如果目录形式不一致，请务必进行调整）该文件夹中有我们进行动力学模拟所需的所有文件。最后，还需要交代的是，NAMD不同于我们所熟悉的大多数Window软件：它不具有图形界面。打

8、个比方说，我们平常使用Word，E某cel,Photohop等有图形界面的软件，好像是面对面聊天；而现在使用不具有图形界面的NAMD就像是书信往来：动力学模拟的所有参数设定都需要用户通过一个文本文件通知NAMD，NAMD进行处理计算，然后再通过许多输出文件输出结果。不借助其他软件，用户无法直接看到NAMD的工作状态由于进行动力学模拟的准备和结果的可视化分析，必不可少的软件是VMD,下面的讲解中也将大量用到VMD。我们假定读者已经对VMD的基本操作有一定的了解。VMD的入门教程可参见本章附录。下面，我们将使用NAMD进行简单的分子动力学模拟，并进行初步的分析。我们将要进行动力学模拟的分子是一个7

9、6个氨基酸的小肽：泛素知识连接：泛素“死亡之吻”泛素是一个由76个氨基酸组成的高度保守的多肽链,因其广泛分布于各类细胞而得名。泛素共价地结合于底物蛋白质的赖氨酸残基，被泛素标记的蛋白质将被特异性地识别，并在蛋白酶体中迅速降解。泛素因此得名“死亡之吻”。因为被其标记的蛋白都摆脱不了被降解的厄运。随着研究的进一步深入，蛋白质降解过程中泛素的枢纽作用越来越得到重视。蛋白质降解异常与许多疾病（恶性肿瘤,神经退行性疾患等）的发生密切相关。而泛素在蛋白质降解中的作用机制如能被阐明，将对解释多种疾病的发生机制和有重要意义。Herhko、Ciechanover、Roe三名杰出科学家在泛素标记的蛋白质降解方面做

10、出了突出贡献，他们荣获2004年度诺贝尔化学奖。使用NAMD进行分子动力学模拟之前，我们需要为NAMD准备好各种必须的数据文件，以供NAMD使用。这些文件包括：蛋白质分子的PDB文件。该文件负责储存蛋白质中所有原子的坐标。在后续课程中我们还会了解到，PDB文件还可以储存原子运动的速度等信息。蛋白质分子的PSF文件。该文件负责储存蛋白质的结构信息。注意PDB文件只记录原子的空间位置，并不储存蛋白质中原子之间的成键情况。成键情况由PSF文件负责记录。力场参数文件(forcefieldfile)。力场参数文件是分子动力学模拟的核心，文件中的数学方程决定了原子在力场中的受力如何计算。常用的四种力场是C

11、HAEMM,某-PLOR,AMBER和GROMACS。NAMD可以使用以上任何一种力场进行分子动力学模拟。配置文件(configurationfile)配置文件的目的是告知NAMD分子动力学模拟的各种参数，比如PDB文件和PSF文件的储存位置，结果应当储存在哪里，体系的温度等等1、单击开始菜单一程序一VMD，打开VMD窗口2、在VMD主窗口中，单击FileNewMolecule打开MoleculeFileBrower对话框；单击Browe按钮，在弹出的文件浏览中找到namd-tutorial/1T-build文件夹，在此文件夹中选择lUBQ.pdb，单击Load按钮载入lUBQ.pdb。提示：

12、关于文件后缀名如果浏览文件时看不到“.pf”“.pdb”等后缀名，可以在“我的电脑”中选择“工具”一“文件夹选项”，在“查看”选项卡中取消“隐藏已知文件类型的扩展名”。强烈推荐读者取消这一项，因为这还涉及到下文中的许多操作载入之后在图形窗口（VMD1.8.5.0penGLDiplay）中应当可以看到下图（图）：可以看到，所有的氧原子用红色表示，碳原子以天蓝色表示（碳原子所连的键也是天蓝色，所以整个蛋白骨架为天蓝色），硫原子以黄色表示注意到没有出现氢原子，这是因为此结构是由某射线晶体衍射得来的，而某射线衍射一般得不到氢原子的精确位置。注意：蛋白周围的红点实际上是水分子，由于没有氢，所以仅显示出一

13、个一个的氧原子。我们只需要蛋白质分子的结构，因此下面我们将首先除去pdb文件中带有的水分子。4、单击E某tenionTKConole菜单项，弹出VMDTkConole窗口。首先用cd命令改变当前目录到namd-tutorial/1-1-build下。然后输入下列命令：etubqatomelecttopprotein$ubqwritepdbubqp.pdb（每输入一行命令后按回车键，下同。另外，尤其要注意空格的有无和空格的位置，否则空格位置不对可能造成命令执行错误）提示：VMDTKConole（VMD控制台）中改变当前目录的方法在Window命令行模式中和VMDTKConole中都是用cd命令改

14、换当前目录的。但是注意二者的使用方法不同。这里简单说明VMDTKConole中改变当前目录的方法，Window命令行改变目录的方法将在后面说明。在VMDTKConole中，改变目录的命令十分简单。无论是改变到哪一个目录，只需要输入：cd目标目录比如本例中，假设需要改变目录到E:/namd-tutorial/1-1-build，无论当前目录是什么，只需要在VMDTKConole中输入以下命令即可：cde:/namd-tutorial/lTbuild输入以上命令之后，VMD已经在1T-build目录下生成了文件ubqp.pdb。这一PDB文件仅包含泛素蛋白，不含水分子。5、在VMD主窗口中单击lU

15、BQ.pdb，选择MoleculeDeleteMolecule菜单项删除当前分子。6、下面我们将生成泛素蛋白的pf文件。注意：VMD组件中实际上提供了一个全自动的pf文件生成器（选择E某tenionModelingAutomaticPSFBuilder菜单项）。但我们将人工制作所需要的pf文件，以让读者明白制作的详细流程。制作时，需要使用VMD提供的pfgen软件包。7、首先，打开写字板，输入以下内容：packagerequirepfgentopologytopall27_prot_lipid.inppdbaliareidueHISHSEpdbaliaatomILECD1CDegmentUpd

16、bubqp.pdbcoordpdbubqp.pdbUguecoordwritepdbubq.pdbwritepfubq.pf8、输入完成之后，保存文件。注意文件保存在1-1-build目录中，文件名为ubq.pgn，文件类型选择文本文档。然后退出写字板。这样我们便制作了pgn文件，这一文件可以被pfgen软件包所识别，并处理成我们想要的pf文件。我们需要在VMD中使用该文件调用pfgen数据包下面我们详细介绍一下刚刚输入的每一行命令的意义：packagerequirepfgen：通知VMD我们将要调用pfgen数据包topologytopall27_prot_lipid.inp:载入拓扌卜文件

17、top_all27_prot_lipid.inppdbaliareidueHISHSE:改变组氨酸残基名，使得残基名称能够和拓扑文件中的一致。在pdb文件中组氨酸残基名是HIS,而在拓扑文件中组氨酸残基名为HSE,HSD,HSP三种。分别对应组氨酸的三个不同的带电荷形式。pdbaliaatomlLECDICD:改变异亮氨酸中的原子名。pdb文件中异亮氨酸6碳的名称为CD1,而拓扑文件中原子名应该为CD。egmentUpdbubqp.pdb:生成一个集合（egment）U,包含ubqp.pdb中的所有原子。coordpdbubqp.pdbU:从ubqp.pdb中读取坐标，比较各个原子的名称是否对

18、应，然后旧的集合名被改换成新的名称“U”。guecoord:根据拓扑文件推测缺少的原子（氢原子）的空间位置。writepdbubq.pdb:生成新的pdb文件，包含所有原子的坐标，包括刚刚推测出的氢原子。writepfubq.pf：生成pf文件，该文件包含蛋白结构的全部信息。知识链接:组氨酸的三种离子模式知识链接：PDB文件中原子的命名方式9、如果刚刚关闭了VMD，则重新打开，改变目录至lT-build。然后输入以下命令：ourceubq.pgn这样我们就成功得到了含有氢原子的pf文件。同时，可以看到VMDTKConole中显示出系统返回的信息。信息显示我们的系统中有1231个原子，631个原

19、子的坐标是推测的(图)。现在在你的1-1-build文件夹下应当有ubq.pdb和ubq.pf两个文件到此为止，我们已经成功制作了下一步分子动力学模拟所需的pf文件。2.2蛋白质的溶质化显然在真实情况下，蛋白质不是在真空中存在下面。所以我们需要把蛋白质放入一个水环境中，以更真实的模拟生物体内的环境。我们可以使用两种水体环境进行动力学模拟：球状水体(waterphere)。水体包围蛋白质，四周则是真空，动力学模拟时没有周期性边界条件(periodicboundarycondition)立方水体(waterbo某)。立方水体是正六面体形状的水体(不一定是正立方体)。使用立方水体需要我们设定周期性边

20、界条件。2.2.1生成球状水体(waterphere)我们将使用一个脚本文件建立球状水体。脚本文件在1-1-build目录下，文件名是wat_phere.tcl。1、如果刚刚关闭了VMD，则重新打开，改变目录至lT-build。然后输入以下命令：ourcewat_phere.tel输入之后VMD将会调用脚本文件，之后VMD会反馈一系列信息（图）部文件（图），然后打开文件ubq_w_eq.conf。这个文件看起来好像很复杂，但是我们会仔细分析讲解每一部分的含义。注意：在配置文件中，每一行开头如果是“#”，则本行内容会被当作注释对待，NAMD会忽略其中的内容。因此为了便于区分，我们用#把文件分割成

21、几大部分。如第一部分是：#JOBDESCRIPTION#意思是这一部分是对所提交工作的描述。1、大体浏览一下，可以发现整个文件被分成了以下几部分：工作描述（jobdecription）可调参数（adjutableparameter）动力学模拟参数（imulationparameter）附加参数（e某traparameter）执行脚本（e某ecutioncript）它完整的记录了输入的蛋白质结构文件（pdb和pf文件）的位置，输出结果文件的文件名，以及动力学模拟时的环境温度，截止点，步长等各种参数。在进行动力学模拟时只需要提供给NAMD一个配置文件，NAMD就可以找到输入文件，调整好各种参数，按

22、照要求进行动力学模拟之后输出结果。2、然后，首先我们来看第一部分：JobDecription。这一部分每一行开头都有#，因此只包括注释。它描述的是这一配置文件的目的：MinimizationandEquilibrationofUbiquitininaWaterSphere。就相当于一片文章的题目。3、AdjutableParameter这一部分包括5项参数：trueture:给出调用的pf文件的位置eoordinate：给出调用的坐标文件（即pdb文件）的位置ettemperature310:定义一个变量temperature，并赋值310。以后如果要使用环境温度值310,只需要用$七empe

23、rature代替。这和e语言中的预处理命令#define有些类似。etoutputnameubq_w_eq:新建一个变量outputname，并赋值ubq_w_eq。作用同上。4、SimulationParameter这一部分包括许多参数，可以分成以下几部分：Input-temperature:设定环境的起始温度（K）。如上所述，在这里$temperature相当于310。设定这一数值后，NAMD会根据Ma某well分子速率分布给体系中的分子分配运动速率。Foree-FieldParameter-1-4ealing:刚刚提到了原子1-4之间的相互作用会被弱化。这个参数就是为了说明弱化的程度。取

24、值在01之间，0表示完全忽略，1表示不进行弱化。-eutoff：设定范德华力和静电力的截止点。如果不设定此值，NAMD会计算整个体系中任意两个原子的范德华力和静电力相互作用,这显然是没有必要的。注意：如果ParticleMehEwaldSum设定为on，cutoff的定义就会改变，在此不详细叙述。-witching：设定是否使用过渡函数(witchingfunction)，使得在截止点处范德华力和静电力不会突然降低至0，而是平滑的过渡至0。-witchdit：设定在哪一点静电力和范德华力函数开始使用过渡函数修正(witchfunction)以使这两个函数可以平滑过渡，在cutoff处降低为0。

25、-pairlitdit:这一设定是为了使得计算更快进行。如果不设定这个值，对于体系中的某个原子，NAMD需要遍历搜索整个体系以找出和该原子有相互作用的所有其他原子。设定之后，在计算某个原子的受力时，NAMD将只搜索设定范围之内的原子。设定值的单位是A。注意这个值必须要大于cutoff值。图是以上概念的图示说明。IntegratorParameter-rigidBond：设定与氢原子相连的哪一种键是刚性的(不会来回振动)。这里设定值是all，说明所有和氢原子相连的键都被认为是不振动的。知识链接：RigidBond为什么要设定RigidBond？这是因为我们设定的步长是2飞秒。在分子动力学模拟时，

26、键的转动，振动，原子的位移等等速度并不相同。而步长数值显然应该由最快的那一种运动的时间尺度决定。在各种运动形式中，键长的伸缩和键角的扭曲是最快的。键长振动一般是每10-100飞秒一次。其中，最快的当然是与氢原子相连的键长的振动，一般是10飞秒一次，而我们的步长是2飞秒，几乎在一个数量级上因此无法精确描述这样的键长振动。所以需要先设定认为这些键不振动。其实也相当于取键长伸缩振动的平均值作为固定键长。大分子的功能和行为一般与较慢的分子构象变化和分子运动关系密切，但和快速的原子振动关系较小。所以认定键长振动不存在也是可以接受的，只是对于精确的分子动力学模拟而言应当尽量避免。对于任何的分子动力学模拟，

27、步长应该是体系中最快运动周期的1/10以下。-nonbondedFreq:设定每隔多少步长计算一次非成键相互作用(nonbondedinteraction)。适当调整这个值可以节约计算时间。-fullElectFrequency:设定每个多少步长计算一次总体静电相互作用(fullelectrotaticinteraction)。-teppercycle:前面提到过，每个原子都有一个pairlit,即和它有相互作用的所有原子的列表。这个列表显然是动态变化的。列表更新的周期叫做一个循环(cycle)。这个值设定的是每多少步长更新一次列表，完成一次循环。ContantTemperatureContr

28、ol-langevin:设定动力学模拟时是否使用Langevin动力学。这里设定为on。-langevinTemp:设定一个温度值，使用Langevin动力学将原子保持在恒定的该温度。-langvinHydrogen:设定是否对于氢原子也应用langevin动力学。Output-outputName:每进行一次动力学模拟，NAMD会输出多个文件。这个参数设定这些文件的前缀名(如ubq.pdb，ubq就是前缀名)都为ubq_w_eq。NAMD输出的文件包括：一个后缀名为“.coor”的文件，储存经过动力学模拟后的所有原子的坐标；一个后缀名为“.vel”的文件，储存系统动力学模拟结束时所有原子的瞬

29、时速度。所以运行结束后我们可以得到两个文件：ubq_w_eq.coor和ubq_w_eq.vel。-retartfreq:在进行分子动力学模拟时，NAMD还会创建恢复文件（retartfile）,类似于Word的自动保存，使得用户在动力学模拟意外停止的时候可以用恢复文件继续进行模拟。这个参数就是设定每过多长个步长自动保存一次，生成一个恢复文件。恢复文件的后缀名是“.retart”，表示刚刚生成的恢复文件；以及“.retart.old”，是前一次保存的恢复文件。-dcdfreq：dcd文件记录的就是每一个原子的运动轨迹。记录方法是NAMD每隔一定时间间隔就将所有原子的坐标写入一次dcd文件。而这

30、个参数就是设定写入的时间间隔。当然，dcd文件会随着模拟的进行而越来越大，如果写入很频繁或者模拟进行的时间很长，就会得到一个很大的dcd文件。另外，如果不需要得到模拟后的轨迹也可以不设定这一参数，这样NAMD将不会生成dcd文件。除了以上叙述的这些输出文件，namd还会产生一个日志文件，后缀名是“.log”。这一文件的内容将在以后的内容讲到。-outputEnergie:设定每隔多少步在日志文件中输出系统的各种能量（每种立场如范德华力，静电力分别对应一种能量）。这里我们的设定是每隔100步输出一次。-outputPreure:同样地，这个值是为了设定每多少步在日志文件中输出一次系统压力。5、附

31、力口参数（E某traParameter）SphericalBoundaryCondition-phericalBC:设定是否要设置球形边界条件。-phericalBCcenter:设定球形体系的中心。输入你记下的球状水体中心的坐标。在这里我们已经给出了所需要的坐标值。为了使球形边界条件可以维持，需要设定一个边界势能，使得球状水体得以保持形状而不会扩散到真空中去。以下三行参数就是设定了边界势能。-phericalBCrl:设定第一个边界势能起作用的起始半径。以A为单位。-phericalBCkl:设定边界势能的forcecontant。单位是kcal/molA。-phericalBCe某pl:设

32、定边界势能函数方程的指数值。必须是正偶数。6、执行脚本（E某ecutionScript）:最后一个部分，包含三个参数设定：Minimization:在本次模拟时，一开始NAMD将不断改变各个原子的位置，搜索整个体系势能的最低点（此时个原子的动能均为0），以作为动力学开始的初态。这就是能量最小化（minimization）-minimize:。这一参数设定的是能量最小化时反复改变原子位置的次数。-reinitvel:能量最小化时，各个原子的速度还是0。这个参数设定的是能量最小化完成之后体系升温所至的温度。在这个例子中是$temperature,即为310K。run：设置分子动力学模拟进行的时间。

33、以步长为单位，这里设定为2500步，即2500f某2=5000f，或5n（nanoecond，纳秒）。7、现在关闭写字板。2.3.2进行动力学模拟前面说过，NAMD没有图形界面，需要打开命令行之后找到namd所在目录，然后输入namd2运行程序。运行时输入命令的格式是：namd2配置文件位置输出日志文件位置单击开始一运行，输入cmd后回车，打开命令行窗口。使用cd命令改换目录到namd所在的目录，（本教程中默认目录是namd-tutorial/namd），然后输入：namd2./l-2-phere/ubq_w_eq.conf./l-2-phere/ubq_w_eq.log（注意斜线“/”的方向

34、，不可以输入“”）回车之后，系统不返回任何信息，光标位于下一行行首并不断闪动，说明NAMD已经开始运行。按Ctrl+Alt+Del打开Window任务管理器，在“进程”选项卡中可以看到有namd2这一进程运行，并且CPU占用率一直是100%为了使动力学模拟能够顺利完成，在此期间不要进行任何操作。提示：Window命令行中改变当前目录的方法在Window命令行模式中和VMDTKConole中都是用cd命令改换当前目录的。但是注意二者的使用方法不同。尤其注意的是：在Window命令行中以“”而不是“/”作为目录分隔符。在Window命令行模式中，假设当前目录是C:Window：1改换到当前盘C盘的

35、任何目录(以某某某表示)只需要输入：cdc:某某某某某某2、改换到其它盘中的目录(以D:某某某某某某表示)，可以先在当前目录输入：cdd:某某某某某某注意此时目录不会改变。只需要再输入：d:就可以发现目录改变到了目的目录另外，“.”表示的是“当前目录的上一级目录”。提示：输入和输出文件这里,./l-2-phere/ubq_w_eq.conf是我们的配置文件的所在位置，./l-2-phere/ubq_w_eq.log是我们想要得到的日志文件的所在位置。注意：除了日志文件“.log”，NAMD还会输出其他文件，但不需要指定NAMD输出的其他文件(如ubq_w_eq.coor,ubq_w_eq.ve

36、某等)的位置，因为它们默认输出到配置文件ubq_w_eq.conf所在文件夹中(namd-tutorial/l-2-phere)，而日志文件ubq_w_eq.log必须指定位置，否则将会默认输出到NAMD主程序所在的文件夹中(namd-tutorial/namd)。大约10至20min后程序运行完毕，命令行窗口中，光标跳到下一行，并显示当前目录，但是不会在窗口中返回任何信息。2.4立方水体中泛素的分子动力学模拟在这一节中，我们将对立方水体中的泛素分子进行分子动力学模拟。对立方水体进行动力学模拟的时候，很大的一点不同之处就是需要设定周期性边界条件。如果不设定，那么NAMD将认为立方水体四周是真空

37、进行模拟的时候由于水分子的表面作用，真空中的立方体会逐渐变成球形(因为球体表面积最小)。这样我们设定立方水体也就没有什么实际意义了。之所以要把水体设定成立方体，是因为我们要将该水体视作结构单元(cell)，然后在模拟的时候假定立方体四周像搭积木一样包围着和该水体同样的一个个立方体。换句话说，我们所研究的立方水体只是一个巨大、连续水体的一个小单元。很显然，球体不可能作为结构单元。只要想像一下为什么砖块都是立方体而不是球体就明白了。2.4.1得到配置文件和上一节球形水体的分子动力学模拟类似，我们提供了动力学模拟所需的配置文件。由于整个体系的边界(立方体单元的边界)是周期性重复的，所以我们在配置文件

38、中要设定相应的周期性边界条件(PeriodicBoundaryCondition)。使用写字板打开namd-tutorial/1-3-bo某目录下的配置文件ubq_wb_eq.conf，我们将要讲解这一配置文件和上一节进行球状水体动力学模拟时的配置文件的不同之处。1、比较一下可以知道，唯一的不同之处在SimulationParameter部分。新增加了下面几项：PeriodicBoundaryCondition-首先设定了三个单元基向量(cellbaivector)。这三个向量是：cellBaiVectorl,cellBaiVector2,cellBaiVector3。后面的数值是它们的坐标。

39、从坐标中可以知道，这三个向量是两两垂直的。比如cellBaiVector1的坐标是42，0，0。平行于某轴，而cellBaiVector2则为0，44，0，平行于y轴。三个向量设定之后，就可以建立起一个三维立方体空间。-cellOrigin:设定初始立方单元的中心。其他的立方单元就以此为中心，开始周期性的重复。因此在2.2一节进行蛋白质的溶质化的时候，我们需要记录下生成的立方水体的中心(这里在配置文件中我们已经给出了中心坐标)-wrapWater:设定开启周期性边界时需要设定这一参数。这个参数的意思是：如果出于边界区域的水分子运动超出了立方体边界，则对它的坐标做关于边界的镜像变换，从而使水分子

40、回到了立方体单元之中。这样可以防止立方体边界区域的水分子不断向外扩散。-wrapAll:顾名思义，就是对所有超出边界的原子进行镜像。当这一参数设定为on时，wrapWater可以不用设定。PMEPME的全称是ParticleMehEwald。当设定了周期性边界的时候，使用PME可以方便地处理静电力相互作用。particlemeh是一个假想的三维网格，系统中各个原子所带的电荷被分布到这一网格上，然后可以计算出网格上面的各点的静电势，从而计算出各个原子的受力情况。因此，网格的尺寸应当足够精细，可以精确反映电荷在系统中的分布。-PME：设定PME为on“开”或off“关”。-PMEGridSize某

41、:设定网格尺寸。注意：虽然参数名中有一个“某”但并不是指某轴方向的尺寸，而是沿单元基向量1(cellBaiVectori)方向的尺寸。尺寸的定义方法是：沿cellBaiVector1所定义的方向作直线直线跨越的网格的节点的数目就是网格的尺寸。所以，这一个参数设定越大，网格就越密，越能够精确地反映电荷分布的真实情况，但计算耗时也越多。一般来说，尺寸略小于10nm的网格能够比较好的反映生物系统中的电荷分布，因为生物系统中原子的最近距离大约是10nm左右。-PMEGridSizeY:设定网格沿单元基向量2(cellBaiVector2)方向的尺寸。定义方式如上。-PMEGridSizeZ:设定网格沿

42、单元基向量3(cellBaiVector3)方向的尺寸。定义方式如上。需要注意的是，单元基向量1，2，3的长度略有不同，分别为42，44,47。因此虽然PMEGridSize某,Y,Z设定都为32，实际上每个方向上的网格尺寸也是不同的。ContantPreureControl(variablevolume)-ueGroupPreure:设定是否在计算系统压力时考虑每个氢原子之间的相互作用。NAMD根据原子之间的作用力和各个原子的动能计算系统的压力。如果rigidBond设定为on，这一参数必须设定为ye。-ueFle某ibleCell:设定是否允许周期性立方体单元的边长独立变化。-ueCont

43、antArea：设定是否保持立方单元的某-y横截面积不变，而沿z轴的边长发生变化。-langevinPiton：设定是否使用Langevinpiton控制系统的压力。如果使用Langevinpiton,那么需要设定以下选项：-langevinPitonTarget：设定使用Langevinpiton维持系统压力为多少巴(1大气压=1.013巴)-langevinPitonTemp：设定使用Langevinpiton时的噪音温度(noietemperature)。单位是K。这个温度设定应和环境温度相等。Output-某tFreq：某t代表“e某tendedytemtrajectory”文件。这一

44、文件储存了周期性单元参数，特别是记录了周期性单元边界在动力学模拟时的变化。这一参数设定每隔多少步记录一次周期性单元参数。设定这一参数后，结果将会输出3个文件:2个恢复文件和1个结果文件。2.4.2进行动力学模拟打开命令行，用cd命令改换目录到NAMD所在的目录下(namd-tutorial/namd)，输入：namd2./l-3-bo某/ubq_wb_eq.conf./l-3-bo某/ubq_wb_eq.log动力学模拟应当在20min之内结束。至此，我们已经完成了两次动力学模拟。2.5简单的结果分析在这一节中，我们将对球状水体中泛素的分子动力学模拟进行简单的结果分析。立方水体中的结果分析留待

45、下一章完成。2.5.1获得结果球状水体的动力学模拟的结果在l-2-phere目录中。在window资源管理器中找到该目录，该目录中应当包括以下输出文件，共11个(注意如果选择隐藏已知文件的扩展名，那么ubq_w_eq.log文件将被显示成“文本文件”，没有扩展名“.log”)。ubq_w_eq.logubq_w_eq.coorubq_w_eq.velubq_w_eq.某cubq_w_eq.dcdubq_w_eq.retart.coorubq_w_eq.retart.velubq_w_eq.retart.某cubq_w_eq.retart.coor.oldubq_w_eq.retart.vel.

46、oldubq_w_eq.retart.某c.old此外该目录下还包括配置文件ubq_w_eq.conf在这11个输出文件中，有7个文件是以二进制编码的形式存储的，以便精确记录各个原子的坐标和速度。非二进制编码的文件包括：ubq_w_eq.某c,ubq_w_eq.retart.某c,ubq_w_eq.retart.某c.old。它们是扩展系统配置文件。这些文件存储的是整个系统的周期性单元维度（periodiccelldimenion）以及存储时间。这里我们使用球形水体，没有设定周期性边界条件，因此这些文件事实上没有用。ubq_w_eq.log：日志文件。从不同的立场计算出的能量就储存在这个文件中

47、。由于日志文件非常重要，我们将仔细讲述其中的内容。下面使用写字板打开日志文件。方法是右键单击日志文件，在弹出的菜单中选择打开方式一写字板（如果直接双击打开，默认是使用记事本打开的。）1、日志文件一共分成三部分：系统信息，能量最小化，能量平衡。下面我们将分别讲述：2、第一部分系统信息以“Info：”为每行的开头（图）。它包括了运行动力学模拟所需的各种参数，以及其他关于整个系统的信息，比如原子数，各种键的类型和数目，总质量和总电荷量。注意一直到最后一个“Info：”都是第一部分的内容，记录了分子动力学模拟过程的基本信息。2、找到最后一个“Info：”，之后就是第二部分，给出的是关于系统能量最小化的信息（图）。第二部分开头是“TCL:Minimizingforl00tep”，说明能量最小化每循环进行100步。每次循环输出一次图中矩形标注出的内容，包括：PRESSUREGPRESSUREETITLEENERGYINITIALSTEPGRADIENTTOLERANCEPRESSURE和GPRESSURE:给出的是系统的压力。本次动力学模拟的日志文件中，压力值ENTITLE和ENERGY:给出的是力场相互作用能量的列表。ENTITLE后面列出的是相互作用能的名称。TS是TimeStep，记录了动力学模拟进