2. 蛋白稳定性核心热力学参数

在蛋白稳定性表征工程中,最常被引用的五个量是折叠自由能 ΔG、熔解温度 Tm、聚集起始温度 Tagg、热容变化 ΔCp 以及可逆性判据。它们并非彼此独立的“指标清单”,而是一组由同一热力学框架串联起来的参数:ΔG 描述天然态与去折叠态之间的能量差,Tm 描述该能量差归零时的温度,ΔCp 决定 ΔG 随温度弯曲的曲率,Tagg 描述不可逆聚集开始主导的温度,可逆性则决定上述平衡量是否还能被热力学解释。理解它们之间的换算与制约关系,才能把一次热稳定性实验转化为可执行的工程决策。

一、ΔG 与 ΔCp:稳定性的“高度”和“曲率”

折叠自由能 ΔG = GU − GF,即去折叠态与天然态的自由能差。按惯例,ΔG > 0 表示天然态更稳定,ΔG = 0 对应 Tm。工程上更关心的是 25 ℃ 下的 ΔG25,因为它接近储存与制剂条件。常用近似式(以 Tm 与 ΔCp 为参数)为:

ΔG(T) ≈ ΔH_m · (1 − T/Tm) − ΔCp · [ (Tm − T) + T·ln(T/Tm) ]
其中 ΔH_m = Tm · ΔS_m,ΔS_m 为熔解熵变

ΔCp 是去折叠态与天然态热容之差,通常为正(去折叠态暴露更多疏水表面)。它有两个工程含义:其一,ΔCp 越大,ΔG(T) 曲线越“弯”,冷变性风险越高;其二,ΔCp 决定 ΔH 随温度的变化率,dΔH/dT = ΔCp。因此仅凭一个 Tm 无法判断 25 ℃ 稳定性——两个 Tm 相同的蛋白,ΔCp 大者可能在低温下 ΔG 更低。示例(合理假设):某抗体 Fab 的 Tm = 72 ℃、ΔHm = 480 kJ/mol、ΔCp = 12 kJ/(mol·K),代入上式可得 ΔG25 ≈ 28 kJ/mol;若把 ΔCp 改为 6 kJ/(mol·K),ΔG25 会升到约 40 kJ/mol。这说明 ΔCp 是连接 Tm 与室温稳定性的关键桥梁。

二、Tm 与 Tagg:平衡熔解与不可逆聚集的分界

Tm 是热变性曲线(常用 CD 222 nm、DSF 内源荧光或 DSC)中点对应的温度,反映平衡态下 50% 分子去折叠的温度。Tagg 则是聚集起始温度,通常由静态光散射(SLS)或动态光散射(DLS)随温度上升的拐点确定。二者关系决定了实验解释方式:

  • Tagg > Tm:聚集发生在去折叠之后,Tm 可作为热力学稳定性指标,ΔG 模型适用。
  • Tagg ≈ Tm:熔解与聚集耦合,表观 Tm 可能被聚集拉低,需谨慎解读。
  • Tagg < Tm:聚集先于明显去折叠发生,Tm 不再是纯热力学量,ΔG 拟合失去意义。

因此,在拿到一个 Tm 之前,应先确认 Tagg 与可逆性。若 Tagg 明显低于 Tm,应优先优化缓冲液、pH 或添加辅料抑制聚集,而不是继续用 ΔG 模型做定量比较。

三、可逆性判据:热力学解释的前提

可逆性判据的核心是:升温去折叠后降温,能否回到与初始一致的天然态。常用判断方式包括:

  1. DSC 第二次升温曲线与第一次是否重合;
  2. CD 或荧光降温复性后信号是否恢复至初始值的 90% 以上;
  3. 复性后 SEC-HPLC 单体峰面积是否与初始一致,聚体是否显著增加;
  4. 浓度依赖性与扫描速率依赖性:不可逆过程常表现为 Tm 随扫描速率升高而升高。

只有可逆(或近似可逆)体系,Tm、ΔH、ΔCp 才能用于计算 ΔG。若不可逆,应改用 Tagg、Tonset、残余单体率等操作性指标,并明确标注其非平衡属性。

热力学参数关系与工程决策链 ΔCp 热容变化 决定 ΔG(T) 曲率 Tm 熔解温度 ΔG = 0 的温度点 ΔHm 熔解焓 dΔH/dT = ΔCp ΔG25 折叠自由能 室温稳定性核心量 ΔG(T) 模型 ΔG = ΔHm(1−T/Tm) − ΔCp[(Tm−T)+T·ln(T/Tm)] 前提:可逆平衡 可逆性判据 DSC 二次扫描重合 复性信号恢复 ≥ 90% SEC 单体峰一致 Tagg 聚集起始温度 SLS / DLS 拐点 不可逆过程主导 工程决策分支 ① Tagg > Tm 且可逆 → 用 ΔG25 排序候选分子 ② Tagg ≈ Tm → 优化缓冲液/pH,再测 Tm ③ Tagg < Tm 或不可逆 → 改用 Tagg、残余单体率   不作为 ΔG 定量依据 ΔG > 0 表示天然态更稳定;ΔG = 0 对应 Tm;ΔCp 越大,ΔG(T) 曲率越强,冷变性风险越高 可逆性是使用 ΔG 模型的前提;不可逆体系应改用 Tagg 等操作性指标

四、参数之间的工程换算与决策

把上述参数用于工程决策时,建议按以下顺序推进:

  1. 先测可逆性:用 DSC 或 CD 做一次升温—降温循环,确认是否可逆。不可逆则直接进入 Tagg 优化路线。
  2. 再测 Tm 与 ΔHm:由 DSC 曲线积分得到 ΔHm,由峰位得到 Tm;若用 DSF,需注意染料可能扰动 Tm。
  3. 估算 ΔCp:可由 DSC 曲线基线差直接得到,也可用 ΔCp ≈ 0.05–0.15 kJ/(mol·K) 每残基的经验值粗估,或通过 ΔH(T) 斜率拟合。
  4. 计算 ΔG25:代入 ΔG(T) 模型,得到室温稳定性排序。
  5. 交叉验证 Tagg:用 SLS/DLS 确认 Tagg 是否高于 Tm;若低于,则 ΔG25 仅作参考。

示例决策表(合理假设数据,用于说明判断逻辑):

候选分子 Tm (℃) ΔHm (kJ/mol) ΔCp (kJ/mol·K) ΔG25 (kJ/mol) Tagg (℃) 可逆性 决策
A7248012≈2878可逆优先推进
B7552018≈2270部分可逆需优化缓冲液
C684008≈3066不可逆改用 Tagg 评估

表中 A 的 Tm 不是最高,但 ΔG25 与 Tagg 组合最优;B 的 Tm 最高,却因 ΔCp 大、Tagg 低而存在聚集风险;C 的 ΔG25 看似不错,但不可逆使其失去热力学意义。这正是“参数必须联合解读”的典型场景。

实操提示:若样品量有限,可先用 DSF 快速筛 Tm 与 Tagg,再用 DSC 对优选分子做 ΔCp 与可逆性确认。DSF 的 Tm 通常比 DSC 高 1–3 ℃,跨方法比较时需注明方法。
常见误区:把 Tm 当作唯一稳定性排序依据。Tm 只反映 ΔG = 0 的位置,不反映 25 ℃ 下的 ΔG 大小;忽略 ΔCp 与可逆性,可能把“高 Tm 但易聚集”的分子误判为最优候选。

综上,ΔG、Tm、Tagg、ΔCp 与可逆性构成一条完整的证据链:可逆性决定模型是否可用,ΔCp 与 ΔHm 决定 ΔG(T) 的形状,Tm 标定 ΔG = 0 的位置,Tagg 则划出热力学描述失效的边界。工程决策应在这条链上逐级确认,而不是孤立地比较某一个数值。