《【优质】三种常用分子模拟软件介绍.docx》由会员分享,可在线阅读,更多相关《【优质】三种常用分子模拟软件介绍.docx(17页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、【优质】三种常用分子模拟软件介绍三种常用分子模拟软件介绍一、NAMDNAMDNAnoscaleMolecularDynamics是用于在大规模并行计算机上快速模拟大分子体系的并行分子动力学代码。NAMD用经历力场,如Amber,CHARMM和Dreiding,通过数值求解运动方程计算原子轨迹。1.软件所能模拟的体系的尺度,如微观,介观或跨尺度等微观。是诸多md软件中并行处理最好的,能够支持几千个cpu运算。在单机上速度也很快。模拟体系常为为10,000-1,000,000个原子。2.软件所属的类型,如MD,DPD,DFT,MC,量化,或穿插等全原子md,有文献上也用它做过cgmd。3.软件能研
2、究的相关领域,使用者的背景最好是?使用的力场有charmm,x-plor,amber等,合适模拟蛋白质,核酸,细胞膜等体系。可以进行团簇和CNT系统的模拟软件原理经典,操作简单。但需要对体系的性质足够了解。4.软件中主要涉及的理论方法范畴经典的md,以及用多种方法计算自由能和SMD模拟。数据分析时候一般很少涉及复杂的热力学和统计热力学的原理,但知道一些最好。5.软件主要包含的处理工具namd是计算部分,本身不能建模和数据分析unix的哲学kiss。但vmd同namd系出同门,已同namd实现无逢链接。vmd的tcl脚本一定要搞懂,别的就不多介绍了。26.与此软件密切相关的软件vmd,及其他数据
3、统计分析软件excel,OOo-calc等足够了NAMD在window环境下的编译安装1.下载NAMD_2.7b2_Win322.解压到任意目录下建议最好直接是C:或D:下3.添加windows的环境变量:右键单击我的电脑-属性-高级-环境变量(在右下角)-在系统的Path变量里添加你NAMD所在文件夹,比方我的%SystemRoot%system32;%SystemRoot%;%SystemRoot%System32Wbem;C:ProgramFilesCommonFilesThunderNetworkKanKanCodecs;C:NAMD_2.7b2_Win32注意:添加的变量名称要和文件
4、夹得名称一致假如文件夹得名称你改为namd,那么变量名称为C:NAMD4.namd2.7需要后面跟conf文件才能够正确运行,并且要在conf文件所在目录执行命令。如:我的命令窗口显示C:DocumentsandSettingsHP因而我的conf文件要放在C:DocumentsandSettingsHP这个文件夹下,然后执行命令C:DocumentsandSettingsHPC:NAMD_2.7b2_Win32namd2da.conf即可。二、GROMACSGROMACS是用于研究生物分子体系的分子动力学程序包。它能够用分子动力学、随机动力学或者途径积分方法模拟溶液或晶体中的任意分子,进行分
5、子能量的最小化,分析构象等。它的模拟程序包包含GROMACS力场(蛋白质、核苷酸、糖等),研究的范围能够包括玻璃和液晶、到聚合物、晶体和生物分子溶液。GROMACS是一个功能强大的分子动力学的模拟软件,其在模拟大量分子系统的牛顿运动方面具有极大的优势。GROMACS支持几乎所有当前流行的分子模拟软件的算法,而且与同类软件相比,它还具有一些特有的优势:(1)GROMACS进行了大量的算法的优化,使其计算功能更强大。例如:在计算矩阵的逆时,算法的内循环会根据本身系统的特点自动选择由C语言或Fortran来编译。GROMACS中对Altivecloops的计算,无论是在Linux还是MacOSX.系
6、统上,它都要比其它软件快3-10倍,而且GROMACS提高计算速度的同时也保证了计算精度。(2)GROMACS具有友好的用户界面,拓扑文件和参数文件都以文档的形式给出。在程序运行经过中,并不用输入脚本注释语言。所有GROMACS的操作都是通过简单的命令行操作进行的。而且运行的经过是分步的,随时能够检查模拟的正确性和可行性,能够减少时间上的浪费。(3)GRMACS操作简单,功能丰富,而且对于初学者来讲易于上手。而且能够通过具体的免费使用手册,用户能够得到更多的信息。(4)在模拟运行的经过中,GROMACS会不断报告用户程序的运算速度和进程。(5)GROMACS具有良好的兼容性。输入文件和输出的轨
7、迹文件的格式都是独立于硬件的。(6)GROMACS能通过二进制文件来写入坐标,这样就提供了一个压缩性很强的轨迹数据存储方法,压缩方式的精度能够由用户来选择。(7)GROMACS还为轨迹分析提供了大量的辅助工具,用户不必再为常规分析编写任何程序。GROMACS还提供了轨迹的可视程序,而且很多可视化工具都能够显示。(8)GROMACS允许并行运算,使用标准的MPI通讯。(9)GROMACS程序包中包括各种常见的蛋白质和核酸的拓扑构造。包括20种标准的氨基酸以及其变异体,4种核苷和4种脱氧核苷,以及糖类和脂类。GROMACS的运行经过,主要由一系列的文件和命令组成。GROMACS一般的模拟经过能够分
8、成下面三个阶段:(1)前处理经过:生成模拟对象的坐标文件、拓扑构造文件以及平衡参数及其外力作用参数等文件。(2)模拟经过:首先要对系统进行能量最小化,避免构造的不合理而在模拟中出现错误;然后是对系统升温经过,先给系统的各个原子以Boltzmen分布初速度,再模拟较短的时间以到达初步的平衡;最后进行真正的分子动力学模拟,即平衡经过。此经过一般时间步长为1fs,运行时间在ns量级,以保证模拟系统尽可能找到势能的最低点。当然,对于其他的操作,如施加外力(模拟AFM加力)需要在平衡之后进行。在MD模拟的经过中,用户能够运用配套的可视化软件,如VMD等随时观测模拟的经过及系统的状态。(3)后处理经过:M
9、D模拟结束后,GROMACS会产生一系列文件,如.pdo文件(受力分析文件)、.trr文件(模拟经过结果文件)、.edr文件(能量文件)等。同时,GROMACS本身还提供了多种分析程序,能够对这些文件进行分析,能够得到分子体系的各种信息。2GROMACS安装流程(1)解压缩fftw,lam-mpi,gromacs源码3tar-zxvffftw-3.1.2.tar.gztarzxvfgromacs-3.3.1.tar.gztar-zxvflam-7.1.3.tar.gz(2)编译lam-mpicdlam-7.1.3./configure-prefix=/home/lam-7.1.3-withou
10、t-fc-with-rsh=ssh-xmakemakeinstall注:-without-fc是不编译mpif77,能够去掉(3)添加mpi环境变量exportPATH=$PATH:/home/lam-7.1.3/bin(appendto.bashrc)(4)编译fftw单双精度版cdfftw-3.1.2./configure-enable-float-enable-mpi-prefix=/home/fftw-3.1.2makemakeinstallmakedistclean./configure-disable-float-enable-mpi-prefix=/home/fftw-3.1.2
11、(3)设置fftw环境变量exportCPPFLAGS=-I/home/fftw-3.1.2/includeexportLDFLAGS=-L/home/fftw-3.1.2/lib(4)编译gromacscdgromacs-3.3.1./configure-prefix=/home/gromacs-3.3.1-enable-mpimakemakeinstallmakedistclean./configure-prefix=/home/gromacs-3.3.1-program-suffix=_d-enable-mpi-disable-float(5)设置gromacs环境变量exportPAT
12、H=$PATH:/home/gromacs-3.3.1/bin(appendto.bashrc)(6)编译gromacs源包里的其它文件可选makecontrib注:这步能够省去update:gromacs-4.0,fftw-3.2.1,lam7.1.4与上面的方法完全一样,只需更换目录即可。模拟步骤Belowispresentedageneralisedprocedureforperformingasimulation.Theexactstepsandprocessesinvolvedwillvarydependingonexactlywhatisbeingattempted.Useasag
13、eneralguideonly!1Clearlyidentifytheproperty/phenomenaofinteresttobestudiedbyperformingthesimulation.2Selecttheappropriatetoolstobeabletoperformthesimulationandobservetheproperty/phenomenaofinterest.Itisimportanttoreadandfamiliariseyourselfwithpublicationsbyotherresearchersonsimilarsystems.Toolsinclu
14、de:-softwaretoperformthesimulationwith,considerationofforcefieldmayinfluencethisdecision.-forcefieldwhichdescribeshowtheatoms/particleswithinthesysteminteractwitheachother.Selectonethatisappropriateforthesystembeingstudiedandtheproperty/phenomenaofinterest.Veryimportantandnon-trivialstep!3Obtain/gen
15、eratetheinitialcoordinatefileforeachmoleculetobeplacedwithinthesystem.4Generatetherawstartingstructureforthesystembyplacingthemoleculeswithinthecoordinatefileasappropriate.Moleculesmaybespecificallyplacedorarrangedrandomly.5Obtain/generatethetopologyfileforthesystem,using(forexample)pdb2gmx,PRODRGor
16、yourfavouritetexteditorinconcertwithchapter5oftheGROMACSManual.6Describeasimulationbox(e.g.usingeditconf)whosesizeisappropriatefortheeventualdensityyouwouldlike,fillitwithsolvent(e.g.usinggenbox),andaddanycounter-ionsneededtoneutralizethesystem(e.g.usinggromppandgenion).Inthesestepsyoumayneedtoedity
17、ourtopologyfiletostaycurrentwithyourcoordinatefile.7Runanenergyminimisationsimulationonthesystem(usinggromppandmdrun).Thisisrequiredtosortoutanybadstartingstructurescausedduringgenerationofthesystem,whichmaycausetheproductionsimulationtocrash.8Selecttheappropriatesimulationparametersfortheequilibrat
18、ionsimulation(definedin.mdpfile).Youneedtobeconsistentwithhowforcefieldwasderived.YoumayneedtosimulateatNVTwithpositionrestraintsonyoursolventand/orsolutetogettheTalmostright,thenrelaxtoNPTtofixthedensity,thenmovefurther(ifneeded)toreachyourproductionsimulationensemble(e.g.NVT,NVE).9Runtheequilibrat
19、ionsimulationforsufficienttimesothatthesystemrelaxessufficientlytoallowtheproductionruntobecommenced(usinggromppandmdrun,theng_energyandtrajectoryvisualisationtools).10Selecttheappropriatesimulationparametersfortheproductionsimulation(definedin.mdpfile),inparticularbecarefulnottore-generatetheveloci
20、ties.Youstillneedtobeconsistentwithhowtheforcefieldwasderivedandhowtomeasuretheproperty/phenomenaofinterest.11Runtheproductionsimulationforsufficienttimesothatproperty/pheno-menaofinterestcanbeobservedinrequireddetail(usinggrompp/tpbconvandmdrun).12Analyse/visualisetheresultingtrajectoryanddatafiles
21、toobtaininformationontheproperty/phenomenaofinterest.4三、AmberAmber是著名的分子动力学软件,用于蛋白质、核酸、糖等生物大分子的计算模拟。Amber也指一种经历力场(empiricalforcefields)。力场和代码是分开的,一些软件中包含amber力场,而其他的力场也包含在此amber的软件中。AMBER提供两部分内容:用于模拟生物分子的一组分子力学力场无版权限制,也用于其它一些模拟程序中;分子模拟程序软件包,包含源代码和演示有版权限制,需要购买。1AMBER主要程序Leap:用于准备分子系统坐标和参数文件,有两个程序:xleap:Xwindows版本的leap,带GUI图形界面tleap:文本界面的LeapAntechamber:用于生成少见小分子力学参数文件的。有的时候一些小分子Leap程序不认识,需要加载其力学参数,这些力学参数文件就要antechamber生成Sander:MD数据产生程序,即MD模拟程序,被称做AMBER的大脑程序。Ptraj:MD模拟轨迹分析程序。