自由能基础知识：热力学基本概念，吉布斯自由能
做药物分子优化，你要是连自由能是什么都说不清楚，那后面的FEP、TI之类的方法基本就是瞎搞。这一章，咱们先把热力学最核心的那几个概念掰扯明白。
我记得刚入行那会儿，带我的老工程师跟我说了一句话，至今都觉得特别在理——「自由能就是分子世界的货币」。你想想看，药物和靶点结合，本质上就是一场「交易」，交易的方向和规模，全看自由能怎么分配。明白了这个，你就能理解为什么我们整天算ΔG算得那么起劲。
1. 系统、环境与状态函数
热力学里首先得分清「系统」和「环境」。系统就是咱们研究的对象，比如一个蛋白质-配体复合物。环境就是系统外面的一切，比如水、离子、膜。做自由能计算时，我最常用的模型就是把蛋白质、配体和一部分水分子包在盒子里，剩下的用周期性边界来处理。
系统分类（按与环境交换情况）
- 孤立系统：不交换物质，也不交换能量。理想模型，实际不存在。
- 封闭系统：不交换物质，但交换能量。实验室最常见，比如密封试管。
- 开放系统：既交换物质也交换能量。活细胞、开放反应池都属于此类。
状态函数这玩意儿你得理解透。它只取决于系统的当前状态，跟怎么走到这个状态没关系。温度、压力、体积、内能、熵、焓、自由能，全算状态函数。我做过一个项目，需要计算配体从水相转移到蛋白结合口袋的自由能差，中间试了三条不同的热力学路径，结果最终ΔG一模一样——这不就印证了状态函数的性质嘛。
2. 焓与熵——自由能的两根支柱
吉布斯自由能公式你们肯定见过：G = H - TS。这里的H是焓，S是熵。我们来拆开讲。
焓（H）：说白了就是系统内部能量加上对外做的体积功。恒压条件下，焓变ΔH就等于吸收或放出的热量。药物结合时，ΔH < 0通常意味着形成了更稳定的弱键（氢键、范德华力）。我见过不少优化案例，配体加个甲基，ΔH变好了，结果ΔS崩了，最后总自由能没捞着好。
熵（S）：系统的混乱度，也可以理解成「可能的微观状态数」的对数。公式是S = kB ln Ω。结合过程往往伴随「熵惩罚」——配体从自由游荡变成被固定，平动熵、转动熵、部分构象熵都没了。我曾经接手过一个先导化合物，结合ΔG只有-6 kcal/mol，操作了半天焓就是提不上去。后来调整了一个柔性链，释放了配体一部分构象熵，ΔG硬生生拉到了-9.5。所以别只盯着焓，熵有时候是更关键的变量。
| 概念 | 物理意义 | 结合时的常见变化 |
|---|---|---|
| 焓变 ΔH | 化学键能量的变化 | 通常 < 0（释放能量） |
| 熵变 ΔS | 体系混乱度的变化 | 通常 < 0（熵惩罚） |
| 自由能变 ΔG | ΔH - TΔS，驱动力的判据 | ΔG < 0 自发结合 |
避坑指南：我曾经在处理一个强疏水体系时，只看ΔH很好（-12 kcal/mol），忽略了ΔS的巨大惩罚（-TΔS 贡献+8 kcal/mol），以为结合很强。结果实验测下来Kd只有μM级别。后来才意识到水分子释放带来的熵增益被我低估了。所以永远要同时看焓和熵两个部分，不要偏科。
3. 吉布斯自由能——平衡的判据
恒温恒压下，吉布斯自由能的减少方向就是自发过程的方向。数学判据很简单：
ΔG < 0：自发过程（体系趋于平衡）
ΔG = 0：平衡状态（净变化为零）
ΔG > 0：非自发（需要外界做功）
在药物优化里，我们最关心的就是配体-蛋白结合过程的ΔGbind。结合自由能 ΔGbind = Gcomplex - (Gprotein + Gligand)。这个值越负，结合越强。一般优化目标就是让ΔG＜-10 kcal/mol（对应nM级别的亲和力）。
一个快速估算方法：室温下（300K），ΔG每变化1.36 kcal/mol，平衡解离常数Kd大约变化10倍。所以你把ΔG优化「出溜」下去3 kcal/mol，亲和力就能提升近1000倍。我个人习惯先算相对自由能（RBFE），快筛几十个类似物，锁定方向再跑绝对自由能。
常用换算关系：
ΔG = -RT ln(Kd)
R = 1.987 cal/(mol·K) ≈ 0.001987 kcal/(mol·K)
在T = 300 K时：
Kd (M) = exp(ΔG / (0.001987 * 300))
反过来：
ΔG (kcal/mol) = -0.592 * log10(Kd) (注意这是常用电化学近似，真正的热力学公式要用ln)4. 从热力学基础到计算实践
学到这里你可能会问：那我们在计算机里怎么算自由能？说白了，分子模拟不可能采样到所有微观状态，所以我们只能靠「热力学积分」或者「自由能微扰」这类手段，把两个状态之间的自由能差算出来。底下这张图概括了热力学基础概念到实际计算方法的核心逻辑：
说回基础——没有扎实的热力学理解，你连什么时候该用TI、什么时候该用FEP都判断不了。焓和熵的拆分尤其重要。我建议你在做任何自由能计算之前，先问自己三个问题：
- 这个结合过程的焓变主要来自哪些相互作用？（氢键？静电？范德华？去溶剂化？）
- 熵惩罚主要来自哪个部分？（配体的构象？蛋白口袋的柔性？水的位移？）
- 温度和浓度对结果有什么影响？（很多实验室条件是室温25°C，但如果你做的是恒温动物体内模拟，温度就得跟着升到37°C）
本章要点速记
- 自由能 = 焓 - 温度 × 熵（G = H - TS）
- ΔG < 0 是自发结合的必要条件
- 药物优化追求更负的ΔGbind（-10 kcal/mol以下为佳）
- 焓和熵需要分开分析，不能只看总和
- 自由能计算的核心是「状态函数」性质——路径无关
这一章讲的都是最底层的热力学概念。你把它吃透了，后面讲自由能微扰、热力学积分的时候才不会觉得那些数学公式是一堆纸老虎。嗯，概念就是这样，越基础的东西，越值得花时间慢慢品。