1
|
Faure L, Mollet B, Liebermeister W, Faulon JL. A neural-mechanistic hybrid approach improving the predictive power of genome-scale metabolic models. Nat Commun 2023; 14:4669. [PMID: 37537192 PMCID: PMC10400647 DOI: 10.1038/s41467-023-40380-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/01/2022] [Accepted: 07/19/2023] [Indexed: 08/05/2023] Open
Abstract
Constraint-based metabolic models have been used for decades to predict the phenotype of microorganisms in different environments. However, quantitative predictions are limited unless labor-intensive measurements of media uptake fluxes are performed. We show how hybrid neural-mechanistic models can serve as an architecture for machine learning providing a way to improve phenotype predictions. We illustrate our hybrid models with growth rate predictions of Escherichia coli and Pseudomonas putida grown in different media and with phenotype predictions of gene knocked-out Escherichia coli mutants. Our neural-mechanistic models systematically outperform constraint-based models and require training set sizes orders of magnitude smaller than classical machine learning methods. Our hybrid approach opens a doorway to enhancing constraint-based modeling: instead of constraining mechanistic models with additional experimental measurements, our hybrid models grasp the power of machine learning while fulfilling mechanistic constrains, thus saving time and resources in typical systems biology or biological engineering projects.
Collapse
Affiliation(s)
- Léon Faure
- MICALIS Institute, INRAE, AgroParisTech, University of Paris-Saclay, 78350, Jouy-en-Josas, France
| | - Bastien Mollet
- Ecole Normale Supérieure of Lyon, 69342, Lyon, France
- UMR MIA, INRAE, AgroParisTech, University of Paris-Saclay, 91120, Palaiseau, France
| | | | - Jean-Loup Faulon
- MICALIS Institute, INRAE, AgroParisTech, University of Paris-Saclay, 78350, Jouy-en-Josas, France.
- Manchester Institute of Biotechnology, University of Manchester, Manchester, M1 7DN, UK.
| |
Collapse
|
2
|
Razavi M, Saberi Fathi SM, Tuszynski JA. The Effect of the Protein Synthesis Entropy Reduction on the Cell Size Regulation and Division Size of Unicellular Organisms. ENTROPY 2022; 24:e24010094. [PMID: 35052120 PMCID: PMC8775074 DOI: 10.3390/e24010094] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Revised: 12/04/2021] [Accepted: 12/28/2021] [Indexed: 12/14/2022]
Abstract
The underlying mechanism determining the size of a particular cell is one of the fundamental unknowns in cell biology. Here, using a new approach that could be used for most of unicellular species, we show that the protein synthesis and cell size are interconnected biophysically and that protein synthesis may be the chief mechanism in establishing size limitations of unicellular organisms. This result is obtained based on the free energy balance equation of protein synthesis and the second law of thermodynamics. Our calculations show that protein synthesis involves a considerable amount of entropy reduction due to polymerization of amino acids depending on the cytoplasmic volume of the cell. The amount of entropy reduction will increase with cell growth and eventually makes the free energy variations of the protein synthesis positive (that is, forbidden thermodynamically). Within the limits of the second law of thermodynamics we propose a framework to estimate the optimal cell size at division.
Collapse
Affiliation(s)
- Mohammad Razavi
- Department of Physics, Faculty of Science, Ferdowsi University of Mashhad, Mashhad 9177948974, Iran;
| | - Seyed Majid Saberi Fathi
- Department of Physics, Faculty of Science, Ferdowsi University of Mashhad, Mashhad 9177948974, Iran;
- Correspondence:
| | - Jack Adam Tuszynski
- Department of Oncology, University of Alberta, Edmonton, AB T6G 1Z2, Canada;
- Department of Physics, University of Alberta, Edmonton, AB T6G 2E1, Canada
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, 10129 Torino, Italy
| |
Collapse
|
3
|
Kishk A, Pacheco MP, Sauter T. DCcov: Repositioning of drugs and drug combinations for SARS-CoV-2 infected lung through constraint-based modeling. iScience 2021; 24:103331. [PMID: 34723158 PMCID: PMC8536485 DOI: 10.1016/j.isci.2021.103331] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Revised: 07/29/2021] [Accepted: 10/19/2021] [Indexed: 12/15/2022] Open
Abstract
The 2019 coronavirus disease (COVID-19) became a worldwide pandemic with currently no approved effective antiviral drug. Flux balance analysis (FBA) is an efficient method to analyze metabolic networks. Here, FBA was applied on human lung cells infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to reposition metabolic drugs and drug combinations against the virus replication within the host tissue. Making use of expression datasets of infected lung tissue, genome-scale COVID-19-specific metabolic models were reconstructed. Then, host-specific essential genes and gene pairs were determined through in silico knockouts that permit reducing the viral biomass production without affecting the host biomass. Key pathways that are associated with COVID-19 severity in lung tissue are related to oxidative stress, ferroptosis, and pyrimidine metabolism. By in silico screening of Food and Drug Administration (FDA)-approved drugs on the putative disease-specific essential genes and gene pairs, 85 drugs and 52 drug combinations were predicted as promising candidates for COVID-19 (https://github.com/sysbiolux/DCcov).
Collapse
Affiliation(s)
- Ali Kishk
- Systems Biology Group, Department of Life Sciences and Medicine, University of Luxembourg, 4367 Esch-sur-Alzette, Luxembourg
| | - Maria Pires Pacheco
- Systems Biology Group, Department of Life Sciences and Medicine, University of Luxembourg, 4367 Esch-sur-Alzette, Luxembourg
| | - Thomas Sauter
- Systems Biology Group, Department of Life Sciences and Medicine, University of Luxembourg, 4367 Esch-sur-Alzette, Luxembourg
| |
Collapse
|
4
|
Panikov NS. Genome-Scale Reconstruction of Microbial Dynamic Phenotype: Successes and Challenges. Microorganisms 2021; 9:2352. [PMID: 34835477 PMCID: PMC8621822 DOI: 10.3390/microorganisms9112352] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 10/18/2021] [Accepted: 10/27/2021] [Indexed: 12/04/2022] Open
Abstract
This review is a part of the SI 'Genome-Scale Modeling of Microorganisms in the Real World'. The goal of GEM is the accurate prediction of the phenotype from its respective genotype under specified environmental conditions. This review focuses on the dynamic phenotype; prediction of the real-life behaviors of microorganisms, such as cell proliferation, dormancy, and mortality; balanced and unbalanced growth; steady-state and transient processes; primary and secondary metabolism; stress responses; etc. Constraint-based metabolic reconstructions were successfully started two decades ago as FBA, followed by more advanced models, but this review starts from the earlier nongenomic predecessors to show that some GEMs inherited the outdated biokinetic frameworks compromising their performances. The most essential deficiencies are: (i) an inadequate account of environmental conditions, such as various degrees of nutrients limitation and other factors shaping phenotypes; (ii) a failure to simulate the adaptive changes of MMCC (MacroMolecular Cell Composition) in response to the fluctuating environment; (iii) the misinterpretation of the SGR (Specific Growth Rate) as either a fixed constant parameter of the model or independent factor affecting the conditional expression of macromolecules; (iv) neglecting stress resistance as an important objective function; and (v) inefficient experimental verification of GEM against simple growth (constant MMCC and SGR) data. Finally, we propose several ways to improve GEMs, such as replacing the outdated Monod equation with the SCM (Synthetic Chemostat Model) that establishes the quantitative relationships between primary and secondary metabolism, growth rate and stress resistance, process kinetics, and cell composition.
Collapse
Affiliation(s)
- Nicolai S Panikov
- Department of Chemistry and Chemical Biology, Northeastern University, 360 Huntington Ave., Boston, MA 02115, USA
| |
Collapse
|
5
|
Intelligent host engineering for metabolic flux optimisation in biotechnology. Biochem J 2021; 478:3685-3721. [PMID: 34673920 PMCID: PMC8589332 DOI: 10.1042/bcj20210535] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2021] [Revised: 09/22/2021] [Accepted: 09/24/2021] [Indexed: 12/13/2022]
Abstract
Optimising the function of a protein of length N amino acids by directed evolution involves navigating a 'search space' of possible sequences of some 20N. Optimising the expression levels of P proteins that materially affect host performance, each of which might also take 20 (logarithmically spaced) values, implies a similar search space of 20P. In this combinatorial sense, then, the problems of directed protein evolution and of host engineering are broadly equivalent. In practice, however, they have different means for avoiding the inevitable difficulties of implementation. The spare capacity exhibited in metabolic networks implies that host engineering may admit substantial increases in flux to targets of interest. Thus, we rehearse the relevant issues for those wishing to understand and exploit those modern genome-wide host engineering tools and thinking that have been designed and developed to optimise fluxes towards desirable products in biotechnological processes, with a focus on microbial systems. The aim throughput is 'making such biology predictable'. Strategies have been aimed at both transcription and translation, especially for regulatory processes that can affect multiple targets. However, because there is a limit on how much protein a cell can produce, increasing kcat in selected targets may be a better strategy than increasing protein expression levels for optimal host engineering.
Collapse
|
6
|
Hala D, Faulkner P, He K, Kamalanathan M, Brink M, Simons K, Apaydin M, Hernout B, Petersen LH, Ivanov I, Qian X. An integrated in vivo and in silico analysis of the metabolism disrupting effects of CPI-613 on embryo-larval zebrafish (Danio rerio). Comp Biochem Physiol C Toxicol Pharmacol 2021; 248:109084. [PMID: 34051378 DOI: 10.1016/j.cbpc.2021.109084] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/09/2021] [Revised: 05/17/2021] [Accepted: 05/19/2021] [Indexed: 01/12/2023]
Abstract
CPI-613 is a mitochondrial metabolism disrupter that inhibits tricarboxylic acid (TCA) cycle activity. The consequences of TCA cycle disruption on various metabolic pathways and overall organismal physiology are not fully known. The present study integrates in vivo experimental data with an in silico stoichiometric metabolism model of zebrafish to study the metabolic pathways perturbed under CPI-613 exposure. Embryo-larval life stages of zebrafish (Danio rerio) were exposed to 1 μM CPI-613 for 20 days. Whole-organism respirometry measurements showed an initial suppression of O2 consumption at Day 5 of exposure, followed by recovery comparable to the solvent control (0.01% DMSO) by Day 20. Comparison of whole-transcriptome RNA-sequencing at Day 5 vs. 20 of exposure showed functional categories related to O2 binding and transport, antioxidant activity, FAD binding, and hemoglobin complexes, to be commonly represented. Metabolic enzyme gene expression changes and O2 consumption rate was used to parametrize two in silico stoichiometric metabolic models representative of Day 5 or 20 of exposure. Computational simulations predicted impaired ATP synthesis, α-ketoglutarate dehydrogenase (KGDH) activity, and fatty acid β-oxidation at Day 5 vs. 20 of exposure. These results show that the targeted disruption of KGDH may also impact oxidative phosphorylation (ATP synthesis) and fatty acid metabolism (β-oxidation), in turn influencing cellular bioenergetics and the observed reduction in whole-organism O2 consumption rate. The results of this study provide an integrated in vivo and in silico framework to study the impacts of metabolic disruption on organismal physiology.
Collapse
Affiliation(s)
- David Hala
- Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, USA; Department of Ecology and Conservation Biology, Texas A&M University, College Station, TX, USA.
| | - Patricia Faulkner
- Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, USA
| | - Kai He
- Department of Electrical & Computer Engineering, Texas A&M University, College Station, TX, USA
| | - Manoj Kamalanathan
- Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, USA
| | - Mikeelee Brink
- Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, USA
| | - Kristina Simons
- Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, USA
| | - Meltem Apaydin
- Department of Electrical & Computer Engineering, Texas A&M University, College Station, TX, USA
| | - Beatrice Hernout
- Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, USA; Institute for a Sustainable Environment, Department of Biology, Clarkson University, Potsdam, NY, USA
| | - Lene H Petersen
- Department of Marine Biology, Texas A&M University at Galveston, Galveston, TX, USA
| | - Ivan Ivanov
- Department of Veterinary Physiology & Pharmacology, Texas A&M University, College Station, TX, USA
| | - Xiaoning Qian
- Department of Electrical & Computer Engineering, Texas A&M University, College Station, TX, USA
| |
Collapse
|
7
|
Palsson BO. Genome‐Scale Models. Metab Eng 2021. [DOI: 10.1002/9783527823468.ch2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
8
|
In Silico Analysis of Functionalized Hydrocarbon Production Using Ehrlich Pathway and Fatty Acid Derivatives in an Endophytic Fungus. J Fungi (Basel) 2021; 7:jof7060435. [PMID: 34072611 PMCID: PMC8228540 DOI: 10.3390/jof7060435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Revised: 05/25/2021] [Accepted: 05/26/2021] [Indexed: 11/17/2022] Open
Abstract
Functionalized hydrocarbons have various ecological and industrial uses, from signaling molecules and antifungal/antibacterial agents to fuels and specialty chemicals. The potential to produce functionalized hydrocarbons using the cellulolytic, endophytic fungus, Ascocoryne sarcoides, was quantified using genome-enabled, stoichiometric modeling. In silico analysis identified available routes to produce these hydrocarbons, including both anabolic- and catabolic-associated strategies, and determined correlations between the type and size of the hydrocarbons and culturing conditions. The analysis quantified the limits of the wild-type metabolic network to produce functionalized hydrocarbons from cellulose-based substrates and identified metabolic engineering targets, including cellobiose phosphorylase (CP) and cytosolic pyruvate dehydrogenase complex (PDHcyt). CP and PDHcyt activity increased the theoretical production limits under anoxic conditions where less energy was extracted from the substrate. The incorporation of both engineering targets resulted in near-complete conservation of substrate electrons in functionalized hydrocarbons. The in silico framework was integrated with in vitro fungal batch growth experiments to support O2 limitation and functionalized hydrocarbon production predictions. The metabolic reconstruction of this endophytic filamentous fungus describes pathways for both specific and general production strategies of 161 functionalized hydrocarbons applicable to many eukaryotic hosts.
Collapse
|
9
|
Seep L, Razaghi-Moghadam Z, Nikoloski Z. Reaction lumping in metabolic networks for application with thermodynamic metabolic flux analysis. Sci Rep 2021; 11:8544. [PMID: 33879809 PMCID: PMC8058346 DOI: 10.1038/s41598-021-87643-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Accepted: 03/26/2021] [Indexed: 11/13/2022] Open
Abstract
Thermodynamic metabolic flux analysis (TMFA) can narrow down the space of steady-state flux distributions, but requires knowledge of the standard Gibbs free energy for the modelled reactions. The latter are often not available due to unknown Gibbs free energy change of formation [Formula: see text], of metabolites. To optimize the usage of data on thermodynamics in constraining a model, reaction lumping has been proposed to eliminate metabolites with unknown [Formula: see text]. However, the lumping procedure has not been formalized nor implemented for systematic identification of lumped reactions. Here, we propose, implement, and test a combined procedure for reaction lumping, applicable to genome-scale metabolic models. It is based on identification of groups of metabolites with unknown [Formula: see text] whose elimination can be conducted independently of the others via: (1) group implementation, aiming to eliminate an entire such group, and, if this is infeasible, (2) a sequential implementation to ensure that a maximal number of metabolites with unknown [Formula: see text] are eliminated. Our comparative analysis with genome-scale metabolic models of Escherichia coli, Bacillus subtilis, and Homo sapiens shows that the combined procedure provides an efficient means for systematic identification of lumped reactions. We also demonstrate that TMFA applied to models with reactions lumped according to the proposed procedure lead to more precise predictions in comparison to the original models. The provided implementation thus ensures the reproducibility of the findings and their application with standard TMFA.
Collapse
Affiliation(s)
- Lea Seep
- Bioinformatics, Institute for Biochemistry and Biology, University of Potsdam, 14476, Potsdam, Germany
| | - Zahra Razaghi-Moghadam
- Bioinformatics, Institute for Biochemistry and Biology, University of Potsdam, 14476, Potsdam, Germany
- Systems Biology and Mathematical Modeling, Max Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany
| | - Zoran Nikoloski
- Bioinformatics, Institute for Biochemistry and Biology, University of Potsdam, 14476, Potsdam, Germany.
- Systems Biology and Mathematical Modeling, Max Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany.
| |
Collapse
|
10
|
Cheng Y, Schlosser P, Hertel J, Sekula P, Oefner PJ, Spiekerkoetter U, Mielke J, Freitag DF, Schmidts M, Kronenberg F, Eckardt KU, Thiele I, Li Y, Köttgen A. Rare genetic variants affecting urine metabolite levels link population variation to inborn errors of metabolism. Nat Commun 2021; 12:964. [PMID: 33574263 PMCID: PMC7878905 DOI: 10.1038/s41467-020-20877-8] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2020] [Accepted: 12/21/2020] [Indexed: 02/07/2023] Open
Abstract
Metabolite levels in urine may provide insights into genetic mechanisms shaping their related pathways. We therefore investigate the cumulative contribution of rare, exonic genetic variants on urine levels of 1487 metabolites and 53,714 metabolite ratios among 4864 GCKD study participants. Here we report the detection of 128 significant associations involving 30 unique genes, 16 of which are known to underlie inborn errors of metabolism. The 30 genes are strongly enriched for shared expression in liver and kidney (odds ratio = 65, p-FDR = 3e-7), with hepatocytes and proximal tubule cells as driving cell types. Use of UK Biobank whole-exome sequencing data links genes to diseases connected to the identified metabolites. In silico constraint-based modeling of gene knockouts in a virtual whole-body, organ-resolved metabolic human correctly predicts the observed direction of metabolite changes, highlighting the potential of linking population genetics to modeling. Our study implicates candidate variants and genes for inborn errors of metabolism.
Collapse
Affiliation(s)
- Yurong Cheng
- Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany
- Faculty of Biology, University of Freiburg, Freiburg, Germany
| | - Pascal Schlosser
- Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany
| | - Johannes Hertel
- School of Medicine, National University of Ireland, Galway, University Road, Galway, Ireland
- University of Greifswald, University Medicine Greifswald, Department of Psychiatry and Psychotherapy, Greifswald, Germany
| | - Peggy Sekula
- Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany
| | - Peter J Oefner
- Institute of Functional Genomics, University of Regensburg, Regensburg, Germany
| | - Ute Spiekerkoetter
- Department of General Pediatrics and Adolescent Medicine, Medical Center and Faculty of Medicine - University of Freiburg, Freiburg, Germany
| | - Johanna Mielke
- Bayer AG, Division Pharmaceuticals, Open Innovation & Digital Technologies, Wuppertal, Germany
| | - Daniel F Freitag
- Bayer AG, Division Pharmaceuticals, Open Innovation & Digital Technologies, Wuppertal, Germany
| | - Miriam Schmidts
- Department of General Pediatrics and Adolescent Medicine, Medical Center and Faculty of Medicine - University of Freiburg, Freiburg, Germany
| | - Florian Kronenberg
- Institute of Genetic Epidemiology, Department of Genetics and Pharmacology, Medical University of Innsbruck, Innsbruck, Austria
| | - Kai-Uwe Eckardt
- Department of Nephrology and Hypertension, University of Erlangen-Nürnberg, Erlangen, Germany
- Department of Nephrology and Medical Intensive Care, Charité - Universitätsmedizin Berlin, Berlin, Germany
| | - Ines Thiele
- School of Medicine, National University of Ireland, Galway, University Road, Galway, Ireland
- Division of Microbiology, National University of Ireland, Galway, University Road, Galway, Ireland
- APC Microbiome Ireland, Galway, Ireland
| | - Yong Li
- Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany
| | - Anna Köttgen
- Institute of Genetic Epidemiology, Faculty of Medicine and Medical Center - University of Freiburg, Freiburg, Germany.
- CIBSS - Centre for Integrative Biological Signalling Studies, Albert-Ludwigs-Universität Freiburg, Freiburg, Germany.
| |
Collapse
|
11
|
Ortega-Quintana FA, Trujillo-Roldán MA, Botero-Castro H, Alvarez H. Modeling the interaction between the central carbon metabolism of Escherichia coli and bioreactor culture media. Biochem Eng J 2020. [DOI: 10.1016/j.bej.2020.107753] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
12
|
Macklin DN, Ahn-Horst TA, Choi H, Ruggero NA, Carrera J, Mason JC, Sun G, Agmon E, DeFelice MM, Maayan I, Lane K, Spangler RK, Gillies TE, Paull ML, Akhter S, Bray SR, Weaver DS, Keseler IM, Karp PD, Morrison JH, Covert MW. Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation. Science 2020; 369:eaav3751. [PMID: 32703847 PMCID: PMC7990026 DOI: 10.1126/science.aav3751] [Citation(s) in RCA: 76] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2018] [Revised: 10/28/2019] [Accepted: 05/26/2020] [Indexed: 12/24/2022]
Abstract
The extensive heterogeneity of biological data poses challenges to analysis and interpretation. Construction of a large-scale mechanistic model of Escherichia coli enabled us to integrate and cross-evaluate a massive, heterogeneous dataset based on measurements reported by various groups over decades. We identified inconsistencies with functional consequences across the data, including that the total output of the ribosomes and RNA polymerases described by data are not sufficient for a cell to reproduce measured doubling times, that measured metabolic parameters are neither fully compatible with each other nor with overall growth, and that essential proteins are absent during the cell cycle-and the cell is robust to this absence. Finally, considering these data as a whole leads to successful predictions of new experimental outcomes, in this case protein half-lives.
Collapse
Affiliation(s)
- Derek N Macklin
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Travis A Ahn-Horst
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Heejo Choi
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Nicholas A Ruggero
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
- Department of Chemical Engineering, Stanford University, Stanford, CA 94305, USA
| | - Javier Carrera
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - John C Mason
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Gwanggyu Sun
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Eran Agmon
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Mialy M DeFelice
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Inbal Maayan
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Keara Lane
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Ryan K Spangler
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Taryn E Gillies
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Morgan L Paull
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | - Sajia Akhter
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | - Samuel R Bray
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA
| | | | | | | | - Jerry H Morrison
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| | - Markus W Covert
- Department of Bioengineering, Stanford University, Stanford, CA 94305, USA.
- Allen Discovery Center at Stanford University, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
13
|
Sigmarsdóttir Þ, McGarrity S, Rolfsson Ó, Yurkovich JT, Sigurjónsson ÓE. Current Status and Future Prospects of Genome-Scale Metabolic Modeling to Optimize the Use of Mesenchymal Stem Cells in Regenerative Medicine. Front Bioeng Biotechnol 2020; 8:239. [PMID: 32296688 PMCID: PMC7136564 DOI: 10.3389/fbioe.2020.00239] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2019] [Accepted: 03/09/2020] [Indexed: 12/15/2022] Open
Abstract
Mesenchymal stem cells are a promising source for externally grown tissue replacements and patient-specific immunomodulatory treatments. This promise has not yet been fulfilled in part due to production scaling issues and the need to maintain the correct phenotype after re-implantation. One aspect of extracorporeal growth that may be manipulated to optimize cell growth and differentiation is metabolism. The metabolism of MSCs changes during and in response to differentiation and immunomodulatory changes. MSC metabolism may be linked to functional differences but how this occurs and influences MSC function remains unclear. Understanding how MSC metabolism relates to cell function is however important as metabolite availability and environmental circumstances in the body may affect the success of implantation. Genome-scale constraint based metabolic modeling can be used as a tool to fill gaps in knowledge of MSC metabolism, acting as a framework to integrate and understand various data types (e.g., genomic, transcriptomic and metabolomic). These approaches have long been used to optimize the growth and productivity of bacterial production systems and are being increasingly used to provide insights into human health research. Production of tissue for implantation using MSCs requires both optimized production of cell mass and the understanding of the patient and phenotype specific metabolic situation. This review considers the current knowledge of MSC metabolism and how it may be optimized along with the current and future uses of genome scale constraint based metabolic modeling to further this aim.
Collapse
Affiliation(s)
- Þóra Sigmarsdóttir
- The Blood Bank, Landspitali – The National University Hospital of Iceland, Reykjavik, Iceland
- School of Science and Engineering, Reykjavik University, Reykjavik, Iceland
| | - Sarah McGarrity
- School of Science and Engineering, Reykjavik University, Reykjavik, Iceland
- Faculty of Medicine, School of Health Sciences, University of Iceland, Reykjavik, Iceland
| | - Óttar Rolfsson
- Faculty of Medicine, School of Health Sciences, University of Iceland, Reykjavik, Iceland
| | | | - Ólafur E. Sigurjónsson
- The Blood Bank, Landspitali – The National University Hospital of Iceland, Reykjavik, Iceland
- School of Science and Engineering, Reykjavik University, Reykjavik, Iceland
| |
Collapse
|
14
|
Yang X, Yuan Q, Luo H, Li F, Mao Y, Zhao X, Du J, Li P, Ju X, Zheng Y, Chen Y, Liu Y, Jiang H, Yao Y, Ma H, Ma Y. Systematic design and in vitro validation of novel one-carbon assimilation pathways. Metab Eng 2019; 56:142-153. [DOI: 10.1016/j.ymben.2019.09.001] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2019] [Revised: 08/17/2019] [Accepted: 09/01/2019] [Indexed: 11/24/2022]
|
15
|
Ghasemi-Kahrizsangi T, Marashi SA, Hosseini Z. Genome-Scale Metabolic Network Models of Bacillus Species Suggest that Model Improvement is Necessary for Biotechnological Applications. IRANIAN JOURNAL OF BIOTECHNOLOGY 2019; 16:e1684. [PMID: 31457023 PMCID: PMC6697824 DOI: 10.15171/ijb.1684] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Revised: 09/07/2017] [Accepted: 09/18/2017] [Indexed: 11/11/2022]
Abstract
Background A genome-scale metabolic network model (GEM) is a mathematical representation of an organism’s metabolism. Today, GEMs are popular tools for computationally simulating the biotechnological processes and for predicting biochemical properties of (engineered) strains. Objectives In the present study, we have evaluated the predictive power of two GEMs, namely iBsu1103 (for Bacillus subtilis 168) and iMZ1055 (for Bacillus megaterium WSH002). Materials and Methods For comparing the predictive power of Bacillus subtilis and Bacillus megaterium GEMs, experimental data were obtained from previous wet-lab studies included in PubMed. By using these data, we set the environmental, stoichiometric and thermodynamic constraints on the models, and FBA is performed to predict the biomass production rate, and the values of other fluxes. For simulating experimental conditions in this study, COBRA toolbox was used. Results By using the wealth of data in the literature, we evaluated the accuracy of in silico simulations of these GEMs. Our results suggest that there are some errors in these two models which make them unreliable for predicting the biochemical capabilities of these species. The inconsistencies between experimental and computational data are even greater where B. subtilis and B. megaterium do not have similar phenotypes. Conclusions Our analysis suggests that literature-based improvement of genome-scale metabolic network models of the two Bacillus species is essential if these models are to be successfully applied in biotechnology and metabolic engineering.
Collapse
Affiliation(s)
| | - Sayed-Amir Marashi
- Department of Biotechnology, College of Science, University of Tehran, Tehran, Iran
| | - Zhaleh Hosseini
- Department of Biotechnology, College of Science, University of Tehran, Tehran, Iran
| |
Collapse
|
16
|
Kim EY, Ashlock D, Yoon SH. Identification of critical connectors in the directed reaction-centric graphs of microbial metabolic networks. BMC Bioinformatics 2019; 20:328. [PMID: 31195955 PMCID: PMC6567475 DOI: 10.1186/s12859-019-2897-z] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2019] [Accepted: 05/13/2019] [Indexed: 12/20/2022] Open
Abstract
BACKGROUND Detection of central nodes in asymmetrically directed biological networks depends on centrality metrics quantifying individual nodes' importance in a network. In topological analyses on metabolic networks, various centrality metrics have been mostly applied to metabolite-centric graphs. However, centrality metrics including those not depending on high connections are largely unexplored for directed reaction-centric graphs. RESULTS We applied directed versions of centrality metrics to directed reaction-centric graphs of microbial metabolic networks. To investigate the local role of a node, we developed a novel metric, cascade number, considering how many nodes are closed off from information flow when a particular node is removed. High modularity and scale-freeness were found in the directed reaction-centric graphs and betweenness centrality tended to belong to densely connected modules. Cascade number and bridging centrality identified cascade subnetworks controlling local information flow and irreplaceable bridging nodes between functional modules, respectively. Reactions highly ranked with bridging centrality and cascade number tended to be essential, compared to reactions that other central metrics detected. CONCLUSIONS We demonstrate that cascade number and bridging centrality are useful to identify key reactions controlling local information flow in directed reaction-centric graphs of microbial metabolic networks. Knowledge about the local flow connectivity and connections between local modules will contribute to understand how metabolic pathways are assembled.
Collapse
Affiliation(s)
- Eun-Youn Kim
- School of Basic Sciences, Hanbat National University, Daejeon, 34158, Republic of Korea
| | - Daniel Ashlock
- Department of Mathematics and Statistics, the University of Guelph, Guelph, Ontario, N1G 2W1, Canada
| | - Sung Ho Yoon
- Department of Bioscience and Biotechnology, Konkuk University, Seoul, 05029, Republic of Korea.
| |
Collapse
|
17
|
Cheng C, O’Brien EJ, McCloskey D, Utrilla J, Olson C, LaCroix RA, Sandberg TE, Feist AM, Palsson BO, King ZA. Laboratory evolution reveals a two-dimensional rate-yield tradeoff in microbial metabolism. PLoS Comput Biol 2019; 15:e1007066. [PMID: 31158228 PMCID: PMC6564042 DOI: 10.1371/journal.pcbi.1007066] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2018] [Revised: 06/13/2019] [Accepted: 05/02/2019] [Indexed: 01/12/2023] Open
Abstract
Growth rate and yield are fundamental features of microbial growth. However, we lack a mechanistic and quantitative understanding of the rate-yield relationship. Studies pairing computational predictions with experiments have shown the importance of maintenance energy and proteome allocation in explaining rate-yield tradeoffs and overflow metabolism. Recently, adaptive evolution experiments of Escherichia coli reveal a phenotypic diversity beyond what has been explained using simple models of growth rate versus yield. Here, we identify a two-dimensional rate-yield tradeoff in adapted E. coli strains where the dimensions are (A) a tradeoff between growth rate and yield and (B) a tradeoff between substrate (glucose) uptake rate and growth yield. We employ a multi-scale modeling approach, combining a previously reported coarse-grained small-scale proteome allocation model with a fine-grained genome-scale model of metabolism and gene expression (ME-model), to develop a quantitative description of the full rate-yield relationship for E. coli K-12 MG1655. The multi-scale analysis resolves the complexity of ME-model which hindered its practical use in proteome complexity analysis, and provides a mechanistic explanation of the two-dimensional tradeoff. Further, the analysis identifies modifications to the P/O ratio and the flux allocation between glycolysis and pentose phosphate pathway (PPP) as potential mechanisms that enable the tradeoff between glucose uptake rate and growth yield. Thus, the rate-yield tradeoffs that govern microbial adaptation to new environments are more complex than previously reported, and they can be understood in mechanistic detail using a multi-scale modeling approach. This study reconciles multiple existing microbial rate-yield tradeoff theories with experimental data. There is great interest in developing quantitative descriptions of the relationship between growth rate and growth yield [1]. However, some reported experiments [2–4] in the literature do not agree with existing theories [5–7]. Specifically, overflow metabolism in E. coli can either be coupled [5, 8] or decoupled [2–4] from growth rate. We found that adaptive laboratory evolution (ALE) experiments of E. coli reveal a two-dimensional rate-yield tradeoff in adapted strains where the dimensions are (i) a tradeoff between growth rate and growth yield, previously reported by [5], and (ii) a tradeoff between substrate uptake rate and growth yield. The appearance of this two-dimensional tradeoff during adaptation suggests that microorganisms adapting to new environments are subject to a more complex set of rate-yield tradeoffs than previously reported [5, 6]. In this study, the two-dimensional rate-yield tradeoff is quantitatively explained through our multi-scale modeling approach, combining a previously reported small-scale proteome allocation model [5] with a genome-scale model of metabolism and gene-expression (ME-model) [9]. The modeling approach is also instrumental to future studies.
Collapse
Affiliation(s)
- Chuankai Cheng
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Edward J. O’Brien
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Douglas McCloskey
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Jose Utrilla
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Connor Olson
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Ryan A. LaCroix
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Troy E. Sandberg
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
| | - Adam M. Feist
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Lyngby, Denmark
| | - Bernhard O. Palsson
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Lyngby, Denmark
- Department of Pediatrics, Univerity of California San Diego, La Jolla, California, United States of America
| | - Zachary A. King
- Department of Bioengineering, University of California San Diego, La Jolla, California, United States of America
- * E-mail:
| |
Collapse
|
18
|
Ajjolli Nagaraja A, Fontaine N, Delsaut M, Charton P, Damour C, Offmann B, Grondin-Perez B, Cadet F. Flux prediction using artificial neural network (ANN) for the upper part of glycolysis. PLoS One 2019; 14:e0216178. [PMID: 31067238 PMCID: PMC6505829 DOI: 10.1371/journal.pone.0216178] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2019] [Accepted: 04/15/2019] [Indexed: 01/08/2023] Open
Abstract
The selection of optimal enzyme concentration in multienzyme cascade reactions for the highest product yield in practice is very expensive and time-consuming process. The modelling of biological pathways is a difficult process because of the complexity of the system. The mathematical modelling of the system using an analytical approach depends on the many parameters of enzymes which rely on tedious and expensive experiments. The artificial neural network (ANN) method has been successively applied in different fields of science to perform complex functions. In this study, ANN models were trained to predict the flux for the upper part of glycolysis as inferred by NADH consumption, using four enzyme concentrations i.e., phosphoglucoisomerase, phosphofructokinase, fructose-bisphosphate-aldolase, triose-phosphate-isomerase. Out of three ANN algorithms, the neuralnet package with two activation functions, “logistic” and “tanh” were implemented. The prediction of the flux was very efficient: RMSE and R2 were 0.847, 0.93 and 0.804, 0.94 respectively for logistic and tanh functions using a cross validation procedure. This study showed that a systemic approach such as ANN could be used for accurate prediction of the flux through the metabolic pathway. This could help to save a lot of time and costs, particularly from an industrial perspective. The R-code is available at: https://github.com/DSIMB/ANN-Glycolysis-Flux-Prediction.
Collapse
Affiliation(s)
- Anamya Ajjolli Nagaraja
- LE2P, Laboratory of Energy, Electronics and Processes EA 4079, Faculty of Sciences and Technology, University of La Reunion, France
| | | | - Mathieu Delsaut
- LE2P, Laboratory of Energy, Electronics and Processes EA 4079, Faculty of Sciences and Technology, University of La Reunion, France
| | - Philippe Charton
- DSIMB, INSERM, UMR S-1134, Laboratory of ExcellenceLABEX GR, Faculty of Sciences and Technology, University of La Reunion & University Paris Diderot, Paris, France
| | - Cedric Damour
- LE2P, Laboratory of Energy, Electronics and Processes EA 4079, Faculty of Sciences and Technology, University of La Reunion, France
| | - Bernard Offmann
- Université de Nantes, Unité Fonctionnalité et Ingénierie des Protéines (UFIP), UMR 6286 CNRS, UFR Sciences et Techniques, chemin de la Houssinière, France
| | - Brigitte Grondin-Perez
- LE2P, Laboratory of Energy, Electronics and Processes EA 4079, Faculty of Sciences and Technology, University of La Reunion, France
| | - Frederic Cadet
- DSIMB, INSERM, UMR S-1134, Laboratory of ExcellenceLABEX GR, Faculty of Sciences and Technology, University of La Reunion & University Paris Diderot, Paris, France
- * E-mail:
| |
Collapse
|
19
|
Metagenome level metabolic network reconstruction analysis reveals the microbiome in the Bogotá River is functionally close to the microbiome in produced water. Ecol Modell 2019. [DOI: 10.1016/j.ecolmodel.2019.02.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
20
|
Yaşar Yildiz S, Nikerel E, Toksoy Öner E. Genome-Scale Metabolic Model of a Microbial Cell Factory ( Brevibacillus thermoruber 423) with Multi-Industry Potentials for Exopolysaccharide Production. OMICS-A JOURNAL OF INTEGRATIVE BIOLOGY 2019; 23:237-246. [PMID: 30932743 DOI: 10.1089/omi.2019.0028] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Brevibacillus thermoruber 423 is a thermophilic bacterium capable of producing high levels of exopolysaccharide (EPS) that has broad applications in nutrition, feed, cosmetics, pharmaceutical, and chemical industries, not to mention in health and bionanotechnology sectors. EPS is a natural, nontoxic, and biodegradable polymer of sugar residues and plays pivotal roles in cell-to-cell interactions, adhesion, biofilm formation, and protection of cell against environmental extremes. This bacterium is a thermophilic EPS producer while exceeding other thermophilic producers by virtue of high level of polymer synthesis. Recently, B. thermoruber 423 was noted for relevance to multiple industry sectors because of its capacity to use xylose, and produce EPS, isoprenoids, ethanol/butanol, lipases, proteases, cellulase, and glucoamylase enzymes as well as its resistance to arsenic. A key step in understanding EPS production with a systems-based approach is the knowledge of microbial genome sequence. To speed biotechnology and industrial applications, this study reports on a genome-scale metabolic model (GSMM) of B. thermoruber 423, constructed using the recently available high-quality genome sequence that we have subsequently validated using physiological data on batch growth and EPS production on seven different carbon sources. The model developed contains 1454 reactions (of which 1127 are assigned an enzyme commission number) and 1410 metabolites from 925 genes. This GSMM offers the promise to enable and accelerate further systems biology and industrial scale studies, not to mention the ability to calculate metabolic flux distribution in large networks and multiomic data integration.
Collapse
Affiliation(s)
- Songül Yaşar Yildiz
- 1 Department of Bioengineering, Istanbul Medeniyet University, Istanbul, Turkey
| | - Emrah Nikerel
- 2 Department of Genetics and Bioengineering, Yeditepe University, Istanbul, Turkey
| | - Ebru Toksoy Öner
- 3 Department of Bioengineering, IBSB, Marmara University, Istanbul, Turkey
| |
Collapse
|
21
|
Kim H, Kim S, Yoon SH. Metabolic network reconstruction and phenome analysis of the industrial microbe, Escherichia coli BL21(DE3). PLoS One 2018; 13:e0204375. [PMID: 30240424 PMCID: PMC6150544 DOI: 10.1371/journal.pone.0204375] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2018] [Accepted: 09/06/2018] [Indexed: 01/25/2023] Open
Abstract
Escherichia coli BL21(DE3) is an industrial model microbe for the mass-production of bioproducts such as biofuels, biorefineries, and recombinant proteins. However, despite its important role in scientific research and biotechnological applications, a high-quality metabolic network model for metabolic engineering is yet to be developed. Here, we present the comprehensive metabolic network model of E. coli BL21(DE3), named iHK1487, based on the latest genome reannotation and phenome analysis. The metabolic model consists of 1,164 unique metabolites, 2,701 metabolic reactions, and 1,487 genes. The model was validated and improved by comparing the simulation results with phenome data from phenotype microarray tests. Previous transcriptome profile data was incorporated during model reconstruction, and flux prediction was simulated using the model. iHK1487 was simulated to explore the metabolic features of BL21(DE3) such as broad spectrum amino acid utilization and enhanced flux through the upper glycolytic pathway and TCA cycle. iHK1487 will contribute to systematic understanding of cellular physiology and metabolism of E. coli BL21(DE3) and highlight its biotechnological applications.
Collapse
Affiliation(s)
- Hanseol Kim
- Department of Bioscience and Biotechnology, Konkuk University, Seoul, Republic of Korea
| | - Sinyeon Kim
- Department of Bioscience and Biotechnology, Konkuk University, Seoul, Republic of Korea
| | - Sung Ho Yoon
- Department of Bioscience and Biotechnology, Konkuk University, Seoul, Republic of Korea
| |
Collapse
|
22
|
Kumar M, Ji B, Babaei P, Das P, Lappa D, Ramakrishnan G, Fox TE, Haque R, Petri WA, Bäckhed F, Nielsen J. Gut microbiota dysbiosis is associated with malnutrition and reduced plasma amino acid levels: Lessons from genome-scale metabolic modeling. Metab Eng 2018; 49:128-142. [PMID: 30075203 PMCID: PMC6871511 DOI: 10.1016/j.ymben.2018.07.018] [Citation(s) in RCA: 50] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2018] [Revised: 07/30/2018] [Accepted: 07/30/2018] [Indexed: 12/31/2022]
Abstract
Malnutrition is a severe non-communicable disease, which is prevalent in children from low-income countries. Recently, a number of metagenomics studies have illustrated associations between the altered gut microbiota and child malnutrition. However, these studies did not examine metabolic functions and interactions between individual species in the gut microbiota during health and malnutrition. Here, we applied genome-scale metabolic modeling to model the gut microbial species, which were selected from healthy and malnourished children from three countries. Our analysis showed reduced metabolite production capabilities in children from two low-income countries compared with a high-income country. Additionally, the models were also used to predict the community-level metabolic potentials of gut microbes and the patterns of pairwise interactions among species. Hereby we found that due to bacterial interactions there may be reduced production of certain amino acids in malnourished children compared with healthy children from the same communities. To gain insight into alterations in the metabolism of malnourished (stunted) children, we also performed targeted plasma metabolic profiling in the first 2 years of life of 25 healthy and 25 stunted children. Plasma metabolic profiling further revealed that stunted children had reduced plasma levels of essential amino acids compared to healthy controls. Our analyses provide a framework for future efforts towards further characterization of gut microbial metabolic capabilities and their contribution to malnutrition.
Collapse
Affiliation(s)
- Manish Kumar
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE41128 Gothenburg, Sweden
| | - Boyang Ji
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE41128 Gothenburg, Sweden
| | - Parizad Babaei
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE41128 Gothenburg, Sweden
| | - Promi Das
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE41128 Gothenburg, Sweden
| | - Dimitra Lappa
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE41128 Gothenburg, Sweden
| | - Girija Ramakrishnan
- Department of Medicine/Division of Infectious Diseases, and University of Virginia, Charlottesville, VA, USA
| | - Todd E Fox
- Department of Pharmacology, University of Virginia, Charlottesville, VA, USA
| | - Rashidul Haque
- International Centre for Diarrheal Disease Research, Dhaka, Bangladesh
| | - William A Petri
- Department of Medicine/Division of Infectious Diseases, and University of Virginia, Charlottesville, VA, USA
| | - Fredrik Bäckhed
- The Wallenberg Laboratory, Department of Molecular and Clinical Medicine, University of Gothenburg, 41345 Gothenburg, Sweden; Novo Nordisk Foundation Center for Basic Metabolic Research, Section for Metabolic Receptology and Enteroendocrinology, Faculty of Health Sciences, University of Copenhagen, 2200 Copenhagen, Denmark
| | - Jens Nielsen
- Department of Biology and Biological Engineering, Chalmers University of Technology, SE41128 Gothenburg, Sweden; Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, DK2800 Lyngby, Denmark.
| |
Collapse
|
23
|
Tajparast M, Frigon D. Predicting the accumulation of storage compounds by Rhodococcus jostii RHA1 in the feast-famine growth cycles using genome-scale flux balance analysis. PLoS One 2018; 13:e0191835. [PMID: 29494607 PMCID: PMC5832212 DOI: 10.1371/journal.pone.0191835] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2017] [Accepted: 01/11/2018] [Indexed: 01/18/2023] Open
Abstract
Feast-famine cycles in biological wastewater resource recovery systems select for bacterial species that accumulate intracellular storage compounds such as poly-β-hydroxybutyrate (PHB), glycogen, and triacylglycerols (TAG). These species survive better the famine phase and resume rapid substrate uptake at the beginning of the feast phase faster than microorganisms unable to accumulate storage. However, ecophysiological conditions favouring the accumulation of either storage compounds remain to be clarified, and predictive capabilities need to be developed to eventually rationally design reactors producing these compounds. Using a genome-scale metabolic modelling approach, the storage metabolism of Rhodococcus jostii RHA1 was investigated for steady-state feast-famine cycles on glucose and acetate as the sole carbon sources. R. jostii RHA1 is capable of accumulating the three storage compounds (PHB, TAG, and glycogen) simultaneously. According to the experimental observations, when glucose was the substrate, feast phase chemical oxygen demand (COD) accumulation was similar for the three storage compounds; when acetate was the substrate, however, PHB accumulation was 3 times higher than TAG accumulation and essentially no glycogen was accumulated. These results were simulated using the genome-scale metabolic model of R. jostii RHA1 (iMT1174) by means of flux balance analysis (FBA) to determine the objective functions capable of predicting these behaviours. Maximization of the growth rate was set as the main objective function, while minimization of total reaction fluxes and minimization of metabolic adjustment (environmental MOMA) were considered as the sub-objective functions. The environmental MOMA sub-objective performed better than the minimization of total reaction fluxes sub-objective function at predicting the mixture of storage compounds accumulated. Additional experiments with 13C-labelled bicarbonate (HCO3−) found that the fluxes through the central metabolism reactions during the feast phases were similar to the ones during the famine phases on acetate due to similarity in the carbon sources in the feast and famine phases (i.e., acetate and poly-β-hydroxybutyrate, respectively); this suggests that the environmental MOMA sub-objective function could be used to analyze successive environmental conditions such as the feast and famine cycles while the metabolically similar carbon sources are taken up by microorganisms.
Collapse
Affiliation(s)
- Mohammad Tajparast
- Microbial Community Engineering Laboratory, Department of Civil Engineering and Applied Mechanics, McGill University, Montreal, Quebec, Canada
| | - Dominic Frigon
- Microbial Community Engineering Laboratory, Department of Civil Engineering and Applied Mechanics, McGill University, Montreal, Quebec, Canada
- * E-mail:
| |
Collapse
|
24
|
Soubeyrand E, Colombié S, Beauvoit B, Dai Z, Cluzet S, Hilbert G, Renaud C, Maneta-Peyret L, Dieuaide-Noubhani M, Mérillon JM, Gibon Y, Delrot S, Gomès E. Constraint-Based Modeling Highlights Cell Energy, Redox Status and α-Ketoglutarate Availability as Metabolic Drivers for Anthocyanin Accumulation in Grape Cells Under Nitrogen Limitation. FRONTIERS IN PLANT SCIENCE 2018; 9:421. [PMID: 29868039 PMCID: PMC5966944 DOI: 10.3389/fpls.2018.00421] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/07/2018] [Accepted: 03/16/2018] [Indexed: 05/18/2023]
Abstract
Anthocyanin biosynthesis is regulated by environmental factors (such as light, temperature, and water availability) and nutrient status (such as carbon, nitrogen, and phosphate nutrition). Previous reports show that low nitrogen availability strongly enhances anthocyanin accumulation in non carbon-limited plant organs or cell suspensions. It has been hypothesized that high carbon-to-nitrogen ratio would lead to an energy excess in plant cells, and that an increase in flavonoid pathway metabolic fluxes would act as an "energy escape valve," helping plant cells to cope with energy and carbon excess. However, this hypothesis has never been tested directly. To this end, we used the grapevine Vitis vinifera L. cultivar Gamay Teinturier (syn. Gamay Freaux or Freaux Tintorier, VIVC #4382) cell suspension line as a model system to study the regulation of anthocyanin accumulation in response to nitrogen supply. The cells were sub-cultured in the presence of either control (25 mM) or low (5 mM) nitrate concentration. Targeted metabolomics and enzyme activity determinations were used to parametrize a constraint-based model describing both the central carbon and nitrogen metabolisms and the flavonoid (phenylpropanoid) pathway connected by the energy (ATP) and reducing power equivalents (NADPH and NADH) cofactors. The flux analysis (2 flux maps generated, for control and low nitrogen in culture medium) clearly showed that in low nitrogen-fed cells all the metabolic fluxes of central metabolism were decreased, whereas fluxes that consume energy and reducing power, were either increased (upper part of glycolysis, shikimate, and flavonoid pathway) or maintained (pentose phosphate pathway). Also, fluxes of flavanone 3β-hydroxylase, flavonol synthase, and anthocyanidin synthase were strongly increased, advocating for a regulation of the flavonoid pathway by alpha-ketoglutarate levels. These results strongly support the hypothesis of anthocyanin biosynthesis acting as an energy escape valve in plant cells, and they open new possibilities to manipulate flavonoid production in plant cells. They do not, however, support a role of anthocyanins as an effective mechanism for coping with carbon excess in high carbon to nitrogen ratio situations in grape cells. Instead, constraint-based modeling output and biomass analysis indicate that carbon excess is dealt with by vacuolar storage of soluble sugars.
Collapse
Affiliation(s)
- Eric Soubeyrand
- UMR 1287 Ecophysiologie et Génomique Fonctionnelle de la Vigne, Université de Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
| | - Sophie Colombié
- UMR 1332 Biologie du Fruit et Pathologie, INRA-Bordeaux, IBVM, Bordeaux, France
| | - Bertrand Beauvoit
- UMR 1332 Biologie du Fruit et Pathologie, INRA-Bordeaux, IBVM, Bordeaux, France
| | - Zhanwu Dai
- UMR 1287 Ecophysiologie et Génomique Fonctionnelle de la Vigne, INRA-Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
| | - Stéphanie Cluzet
- EA 3675 GESVAB, Université de Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
| | - Ghislaine Hilbert
- UMR 1287 Ecophysiologie et Génomique Fonctionnelle de la Vigne, INRA-Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
| | - Christel Renaud
- UMR 1287 Ecophysiologie et Génomique Fonctionnelle de la Vigne, INRA-Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
| | - Lilly Maneta-Peyret
- UMR 5200 Laboratoire de Biogenèse Membranaire, Université de Bordeaux, Bordeaux, France
| | | | - Jean-Michel Mérillon
- EA 3675 GESVAB, Université de Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
| | - Yves Gibon
- UMR 1332 Biologie du Fruit et Pathologie, INRA-Bordeaux, IBVM, Bordeaux, France
| | - Serge Delrot
- UMR 1287 Ecophysiologie et Génomique Fonctionnelle de la Vigne, Université de Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
| | - Eric Gomès
- UMR 1287 Ecophysiologie et Génomique Fonctionnelle de la Vigne, Université de Bordeaux, Institut des Sciences de la Vigne et du Vin, Bordeaux, France
- *Correspondence: Eric Gomès,
| |
Collapse
|
25
|
Fyson N, King J, Belcher T, Preston A, Colijn C. A curated genome-scale metabolic model of Bordetella pertussis metabolism. PLoS Comput Biol 2017; 13:e1005639. [PMID: 28715411 PMCID: PMC5553986 DOI: 10.1371/journal.pcbi.1005639] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2016] [Revised: 08/11/2017] [Accepted: 06/15/2017] [Indexed: 11/30/2022] Open
Abstract
The Gram-negative bacterium Bordetella pertussis is the causative agent of whooping cough, a serious respiratory infection causing hundreds of thousands of deaths annually worldwide. There are effective vaccines, but their production requires growing large quantities of B. pertussis. Unfortunately, B. pertussis has relatively slow growth in culture, with low biomass yields and variable growth characteristics. B. pertussis also requires a relatively expensive growth medium. We present a new, curated flux balance analysis-based model of B. pertussis metabolism. We enhance the model with an experimentally-determined biomass objective function, and we perform extensive manual curation. We test the model’s predictions with a genome-wide screen for essential genes using a transposon-directed insertional sequencing (TraDIS) approach. We test its predictions of growth for different carbon sources in the medium. The model predicts essentiality with an accuracy of 83% and correctly predicts improvements in growth under increased glutamate:fumarate ratios. We provide the model in SBML format, along with gene essentiality predictions. Metabolic flux models have been used to understand how organisms adapt their metabolism under different growth conditions, and are finding increasing application in synthetic biology and biotechnology. One barrier to progress in this field is the construction and curation of metabolic flux models for new organisms. Here we present a curated genome-scale metabolic flux model for Bordetella pertussis, the causative agent of whooping cough. Producing vaccines against whooping cough requires growing B. pertussis in large volumes. However, its growth is relatively slow, final yields of biomass are relatively low and growth characteristics can be variable. Understanding B. pertussis metabolism has applications to improving vaccine production, as well as in understanding the basic biology of this organism.
Collapse
Affiliation(s)
- Nick Fyson
- Department of Mathematics, Imperial College, London, UK
| | - Jerry King
- The Milner Centre for Evolution and Department of Biology and Biochemistry, University of Bath, Bath, UK
| | - Thomas Belcher
- The Milner Centre for Evolution and Department of Biology and Biochemistry, University of Bath, Bath, UK
| | - Andrew Preston
- The Milner Centre for Evolution and Department of Biology and Biochemistry, University of Bath, Bath, UK
| | - Caroline Colijn
- Department of Mathematics, Imperial College, London, UK
- * E-mail:
| |
Collapse
|
26
|
Cook DJ, Nielsen J. Genome-scale metabolic models applied to human health and disease. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2017. [DOI: 10.1002/wsbm.1393] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Daniel J Cook
- Department of Biology and Biological Engineering; Chalmers University of Technology; Gothenburg Sweden
| | - Jens Nielsen
- Department of Biology and Biological Engineering; Chalmers University of Technology; Gothenburg Sweden
| |
Collapse
|
27
|
Salinas D, Minor CA, Carlson RP, McCutchen CN, Mumey BM, June RK. Combining Targeted Metabolomic Data with a Model of Glucose Metabolism: Toward Progress in Chondrocyte Mechanotransduction. PLoS One 2017; 12:e0168326. [PMID: 28056047 PMCID: PMC5215894 DOI: 10.1371/journal.pone.0168326] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Accepted: 11/30/2016] [Indexed: 11/19/2022] Open
Abstract
Osteoarthritis is a debilitating disease likely involving altered metabolism of the chondrocytes in articular cartilage. Chondrocytes can respond metabolically to mechanical loads via cellular mechanotransduction, and metabolic changes are significant because they produce the precursors to the tissue matrix necessary for cartilage health. However, a comprehensive understanding of how energy metabolism changes with loading remains elusive. To improve our understanding of chondrocyte mechanotransduction, we developed a computational model to calculate the rate of reactions (i.e. flux) across multiple components of central energy metabolism based on experimental data. We calculated average reaction flux profiles of central metabolism for SW1353 human chondrocytes subjected to dynamic compression for 30 minutes. The profiles were obtained solving a bounded variable linear least squares problem, representing the stoichiometry of human central energy metabolism. Compression synchronized chondrocyte energy metabolism. These data are consistent with dynamic compression inducing early time changes in central energy metabolism geared towards more active protein synthesis. Furthermore, this analysis demonstrates the utility of combining targeted metabolomic data with a computational model to enable rapid analysis of cellular energy utilization.
Collapse
Affiliation(s)
- Daniel Salinas
- Computer Science, Montana State University, Bozeman, MT United States of America
| | - Cody A. Minor
- Mathematics, Montana State University, Bozeman, MT United States of America
| | - Ross P. Carlson
- Chemical & Biological Engineering, Montana State University, Bozeman, MT United States of America
| | - Carley N. McCutchen
- Mechanical & Industrial Engineering, Montana State University, Bozeman, MT United States of America
| | - Brendan M. Mumey
- Computer Science, Montana State University, Bozeman, MT United States of America
| | - Ronald K. June
- Mechanical & Industrial Engineering, Montana State University, Bozeman, MT United States of America
- Department of Cell Biology & Neurosciences, Montana State University, Bozeman, MT United States of America
- Department of Orthopaedics and Sports Medicine, University of Washington, Seattle, WA United States of America
- * E-mail:
| |
Collapse
|
28
|
Hunt KA, Jennings RD, Inskeep WP, Carlson RP. Stoichiometric modelling of assimilatory and dissimilatory biomass utilisation in a microbial community. Environ Microbiol 2016; 18:4946-4960. [PMID: 27387069 PMCID: PMC5629010 DOI: 10.1111/1462-2920.13444] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2016] [Accepted: 06/30/2016] [Indexed: 11/26/2022]
Abstract
Assimilatory and dissimilatory utilisation of autotroph biomass by heterotrophs is a fundamental mechanism for the transfer of nutrients and energy across trophic levels. Metagenome data from a tractable, thermoacidophilic microbial community in Yellowstone National Park was used to build an in silico model to study heterotrophic utilisation of autotroph biomass using elementary flux mode analysis and flux balance analysis. Assimilatory and dissimilatory biomass utilisation was investigated using 29 forms of biomass-derived dissolved organic carbon (DOC) including individual monomer pools, individual macromolecular pools and aggregate biomass. The simulations identified ecologically competitive strategies for utilizing DOC under conditions of varying electron donor, electron acceptor or enzyme limitation. The simulated growth environment affected which form of DOC was the most competitive use of nutrients; for instance, oxygen limitation favoured utilisation of less reduced and fermentable DOC while carbon-limited environments favoured more reduced DOC. Additionally, metabolism was studied considering two encompassing metabolic strategies: simultaneous versus sequential use of DOC. Results of this study bound the transfer of nutrients and energy through microbial food webs, providing a quantitative foundation relevant to most microbial ecosystems.
Collapse
Affiliation(s)
- Kristopher A. Hunt
- Center for Biofilm Engineering, Montana State University, Bozeman, MT, USA
- Department of Chemical and Biological Engineering, Montana State University, Bozeman, MT, USA
- Thermal Biology Institute, Montana State University, Bozeman, MT, USA
| | - Ryan deM. Jennings
- Thermal Biology Institute, Montana State University, Bozeman, MT, USA
- Department of Land Resources and Environmental Sciences, Montana State University, Bozeman, MT, USA
| | - William P. Inskeep
- Thermal Biology Institute, Montana State University, Bozeman, MT, USA
- Department of Land Resources and Environmental Sciences, Montana State University, Bozeman, MT, USA
| | - Ross P. Carlson
- Center for Biofilm Engineering, Montana State University, Bozeman, MT, USA
- Department of Chemical and Biological Engineering, Montana State University, Bozeman, MT, USA
- Thermal Biology Institute, Montana State University, Bozeman, MT, USA
| |
Collapse
|
29
|
Lechner A, Brunk E, Keasling JD. The Need for Integrated Approaches in Metabolic Engineering. Cold Spring Harb Perspect Biol 2016; 8:cshperspect.a023903. [PMID: 27527588 DOI: 10.1101/cshperspect.a023903] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
This review highlights state-of-the-art procedures for heterologous small-molecule biosynthesis, the associated bottlenecks, and new strategies that have the potential to accelerate future accomplishments in metabolic engineering. We emphasize that a combination of different approaches over multiple time and size scales must be considered for successful pathway engineering in a heterologous host. We have classified these optimization procedures based on the "system" that is being manipulated: transcriptome, translatome, proteome, or reactome. By bridging multiple disciplines, including molecular biology, biochemistry, biophysics, and computational sciences, we can create an integral framework for the discovery and implementation of novel biosynthetic production routes.
Collapse
Affiliation(s)
- Anna Lechner
- Joint Bioenergy Institute (JBEI), Emeryville, California 94608.,Department of Chemical & Biomolecular Engineering, Department of Bioengineering, University of California, Berkeley, California 94720
| | - Elizabeth Brunk
- Department of Bioengineering, University of California, San Diego, California 92093
| | - Jay D Keasling
- Joint Bioenergy Institute (JBEI), Emeryville, California 94608.,Department of Chemical & Biomolecular Engineering, Department of Bioengineering, University of California, Berkeley, California 94720.,Physical Biosciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720
| |
Collapse
|
30
|
Dersch LM, Beckers V, Wittmann C. Green pathways: Metabolic network analysis of plant systems. Metab Eng 2016; 34:1-24. [DOI: 10.1016/j.ymben.2015.12.001] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2015] [Revised: 11/30/2015] [Accepted: 12/01/2015] [Indexed: 12/18/2022]
|
31
|
Jamshidi N, Raghunathan A. Cell scale host-pathogen modeling: another branch in the evolution of constraint-based methods. Front Microbiol 2015; 6:1032. [PMID: 26500611 PMCID: PMC4594423 DOI: 10.3389/fmicb.2015.01032] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2015] [Accepted: 09/11/2015] [Indexed: 12/12/2022] Open
Abstract
Constraint-based models have become popular methods for systems biology as they enable the integration of complex, disparate datasets in a biologically cohesive framework that also supports the description of biological processes in terms of basic physicochemical constraints and relationships. The scope, scale, and application of genome scale models have grown from single cell bacteria to multi-cellular interaction modeling; host-pathogen modeling represents one of these examples at the current horizon of constraint-based methods. There are now a small number of examples of host-pathogen constraint-based models in the literature, however there has not yet been a definitive description of the methodology required for the functional integration of genome scale models in order to generate simulation capable host-pathogen models. Herein we outline a systematic procedure to produce functional host-pathogen models, highlighting steps which require debugging and iterative revisions in order to successfully build a functional model. The construction of such models will enable the exploration of host-pathogen interactions by leveraging the growing wealth of omic data in order to better understand mechanism of infection and identify novel therapeutic strategies.
Collapse
Affiliation(s)
- Neema Jamshidi
- Institute of Engineering in Medicine, University of California San Diego, La Jolla, CA, USA ; Department of Radiological Sciences, University of California, Los Angeles Los Angeles, CA, USA
| | - Anu Raghunathan
- Chemical Engineering Division, National Chemical Laboratory Pune, India
| |
Collapse
|
32
|
Gene knockout identification for metabolite production improvement using a hybrid of genetic ant colony optimization and flux balance analysis. BIOTECHNOL BIOPROC E 2015. [DOI: 10.1007/s12257-015-0276-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
33
|
Vaitheesvaran B, Xu J, Yee J, Q-Y L, Go VL, Xiao GG, Lee WN. The Warburg effect: a balance of flux analysis. Metabolomics 2015; 11:787-796. [PMID: 26207106 PMCID: PMC4507278 DOI: 10.1007/s11306-014-0760-9] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Cancer metabolism is characterized by increased macromolecular syntheses through coordinated increases in energy and substrate metabolism. The observation that cancer cells produce lactate in an environment of oxygen sufficiency (aerobic glycolysis) is a central theme of cancer metabolism known as the Warburg effect. Aerobic glycolysis in cancer metabolism is accompanied by increased pentose cycle and anaplerotic activities producing energy and substrates for macromolecular synthesis. How these processes are coordinated is poorly understood. Recent advances have focused on molecular regulation of cancer metabolism by oncogenes and tumor suppressor genes which regulate numerous enzymatic steps of central glucose metabolism. In the past decade, new insights in cancer metabolism have emerged through the application of stable isotopes particularly from 13C carbon tracing. Such studies have provided new evidence for system-wide changes in cancer metabolism in response to chemotherapy. Interestingly, experiments using metabolic inhibitors on individual biochemical pathways all demonstrate similar system-wide effects on cancer metabolism as in targeted therapies. Since biochemical reactions in the Warburg effect place competing demands on available precursors, high energy phosphates and reducing equivalents, the cancer metabolic system must fulfill the condition of balance of flux (homeostasis). In this review, the functions of the pentose cycle and of the tricarboxylic acid (TCA) cycle in cancer metabolism are analyzed from the balance of flux point of view. Anticancer treatments that target molecular signaling pathways or inhibit metabolism alter the invasive or proliferative behavior of the cancer cells by their effects on the balance of flux (homeostasis) of the cancer metabolic phenotype.
Collapse
Affiliation(s)
- B Vaitheesvaran
- Department of Medicine, Diabetes Center, Stable Isotope and
Metabolomics Core Facility, Albert Einstein College of Medicine Diabetes Center,
Bronx, New York, USA
| | - J Xu
- Department of Pathology, University of Southern California, Los
Angeles, Caligornia, USA
| | - J Yee
- Department of Pediatrics, Division of Endocrinology and Metabolism,
University of California, Los Angeles, California, USA
| | - Lu Q-Y
- Department of Medicine, University of California, Los Angeles, CA,
USA
| | - VL Go
- Department of Medicine, University of California, Los Angeles, CA,
USA
| | - G G Xiao
- Functional Genomics/Proteomics Laboratories Creighton University
medical Center, Nebraska, and School of Pharmaceutical Science and Technology at
Dalian University of Technology, Dalian, China
| | - WN Lee
- LA Biomedical Research Institute, Torrance, CA, USA and Department
of Pediatrics, Division of Endocrinology and Metabolism, University of California,
Los Angeles, California USA
| |
Collapse
|
34
|
O'Brien EJ, Monk JM, Palsson BO. Using Genome-scale Models to Predict Biological Capabilities. Cell 2015; 161:971-987. [PMID: 26000478 DOI: 10.1016/j.cell.2015.05.019] [Citation(s) in RCA: 439] [Impact Index Per Article: 48.8] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2014] [Indexed: 11/29/2022]
Abstract
Constraint-based reconstruction and analysis (COBRA) methods at the genome scale have been under development since the first whole-genome sequences appeared in the mid-1990s. A few years ago, this approach began to demonstrate the ability to predict a range of cellular functions, including cellular growth capabilities on various substrates and the effect of gene knockouts at the genome scale. Thus, much interest has developed in understanding and applying these methods to areas such as metabolic engineering, antibiotic design, and organismal and enzyme evolution. This Primer will get you started.
Collapse
Affiliation(s)
- Edward J O'Brien
- Department of Bioengineering, University of California, San Diego, La Jolla, CA 92093, USA; Bioinformatics and Systems Biology Program, University of California, San Diego, La Jolla, CA 92093, USA
| | - Jonathan M Monk
- Department of Bioengineering, University of California, San Diego, La Jolla, CA 92093, USA; Department of NanoEngineering, University of California, San Diego, La Jolla, CA 92093, USA
| | - Bernhard O Palsson
- Department of Bioengineering, University of California, San Diego, La Jolla, CA 92093, USA; Department of Pediatrics, University of California, San Diego, La Jolla, CA 92093, USA; Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Lyngby 2800, Denmark.
| |
Collapse
|
35
|
Krumholz EW, Libourel IGL. Sequence-based Network Completion Reveals the Integrality of Missing Reactions in Metabolic Networks. J Biol Chem 2015; 290:19197-207. [PMID: 26041773 DOI: 10.1074/jbc.m114.634121] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2014] [Indexed: 11/06/2022] Open
Abstract
Genome-scale metabolic models are central in connecting genotypes to metabolic phenotypes. However, even for well studied organisms, such as Escherichia coli, draft networks do not contain a complete biochemical network. Missing reactions are referred to as gaps. These gaps need to be filled to enable functional analysis, and gap-filling choices influence model predictions. To investigate whether functional networks existed where all gap-filling reactions were supported by sequence similarity to annotated enzymes, four draft networks were supplemented with all reactions from the Model SEED database for which minimal sequence similarity was found in their genomes. Quadratic programming revealed that the number of reactions that could partake in a gap-filling solution was vast: 3,270 in the case of E. coli, where 72% of the metabolites in the draft network could connect a gap-filling solution. Nonetheless, no network could be completed without the inclusion of orphaned enzymes, suggesting that parts of the biochemistry integral to biomass precursor formation are uncharacterized. However, many gap-filling reactions were well determined, and the resulting networks showed improved prediction of gene essentiality compared with networks generated through canonical gap filling. In addition, gene essentiality predictions that were sensitive to poorly determined gap-filling reactions were of poor quality, suggesting that damage to the network structure resulting from the inclusion of erroneous gap-filling reactions may be predictable.
Collapse
Affiliation(s)
| | - Igor G L Libourel
- From the Department of Plant Biology and the Biotechnology Institute, University of Minnesota, Saint Paul, Minnesota 55108
| |
Collapse
|
36
|
Monk J, Nogales J, Palsson BO. Optimizing genome-scale network reconstructions. Nat Biotechnol 2015; 32:447-52. [PMID: 24811519 DOI: 10.1038/nbt.2870] [Citation(s) in RCA: 138] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Jonathan Monk
- 1] Department of Bioengineering, University of California, San Diego, La Jolla, California, USA. [2]
| | - Juan Nogales
- 1] Department of Bioengineering, University of California, San Diego, La Jolla, California, USA. [2] Department of Environmental Biology, Centro de Investigaciones Biológicas, CSIC, Madrid, Spain. [3]
| | - Bernhard O Palsson
- Department of Bioengineering, University of California, San Diego, La Jolla, California, USA
| |
Collapse
|
37
|
O'Brien EJ, Palsson BO. Computing the functional proteome: recent progress and future prospects for genome-scale models. Curr Opin Biotechnol 2015; 34:125-34. [PMID: 25576845 DOI: 10.1016/j.copbio.2014.12.017] [Citation(s) in RCA: 56] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2014] [Revised: 12/16/2014] [Accepted: 12/17/2014] [Indexed: 11/18/2022]
Abstract
Constraint-based models enable the computation of feasible, optimal, and realized biological phenotypes from reaction network reconstructions and constraints on their operation. To date, stoichiometric reconstructions have largely focused on metabolism, resulting in genome-scale metabolic models (M-Models). Recent expansions in network content to encompass proteome synthesis have resulted in models of metabolism and protein expression (ME-Models). ME-Models advance the predictions possible with constraint-based models from network flux states to the spatially resolved molecular composition of a cell. Specifically, ME-Models enable the prediction of transcriptome and proteome allocation and limitations, and basal expression states and regulatory needs. Continued expansion in reconstruction content and constraints will result in an increasingly refined representation of cellular composition and behavior.
Collapse
Affiliation(s)
- Edward J O'Brien
- Bioinformatics and Systems Biology Program, University of California, San Diego, United States; Department of Bioengineering, University of California, San Diego, United States
| | - Bernhard O Palsson
- Bioinformatics and Systems Biology Program, University of California, San Diego, United States; Department of Bioengineering, University of California, San Diego, United States; Department of Pediatrics, University of California, San Diego, United States; Novo Nordisk Center for Biosustainability, The Danish Technical University, Denmark.
| |
Collapse
|
38
|
Perturbation Experiments: Approaches for Metabolic Pathway Analysis in Bioreactors. ADVANCES IN BIOCHEMICAL ENGINEERING/BIOTECHNOLOGY 2015; 152:91-136. [PMID: 25981857 DOI: 10.1007/10_2015_326] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
In the last decades, targeted metabolic engineering of microbial cells has become one of the major tools in bioprocess design and optimization. For successful application, a detailed knowledge is necessary about the relevant metabolic pathways and their regulation inside the cells. Since in vitro experiments cannot display process conditions and behavior properly, process data about the cells' metabolic state have to be collected in vivo. For this purpose, special techniques and methods are necessary. Therefore, most techniques enabling in vivo characterization of metabolic pathways rely on perturbation experiments, which can be divided into dynamic and steady-state approaches. To avoid any process disturbance, approaches which enable perturbation of cell metabolism in parallel to the continuing production process are reasonable. Furthermore, the fast dynamics of microbial production processes amplifies the need of parallelized data generation. These points motivate the development of a parallelized approach for multiple metabolic perturbation experiments outside the operating production reactor. An appropriate approach for in vivo characterization of metabolic pathways is presented and applied exemplarily to a microbial L-phenylalanine production process on a 15 L-scale.
Collapse
|
39
|
Colombié S, Nazaret C, Bénard C, Biais B, Mengin V, Solé M, Fouillen L, Dieuaide-Noubhani M, Mazat JP, Beauvoit B, Gibon Y. Modelling central metabolic fluxes by constraint-based optimization reveals metabolic reprogramming of developing Solanum lycopersicum (tomato) fruit. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2015; 81:24-39. [PMID: 25279440 PMCID: PMC4309433 DOI: 10.1111/tpj.12685] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2014] [Revised: 09/19/2014] [Accepted: 09/19/2014] [Indexed: 05/18/2023]
Abstract
Modelling of metabolic networks is a powerful tool to analyse the behaviour of developing plant organs, including fruits. Guided by our current understanding of heterotrophic metabolism of plant cells, a medium-scale stoichiometric model, including the balance of co-factors and energy, was constructed in order to describe metabolic shifts that occur through the nine sequential stages of Solanum lycopersicum (tomato) fruit development. The measured concentrations of the main biomass components and the accumulated metabolites in the pericarp, determined at each stage, were fitted in order to calculate, by derivation, the corresponding external fluxes. They were used as constraints to solve the model by minimizing the internal fluxes. The distribution of the calculated fluxes of central metabolism were then analysed and compared with known metabolic behaviours. For instance, the partition of the main metabolic pathways (glycolysis, pentose phosphate pathway, etc.) was relevant throughout fruit development. We also predicted a valid import of carbon and nitrogen by the fruit, as well as a consistent CO2 release. Interestingly, the energetic balance indicates that excess ATP is dissipated just before the onset of ripening, supporting the concept of the climacteric crisis. Finally, the apparent contradiction between calculated fluxes with low values compared with measured enzyme capacities suggest a complex reprogramming of the metabolic machinery during fruit development. With a powerful set of experimental data and an accurate definition of the metabolic system, this work provides important insight into the metabolic and physiological requirements of the developing tomato fruits.
Collapse
Affiliation(s)
- Sophie Colombié
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
- *For correspondence (e-mail )
| | - Christine Nazaret
- Institut de Mathématiques de Bordeaux, ENSTBB-Institut Polytechnique de Bordeaux351 Cours de la Liberation, Talence, F-33400, France
| | - Camille Bénard
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
| | - Benoît Biais
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
| | - Virginie Mengin
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
| | - Marion Solé
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
| | - Laëtitia Fouillen
- CNRS, UMR 5200Laboratoire de Biogenèse Membranaire, Villenave D'Ornon, F-33883, France
- Univ. Bordeaux146 rue Léo-Saignat, Bordeaux Cedex, F-33076, France
| | - Martine Dieuaide-Noubhani
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
- Univ. Bordeaux146 rue Léo-Saignat, Bordeaux Cedex, F-33076, France
| | - Jean-Pierre Mazat
- Univ. Bordeaux146 rue Léo-Saignat, Bordeaux Cedex, F-33076, France
- IBGC-CNRS1 rue Camille Saint-Saëns, Bordeaux Cedex, F-33077, France
| | - Bertrand Beauvoit
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
- Univ. Bordeaux146 rue Léo-Saignat, Bordeaux Cedex, F-33076, France
| | - Yves Gibon
- INRAUMR 1332 Biologie du Fruit et Pathologie, Villenave d'Ornon, F-33883, France
| |
Collapse
|
40
|
Shi C, Yin J, Liu Z, Wu JX, Zhao Q, Ren J, Yao N. A systematic simulation of the effect of salicylic acid on sphingolipid metabolism. FRONTIERS IN PLANT SCIENCE 2015; 6:186. [PMID: 25859253 PMCID: PMC4373270 DOI: 10.3389/fpls.2015.00186] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/20/2014] [Accepted: 03/08/2015] [Indexed: 05/04/2023]
Abstract
The phytohormone salicylic acid (SA) affects plant development and defense responses. Recent studies revealed that SA also participates in the regulation of sphingolipid metabolism, but the details of this regulation remain to beexplored. Here, we use in silico Flux Balance Analysis (FBA) with published microarray data to construct a whole-cell simulation model, including 23 pathways, 259 reactions, and 172 metabolites, to predict the alterations in flux of major sphingolipid species after treatment with exogenous SA. This model predicts significant changes in fluxes of certain sphingolipid species after SA treatment, changes that likely trigger downstream physiological and phenotypic effects. To validate the simulation, we used (15)N-labeled metabolic turnover analysis to measure sphingolipid contents and turnover rate in Arabidopsis thaliana seedlings treated with SA or the SA analog benzothiadiazole (BTH). The results show that both SA and BTH affect sphingolipid metabolism, altering the concentrations of certain species and also changing the optimal flux distribution and turnover rate of sphingolipids. Our strategy allows us to estimate sphingolipid fluxes on a short time scale and gives us a systemic view of the effect of SA on sphingolipid homeostasis.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Nan Yao
- *Correspondence: Nan Yao, State Key Laboratory of Biocontrol, Guangdong Key Laboratory of Plant Resources, Department of Biological Science and Technology, School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, China
| |
Collapse
|
41
|
Goyal N, Widiastuti H, Karimi IA, Zhou Z. A genome-scale metabolic model of Methanococcus maripaludis S2 for CO2 capture and conversion to methane. MOLECULAR BIOSYSTEMS 2014; 10:1043-54. [PMID: 24553424 DOI: 10.1039/c3mb70421a] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Methane is a major energy source for heating and electricity. Its production by methanogenic bacteria is widely known in nature. M. maripaludis S2 is a fully sequenced hydrogenotrophic methanogen and an excellent laboratory strain with robust genetic tools. However, a quantitative systems biology model to complement these tools is absent in the literature. To understand and enhance its methanogenesis from CO2, this work presents the first constraint-based genome-scale metabolic model (iMM518). It comprises 570 reactions, 556 distinct metabolites, and 518 genes along with gene-protein-reaction (GPR) associations, and covers 30% of open reading frames (ORFs). The model was validated using biomass growth data and experimental phenotypic studies from the literature. Its comparison with the in silico models of Methanosarcina barkeri, Methanosarcina acetivorans, and Sulfolobus solfataricus P2 shows M. maripaludis S2 to be a better organism for producing methane. Using the model, genes essential for growth were identified, and the efficacies of alternative carbon, hydrogen and nitrogen sources were studied. The model can predict the effects of reengineering M. maripaludis S2 to guide or expedite experimental efforts.
Collapse
Affiliation(s)
- Nishu Goyal
- Department of Chemical and Biomolecular Engineering, National University of Singapore, 4 Engineering Drive 4, Singapore 117576.
| | | | | | | |
Collapse
|
42
|
Song HS, Reifman J, Wallqvist A. Prediction of metabolic flux distribution from gene expression data based on the flux minimization principle. PLoS One 2014; 9:e112524. [PMID: 25397773 PMCID: PMC4232356 DOI: 10.1371/journal.pone.0112524] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2014] [Accepted: 10/08/2014] [Indexed: 11/30/2022] Open
Abstract
Prediction of possible flux distributions in a metabolic network provides detailed phenotypic information that links metabolism to cellular physiology. To estimate metabolic steady-state fluxes, the most common approach is to solve a set of macroscopic mass balance equations subjected to stoichiometric constraints while attempting to optimize an assumed optimal objective function. This assumption is justifiable in specific cases but may be invalid when tested across different conditions, cell populations, or other organisms. With an aim to providing a more consistent and reliable prediction of flux distributions over a wide range of conditions, in this article we propose a framework that uses the flux minimization principle to predict active metabolic pathways from mRNA expression data. The proposed algorithm minimizes a weighted sum of flux magnitudes, while biomass production can be bounded to fit an ample range from very low to very high values according to the analyzed context. We have formulated the flux weights as a function of the corresponding enzyme reaction's gene expression value, enabling the creation of context-specific fluxes based on a generic metabolic network. In case studies of wild-type Saccharomyces cerevisiae, and wild-type and mutant Escherichia coli strains, our method achieved high prediction accuracy, as gauged by correlation coefficients and sums of squared error, with respect to the experimentally measured values. In contrast to other approaches, our method was able to provide quantitative predictions for both model organisms under a variety of conditions. Our approach requires no prior knowledge or assumption of a context-specific metabolic functionality and does not require trial-and-error parameter adjustments. Thus, our framework is of general applicability for modeling the transcription-dependent metabolism of bacteria and yeasts.
Collapse
Affiliation(s)
- Hyun-Seob Song
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, U.S. Army Medical Research and Materiel Command, Fort Detrick, Maryland, United States of America
| | - Jaques Reifman
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, U.S. Army Medical Research and Materiel Command, Fort Detrick, Maryland, United States of America
| | - Anders Wallqvist
- Department of Defense Biotechnology High Performance Computing Software Applications Institute, Telemedicine and Advanced Technology Research Center, U.S. Army Medical Research and Materiel Command, Fort Detrick, Maryland, United States of America
- * E-mail:
| |
Collapse
|
43
|
Waldherr S, Oyarzún DA, Bockmayr A. Dynamic optimization of metabolic networks coupled with gene expression. J Theor Biol 2014; 365:469-85. [PMID: 25451533 DOI: 10.1016/j.jtbi.2014.10.035] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2014] [Revised: 10/22/2014] [Accepted: 10/27/2014] [Indexed: 11/24/2022]
Abstract
The regulation of metabolic activity by tuning enzyme expression levels is crucial to sustain cellular growth in changing environments. Metabolic networks are often studied at steady state using constraint-based models and optimization techniques. However, metabolic adaptations driven by changes in gene expression cannot be analyzed by steady state models, as these do not account for temporal changes in biomass composition. Here we present a dynamic optimization framework that integrates the metabolic network with the dynamics of biomass production and composition. An approximation by a timescale separation leads to a coupled model of quasi-steady state constraints on the metabolic reactions, and differential equations for the substrate concentrations and biomass composition. We propose a dynamic optimization approach to determine reaction fluxes for this model, explicitly taking into account enzyme production costs and enzymatic capacity. In contrast to the established dynamic flux balance analysis, our approach allows predicting dynamic changes in both the metabolic fluxes and the biomass composition during metabolic adaptations. Discretization of the optimization problems leads to a linear program that can be efficiently solved. We applied our algorithm in two case studies: a minimal nutrient uptake network, and an abstraction of core metabolic processes in bacteria. In the minimal model, we show that the optimized uptake rates reproduce the empirical Monod growth for bacterial cultures. For the network of core metabolic processes, the dynamic optimization algorithm predicted commonly observed metabolic adaptations, such as a diauxic switch with a preference ranking for different nutrients, re-utilization of waste products after depletion of the original substrate, and metabolic adaptation to an impending nutrient depletion. These examples illustrate how dynamic adaptations of enzyme expression can be predicted solely from an optimization principle.
Collapse
Affiliation(s)
- Steffen Waldherr
- Institute for Automation Engineering, Otto von Guericke University Magdeburg, Universitätsplatz 2, 39106 Magdeburg, Germany.
| | - Diego A Oyarzún
- Department of Mathematics, Imperial College London, SW7 2AZ London, United Kingdom
| | - Alexander Bockmayr
- DFG Research Center Matheon, Freie Universität Berlin, Arnimallee 6, 14195 Berlin, Germany
| |
Collapse
|
44
|
Manor O, Levy R, Borenstein E. Mapping the inner workings of the microbiome: genomic- and metagenomic-based study of metabolism and metabolic interactions in the human microbiome. Cell Metab 2014; 20:742-752. [PMID: 25176148 PMCID: PMC4252837 DOI: 10.1016/j.cmet.2014.07.021] [Citation(s) in RCA: 66] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
The human gut microbiome is a major contributor to human metabolism and health, yet the metabolic processes that are carried out by various community members, the way these members interact with each other and with the host, and the impact of such interactions on the overall metabolic machinery of the microbiome have not yet been mapped. Here, we discuss recent efforts to study the metabolic inner workings of this complex ecosystem. We will specifically highlight two interrelated lines of work, the first aiming to deconvolve the microbiome and to characterize the metabolic capacity of various microbiome species and the second aiming to utilize computational modeling to infer and study metabolic interactions between these species.
Collapse
Affiliation(s)
- Ohad Manor
- Department of Genome Sciences, University of Washington, Seattle, WA 98102, USA
| | - Roie Levy
- Department of Genome Sciences, University of Washington, Seattle, WA 98102, USA
| | - Elhanan Borenstein
- Department of Genome Sciences, University of Washington, Seattle, WA 98102, USA; Department of Computer Science and Engineering, University of Washington, Seattle, WA 98102, USA; Santa Fe Institute, Santa Fe, NM 87501, USA.
| |
Collapse
|
45
|
Benedict MN, Mundy MB, Henry CS, Chia N, Price ND. Likelihood-based gene annotations for gap filling and quality assessment in genome-scale metabolic models. PLoS Comput Biol 2014; 10:e1003882. [PMID: 25329157 PMCID: PMC4199484 DOI: 10.1371/journal.pcbi.1003882] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2014] [Accepted: 08/25/2014] [Indexed: 12/27/2022] Open
Abstract
Genome-scale metabolic models provide a powerful means to harness information from genomes to deepen biological insights. With exponentially increasing sequencing capacity, there is an enormous need for automated reconstruction techniques that can provide more accurate models in a short time frame. Current methods for automated metabolic network reconstruction rely on gene and reaction annotations to build draft metabolic networks and algorithms to fill gaps in these networks. However, automated reconstruction is hampered by database inconsistencies, incorrect annotations, and gap filling largely without considering genomic information. Here we develop an approach for applying genomic information to predict alternative functions for genes and estimate their likelihoods from sequence homology. We show that computed likelihood values were significantly higher for annotations found in manually curated metabolic networks than those that were not. We then apply these alternative functional predictions to estimate reaction likelihoods, which are used in a new gap filling approach called likelihood-based gap filling to predict more genomically consistent solutions. To validate the likelihood-based gap filling approach, we applied it to models where essential pathways were removed, finding that likelihood-based gap filling identified more biologically relevant solutions than parsimony-based gap filling approaches. We also demonstrate that models gap filled using likelihood-based gap filling provide greater coverage and genomic consistency with metabolic gene functions compared to parsimony-based approaches. Interestingly, despite these findings, we found that likelihoods did not significantly affect consistency of gap filled models with Biolog and knockout lethality data. This indicates that the phenotype data alone cannot necessarily be used to discriminate between alternative solutions for gap filling and therefore, that the use of other information is necessary to obtain a more accurate network. All described workflows are implemented as part of the DOE Systems Biology Knowledgebase (KBase) and are publicly available via API or command-line web interface.
Collapse
Affiliation(s)
- Matthew N. Benedict
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
| | - Michael B. Mundy
- Center for Individualized Medicine, Mayo Clinic, Rochester, Minnesota, United States of America
| | - Christopher S. Henry
- Mathematics and Computer Science Division, Argonne National Laboratory, Lemont, Illinois, United States of America
| | - Nicholas Chia
- Center for Individualized Medicine, Mayo Clinic, Rochester, Minnesota, United States of America
- Department of Surgery, Mayo Clinic, Rochester, Minnesota, United States of America
- Department of Physiology and Bioengineering, Mayo Clinic, Rochester, Minnesota, United States of America
- * E-mail: (NC); (NDP)
| | - Nathan D. Price
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois, United States of America
- Institute for Systems Biology, Seattle, Washington, United States of America
- * E-mail: (NC); (NDP)
| |
Collapse
|
46
|
Khodayari A, Zomorrodi AR, Liao JC, Maranas CD. A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data. Metab Eng 2014; 25:50-62. [DOI: 10.1016/j.ymben.2014.05.014] [Citation(s) in RCA: 145] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2014] [Revised: 04/17/2014] [Accepted: 05/28/2014] [Indexed: 01/27/2023]
|
47
|
Chiu HC, Levy R, Borenstein E. Emergent biosynthetic capacity in simple microbial communities. PLoS Comput Biol 2014; 10:e1003695. [PMID: 24992662 PMCID: PMC4084645 DOI: 10.1371/journal.pcbi.1003695] [Citation(s) in RCA: 69] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2013] [Accepted: 05/16/2014] [Indexed: 12/22/2022] Open
Abstract
Microbes have an astonishing capacity to transform their environments. Yet, the metabolic capacity of a single species is limited and the vast majority of microorganisms form complex communities and join forces to exhibit capabilities far exceeding those achieved by any single species. Such enhanced metabolic capacities represent a promising route to many medical, environmental, and industrial applications and call for the development of a predictive, systems-level understanding of synergistic microbial capacity. Here we present a comprehensive computational framework, integrating high-quality metabolic models of multiple species, temporal dynamics, and flux variability analysis, to study the metabolic capacity and dynamics of simple two-species microbial ecosystems. We specifically focus on detecting emergent biosynthetic capacity--instances in which a community growing on some medium produces and secretes metabolites that are not secreted by any member species when growing in isolation on that same medium. Using this framework to model a large collection of two-species communities on multiple media, we demonstrate that emergent biosynthetic capacity is highly prevalent. We identify commonly observed emergent metabolites and metabolic reprogramming patterns, characterizing typical mechanisms of emergent capacity. We further find that emergent secretion tends to occur in two waves, the first as soon as the two organisms are introduced, and the second when the medium is depleted and nutrients become limited. Finally, aiming to identify global community determinants of emergent capacity, we find a marked association between the level of emergent biosynthetic capacity and the functional/phylogenetic distance between community members. Specifically, we demonstrate a "Goldilocks" principle, where high levels of emergent capacity are observed when the species comprising the community are functionally neither too close, nor too distant. Taken together, our results demonstrate the potential to design and engineer synthetic communities capable of novel metabolic activities and point to promising future directions in environmental and clinical bioengineering.
Collapse
Affiliation(s)
- Hsuan-Chao Chiu
- Department of Genome Sciences, University of Washington, Seattle, Washington, United States of America
| | - Roie Levy
- Department of Genome Sciences, University of Washington, Seattle, Washington, United States of America
| | - Elhanan Borenstein
- Department of Genome Sciences, University of Washington, Seattle, Washington, United States of America
- Department of Computer Science and Engineering, University of Washington, Seattle, Washington, United States of America
- Santa Fe Institute, Santa Fe, New Mexico, United States of America
- * E-mail:
| |
Collapse
|
48
|
Stitt M, Gibon Y. Why measure enzyme activities in the era of systems biology? TRENDS IN PLANT SCIENCE 2014; 19:256-65. [PMID: 24332227 DOI: 10.1016/j.tplants.2013.11.003] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/26/2013] [Revised: 11/05/2013] [Accepted: 11/08/2013] [Indexed: 05/22/2023]
Abstract
Information about the abundance and biological activities of proteins is essential to reveal how genes affect phenotypes. Over the past decade, mass spectrometry (MS)-based proteomics has revolutionized the identification and quantification of proteins, and the detection of post-translational modifications. Interpretation of proteomics data depends on information about the biological activities of proteins, which has created a bottleneck in research. This review focuses on enzymes in central metabolism. We examine the methods used for measuring enzyme activities, and discuss how these methods provide information about the kinetic and regulatory properties of enzymes, their turnover, and how this information can be integrated into metabolic models. We also discuss how robotized assays could enable the genetic networks that control enzyme abundance to be analyzed.
Collapse
Affiliation(s)
- Mark Stitt
- Max Planck Institute of Molecular Plant Physiology, Am Muehlenberg 1, 14476 Potsdam-Golm, Germany.
| | - Yves Gibon
- INRA, University of Bordeaux, UMR 1332 Fruit Biology and Pathology, F-33883 Villenave d'Ornon, France
| |
Collapse
|
49
|
Rezola A, Pey J, Tobalina L, Rubio A, Beasley JE, Planes FJ. Advances in network-based metabolic pathway analysis and gene expression data integration. Brief Bioinform 2014; 16:265-79. [DOI: 10.1093/bib/bbu009] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Xi Y, Zhao Y, Wang L, Wang F. Comparison on extreme pathways reveals nature of different biological processes. BMC SYSTEMS BIOLOGY 2014; 8 Suppl 1:S10. [PMID: 24565046 PMCID: PMC4080357 DOI: 10.1186/1752-0509-8-s1-s10] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
Background Constraint-based reconstruction and analysis (COBRA) is used for modeling genome-scale metabolic networks (MNs). In a COBRA model, extreme pathways (ExPas) are the edges of its conical solution space, which is formed by all viable steady-state flux distributions. ExPa analysis has been successfully applied to MNs to reveal their phenotypic capabilities and properties. Recently, the COBRA framework has been extended to transcriptional regulatory networks (TRNs) and transcriptional and translational networks (TTNs), so efforts are needed to determine whether ExPa analysis is also effective on these two types of networks. Results In this paper, the ExPas resulting from the COBRA models of E.coli's MN, TRN and TTN were horizontally compared from 5 aspects: (1) Total number and the ratio of their amount to reaction amount; (2) Length distribution; (3) Reaction participation; (4) Correlated reaction sets (CoSets); (5) interconnectivity degree. Significant discrepancies in above properties were observed during the comparison, which reveals the biological natures of different biological processes. Besides, by demonstrating the application of ExPa analysis on E.coli, we provide a practical guidance of an improved approach to compute ExPas on COBRA models of TRNs. Conclusions ExPas of E.coli's MN, TRN and TTN have different properties, which are strongly connected with various biological natures of biochemical networks, such as topological structure, specificity and redundancy. Our study shows that ExPas are biologically meaningful on the newborn models and suggests the effectiveness of ExPa analysis on them.
Collapse
|