Gromacs——教程学习(6)

  1. 谈谈怎么判断分子动力学模拟是否达到了平衡

在计算RMSD之前必须先通过最小二乘法将各帧结构相对于参考结构进行最大程度叠合,从而消除体系的整体运动而令RMSD只体现生物分子内部结构的变化,这称为align或者least squares fit。

需要注意的是,蛋白质骨架的RMSD曲线并不总能反映出蛋白结构特征变化的全貌。一方面,RMSD的计算涉及到一大堆原子的平均,若一个很大的蛋白某个局部出现一定程度的显著的、不可忽视的结构变化,在RMSD计算过程中它会被很大程度淹没掉。

另一方面,骨架的RMSD不能反映侧链上出现的一些关键变化。例如上面这篇文章发现,在423 ps的时候磷酸化酪氨酸574残基侧链出现了显著的翻转,这使得被研究的酪氨酸激酶活性位点打开,这在分子生物学角度是一个很重要的现象,而从alpha碳的RMSD曲线上则完全体现不出来这一点。文中图4(A)绘制了这个残基侧链二面角随时间的变化,这才让这关键信息被展现出来。

如果蛋白质的RMSD曲线在已经跑的模拟时间内还没有稳定在一个值附近且持续了较长时间(有零点零几nm的上下波动是可以接受的),而你又需要蛋白质已达到平衡的结构,通常就需要延长模拟时间续跑,直到RMSD基本平稳。

有些人试图用衡量小分子液体的方式衡量生物分子模拟体系是否达到稳定,即还用温度、密度、势能等参数来衡量,这明显不行。因为这些量的收敛速度远快于蛋白质的RMSD的收敛,在这些指标已平稳的时候蛋白质部分通常还明显没达到平衡。虽然蛋白质整体结构的进一步变化必然会对势能等因素有多多少少的影响,但在势能曲线上基本体现不出来。通常模拟蛋白质都是在显式水模型下模拟的,势能绝大部分都是水所贡献的,显然蛋白质结构的改变对总势能的影响通常会被淹没;

下面两张图是计算化学公社论坛上我答疑时看到的,可见中途RMSD都出现过剧烈的变化,完全超越了平衡状态下的正常波动的程度。这种突变必须给予明确的解释,要不然审稿人看见肯定会问你是怎么回事。始终要记住,RMSD曲线是轨迹变化的定量反映,想解释抽象的RMSD曲线的突然变化是什么导致的,当然要去用VMD程序看轨迹动画才能充分弄清楚原因。绝对不可能RMSD曲线出现厉害的突变时在轨迹动画上却看不出任何端倪。

在这里插入图片描述在这里插入图片描述
上图这种情况只要在VMD里观看轨迹动画自然能明白原因。就算不看轨迹也100%能断定肯定是因为周期性方式记录轨迹,而被计算RMSD的对象在模拟过程中越过了盒子所致。当某一帧突然体系被弄到了盒子另一头,RMSD自然会出现瞬间的大幅变化。解决这种RMSD曲线不连续性的方法很简单,消除被计算RMSD的对象的整体运动,并保持此对象的完整性就完了。

比如在GROMACS里,可以用gmx trjconv结合-pbc mol使得转换出的新轨迹中的分子保持完整,计算RMSD前再经过叠合来消除蛋白质的整体运动,其RMSD曲线就不可能突变。还有一种情况是被计算RMSD的对象是由多个分子组成的,比如有多条链的蛋白、蛋白+配体复合物、含有两条链的核酸,这种情况可能模拟中途只有其中一部分跑出盒子而另一部分还在盒子内,光是用-pbc mol保持单分子完整的话此对象的不同部分还是分家的,对这种情况需要用gmx trjconv结合-pbc cluster让指定的组在轨迹中一直保持完整,组的定义就对应于你要计算的RMSD的对象。顺带一提,对生物分子的模拟,我一般建议让程序在模拟过程中就一直对生物分子消除平动转动,对于GROMACS跑蛋白质来说就是设comm-grps = protein和comm-mode = angular,这样就避免了生物分子中途跑出盒子的问题。

不同的生物分子体系的RMSD达到平稳时RMSD的数值会有一定不同。柔性越大的体系,平衡时RMSD越大。比如DNA体系骨架柔性比一般较刚性的蛋白质要大,这会体现在达到平衡时RMSD相对较大上。

通过RMSD考察是否平衡不是必须得用初始帧做参考结构。J. Chem. Phys., 109, 10115 (1998)认为以初始结构当参考帧来观察RMSD会高估达到平衡的时间,因此会把一部分其实已经达到平衡的轨迹误认为还没平衡,导致浪费了一些轨迹。作者的观点是相对于初始结构,模拟开始后RMSD的增加由两部分组成:(1)非平衡态到达平衡态 (2)平衡态时在相空间采样。实际上(1)满足了的时候本质上生物分子就已经平衡了,这比起以初始结构为参考帧计算的RMSD曲线达到平稳时更早。为了判断什么时候就已经满足了(1),作者建议以已经平衡了的结构作为参考帧来计算RMSD。下图是文中的图,作者取500、600、700、800、900 ps分别作为参考帧绘制了RMSD曲线,并假定这些时刻体系都已经达到了平衡。由图可见,每条曲线在参考帧附近都是迅速上升,这体现出已平衡的结构在相空间采样期间由于热运动自然而然地相对于参考结构发生的结构改变;而在图的前200 ps,这几个曲线都明显上升,因此体现出前200 ps是确确实实结构没有达到平衡、有整体性变化的阶段。
在这里插入图片描述
所以,像温度、势能、密度那些很容易达到平稳的量,只是体系达到平衡的必要而非充分条件,不能只拿这些量说事,除非是本文第1节说的均匀的小分子液体之类的简单体系。
势能只是能量角度上的一个标准,但凡有可能,最好还是再结合一个结构方面的标准来进一步判断平衡,明显会更有说服力。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/bicheng/4822.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

国产Sora诞生!清华团队发布Vidu大模型,可直接生成16秒视频

大模型之争已从单模态转向多模态。 4月27日,在2024中关村论坛年会未来人工智能先锋论坛上,清华大学联合北京生数科技有限公司正式发布了文生视频大模型——Vidu。 在会议上,清华大学人工智能研究院副院长、生数科技首席科学家朱军对外展示了…

stack,queue的模拟实现以及优先级队列

这篇博客用来记录stack&#xff0c;queue的学习。 stack的模拟实现 stack的模拟实现比较简单&#xff0c;先上代码 #pragma once #include<vector> #include<list> #include<deque> #include<iostream> using std::deque; using namespace std;name…

PyTorch中Torch.arange()函数详解

函数原型 arange(start0, end, step1, *, outNone, dtypeNone, layouttorch.strided, deviceNone, requires_gradFalse) -> Tensor 用法 返回大小为的一维张量&#xff0c;其值介于区间 为步长等间隔取值 参数说明 参数类型说明startNumber起始值&#xff0c;默认值&…

【Python的魅力】:利用Pygame实现游戏坦克大战——含完整源码

文章目录 一、游戏运行效果二、代码实现2.1 项目搭建2.2 加载我方坦克2.3 加载敌方坦克2.4 添加爆炸效果2.5 坦克大战之音效处理 三、完整代码 一、游戏运行效果 二、代码实现 坦克大战游戏 2.1 项目搭建 本游戏主要分为两个对象&#xff0c;分别是我方坦克和敌方坦克。用户可…

C++ 中的隐式类型转换

基本数据类型的转换(不是重点) int i 42; double d i; // int 自动转换为 double 类类型转换: 如果一个类定义了接受单一参数的构造函数&#xff08;且该构造函数未被声明为 explicit&#xff09;&#xff0c;那么该构造函数可以被用来进行隐式转换,可以理解为int被隐式转…

告别盲目投放,Xinstall让App广告效果一目了然

在移动互联网时代&#xff0c;App广告的投放已经成为企业营销的重要手段。然而&#xff0c;广告主在投放广告后&#xff0c;如何有效追踪广告效果、衡量广告ROI&#xff0c;一直是困扰他们的难题。今天&#xff0c;我们就来聊聊Xinstall这家国内专业的App全渠道统计服务商&…

计算一个结构的平方

在行列可自由变换的平面上计算3a1*3a1 得到的结构应该有9个点&#xff0c;并且点的分布有3a1的特征 这个9点结构可以化简成 1 1 1 1 1 1 1 1 1 显然这个结构可以再次平方&#xff0c;得到 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 …

python:reportlab 生成pdf:基本用法。

1.首先&#xff0c;打开cmd&#xff0c;安装reportlab pip install -i https://pypi.tuna.tsinghua.edu.cn/simple reportlab #从清华镜像安装更快 然后就可以使用其基本用法。 from reportlab.lib.pagesizes import letter from reportlab.pdfgen import canvasdef genera…

Day43 JDK1.8新特性 下 接口的类优先原则和接口冲突,日期时间组件,重复注解

Day43 JDK1.8新特性 1.接口的默认方法与静态方法 从JDK1.8开始&#xff0c;接口中可以有默认方法&#xff0c;既default修饰的方法&#xff0c;此方法可以让接口的实现类所调用&#xff0c;而接 口中的静态方法直接用接口名调用即可 public class Test1 {Testpublic void test…

小程序评分/关键词/UV优化助力小程序登顶

随着小程序市场的日益繁荣&#xff0c;小程序搜索排名优化成为了众多开发者关注的焦点。小程序搜索排名被很多因素影响着&#xff0c;关键词、评分还有uv&#xff08;授权&#xff09;等。在本文小柚和各位老板分享如何有效优化小程序搜索排名的经验。 一、关键词策略 关键词是…

Python | Leetcode Python题解之第55题跳跃游戏

题目&#xff1a; 题解&#xff1a; class Solution:def canJump(self, nums: List[int]) -> bool:n, rightmost len(nums), 0for i in range(n):if i < rightmost:rightmost max(rightmost, i nums[i])if rightmost > n - 1:return Truereturn False

【c++】cpp运算符重载

目录 &#xff08;1&#xff09;什么是运算符重载 &#xff08;2&#xff09;运算符重载的本质是函数调用 &#xff08;3&#xff09;可以与不可以重载的运算符 &#xff08;4&#xff09;单目运算符与双目符重载区别 &#xff08;5&#xff09;双目运算符重载举例 重载操…

Web3的可持续性:构建环境友好的去中心化系统

引言 随着全球对可持续发展和环境问题的日益关注&#xff0c;Web3技术作为一种新型的互联网模式&#xff0c;也开始受到社区和开发者的关注。但很少有人关注到Web3对环境可持续性的潜在影响。本文将探讨Web3如何构建一个环境友好的去中心化系统&#xff0c;以及这如何促进一个…

偏微分方程算法之五点菱形差分法

目录 一、研究目标 二、理论推导 三、算例实现 四、结论 一、研究目标 上个专栏我们介绍了双曲型偏微分方程的主要算法及实现。从今天开始&#xff0c;我们在新的专栏介绍另一种形式偏微分方程-椭圆型的解法。 研究目标选取经典的二维椭圆型方程&#xff08;也称泊松Poisso…

马斯克突击访华;谷歌 Python 基础团队全数被裁;丨 RTE 开发者日报 Vol.195

开发者朋友们大家好&#xff1a; 这里是 「RTE 开发者日报」 &#xff0c;每天和大家一起看新闻、聊八卦。我们的社区编辑团队会整理分享 RTE&#xff08;Real Time Engagement&#xff09; 领域内「有话题的 新闻 」、「有态度的 观点 」、「有意思的 数据 」、「有思考的 文…

python基础学习之写入csv文件

前言 在Python编程中&#xff0c;经常会遇到要将数据存到csv文件中&#xff0c;今天来详细讲一下。 一.功能目的 将以下数据存到csv文件中。 data [ {name: Alice, age: 25, city: New York}, {name: Bob, age: 30, city: Los Angeles}, {name: Charlie, age: 35, city:…

共享模型之无锁——Unsafe

文章目录 概述Unsafe CAS 操作 名字虽然叫Unsafe,但并不是线程不安全&#xff0c;而是因为他会操作内存&#xff0c;操作线程&#xff0c;不建议开发人员使用。 概述 Unsafe 对象提供了非常底层的&#xff0c;操作内存、线程的方法&#xff0c;Unsafe 对象不能直接调用&#x…

群狼调研(长沙商业咨询)广告效果测评的关键指标

广告效果测评涉及多个关键指标&#xff0c;以下是其中一些常见的指标&#xff1a; 1.广告到达&#xff08;Ad Reach&#xff09;&#xff1a;衡量广告在目标受众中的覆盖范围和到达程度。它可以包括广告曝光的人数、频次、覆盖率等方面的评估。 2.广告认知&#xff08;Ad Aware…

基于Python和Selenium的BOSS直聘Python岗位数据分析系统的设计与实现

基于Python和Selenium的BOSS直聘Python岗位数据分析系统的设计与实现 Design and Implementation of Python-based Selenium-powered BOSS Direct Recruitment Python Job Data Analysis System 完整下载链接:基于Python和Selenium的BOSS直聘Python岗位数据分析系统的设计与实…

H5 录音功能

Recorder: html5 js 录音 mp3 wav ogg webm amr g711a g711u 格式&#xff0c;支持pc和Android、iOS部分浏览器、Hybrid App&#xff08;提供Android iOS App源码&#xff09;、微信&#xff0c;提供ASR语音识别转文字 H5版语音通话聊天示例 DTMF编码解码 git 地址&#xff1a…