大分子分子动力学模拟

本项目对一个大分子体系（约 50,000 原子）做了 100 ns 的显式溶剂全原子分子动力学模拟，通过温度、RMSD 与势能三条曲线确认体系是否真正达到平衡。这张图展示的是平衡阶段的三个关键判据，是判断模拟结果能否用于后续分析的前提。

三张曲线图，怎么看算”稳了”

这三张图是判断模拟”有没有跑稳”的体检报告。第一张温度在 300 K 上下晃，说明热浴正常；第二张 RMSD 先涨后平，说明结构不再乱漂；第三张势能先降后平，说明体系找到了低能态。三张一起看，约 400 ps 之后全进平台，这时候再拿轨迹做统计才靠谱。只要有一张还在飘，都不能算平衡。

计算逻辑

1. 体系构建：大分子 + 溶剂盒 + 离子中和，总原子数约 50,000。
2. 能量最小化：消除初始结构冲突。
3. 升温与平衡：NVT 下从 0 K 升温至 300 K，再切 NPT 调整密度。
4. 生产模拟：NVT 或 NPT 下跑 100 ns，步长 2 fs，氢键约束用 SHAKE/LINCS。
5. 平衡判据：温度、RMSD、势能三条曲线同时稳定。

关键参数

体系构建清单（可直接照做）

本项目从 PDB 到可跑模拟的完整准备链路：

6. 结构清洗：用 `pdb4amber` 去掉水分子/配体，补缺失残基，修正非标准残基。
7. 质子化状态：按目标 pH（通常 0–7.4）用 `H++` 或 `propKa` 确定组氨酸/天冬氨酸/谷氨酸电荷。
8. 溶剂化：`tleap` 加 TIP3P 水盒，边缘离蛋白 ≥ 10 Å；加 Na⁺/Cl⁻ 至 15 M 中性盐。
9. 力场装配：`leaprc.protein.ff19SB` + `leaprc.water.tip3p` + `leaprc.gaff2`（配体）。
10. 能量最小化：先约束蛋白跑 2000 步，再放开跑 5000 步，消除原子重叠。

这一步没做扎实，后面 100 ns 跑得再久也是白费。

结果判读

• 温度曲线 (a)：前 100 ps 波动较大，之后围绕 300 K 小幅震荡，说明热浴已稳定。
• RMSD 曲线 (b)：从 0 快速上升到约 5 Å 后进入平台，表示结构已弛豫到平衡构象族。
• 势能曲线 (c)：前 200 ps 快速下降，随后进入噪声平台，说明体系已找到低能构象。

可借鉴经验

• 大分子模拟不要只看 RMSD，势能不降、温度不稳都说明没平衡。
• 100 ns 对大蛋白构象变化可能只是”快照”，如果关心慢过程，要补副本交换或元动力学。
• 溶剂化层对大分子稳定性影响很大，构建体系时水盒边缘至少离蛋白 10 Å。

平衡判据量化门槛

本项目对”平衡”不靠肉眼，设了三条可量化门槛，三条同时满足才切生产轨迹做统计：

补充：若关心配体结合或别构这样的慢过程，100 ns 往往不够，本项目会先跑 3–5 个独立副本（replica）取并集，或上副本交换 MD。

局限

经典 MD 无法描述成键/断键；力场参数缺失的配体或金属中心需用 QM/MM 或专用参数。更多分子动力学案例可参考本站的分子动力学计算栏目。整体服务见科研学术网首页。