手机版

DS 分子动力学模拟

发布时间：2026-08-20 来源：科研学术网

小中大

字号：

这张图是本项目用 Discovery Studio（DS）跑一个约 5 万原子蛋白体系的 MD 结果。左图是三条独立副本（replica）的骨架 RMSD，100 ns 后都收敛到约 2.6 Å 的平台；右图是每残基 RMSF，两个红色高亮区域（约 40–50 残基、105–115 残基）柔性明显高于其他区域，对应蛋白表面的 flexible loop。两图合起来说明：体系整体稳定，但局部 loop 有明显运动，做药物设计时不能把它们当刚性口袋处理。

左图 RMSD 在说什么

三条曲线一开始快速上升，前 20–30 ns 属于平衡调整期；之后进入平台，波动在 2.4–2.7 Å 之间。τ 值（12 ns、18 ns、26 ns）是各自副本的弛豫时间，越短说明平衡越快。三个副本最终都收敛到同一水平，说明采样起点独立但殊途同归，结果可重复。

本项目判稳的量化门槛：

| 判据 | 标准 | 本案例 |

|——|——|——|

| RMSD 平台 | 后 50 ns 漂移 < 0.5 Å | 达标 |

| 副本一致性 | 三个副本 RMSD 均值差 < 0.5 Å | 达标 |

| 能量漂移 | 势能标准差 < 1% 均值 | 达标 |

右图 RMSF 在说什么

RMSF 是残基层面的平均涨落。刚性二级结构（α-螺旋、β-折叠）的 RMSF 通常在 0.5–1.5 Å；loop 区 2–3 Å 很常见。本案例两个 loop region 的 RMSF 超过 2.5 Å，说明它们是构象开关或配体进出口。做分子对接时如果 binding pocket 毗邻这些 loop，必须上诱导契合或 MD 采样，不能 rigid docking。

Discovery Studio MD 设置要点

| 设置项 | 本项目取值 | 说明 |

|——–|——|——|

| 力场 | CHARMm / AMBER ff99SB-ILDN | 蛋白体系用 CHARMm，核酸用 AMBER |

| 溶剂化 | TIP3P 水盒，边缘 ≥ 10 Å | 防止蛋白镜像相互作用 |

| 离子强度 | 0.15 M NaCl | 中和体系并模拟生理盐 |

| 时间步长 | 2 fs | 用 SHAKE/LINCS 约束氢键 |

| 温度 | 300 K，Nose-Hoover | 副本独立初始速度 |

| 输出帧 | 每 10 ps | 共 10,000 帧 |

平衡阶段怎么切

本项目典型流程：

能量最小化 5000 步，去掉原子重叠。
约束蛋白重原子，NVT 升温到 300 K，跑 100 ps。
放开约束，NPT 平衡 1 ns，让水盒密度稳定。
生产跑 100 ns（每个副本），三个副本独立采样。

轨迹分析做什么

除了 RMSD/RMSF，本项目还会做：

– **二级结构演化**：看 α-螺旋/β-折叠是否坍塌。

– **氢键占有率**：配体与关键残基氢键在 MD 中是否 > 50%。

– **回转半径 Rg**：判断整体紧致程度。

– **PCA**：找主要运动模式，定位别构位点。

常见坑

– **单副本陷阱**：一个副本收敛可能是偶然，三个副本一致才算稳。

– **loop 没平衡**：loop 运动慢，100 ns 不够，需加副本或上增强采样。

– **质子化状态错**：Asp/Glu/His 电荷不对，loop 行为全变。

DS 分子动力学模拟是验证对接结果、研究构象动态的核心工具。更多蛋白体系 MD 建模与轨迹分析方法可参考科研学术网 [DS分子动力学模拟](https://www.keyanxueshu.com/category/md/) 栏目，完整流程见首页 https://www.keyanxueshu.com/。