用途
分子动力学（MD）模拟通过积分牛顿运动方程，在计算机上模拟分子体系随时间的演化。该技能覆盖两个互补工具：
- OpenMM（https://openmm.org/）：支持 GPU 的高性能 MD 模拟引擎，提供 Python API 和灵活的力场支持。
- MDAnalysis（https://mdanalysis.org/）：用于读取、写入和分析来自所有主流模拟软件包的 MD 轨迹的 Python 库。
适用场景
- 蛋白稳定性分析：突变如何影响蛋白动力学
- 药物结合模拟：表征配体的结合模式与停留时间
- 构象采样：探索蛋白柔性与构象变化
- 蛋白-蛋白相互作用：建模界面动力学与结合能量学
- RMSD/RMSF 分析：从参考结构量化结构涨落
- 自由能估计：计算结合自由能或构象自由能
- 膜模拟：在脂双层中建模蛋白
- 内在无序蛋白：研究 IDR 构象系综
安装
conda install -c conda-forge openmm mdanalysis nglview
# or
pip install openmm mdanalysis
主要能力与使用方式
1. 体系准备（OpenMM）
从 PDB 文件加载结构，加载力场（如 amber14-all.xml 与 amber14/tip3pfb.xml），通过 Modeller 加氢并添加溶剂盒（10 Å padding、150 mM NaCl），再创建 System，采用 PME 处理长程静电、1.0 nm 非键截断、HBonds 约束、刚性水模型。
2. 能量最小化
使用 LangevinMiddleIntegrator，优先选择 CUDA，其次 OpenCL，最后回退 CPU。检查初始能量后执行 minimizeEnergy，并保存最小化结构。
3. NVT 平衡
设置目标温度速度，添加 StateDataReporter 日志与 DCDReporter 轨迹，运行指定步数（如 50000 步 × 2 fs = 100 ps）。
4. NPT 平衡与生产
添加 MonteCarloBarostat 控制压力，更新报告器（含密度），添加 CheckpointReporter 保存检查点，运行生产模拟（如 500000 步 × 2 fs = 1 ns）。
5. 轨迹分析（MDAnalysis）
- 加载轨迹：Universe(topology_file, trajectory_file)，输出原子数、帧数与时间范围。
- RMSD 分析：AlignTraj 对齐后计算 RMSD，并绘制随时间变化曲线。
- RMSF 分析：按残基计算柔性，返回残基编号与 RMSF 值。
- 蛋白-配体接触：按距离截断（如 4.5 Å）逐帧统计接触残基。
6. 力场选择指南
- 标准蛋白：AMBER14（amber14-all.xml），水模型 TIP3P-FB
- 蛋白 + 小分子：AMBER14 + GAFF2，TIP3P-FB
- 膜蛋白：CHARMM36m，TIP3P
- 核酸：AMBER99-bsc1 或 AMBER14，TIP3P
- 无序蛋白：ff19SB 或 CHARMM36m，TIP3P
7. 体系准备工具
- PDBFixer：修复缺失残基、非标准残基、缺失原子，加氢并标准化，可去除水/配体。
- GAFF2 小分子参数化：通过 OpenFF Toolkit（pip install openff-toolkit）从 SMILES 生成分子、构象与 Interchange 参数。
最佳实践
- MD 前务必最小化：原始 PDB 结构存在空间冲突
- 生产前先平衡：NVT（50–100 ps）→ NPT（100–500 ps）→ 生产
- 使用 GPU：模拟在 GPU（CUDA/OpenCL）上快 10–100 倍
- 2 fs 时间步配合 HBonds 约束为标准做法；使用 HMR（氢质量重分配）时可用 4 fs
- 仅分析平衡后的轨迹：丢弃前 20–50% 作为平衡段
- 保存检查点：MD 运行可能失败，检查点可重启
- 溶剂化体系需周期性边界条件
- 带电体系静电采用 PME，比截断方法更准确
限制
资料未提供该技能的具体收费、价格、平台支持范围、性能基准或效果保证等信息。