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

Technische Samenvatting: Audit van de Haldane-consistentie in reversibele enzymkinetiek

Probleemstelling
Reversibele enzymkinetische constanten vormen een kwantitatieve intersectie tussen enzymkinetiek en biochemische thermodynamica. Voor een reversibele reactie in een gesloten systeem dicteert microscopische reversibiliteit dat de kinetische constanten aan de Haldane-relatie moeten voldoen, die de impliciete evenwichtsconstante ($K'_{eq,kin}$) in de reactiewet vastlegt. Echter, publieke databases (bijv. SABIO-RK, BREBREDA, TECRDB) bevatten vaak heterogene gegevens waarbij kinetische parameters en thermodynamische evenwichtsconstanten gerapporteerd worden voor verschillende enzymvormen, organismen of condities, of waar de voorwaartse en achterwaartse kinetiek niet in dezelfde studie zijn gerapporteerd. Bijgevolg is er een gebrek aan gecureerde, enkelvoudige "tweezijdige" records die een onafhankelijke audit mogelijk maken om te controleren of gerapporteerde kinetische constanten thermodynamisch consistent zijn met gemeten of geschatte evenwichtsconstanten. Bovs staat er ook geen gestandaardiseerde, gekalibreerdeerde, richting-symmetrische rapportatieschaal beschikbaar om de discrepantie tussen deze twee onafhankelijke vensters op dezelfde reactie te kwantificeren.

Methodologie
De auteurs stellen een workflow voor die gecentreerd is rond een specifieke score-metriek en een rigoureus cureerprotocol:

1. De Haldane-consistentiescore ($C_{Haldane}$):
   Het artikel maakt gebruik van de reciproke kostenfunctie gedefinieerd door Washburn en Zlatanović [1] als een richting-symmetrische audit-schaal:
   $$C_{Haldane} = J(x) = \frac{1}{2}(x + x^{-1}) - 1 = \cosh(\ln x) - 1$$
   waarbij $x = K'_{eq,kin} / K'_{eq,thermo}$.

   • Eigenschappen: De score is nul bij perfecte overeenstemming ($x=1$), symmetrisch voor reciproke fouten ($J(x) = J(1/x)$), en codeert de vrije-energie-discrepantie ($\Delta\Delta G$) in $RT$-eenheden. Het is strikt toenemend in $|\ln x|$, waardoor de ordening van records op basis van de magnitude van de fold-fout behouden blijft.
   • Axiomatische basis: De functionele vorm wordt uniek gekenmerkt door vijf axioma's (reciproque symmetrie, normalisatie, een specifieke compositiewet, continuïteit en eenheidscalibratie), wat het onderscheidt van eenvoudigere metrieken zoals de gekwadrateerde log-fout.
   • Rapporteringsbanden: Continue scores worden in kaart gebracht naar beschrijvende fold-discrepantiebanden (bijv. "binnen een tweevoudige fout" komt overeen met $C_{Haldane} \le 0.25$).

2. Cureerprotocol:
   Een reproduceerbare workflow werd opgezet om "tweezijdige" records te oogsten:

   • Inclusiecriteria: Records moeten reversibele uni-uni (of pseudo-uni-uni met vaste cosubstraten) zijn of een specifieke multi-substraat Haldane-relatie volgen. Cruciaal is dat de voorwaartse en achterwaartse kinetische constanten afkomstig moeten zijn van dezelfde enzymvorm, hetzelfde organisme en binnen dezelfde assay-context binnen één enkele studie, gepaard met een conditie-gematchte thermodynamische comparator (uit TECRDB of eQuilibrator).
   • Exclusiecriteria: Thermodynamisch beperkte globale fits worden uitgesloten omdat zij de Haldane-relatie bij constructie bevredigen en geen onafhankelijke audit bieden. Records met ontbrekende achterwaartse constanten, ambigue identiteiten of mismatchende condities worden verworpen.
   • Mechanisme-afhandeling: Voor bi-bi reacties (ordered, rapid-equilibrium random, ping-pong) worden mechanisme-specifieke Haldane-relaties toegepast om $K'_{eq,kin}$ te berekenen, terwijl de scoring-metriek ongewijzigd blijft.

3. Benchmarking-strategie:
   Aangezien echte gecureerde records geen "ground-truth" foutlabels bevatten, werd een semi-synthetische gelabelde benchmark geconstrueerd. Vertrekkend vanuit 29 "consistente" seed-records (binnen tweevoud), injecteerden de auteurs zes verschillende foutmodi (richting-swaps, eenheids-transcriptiefouten, proton-mismatches, conditie-mismatches, isoform-mismatches en mechanisme/formule-mismatches) om 1.885 gelabelde "inconsistente" gevallen te genereren. Dit maakt het mogelijk om de operationele kenmerken (bijv. AUC) te schatten, voorwaardelijk op de geïnjecteerde taxonomische fouten.

Belangrijkste bijdragen

• Een gekalibreerde rapportageschaal: Het artikel stelt $C_{Haldane}$ niet op als een nieuwe rangschikkingscriteria, maar als een gekalibreerde, richting-symmetrische schaal die zowel de fold-fout als de vrije-energie-discrepantie tegelijkertijd rapporteert, inclusief een gesloten vorm voor de inverse.
• Een gecureerde tweezijdige backbone: De auteurs hebben een dataset samengesteld van 21 geauditeerde single-study records (16 uni-uni, 5 bi-bi) die voldoen aan strikte inclusiecriteria. Dit is het eerste publieke corpus dat aan deze specifieke standaard is gecureerd.
• Een semi-synthetische benchmark: Een gelabelde dataset die de detecteerbaarheid van specifieke foutmodi kwantificeert met behulp van de voorgestelde score, wat de eerste operationele kenmerken biedt voor Haldane-consistentie auditing.
• Reproduceerbaarheid: Alle data, code, protocollen en benchmark-generatoren zijn gearchiveerd met SHA-256 checksums om exacte reproductie te garanderen.

Resultaten

• Demonstratieset: Toegepast op een kleine set handmatig gecureerde records (inclusief fosfoglucose-isomerase, triosefosfaat-isomerase en fumaraat), identificeerde de score succesvol buffer-afhankelijke discrepanties in fumaraat (niet-fosfaatbuffers vertoonden 3,8–4,5-voudige fouten) en bevestigde reciprociteit in racemase-controles.
• De Backbone (21 Records):
  • Distributie: 18 van de 21 records vielen binnen de "binnen-tweevoudige" band ($C_{Haldane} \le 0.25$).
  • Geflagde Records: Drie records werden gemarkeerd als inconsistenties (overschrijding van de tweevoudige grens): Pea phosphoriboisomerase (3,09-voudig), Clostridium scindens D-psicose 3-epimerase (2,88-voudig) en Sinorhizobium D-tagatose 3-epimerase (3,96-voudig).
  • Onafhankelijkheid: Acht records waren "genuin onafhankelijk" (kinetiek past zonder thermodynamische priors vs. aparte evenwichtsmetingen); alle acht waren binnen de tweevoudige band.
  • Chemische reikwijdte: De backbone is sterk geconcentreerd in koolhydraat-isomerases en epimerases; daarom worden de resultaten beschreven als "observaties binnen de familie" voor deze specifieke chemie.
• Benchmark Prestaties:
  • Als continue classifier bereikte de score een Area Under the ROC Curve (AUC) van 0,784 (95% bootstrap CI: 0,725–0,838).
  • Detecteerbaarheid: Grove fouten (eenheidstranscriptie, proton-mismatches) werden gedetecteerd met 100% sensitiviteit. De minst detecteerbare modi waren richting-swaps (wanneer $K'_{eq} \approx 1$) en conditie-mismatches (verkeerde temperatuur), met respectievelijk sensitiviteiten van 0,76 en 0,11 bij de tweevoudige cut-off.
  • De AUC is invariant onder monotone herschaling van $|\ln x|$, wat bevestigt dat het de discriminatie van fold-foutinformatie meet in plaats van een uniek voordeel van de specifieke $C_{Haldane}$ functionele vorm.

Betekenis en Claims
Het artikel claimt bescheiden een reproduceerbare workflow en een beschrijvende empirische referentieverdeling te bieden voor een specifieke, gecureerde steekproef van enzymkinetiek.

• Geen algemene populatieschatting: De auteurs stellen expliciet dat de resultaten (18/21 binnen tweevoud) nog niet kunnen worden geëxtrapoleerd naar alle reversibele enzymkinetiek, aangezien de steekproef chemisch smal is (koolhydraat-isomerases/epimerases) en geen ground-truth foutlabels heeft voor de echte records.
• Haalbaarheidsdemonstratie: De studie dient als een haalbaarheidsdemonstratie dat een tweezijdige audit mogelijk is en dat een aanzienlijk deel van de goed gekarakteriseerde records in deze specifieke klasse thermodynamisch consistent is.
• Curatiedoelen: De drie geflagde records worden niet gepresenteerd als gevalideerde fouten, maar als "curatie-doelen" die verdere inspectie vereisen, aangezien de discrepantie kan voortkomen uit niet-gemodelleerde secundaire kenmerken (bijv. productinhibitie) of problemen met de assemblage van gegevens.
• Operationele bruikbaarheid: De primaire bijdrage is de instelling van een gestandaardiseerde, gekalibreerde schaal ($C_{Haldane}$) en een protocol dat kan worden toegepast op toekomstige data-oogsten zonder methodologische wijzigingen, wat de groei van een grotere, diversere audit-corpus faciliteert.