【摘要】 掌握GROMACS中的RMSD分析流程，理解参考结构、质量权重、Protein组选择和RMSD平台期判断。

科研人必备干货｜全套科研知识库

常用科研资料汇总，点击领取即可获取。

 

RMSD（root mean square deviation，均方根偏差）是分子动力学后处理中最常用的结构变化指标之一。它通过比较轨迹中每一帧结构与参考结构的差异，帮助研究者判断蛋白质骨架是否趋于稳定，以及模拟过程中结构偏移的幅度。

 

RMSD到底测量什么

RMSD回答的是一个具体问题：轨迹当前帧与参考结构经过对应原子匹配后，平均偏移了多少。数值越小，说明当前结构与参考结构越接近；曲线是否趋于平稳，则比某个孤立数值更适合用于观察结构变化趋势。

RMSD常用于：

  • 判断蛋白质主链是否在模拟过程中趋于稳定；

  • 比较不同阶段的构象偏移；

  • 识别明显的结构跃迁或持续漂移；

  • 为后续RMSF、氢键和构象聚类分析提供基础判断。

 

RMSD公式与质量权重

RMSD本质上是当前结构与参考结构之间均方偏差的平方根。实际计算时，需要先确定原子选择、参考结构、是否进行拟合以及是否采用质量权重。

图：RMSD定义公式。

质量权重版本会按照原子质量调整不同原子对结果的贡献。如果不希望使用质量权重，可在GROMACS命令中加入 -nomw。

图：质量权重RMSD的表达形式。

 

GROMACS标准操作流程

1. 准备生产轨迹

在进行RMSD分析前，应确认生产模拟已经完成，并检查轨迹、结构输入和拓扑文件是否对应。来源教程使用了 prod.xtc 和 prod.tpr 作为示例。

2. 运行gmx rms

gmx rms -f prod.xtc -s prod.tpr -o rmsd.xvg

命令参数含义如下：

参数

含义

gmx rms

调用RMSD分析模块

-f prod.xtc

指定压缩轨迹文件

-s prod.tpr

指定模拟输入结构

-o rmsd.xvg

输出RMSD数据文件

在交互选择中，通常需要分别指定拟合组和RMSD计算组。蛋白质体系的基础分析可以选择 Protein，但对于膜蛋白、配体或核酸体系，应根据研究问题选择更具体的原子组。

3. 绘制RMSD曲线

xmgrace rmsd.xvg

也可以使用Python、Origin或其他绘图工具读取XVG文件并统一图表样式。

 

如何判断曲线是否趋于稳定

RMSD曲线整体上升，通常意味着体系仍在经历构象变化或从初始结构中偏离；曲线在一定时间后围绕平均值波动，说明结构变化趋于稳定；出现持续漂移、多次跳变或多个平台，则需要结合轨迹可视化、RMSF、聚类和能量指标进一步分析。

图：RMSD曲线用于观察结构偏移与平台期。

重要的是，RMSD平台期只能说明所选原子组的整体偏移趋于稳定，不代表体系所有性质都已经收敛。例如，局部侧链、配体姿态、溶剂结构和氢键网络可能仍在变化。

 

2 ns是否足够

2 ns可以作为教程演示量级，但不能直接作为所有蛋白质体系的模拟标准。是否足够取决于体系大小、初始构象、温度和压力耦合、目标性质以及多个独立重复轨迹的一致性。对于论文研究，应在RMSD趋稳后继续采样，并结合其他收敛指标判断统计可靠性。

 

常见问题

RMSD和MSD是同一个指标吗？

不是。RMSD通常描述结构相对参考结构的偏移，常用于蛋白质构象分析；MSD则常用于描述粒子或离子随时间的位移，并可进一步估算扩散行为。两者都含有“平方位移”思想，但研究对象和解释方式不同。

RMSD曲线变平就说明模拟完成了吗？

不一定。RMSD变平说明选定原子组的整体结构偏移趋于稳定，但还需要结合RMSF、能量、半径、氢键、构象聚类或目标性质的时间序列进行判断。

 

结语

GROMACS RMSD分析的关键不只是运行一条命令，而是明确参考结构、原子选择、拟合方式、质量权重和收敛判断标准。将RMSD曲线与轨迹观察及其他后处理指标结合，才能把“曲线变平”转化为更可靠的结构稳定性结论。