Towards kinetic modeling of genome-scale metabolic networks without sacrificing stoichiometric, thermodynamic and physiology constraints

Mathematical modeling is an essential tool for a comprehensive understanding of cell metabolism and its interactions with the environmental and process conditions. Recent developments in the construction and analysis of stoichiometric models made it possible to define limits on steady-state metabolic behavior using flux balance analysis. However, detailed information about enzyme kinetics and enzyme regulation is needed to formulate kinetic models that can accurately capture the dynamic metabolic responses. The use of mechanistic enzyme kinetics is a difficult task due to uncertainty in the kinetic properties of enzymes. Therefore, the majority of recent works consider only the mass action kinetics for the reactions in the metabolic networks. In this work, we applied the ORACLE framework and constructed a large-scale, mechanistic kinetic model of optimally grown E. coli. We investigated the complex interplay between stoichiometry, thermodynamics, and kinetics in determining the flexibility and capabilities of metabolism. Our results indicate that enzyme saturation is a necessary consideration in modeling metabolic networks and it extends the feasible ranges of the metabolic fluxes and metabolite concentrations. Our results further suggest that the enzymes in metabolic networks have evolved to function at different saturation states to ensure greater flexibility and robustness of the cellular metabolism.

Identification of Feasible Metabolic Fluxes and Metabolite Concentrations using Large-scale Kinetic Models

Anirikh Chakrabarti, Vassily Hatzimanikatis, Ljubisa Miskovic, Keng Cher Soh

The constraints imposed during modeling must satisfy biologically representative phenotypes of the studied organism. Simultaneously, identification and analysis of these constraints enhances our understanding of the evolution/operational paradigms of the organism. It was postulated by Varma and Palsson[1] that it is possible to define limits on metabolic behavior using flux balance analysis, but in order to accurately capture the metabolic responses, detailed information about enzyme kinetics and their regulation is needed. Since development of mechanistic kinetic models is a difficult task due to uncertainty in kinetic properties of enzymes, a substantial number of recent works consider only the mass action (MA) term in their model formulation. As kinetics is one of crucial factors in governing the metabolic capabilities of a cell, i.e. realizable metabolic flux and concentration states, considering only the mass action term does not necessarily provide a realistic description of the feasible space of fluxes and concentrations. In this work, using the ORACLE[2] framework, we constructed a large-scale mechanistic kinetic model of optimally grown E. coli that considers the enzyme saturations as observed in biological systems. Using this model, we performed an analysis of the complex interplay between stoichiometry, thermodynamics, and kinetics in determining flexibility and capabilities of metabolic networks. Our analysis indicates that enzyme saturation is an important and necessary consideration in modeling metabolic networks. Extended ranges of feasibility, both in the space of metabolic fluxes and metabolite concentrations, of kinetic models involving the enzyme saturation suggests that the enzymes in metabolic networks have evolved to function at different saturation states so as to ensure higher flexibility and robustness of the cell.

2013

Mapping and Analysis of the Complexity of Lipid Metabolism for Applications in Health and Sustainability

Vassily Hatzimanikatis, Sofia Tsouka

Metabolism is a network of biotransformations that sustain life in cells. Over the past decades, metabolism has been studied in multiple contexts, ranging from medicine to industrial manufacturing. In particular, lipid metabolism has been linked to various physiopathologies, and its study could lead to new diagnostics and treatments. Metabolic networks are extremely complex and a systems approach is essential in deciphering the different levels of cellular organization and regulation that an organism possesses. Genome-scale metabolic models (GEMs) encompass all the available information for an organism, though their large size and complexity introduce a lot of uncertainty in the accuracy of predictions. GEM reductions are usually done in an ad hoc manner to produce context-specific models, and cannot serve multiple studies. The study of regulatory mechanisms requires the development of kinetic models that can accurately capture dynamic responses to perturbations. To this end, constraint-based models can be enhanced through the integration of kinetic information. Constructing consistent large-scale kinetic models is a challenging endeavor, since detailed knowledge about regulatory mechanisms and kinetic parameter values is scarce.Examining the mechanisms of enzymatic and metabolic restructuring in states of mutation could increase our fundamental understanding of how diseases evolve, and facilitate phenotype mapping. Metabolic Control Analysis (MCA) has been a well-established tool for the prediction and evaluation of genetic modification strategies. However, it fails to account for the physiological limitations of an organism and can often lead to unrealistic predictions.In this thesis, we developed computational models, tools, and methodologies to facilitate the study of lipids and their regulatory mechanisms. Firstly, we constructed and curated a metabolic model that focuses on lipid metabolism, through the integration of detailed lipid pathways into a GEM of S. cerevisiae. This model was then systematically reduced around these pathways to provide a more manageable model size for complex studies. We show that this model is as consistent and inclusive as other yeast GEMs, and can be used as a scaffold for integrating lipidomics data to improve predictions in studies of lipid-related biological functions. Secondly, we used this model as a basis to build a large-scale kinetic model. Enzymes in the lipid metabolism are typically promiscuous and multifunctional, giving rise to enzymatic coupling. To accurately capture these properties, we assigned suitable kinetic rate expressions to the reactions in the network. We generated populations of kinetic parameter sets through a sampling-based workflow and we demonstrate how the consideration of enzymatic coupling is essential for factual predictions. Thirdly, we developed a constraint-based formulation that utilizes MCA-based control coefficients for the consistent derivation of metabolic engineering strategies. We show how the parametrization and introduction of biological constraints enhances this formulation in comparison to the classical MCA approach, and we highlight its ability to generate alternative optimal strategies. Fourthly, we defined and introduced additional mathematical objectives and constraints to this formulation to enable the mapping of enzymes that are responsible for different phenotypes. Concluding, we discuss the contribution and potential applications of this thesis.

EPFL2020

Large-scale dynamic models of metabolic networks

Anirikh Chakrabarti, Vassily Hatzimanikatis, Ljubisa Miskovic, Keng Cher Soh

The development and application of methods quantifying the metabolic fluxes such as Flux Balance Analysis (FBA) has been one of the driving forces behind the successful growth of metabolic engineering. However, applications of FBA have been limited by the fact that FBA does not take into account kinetic properties of the network and therefore it cannot be used to identify rate-limiting steps and comprehend time course evolutions of the system. Dynamic mathematical descriptions of the metabolism offer a large advantage compared to constraint-based stoichiometric models, but unfortunately their development comes with inevitable difficulties due to: (i) structural uncertainties, such as incomplete knowledge about stoichiometry or about kinetic laws of the enzymes, and (ii) quantitative/parametric uncertainties such as lack of knowledge concerning kinetic parameters. Recent developments and vast resources of curated genome scale metabolic networks address to a great extent the issues around stoichiometric uncertainty. However, the knowledge about kinetic rate laws and in particular their parameters is to these days still limited. In this contribution, starting from large-scale stoichiometric models we use the ORACLE (Optimization and Risk Analysis of Complex Living Entities) framework that integrates available information in a set of stable log-linear kinetic models sharing the same steady state. These models are used to compute kinetic parameters of the enzymatic mechanisms in the metabolic network, e.g. the maximal velocities, Vmax, and Michaelis constants, Km. Using these kinetic parameters, we systematically develop populations of stable dynamic models having the same steady-state as the log-linear ones. The estimated parameters are comparable to the experimental information as seen in BRENDA and other databases. These non-linear estimations about the stable state can be therefore used to analyze properties of the system upon large perturbations and investigate time course evolutions in and around this steady state. We demonstrate the capabilities of the proposed approach by building a dynamic E. coli core model that includes ca. 200 metabolites and more than 400 reactions. We discuss the strengths and limitations of this approach and possible avenues for development.

2012

Large-scale dynamic models of metabolic networks

Anirikh Chakrabarti, Vassily Hatzimanikatis, Ljubisa Miskovic, Keng Cher Soh

2012

Identification of Feasible Metabolic Fluxes and Metabolite Concentrations using Large-scale Kinetic Models

Anirikh Chakrabarti, Vassily Hatzimanikatis, Ljubisa Miskovic, Keng Cher Soh

2013

Mapping and Analysis of the Complexity of Lipid Metabolism for Applications in Health and Sustainability

Vassily Hatzimanikatis, Sofia Tsouka

EPFL2020