2. 蛋白质稳定性基础
要预测一个突变会不会让蛋白“变不稳”，必须先明确“稳定”在热力学上到底指什么。蛋白质的天然态（N）与去折叠态（U）之间是一个可逆平衡：N ⇌ U。在给定温度、pH、缓冲液和离子强度下，这个平衡由折叠自由能 ΔGfold = GU − GN 描述。注意符号约定：ΔGfold 为正值表示天然态更稳定，蛋白主要处于折叠状态；ΔGfold 为负值表示去折叠态占优，蛋白趋于解折叠。工程上常说的“稳定性”通常指 ΔGfold 的大小，典型球状蛋白在 25 ℃ 下的 ΔGfold 约为 5～15 kcal/mol，这个数值并不大，相当于十几个氢键或疏水相互作用的净贡献，因此单个突变就足以显著改变平衡。
ΔGfold 与平衡常数 KU = [U]/[N] 的关系为 ΔGfold = −RT ln KU，其中 R = 1.987×10−3 kcal·mol−1·K−1，T 为绝对温度。在 298 K 时，RT ≈ 0.592 kcal/mol。这意味着 ΔGfold 每变化 1 kcal/mol，KU 约变化 e1/0.592 ≈ 5.4 倍。例如某蛋白 ΔGfold = 8.0 kcal/mol，则 KU ≈ e−8.0/0.592 ≈ 1.4×10−6，即约百万分之一的分子处于去折叠态；若突变使其 ΔGfold 降至 6.0 kcal/mol，KU 升至约 4.0×10−5，去折叠态比例增加近 30 倍。这就是为什么稳定性预测必须做到 1 kcal/mol 量级的精度才有实用意义。
突变效应用 ΔΔG 表示，定义为突变体与野生型折叠自由能之差：ΔΔGmut = ΔGfold(mut) − ΔGfold(WT)。按此定义，ΔΔG > 0 表示突变体更不稳定（去稳定突变），ΔΔG < 0 表示突变体更稳定（稳定化突变），ΔΔG ≈ 0 表示中性突变。这一符号约定在文献和预测工具中并不完全统一，有些工具输出 ΔΔGfold = ΔGfold(WT) − ΔGfold(mut)，此时正值代表稳定化。工程实践中必须首先确认所用工具或数据的符号方向，否则会把稳定化突变误判为去稳定突变。本课程统一采用 ΔΔG = ΔGfold(mut) − ΔGfold(WT)，正值去稳定。
从结构角度看，突变通过改变分子内相互作用来影响 ΔGfold。埋藏在疏水核心的突变若把大侧链换成小侧链，会留下空腔并损失范德华接触，通常 ΔΔG 为 +1～+4 kcal/mol；把带电残基引入疏水核心，去溶剂化惩罚极大，ΔΔG 可达 +5 kcal/mol 以上；而表面电荷突变通过优化静电相互作用或引入盐桥，可能产生 −0.3～−1.5 kcal/mol 的稳定化效应。此外，突变还可能改变去折叠态的能量：例如在柔性 loop 区引入脯氨酸会降低去折叠态熵，从而稳定天然态，表现为负 ΔΔG。理解这些结构—能量对应关系，是后续使用 Rosetta、FoldX 或机器学习预测器的基础。
实验上测定 ΔΔG 的经典方法是化学变性剂滴定（如 GdmCl 或尿素）结合圆二色或荧光监测，通过拟合两态模型得到 Cm 和 m 值，再外推至零变性剂浓度得到 ΔGfold(H2O)。也可用热变性 DSC 或 CD 熔解曲线获得 Tm 和 ΔHm，再计算 ΔGfold。需要强调的是，ΔΔG 是状态函数差值，原则上与路径无关，但实验测定值依赖于温度、pH 和缓冲液条件，因此比较突变体与野生型时必须在同一条件下进行。若突变显著改变 m 值或熔解曲线形状，说明两态假设可能不成立，此时简单相减得到的 ΔΔG 需要谨慎解释。
在预测工程中，ΔΔG 是核心标签。训练数据多来自 ProThermDB 等数据库中的实验 ΔΔG，预测器输入为野生型结构、突变位点和突变类型，输出为 ΔΔG 估计值。评估指标常用 Pearson 相关系数、RMSE 和符号准确率。一个常见误区是只关注相关系数而忽略符号准确率：若模型把去稳定突变预测为稳定化，即使数值接近，也会误导后续实验设计。因此，实际筛选稳定化突变时，通常设定 ΔΔG < −0.5 kcal/mol 作为候选阈值，并结合结构检查排除破坏活性位点或二硫键的突变。理解 ΔGfold 与 ΔΔG 的热力学定义和符号约定，是贯穿整门课程的第一块基石。