1. FEP 计算概述：到底什么是 FEP？

大家好，欢迎来到这门课。我是你们这趟旅程的向导，一个在计算化学和药物设计里摸爬滚打了好些年的老工程师。

咱们开门见山。FEP，全称是 Free Energy Perturbation，自由能微扰。听起来有点玄乎对吧？

其实说白了，它就是一套计算化学方法。用来干嘛呢？用来算两个分子在结合时，自由能到底差多少。

你想想看，药物分子和靶点蛋白结合，就像钥匙插进锁孔。结合得好不好，顺不顺，关键就看这个“自由能变化”。FEP 就是那个能提前告诉你“这把钥匙行不行”的计算工具。

1.1 物理图像：从“态”到“态”的能量差

FEP 的数学基础其实就一个核心公式：

ΔG(B←A) = -kBT · ln〈exp[-(HB - HA)/kBT]〉A

嗯，看起来有点吓人。但别怕，咱们不用背公式。

它的物理含义很简单：我们从状态 A（比如老的配体）慢慢“微扰”到状态 B（新的配体），去算这中间自由能的变化。

为什么会用“微扰”这个词？因为 FEP 不会让你一步跨过去。它会在 A 和 B 之间，插入好多好多“中间态”。每个中间态都是一个混合体，既有 A 的影子，又有 B 的样子。

我记得自己刚开始接触 FEP 时，总觉得这东西有点“玄学”。直到我亲手跑了一个最简单的甲烷变乙烷的体系，看到算出来的自由能差和实验值几乎一模一样——那种感觉，嗯，就像第一次看到实验数据跟自己算的吻合一样，挺激动的。

1.2 FEP 在药物设计中的角色

谈到药物设计，大家更熟悉的可能是对接（Docking）、分子动力学模拟（MD）。那 FEP 在这中间扮演什么角色呢？

我个人习惯把计算工具按“精度”和“速度”排个序：

• 对接（Docking）：几秒钟筛几十万分子，但能量精度差，经常“假阳性”。
• 分子力学（MM/GBSA）：几小时算上百分子，精度比对接好一点，但只算“终态”能量。
• FEP：几小时算一对分子。精度能到 1 kcal/mol 以内。这才是“金牌标准”。

说白了，FEP 不是用来做“大海捞针”的虚拟筛选的。它是用来做“最后一公里”的精准优化的。

我在项目里遇到过不止一次：通过对接找到了一堆活性化合物，但怎么改都提不高活性。用 FEP 一算，才发现原来我们一直在朝着错误的方向做替换。嗯，那之后，我们就把 FEP 加进了核心流程里。

1.3 它在药物设计中的五大应用场景

我把 FEP 最常见的用途整理成了一个表格。这个表你在项目里随时可以拿来参考：

应用场景	具体干什么	价值
先导化合物优化	比较一系列类似物的结合自由能	快速筛选出活性最好的骨架
官能团替换评估	甲基变乙基、苯环变吡啶等微小改动	避免合成“越改越差”的化合物
选择性预测	计算化合物对野生型 vs 突变型的亲和力差异	解决耐药性问题
手性/构象偏好	判断哪个手性中心活性和选择性更好	指导不对称合成方向
代谢稳定性预估	间接预测哪个位置容易被代谢酶攻击	避免合成后才发现代谢太快

你可能会问：FEP 这么准，那以后是不是不用做实验了？

我的答案是：千万别这么想。FEP 是给实验指方向的，不是替代实验的。我在实际项目中踩过一个坑：某个体系 FEP 跑得特别漂亮，结果合成出来后活性降低了 10 倍。后来查了好久才发现，是我们算的那套力场参数，根本不准。

1.4 一张图看懂 FEP 的核心逻辑

说了这么多，我们来画一张图总结一下。这张图展示了从“物理原理”到“实际应用”的完整脉络。

1.5 版本控制？为什么 FEP 需要它？

嗯，说到这，你可能会想：FEP 听起来就是个软件工具，为什么还要搞版本控制？

我告诉你一个真实的教训吧。之前我管理过一个 FEP 项目，组里三个人同时算同一个体系的类似物。结果一个人用的力场是 GAFF 2.1，另一个人用的正好是 2.11——两个版本算出来的 ΔΔG 差了 0.4 kcal/mol。

你说这两个结果哪个对？没人说得清。因为没有版本记录。

这就是为什么我要把“版本控制”和“FEP 计算”放在一起讲。在后续的章节里，我们会一步步搭建一个可重现、可追溯的 FEP 计算管线。