分子对接结合能

分子对接结合能，是每一个做完对接、想把这个”结合有多强”从排名数字变成可解释能量的人，必须跨的一道坎。我刚学对接那会儿，兴冲冲拿对接分当结合能写进报告，说”这个分子结合能 -9 kcal/mol”，结果被合作者一句话问住：”你这单位哪来的？对接分是无量纲的。”后来才懂，分子对接结合能难的从来不是读个数，而是你分不分的清”对接分”和”结合自由能”根本不是一回事。

一、背景：为什么对接分不能当结合能分子对接结合能

对接打分是为了”区分好坏姿态”设计经验函数，它和真实结合自由能（实验上 ΔG = RT ln K_d）差着一个数量级，还无量纲。做苗头优化的人要算相对结合能排取代基优先级，做机理的人要解释为什么某个突变削弱结合，都需要真正的能量，而不是排名。可以说，分子对接结合能是你把”排第几”升级成”强多少”的关键。很多新手把对接分当 kcal/mol 报，等到别人拿实验 K_d 一对，差出几倍才慌。更现实的是，只排名不定量，你没法指导结构优化——不知道哪个相互作用贡献最大，改起来就是瞎蒙。

二、核心原理：从经验分到自由能

要把结合量化，得从对接分走到结合自由能 ΔG_bind。MM-GBSA 是最常用的端点到端点方法：在对接 pose 基础上，对复合物、受体、配体分别做短时 MD 采样，用分子力学（键/角/二面+静电+范德华）减去广义玻恩（GB）溶剂化能，得到 ΔG = <E_complex − E_receptor − E_ligand> + 溶剂化项。它比对接分准得多，能给出 kcal/mol 量级的相对值。更准的是 FEP（自由能微扰）和 TI（热力学积分），用炼金术路径逐点算 ΔG，误差可到 1 kcal/mol 内，但只适合小改（如逐个取代基）。关键认知是：对接分用于排序，MM-GBSA 用于相对定量，FEP 用于高精度取代基优化，三者精度递增、代价也递增。还有一个点：熵，MM-GBSA 常用 Normal Mode 估算构象熵，但噪声大，相对能比较时常忽略熵只比焓。

三、关键技术要点：轨迹、熵与分解

MM-GBSA 的质量全在采样：pose 要先做约束 MD 松弛再采样，frames 取 50–200 帧取平均，太少噪声大。隐性溶剂的 GB 模型（如 GBneck2）要和力场匹配，否则静电项失真。熵项（Normal Mode 或准谐波）对大体系噪声大，我一般只比较焓差 ΔH 用于相对排序，除非要 absolute ΔG。关键的增值项是每残基能量分解（per-residue decomposition），能直接告诉你哪个残基贡献最大——这是指导结构优化最实用的输出。我早年优化一个激酶苗头，对接分看不出门道，MM-GBSA 分解一出，发现贡献全在一个水介导的氢键上，改成直接氢键后活性涨了 5 倍。另一个常被忽略的点：ΔG 的误差棒，我交付一定给均值±标准差，单点 MM-GBSA 数值不可信。

四、实操流程：从 pose 到能量结论分子对接结合能

第一步，取对接 top pose，做约束 MD 松弛消除初始应变；第二步，生产采样 MD（ns 级），取若干帧；第三步，对每帧算 MM-GBSA 能量，求平均和标准差；第四步，做 per-residue 分解，找关键贡献残基；第五步，若要做取代基优化，挑关键位点用 FEP 逐点算 ΔΔG。我习惯用 MM-GBSA 先排相对强弱、定位关键残基，再对少数候选上 FEP 精修。对蛋白-核酸体系还要确认离子在能量里的处理一致。对结合模式不确定的 pose，我会先 MD 验证稳定再做能量，避免给假 pose 算出一个漂亮但假的 ΔG。

我把对接分当 kcal/mol 报的错，早年犯过一次就再没犯。对接分无量纲只用于排序，真要能量走 MM-GBSA，要精度上 FEP。凡是指导结构优化的能量，必须过 MM-GBSA 甚至 FEP，对接分只用于初筛不用于排名结论。把排名和定量混为一谈，改结构就是瞎蒙。

per-residue 能量分解是我最常用的指导工具。单看总 ΔG 不知道改哪，分解一出哪个残基贡献最大、加什么相互作用立刻清晰。这是对接分永远给不了的指导价值，也是结合能计算最实用的产出。我交付结合能，分解表往往比总能量更被合作者看重，因为它直接指向下一步怎么改分子。

我交付任何 MM-GBSA 结果都附均值±标准差。单点能量不可信，采样帧少噪声大，没有误差棒的数值容易被当成精确结论。相对能比较时我常忽略熵项只比焓，但会注明；要绝对 ΔG 才上 FEP。把不确定度量出来，比给一个漂亮但虚的数负责任得多。

五、常见问题

Q1：对接分能当 kcal/mol 吗？不能，它无量纲只用于排序。 Q2：MM-GBSA 噪声大？采样帧太少，加到 100–200 帧取平均。 Q3：相对能差不对？pose 没松弛或 GB 模型错，重做采样。 Q4：哪个残基关键？做 per-residue 分解，直接定位贡献最大者。 Q5：要绝对结合能？上 FEP/TI，但只适合小改且机时高。 Q6：熵项忽高忽低？Normal Mode 噪声大，相对比较可省熵。 Q7：单点数值能报吗？不能，必须给均值±标准差。

讲一个具体的坑：有次给一个系列苗头算结合能，我直接拿对接分排序就让合作者改结构，改了三轮活性纹丝不动。回头用 MM-GBSA 一算，发现排名第二的分子实际 ΔG 比第一还低 1.5 kcal/mol，排序被对接分误导了。从那以后我固定一条规矩：凡是指导结构优化的能量，必须走 MM-GBSA 甚至 FEP，对接分只用于初筛不用于排序结论。另一个常被忽略的点：能量分解，单看总 ΔG 不知道改哪，per-residue 一出来改哪个残基、加什么相互作用立刻清晰，这是对接分永远给不了的指导价值。

六、复盘总结

回过头看，分子对接结合能最容易被当成”读个对接分”的偷懒活，但它真正的价值在把排名升级成可解释的能量：对接分排序、MM-GBSA 定量、FEP 精修、分解定位。我现在的习惯是用对接分初筛、MM-GBSA 排相对强弱并分解关键残基、关键候选上 FEP，四步少一步都不敢拿能量去指导结构优化。结合能不是对接分换个单位，它是从排序到机理的一次跃迁。把采样、分解、误差三件事钉死，算出来的每个 kcal/mol 才敢拿去改分子。