【摘要】 从RMSD平台期、持续漂移、多峰曲线到参考结构选择，系统讲解分子动力学模拟中RMSD结果的判断和常见问题。

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

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

 

RMSD曲线不是简单的“越低越好”。正确的解读要同时考虑参考结构、拟合组、计算组、时间窗口和体系本身的构象变化。曲线趋于平台，通常意味着整体结构偏移趋于稳定；持续上升或反复跳变，则需要进一步排查。

 

四类常见曲线形态

1. 快速上升后进入平台

这通常表示体系从初始构象调整到新的平衡附近。若平台期波动范围稳定，并且能量、温度和其他结构指标也合理，可以将其作为进一步分析的候选区间。

2. 全程持续上升

可能原因包括体系仍在发生构象转变、初始结构与模拟条件不匹配、轨迹预处理不正确，或所选原子组包含了持续移动的柔性区域。

3. 多次跳变并形成多个平台

这可能对应不同构象状态之间的转换。此时不应简单删除跳变区间，而应结合轨迹可视化、聚类分析和功能状态判断其物理意义。

4. 高频大幅波动

需要检查体系是否发生PBC跳跃、拟合组是否过小、轨迹采样是否不足，以及单位和时间轴是否正确。

图：RMSD平台期和结构漂移需要结合曲线整体形态判断。

 

参考结构怎么选

常见参考结构包括第一帧、能量最小化后的结构、平衡阶段末帧或实验结构。不同参考结构会改变RMSD的绝对数值，因此比较不同模拟体系时，应尽量保持参考结构和原子选择一致。

如果研究关注的是蛋白质主链稳定性，可以使用主链或Cα作为计算组；如果关注配体结合姿态或特定结构域，则应设置针对性原子组，并明确是否先对整体蛋白进行拟合。

图：RMSD计算依赖当前结构与参考结构之间的对应关系。

 

质量权重什么时候使用

质量权重会让不同质量原子对RMSD的贡献不同。对于常规蛋白质主链分析，是否使用质量权重应在方法部分明确说明，并在不同体系比较时保持一致。如果需要不考虑质量权重，可在GROMACS命令中加入 -nomw。

 

一个可复核的GROMACS命令

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

建议在记录中同时保存：

  • 拟合组和RMSD计算组；

  • 参考结构或参考帧；

  • 是否使用质量权重；

  • 轨迹起止时间和步长；

  • GROMACS版本；

  • 轨迹是否经过PBC处理和居中。

 

RMSD不能单独证明什么

RMSD不能单独证明模拟已经完全平衡，也不能直接证明蛋白质功能、结合自由能或实验活性。更稳妥的做法是将RMSD与RMSF、回转半径、氢键、二级结构、能量和目标性质的时间序列结合起来。

 

2 ns演示与论文级模拟的区别

教程中的2 ns有助于展示命令和曲线，但论文级模拟所需时间取决于体系和研究问题。对于构象变化较慢的体系，应在曲线趋稳后继续采样，并通过独立重复、分块平均或相关时间分析评估统计可靠性。

图：模拟前的蛋白质初始模型。

图：演示体系完成2 ns模拟后的结构示意。

 

FAQ

RMSD平台期应该从哪一帧开始统计？

没有适用于所有体系的固定帧数。应先观察曲线，再结合其他指标确定平台区间，并在方法中说明排除的预平衡阶段和最终统计窗口。

RMSD异常时最先检查什么？

优先检查轨迹是否发生PBC跳跃、结构是否正确居中、参考结构是否匹配、拟合组是否合理，以及时间单位和输出文件是否正确。确认这些基础问题后，再判断异常是否来自真实构象变化。

 

结语

RMSD的价值在于把结构变化转化为可量化的时间序列，但它必须放在完整的模拟质量控制流程中解释。命令、原子组、参考结构和统计窗口记录得越清楚，结果越容易复核，也越适合进入论文方法和补充信息。