Modeling genome-wide enzyme evolution predicts strong epistasis underlying catalytic turnover rates

David Heckmann, Daniel C. Zielinski, Bernhard O. Palsson*

*Corresponding author for this work

Research output: Contribution to journalJournal articleResearchpeer-review

199 Downloads (Pure)


Systems biology describes cellular phenotypes as properties that emerge from the complex interactions of individual system components. Little is known about how these interactions have affected the evolution of metabolic enzymes. Here, we combine genome-scale metabolic modeling with population genetics models to simulate the evolution of enzyme turnover numbers (kcats) from a theoretical ancestor with inefficient enzymes. This systems view of biochemical evolution reveals strong epistatic interactions between metabolic genes that shape evolutionary trajectories and influence the magnitude of evolved kcats. Diminishing returns epistasis prevents enzymes from developing higher kcats in all reactions and keeps the organism far from the potential fitness optimum. Multifunctional enzymes cause synergistic epistasis that slows down adaptation. The resulting fitness landscape allows kcat evolution to be convergent. Predicted kcat parameters show a significant correlation with experimental data, validating our modeling approach. Our analysis reveals how evolutionary forces shape modern kcatS and the whole of metabolism.
Original languageEnglish
Article number5270
JournalNature Communications
Number of pages9
Publication statusPublished - 2018

Fingerprint Dive into the research topics of 'Modeling genome-wide enzyme evolution predicts strong epistasis underlying catalytic turnover rates'. Together they form a unique fingerprint.

Cite this