DP Technology RiDYMO：分子动力学平台案例

RiDYMO 是深势科技（DP Technology）的分子动力学平台，特色在增强采样——用改进的采样方法探索常规 MD 到不了的构象空间，用于发现隐蔽口袋（cryptic pocket）和研究难成药靶点。

常规 MD 的采样困境

• 生物学相关的构象变化（口袋开合、别构转换、配体解离）常发生在微秒到毫秒时间尺度；
• 常规 MD 单条轨迹通常只能跑到纳秒到微秒；
• 体系被困在能量势阱中，跨越能垒的事件在可及时间内根本不发生；
• 结果是：你跑了 100 ns，看到蛋白在原地抖动，什么也没发生，然后误以为「这个口袋不会打开」。

增强采样的主要方法

方法	思路	特点
Metadynamics	在已访问区域堆积偏置势，「填平」势阱	需预定义集体变量（CV）
Umbrella Sampling	沿预设路径加约束分窗采样	需知道路径，可算自由能剖面
Replica Exchange（REMD）	多温度副本交换，高温帮助翻越能垒	不需 CV，但计算量大
Accelerated MD / GaMD	整体降低能垒	不需 CV，易上手
Adaptive Sampling	迭代地从欠采样区域重启轨迹	适合分布式算力
RiD（Reinforced Dynamics）	用机器学习自适应地构造偏置	减少对人工选 CV 的依赖

「集体变量怎么选」是增强采样最大的实践难点：Metadynamics 需要你事先指定「沿哪个方向加速探索」，选错了就白算。这正是机器学习辅助方法（如 RiD 类思路）试图解决的问题。

隐蔽口袋：为什么值得找

• 很多靶点在晶体结构中「看起来无处下手」——表面平坦、没有明显凹陷；
• 但蛋白是动态的，某些口袋只在特定构象下短暂出现，晶体结构捕捉不到；
• 找到这类隐蔽口袋，可能把「不可成药」的靶点变成可成药——KRAS G12C 的 switch II 口袋是最著名的例子，它在常规结构中并不明显，却成为多个已上市药物的结合位点。
• 别构位点也常属于这一类。

开源实现路径

# OpenMM + PLUMED 是最主流的开源组合
mamba install -c conda-forge openmm plumed py-plumed openmm-plumed

# plumed.dat —— Metadynamics 示例
# 用两个残基间距离作为集体变量（探索口袋开合）
d1: DISTANCE ATOMS=1234,5678

METAD ...
  ARG=d1
  PACE=500 HEIGHT=1.2 SIGMA=0.05
  BIASFACTOR=15          # well-tempered，保证收敛
  FILE=HILLS
... METAD

PRINT ARG=d1,metad.bias FILE=COLVAR STRIDE=100

from openmm.app import *
from openmm import *
from openmmplumed import PlumedForce

system.addForce(PlumedForce(open("plumed.dat").read()))
# 其余与常规 OpenMM 流程相同（见 200）

其他开源选项：GaMD（OpenMM/AMBER 内置，不需选 CV，最易上手）、PySAGES、SSAGES；口袋检测用 fpocket 逐帧扫描轨迹，统计哪些口袋在多少比例的帧中出现。

# 用 fpocket 扫轨迹找隐蔽口袋的思路
for i in $(seq 0 100 10000); do
  # 从轨迹抽帧
  python extract_frame.py --traj traj.dcd --top top.pdb --frame $i --out f$i.pdb
  fpocket -f f$i.pdb
done
# 统计各口袋出现频率与体积分布，
# 在晶体结构中不存在、但在轨迹中反复出现的口袋 = 隐蔽口袋候选

使用注意

• 增强采样不创造真相：它加速了采样，但如果力场不准或体系设置有误，只会更快地得到错误结果。
• 收敛性必须检验：Metadynamics 要看自由能面是否收敛（偏置势增长是否趋缓），不收敛的结果没有意义。
• 找到的口袋要验证可成药性：能打开不等于能结合分子。需要评估口袋体积、疏水性、是否有足够的相互作用位点。
• 算力投入大：增强采样的成本远高于常规 MD，要有预算准备。
• 最终仍需实验验证：隐蔽口袋的存在需要通过片段筛选、共晶结构等实验手段确证。

延伸资源

• 模拟引擎：200《OpenMM》；结构动力学：209《ProDy》；
• 公司：316《DP Technology 深势科技》；口袋与别构概念见「结构与模拟」模块。