分享CADD之分子动力学的简介(下篇)
在上篇中,我们了解了分子动力学的基础和基本原理,接下来我们来看一下MD模拟是如何保证准确性的,以及它的主要流程。
MD的准确性如何?
通常一种理论方法的出现,会伴随各种声音的出现,MD模拟也毫无例外,有支持派,当然也有质疑其准确性的。那它对分子结构和原子间相互作用描述的到底准确吗?
在上篇当中提到的范德华相互作用,其是通过LJ势进行描述的,该势中除了有与位置坐标相关的r之外,LJ势中还含有两个参数:σ和ε。以AMBER99SB力场为例,在该力场下脂肪族中的C和羰基中的C,这两个参数的值会不同。当然在不同的力场下,同种原子的这些参数也会有差别。那么这类参数是如何得到的呢?这些参数是拟合自实验数据或量子化学的结果,是通过选定初始参数,做分子模拟,然后看模拟结果和实验值的差别,重新调整参数,继续模拟,直到达到收敛标准。这一过程属于经验描述 ,都是为了和实验数据相符合,这也就是为什么会出现四点水模型和五点水模型的原因。而且各种力场也会根据新的实验数据去做参数优化,使其更加准确。现如今已发展了很多分子力场,比如生物模拟 常用的AMBER , CHARMM , OPLS , GROMOS ,材料领域 常用的CFF , MMFF , COMPASS 等等。
图1是与MD模拟相关的文献发表量(从1977-2017年),并且这些文章发表在顶级期刊。这表明MD模拟在近年来已经变得越来越常见,也逐渐受到了认可。
图1 与MD模拟相关的文献量
MD如何应用?
应用软件
随着计算机性能的发展,能执行MD模拟的软件也逐渐增多。现如今比较主流的MD模拟软件(见表 1),如AMBER ,GROMACS ,CHARMM 等。当然也有一些其他可以进行MD模拟的计算软件,但由于只是该软件的一个模块,计算十分耗时,并不推荐。
表 1 几款常用的MD模拟软件
对于初学者而言,GROMACS是相对来说较为友好的软件。这款软件是开源的,也支持GPU加速,安装简单,很容易上手。鉴于GROMACS官网上已经有很多分子动力学相关教程(包括蛋白、蛋白-小分子、蛋白-磷脂分子层、虚拟位点等),我们这里就不重复介绍已有的蛋白模拟教程,小伙伴们可以在文末找到教程网址,自行操作。
主要过程
无论是哪一种软件,它们的步骤都是相类似的,大致包含以下过程(图2):
图2 主要过程
初始化
初始化是分子动力学里重要的一步,它会读取模型的参数,模拟控制参数。如增加非重原子的缺失;确定溶剂;确定周期性边界条件盒子;确定力场等。(选择恰当的力场是保证体系结果合理的前提,因此力场选择需谨慎! )
能量优化
通常有些结构模型会存在原子重叠(图3),键长和键角的严重扭曲等一些不合理的状态,若以这种状态直接进行分子动力学模拟,得到的数据也是不准确的。为了解决这一类问题,通常会对模拟体系进行能量优化,优化其模型结构。一般会优化几千到几万步。
图3 原子重叠
平衡
平衡分子动力学模拟,总是在一定的系综下进行的,系综 大致可分为以下几类:(系综:代表一定条件下,一个体系的大量可能状态的集合 )
微正则系统: 系统原子数N,体积V,能量E保持不变。又称为 NVE系综。微正则系综里的每个体系具有同等的能量。
正则系统 : 系统原子数N,体积V,温度T保持不变,且总动量保持不变。又称为NVT系综。该系综里的各体系可以和其他体系进行交换能量。
等温等压系统: 系统原子数N,压强P,温度T保持不变,又称为NPT系综。在等温等压系综下,各体系可以和其他体系交换能量和体积,但系综内各个体系有相同的温度和压强。
等压等焓系统: 系统原子数N,压强P,焓值H=E+PV保持不变。在模拟中较少见。
这一步是为了保证体系的稳定,否则在执行后面的操作时,体系会发生崩溃。
数据产出: 前几个步骤都是准备工作,为的就是最后的数据产出这一步。数据产出的这一过程就是之前提到的基本原理中力场和势能的应用,要在整个系统平衡后进行。该过程可记录下体系中粒子随时间的变化的坐标、速度和能量。为了可观测到研究体系的性质,并保证此过程具有可重复性,模拟的时间(步长)一定要够久。
当完成以上步骤后,就可以对体系的轨迹进行分析处理,观察体系具有什么现象或性质。
总结
随着算力不断地提升,计算小体系下的毫秒级别和百万原子体系已经可以实现。而且MD模拟也衍生出来更多分支,如增强采样的RMED;适用于大分子体系的粗粒化模拟;提升计算速度的DMD等。MD模拟的参数是由经验参数组成的;模拟可以提供详细的构象分布和时间序列,是对实验的一种补充。所以说这是实验支撑理论,理论补充实验的一个过程。当然经验势函数也是有局限性的,针对不同体系的准确性会有差异。为了提升势函数的精确性及计算效率性,研究人员也在不断改进完善算法。
希望该文能让初学者对分子动力学模拟有初步的了解,欢迎大家一起讨论CADD的相关内容!
GROMACS教程网址 http://www.mdtutorials.com/gmx/
参考文献
[1]. Hollingsworth SA, Dror RO. Molecular Dynamics Simulation for All[J]. Neuron. 2018 Sep 19;99(6):1129-1143.
[2]. 李文飞, 张建, 王骏, 王炜. 生物大分子多尺度理论和计算方法[J]. 物理学报, 2015, 64(9): 098701.
[3]. Wang Y, Liu T, Xie J, Cheng M, Sun L, Zhang S, Xin J, Zhang N. A review on application of molecular simulation technology in food molecules interaction[J]. Curr Res Food Sci. 2022 Oct 11;5:1873-1881.
[4]. 王存新. 蛋白质模拟:原理,发展和应用[M]. 科学出版社, 2016.
[5]. 刘冠辰.分子动力学模拟及在生物大分子模拟领域的应用[J].吉林化工学院学报,2015,32(11):112-116.
[6]. 刘远超. 分子力场[EB/OL].https://blog.sciencenet.cn/blog-3330179-1197033.html.
[7]. Yunxiang Sun, Aleksandr Kakinen, Yanting Xing, Pouya Faridi, Aparna Nandakumar, Anthony W. Purcell, Thomas P. Davis, Pu Chun Ke and Feng Ding. Amyloid Self-Assembly of hIAPP8-20 via the Accumulation of Helical Oligomers, α-Helix to β-Sheet Transition, and Formation of β-Barrel Intermediates[J]. Small, 2019: e1805166.
[8]. Hess Berk, Kutzner Carsten, Spoel David Van Der, et al. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation [J]. Journal of chemical theory and computation, 2008,4(3): 435-447.
人生的意义也许永远没有答案,但也要尽情感受这种没有答案的人生一诗一文每一天的原创诗歌和梦想,不要再和我擦肩而过无论世界如何改变,是繁荣或又要握紧长剑,我在,你便只管往前。云兮人不应该是插在花瓶里供人观赏的静物,而是蔓延在草原上随风起舞的韵律
感情总是不顺人际关系不好,是因为过去的伤每个孩子都需要一个地方,在这里,他感到安全,并受周围的人喜爱。每个人都向往一个地方,在这里,他可以彻底放松,完全成为他自己。理想情况下,这个地方应该是和父母一起生活的原生家庭。在原
气温大跳水!江西人,你的降温套餐即将发货本周的晴好天气余额已经不多了一波冷空气已经蓄势待发本周末即将到货全国天气本周大范围雨雪即将上线11日至15日西北地区东部华北东北地区大部黄淮江淮江汉西南地区东部江南北部等地将有小到
网易元宇宙业务又添新动态,合作央视网,定制春晚元宇宙会场今日,网易在元宇宙业务上又添新动态,其宣布将和央视网合作,为网络春晚观众定制首个元宇宙会场,来实现让上万人实时互动。近一月来,网易还与茅台合作推出元宇宙相关的APP,投资了张艺谋联
东威科技业绩快报2022年净利同比增34。3新能源镀膜设备收入同比增1462东威科技业绩快报2022年净利同比增34。3新能源镀膜设备收入同比增1462科创板日报11日讯,东威科技公布业绩快报,公司2022年度实现营业总收入10。18亿元,同比增长26。5
2023渴望已久,向我们的星辰大海进军吧总有一些什么值得我们去热爱总有一些什么值得我们去期待总有一些什么值得我们去相信总有一些什么值得我们去坚守总有一些什么值得我们去奔赴总有一些什么值得我们去成就2023年,日子是全新的
浅谈我的2022年(原创首发)我的生活也是头条在过去的2022年,世界上发生了众多的事情,最让世人瞩目的是新冠疫情和俄乌战争。这些在各种媒体上都可以看到,我不多说。我只说我在过去的一年里取得的微不足道的进步。四
与家父之不相见,已逾二十多个年头了,但近来却屡屡在梦中念起他家父相里与家父之不相见,已逾二十多个年头了。但,近来却屡屡在梦中念起他。与此同时,梦幻中出现的,还有我的母亲。此乃今年以来,总在脑际萦绕常在心头纠结的事了。有关母亲的事,我曾在早前
2023最好的打开方式1重启自己新的一年已然到来。好事坏事,皆成往事,昨年的成就和遗憾都一并归零。请重启自己,然后大步向前。人生的主动权,永远把握在那些敢于重新出发的人手中。正如阅读书籍,只有一页一页地
雪落古凉州今夜的凉州,又落雪了。三分羞涩,七分自觉。而前些日子款款来过的一场初雪,没有今夜这般酣畅快意。那场初雪,怯怯地来,轻如柳絮,飘洒了半宿,便依依地走了。对凉州而言,一花一世界,一雪一
金句摘抄1。多难殷忧兴国运,动心忍性希前哲。2。生逢盛世,我们享受着新时代的和平阳光,书写着新时代的奋斗华章。3。雷霆万钧壮志不改,沧海横流方显英雄气概。4。真相不是一块橡皮泥,可以随意揉