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

Resumo Técnico: Auditoria da Consistência de Haldane em Cinética Enzimática Reversível

Definição do Problema
As constantes cinéticas de enzimas reversíveis fornecem uma interseção quantitativa entre a cinética enzimática e a termodinâmica bioquímica. Para uma reação reversível em um sistema fechado, a reversibilidade microscópica dita que as constantes cinéticas devem satisfazer a relação de Haldane, que fixa a constante de equilíbrio aparente ($K'_{eq,kin}$) implícida pela lei de velocidade. No entanto, bancos de dados públicos (ex: SABIO-RK, BRENDA, TECRDB) frequentemente contêm dados heterogêneos onde parâmetros cinéticos e constantes de equilíbrio termodinâmico são relatados para diferentes formas de enzimas, organismos ou condições, ou onde a cinética direta e reversa não são relatadas no mesmo estudo. Consequentemente, há uma falta de registros "bilaterais" curados de um único estudo que permitam uma auditoria independente de se os constantes cinéticos relatados são termodinamicamente consistentes com as constantes de equilíbrio medidas ou estimadas. Além disso, não existe uma escala de relato padronizada, calibrada e simétrica em relação à direção para quantificar a discrepância entre essas duas janelas independentes sobre a mesma reação.

Metodologia
Os autores propõem um fluxo de trabalho centrado em uma métrica de pontuação específica e um protocolo de curadoria rigoroso:

1. O Score de Consistência de Haldane ($C_{Haldane}$):
   O artigo utiliza a função de custo recíproca definida por Washburn e Zlatanović [1] como uma escala de auditoria simétrica em relação à direção:
   $$C_{Haldane} = J(x) = \frac{1}{2}(x + x^{-1}) - 1 = \cosh(\ln x) - 1$$
   onde $x = K'_{eq,kin} / K'_{eq,thermo}$.

   • Propriedades: O score é zero em concordância perfeita ($x=1$), simétrico para erros recíprocos ($J(x) = J(1/x)$) e codifica a discrepância de energia livre ($\Delta\Delta G$) em unidades de $RT$. É estritamente crescente em $|\ln x|$, preservando a ordenação dos registros pelo magnitude do erro de dobra (fold-error).
   • Base Axiomática: A forma funcional é unicamente caracterizada por cinco axiomas (simetria recíproca, normalização, uma lei de composição específica, continuidade e calibração de unidade), distinguindo-se de métricas mais simples como o erro logarítmico ao quadrado.
   • Bandas de Relato: Scores contínuos são mapeados para bandas descritivas de discrepância de dobra (ex: "dentro de duas vezes" corresponde a $C_{Haldane} \le 0,25$).

2. Protocolo de Curadoria:
   Um fluxo de trabalho reprodutível foi estabelecido para colher registros "bilaterais":

   • Critérios de Inclusão: Registros devem ser uni-uni reversíveis (ou pseudo-uni-uni com cossubstratos fixos) ou seguir uma relação de Haldane específica para múltiplos substratos. Crucialmente, as constantes cinéticas direta e reversa devem provir da mesma forma de enzima, organismo e contexto de ensaio dentro de um único estudo, pareadas com um comparador termodinâmico de condições correspondentes (do TECRDB ou eQuilibrator).
   • Critérios de Exclusão: Ajustes globais termodinamicamente constrangidos são excluídos porque satisfazem a relação de Haldane por construção e não oferecem auditoria independente. Registros com constantes reversas ausentes, identidades ambíguas ou condições incompatíveis são descartados.
   • Tratamento de Mecanismo: Para reações bi-bi (ordenadas, de equilíbrio rápido aleatório, ping-pong), relações de Haldane específicas de mecanismo são aplicadas para computar $K'_{eq,kin}$, enquanto a métrica de pontuação permanece inalterada.

3. Estratégia de Benchmarking:
   Como os registros curados reais carecem de rótulos de erro de verdade fundamental (ground-truth), um benchmark rotulado semi-sintético foi construído. Partindo de 29 registros "consistentes" semente (dentro de duas vezes), os autores injetaram seis modos distintos de erro (trocas de direção, erros de transcrição de unidade, incompatibilidade de prótons, incompatibilidade de condições, incompatibilidade de isoformas e incompatibilidade de mecanismo/fórmula) para gerar 1.885 casos "inconsistentes" rotulados. Isso permite a estimativa das características operacionais (ex: AUC) condicional à taxonomia de erro injetada.

Principais Contribuições

• Uma Escala de Relato Calibrada: O artigo estabelece $C_{Haldane}$ não como um novo critério de classificação, mas como uma escala calibrada e simétrica em relação à direção que reporta simultaneamente o erro de dobra e a discrepância de energia livre, possuindo uma inversa de forma fechada.
• Um Backbone Bilateral Curado: Os autores montaram um conjunto de 21 registros auditados de um único estudo (16 uni-uni, 5 bi-bi) que atendem aos rigorosos critérios de inclusão. Este é o primeiro corpus público curado para este padrão específico.
• Um Benchmark Semi-Sintético: Um conjunto de dados rotulado quantificando a detectabilidade de modos de erro específicos usando o score proposto, fornecendo as primeiras características operacionais para a auditoria de consistência de Haldane.
• Reprodutibilidade: Todos os dados, códigos, protocolos e geradores de benchmark foram arquivados com checksums SHA-256 para garantir a reprodução exata dos resultados.

Resultos

• Conjunto de Demonstração: Aplicado a um pequeno conjunto de registros curados manualmente (incluindo fosfoglicose isomerase, triosefosfato isomerase e fumarase), o score identificou com sucesso discrepâncias dependentes de tampão na fumarase (tampões não-fosfato mostraram erros de 3,8–4,5 vezes) e confirmou a reciprocidade em controles de racemase.
• O Backbone (21 Registros):
  • Distribuição: 18 dos 21 registros caíram dentro da banda "dentro de duas vezes" ($C_{Haldane} \le 0,25$).
  • Registros Sinalizados: Três registros foram sinalizados como inconsistências (excedendo o limite de duas vezes): fosforiboisomerase de ervilha (3,09 vezes), D-psicose 3-epimerase de Clostridium scindens (2,88 vezes) e D-tagatose 3-epimerase de Sinorhizobium (3,96 vezes).
  • Independência: Oito registros eram "genuinamente independentes" (ajuste cinético sem priors termodinâmicos vs. medições de equilíbrio separadas); todos os oito estavam dentro da banda de duas vezes.
  • Escopo Químico: O backbone está fortemente concentrado em isomerases e epimerases de carboidratos; assim, os resultados são descritos como "observações dentro da mesma família" para esta química específica.
• Desempenho do Benchmark:
  • Como classificador contínuo, o score alcançou uma Área Sob a Curva ROC (AUC) de 0,784 (IC bootstrap 95%: 0,725–0,838).
  • Detectabilidade: Erros grosseiros (transcrição de unidade, incompatibilidade de prótons) foram detectados com 100% de sensibilidade. Os modos menos detectáveis foram trocas de direção (quando $K'_{eq} \approx 1$) e incompatibilidades de condições (temperatura errada), com sensibilidades de 0,76 e 0,11 respectivamente no corte de duas vezes.
  • A AUC é invariante sob o escalonamento monotônico de $|\ln x|$, confirmando que ela mede a discriminabilidade da informação de erro de dobra em vez de uma vantagem única da forma funcional $C_{Haldane}$.

Significância e Alegações
O artigo alega modestamente fornecer um fluxo de trabalho reprodutível e uma distribuição de referência empírica descritiva para um exemplo específico e curado de cinética enzimática.

• Não é uma Estimativa de População Geral: Os autores declaram explicitamente que os resultados (18/21 dentro de duas vezes) ainda não podem ser extrapolados para toda a cinética enzimática reversível, pois a amostra é quimicamente estreita (isomerases/epimerases de carboidratos) e carece de rótulos de erro de verdade fundamental para os registros reais.
• Demonstração de Viabilidade: O estudo serve como uma demonstração de viabilidade de que uma auditoria bilateral é possível e que a grande maioria dos registros bem caracterizados nesta classe específica é termodinamicamente consistente.
• Alvos de Curadoria: Os três registros sinalizados são apresentados não como erros validados, mas como "alvos de curadoria" que requerem maior escrutínio, pois a discrepância pode originar-se de características secundárias não modeladas (ex: inibição pelo produto) ou problemas de montagem de dados.
• Utilidade Operacional: A principal contribuição é o estabelecimento de uma escala padronizada e calibrada ($C_{Haldane}$) e um protocolo que pode ser aplicado a futuras coletas de dados sem mudanças metodológicas, facilitando o crescimento de um corpus de auditoria maior e mais diverso.