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

Resumen Técnico: Auditoría de la Consistencia de Haldane en la Cinética Enzimática Reversible

Planteamiento del Problema
Las constantes cinéticas de enzimas reversibles proporcionan una intersección cuantitativa entre la cinética enzimática y la termodinámica bioquímica. Para una reacción reversible en un sistema cerrado, la reversibilidad microscópica dicta que las constantes cinéticas deben satisfacer la relación de Haldane, la cual fija la constante de equilibrio aparente ($K'_{eq,kin}$) implicada por la ley de velocidad. Sin embargo, las bases de datos públicas (por ejemplo, SABIO-RK, BRENDA, TECRDB) suelen contener datos heterogéneos donde los parámetros cinéticos y las constantes de equilibrio termodinámicas se reportan para diferentes formas de la enzima, organismos o condiciones, o donde las cinéticas directa e inversa no se reportan en el mismo estudio. En consecuencia, existe una falta de registros "de dos lados" curados de un solo estudio que permitan una auditoría independiente de si las constantes cinéticas reportadas son termodinámicamente consistentes con las constantes de equilibrio medidas o estimadas. Además, no existe una escala de reporte estandarizada, calibrada y simétrica en dirección para cuantificar la discrepancia entre estas dos ventanas independientes sobre la misma reacción.

Metodología
Los autores proponen un flujo de trabajo centrado en una métrica de puntuación específica y un protocolo de curación riguroso:

1. La Puntuación de Consistencia de Haldane ($C_{Haldane}$):
   El artículo utiliza la función de costo recíproco definida por Washburn y Zlatanović [1] como una escala de auditoría simétrica en dirección:
   $$C_{Haldane} = J(x) = \frac{1}{2}(x + x^{-1}) - 1 = \cosh(\ln x) - 1$$
   donde $x = K'_{eq,kin} / K'_{eq,thermo}$.

   • Propiedades: La puntuación es cero en un acuerdo perfecto ($x=1$), es simétrica para errores recíprocos ($J(x) = J(1/x)$) y codifica la discrepancia de energía libre ($\Delta\Delta G$) en unidades de $RT$. Es estrictamente creciente en $|\ln x|$, preservando el orden de los registros por magnitud de error de veces (fold-error).
   • Base Axiomática: La forma funcional está caracterizada unívocamente por cinco axiomas (simetría recíproca, normalización, una ley de composición específica, continuidad y calibración de unidad), lo que la distingue de métricas más simples como el error de logaritmo al cuadrado.
   • Bandas de Reporte: Las puntuaciones continuas se mapean a bandas descriptivas de discrepancia de veces (por ejemplo, "dentro de un factor de dos" corresponde a $C_{Haldane} \le 0.25$).

2. Protocolo de Curación:
   Se estableció un flujo de trabajo reproducible para recolectar registros "de dos lados":

   • Criterios de Inclusión: Los registros deben ser uni-uni reversibles (o pseudo-uni-uni con cosustratos fijos) o seguir una relación de Haldane específica para multi-sustratos. Crucialmente, las constantes cinéticas directa e inversa deben provenir de la misma forma de enzima, organismo y contexto de ensayo dentro de un mismo estudio, emparejadas con un comparador termodinámico de condiciones coincidentes (de TECRDB o eQuilibrator).
   • Criterios de Exclusión: Se excluyen los ajustes globales termodinámicamente restringidos porque satisfacen la relación de Haldane por construcción y no ofrecen una auditoría independiente. Se descartan los registros con constantes inversas faltantes, identidades ambiguas o condiciones desajustadas.
   • Manejo de Mecanismos: Para reacciones bi-bi (ordenadas, de equilibrio rápido aleatorio, ping-pong), se aplican relaciones de Haldane específicas del mecanismo para calcular $K'_{eq,kin}$, mientras que la métrica de puntuación permanece sin cambios.

3. Estrategia de Benchmarking:
   Dado que los registros curados reales carecen de etiquetas de error de verdad fundamental (ground-truth), se construyó un benchmark semi-sintético etiquetado. Partiendo de 29 registros "consistentes" semilla (dentro de un factor de dos), los autores inyectaron seis modos de error distintos (intercambios de dirección, errores de transcripción de unidades, desajustes de protones, desajustes de condiciones, desajustes de isoformas y desajustes de mecanismo/fórmula) para generar 1,885 casos "inconsistentes" etiquetados. Esto permite la estimación de las características operativas (por ejemplo, AUC) condicional a la taxonomía del error inyectado.

Contribuciones Clave

• Una Escala de Reporte Calibrada: El artículo establece $C_{Haldane}$ no como un nuevo criterio de clasificación, sino como una escala calibrada y simétrica en dirección que reporta simultáneamente el error de veces y la discrepancia de energía libre, con una inversa de forma cerrada.
• Un Esqueleto Curado de Dos Lados: Los autores ensamblaron un conjunto de datos de 21 registros auditados de un solo estudio (16 uni-uni, 5 bi-bi) que cumplen con estrictos criterios de inclusión. Este es el primer corpus público curado para este estándar específico.
• Un Benchmark Semi-Sintético: Un conjunto de datos etiquetado que cuantifica la detectabilidad de modos de error específicos utilizando la puntuación propuesta, proporcionando las primeras características operativas para la auditoría de la consistencia de Haldane.
• Reproducibilidad: Todos los datos, códigos, protocolos y generadores de benchmark están archivados con sumas de verificación SHA-256 para asegurar la reproducción exacta de los resultados.

Resultados

• Conjunto de Demostración: Aplicado a un pequeño conjunto de registros curados manualmente (incluyendo fosfoglucosa isomerasa, triosa fosfato isomerasa y fumarasa), la puntuación identificó con éxito discrepancias dependientes del buffer en la fumarasa (los buffers sin fosfato mostraron errores de 3.8–4.5 veces) y confirmó la reciprocidad en los controles de la racemasa.
• El Esqueleto (21 Registros):
  • Distribución: 18 de los 21 registros se situaron dentro de la banda de "dentro de un factor de dos" ($C_{Haldane} \le 0.25$).
  • Registros Marcados: Tres registros fueron señalados como inconsistencias (excediendo el límite de dos veces): la fosforoboisomerasa de guisante (3.09 veces), la D-psicosa 3-epimerasa de Clostridium scindens (2.88 veces) y la D-tagatosa 3-epimerasa de Sinorhizobium (3.96 veces).
  • Independencia: Ocho registros eran "genuinamente independientes" (la cinética se ajustaba sin prioris termodinámicas vs. mediciones de equilibrio separadas); los ocho estaban dentro de la banda de un factor de dos.
  • Alcance Químico: El esqueleto está fuertemente concentrado en isomerasas y epimerasas de carbohidratos; por lo tanto, los resultados se describen como "observaciones dentro de la familia" para esta química específica.
• Desempeño del Benchmark:
  • Como clasificador continuo, la puntuación alcanzó un Área Bajo la Curva ROC (AUC) de 0.784 (IC de 95% bootstrap: 0.725–0.838).
  • Detectabilidad: Los errores groseros (transcripción de unidades, desajustes de protones) fueron detectados con un 100% de sensibilidad. Los modos menos detectables fueron los intercambios de dirección (cuando $K'_{eq} \approx 1$) y los desajustes de condiciones (temperatura incorrecta), con sensibilidades de 0.76 y 0.11 respectivamente en el corte de dos veces.
  • El AUC es invariante bajo el reescalado monotónico de $|\ln x|$, confirmando que mide la discriminabilidad de la información del error de veces en lugar de una ventaja única de la forma funcional $C_{Haldane}$.

Significancia y Reivindicaciones
El artículo reivindica modestamente proporcionar un flujo de trabajo reproducible y una distribución de referencia empírica descriptiva para una muestra específica y curada de cinética enzimática.

• No es una Estimación de la Población General: Los autores declaran explícitamente que los resultados (18/21 dentro de un factor de dos) no pueden extrapolarse aún a todas las cinéticas enzimáticas reversibles, ya que la muestra es químicamente estrecha (isomerasas/epimerasas de carbohidratos) y carece de etiquetas de verdad fundamental para los registros reales.
• Demostración de Viabilidad: El estudio sirve como una demostración de viabilidad de que una auditoría de dos lados es posible y que una mayoría significativa de los registros bien caracterizados en esta clase específica son termodinámicamente consistentes.
• Objetivos de Curación: Los tres registros marcados se presentan no como errores validados, sino como "objetivos de curación" que requieren mayor escrutinio, ya que la discrepancia podría derivarse de características secundarias no modeladas (por ejemplo, inhibición por producto) o problemas de ensamblaje de datos.
• Utilidad Operacional: La contribución principal es el establecimiento de una escala estandarizada y calibrada ($C_{Haldane}$) y un protocolo que puede aplicarse a futuras recolecciones de datos sin cambios metodológicos, facilitando el crecimiento de un corpus de auditoría más grande y diverso.