用途：分子动力学（MD）模拟计算分子体系随时间的演化，通过积分牛顿运动方程实现。该技能覆盖两个互补工具：OpenMM（https://openmm.org/）为支持 GPU 的高性能 MD 模拟引擎，提供 Python API 和灵活力场支持；MDAnalysis（https://mdanalysis.org/）为读取、写入和分析主流模拟包轨迹的 Python 库。
适用场景：蛋白质稳定性分析（如突变对动力学的影响）、药物结合模拟（表征配体结合模式与停留时间）、构象采样、蛋白-蛋白相互作用、RMSD/RMSF 分析、自由能估算、膜模拟、固有无序蛋白研究等。
安装： conda install -c conda-forge openmm mdanalysis nglview 或 pip install openmm mdanalysis
主要能力与用法：
- 体系准备：从 PDB 读取拓扑与坐标，加载力场（如 amber14-all.xml + amber14/tip3pfb.xml），加氢并加溶剂盒（10 Å padding、150 mM NaCl），创建 System（PME 长程静电、1.0 nm 截断、HBonds 约束、刚性水）。
- 能量最小化：用 LangevinMiddleIntegrator 与 Simulation，优先 CUDA、回退 OpenCL/CPU，调用 simulation.minimizeEnergy，输出 minimized.pdb。
- NVT 平衡：setVelocitiesToTemperature，配合 StateDataReporter、DCDReporter 记录。
- NPT 生产：添加 MonteCarloBarostat，重初始化上下文，记录密度、轨迹并保存 CheckpointReporter 检查点。
- 轨迹分析：用 MDAnalysis 加载拓扑与轨迹，进行 RMSD（AlignTraj + rms.RMSD）、RMSF（逐残基柔性）、蛋白-配体接触（4.5 Å 截断）计算与绘图。
力场选择：标准蛋白用 AMBER14 + TIP3P-FB；含小分子用 AMBER14 + GAFF2；膜蛋白用 CHARMM36m + TIP3P；核酸用 AMBER99-bsc1 或 AMBER14；无序蛋白用 ff19SB 或 CHARMM36m。
体系准备工具：PDBFixer 修复缺失残基/原子、去杂原子、补氢（按 pH）；小分子可用 OpenFF Toolkit 生成参数（openff-2.0.0.offxml）。
最佳实践：模拟前必须最小化；平衡顺序为 NVT（50–100 ps）→ NPT（100–500 ps）→ 生产；建议使用 GPU；配合 HBonds 约束使用 2 fs 时间步，HMR 可用 4 fs；分析时舍弃前 20–50% 平衡段；定期保存检查点；溶剂化体系需周期性边界条件；带电荷体系静电用 PME 更准确。
限制：资料未提供具体性能基准、收费与平台支持细节；模拟精度与稳定性取决于力场、体系准备和参数选择，需要具备 MD 基础与计算资源。
参考资料：OpenMM 文档 https://openmm.org/documentation.html；MDAnalysis 用户指南 https://docs.mdanalysis.org/；替代引擎 GROMACS https://manual.gromacs.org/、NAMD https://www.ks.uiuc.edu/Research/namd/；CHARMM-GUI https://charmm-gui.org/；AmberTools https://ambermd.org/AmberTools.php。