Simultaneous application of enzyme and thermodynamic constraints to metabolic models using an updated Python implementation of GECKO

Abstract

GEnome-scale Metabolic (GEM) models are knowledge bases of the reactions and metabolites of a particular organism. These GEM models allow for the simulation of the metabolism - e.g. calculating growth and production yields - based on the stoichiometry, reaction directionality and uptake rates of the metabolic network. Over the years, several extensions have been added to take into account other actors in metabolism, going beyond pure stoichiometry. One such extension is enzyme-constraint models, which enable the integration of kinetic data and proteomics data into GEM models. Given its relatively recent formulation, there are still challenges in standardization and data reconciliation between the model and the experimental measurements. In this work, we present geckopy 3.0 (Genome-scale model Enzyme Constraints, using Kinetics and Omics in python), an actualization from scratch of the previous python implementation of the same name. This update tackles the aforementioned challenges, in an effort to reach maturity in enzyme-constraint modeling. With the new geckopy, proteins are typed in the SBML document, taking advantage of the SBML Groups extension, in compliance with community standards. Additionally, a suite of relaxation algorithms - in the form of linear and mixed-integer linear programming problems - has been added to facilitate reconciliation of raw proteomics data with the metabolic model. Several functionalities to integrate experimental data were implemented, including an interface layer with pytfa for the usage of thermodynamics and metabolomics constraints. Finally, the relaxation algorithms were benchmarked against public proteomics datasets in Escherichia coli for different conditions, revealing targets for improving the enzyme constrained model and/or the proteomics pipeline.