分子动力学结果分析

分子动力学结果分析是把一段轨迹翻译成可发表的物理解释的关键环节。本项目在交付 MD 模拟后， routinely 输出 RMSD、RMSF、氢键、溶剂可及面积与自由能分解等成套指标。下文以一条 100 ns 蛋白-配体轨迹为例，说明各指标该看什么、怎么判读、常见陷阱在哪里。

先确认轨迹是否”收敛”

拿到轨迹第一件事不是看结合模式，而是确认体系是否平衡。本项目习惯先画 Cα-RMSD：蛋白骨架 RMSD 在 20 ns 后波动进入平台（通常 < 3 Å），才认为平衡区可用。如果 100 ns 仍在单调漂移，说明平衡时间不足或约束设置有问题，后续所有统计量都不可信。

RMSF 定位柔性区域

RMSF（均方根涨落） 按残基给出波动幅度。loop 区与末端天然偏高，但若活性位点附近出现意料外的高峰，往往提示该区域未被正确约束或力场对局部构象描述有偏差。本项目曾遇到一个激酶 activation loop 在 RMSF 上异常活跃，回头发现是初始结构缺失了该 loop 的二面角约束。

氢键与水桥

配体结合稳定性不能只看对接 pose。本项目统计 氢键占有率：占有率 > 30% 的氢键才视为有效相互作用；低于此值的”瞬时氢键”在快照里好看，但对结合贡献有限。对极性口袋，还要看水桥（water bridge）是否持续存在，它常常决定了 MM-PBSA 中极性溶剂项的正负。

溶剂可及面积（SASA）与聚类

配体结合后蛋白 SASA 下降，是结合事件发生的直接证据。本项目通常再对配体构象做聚类（基于 RMSD 或二面角），看 top 簇是否占轨迹主体——如果 top 簇占比 < 50%，说明配体在口袋里反复翻转换姿态，结合并不专一。

自由能分解

MM-PBSA 单轨迹法对每帧做能量分解到残基，得到 per-residue 贡献图。本项目重点看贡献 > 1 kcal/mol 的残基是否与已知突变热点吻合；若不吻合，要警惕极性项被溶剂模型高估。对结合自由能差异 < 2 kcal/mol 的两个配体，结论应谨慎，建议补做热力学积分（TI）或 FEP 交叉验证。

一条 100 ns 轨迹的判读实例

某蛋白-小分子体系：Cα-RMSD 在 15 ns 平台期 2.1 Å；配体位点氢键占有率 42%/38%（双氢键）；top 簇占比 71%；MM-PBSA ΔG = −28.4 kcal/mol，其中三个疏水残基合计贡献 −11.2 kcal/mol。结论：结合稳定且由疏水主导，与湿实验 IC50 趋势一致。

分析报告的边界

需要说明：RMSD/RMSF 描述”结构波动”，不直接等于”功能状态”。柔性蛋白的多个亚稳态可能都需要采样，单条 100 ns 轨迹未必覆盖全貌，本项目对关键结论会建议补长程或副本（replica）采样。分子动力学分析的技术细节与工具链可进一步查阅本站的分子动力学计算栏目。整体服务范围见科研学术网首页。