📚 课程目录第2章结合自由能的热力学与统计力学基础

2. 结合自由能的热力学与统计力学基础

在蛋白–配体结合自由能计算中，我们真正想回答的问题并不是“分子之间有没有相互作用”，而是“在给定温度、浓度和标准态下，配体从自由溶液进入蛋白结合位点，体系自由能改变了多少”。这个量就是结合自由能 ΔG_bind。它决定了结合亲和力，也决定了我们后续用何种方法、以何种精度去计算。

一、结合自由能的定义与标准态

设蛋白 P 与配体 L 形成复合物 PL：

P + L ⇌ PL

结合自由能定义为该过程的自由能变化：

ΔG_bind = G(PL) − G(P) − G(L)

与之对应的是结合常数 K_a 与解离常数 K_d：

ΔG_bind = −RT ln K_a = RT ln K_d

其中 R 为气体常数（8.314 J·mol⁻¹·K⁻¹），T 为绝对温度。在 298.15 K 下，ΔG_bind 每变化约 1.36 kcal/mol，K_d 变化一个数量级；每变化约 5.7 kJ/mol，K_d 变化约 10 倍。这是判断计算精度是否“够用”的第一把尺子。

必须强调标准态。K_a 的量纲依赖于浓度标准态，通常取 1 mol/L（即 1 M）。因此 ΔG_bind 是相对于 1 M 标准态的自由能变化，而不是任意浓度下的实际自由能。若在计算中使用了不同的标准浓度，或没有正确处理平动、转动自由能中的标准态项，结果会出现系统性偏差。工程上常见的做法是：在 MM/PBSA 或 LIE 类方法中显式加入标准态校正项，或在热力学积分中通过约束/释放步骤保证参考态一致。

二、焓熵分解：ΔG = ΔH − TΔS

结合自由能可分解为焓与熵两部分：

ΔG_bind = ΔH − TΔS

ΔH 主要来自蛋白–配体之间的非键相互作用（静电、范德华、氢键）以及结合诱导的构象能变化；−TΔS 则包含配体与蛋白去溶剂化、平动/转动自由度损失、构象熵变化以及结合后水分子重排等贡献。焓熵分解的价值在于：它解释了为什么两个亲和力相近的配体可能具有完全不同的物理机制——一个靠强氢键网络，另一个靠疏水埋藏与熵驱动。

但在实际计算中，ΔH 与 −TΔS 往往都是大数相减得到小数，误差放大严重。例如 ΔH = −40 kcal/mol、−TΔS = +32 kcal/mol，得到 ΔG = −8 kcal/mol；若 ΔH 有 2 kcal/mol 误差，ΔG 的相对误差就非常显著。因此，除非方法对熵有可靠估计，否则直接比较 ΔH 或 TΔS 的绝对值意义有限，更稳妥的是比较同一系列配体之间的相对值。

工程经验：在先导化合物优化阶段，优先关注 ΔG 的相对排序；只有在需要解释机制或做骨架跃迁时，才深入分析焓熵补偿。

三、配分函数与系综平均

统计力学把宏观自由能与微观状态联系起来。对于正则系综（NVT），配分函数为：

Q = Σ_i exp(−E_i / k_B T)

自由能与配分函数的关系为：

A = −k_B T ln Q

在恒压系综（NPT）下对应的是吉布斯自由能 G = −k_B T ln Δ，其中 Δ 为等温等压配分函数。任何可观测量的系综平均为：

⟨A⟩ = Σ_i A_i exp(−E_i / k_B T) / Q

这解释了为什么自由能计算不能只看单一构象：它本质上是对所有相关微观状态的加权平均。分子动力学采样越充分，越接近真实系综平均；采样不足时，得到的只是局部构象的平均，误差不可控。

配分函数的分解也直接对应计算中的自由度处理：平动、转动、振动、构象。平动和转动有解析表达式，构象部分必须靠采样。结合过程中，配体从自由态到结合态损失了平动和转动自由度，这部分熵损失在标准态下约为 6–8 kcal/mol 量级（取决于分子量和标准态），是结合自由能中不可忽略的负贡献。

四、计算量与实际可观测量的对应关系

实际实验可观测的是 K_d、IC₅₀、ΔG_bind，而计算直接得到的是能量、力、构象轨迹。二者之间的桥梁是系综平均与标准态校正。常见对应关系如下：

计算量	物理含义	对应实验量
ΔG_bind	标准态结合自由能	K_d、IC₅₀（需换算）
ΔH	结合焓变	ITC 直接测得
−TΔS	结合熵变贡献	ITC 间接得到
⟨E_vdW⟩、⟨E_elec⟩	系综平均非键能	无直接实验对应，用于机制解释
RMSD、氢键占有率	构象稳定性与相互作用持续性	晶体结构、NMR 等结构信息

需要特别注意的是，IC₅₀ 并不等于 K_d，它受实验条件（底物浓度、酶浓度）影响。若要用计算值直接对比 IC₅₀，必须先用 Cheng–Prusoff 关系换算，否则相关性会被系统误差掩盖。

五、从热力学到工程判断

理解上述基础后，面对一个具体计算任务时，应先明确三件事：第一，目标量是绝对 ΔG 还是相对 ΔΔG；第二，标准态和参考态是否一致；第三，熵的贡献是否被合理处理。若只做同系列配体排序，MM/PBSA 或 LIE 可能足够；若要做绝对亲和力预测或骨架差异较大的比较，则需要自由能扰动（FEP）或热力学积分（TI）这类严格方法。

下面的示意图把“热力学定义—标准态—焓熵分解—配分函数—系综平均—可观测量”之间的对应关系集中呈现，便于在工程实践中快速定位误差来源。

把这张关系图记在脑中，后续选择 MM/PBSA、LIE、FEP 或 TI 时，就不会只凭“方法名气”做决定，而是能根据目标量、标准态、熵处理和采样成本做出有依据的工程取舍。