大分子分子动力学模拟：三链蛋白体系骨架涨落与溶剂可及表面积演化

本项目对一条 158 残基蛋白进行 1 μs 全原子分子动力学模拟，力场 AMBER ff14SB，水模型 TIP3P，三个独立 500 ns 副本。分析骨架 RMSD、逐残基 RMSF、回旋半径 Rg 与溶剂可及表面积（SASA）随时间演化，评估体系收敛性。

案例分析

三个副本的骨架 RMSD 在 100–200 ns 后进入稳定平台，平均约 2.0 Å，说明体系已达到平衡。两个柔性 loop（loop1、loop2）区域的 RMSF 超过 1.02 Å 阈值，是结合口袋动态变化的主要来源；核心 β-barrel 区域 RMSF 低于 0.7 Å，结构稳定。回旋半径 Rg 稳定在 21.4 Å 左右，SASA 约 13800–14500 Ų，压实约 0.6 Å，埋藏面积 640 Ų，二级结构含量 41%。

模型细节、验证与工程意义

模拟在 AMBER 22 中运行，力场 ff14SB + TIP3P，SHAKE 约束氢键，PME 处理静电，Langevin 控温。三个独立副本从不同初始速度出发，每体系 500 ns。RMSF 用 cpptraj 计算，二级结构含量用 DSSP。验证时核心 β-barrel RMSF 与 B-factor 折算值趋势一致；Rg 与 SAXS 估算值偏差 <5%。

总结教训

第一个教训是 1 μs 模拟对中小型蛋白通常足够收敛，但柔性 loop 可能仍未完全采样，需结合副本交换或增强采样。第二个教训是 RMSF 阈值 1.0 Å 可作为 loop 柔性的经验界限，超过此值提示药物设计时需考虑诱导契合。第三个教训是三个独立副本一致性比单条长轨迹更能判断收敛，副本间 RMSD 差异应 <0.5 Å。

项目经验

1. 蛋白全原子 MD：AMBER ff14SB + TIP3P。
2. 总模拟时长 1 μs，三个独立 500 ns 副本。
3. 骨架 RMSD 约 2.0 Å，100–200 ns 后达稳态。
4. 两个柔性 loop 区域 RMSF >1.02 Å。
5. 核心 β-barrel 区域 RMSF <0.7 Å。
6. 回旋半径 Rg ≈ 21.4 Å。
7. SASA 约 13800–14500 Ų。
8. 体系压实约 0.6 Å，埋藏面积 640 Ų。
9. 二级结构含量 41%。
10. 1 μs 对核心区域足够收敛，loop 可能需增强采样。
11. RMSF >1.0 Å 提示 loop 柔性，药物设计需考虑诱导契合。
12. SHAKE 约束氢键，积分步长 2 fs。
13. PME 处理长程静电，cutoff 10 Å。
14. 本项目方法可迁移到蛋白-配体结合、变构机制研究。

补充经验

1. 增强采样（metadynamics、REST）可加速 loop 构象采样。
2. MM-GBSA 结合自由能计算建议取 100 ns 以上平衡段。
3. loop 区域若 RMSF 持续上升，可能预示部分去折叠，需延长采样。
4. pH replica exchange 可研究质子化态对柔性的影响。

复现参数速查（案例图硬参数）

更多相关仿真经验，可参考本站的 MD 仿真案例栏目。