1、自由能估计基础：什么是自由能，为什么要估计自由能，热力学原理简介

咱们先聊点最根本的东西。搞药物筛选的人，嘴上天天挂着“结合自由能”，但自由能到底是什么？说白了，它就是个系统的“身份证”。系统在某个状态下的稳定性，全看它的自由能高低。自由能越低，状态越稳定，就像一块石头滚到谷底，你再想推它上去就得费力气。

我刚开始做药物筛选的时候，总对这东西感觉抽象。直到有一次，我拿一个看起来不错的分子去做活性测试，结果完全没活性。后来一算，才发现它的结合自由能竟然比参考分子还高。嗯，从那以后，我牢牢记住了——自由能不是纸上谈兵，它直接决定了分子的死活。

1.1 自由能是什么——做药物筛选的人必须懂的概念

自由能（Free Energy）这名字听起来高大上，其实可以拆两层理解：

• 吉布斯自由能（G）——最常用，对应恒定温度和压力下的过程。比如在人体内（37°C、1 atm）发生的药物结合反应，就看它。
• 亥姆霍兹自由能（A）——对应恒定温度和体积，在溶液体系里用得少一些，但搞清楚概念也有用。

为什么“自由”这两个字特别重要？你想想看，系统里其实大部分能量都锁在热运动里，不能用来做功。只有自由能这笔钱，才是随时可以花的“现钞”。一个药物分子能不能从水中挤出身来、钻进靶蛋白的口袋里，就看它在不同环境下的自由能差额。

1.2 为什么非要估计自由能？——关于结合亲和力

我自己入行时经常问：给定10万个候选化合物，怎么挑出那1-2个值得做实验的？对接分数显然不够可靠。真正的“铁指标”是结合自由能（ΔG_bind）。

ΔG_bind 越负，结合越强。一般规则如下：

• ΔG < -10 kcal/mol → 很强的结合（IC50 通常在 nM 级别）
• -8 kcal/mol < ΔG < -6 kcal/mol → 中等结合
• ΔG > -5 kcal/mol → 很弱，不太值得做候选

我刚开始做药物筛选的时候，觉得直接看对接分数就行，结果踩了坑。后来我手头的一个项目，用了 FEP 重新评估7个分子，3个被筛掉（因为它们结合自由能太高），剩下4个里2个在实验里果然有 nM 级活性。

1.3 热力学基础：自由能从哪来？——用统计力学的语言说清楚

自由能不是凭空冒出来的。它来自于统计力学，核心公式长这样：

G = -k_B T ln Q

其中：
G  = 吉布斯自由能
k_B = 玻尔兹曼常数
T  = 绝对温度
Q  = 配分函数（系统所有可能状态的能量加权和）

好，公式归公式，实际算的时候我们没法直接算 Q——因为状态空间太大了，接近无限。那我们就用采样。说白了，把所有可能的构象都跑一遍，然后统计平均。这就是为什么分子动力学模拟（MD）是自由能估计的标配。

我记得第一次跑自由能微扰的时候，温度指定错了，结果跑了三天才发现系统热化有问题——浪费时间啊。后来就养成了习惯：模拟之前先跑5次不同随机种子的平衡，确认体系稳定了，再开正式的ΔG计算。

这里有件很多人不知道的事：ΔG 不是单独通过MD直接测的，而是通过系综平均、逐步微扰或热力学积分算出来的。它是一个统计物理量，不是瞬时的东西。所以采样不够，结果就漂移得离谱。

1.4 自由能估计方法概览——不同场景用不同的武器

方法很多，但核心就几类。我个人习惯把它分成三组：

• 方法A —— 热力学积分（TI）/ 自由能微扰（FEP）：精度高，适合成熟前期到先导优化阶段。算力贵、耗时长。
• 方法B —— MM-PBSA/GBSA：精度中等，速度快，适合批量过滤。我当时拿它先筛100个分子，选出10个再扔给FEP强化。
• 方法C —— 能量打分函数（AutoDock Vina 等）：快如闪电，但从严格意义上来说这不是自由能估计。用来当第一道筛子还行。

1.5 本章小结框架 —— 一张图看懂

下面这张图把第一章的知识体系串起来了。你看一下大概就有了整体把握：