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

Résumé technique : Audit de la cohérence de Haldane dans la cinétique enzymatique réversible

Énoncé du problème
Les constantes cinétiques des enzymes réversibles constituent une intersection quantitative entre la cinétique enzymatique et la thermodynamique biochimique. Pour une réaction réversible dans un système fermé, la réversibilité microscopique dicte que les constantes cinétiques doivent satisfaire la relation de Haldane, laquelle fixe la constante d'équilibre apparente ($K'_{eq,kin}$) impliquée par la loi de vitesse. Cependant, les bases de données publiques (par exemple, SABIO-RK, BREBREDA, TECRDB) contiennent souvent des données hétérogènes où les paramètres cinétiques et les constantes d'équilibre thermodynamiques sont rapportés pour des formes enzymatiques, des organismes ou des conditions différents, ou bien où les cinétiques directe et inverse ne sont pas rapportées dans la même étude. Par conséquent, il existe un manque de dossiers « bilatéraux » curés issus d'une étude unique qui permettent un audit indépendant pour vérifier si les constantes cinétiques rapportées sont thermodynamiquement cohérentes avec les constantes d'équilibre mesurées ou estimées. De plus, il n'existe pas d'échelle de rapport standardisée, calibrée et symétrique par rapport à la direction, pour quantifier l'écart entre ces deux fenêtres indépendantes sur la même réaction.

Méthodologie
Les auteurs proposent un flux de travail centré sur une métrique de score spécifique et un protocole de curation rigoureux :

1. Le score de cohérence de Haldane ($C_{Haldane}$) :
   L'article utilise la fonction de coût réciproque définie par Washburn et Zlatanović [1] comme une échelle d'audit symétrique par rapport à la direction :
   $$C_{Haldane} = J(x) = \frac{1}{2}(x + x^{-1}) - 1 = \cosh(\ln x) - 1$$
   où $x = K'_{eq,kin} / K'_{eq,thermo}$.

   • Propriétés : Le score est nul en cas d'accord parfait ($x=1$), symétrique pour les erreurs réciproques ($J(x) = J(1/x)$), et encode la divergence d'énergie libre ($\Delta\Delta G$) en unités de $RT$. Il est strictement croissant en $|\ln x|$, préservant l'ordonnancement des enregistrements par magnitude d'erreur de facteur.
   • Base axiomatique : La forme fonctionnelle est caractérisée de manière unique par cinq axiomes (symétrie réciproque, normalisation, une loi de composition spécifique, continuité et étalonnage unitaire), se distinguant ainsi de métriques plus simples comme l'erreur logarithmique au carré.
   • Bandes de rapport : Les scores continus sont mappés sur des bandes de divergence de facteur descriptives (par exemple, « dans un facteur deux » correspond à $C_{Haldane} \le 0,25$).

2. Protocole de curation :
   Un flux de travail reproductible a été établi pour récolter des enregistrements « bilatéraux » :

   • Critères d'inclusion : Les enregistrements doivent être de type uni-uni réversible (ou pseudo-uni-uni avec cosubstrats fixes) ou suivre une relation de Haldane multi-substrat spécifique. Crucialement, les constantes cinétiques directe et inverse doivent provenir de la même forme enzymatique, du même organisme et du même contexte de dosage au sein d'une même étude, appariées à un comparateur thermodynamique correspondant aux conditions (provenant de TECRDB ou eQuilibrator).
   • Critères d'exclusion : Les ajustements globaux contraints par la thermodynamique sont exclus car ils satisfont la relation de Haldane par construction et n'offrent aucun audit indépendant. Les enregistrements avec des constantes inverses manquantes, des identités ambiguës ou des conditions discordantes sont écartés.
   • Gestion des mécanismes : Pour les réactions bi-bi (ordre, équilibre rapide, ping-pong), des relations de Haldane spécifiques au mécanisme sont appliquées pour calculer $K'_{eq,kin}$, tandis que la métrique de score reste inchangée.

3. Stratégie de benchmarking :
   Puisque les enregistrements curés réels manquent d'étiquettes d'erreur de vérité de terrain (ground-truth), un benchmark semi-synthétique étiqueté a été construit. En partant de 29 enregistrements « cohérents » sources (dans un facteur deux), les auteurs ont injecté six modes d'erreur distincts (inversions de direction, erreurs de transcription d'unités, discordances de protons, discordances de conditions, discordances d'isoformes, et discordances de mécanisme/formule) pour générer 1 885 cas « incohérents » étiquetés. Cela permet d'estimer les caractéristiques de fonctionnement (par exemple, l'AUC) conditionnellement à la taxonomie de l'erreur injectée.

Contributions clés

• Une échelle de rapport calibrée : L'article établit $C_{Haldane}$ non pas comme un nouveau critère de classement, mais comme une échelle calibrée, symétrique par rapport à la direction, qui rapporte simultanément le facteur d'erreur et la divergence d'énergie libre, avec une inverse de forme fermée.
• Un squelette bilatéral curé : Les auteurs ont assemblé un ensemble de 21 enregistrements audités de l'étude unique (16 uni-uni, 5 bi-bi) répondant à des critères d'inclusion stricts. C'est le premier corpus public curé selon ce standard spécifique.
• Un benchmark semi-synthétique : Un ensemble de données étiquetées quantifiant la détectabilité de modes d'erreur spécifiques en utilisant le score proposé, fournissant les premières caractéristiques de fonctionnement pour l'audit de la cohérence de Haldane.
• Reproductibilité : Toutes les données, le code, les protocoles et les générateurs de benchmark sont archivés avec des sommes de contrôle SHA-256 pour garantir une reproduction exacte des résultats.

Résultats

• Ensemble de démonstration : Appliqué à un petit ensemble d'enregistrements curés manuellement (incluant la phosphoglucose isomérase, la triosephosphate isomérase et la fumarase), le score a identifié avec succès des divergences dépendantes du tampon (les tampons non-phosphates montraient des erreurs de facteur 3,8 à 4,5) et a confirmé la réciprocité dans les contrôles de la racémase.
• Le squelette (21 enregistrements) :
  • Distribution : 18 des 21 enregistrements se situaient dans la bande « dans un facteur deux » ($C_{Haldane} \le 0,25$).
  • Enregistrements signalés : Trois enregistrements ont été signalés comme incohérents (dépassant la limite du facteur deux) : la phospho-ribose isomérase de pois (3,09 fois), la D-psicose 3-épimérase de Clostridium scindens (2,88 fois) et la D-tagatose 3-épimérase de Sinorhizobium (3,96 fois).
  • Indépendance : Huit enregistrements étaient « véritablement indépendants » (la cinétique a été ajustée sans priors thermodynamiques vs des mesures d'équilibre séparées) ; tous les huit étaient dans la bande du facteur deux.
  • Portée chimique : Le squelette est fortement concentré dans les isomérases et épimérases de glucides ; ainsi, les résultats sont décrits comme des « observations intra-famille » pour cette chimie spécifique.
• Performance du benchmark :
  • En tant que classificateur continu, le score a atteint une Aire Sous la Courbe ROC (AUC) de 0,784 (IC bootstrap 95 % : 0,725–0,838).
  • Détectabilité : Les erreurs grossières (transcription d'unités, discordances de protons) ont été détectées avec une sensibilité de 100 %. Les modes les moins détectables étaient les inversions de direction (lorsque $K'_{eq} \approx 1$) et les discordances de conditions (mauvaise température), avec des sensibilités respectives de 0,76 et 0,11 au seuil du facteur deux.
  • L'AUC est invariante sous un rescalage monotone de $|\ln x|$, confirmant qu'elle mesure la discriminabilité de l'information de facteur d'erreur plutôt qu'un avantage unique de la forme fonctionnelle $C_{Haldane}$.

Signification et affirmations
L'article affirme modestement fournir un flux de travail reproductible et une distribution de référence empirique descriptive pour un échantillon curé spécifique de la cinétique enzymatique.

• Pas une estimation de la population générale : Les auteurs précisent explicitement que les résultats (18/21 dans un facteur deux) ne peuvent pas encore être extrapolés à toute la cinétique enzymatique réversible, car l'échantillon est chimiquement étroit (isomérases/épimérases de glucides) et manque d'étiquettes d'erreur de vérité de terrain pour les enregistrements réels.
• Démonstration de faisabilité : L'étude sert de démonstration de faisabilité prouvant qu'un audit bilatéral est possible et qu'une majorité significative d'enregistrements bien caractérisés dans cette classe spécifique sont thermodynamiquement cohérents.
• Cibles de curation : Les trois enregistrements signalés sont présentés non pas comme des erreurs validées, mais comme des « cibles de curation » nécessitant un examen plus approfondi, car la divergence pourrait provenir de caractéristiques secondaires non modélisées (par exemple, l'inhibition par le produit) ou de problèmes d'assemblage de données.
• Utilité opérationnelle : La principale contribution est l'établissement d'une échelle standardisée et calibrée ($C_{Haldane}$) et d'un protocole qui peut être appliqué à de futures récoltes de données sans changements méthodologiques, facilitant la croissance d'un corpus d'audit plus large et plus diversifié.