|||
以下所有内容均属于个人学习过程中的总结,如有错误,欢迎批评指正!
用Amber进行能量分解和计算结合自由能
First release:2018-01-13 Last update: 2018-03-17
Amber是一款适用于生物大分子的动力学模拟的软件,其官方网站是http://ambermd.org/index.html,目前最新版本为Amber17,Amber软件包主要包括2个部分:Amber Tools和Amber,其中Amber Tools可以免费在Amber官网下载,Tools中包含了Amber几乎所有的模块,Sander、tleap、MMPBSA等最核心的内容都可以免费使用。而另外的Amber则是唯一收费的部分,该部分主要包括了PMEMD以及GPU加速 的功能,简而言之, 如果你不追求计算速度,完全可以用免费的Amber(速度是相当慢,16线程的机器,跑2~3万原子的体系,一天大概是2~3 ns),而Amber加速之后,速度提升10倍以上。
Amber的能量分解部分是非常重要的模块,即MMPBSA,能量分解其实就是计算配体和受体之间结合的亲密程度,主要思想是将驱动二者结合的能量分解到每一个氨基酸残基上,我们可以清楚的看到每个氨基酸对配体结合的具体能量贡献,包括VDW、溶剂化能、静电能等。下面笔者简单介绍一下,用MMPBSA.py快速计算结合自由能的方法,希望对正在学习Amber的新手朋友带来帮助。
MMPBSA.py的计算过程只有一步,以前的MMPBSA.pl脚本需要多步操作,非常繁琐。输入文件mmpbsa.in内容如下:
Input file for running entropy calculations using NMode
&general
startframe=395, endframe=490, interval=5, keep_files=2, ########## startframe=100, endframe=200分别表示用于计算能量分解的起始构象和终止构象;interval表示每隔5个构象提取一个,一共是20个构象;keep_files关键词是保存临时文件,0为不保存,1为保存所有临时文件,2为保存临时文件并在计算成功结束后自动删除。
receptor_mask=:185-404, ligand_mask=:1-184, ########## receptor_mask=:2-347, ligand_mask=:1-1,分别表示定义的受体和配体的序号,这个需要必须与动力学模拟时的序号一致(在计算蛋白和蛋白的相互作用时,可以将ligand_mask的值改为对应的残基序号即可)
/
&gb
igb = 5, saltcon = 0.15, ########## igb:Generalized Born method to use,可以设置为1,2,5,7,默认值为5 ;saltcon为盐浓度,单位mol/L
/
&pb
istrng = 0.15, ########## istrng表示离子强度,在PBSA计算中为mmol/L
/
以上参数设置参考自Amber说明书,MMPBSA.py的参数设置还有很多,但上述参数基本上能满足大多数需求,具体的参数需要根据不同的实验条件来进行设置和修改,请参考Amber说明书中的31. MMPBSA.py,page 660
上述mmpbsa.in文件为参数输入文件,计算结合自由能还需要三个拓扑文件,complex.prmtop,receptor.prmtop和ligand.prmtop,这些拓扑文件在做动力学模拟之前就已经做好,以及最终的轨迹文件(建议将轨迹文件采用cpptraj工具去除水分子和其他离子),md1.mdcrd
准备好上述文件后,执行命令:MMPBSA.py -O -i mmpbsa1.in -o mmpbsa.dat -cp complex.prmtop -rp receptor.prmtop -lp ligand.prmtop -y md1.mdcrd > progress.log,就可以进行计算了。结果文件中会给出三个体系(complex、receptor和ligand)的各项能量,以及总能量。结合自由能=receptor+ligand-complex
如果你对Amber分子动力学模拟或者量子化学感兴趣,或者对本文有什么问题,欢迎加入我们的交流群:
qq群 580744615
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-10 06:29
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社