MétaCan
Menu
Back to cohort
Record W4414259289 · doi:10.1101/2025.09.05.674464

Free energy perturbations in enzyme kinetic models reveal cryptic epistasis

2025· preprint· en· W4414259289 on OpenAlexafffund
Karol Buda, Nobuhiko Tokuriki

Bibliographic record

VenuebioRxiv (Cold Spring Harbor Laboratory) · 2025
Typepreprint
Languageen
FieldBiochemistry, Genetics and Molecular Biology
TopicProtein Structure and Dynamics
Canadian institutionsCanada's Michael Smith Genome Sciences CentreUniversity of British Columbia
FundersNatural Sciences and Engineering Research Council of Canada
KeywordsEpistasisKinetic energyParameterized complexityEnzymeMutant

Abstract

fetched live from OpenAlex

Abstract Epistasis—the context-dependence of mutational effects—is a key driver of protein evolution, influencing adaptive pathways and functional diversity. While specific epistasis arises from direct physical interactions between mutations, non-specific epistasis emerges when a non-linear mapping links a protein’s biophysical properties to its function. Enzyme kinetic parameters are assumed to be devoid of non-specific epistasis when measured in vitro , enabling direct connection of epistasis to the enzyme’s structural features. Here, we show that this assumption is incorrect: enzyme catalytic parameters like k cat and K M inherently exhibit non-specific epistasis due to the multi-state nature of the catalytic cycle. Using in silico enzyme models, parameterized by free energies of ground and transition states, we simulated 1000 “mutations” or perturbations to the sub-state free energies within the kinetic ensemble. We then combined these mutations, creating one million double mutants with strictly additive free energy effects. Despite the absence of explicit mutational interactions, we observed substantial epistasis in catalytic parameters; its prevalence and complexity increasing with the number of kinetic states in the mechanism. We derived analytical conditions for the emergence of this form of epistasis in a simple kinetic model, demonstrating that non-specific epistasis depends on the relative values of key microscopic rate constants. Finally, we validated our framework by reanalyzing kinetic data for double mutants in Bacillus cereus β-lactamase I and found that reported specific epistasis in catalytic efficiency was substantially stronger than previously inferred, altering mechanistic interpretations. Our results identify an intrinsic, previously unknown source of epistasis that can distort both the magnitude and sign of mutational effects in enzyme kinetics. We provide theoretical and computational tools for recognizing and correcting for this form of non-specific epistasis, enabling accurate mechanistic inference from kinetic data and improving our understanding of the links between epistasis, structure-function relationships, enzyme evolution, and protein design. Author Summary Enzymes help organisms convert reactants to products through a series of steps; each associated with an energy that dictates how well the enzyme catalyzes a reaction. Enzymes evolve to become more efficient, or catalyze new reactions, through mutations that change the free energies of these steps. Sometimes, the effect of one mutation depends on the presence of another, a phenomenon called epistasis. Epistasis is typically studied by measuring the effect of mutations on standard enzyme parameters under the assumption that changes in these values reflect structural interactions between mutations in the protein. Our study shows that this assumption is misleading. Even when mutations act independently on the energies of steps in an enzyme’s reaction, when combined they can create epistasis. This phenomenon arises from the complex, non-linear relationships between the parameters that define each step in the enzyme reaction and the measurements we obtain during experiments probing enzyme functions. Using computational simulations, we mathematically derive the necessary conditions for this form of epistasis, demonstrate that epistasis increases in prevalence as the enzyme reaction becomes more complex, and apply our model to published experimental data. Our findings urge that researchers should account for these effects before drawing structural conclusions from epistasis.

Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.

How this classification was reachedexpand

Full frame machine prediction

Teacher imitation

Not calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.

metaresearch head score (Codex)0.001
metaresearch head score (Gemma)0.003
Version: metacan-v3-hybrid-931329e0061cValidation status: machine_predicted_unvalidated
Candidate categoriesnone
Consensus categoriesnone
DomainCandidate signal: none · Consensus signal: none
Study designCandidate signal: Simulation or modeling · Consensus signal: Simulation or modeling
GenreCandidate signal: Empirical · Consensus signal: Empirical
Teacher disagreement score0.005
Threshold uncertainty score0.010

Distilled classifier scores by category (both heads)

CategoryCodexGemma
Metaresearch0.0010.003
Meta-epidemiology (narrow)0.0010.000
Meta-epidemiology (broad)0.0010.001
Bibliometrics0.0010.000
Science and technology studies0.0000.001
Scholarly communication0.0010.001
Open science0.0000.001
Research integrity0.0010.001
Insufficient payload (model declined to judge)0.0010.000

Machine scores (provisional)

The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.

Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.

Opus teacher head0.009
GPT teacher head0.211
Teacher spread0.202 · how far apart the two teachers sit on this one work
Validation statusscore_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from it

Classification

machine, unvalidated

Machine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.

The models applied no category: nothing in the taxonomy fit this work.
Study designSimulation or modeling
Domainnot available
GenreEmpirical

How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".

Quick stats

Citations2
Published2025
Admission routes2
Has abstractyes

Explore more

Same venuebioRxiv (Cold Spring Harbor Laboratory)Same topicProtein Structure and DynamicsFrench-language works237,207