本项目对一个大分子体系（约 50,000 原子）做了 100 ns 的显式溶剂全原子分子动力学模拟，通过温度、RMSD 与势能三条曲线确认体系是否真正达到平衡。这张图展示的是平衡阶段的三个关键判据，是判断模拟结果能否用于后续分析的前提。
三张曲线图，怎么看算”稳了”
这三张图是判断模拟”有没有跑稳”的体检报告。第一张温度在 300 K 上下晃，说明热浴正常；第二张 RMSD 先涨后平，说明结构不再乱漂；第三张势能先降后平，说明体系找到了低能态。三张一起看，约 400 ps 之后全进平台，这时候再拿轨迹做统计才靠谱。只要有一张还在飘，都不能算平衡。
计算逻辑
- 体系构建：大分子 + 溶剂盒 + 离子中和，总原子数约 50,000。
- 能量最小化：消除初始结构冲突。
- 升温与平衡：NVT 下从 0 K 升温至 300 K，再切 NPT 调整密度。
- 生产模拟：NVT 或 NPT 下跑 100 ns，步长 2 fs，氢键约束用 SHAKE/LINCS。
- 平衡判据：温度、RMSD、势能三条曲线同时稳定。
关键参数
| 参数 |
取值 |
说明 |
| 力场 |
AMBER ff19SB + TIP3P |
蛋白-水体系经典组合 |
| 温度 |
300 K |
Nose-Hoover 热浴 |
| 时间步长 |
2 fs |
含 H 约束 |
| 模拟时长 |
100 ns |
大分子构象采样 |
| 输出帧 |
每 10 ps 一帧 |
共 10,000 帧用于分析 |
体系构建清单（可直接照做）
本项目从 PDB 到可跑模拟的完整准备链路：
- 结构清洗：用 `pdb4amber` 去掉水分子/配体，补缺失残基，修正非标准残基。
- 质子化状态：按目标 pH（通常 0–7.4）用 `H++` 或 `propKa` 确定组氨酸/天冬氨酸/谷氨酸电荷。
- 溶剂化：`tleap` 加 TIP3P 水盒，边缘离蛋白 ≥ 10 Å；加 Na⁺/Cl⁻ 至 15 M 中性盐。
- 力场装配：`leaprc.protein.ff19SB` + `leaprc.water.tip3p` + `leaprc.gaff2`（配体）。
- 能量最小化：先约束蛋白跑 2000 步，再放开跑 5000 步，消除原子重叠。
这一步没做扎实，后面 100 ns 跑得再久也是白费。
结果判读
- 温度曲线 (a)：前 100 ps 波动较大，之后围绕 300 K 小幅震荡，说明热浴已稳定。
- RMSD 曲线 (b)：从 0 快速上升到约 5 Å 后进入平台，表示结构已弛豫到平衡构象族。
- 势能曲线 (c)：前 200 ps 快速下降，随后进入噪声平台，说明体系已找到低能构象。
可借鉴经验
- 大分子模拟不要只看 RMSD，势能不降、温度不稳都说明没平衡。
- 100 ns 对大蛋白构象变化可能只是”快照”，如果关心慢过程，要补副本交换或元动力学。
- 溶剂化层对大分子稳定性影响很大，构建体系时水盒边缘至少离蛋白 10 Å。
平衡判据量化门槛
本项目对”平衡”不靠肉眼，设了三条可量化门槛，三条同时满足才切生产轨迹做统计：
| 判据 |
稳定标准 |
本项目结果 |
| 温度波动 |
±2% 目标值且无明显漂移 |
300 K ± 4 K |
| RMSD 平台 |
最后 20 ns 漂移 < 0.5 Å |
平台区 < 0.3 Å |
| 势能波动 |
最后 20 ns 标准差 < 1% 均值 |
0.6% |
补充：若关心配体结合或别构这样的慢过程，100 ns 往往不够，本项目会先跑 3–5 个独立副本（replica）取并集，或上副本交换 MD。
局限
经典 MD 无法描述成键/断键；力场参数缺失的配体或金属中心需用 QM/MM 或专用参数。更多分子动力学案例可参考本站的分子动力学计算栏目。整体服务见科研学术网首页。
