Auditing Haldane Consistency in Reversible Enzyme Kinetics: A Curated Two-Sided Backbone and a Labeled Fold-Error Benchmark

技术摘要：审计可逆酶动力学中的 Haldane 一致性

问题陈述
可逆酶动力学常数提供了酶动力学与生物化学热力学之间的定量交集。对于封闭系统中的可逆反应，微观可逆性要求动力学常数必须满足 Haldane 关系，该关系固定了速率方程所暗示的表观平衡常数 ($K'_{eq,kin}$)。然而，公共数据库（如 SABIO-RK、BRENDA、TECRDB）通常包含异质数据，其中动力学参数和热力学平衡常数是针对不同的酶形式、生物体或实验条件报告的，或者正向和反向动力学并非在同一项研究中报告。因此，目前缺乏经过精细整理的、单研究来源的“双侧”记录，无法对报告的动力学常数是否与测定或估算的热力学平衡常数保持热力学一致进行独立审计。此外，也缺乏一种标准化的、经过校准的、方向对称的报告尺度，用于量化这两个观察同一反应的不同窗口之间的差异。

方法论
作者提出了一个以特定评分指标和严格整理协议为核心的工作流：

1. Haldane 一致性评分 ($C_{Haldane}$):
   论文采用了由 Washburn 和 Zlatanović [1] 定义的倒数代价函数作为方向对称的审计尺度：
   $$C_{Haldane} = J(x) = \frac{1}{2}(x + x^{-1}) - 1 = \cosh(\ln x) - 1$$
   其中 $x = K'_{eq,kin} / K'_{eq,thermo}$。

   • 属性： 该评分在完全一致时 ($x=1$) 为零，对于倒数误差具有对称性 ($J(x) = J(1/x)$)，并且以 $RT$单位编码了自由能差异 ($\Delta\Delta G$)。它在$|\ln x|$ 方面是严格递增的，保留了按倍数误差大小排序的顺序。
   • 公理基础： 该函数形式由五个公理唯一确定（倒数对称性、归一化、特定的组合律、连续性和单位校准），这使其区别于更简单的平方对数误差度量。
   • 报告区间： 连续评分被映射到描述性的倍数差异区间（例如，“在两倍以内”对应于 $C_{Haldane} \le 0.25$）。

2. 整理协议：
   建立了一个用于获取“双侧”记录的可重复工作流：

   • 纳入标准： 记录必须是可逆的 uni-uni（或在固定共底物情况下的伪 uni-uni）或遵循特定的多底物 Haldane 关系。至关重要的是，正向和反向动力学常数必须来自同一种酶形式、同一种生物体以及同一种实验检测背景，并与匹配条件的参数进行对比（来自 TECRDB 或 eQuilibrator）。
   • 排除标准： 排除受热力学约束的全局拟合，因为它们通过构造方式满足 Haldane 关系，无法提供独立的审计。缺失反向常数、身份不明或条件不匹配的记录将被剔除。
   • 机制处理： 对于 bi-bi 反应（有序、快速平衡随机、乒乓机制），应用特定机制的 Haldane 关系来计算 $K'_{eq,kin}$，而评分指标保持不变。

3. 基准测试策略：
   由于真实的整理记录缺乏地面真值（ground-truth）误差标签，本文构建了一个半合成的有标签基准。研究从 29 个“一致”种子记录（在两倍以内）开始，注入了六种不同的误差模式（方向颠倒、单位转录错误、质子不匹配、条件不匹配、异构体不匹配、机制/公式不匹配）以生成 1,885 个有标签的“不一致”案例。这使得可以根据注入的误差分类法来估计操作特性（如 AUC）。

核心贡献

• 一个校准的报告尺度： 论文确立了 $C_{Haldane}$ 不仅仅是一个新的排序标准，而是一个校准过的、方向对称的尺度，能够同时报告倍数误差和自由能差异，并具有闭式反函数。
• 一个整理的“双侧骨架”： 作者组建了一个包含 21 条经过审计的单研究记录（16 条 uni-uni，5 条 bi-bi）的数据集。这是第一个专门针对此标准进行整理的公开语料库。
• 一个半合成基准： 一个量化了使用该评分检测特定误差模式能力的有标签数据集，提供了关于 Haldane 一致性审计的操作特性。
• 可重复性： 所有数据、代码、协议和基准生成器均附带 SHA-256 校验和存档，以确保结果的精确复现。

结果

• 演示集： 应用于一组手动整理的记录（包括磷酸葡萄糖异构酶、磷酸三糖异构酶和延胡索酸酶），该评分成功识别了延胡索酸酶的缓冲液依赖性差异（非磷酸盐缓冲液显示出 3.8–4.5 倍的误差），并确认了消旋酶对照组的倒数关系。
• 骨架数据（21 条记录）：
  • 分布： 21 条记录中有 18 条落在“两倍以内”区间（$C_{Haldane} \le 0.25$）。
  • 标记记录： 三条记录被标记为不一致（超过两倍边界）：豌豆磷酸核糖异构酶 (3.09 倍)、Clostridium scindens D-甘露醇 3-差向异构酶 (2.88 倍) 以及 Sinorhizobium D-塔格图糖 3-差向异构酶 (3.96 倍)。
  • 独立性： 八条记录是“真正独立的”（动力学拟合不含热力学先验 vs. 独立的平衡测量）；所有八条均在两倍以内区间。
  • 化学范围： 骨架数据高度集中在碳水化合物异构酶和差向异构酶；因此，结果被描述为针对该特定化学类型的“家族内观察”。
• 基准性能：
  • 作为一个连续分类器，该评分实现的曲线下面积 (AUC) 为 0.784（95% Bootstrap 置信区间：0.725–0.838）。
  • 可检测性： 粗大错误（单位转录、质子不匹配）以 100% 的敏感度被检测到。最难检测的模式是方向颠倒（当 $K'_{eq} \approx 1$ 时）和条件不匹配（错误的温度），在两倍切点下的敏感度分别为 0.76 和 0.11。
  • AUC 在 $|\ln x|$ 的单调重缩放下保持不变，证实了它衡量的是倍数误差信息的判别能力，而非 $C_{Haldane}$ 特定函数形式带来的独特优势。

意义与主张
论文谦虚地声称其提供了一个可重复的工作流和一个针对特定整理样本的描述性经验分布参考。

• 并非全人群估计： 作者明确指出，由于样本在化学性质上较为狭窄（碳水化合物异构酶/差向异构酶）且缺乏真实记录的地面真值误差标签，因此目前的结果不能推及到所有可逆酶动力学。
• 可行性演示： 本研究作为一项可行性演示，证明了进行“双侧审计”是可能的，并且在该特定类别中，绝大多数特征良好的记录在热力学上是一致的。
• 整理目标： 被标记的三条记录并非被视为已证实的错误，而是作为“整理目标”需要进一步审查，因为差异可能源于未建模的次要特征（如产物抑制）或数据组装问题。
• 操作效用： 本文的主要贡献在于建立了一个标准化的、校准过的尺度（$C_{Haldane}$）以及一个可以在未来数据采集中无需改变方法论即可应用的协议，从而促进更大规模、更多样化审计语料库的增长。