1
|
Narayanan B, Weilandt D, Masid M, Miskovic L, Hatzimanikatis V. Rational strain design with minimal phenotype perturbation. Nat Commun 2024; 15:723. [PMID: 38267425 PMCID: PMC10808392 DOI: 10.1038/s41467-024-44831-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 01/08/2024] [Indexed: 01/26/2024] Open
Abstract
Devising genetic interventions for desired cellular phenotypes remains challenging regarding time and resources. Kinetic models can accelerate this task by simulating metabolic responses to genetic perturbations. However, exhaustive design evaluations with kinetic models are computationally impractical, especially when targeting multiple enzymes. Here, we introduce a framework for efficiently scouting the design space while respecting cellular physiological requirements. The framework employs mixed-integer linear programming and nonlinear simulations with large-scale nonlinear kinetic models to devise genetic interventions while accounting for the network effects of these perturbations. Importantly, it ensures the engineered strain's robustness by maintaining its phenotype close to that of the reference strain. The framework, applied to improve the anthranilate production in E. coli, devises designs for experimental implementation, including eight previously experimentally validated targets. We expect this framework to play a crucial role in future design-build-test-learn cycles, significantly expediting the strain design compared to exhaustive design enumeration.
Collapse
Affiliation(s)
- Bharath Narayanan
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
- Department of Oncology, University of Cambridge, Cambridge, CB2 0XZ, UK
| | - Daniel Weilandt
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, 08544, USA
| | - Maria Masid
- Ludwig Institute for Cancer Research, Department of Oncology, University of Lausanne, and Lausanne University Hospital (CHUV), Lausanne, Switzerland
| | - Ljubisa Miskovic
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland.
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland.
| |
Collapse
|
2
|
Marín-Hernández Á, Saavedra E. Metabolic control analysis as a strategy to identify therapeutic targets, the case of cancer glycolysis. Biosystems 2023; 231:104986. [PMID: 37506818 DOI: 10.1016/j.biosystems.2023.104986] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 07/23/2023] [Accepted: 07/25/2023] [Indexed: 07/30/2023]
Abstract
The use of kinetic modeling and metabolic control analysis (MCA) to identify possible therapeutic targets and to investigate the controlling and regulatory mechanisms in cancer glycolysis is here reviewed. The glycolytic pathway has been considered a target to decrease cancer cell growth; however, its occurrence in normal cells makes it difficult to design therapeutic strategies that target this pathway in pathological cells. Notwithstanding, the over-expression of all enzymes and transporters, as well as the expression of isoenzymes with different kinetic and regulatory properties in cancer cells, suggested a different distribution of the control of glycolytic flux than that observed in normal cells. Kinetic models of glycolysis are constructed with enzyme kinetics experimental data, validated with the steady-state metabolite concentrations and glycolytic fluxes; applying MCA, permitted us to identify the steps with the highest control of glycolysis in cancer cells, but low control in normal cells. The cancer glycolysis main controlling steps under several metabolic conditions were: glucose transport, hexokinase and hexose-6-phosphate isomerase (HPI); whereas in normal cells were: the first two and phosphofructokinase-1. HPI is the best therapeutic target because it exerts high control in cancer glycolytic flux, but not in normal cells. Furthermore, kinetic modeling also contributed to identifying new feed-back and feed-forward regulatory loops in cancer cells glycolysis, and to understanding the mode of metabolic action of glycolytic inhibitors. Thus, MCA and metabolic modeling allowed to propose new strategies for inhibiting glycolysis in cancer cells.
Collapse
Affiliation(s)
- Álvaro Marín-Hernández
- Departamento de Bioquímica, Instituto Nacional de Cardiología Ignacio Chávez, Mexico City, 14080, Mexico.
| | - Emma Saavedra
- Departamento de Bioquímica, Instituto Nacional de Cardiología Ignacio Chávez, Mexico City, 14080, Mexico
| |
Collapse
|
3
|
Sahin A, Weilandt DR, Hatzimanikatis V. Optimal enzyme utilization suggests that concentrations and thermodynamics determine binding mechanisms and enzyme saturations. Nat Commun 2023; 14:2618. [PMID: 37147292 PMCID: PMC10162984 DOI: 10.1038/s41467-023-38159-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Accepted: 04/19/2023] [Indexed: 05/07/2023] Open
Abstract
Deciphering the metabolic functions of organisms requires understanding the dynamic responses of living cells upon genetic and environmental perturbations, which in turn can be inferred from enzymatic activity. In this work, we investigate the optimal modes of operation for enzymes in terms of the evolutionary pressure driving them toward increased catalytic efficiency. We develop a framework using a mixed-integer formulation to assess the distribution of thermodynamic forces and enzyme states, providing detailed insights into the enzymatic mode of operation. We use this framework to explore Michaelis-Menten and random-ordered multi-substrate mechanisms. We show that optimal enzyme utilization is achieved by unique or alternative operating modes dependent on reactant concentrations. We find that in a bimolecular enzyme reaction, the random mechanism is optimal over any other ordered mechanism under physiological conditions. Our framework can investigate the optimal catalytic properties of complex enzyme mechanisms. It can further guide the directed evolution of enzymes and fill in the knowledge gaps in enzyme kinetics.
Collapse
Affiliation(s)
- Asli Sahin
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), 1015, Lausanne, Switzerland
| | - Daniel R Weilandt
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), 1015, Lausanne, Switzerland
- Department of Chemistry and Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ, USA
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), 1015, Lausanne, Switzerland.
| |
Collapse
|
4
|
Choudhury S, Moret M, Salvy P, Weilandt D, Hatzimanikatis V, Miskovic L. Reconstructing Kinetic Models for Dynamical Studies of Metabolism using Generative Adversarial Networks. NAT MACH INTELL 2022; 4:710-719. [PMID: 37790987 PMCID: PMC10543203 DOI: 10.1038/s42256-022-00519-y] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Accepted: 07/11/2022] [Indexed: 11/09/2022]
Abstract
Kinetic models of metabolism relate metabolic fluxes, metabolite concentrations and enzyme levels through mechanistic relations, rendering them essential for understanding, predicting and optimizing the behaviour of living organisms. However, due to the lack of kinetic data, traditional kinetic modelling often yields only a few or no kinetic models with desirable dynamical properties, making the analysis unreliable and computationally inefficient. We present REKINDLE (Reconstruction of Kinetic Models using Deep Learning), a deep-learning-based framework for efficiently generating kinetic models with dynamic properties matching the ones observed in cells. We showcase REKINDLE's capabilities to navigate through the physiological states of metabolism using small numbers of data with significantly lower computational requirements. The results show that data-driven neural networks assimilate implicit kinetic knowledge and structure of metabolic networks and generate kinetic models with tailored properties and statistical diversity. We anticipate that our framework will advance our understanding of metabolism and accelerate future research in biotechnology and health.
Collapse
Affiliation(s)
- Subham Choudhury
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Michael Moret
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Pierre Salvy
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Present Address: Cambrium GmBH, Berlin, Germany
| | - Daniel Weilandt
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
- Present Address: Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ USA
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Ljubisa Miskovic
- Laboratory of Computational Systems Biology (LCSB), Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|
5
|
Hurbain J, Thommen Q, Anquez F, Pfeuty B. Quantitative modeling of pentose phosphate pathway response to oxidative stress reveals a cooperative regulatory strategy. iScience 2022; 25:104681. [PMID: 35856027 PMCID: PMC9287814 DOI: 10.1016/j.isci.2022.104681] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Revised: 05/12/2022] [Accepted: 06/23/2022] [Indexed: 01/22/2023] Open
Abstract
Living cells use signaling and regulatory mechanisms to adapt to environmental stresses. Adaptation to oxidative stress involves the regulation of many enzymes in both glycolysis and pentose phosphate pathways (PPP), so as to support PPP-driven NADPH recycling for antioxidant defense. The underlying regulatory logic is investigated by developing a kinetic modeling approach fueled with metabolomics and 13C-fluxomics datasets from human fibroblast cells. Bayesian parameter estimation and phenotypic analysis of models highlight complementary roles for several metabolite-enzyme regulations. Specifically, carbon flux rerouting into PPP involves a tight coordination between the upregulation of G6PD activity concomitant to a decreased NADPH/NADP+ ratio and the differential control of downward and upward glycolytic fluxes through the joint inhibition of PGI and GAPD enzymes. Such functional interplay between distinct regulatory feedbacks promotes efficient detoxification and homeostasis response over a broad range of stress level, but can also explain paradoxical pertubation phenotypes for instance reported for 6PGD modulation in mammalian cells.
Collapse
Affiliation(s)
- Julien Hurbain
- CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, University of Lille, 59000 Lille, France
| | - Quentin Thommen
- CNRS, Inserm, CHU Lille, Institut Pasteur de Lille, UMR9020-U1277 - CANTHER - Cancer Heterogeneity Plasticity and Resistance to Therapies, University of Lille, 59000 Lille, France
| | - Francois Anquez
- CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, University of Lille, 59000 Lille, France
| | - Benjamin Pfeuty
- CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, University of Lille, 59000 Lille, France
| |
Collapse
|
6
|
Ramos JRC, Bissinger T, Genzel Y, Reichl U. Impact of Influenza A Virus Infection on Growth and Metabolism of Suspension MDCK Cells Using a Dynamic Model. Metabolites 2022; 12:metabo12030239. [PMID: 35323683 PMCID: PMC8950586 DOI: 10.3390/metabo12030239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 03/04/2022] [Accepted: 03/09/2022] [Indexed: 11/21/2022] Open
Abstract
Cell cultured-based influenza virus production is a viable option for vaccine manufacturing. In order to achieve a high concentration of viable cells, is requirement to have not only optimal process conditions, but also an active metabolism capable of intracellular synthesis of viral components. Experimental metabolic data collected in such processes are complex and difficult to interpret, for which mathematical models are an appropriate way to simulate and analyze the complex and dynamic interaction between the virus and its host cell. A dynamic model with 35 states was developed in this study to describe growth, metabolism, and influenza A virus production in shake flask cultivations of suspension Madin-Darby Canine Kidney (MDCK) cells. It considers cell growth (concentration of viable cells, mean cell diameters, volume of viable cells), concentrations of key metabolites both at the intracellular and extracellular level and virus titers. Using one set of parameters, the model accurately simulates the dynamics of mock-infected cells and correctly predicts the overall dynamics of virus-infected cells for up to 60 h post infection (hpi). The model clearly suggests that most changes observed after infection are related to cessation of cell growth and the subsequent transition to apoptosis and cell death. However, predictions do not cover late phases of infection, particularly for the extracellular concentrations of glutamate and ammonium after about 12 hpi. Results obtained from additional in silico studies performed indicated that amino acid degradation by extracellular enzymes resulting from cell lysis during late infection stages may contribute to this observed discrepancy.
Collapse
Affiliation(s)
- João Rodrigues Correia Ramos
- Bioprocess Engineering, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany; (T.B.); (Y.G.); (U.R.)
- Correspondence:
| | - Thomas Bissinger
- Bioprocess Engineering, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany; (T.B.); (Y.G.); (U.R.)
| | - Yvonne Genzel
- Bioprocess Engineering, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany; (T.B.); (Y.G.); (U.R.)
| | - Udo Reichl
- Bioprocess Engineering, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany; (T.B.); (Y.G.); (U.R.)
- Institute of Process Engineering, Faculty of Process & Systems Engineering, Otto-von-Guericke University, Universitätsplatz 2, 39106 Magdeburg, Germany
| |
Collapse
|
7
|
Schölzel C, Blesius V, Ernst G, Dominik A. Characteristics of mathematical modeling languages that facilitate model reuse in systems biology: a software engineering perspective. NPJ Syst Biol Appl 2021; 7:27. [PMID: 34083542 PMCID: PMC8175692 DOI: 10.1038/s41540-021-00182-w] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2019] [Accepted: 04/19/2021] [Indexed: 02/06/2023] Open
Abstract
Reuse of mathematical models becomes increasingly important in systems biology as research moves toward large, multi-scale models composed of heterogeneous subcomponents. Currently, many models are not easily reusable due to inflexible or confusing code, inappropriate languages, or insufficient documentation. Best practice suggestions rarely cover such low-level design aspects. This gap could be filled by software engineering, which addresses those same issues for software reuse. We show that languages can facilitate reusability by being modular, human-readable, hybrid (i.e., supporting multiple formalisms), open, declarative, and by supporting the graphical representation of models. Modelers should not only use such a language, but be aware of the features that make it desirable and know how to apply them effectively. For this reason, we compare existing suitable languages in detail and demonstrate their benefits for a modular model of the human cardiac conduction system written in Modelica.
Collapse
Affiliation(s)
- Christopher Schölzel
- Technische Hochschule Mittelhessen - University of Applied Sciences, Giessen, Germany.
| | - Valeria Blesius
- Technische Hochschule Mittelhessen - University of Applied Sciences, Giessen, Germany
| | - Gernot Ernst
- Vestre Viken Hospital Trust, Kongsberg, Norway
- University of Oslo, Oslo, Norway
| | - Andreas Dominik
- Technische Hochschule Mittelhessen - University of Applied Sciences, Giessen, Germany
| |
Collapse
|
8
|
Constraint-based metabolic control analysis for rational strain engineering. Metab Eng 2021; 66:191-203. [PMID: 33895366 DOI: 10.1016/j.ymben.2021.03.003] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Revised: 02/10/2021] [Accepted: 03/02/2021] [Indexed: 11/20/2022]
Abstract
The advancements in genome editing techniques over the past years have rekindled interest in rational metabolic engineering strategies. While Metabolic Control Analysis (MCA) is a well-established method for quantifying the effects of metabolic engineering interventions on flows in metabolic networks and metabolite concentrations, it does not consider the physiological limitations of the cellular environment and metabolic engineering design constraints. We report here a constraint-based framework, Network Response Analysis (NRA), for rational genetic strain design. NRA is cast as a Mixed-Integer Linear Programming problem that integrates MCA, Thermodynamically-based Flux Analysis (TFA), biologically relevant constraints, as well as genome editing restrictions into a comprehensive platform for identifying metabolic engineering targets. We show that the NRA formulation and its core constraints are equivalent to the ones of Flux Balance Analysis (FBA) and TFA, which allows it to be used for a wide range of optimization criteria and with various physiological constraints. We also show how the parametrization and introduction of biological constraints enhance the NRA formulation compared to the classical MCA approach, and we demonstrate its features and its ability to generate multiple alternative optimal strategies given several user-defined boundaries and objectives. In summary, NRA is a sophisticated alternative to classical MCA for rational metabolic engineering that accommodates the incorporation of physiological data at metabolic flux, metabolite concentration, and enzyme expression levels.
Collapse
|
9
|
Hameri T, Fengos G, Hatzimanikatis V. The effects of model complexity and size on metabolic flux distribution and control: case study in Escherichia coli. BMC Bioinformatics 2021; 22:134. [PMID: 33743594 PMCID: PMC7981984 DOI: 10.1186/s12859-021-04066-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2020] [Accepted: 03/08/2021] [Indexed: 12/31/2022] Open
Abstract
Background Significant efforts have been made in building large-scale kinetic models of cellular metabolism in the past two decades. However, most kinetic models published to date, remain focused around central carbon pathways or are built around ad hoc reduced models without clear justification on their derivation and usage. Systematic algorithms exist for reducing genome-scale metabolic reconstructions to build thermodynamically feasible and consistently reduced stoichiometric models. However, it is important to study how network complexity affects conclusions derived from large-scale kinetic models built around consistently reduced models before we can apply them to study biological systems. Results We reduced the iJO1366 Escherichia Coli genome-scale metabolic reconstruction systematically to build three stoichiometric models of different size. Since the reduced models are expansions around the core subsystems for which the reduction was performed, the models are nested. We present a method for scaling up the flux profile and the concentration vector reference steady-states from the smallest model to the larger ones, whilst preserving maximum equivalency. Populations of kinetic models, preserving similarity in kinetic parameters, were built around the reference steady-states and their metabolic sensitivity coefficients (MSCs) were computed. The MSCs were sensitive to the model complexity. We proposed a metric for measuring the sensitivity of MSCs to these structural changes. Conclusions We proposed for the first time a workflow for scaling up the size of kinetic models while preserving equivalency between the kinetic models. Using this workflow, we demonstrate that model complexity in terms of networks size has significant impact on sensitivity characteristics of kinetic models. Therefore, it is essential to account for the effects of network complexity when constructing kinetic models. The presented metric for measuring MSC sensitivity to structural changes can guide modelers and experimentalists in improving model quality and guide synthetic biology and metabolic engineering. Our proposed workflow enables the testing of the suitability of a kinetic model for answering certain study-specific questions. We argue that the model-based metabolic design targets that are common across models of different size are of higher confidence, while those that are different could be the objective of investigations for model improvement. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04066-y.
Collapse
Affiliation(s)
- Tuure Hameri
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), 1015, Lausanne, Switzerland
| | - Georgios Fengos
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), 1015, Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), 1015, Lausanne, Switzerland.
| |
Collapse
|
10
|
Masid M, Ataman M, Hatzimanikatis V. Analysis of human metabolism by reducing the complexity of the genome-scale models using redHUMAN. Nat Commun 2020; 11:2821. [PMID: 32499584 PMCID: PMC7272419 DOI: 10.1038/s41467-020-16549-2] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Accepted: 05/07/2020] [Indexed: 01/31/2023] Open
Abstract
Altered metabolism is associated with many human diseases. Human genome-scale metabolic models (GEMs) were reconstructed within systems biology to study the biochemistry occurring in human cells. However, the complexity of these networks hinders a consistent and concise physiological representation. We present here redHUMAN, a workflow for reconstructing reduced models that focus on parts of the metabolism relevant to a specific physiology using the recently established methods redGEM and lumpGEM. The reductions include the thermodynamic properties of compounds and reactions guaranteeing the consistency of predictions with the bioenergetics of the cell. We introduce a method (redGEMX) to incorporate the pathways used by cells to adapt to the medium. We provide the thermodynamic curation of the human GEMs Recon2 and Recon3D and we apply the redHUMAN workflow to derive leukemia-specific reduced models. The reduced models are powerful platforms for studying metabolic differences between phenotypes, such as diseased and healthy cells.
Collapse
Affiliation(s)
- Maria Masid
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
| | - Meric Ataman
- Computational and Systems Biology, Biozentrum, University of Basel, Basel, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland.
| |
Collapse
|
11
|
Tiwary BK. Computational medicine: quantitative modeling of complex diseases. Brief Bioinform 2020; 21:429-440. [PMID: 30698665 DOI: 10.1093/bib/bbz005] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2018] [Revised: 12/21/2018] [Accepted: 12/26/2018] [Indexed: 12/18/2022] Open
Abstract
Biological complex systems are composed of numerous components that interact within and across different scales. The ever-increasing generation of high-throughput biomedical data has given us an opportunity to develop a quantitative model of nonlinear biological systems having implications in health and diseases. Multidimensional molecular data can be modeled using various statistical methods at different scales of biological organization, such as genome, transcriptome and proteome. I will discuss recent advances in the application of computational medicine in complex diseases such as network-based studies, genome-scale metabolic modeling, kinetic modeling and support vector machines with specific examples in the field of cancer, psychiatric disorders and type 2 diabetes. The recent advances in translating these computational models in diagnosis and identification of drug targets of complex diseases are discussed, as well as the challenges researchers and clinicians are facing in taking computational medicine from the bench to bedside.
Collapse
Affiliation(s)
- Basant K Tiwary
- Centre for Bioinformatics, School of Life Sciences, Pondicherry University, Pondicherry, India
| |
Collapse
|
12
|
Tokic M, Hatzimanikatis V, Miskovic L. Large-scale kinetic metabolic models of Pseudomonas putida KT2440 for consistent design of metabolic engineering strategies. BIOTECHNOLOGY FOR BIOFUELS 2020; 13:33. [PMID: 32140178 PMCID: PMC7048048 DOI: 10.1186/s13068-020-1665-7] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Accepted: 01/22/2020] [Indexed: 05/15/2023]
Abstract
BACKGROUND Pseudomonas putida is a promising candidate for the industrial production of biofuels and biochemicals because of its high tolerance to toxic compounds and its ability to grow on a wide variety of substrates. Engineering this organism for improved performances and predicting metabolic responses upon genetic perturbations requires reliable descriptions of its metabolism in the form of stoichiometric and kinetic models. RESULTS In this work, we developed kinetic models of P. putida to predict the metabolic phenotypes and design metabolic engineering interventions for the production of biochemicals. The developed kinetic models contain 775 reactions and 245 metabolites. Furthermore, we introduce here a novel set of constraints within thermodynamics-based flux analysis that allow for considering concentrations of metabolites that exist in several compartments as separate entities. We started by a gap-filling and thermodynamic curation of iJN1411, the genome-scale model of P. putida KT2440. We then systematically reduced the curated iJN1411 model, and we created three core stoichiometric models of different complexity that describe the central carbon metabolism of P. putida. Using the medium complexity core model as a scaffold, we generated populations of large-scale kinetic models for two studies. In the first study, the developed kinetic models successfully captured the experimentally observed metabolic responses to several single-gene knockouts of a wild-type strain of P. putida KT2440 growing on glucose. In the second study, we used the developed models to propose metabolic engineering interventions for improved robustness of this organism to the stress condition of increased ATP demand. CONCLUSIONS The study demonstrates the potential and predictive capabilities of the kinetic models that allow for rational design and optimization of recombinant P. putida strains for improved production of biofuels and biochemicals. The curated genome-scale model of P. putida together with the developed large-scale stoichiometric and kinetic models represents a significant resource for researchers in industry and academia.
Collapse
Affiliation(s)
- Milenko Tokic
- Laboratory of Computational Systems Biotechnology (LCSB), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology (LCSB), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Ljubisa Miskovic
- Laboratory of Computational Systems Biotechnology (LCSB), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| |
Collapse
|
13
|
Miskovic L, Béal J, Moret M, Hatzimanikatis V. Uncertainty reduction in biochemical kinetic models: Enforcing desired model properties. PLoS Comput Biol 2019; 15:e1007242. [PMID: 31430276 PMCID: PMC6716680 DOI: 10.1371/journal.pcbi.1007242] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2019] [Revised: 08/30/2019] [Accepted: 07/03/2019] [Indexed: 11/18/2022] Open
Abstract
A persistent obstacle for constructing kinetic models of metabolism is uncertainty in the kinetic properties of enzymes. Currently, available methods for building kinetic models can cope indirectly with uncertainties by integrating data from different biological levels and origins into models. In this study, we use the recently proposed computational approach iSCHRUNK (in Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models), which combines Monte Carlo parameter sampling methods and machine learning techniques, in the context of Bayesian inference. Monte Carlo parameter sampling methods allow us to exploit synergies between different data sources and generate a population of kinetic models that are consistent with the available data and physicochemical laws. The machine learning allows us to data-mine the a priori generated kinetic parameters together with the integrated datasets and derive posterior distributions of kinetic parameters consistent with the observed physiology. In this work, we used iSCHRUNK to address a design question: can we identify which are the kinetic parameters and what are their values that give rise to a desired metabolic behavior? Such information is important for a wide variety of studies ranging from biotechnology to medicine. To illustrate the proposed methodology, we performed Metabolic Control Analysis, computed the flux control coefficients of the xylose uptake (XTR), and identified parameters that ensure a rate improvement of XTR in a glucose-xylose co-utilizing S. cerevisiae strain. Our results indicate that only three kinetic parameters need to be accurately characterized to describe the studied physiology, and ultimately to design and control the desired responses of the metabolism. This framework paves the way for a new generation of methods that will systematically integrate the wealth of available omics data and efficiently extract the information necessary for metabolic engineering and synthetic biology decisions. Kinetic models are the most promising tool for understanding the complex dynamic behavior of living cells. The primary goal of kinetic models is to capture the properties of the metabolic networks as a whole, and thus we need large-scale models for dependable in silico analyses of metabolism. However, uncertainty in kinetic parameters impedes the development of kinetic models, and uncertainty levels increase with the model size. Tools that will address the issues with parameter uncertainty and that will be able to reduce the uncertainty propagation through the system are therefore needed. In this work, we applied a method called iSCHRUNK that combines parameter sampling and machine learning techniques to characterize the uncertainties and uncover intricate relationships between the parameters of kinetic models and the responses of the metabolic network. The proposed method allowed us to identify a small number of parameters that determine the responses in the network regardless of the values of other parameters. As a consequence, in future studies of metabolism, it will be sufficient to explore a reduced kinetic space, and more comprehensive analyses of large-scale and genome-scale metabolic networks will be computationally tractable.
Collapse
Affiliation(s)
- Ljubisa Miskovic
- Laboratory of Computational Systems Biology (LCSB), EPFL, CH, Lausanne, Switzerland
| | - Jonas Béal
- Master's Program in Life Sciences and Technology, EPFL, CH, Lausanne, Switzerland
| | - Michael Moret
- Master's Program in Life Sciences and Technology, EPFL, CH, Lausanne, Switzerland
| | | |
Collapse
|
14
|
Miskovic L, Tokic M, Savoglidis G, Hatzimanikatis V. Control Theory Concepts for Modeling Uncertainty in Enzyme Kinetics of Biochemical Networks. Ind Eng Chem Res 2019. [DOI: 10.1021/acs.iecr.9b00818] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Ljubisa Miskovic
- Laboratory of Computational Systems Biology (LCSB), EPFL, CH-1015 Lausanne, Switzerland
| | - Milenko Tokic
- Laboratory of Computational Systems Biology (LCSB), EPFL, CH-1015 Lausanne, Switzerland
| | - Georgios Savoglidis
- Laboratory of Computational Systems Biology (LCSB), EPFL, CH-1015 Lausanne, Switzerland
| | | |
Collapse
|
15
|
Maeda K, Westerhoff HV, Kurata H, Boogerd FC. Ranking network mechanisms by how they fit diverse experiments and deciding on E. coli's ammonium transport and assimilation network. NPJ Syst Biol Appl 2019; 5:14. [PMID: 30993002 PMCID: PMC6461619 DOI: 10.1038/s41540-019-0091-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2018] [Accepted: 03/12/2019] [Indexed: 11/17/2022] Open
Abstract
The complex ammonium transport and assimilation network of E. coli involves the ammonium transporter AmtB, the regulatory proteins GlnK and GlnB, and the central N-assimilating enzymes together with their highly complex interactions. The engineering and modelling of such a complex network seem impossible because functioning depends critically on a gamut of data known at patchy accuracy. We developed a way out of this predicament, which employs: (i) a constrained optimization-based technology for the simultaneous fitting of models to heterogeneous experimental data sets gathered through diverse experimental set-ups, (ii) a 'rubber band method' to deal with different degrees of uncertainty, both in experimentally determined or estimated parameter values and in measured transient or steady-state variables (training data sets), (iii) integration of human expertise to decide on accuracies of both parameters and variables, (iv) massive computation employing a fast algorithm and a supercomputer, (v) an objective way of quantifying the plausibility of models, which makes it possible to decide which model is the best and how much better that model is than the others. We applied the new technology to the ammonium transport and assimilation network, integrating recent and older data of various accuracies, from different expert laboratories. The kinetic model objectively ranked best, has E. coli's AmtB as an active transporter of ammonia to be assimilated with GlnK minimizing the futile cycling that is an inevitable consequence of intracellular ammonium accumulation. It is 130 times better than a model with facilitated passive transport of ammonia.
Collapse
Affiliation(s)
- Kazuhiro Maeda
- Frontier Research Academy for Young Researchers, Kyushu Institute of Technology, Kitakyushu, Fukuoka, Japan
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka Japan
| | - Hans V. Westerhoff
- Department of Molecular Cell Biology, Faculty of Science, VU University Amsterdam, O|2 building, Amsterdam, Netherlands
- Manchester Centre for Integrative Systems Biology, Manchester Interdisciplinary Biocentre, School of Chemical Engineering and Analytical Science, The University of Manchester, Manchester, UK
- Synthetic Systems Biology and Nuclear Organization, Swammerdam Institute for Life Sciences, University of Amsterdam, Amsterdam, Netherlands
| | - Hiroyuki Kurata
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka Japan
- Biomedical Informatics R&D Center, Kyushu Institute of Technology, Iizuka, Fukuoka Japan
| | - Fred C. Boogerd
- Department of Molecular Cell Biology, Faculty of Science, VU University Amsterdam, O|2 building, Amsterdam, Netherlands
| |
Collapse
|
16
|
Hameri T, Fengos G, Ataman M, Miskovic L, Hatzimanikatis V. Kinetic models of metabolism that consider alternative steady-state solutions of intracellular fluxes and concentrations. Metab Eng 2018; 52:29-41. [PMID: 30455161 DOI: 10.1016/j.ymben.2018.10.005] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2018] [Revised: 10/07/2018] [Accepted: 10/22/2018] [Indexed: 11/16/2022]
Abstract
Large-scale kinetic models are used for designing, predicting, and understanding the metabolic responses of living cells. Kinetic models are particularly attractive for the biosynthesis of target molecules in cells as they are typically better than other types of models at capturing the complex cellular biochemistry. Using simpler stoichiometric models as scaffolds, kinetic models are built around a steady-state flux profile and a metabolite concentration vector that are typically determined via optimization. However, as the underlying optimization problem is underdetermined, even after incorporating available experimental omics data, one cannot uniquely determine the operational configuration in terms of metabolic fluxes and metabolite concentrations. As a result, some reactions can operate in either the forward or reverse direction while still agreeing with the observed physiology. Here, we analyze how the underlying uncertainty in intracellular fluxes and concentrations affects predictions of constructed kinetic models and their design in metabolic engineering and systems biology studies. To this end, we integrated the omics data of optimally grown Escherichia coli into a stoichiometric model and constructed populations of non-linear large-scale kinetic models of alternative steady-state solutions consistent with the physiology of the E. coli aerobic metabolism. We performed metabolic control analysis (MCA) on these models, highlighting that MCA-based metabolic engineering decisions are strongly affected by the selected steady state and appear to be more sensitive to concentration values rather than flux values. To incorporate this into future studies, we propose a workflow for moving towards more reliable and robust predictions that are consistent with all alternative steady-state solutions. This workflow can be applied to all kinetic models to improve the consistency and accuracy of their predictions. Additionally, we show that, irrespective of the alternative steady-state solution, increased activity of phosphofructokinase and decreased ATP maintenance requirements would improve cellular growth of optimally grown E. coli.
Collapse
Affiliation(s)
- Tuure Hameri
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Georgios Fengos
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Meric Ataman
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Ljubisa Miskovic
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland.
| |
Collapse
|
17
|
Tokic M, Hadadi N, Ataman M, Neves D, Ebert BE, Blank LM, Miskovic L, Hatzimanikatis V. Discovery and Evaluation of Biosynthetic Pathways for the Production of Five Methyl Ethyl Ketone Precursors. ACS Synth Biol 2018; 7:1858-1873. [PMID: 30021444 DOI: 10.1021/acssynbio.8b00049] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The limited supply of fossil fuels and the establishment of new environmental policies shifted research in industry and academia toward sustainable production of the second generation of biofuels, with methyl ethyl ketone (MEK) being one promising fuel candidate. MEK is a commercially valuable petrochemical with an extensive application as a solvent. However, as of today, a sustainable and economically viable production of MEK has not yet been achieved despite several attempts of introducing biosynthetic pathways in industrial microorganisms. We used BNICE.ch as a retrobiosynthesis tool to discover all novel pathways around MEK. Out of 1325 identified compounds connecting to MEK with one reaction step, we selected 3-oxopentanoate, but-3-en-2-one, but-1-en-2-olate, butylamine, and 2-hydroxy-2-methylbutanenitrile for further study. We reconstructed 3 679 610 novel biosynthetic pathways toward these 5 compounds. We then embedded these pathways into the genome-scale model of E. coli, and a set of 18 622 were found to be the most biologically feasible ones on the basis of thermodynamics and their yields. For each novel reaction in the viable pathways, we proposed the most similar KEGG reactions, with their gene and protein sequences, as candidates for either a direct experimental implementation or as a basis for enzyme engineering. Through pathway similarity analysis we classified the pathways and identified the enzymes and precursors that were indispensable for the production of the target molecules. These retrobiosynthesis studies demonstrate the potential of BNICE.ch for discovery, systematic evaluation, and analysis of novel pathways in synthetic biology and metabolic engineering studies.
Collapse
Affiliation(s)
- Milenko Tokic
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Noushin Hadadi
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Meric Ataman
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Dário Neves
- Institute of Applied Microbiology (iAMB), Aachen Biology and Biotechnology (ABBt), RWTH Aachen University, D-52056 Aachen, Germany
| | - Birgitta E. Ebert
- Institute of Applied Microbiology (iAMB), Aachen Biology and Biotechnology (ABBt), RWTH Aachen University, D-52056 Aachen, Germany
| | - Lars M. Blank
- Institute of Applied Microbiology (iAMB), Aachen Biology and Biotechnology (ABBt), RWTH Aachen University, D-52056 Aachen, Germany
| | - Ljubisa Miskovic
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland
| |
Collapse
|
18
|
Yurkovich JT, Alcantar MA, Haiman ZB, Palsson BO. Network-level allosteric effects are elucidated by detailing how ligand-binding events modulate utilization of catalytic potentials. PLoS Comput Biol 2018; 14:e1006356. [PMID: 30086174 PMCID: PMC6097697 DOI: 10.1371/journal.pcbi.1006356] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Revised: 08/17/2018] [Accepted: 07/09/2018] [Indexed: 01/08/2023] Open
Abstract
Allosteric regulation has traditionally been described by mathematically-complex allosteric rate laws in the form of ratios of polynomials derived from the application of simplifying kinetic assumptions. Alternatively, an approach that explicitly describes all known ligand-binding events requires no simplifying assumptions while allowing for the computation of enzymatic states. Here, we employ such a modeling approach to examine the "catalytic potential" of an enzyme-an enzyme's capacity to catalyze a biochemical reaction. The catalytic potential is the fundamental result of multiple ligand-binding events that represents a "tug of war" among the various regulators and substrates within the network. This formalism allows for the assessment of interacting allosteric enzymes and development of a network-level understanding of regulation. We first define the catalytic potential and use it to characterize the response of three key kinases (hexokinase, phosphofructokinase, and pyruvate kinase) in human red blood cell glycolysis to perturbations in ATP utilization. Next, we examine the sensitivity of the catalytic potential by using existing personalized models, finding that the catalytic potential allows for the identification of subtle but important differences in how individuals respond to such perturbations. Finally, we explore how the catalytic potential can help to elucidate how enzymes work in tandem to maintain a homeostatic state. Taken together, this work provides an interpretation and visualization of the dynamic interactions and network-level effects of interacting allosteric enzymes.
Collapse
Affiliation(s)
- James T. Yurkovich
- Department of Bioengineering, University of California, San Diego, La Jolla, CA, USA
- Bioinformatics and Systems Biology Program, University of California, San Diego, La Jolla, CA, USA
| | - Miguel A. Alcantar
- Department of Bioengineering, University of California, San Diego, La Jolla, CA, USA
| | - Zachary B. Haiman
- Department of Bioengineering, University of California, San Diego, La Jolla, CA, USA
| | - Bernhard O. Palsson
- Department of Bioengineering, University of California, San Diego, La Jolla, CA, USA
- Bioinformatics and Systems Biology Program, University of California, San Diego, La Jolla, CA, USA
- Department of Pediatrics, University of California, San Diego, La Jolla, CA, USA
| |
Collapse
|
19
|
Ataman M, Hernandez Gardiol DF, Fengos G, Hatzimanikatis V. redGEM: Systematic reduction and analysis of genome-scale metabolic reconstructions for development of consistent core metabolic models. PLoS Comput Biol 2017; 13:e1005444. [PMID: 28727725 PMCID: PMC5519011 DOI: 10.1371/journal.pcbi.1005444] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Accepted: 03/01/2017] [Indexed: 11/18/2022] Open
Abstract
Genome-scale metabolic reconstructions have proven to be valuable resources in enhancing our understanding of metabolic networks as they encapsulate all known metabolic capabilities of the organisms from genes to proteins to their functions. However the complexity of these large metabolic networks often hinders their utility in various practical applications. Although reduced models are commonly used for modeling and in integrating experimental data, they are often inconsistent across different studies and laboratories due to different criteria and detail, which can compromise transferability of the findings and also integration of experimental data from different groups. In this study, we have developed a systematic semi-automatic approach to reduce genome-scale models into core models in a consistent and logical manner focusing on the central metabolism or subsystems of interest. The method minimizes the loss of information using an approach that combines graph-based search and optimization methods. The resulting core models are shown to be able to capture key properties of the genome-scale models and preserve consistency in terms of biomass and by-product yields, flux and concentration variability and gene essentiality. The development of these “consistently-reduced” models will help to clarify and facilitate integration of different experimental data to draw new understanding that can be directly extendable to genome-scale models. Reduced models are used commonly to understand the metabolism of organisms and to integrate experimental data for many different studies such as physiology, fluxomics and metabolomics. Without consistent or clear criteria on how these reduced models are actually developed, it is difficult to ensure that they reflect the detailed knowledge that is kept in genome scale metabolic network models (GEMs). The redGEM algorithm presented here allows us to systematically develop consistently reduced metabolic models from their genome-scale counterparts. We applied redGEM for the construction of a core model for E. coli central carbon metabolism. We constructed the core model irJO1366 based on the latest genome-scale E. coli metabolic reconstruction (iJO1366). irJO1366 contains the central carbon pathways and other immediate pathways that must be connected to them for consistency with the iJO1366. irJO1366 can be used to understand metabolism of the organism and also to provide guidance for metabolic engineering purposes. The algorithm is also designed to be modular so that heterologous reactions or pathways can be appended to the core model akin to a “plug-and-play”, synthetic biology approach. The algorithm is applicable to any compartmentalized or non-compartmentalized GEM.
Collapse
Affiliation(s)
- Meric Ataman
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH, Lausanne, Switzerland
| | - Daniel F. Hernandez Gardiol
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH, Lausanne, Switzerland
| | - Georgios Fengos
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH, Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH, Lausanne, Switzerland
- * E-mail:
| |
Collapse
|
20
|
Ataman M, Hatzimanikatis V. lumpGEM: Systematic generation of subnetworks and elementally balanced lumped reactions for the biosynthesis of target metabolites. PLoS Comput Biol 2017; 13:e1005513. [PMID: 28727789 PMCID: PMC5519008 DOI: 10.1371/journal.pcbi.1005513] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2016] [Accepted: 03/31/2017] [Indexed: 01/18/2023] Open
Abstract
In the post-genomic era, Genome-scale metabolic networks (GEMs) have emerged as invaluable tools to understand metabolic capabilities of organisms. Different parts of these metabolic networks are defined as subsystems/pathways, which are sets of functional roles to implement a specific biological process or structural complex, such as glycolysis and TCA cycle. Subsystem/pathway definition is also employed to delineate the biosynthetic routes that produce biomass building blocks. In databases, such as MetaCyc and SEED, these representations are composed of linear routes from precursors to target biomass building blocks. However, this approach cannot capture the nested, complex nature of GEMs. Here we implemented an algorithm, lumpGEM, which generates biosynthetic subnetworks composed of reactions that can synthesize a target metabolite from a set of defined core precursor metabolites. lumpGEM captures balanced subnetworks, which account for the fate of all metabolites along the synthesis routes, thus encapsulating reactions from various subsystems/pathways to balance these metabolites in the metabolic network. Moreover, lumpGEM collapses these subnetworks into elementally balanced lumped reactions that specify the cost of all precursor metabolites and cofactors. It also generates alternative subnetworks and lumped reactions for the same metabolite, accounting for the flexibility of organisms. lumpGEM is applicable to any GEM and any target metabolite defined in the network. Lumped reactions generated by lumpGEM can be also used to generate properly balanced reduced core metabolic models. Stoichiometric models have been used in the area of metabolic engineering and systems biology for many decades. The early examples of these models include simplified ad hoc built metabolic pathways, and biomass compositions. The development of genome scale models (GEMs) brought a standard to metabolic network modeling. However, the vast amount of detailed biochemistry in GEMs makes it necessary to develop methods to manage the complexity in them. In this study, we developed lumpGEM, a tool that can systematically identify subnetworks from metabolic networks that can perform certain tasks, such as biosynthesis of a biomass building block and any other target metabolite. By generating alternative subnetworks, lumpGEM also accounts for the redundancy in metabolic networks. We applied lumpGEM on latest E. coli GEM iJO1366 and identified subnetworks/lumped reactions for every biomass building block defined in its biomass formulation. We also compared the results from lumpGEM with existing knowledge in the literature. The lumped reactions generated by lumpGEM can be used to generate consistently reduced metabolic network models.
Collapse
Affiliation(s)
- Meric Ataman
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology, Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland
- * E-mail:
| |
Collapse
|
21
|
Miskovic L, Alff-Tuomala S, Soh KC, Barth D, Salusjärvi L, Pitkänen JP, Ruohonen L, Penttilä M, Hatzimanikatis V. A design-build-test cycle using modeling and experiments reveals interdependencies between upper glycolysis and xylose uptake in recombinant S. cerevisiae and improves predictive capabilities of large-scale kinetic models. BIOTECHNOLOGY FOR BIOFUELS 2017; 10:166. [PMID: 28674555 PMCID: PMC5485749 DOI: 10.1186/s13068-017-0838-5] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2016] [Accepted: 06/06/2017] [Indexed: 05/28/2023]
Abstract
BACKGROUND Recent advancements in omics measurement technologies have led to an ever-increasing amount of available experimental data that necessitate systems-oriented methodologies for efficient and systematic integration of data into consistent large-scale kinetic models. These models can help us to uncover new insights into cellular physiology and also to assist in the rational design of bioreactor or fermentation processes. Optimization and Risk Analysis of Complex Living Entities (ORACLE) framework for the construction of large-scale kinetic models can be used as guidance for formulating alternative metabolic engineering strategies. RESULTS We used ORACLE in a metabolic engineering problem: improvement of the xylose uptake rate during mixed glucose-xylose consumption in a recombinant Saccharomyces cerevisiae strain. Using the data from bioreactor fermentations, we characterized network flux and concentration profiles representing possible physiological states of the analyzed strain. We then identified enzymes that could lead to improved flux through xylose transporters (XTR). For some of the identified enzymes, including hexokinase (HXK), we could not deduce if their control over XTR was positive or negative. We thus performed a follow-up experiment, and we found out that HXK2 deletion improves xylose uptake rate. The data from the performed experiments were then used to prune the kinetic models, and the predictions of the pruned population of kinetic models were in agreement with the experimental data collected on the HXK2-deficient S. cerevisiae strain. CONCLUSIONS We present a design-build-test cycle composed of modeling efforts and experiments with a glucose-xylose co-utilizing recombinant S. cerevisiae and its HXK2-deficient mutant that allowed us to uncover interdependencies between upper glycolysis and xylose uptake pathway. Through this cycle, we also obtained kinetic models with improved prediction capabilities. The present study demonstrates the potential of integrated "modeling and experiments" systems biology approaches that can be applied for diverse applications ranging from biotechnology to drug discovery.
Collapse
Affiliation(s)
- Ljubisa Miskovic
- Ecole Polytechnique Federale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | | | - Keng Cher Soh
- Ecole Polytechnique Federale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Dorothee Barth
- VTT Technical Research Centre of Finland Ltd, Espoo, Finland
| | | | | | - Laura Ruohonen
- VTT Technical Research Centre of Finland Ltd, Espoo, Finland
| | - Merja Penttilä
- VTT Technical Research Centre of Finland Ltd, Espoo, Finland
| | | |
Collapse
|
22
|
Gábor A, Villaverde AF, Banga JR. Parameter identifiability analysis and visualization in large-scale kinetic models of biosystems. BMC SYSTEMS BIOLOGY 2017; 11:54. [PMID: 28476119 PMCID: PMC5420165 DOI: 10.1186/s12918-017-0428-y] [Citation(s) in RCA: 62] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/04/2017] [Accepted: 04/25/2017] [Indexed: 01/13/2023]
Abstract
Background Kinetic models of biochemical systems usually consist of ordinary differential equations that have many unknown parameters. Some of these parameters are often practically unidentifiable, that is, their values cannot be uniquely determined from the available data. Possible causes are lack of influence on the measured outputs, interdependence among parameters, and poor data quality. Uncorrelated parameters can be seen as the key tuning knobs of a predictive model. Therefore, before attempting to perform parameter estimation (model calibration) it is important to characterize the subset(s) of identifiable parameters and their interplay. Once this is achieved, it is still necessary to perform parameter estimation, which poses additional challenges. Methods We present a methodology that (i) detects high-order relationships among parameters, and (ii) visualizes the results to facilitate further analysis. We use a collinearity index to quantify the correlation between parameters in a group in a computationally efficient way. Then we apply integer optimization to find the largest groups of uncorrelated parameters. We also use the collinearity index to identify small groups of highly correlated parameters. The results files can be visualized using Cytoscape, showing the identifiable and non-identifiable groups of parameters together with the model structure in the same graph. Results Our contributions alleviate the difficulties that appear at different stages of the identifiability analysis and parameter estimation process. We show how to combine global optimization and regularization techniques for calibrating medium and large scale biological models with moderate computation times. Then we evaluate the practical identifiability of the estimated parameters using the proposed methodology. The identifiability analysis techniques are implemented as a MATLAB toolbox called VisId, which is freely available as open source from GitHub (https://github.com/gabora/visid). Conclusions Our approach is geared towards scalability. It enables the practical identifiability analysis of dynamic models of large size, and accelerates their calibration. The visualization tool allows modellers to detect parts that are problematic and need refinement or reformulation, and provides experimentalists with information that can be helpful in the design of new experiments. Electronic supplementary material The online version of this article (doi:10.1186/s12918-017-0428-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Attila Gábor
- BioProcess Engineering Group, IIM-CSIC, Eduardo Cabello 6, Vigo, 36208, Spain.,JRC-COMBINE, RWTH Aachen University, Photonics Cluster, Level 4, Campus-Boulevard 79, Aachen, 52074, Germany
| | | | - Julio R Banga
- BioProcess Engineering Group, IIM-CSIC, Eduardo Cabello 6, Vigo, 36208, Spain.
| |
Collapse
|
23
|
Chiappino-Pepe A, Pandey V, Ataman M, Hatzimanikatis V. Integration of metabolic, regulatory and signaling networks towards analysis of perturbation and dynamic responses. ACTA ACUST UNITED AC 2017. [DOI: 10.1016/j.coisb.2017.01.007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
24
|
Penas DR, González P, Egea JA, Doallo R, Banga JR. Parameter estimation in large-scale systems biology models: a parallel and self-adaptive cooperative strategy. BMC Bioinformatics 2017; 18:52. [PMID: 28109249 PMCID: PMC5251293 DOI: 10.1186/s12859-016-1452-4] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2016] [Accepted: 12/24/2016] [Indexed: 12/02/2022] Open
Abstract
Background The development of large-scale kinetic models is one of the current key issues in computational systems biology and bioinformatics. Here we consider the problem of parameter estimation in nonlinear dynamic models. Global optimization methods can be used to solve this type of problems but the associated computational cost is very large. Moreover, many of these methods need the tuning of a number of adjustable search parameters, requiring a number of initial exploratory runs and therefore further increasing the computation times. Here we present a novel parallel method, self-adaptive cooperative enhanced scatter search (saCeSS), to accelerate the solution of this class of problems. The method is based on the scatter search optimization metaheuristic and incorporates several key new mechanisms: (i) asynchronous cooperation between parallel processes, (ii) coarse and fine-grained parallelism, and (iii) self-tuning strategies. Results The performance and robustness of saCeSS is illustrated by solving a set of challenging parameter estimation problems, including medium and large-scale kinetic models of the bacterium E. coli, bakerés yeast S. cerevisiae, the vinegar fly D. melanogaster, Chinese Hamster Ovary cells, and a generic signal transduction network. The results consistently show that saCeSS is a robust and efficient method, allowing very significant reduction of computation times with respect to several previous state of the art methods (from days to minutes, in several cases) even when only a small number of processors is used. Conclusions The new parallel cooperative method presented here allows the solution of medium and large scale parameter estimation problems in reasonable computation times and with small hardware requirements. Further, the method includes self-tuning mechanisms which facilitate its use by non-experts. We believe that this new method can play a key role in the development of large-scale and even whole-cell dynamic models. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1452-4) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- David R Penas
- BioProcess Engineering Group, IIM-CSIC, Eduardo Cabello 6, Vigo, 36208, Spain
| | - Patricia González
- Computer Architecture Group, Universidade da Coruña, Campus de Elviña s/n, Coruña, 15071 A, Spain
| | - Jose A Egea
- Department of Applied Mathematics and Statistics, Universidad Politécnica de Cartagena, c/ Dr. Fleming s/n, Cartagena, 30202, Spain
| | - Ramón Doallo
- Computer Architecture Group, Universidade da Coruña, Campus de Elviña s/n, Coruña, 15071 A, Spain
| | - Julio R Banga
- BioProcess Engineering Group, IIM-CSIC, Eduardo Cabello 6, Vigo, 36208, Spain.
| |
Collapse
|
25
|
Kiparissides A, Hatzimanikatis V. Thermodynamics-based Metabolite Sensitivity Analysis in metabolic networks. Metab Eng 2016; 39:117-127. [PMID: 27845184 DOI: 10.1016/j.ymben.2016.11.006] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2015] [Revised: 11/07/2016] [Accepted: 11/10/2016] [Indexed: 11/29/2022]
Abstract
The increasing availability of large metabolomics datasets enhances the need for computational methodologies that can organize the data in a way that can lead to the inference of meaningful relationships. Knowledge of the metabolic state of a cell and how it responds to various stimuli and extracellular conditions can offer significant insight in the regulatory functions and how to manipulate them. Constraint based methods, such as Flux Balance Analysis (FBA) and Thermodynamics-based flux analysis (TFA), are commonly used to estimate the flow of metabolites through genome-wide metabolic networks, making it possible to identify the ranges of flux values that are consistent with the studied physiological and thermodynamic conditions. However, unless key intracellular fluxes and metabolite concentrations are known, constraint-based models lead to underdetermined problem formulations. This lack of information propagates as uncertainty in the estimation of fluxes and basic reaction properties such as the determination of reaction directionalities. Therefore, knowledge of which metabolites, if measured, would contribute the most to reducing this uncertainty can significantly improve our ability to define the internal state of the cell. In the present work we combine constraint based modeling, Design of Experiments (DoE) and Global Sensitivity Analysis (GSA) into the Thermodynamics-based Metabolite Sensitivity Analysis (TMSA) method. TMSA ranks metabolites comprising a metabolic network based on their ability to constrain the gamut of possible solutions to a limited, thermodynamically consistent set of internal states. TMSA is modular and can be applied to a single reaction, a metabolic pathway or an entire metabolic network. This is, to our knowledge, the first attempt to use metabolic modeling in order to provide a significance ranking of metabolites to guide experimental measurements.
Collapse
Affiliation(s)
- A Kiparissides
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland; Department of Biochemical Engineering, University College London, WC1E 6BT, London, UK
| | - V Hatzimanikatis
- Laboratory of Computational Systems Biotechnology, École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland.
| |
Collapse
|
26
|
Andreozzi S, Chakrabarti A, Soh KC, Burgard A, Yang TH, Van Dien S, Miskovic L, Hatzimanikatis V. Identification of metabolic engineering targets for the enhancement of 1,4-butanediol production in recombinant E. coli using large-scale kinetic models. Metab Eng 2016; 35:148-159. [PMID: 26855240 DOI: 10.1016/j.ymben.2016.01.009] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2015] [Revised: 12/16/2015] [Accepted: 01/29/2016] [Indexed: 11/24/2022]
Abstract
Rational metabolic engineering methods are increasingly employed in designing the commercially viable processes for the production of chemicals relevant to pharmaceutical, biotechnology, and food and beverage industries. With the growing availability of omics data and of methodologies capable to integrate the available data into models, mathematical modeling and computational analysis are becoming important in designing recombinant cellular organisms and optimizing cell performance with respect to desired criteria. In this contribution, we used the computational framework ORACLE (Optimization and Risk Analysis of Complex Living Entities) to analyze the physiology of recombinant Escherichia coli producing 1,4-butanediol (BDO) and to identify potential strategies for improved production of BDO. The framework allowed us to integrate data across multiple levels and to construct a population of large-scale kinetic models despite the lack of available information about kinetic properties of every enzyme in the metabolic pathways. We analyzed these models and we found that the enzymes that primarily control the fluxes leading to BDO production are part of central glycolysis, the lower branch of tricarboxylic acid (TCA) cycle and the novel BDO production route. Interestingly, among the enzymes between the glucose uptake and the BDO pathway, the enzymes belonging to the lower branch of TCA cycle have been identified as the most important for improving BDO production and yield. We also quantified the effects of changes of the target enzymes on other intracellular states like energy charge, cofactor levels, redox state, cellular growth, and byproduct formation. Independent earlier experiments on this strain confirmed that the computationally obtained conclusions are consistent with the experimentally tested designs, and the findings of the present studies can provide guidance for future work on strain improvement. Overall, these studies demonstrate the potential and effectiveness of ORACLE for the accelerated design of microbial cell factories.
Collapse
Affiliation(s)
- Stefano Andreozzi
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland; Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Anirikh Chakrabarti
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland; Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Keng Cher Soh
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland; Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | | | | | | | - Ljubisa Miskovic
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland; Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland
| | - Vassily Hatzimanikatis
- Laboratory of Computational Systems Biotechnology (LCSB), Swiss Federal Institute of Technology (EPFL), CH-1015 Lausanne, Switzerland; Swiss Institute of Bioinformatics, CH-1015 Lausanne, Switzerland.
| |
Collapse
|
27
|
Zhang C, Hua Q. Applications of Genome-Scale Metabolic Models in Biotechnology and Systems Medicine. Front Physiol 2016; 6:413. [PMID: 26779040 PMCID: PMC4703781 DOI: 10.3389/fphys.2015.00413] [Citation(s) in RCA: 80] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2015] [Accepted: 12/15/2015] [Indexed: 12/21/2022] Open
Abstract
Genome-scale metabolic models (GEMs) have become a popular tool for systems biology, and they have been used in many fields such as industrial biotechnology and systems medicine. Since more and more studies are being conducted using GEMs, they have recently received considerable attention. In this review, we introduce the basic concept of GEMs and provide an overview of their applications in biotechnology, systems medicine, and some other fields. In addition, we describe the general principle of the applications and analyses built on GEMs. The purpose of this review is to introduce the application of GEMs in biological analysis and to promote its wider use by biologists.
Collapse
Affiliation(s)
- Cheng Zhang
- State Key Laboratory of Bioreactor Engineering, East China University of Science and TechnologyShanghai, China
| | - Qiang Hua
- State Key Laboratory of Bioreactor Engineering, East China University of Science and TechnologyShanghai, China
- Shanghai Collaborative Innovation Center for Biomanufacturing TechnologyShanghai, China
| |
Collapse
|