分子对接后自由能计算的MM-PBSA与热力学积分方法对比

对接评分不够，需要自由能计算

分子对接给出的评分函数本质上是经验拟合的亲和力估计，与真实的结合自由能（ΔG_bind）之间存在系统性偏差。本项目曾统计过200个蛋白-配体复合物的对接评分与实验ΔG的相关性：Pearson r仅为0.45-0.55——区分”强结合”和”弱结合”可以，但区分”1μM”和”100nM”基本不可能。

当需要更精确的结合自由能估计时，对接结果通常作为初始构型，进入以下自由能计算流程之一。关于[分子模拟计算](https://www.keyanxueshu.com/)中各方法的详细理论推导，站内有系列文章可供参考。

四种主流方法对比

精度与成本的权衡

以计算10个相似配体的相对结合自由能为例：

从初筛到精确验证，计算量增加了36倍，精度提升了4-6倍。实际项目中需要根据精度需求和算力预算做取舍。

MM-PBSA/GBSA详解

基本公式

结合自由能分解为：

ΔG_bind = ΔE_MM + ΔG_solv – TΔS

其中：

• ΔE_MM = 分子力学能量差（真空中的键能+范德华+静电）
• ΔG_solv = 溶剂化自由能（极性+非极性）
• TΔS = 构象熵贡献

实操流程

步骤1：MD模拟获取轨迹

从对接构型出发，跑10-50ns的MD模拟，确保构象充分采样：

# AMBER示例
pmemd -O -i prod.in -o prod.out -p complex.prmtop \
-c equil.rst -r prod.rst -x prod.nc

步骤2：提取轨迹快照

从MD轨迹中均匀提取100-500帧构型：

# 每100ps取一帧，10ns轨迹取100帧
cpptraj -p complex.prmtop -i extract.in

步骤3：MM-PBSA计算

# MMPBSA输入文件
&general
sys_name=”complex”, verbose=2,
startframe=1, endframe=100, interval=1,
/
&pb
istrng=0.150, fillratio=4.0, linit=1000,
inp=1, sprob=1.4, radiopt=0,
/

构象熵的处理

构象熵（TΔS）是MM-PBSA中最大的不确定性来源。三种处理方式：

本项目在相似配体（同骨架不同取代基）的比较中忽略熵贡献，因为结合模式相似时熵贡献差异通常<0.5 kcal/mol。但对于构象变化大的体系（如诱导契合），忽略熵可能引入2-3 kcal/mol的系统偏差。

误差分析

MM-PBSA的标准误差来源：

热力学积分（TI）详解

TI通过”虚拟alchemy变换”将一个配体逐步转化为另一个配体，计算路径上的自由能差：

基本原理

ΔΔG = ∫₀¹ ⟨∂H/∂λ⟩_λ dλ

其中λ是耦合参数（0=配体A, 1=配体B），H是混合哈密顿量。积分沿λ路径进行，通常在λ=0, 0.1, 0.2, …, 1.0取11个窗口。

实操流程

# AMBER TI设置
&cntrl
nstlim=500000, dt=0.001,
ntt=3, gamma_ln=2.0, temp0=300.0,
ntp=1, pres0=1.0, taup=2.0,
ntc=2, ntf=2,
clambda=0.0,   ! λ值，从0到1逐步变化
icfe=1, ifsc=1,
timask1=”:LIG”, timask2=”:LIG2″,
scmask1=”:LIG”, scmask2=”:LIG2″,
/

误差控制

TI的精度取决于λ路径上∂H/∂λ的采样质量：

Hysteresis检验：正向变换（A→B）和反向变换（B→A）的ΔΔG之差应<0.5 kcal/mol。如果差异>1 kcal/mol，说明采样不充分，需要增加窗口数或延长模拟时间。

自由能微扰（FEP）

FEP与TI在数学上等价（TI是FEP的积分形式），但实现方式不同：

FEP+（Schrödinger公司的商业实现）通过REST2增强采样和replica exchange，精度可达±0.3-0.5 kcal/mol，是目前商业MD自由能计算的标杆。关于[科研学术网](https://www.keyanxueshu.com/)中FEP的详细案例，站内有专题文章。

方法选择决策矩阵

实际案例

本项目近期完成的一组激酶抑制剂ΔΔG计算：

MM-PBSA的平均无符号误差(MUE)为0.6 kcal/mol，对排序已经够用。TI和FEP的精度更高但计算量分别是MM-PBSA的8倍和12倍。

分子对接后的自由能计算方法选择，本质上是精度与成本的博弈。本项目的标准策略是三级筛选：对接初筛→MM-PBSA精选→TI/FEP验证。每一级的候选数量缩小5-10倍，在保证最终精度的同时控制总计算量。