FEP、RBFE 与 ABFE：自由能计算到底在算什么

自由能微扰（FEP）是基于严格统计力学的自由能计算方法。它是目前唯一能在药物设计中提供可靠定量预测的计算方法——但适用条件严格，成本高。理解 RBFE 与 ABFE 的区别，是用好它的前提。

RBFE 与 ABFE

	RBFE（相对结合自由能）	ABFE（绝对结合自由能）
计算的量	ΔΔG（A 与 B 的差）	ΔG（单个分子的绝对值）
适用	同系列类似物	任意分子
成本	较低	高 3～10 倍
可靠性	较高（误差常 < 1.5 kcal/mol）	较低
成熟度	生产可用	仍在发展
典型用途	先导优化中的排序	苗头验证、跨骨架比较

RBFE 是实践中的主力：先导优化阶段的核心问题正是「这几个类似物哪个更强」，而这正是 RBFE 擅长的。

炼金术路径：核心思想

# 问题：直接算结合自由能需要模拟配体进出口袋 ——
#      时间尺度太长，采样不可行
#
# 解决：利用自由能是【状态函数】的性质 ——
#      只与始末状态有关，与路径无关
#      → 可以走一条【物理上不存在但计算上方便】的路径
#
# RBFE 的热力学循环：
#
#          ΔG_bind(A)
#   A + P --------------→ A·P
#     |                    |
#     | ΔG_solv(A→B)       | ΔG_complex(A→B)
#     ↓                    ↓
#   B + P --------------→ B·P
#          ΔG_bind(B)
#
#   ΔΔG = ΔG_bind(B) - ΔG_bind(A)
#       = ΔG_complex(A→B) - ΔG_solv(A→B)
#
# 关键：
#   我们不模拟结合过程（竖线是「炼金术变化」），
#   而是模拟【把分子 A 变成分子 B】这个非物理过程
#   —— 一次在水中，一次在蛋白口袋中
#
# 「炼金术」的含义：
#   通过耦合参数 λ 从 0 到 1，
#   逐渐把 A 的原子「变成」B 的原子
#   （改变原子类型、电荷、消失/出现原子）
#
# 实现：
#   在 λ = 0, 0.1, 0.2, ..., 1.0 各跑一个模拟窗口
#   用 BAR 或 MBAR 方法从各窗口的能量差估计自由能
#
# ABFE 的路径：
#   直接把配体「消失」（λ: 1 → 0）
#   → 需要额外的限制势防止配体在消失过程中漂走
#   → 再校正限制势的贡献
#   → 步骤更多，误差来源更多

成功的必要条件

可靠的结合模式。这是第一位的。FEP 假设所有类似物以相同方式结合。如果结合模式判断错了，计算再精确也没用——最好有共晶结构；
结构差异要小：适合改一个取代基、加一个原子、换一个杂原子。骨架大改超出 RBFE 的可靠范围；
好的力场参数：配体参数化的质量直接决定结果（见 202《OpenFF Toolkit》）。非常规基团、金属配位、极化效应强的体系要格外小心；
充分采样：口袋中的水、侧链构象都需要充分弛豫。柔性大的体系需要更长模拟；
正确的质子化态与互变异构体：这些错了，结果就错了；
净电荷变化要特别处理：如果 A 与 B 的净电荷不同，需要额外的有限尺寸校正，否则误差很大。

准确度的现实

条件	典型 RMSE
理想条件（好结构、小改动、刚性口袋）	约 1 kcal/mol
一般条件	1～2 kcal/mol
困难体系（柔性、带电变化、构象不确定）	> 2 kcal/mol
（对比）对接分数	> 2～3 kcal/mol 且相关性差

1 kcal/mol 约对应 5 倍的活性差异。这意味着 FEP 能可靠区分「10 nM 与 100 nM」，但不能可靠区分「10 nM 与 20 nM」——理解这个分辨率，才能设定正确的期望。

微扰图的设计

# RBFE 计算的是「配对」的自由能差，
# 因此需要设计一张【微扰图】（perturbation map）：
# 节点 = 分子，边 = 要计算的转换
#
# 设计原则：
#
# 1) 每条边的变化要小
#    改动越大，采样越难，误差越大
#
# 2) 【图要有闭合环路】
#    环路上的 ΔΔG 之和应该为 0
#    → 偏离 0 的程度就是【自洽性误差】
#    → 这是【检验计算可靠性的关键手段】
#
# 3) 星形 vs 网状
#    星形（都连到参考分子）：边少，成本低，但无环路检验
#    网状（多连接）：成本高，但有冗余与自洽检验
#    → 【建议至少保证有环路】
#
# 4) 有实验值的分子作为锚点
#    有了锚点，相对值可以转换成绝对值
#
# 5) 避免电荷变化的边
#    如果必须有，尽量安排在图的边缘并单独处理
#
# 工具：
#   Lomap（自动生成微扰图，见 130）
#   konnektor
#
# 环路闭合误差的判读：
#   < 0.5 kcal/mol   良好
#   0.5~1.0          可接受
#   > 1.0            该区域的计算不可信

软件选择

软件	说明
OpenFE	开源，工作流完整（见 130《OpenFE 做 RBFE》）
FEP+（Schrödinger）	商业，最成熟，验证最充分
Amber TI	成熟，需要较多手工配置
GROMACS + PMX	开源组合
BioSimSpace	统一接口，可切换后端
Flare FEP（Cresset）	商业，界面友好

成本与项目中的定位

单条边的成本：约 5～20 GPU 小时（取决于体系大小与 λ 窗口数）；
一个典型的微扰图：20～50 条边 → 数百 GPU 小时 → 单卡需数周，多卡集群数天；
因此 FEP 只用于关键决策：合成前对最终候选的排序，而非大规模筛选；
价值判断：一个化合物的合成 + 测试成本通常远高于 FEP 计算——如果 FEP 能避免合成 5 个不该合成的分子，它就是划算的；
正确的期望：FEP 用于提高合成的命中率，而非取代实验。

关键要点

RBFE 是实践主力，ABFE 成本高数倍且仍在发展；
可靠的结合模式是第一前提——模式错了，计算再精确也没用；
微扰图必须有闭合环路，环路误差是检验可靠性的关键手段；
典型精度约 1 kcal/mol（约 5 倍活性差异）——能区分 10 nM 与 100 nM，不能区分 10 与 20 nM。

延伸资源

OpenFE 实战：130《OpenFE 做 RBFE》；MM/GBSA：126《MM/GBSA》；MD 入门：123《分子动力学 MD 入门》；
力场参数化：202《OpenFF Toolkit》；先导优化：405《Lead Optimization》。