1
|
Theorell A, Jadebeck JF, Wiechert W, McFadden J, Nöh K. Rethinking 13C-metabolic flux analysis - The Bayesian way of flux inference. Metab Eng 2024; 83:137-149. [PMID: 38582144 DOI: 10.1016/j.ymben.2024.03.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2023] [Revised: 03/22/2024] [Accepted: 03/23/2024] [Indexed: 04/08/2024]
Abstract
Metabolic reaction rates (fluxes) play a crucial role in comprehending cellular phenotypes and are essential in areas such as metabolic engineering, biotechnology, and biomedical research. The state-of-the-art technique for estimating fluxes is metabolic flux analysis using isotopic labelling (13C-MFA), which uses a dataset-model combination to determine the fluxes. Bayesian statistical methods are gaining popularity in the field of life sciences, but the use of 13C-MFA is still dominated by conventional best-fit approaches. The slow take-up of Bayesian approaches is, at least partly, due to the unfamiliarity of Bayesian methods to metabolic engineering researchers. To address this unfamiliarity, we here outline similarities and differences between the two approaches and highlight particular advantages of the Bayesian way of flux analysis. With a real-life example, re-analysing a moderately informative labelling dataset of E. coli, we identify situations in which Bayesian methods are advantageous and more informative, pointing to potential pitfalls of current 13C-MFA evaluation approaches. We propose the use of Bayesian model averaging (BMA) for flux inference as a means of overcoming the problem of model uncertainty through its tendency to assign low probabilities to both, models that are unsupported by data, and models that are overly complex. In this capacity, BMA resembles a tempered Ockham's razor. With the tempered razor as a guide, BMA-based 13C-MFA alleviates the problem of model selection uncertainty and is thereby capable of becoming a game changer for metabolic engineering by uncovering new insights and inspiring novel approaches.
Collapse
Affiliation(s)
- Axel Theorell
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany
| | - Johann F Jadebeck
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; Computational Systems Biotechnology (AVT.CSB), RWTH Aachen University, 52062 Aachen, Germany
| | - Wolfgang Wiechert
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany; Computational Systems Biotechnology (AVT.CSB), RWTH Aachen University, 52062 Aachen, Germany
| | - Johnjoe McFadden
- Department of Microbial and Cellular Sciences, University of Surrey, GU2 7XH Guildford, United Kingdom
| | - Katharina Nöh
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, 52425 Jülich, Germany.
| |
Collapse
|
2
|
Backman TWH, Schenk C, Radivojevic T, Ando D, Singh J, Czajka JJ, Costello Z, Keasling JD, Tang Y, Akhmatskaya E, Garcia Martin H. BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale. PLoS Comput Biol 2023; 19:e1011111. [PMID: 37948450 PMCID: PMC10664898 DOI: 10.1371/journal.pcbi.1011111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Revised: 11/22/2023] [Accepted: 09/27/2023] [Indexed: 11/12/2023] Open
Abstract
Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13C Metabolic Flux Analysis (13C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in "non-gaussian" situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P-13C MOMA and P-13C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.
Collapse
Affiliation(s)
- Tyler W. H. Backman
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- Biofuels and Bioproducts Division, Joint BioEnergy Institute, Emeryville, California, United States of America
| | - Christina Schenk
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- BCAM, Basque Center for Applied Mathematics, Bilbao, Spain
- DOE Agile BioFoundry, Emeryville, California, United States of America
| | - Tijana Radivojevic
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- Biofuels and Bioproducts Division, Joint BioEnergy Institute, Emeryville, California, United States of America
- DOE Agile BioFoundry, Emeryville, California, United States of America
| | - David Ando
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- Biofuels and Bioproducts Division, Joint BioEnergy Institute, Emeryville, California, United States of America
| | - Jahnavi Singh
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, California, United States of America
| | - Jeffrey J. Czajka
- Department of Energy, Environmental and Chemical Engineering, Washington University in St. Louis, St. Louis, Missouri, United States of America
| | - Zak Costello
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- Biofuels and Bioproducts Division, Joint BioEnergy Institute, Emeryville, California, United States of America
- DOE Agile BioFoundry, Emeryville, California, United States of America
| | - Jay D. Keasling
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- Biofuels and Bioproducts Division, Joint BioEnergy Institute, Emeryville, California, United States of America
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California, United States of America
- Department of Bioengineering, University of California, Berkeley, California, United States of America
- QB3 Institute, University of California, Berkeley, California, United States of America
- Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, Copenhagen, Denmark
- Center for Synthetic Biochemistry, Institute for Synthetic Biology, Shenzhen Institutes for Advanced Technologies, Shenzhen, China
| | - Yinjie Tang
- Department of Energy, Environmental and Chemical Engineering, Washington University in St. Louis, St. Louis, Missouri, United States of America
| | - Elena Akhmatskaya
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- BCAM, Basque Center for Applied Mathematics, Bilbao, Spain
- IKERBASQUE, Basque Foundation for Science, Bilbao, Spain
| | - Hector Garcia Martin
- Biological Systems and Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, California, United States of America
- Biofuels and Bioproducts Division, Joint BioEnergy Institute, Emeryville, California, United States of America
- BCAM, Basque Center for Applied Mathematics, Bilbao, Spain
- DOE Agile BioFoundry, Emeryville, California, United States of America
| |
Collapse
|
3
|
Copeland CA, Olenchock BA, Ziehr D, McGarrity S, Leahy K, Young JD, Loscalzo J, Oldham WM. MYC overrides HIF-1α to regulate proliferating primary cell metabolism in hypoxia. eLife 2023; 12:e82597. [PMID: 37428010 DOI: 10.7554/elife.82597] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2022] [Accepted: 06/27/2023] [Indexed: 07/11/2023] Open
Abstract
Hypoxia requires metabolic adaptations to sustain energetically demanding cellular activities. While the metabolic consequences of hypoxia have been studied extensively in cancer cell models, comparatively little is known about how primary cell metabolism responds to hypoxia. Thus, we developed metabolic flux models for human lung fibroblast and pulmonary artery smooth muscle cells proliferating in hypoxia. Unexpectedly, we found that hypoxia decreased glycolysis despite activation of hypoxia-inducible factor 1α (HIF-1α) and increased glycolytic enzyme expression. While HIF-1α activation in normoxia by prolyl hydroxylase (PHD) inhibition did increase glycolysis, hypoxia blocked this effect. Multi-omic profiling revealed distinct molecular responses to hypoxia and PHD inhibition, and suggested a critical role for MYC in modulating HIF-1α responses to hypoxia. Consistent with this hypothesis, MYC knockdown in hypoxia increased glycolysis and MYC over-expression in normoxia decreased glycolysis stimulated by PHD inhibition. These data suggest that MYC signaling in hypoxia uncouples an increase in HIF-dependent glycolytic gene transcription from glycolytic flux.
Collapse
Affiliation(s)
- Courtney A Copeland
- Department of Medicine, Brigham and Women's Hospital, Boston, United States
- Department of Medicine, Harvard Medical School, Boston, United States
| | - Benjamin A Olenchock
- Department of Medicine, Brigham and Women's Hospital, Boston, United States
- Department of Medicine, Harvard Medical School, Boston, United States
| | - David Ziehr
- Department of Medicine, Brigham and Women's Hospital, Boston, United States
- Department of Medicine, Harvard Medical School, Boston, United States
- Department of Medicine, Massachusetts General Hospital, Boston, United States
| | - Sarah McGarrity
- Department of Medicine, Brigham and Women's Hospital, Boston, United States
- Department of Medicine, Harvard Medical School, Boston, United States
- Center for Systems Biology, School of Health Sciences, University of Iceland, Reykjavik, Iceland
| | - Kevin Leahy
- Department of Medicine, Brigham and Women's Hospital, Boston, United States
- Department of Medicine, Harvard Medical School, Boston, United States
| | - Jamey D Young
- Departments of Chemical & Biomolecular Engineering and Molecular Physiology & Biophysics, Vanderbilt University, Nashville, United States
| | - Joseph Loscalzo
- Department of Medicine, Brigham and Women's Hospital, Boston, United States
- Department of Medicine, Harvard Medical School, Boston, United States
| | - William M Oldham
- Department of Medicine, Brigham and Women's Hospital, Boston, United States
- Department of Medicine, Harvard Medical School, Boston, United States
| |
Collapse
|
4
|
Moiz B, Sriram G, Clyne AM. Interpreting metabolic complexity via isotope-assisted metabolic flux analysis. Trends Biochem Sci 2023; 48:553-567. [PMID: 36863894 PMCID: PMC10182253 DOI: 10.1016/j.tibs.2023.02.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2022] [Revised: 01/22/2023] [Accepted: 02/03/2023] [Indexed: 03/04/2023]
Abstract
Isotope-assisted metabolic flux analysis (iMFA) is a powerful method to mathematically determine the metabolic fluxome from experimental isotope labeling data and a metabolic network model. While iMFA was originally developed for industrial biotechnological applications, it is increasingly used to analyze eukaryotic cell metabolism in physiological and pathological states. In this review, we explain how iMFA estimates the intracellular fluxome, including data and network model (inputs), the optimization-based data fitting (process), and the flux map (output). We then describe how iMFA enables analysis of metabolic complexities and discovery of metabolic pathways. Our goal is to expand the use of iMFA in metabolism research, which is essential to maximizing the impact of metabolic experiments and continuing to advance iMFA and biocomputational techniques.
Collapse
Affiliation(s)
- Bilal Moiz
- Department of Bioengineering, University of Maryland, College Park, MD 20742, USA
| | - Ganesh Sriram
- Department of Chemical and Biomolecular Engineering, University of Maryland, College Park, MD 20742, USA
| | - Alisa Morss Clyne
- Department of Bioengineering, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
5
|
Borah Slater K, Beyß M, Xu Y, Barber J, Costa C, Newcombe J, Theorell A, Bailey MJ, Beste DJV, McFadden J, Nöh K. One-shot 13 C 15 N-metabolic flux analysis for simultaneous quantification of carbon and nitrogen flux. Mol Syst Biol 2023; 19:e11099. [PMID: 36705093 PMCID: PMC9996240 DOI: 10.15252/msb.202211099] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2022] [Revised: 01/06/2023] [Accepted: 01/13/2023] [Indexed: 01/28/2023] Open
Abstract
Metabolic flux is the final output of cellular regulation and has been extensively studied for carbon but much less is known about nitrogen, which is another important building block for living organisms. For the tuberculosis pathogen, this is particularly important in informing the development of effective drugs targeting the pathogen's metabolism. Here we performed 13 C15 N dual isotopic labeling of Mycobacterium bovis BCG steady state cultures, quantified intracellular carbon and nitrogen fluxes and inferred reaction bidirectionalities. This was achieved by model scope extension and refinement, implemented in a multi-atom transition model, within the statistical framework of Bayesian model averaging (BMA). Using BMA-based 13 C15 N-metabolic flux analysis, we jointly resolve carbon and nitrogen fluxes quantitatively. We provide the first nitrogen flux distributions for amino acid and nucleotide biosynthesis in mycobacteria and establish glutamate as the central node for nitrogen metabolism. We improved resolution of the notoriously elusive anaplerotic node in central carbon metabolism and revealed possible operation modes. Our study provides a powerful and statistically rigorous platform to simultaneously infer carbon and nitrogen metabolism in any biological system.
Collapse
Affiliation(s)
| | - Martin Beyß
- Forschungszentrum Jülich GmbH, Institute of Bio- and Geosciences, IBG-1: Biotechnology, Jülich, Germany.,Computational Systems Biotechnology, RWTH Aachen University, Aachen, Germany
| | - Ye Xu
- Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Jim Barber
- Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Catia Costa
- Faculty of Engineering and Physical Sciences, University of Surrey, Guildford, UK
| | - Jane Newcombe
- Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Axel Theorell
- Forschungszentrum Jülich GmbH, Institute of Bio- and Geosciences, IBG-1: Biotechnology, Jülich, Germany
| | - Melanie J Bailey
- Faculty of Engineering and Physical Sciences, University of Surrey, Guildford, UK
| | - Dany J V Beste
- Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Johnjoe McFadden
- Faculty of Health and Medical Sciences, University of Surrey, Guildford, UK
| | - Katharina Nöh
- Forschungszentrum Jülich GmbH, Institute of Bio- and Geosciences, IBG-1: Biotechnology, Jülich, Germany
| |
Collapse
|
6
|
Khana DB, Callaghan MM, Amador-Noguez D. Novel computational and experimental approaches for investigating the thermodynamics of metabolic networks. Curr Opin Microbiol 2021; 66:21-31. [PMID: 34974376 DOI: 10.1016/j.mib.2021.11.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Revised: 11/08/2021] [Accepted: 11/10/2021] [Indexed: 11/18/2022]
Abstract
Thermodynamic analysis of metabolic networks has emerged as a useful new tool for pathway design and metabolic engineering. Understanding the relationship between the thermodynamic driving force of biochemical reactions and metabolic flux has generated new insights regarding the design principles of microbial carbon metabolism. This review summarizes the various lessons that can be obtained from the thermodynamic analysis of metabolic pathways, illustrates concepts of computational thermodynamic tools, and highlights recent applications of thermodynamic analysis to pathway design in industrially relevant microbes.
Collapse
Affiliation(s)
- Daven B Khana
- Department of Bacteriology, University of Wisconsin-Madison, Madison, WI, USA; Great Lakes Bioenergy Research Center, University of Wisconsin-Madison, Madison, WI, USA; Microbiology Doctoral Training Program, University of Wisconsin-Madison, Madison, WI, USA
| | - Melanie M Callaghan
- Department of Bacteriology, University of Wisconsin-Madison, Madison, WI, USA; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Daniel Amador-Noguez
- Department of Bacteriology, University of Wisconsin-Madison, Madison, WI, USA; Great Lakes Bioenergy Research Center, University of Wisconsin-Madison, Madison, WI, USA; Center for Bioenergy Innovation, Oak Ridge National Laboratory, Oak Ridge, TN, USA.
| |
Collapse
|
7
|
de Oliveira RD, Novello V, da Silva LF, Gomez JGC, Le Roux GAC. Glucose metabolism in Pseudomonas aeruginosa is cyclic when producing Polyhydroxyalkanoates and Rhamnolipids. J Biotechnol 2021; 342:54-63. [PMID: 34687809 DOI: 10.1016/j.jbiotec.2021.10.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Revised: 10/11/2021] [Accepted: 10/14/2021] [Indexed: 10/20/2022]
Abstract
Pseudomonas aeruginosa is an important chassis for production of polyhydroxyalkanoates (PHA) and rhamnolipids (RHL). Advances in the understanding of the biosynthesis metabolism of these biocompounds are crucial for increasing yield. 13C-Metabolic Flux Ratio Analysis (13C-MFA) is a technique to estimate in vivo metabolic fluxes ratios. PHA and RHL are essentially non-growth associated products of biotechnological interest and both contain hydroxyalkanoates (HAs), whose labeling patterns could be accessed by GC-MS. In this study, to reveal the relative contributions of the Entner-Doudoroff (ED) pathway and the non-oxidative Pentose Phosphate (PP) pathway to PHA and RHL production, 13C-MFA was performed in Pseudomonas aeruginosa LFM634 when supplied with labeled glucose. This bacterial strain lacks both functional EMP and the oxidative PP branch. Labeling patterns in HAs were measured. Experiments with [U-13C] glucose indicated a low flux though PP pathway. An optimal design of labeling experiment showed that [6-13C] glucose would be the best substrate to enable an estimation of the ED flux with high accuracy. Results of experiments performed with this isotope indicated that about two-thirds of glyceraldehyde 3-phosphate is recycled through a cyclic ED architecture, suggesting that P. aeruginosa utilizes that cycle to regulate the NADPH/Acetyl-CoA ratio for PHA and RHL biosynthesis.
Collapse
Affiliation(s)
| | - Vânia Novello
- Department of Microbiology, Institute of Biomedical Sciences, University of São Paulo, Brazil
| | | | | | | |
Collapse
|
8
|
Wiechert W, Nöh K. Quantitative Metabolic Flux Analysis Based on Isotope Labeling. Metab Eng 2021. [DOI: 10.1002/9783527823468.ch3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
|
9
|
Wegener G, Gropp J, Taubner H, Halevy I, Elvert M. Sulfate-dependent reversibility of intracellular reactions explains the opposing isotope effects in the anaerobic oxidation of methane. SCIENCE ADVANCES 2021; 7:7/19/eabe4939. [PMID: 33952515 PMCID: PMC8099194 DOI: 10.1126/sciadv.abe4939] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Accepted: 03/17/2021] [Indexed: 06/01/2023]
Abstract
The anaerobic oxidation of methane (AOM) is performed by methanotrophic archaea (ANME) in distinct sulfate-methane interfaces of marine sediments. In these interfaces, AOM often appears to deplete methane in the heavy isotopes toward isotopic compositions similar to methanogenesis. Here, we shed light on this effect and its physiological underpinnings using a thermophilic ANME-1-dominated culture. At high sulfate concentrations, residual methane is enriched in both 13C and 2H (13α = 1.016 and 2α = 1.155), as observed previously. In contrast, at low sulfate concentrations, the residual methane is substantially depleted in 13C (13α = 0.977) and, to a lesser extent, in 2H. Using a biochemical-isotopic model, we explain the sulfate dependence of the net isotopic fractionation through the thermodynamic drive of the involved intracellular reactions. Our findings relate these isotopic patterns to the physiology and environment of the ANME, thereby explaining a commonly observed isotopic enigma.
Collapse
Affiliation(s)
- Gunter Wegener
- Max Planck Institute for Marine Microbiology, 28359 Bremen, Germany.
- MARUM, Center for Marine Environmental Sciences, University of Bremen, 28359 Bremen, Germany
| | - Jonathan Gropp
- Department of Earth and Planetary Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel.
| | - Heidi Taubner
- MARUM, Center for Marine Environmental Sciences, University of Bremen, 28359 Bremen, Germany
- Faculty of Geosciences, University of Bremen, 28359 Bremen, Germany
| | - Itay Halevy
- Department of Earth and Planetary Sciences, Weizmann Institute of Science, Rehovot 7610001, Israel
| | - Marcus Elvert
- MARUM, Center for Marine Environmental Sciences, University of Bremen, 28359 Bremen, Germany
- Faculty of Geosciences, University of Bremen, 28359 Bremen, Germany
| |
Collapse
|
10
|
Theorell A, Nöh K. Reversible jump MCMC for multi-model inference in Metabolic Flux Analysis. Bioinformatics 2019; 36:232-240. [DOI: 10.1093/bioinformatics/btz500] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2018] [Revised: 03/15/2019] [Accepted: 06/12/2019] [Indexed: 02/06/2023] Open
Abstract
Abstract
Motivation
The validity of model based inference, as used in systems biology, depends on the underlying model formulation. Often, a vast number of competing models is available, that are built on different assumptions, all consistent with the existing knowledge about the studied biological phenomenon. As a remedy for this, Bayesian Model Averaging (BMA) facilitates parameter and structural inferences based on multiple models simultaneously. However, in fields where a vast number of alternative, high-dimensional and non-linear models are involved, the BMA-based inference task is computationally very challenging.
Results
Here we use BMA in the complex setting of Metabolic Flux Analysis (MFA) to infer whether potentially reversible reactions proceed uni- or bidirectionally, using 13C labeling data and metabolic networks. BMA is applied on a large set of candidate models with differing directionality settings, using a tailored multi-model Markov Chain Monte Carlo (MCMC) approach. The applicability of our algorithm is shown by inferring the in vivo probability of reaction bidirectionalities in a realistic network setup, thereby extending the scope of 13C MFA from parameter to structural inference.
Supplementary information
Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Axel Theorell
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
| | - Katharina Nöh
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich 52428, Germany
| |
Collapse
|
11
|
Beyß M, Azzouzi S, Weitzel M, Wiechert W, Nöh K. The Design of FluxML: A Universal Modeling Language for 13C Metabolic Flux Analysis. Front Microbiol 2019; 10:1022. [PMID: 31178829 PMCID: PMC6543931 DOI: 10.3389/fmicb.2019.01022] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2018] [Accepted: 04/24/2019] [Indexed: 12/16/2022] Open
Abstract
13C metabolic flux analysis (MFA) is the method of choice when a detailed inference of intracellular metabolic fluxes in living organisms under metabolic quasi-steady state conditions is desired. Being continuously developed since two decades, the technology made major contributions to the quantitative characterization of organisms in all fields of biotechnology and health-related research. 13C MFA, however, stands out from other "-omics sciences," in that it requires not only experimental-analytical data, but also mathematical models and a computational toolset to infer the quantities of interest, i.e., the metabolic fluxes. At present, these models cannot be conveniently exchanged between different labs. Here, we present the implementation-independent model description language FluxML for specifying 13C MFA models. The core of FluxML captures the metabolic reaction network together with atom mappings, constraints on the model parameters, and the wealth of data configurations. In particular, we describe the governing design processes that shaped the FluxML language. We demonstrate the utility of FluxML to represent many contemporary experimental-analytical requirements in the field of 13C MFA. The major aim of FluxML is to offer a sound, open, and future-proof language to unambiguously express and conserve all the necessary information for model re-use, exchange, and comparison. Along with FluxML, several powerful computational tools are supplied for easy handling, but also to maintain a maximum of flexibility. Altogether, the FluxML collection is an "all-around carefree package" for 13C MFA modelers. We believe that FluxML improves scientific productivity as well as transparency and therewith contributes to the efficiency and reproducibility of computational modeling efforts in the field of 13C MFA.
Collapse
Affiliation(s)
- Martin Beyß
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Salah Azzouzi
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Michael Weitzel
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Wolfgang Wiechert
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich, Germany.,Computational Systems Biotechnology (AVT.CSB), RWTH Aachen University, Aachen, Germany
| | - Katharina Nöh
- Institute of Bio- and Geosciences, IBG-1: Biotechnology, Forschungszentrum Jülich GmbH, Jülich, Germany
| |
Collapse
|
12
|
Alkan HF, Walter KE, Luengo A, Madreiter-Sokolowski CT, Stryeck S, Lau AN, Al-Zoughbi W, Lewis CA, Thomas CJ, Hoefler G, Graier WF, Madl T, Vander Heiden MG, Bogner-Strauss JG. Cytosolic Aspartate Availability Determines Cell Survival When Glutamine Is Limiting. Cell Metab 2018; 28:706-720.e6. [PMID: 30122555 PMCID: PMC6390946 DOI: 10.1016/j.cmet.2018.07.021] [Citation(s) in RCA: 112] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/17/2017] [Revised: 05/29/2018] [Accepted: 07/29/2018] [Indexed: 12/23/2022]
Abstract
Mitochondrial function is important for aspartate biosynthesis in proliferating cells. Here, we show that mitochondrial aspartate export via the aspartate-glutamate carrier 1 (AGC1) supports cell proliferation and cellular redox homeostasis. Insufficient cytosolic aspartate delivery leads to cell death when TCA cycle carbon is reduced following glutamine withdrawal and/or glutaminase inhibition. Moreover, loss of AGC1 reduces allograft tumor growth that is further compromised by treatment with the glutaminase inhibitor CB-839. Together, these findings argue that mitochondrial aspartate export sustains cell survival in low-glutamine environments and AGC1 inhibition can synergize with glutaminase inhibition to limit tumor growth.
Collapse
Affiliation(s)
- H Furkan Alkan
- Institute of Biochemistry, Graz University of Technology, Humboldtstrasse 46/III, 8010 Graz, Austria; The Koch Institute for Integrative Cancer Research and Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Katharina E Walter
- Institute of Biochemistry, Graz University of Technology, Humboldtstrasse 46/III, 8010 Graz, Austria
| | - Alba Luengo
- The Koch Institute for Integrative Cancer Research and Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Corina T Madreiter-Sokolowski
- Gottfried Schatz Research Center, Molecular Biology and Biochemistry, Medical University of Graz, Neue Stiftingtalstrasse 6/6, A-8010 Graz, Austria
| | - Sarah Stryeck
- Gottfried Schatz Research Center, Molecular Biology and Biochemistry, Medical University of Graz, Neue Stiftingtalstrasse 6/6, A-8010 Graz, Austria
| | - Allison N Lau
- The Koch Institute for Integrative Cancer Research and Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Wael Al-Zoughbi
- Diagnostic and Research Institute of Pathology, Medical University of Graz, Neue Stiftingtalstraße 6, A-8010 Graz, Austria
| | - Caroline A Lewis
- Whitehead Institute for Biomedical Research, 455 Main Street, Cambridge, MA 02142, USA
| | - Craig J Thomas
- Division of Preclinical Innovation, National Center for Advancing Translational Sciences, National Institutes of Health, Bethesda, MD 20892, USA; Lymphoid Malignancies Branch, National Cancer Institute, Bethesda, MD 20892, USA
| | - Gerald Hoefler
- Diagnostic and Research Institute of Pathology, Medical University of Graz, Neue Stiftingtalstraße 6, A-8010 Graz, Austria; BioTechMed-Graz, Graz, Austria
| | - Wolfgang F Graier
- Gottfried Schatz Research Center, Molecular Biology and Biochemistry, Medical University of Graz, Neue Stiftingtalstrasse 6/6, A-8010 Graz, Austria; BioTechMed-Graz, Graz, Austria
| | - Tobias Madl
- Gottfried Schatz Research Center, Molecular Biology and Biochemistry, Medical University of Graz, Neue Stiftingtalstrasse 6/6, A-8010 Graz, Austria; BioTechMed-Graz, Graz, Austria
| | - Matthew G Vander Heiden
- The Koch Institute for Integrative Cancer Research and Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA; Dana-Farber Cancer Institute, Boston, MA 02115, USA.
| | - Juliane G Bogner-Strauss
- Institute of Biochemistry, Graz University of Technology, Humboldtstrasse 46/III, 8010 Graz, Austria; BioTechMed-Graz, Graz, Austria.
| |
Collapse
|
13
|
Nagai H, Masuda A, Toya Y, Matsuda F, Shimizu H. Metabolic engineering of mevalonate-producing Escherichia coli strains based on thermodynamic analysis. Metab Eng 2018; 47:1-9. [DOI: 10.1016/j.ymben.2018.02.012] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2017] [Revised: 12/07/2017] [Accepted: 02/25/2018] [Indexed: 01/07/2023]
|
14
|
Golubeva LI, Shupletsov MS, Mashko SV. Metabolic Flux Analysis Using 13C Isotopes (13C-MFA). 1. Experimental Basis of the Method and the Present State of Investigations. APPL BIOCHEM MICRO+ 2018. [DOI: 10.1134/s0003683817070031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
15
|
Cheah YE, Young JD. Isotopically nonstationary metabolic flux analysis (INST-MFA): putting theory into practice. Curr Opin Biotechnol 2018. [PMID: 29522915 DOI: 10.1016/j.copbio.2018.02.013] [Citation(s) in RCA: 48] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Typically, 13C flux analysis relies on assumptions of both metabolic and isotopic steady state. If metabolism is steady but isotope labeling is not allowed to fully equilibrate, isotopically nonstationary metabolic flux analysis (INST-MFA) can be used to estimate fluxes. This requires solution of differential equations that describe the time-dependent labeling of network metabolites, while iteratively adjusting the flux and pool size parameters to match the transient labeling measurements. INST-MFA holds a number of unique advantages over approaches that rely solely upon steady-state isotope enrichments. First, INST-MFA can be applied to estimate fluxes in autotrophic systems, which consume only single-carbon substrates. Second, INST-MFA is ideally suited to systems that label slowly due to the presence of large intermediate pools or pathway bottlenecks. Finally, INST-MFA provides increased measurement sensitivity to estimate reversible exchange fluxes and metabolite pool sizes, which represents a potential framework for integrating metabolomic analysis with 13C flux analysis. This review highlights the unique capabilities of INST-MFA, describes newly available software tools that automate INST-MFA calculations, presents several practical examples of recent INST-MFA applications, and discusses the technical challenges that lie ahead.
Collapse
Affiliation(s)
- Yi Ern Cheah
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, PMB 351604, Nashville, TN 37235-1604, USA
| | - Jamey D Young
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, PMB 351604, Nashville, TN 37235-1604, USA; Department of Molecular Physiology and Biophysics, Vanderbilt University, PMB 351604, Nashville, TN 37235-1604, USA.
| |
Collapse
|
16
|
Badur MG, Metallo CM. Reverse engineering the cancer metabolic network using flux analysis to understand drivers of human disease. Metab Eng 2018; 45:95-108. [PMID: 29199104 PMCID: PMC5927620 DOI: 10.1016/j.ymben.2017.11.013] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2017] [Revised: 10/11/2017] [Accepted: 11/29/2017] [Indexed: 12/16/2022]
Abstract
Metabolic dysfunction has reemerged as an essential hallmark of tumorigenesis, and metabolic phenotypes are increasingly being integrated into pre-clinical models of disease. The complexity of these metabolic networks requires systems-level interrogation, and metabolic flux analysis (MFA) with stable isotope tracing present a suitable conceptual framework for such systems. Here we review efforts to elucidate mechanisms through which metabolism influences tumor growth and survival, with an emphasis on applications using stable isotope tracing and MFA. Through these approaches researchers can now quantify pathway fluxes in various in vitro and in vivo contexts to provide mechanistic insights at molecular and physiological scales respectively. Knowledge and discoveries in cancer models are paving the way toward applications in other biological contexts and disease models. In turn, MFA approaches will increasingly help to uncover new therapeutic opportunities that enhance human health.
Collapse
Affiliation(s)
- Mehmet G Badur
- Department of Bioengineering, University of California, San Diego, La Jolla, USA
| | - Christian M Metallo
- Department of Bioengineering, University of California, San Diego, La Jolla, USA; Moores Cancer Center, University of California, San Diego, La Jolla, USA; Diabetes and Endocrinology Research Center, University of California, San Diego, La Jolla, USA; Institute of Engineering in Medicine, University of California, San Diego, La Jolla, USA.
| |
Collapse
|
17
|
Metabolite concentrations, fluxes and free energies imply efficient enzyme usage. Nat Chem Biol 2016; 12:482-9. [PMID: 27159581 PMCID: PMC4912430 DOI: 10.1038/nchembio.2077] [Citation(s) in RCA: 265] [Impact Index Per Article: 33.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2015] [Accepted: 03/10/2016] [Indexed: 01/11/2023]
Abstract
In metabolism, available free energy is limited and must be divided across pathway steps to maintain ΔG negative throughout. For each reaction, ΔG is log-proportional both to a concentration ratio (reaction quotient-to-equilibrium constant) and to a flux ratio (backward-to-forward flux). Here we use isotope labeling to measure absolute metabolite concentrations and fluxes in Escherichia coli, yeast, and a mammalian cell line. We then integrate this information to obtain a unified set of concentrations and ΔG for each organism. In glycolysis, we find that free energy is partitioned so as to mitigate unproductive backward fluxes associated with ΔG near zero. Across metabolism, we observe that absolute metabolite concentrations and ΔG are substantially conserved, and that most substrate (but not inhibitor) concentrations exceed the associated enzyme binding site affinity. The observed conservation of metabolite concentrations is consistent with an evolutionary drive to utilize enzymes efficiently given thermodynamic and osmotic constraints.
Collapse
|
18
|
Fluxes through plant metabolic networks: measurements, predictions, insights and challenges. Biochem J 2015; 465:27-38. [PMID: 25631681 DOI: 10.1042/bj20140984] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Although the flows of material through metabolic networks are central to cell function, they are not easy to measure other than at the level of inputs and outputs. This is particularly true in plant cells, where the network spans multiple subcellular compartments and where the network may function either heterotrophically or photoautotrophically. For many years, kinetic modelling of pathways provided the only method for describing the operation of fragments of the network. However, more recently, it has become possible to map the fluxes in central carbon metabolism using the stable isotope labelling techniques of metabolic flux analysis (MFA), and to predict intracellular fluxes using constraints-based modelling procedures such as flux balance analysis (FBA). These approaches were originally developed for the analysis of microbial metabolism, but over the last decade, they have been adapted for the more demanding analysis of plant metabolic networks. Here, the principal features of MFA and FBA as applied to plants are outlined, followed by a discussion of the insights that have been gained into plant metabolic networks through the application of these time-consuming and non-trivial methods. The discussion focuses on how a system-wide view of plant metabolism has increased our understanding of network structure, metabolic perturbations and the provision of reducing power and energy for cell function. Current methodological challenges that limit the scope of plant MFA are discussed and particular emphasis is placed on the importance of developing methods for cell-specific MFA.
Collapse
|
19
|
Lien SK, Niedenführ S, Sletta H, Nöh K, Bruheim P. Fluxome study of Pseudomonas fluorescens reveals major reorganisation of carbon flux through central metabolic pathways in response to inactivation of the anti-sigma factor MucA. BMC SYSTEMS BIOLOGY 2015; 9:6. [PMID: 25889900 PMCID: PMC4351692 DOI: 10.1186/s12918-015-0148-0] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/28/2014] [Accepted: 01/27/2015] [Indexed: 11/25/2022]
Abstract
Background The bacterium Pseudomonas fluorescens switches to an alginate-producing phenotype when the pleiotropic anti-sigma factor MucA is inactivated. The inactivation is accompanied by an increased biomass yield on carbon sources when grown under nitrogen-limited chemostat conditions. A previous metabolome study showed significant changes in the intracellular metabolite concentrations, especially of the nucleotides, in mucA deletion mutants compared to the wild-type. In this study, the P. fluorescens SBW25 wild-type and an alginate non-producing mucA- ΔalgC double-knockout mutant are investigated through model-based 13C-metabolic flux analysis (13C-MFA) to explore the physiological consequences of MucA inactivation at the metabolic flux level. Intracellular metabolite extracts from three carbon labelling experiments using fructose as the sole carbon source are analysed for 13C-label incorporation in primary metabolites by gas and liquid chromatography tandem mass spectrometry. Results From mass isotopomer distribution datasets, absolute intracellular metabolic reaction rates for the wild type and the mutant are determined, revealing extensive reorganisation of carbon flux through central metabolic pathways in response to MucA inactivation. The carbon flux through the Entner-Doudoroff pathway was reduced in the mucA- ΔalgC mutant, while flux through the pentose phosphate pathway was increased. Our findings also indicated flexibility of the anaplerotic reactions through down-regulation of the pyruvate shunt in the mucA- ΔalgC mutant and up-regulation of the glyoxylate shunt. Conclusions Absolute metabolic fluxes and metabolite levels give detailed, integrated insight into the physiology of this industrially, medically and agriculturally important bacterial species and suggest that the most efficient way of using a mucA- mutant as a cell factory for alginate production would be to use non-growing conditions and nitrogen deprivation. Electronic supplementary material The online version of this article (doi:10.1186/s12918-015-0148-0) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Stina K Lien
- Department of Biotechnology, Norwegian University of Science and Technology, Sem Sælands vei 6/8, N-7491, Trondheim, Norway.
| | - Sebastian Niedenführ
- Institute of Bio- and Geosciences IBG-1: Biotechnology, Forschungszentrum Jülich, D-52425, Jülich, Germany.
| | - Håvard Sletta
- Department of Bioprocess technology, SINTEF Materials and Chemistry, Sem Sælands vei 2a, N-7465, Trondheim, Norway.
| | - Katharina Nöh
- Institute of Bio- and Geosciences IBG-1: Biotechnology, Forschungszentrum Jülich, D-52425, Jülich, Germany.
| | - Per Bruheim
- Department of Biotechnology, Norwegian University of Science and Technology, Sem Sælands vei 6/8, N-7491, Trondheim, Norway.
| |
Collapse
|
20
|
Heise R, Fernie AR, Stitt M, Nikoloski Z. Pool size measurements facilitate the determination of fluxes at branching points in non-stationary metabolic flux analysis: the case of Arabidopsis thaliana. FRONTIERS IN PLANT SCIENCE 2015; 6:386. [PMID: 26082786 PMCID: PMC4451360 DOI: 10.3389/fpls.2015.00386] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/17/2014] [Accepted: 05/14/2015] [Indexed: 05/08/2023]
Abstract
Pool size measurements are important for the estimation of absolute intracellular fluxes in particular scenarios based on data from heavy carbon isotope experiments. Recently, steady-state fluxes estimates were obtained for central carbon metabolism in an intact illuminated rosette of Arabidopsis thaliana grown photoautotrophically (Szecowka et al., 2013; Heise et al., 2014). Fluxes were estimated therein by integrating mass-spectrometric data of the dynamics of the unlabeled metabolic fraction, data on metabolic pool sizes, partitioning of metabolic pools between cellular compartments and estimates of photosynthetically inactive pools, with a simplified model of plant central carbon metabolism. However, the fluxes were determined by treating the pool sizes as fixed parameters. Here we investigated whether and, if so, to what extent the treatment of pool sizes as parameters to be optimized in three scenarios may affect the flux estimates. The results are discussed in terms of benchmark values for canonical pathways and reactions, including starch and sucrose synthesis as well as the ribulose-1,5-bisphosphate carboxylation and oxygenation reactions. In addition, we discuss pathways emerging from a divergent branch point for which pool sizes are required for flux estimation, irrespective of the computational approach used for the simulation of the observable labeling pattern. Therefore, our findings indicate the necessity for development of techniques for accurate pool size measurements to improve the quality of flux estimates from non-stationary flux estimates in intact plant cells in the absence of alternative flux measurements.
Collapse
Affiliation(s)
- Robert Heise
- Systems Biology and Mathematical Modeling Group, Max Planck Institute of Molecular Plant PhysiologyPotsdam, Germany
| | - Alisdair R. Fernie
- Central Metabolism Group, Max Planck Institute of Molecular Plant PhysiologyPotsdam, Germany
| | - Mark Stitt
- System Regulation Group, Max Planck Institute of Molecular Plant PhysiologyPotsdam, Germany
| | - Zoran Nikoloski
- Systems Biology and Mathematical Modeling Group, Max Planck Institute of Molecular Plant PhysiologyPotsdam, Germany
- *Correspondence: Zoran Nikoloski, Systems Biology and Mathematical Modeling Group, Max Planck Institute of Molecular Plant Physiology, Am Mühlenberg 1, 14476 Potsdam, Germany
| |
Collapse
|
21
|
Millard P, Massou S, Portais JC, Létisse F. Isotopic Studies of Metabolic Systems by Mass Spectrometry: Using Pascal’s Triangle To Produce Biological Standards with Fully Controlled Labeling Patterns. Anal Chem 2014; 86:10288-95. [DOI: 10.1021/ac502490g] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Pierre Millard
- Université de Toulouse; INSA, UPS, INP; LISBP, 135 Avenue de Rangueil, 31077 Toulouse, France
- INRA, UMR792, Ingénierie des Systèmes Biologiques et des Procédés, 31400 Toulouse, France
- CNRS, UMR5504, 31400 Toulouse, France
| | - Stéphane Massou
- Université de Toulouse; INSA, UPS, INP; LISBP, 135 Avenue de Rangueil, 31077 Toulouse, France
- INRA, UMR792, Ingénierie des Systèmes Biologiques et des Procédés, 31400 Toulouse, France
- CNRS, UMR5504, 31400 Toulouse, France
| | - Jean-Charles Portais
- Université de Toulouse; INSA, UPS, INP; LISBP, 135 Avenue de Rangueil, 31077 Toulouse, France
- INRA, UMR792, Ingénierie des Systèmes Biologiques et des Procédés, 31400 Toulouse, France
- CNRS, UMR5504, 31400 Toulouse, France
| | - Fabien Létisse
- Université de Toulouse; INSA, UPS, INP; LISBP, 135 Avenue de Rangueil, 31077 Toulouse, France
- INRA, UMR792, Ingénierie des Systèmes Biologiques et des Procédés, 31400 Toulouse, France
- CNRS, UMR5504, 31400 Toulouse, France
| |
Collapse
|
22
|
Wiechert W, Nöh K. Isotopically non-stationary metabolic flux analysis: complex yet highly informative. Curr Opin Biotechnol 2013; 24:979-86. [DOI: 10.1016/j.copbio.2013.03.024] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2012] [Revised: 03/28/2013] [Accepted: 03/30/2013] [Indexed: 12/16/2022]
|
23
|
Través PG, de Atauri P, Marín S, Pimentel-Santillana M, Rodríguez-Prados JC, Marín de Mas I, Selivanov VA, Martín-Sanz P, Boscá L, Cascante M. Relevance of the MEK/ERK signaling pathway in the metabolism of activated macrophages: a metabolomic approach. JOURNAL OF IMMUNOLOGY (BALTIMORE, MD. : 1950) 2012; 188:1402-10. [PMID: 22190182 DOI: 10.4049/jimmunol.1101781] [Citation(s) in RCA: 57] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The activation of immune cells in response to a pathogen involves a succession of signaling events leading to gene and protein expression, which requires metabolic changes to match the energy demands. The metabolic profile associated with the MAPK cascade (ERK1/2, p38, and JNK) in macrophages was studied, and the effect of its inhibition on the specific metabolic pattern of LPS stimulation was characterized. A [1,2-[(13)C](2)]glucose tracer-based metabolomic approach was used to examine the metabolic flux distribution in these cells after MEK/ERK inhibition. Bioinformatic tools were used to analyze changes in mass isotopomer distribution and changes in glucose and glutamine consumption and lactate production in basal and LPS-stimulated conditions in the presence and absence of the selective inhibitor of the MEK/ERK cascade, PD325901. Results showed that PD325901-mediated ERK1/2 inhibition significantly decreased glucose consumption and lactate production but did not affect glutamine consumption. These changes were accompanied by a decrease in the glycolytic flux, consistent with the observed decrease in fructose-2,6-bisphosphate concentration. The oxidative and nonoxidative pentose phosphate pathways and the ratio between them also decreased. However, tricarboxylic acid cycle flux did not change significantly. LPS activation led to the opposite responses, although all of these were suppressed by PD325901. However, LPS also induced a small decrease in pentose phosphate pathway fluxes and an increase in glutamine consumption that were not affected by PD325901. We concluded that inhibition of the MEK/ERK cascade interferes with central metabolism, and this cross-talk between signal transduction and metabolism also occurs in the presence of LPS.
Collapse
Affiliation(s)
- Paqui G Través
- Instituto de Investigaciones Biomédicas Alberto Sols, Consejo Superior de Investigaciones Científicas-Universidad Autónoma de Madrid, 28029 Madrid, Spain
| | | | | | | | | | | | | | | | | | | |
Collapse
|
24
|
The benefits of being transient: isotope-based metabolic flux analysis at the short time scale. Appl Microbiol Biotechnol 2011; 91:1247-65. [DOI: 10.1007/s00253-011-3390-4] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2011] [Revised: 05/15/2011] [Accepted: 05/16/2011] [Indexed: 12/24/2022]
|
25
|
An in vivo data-driven framework for classification and quantification of enzyme kinetics and determination of apparent thermodynamic data. Metab Eng 2011; 13:294-306. [DOI: 10.1016/j.ymben.2011.02.005] [Citation(s) in RCA: 80] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2010] [Revised: 01/10/2011] [Accepted: 02/15/2011] [Indexed: 01/21/2023]
|
26
|
Thermodynamic calculations for biochemical transport and reaction processes in metabolic networks. Biophys J 2011; 99:3139-44. [PMID: 21081060 DOI: 10.1016/j.bpj.2010.09.043] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2010] [Revised: 08/17/2010] [Accepted: 09/09/2010] [Indexed: 11/22/2022] Open
Abstract
Thermodynamic analysis of metabolic networks has recently generated increasing interest for its ability to add constraints on metabolic network operation, and to combine metabolic fluxes and metabolite measurements in a mechanistic manner. Concepts for the calculation of the change in Gibbs energy of biochemical reactions have long been established. However, a concept for incorporation of cross-membrane transport in these calculations is still missing, although the theory for calculating thermodynamic properties of transport processes is long known. Here, we have developed two equivalent equations to calculate the change in Gibbs energy of combined transport and reaction processes based on two different ways of treating biochemical thermodynamics. We illustrate the need for these equations by showing that in some cases there is a significant difference between the proposed correct calculation and using an approximative method. With the developed equations, thermodynamic analysis of metabolic networks spanning over multiple physical compartments can now be correctly described.
Collapse
|
27
|
Integrated stoichiometric, thermodynamic and kinetic modelling of steady state metabolism. J Theor Biol 2010; 264:683-92. [PMID: 20230840 DOI: 10.1016/j.jtbi.2010.02.044] [Citation(s) in RCA: 36] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2009] [Revised: 01/15/2010] [Accepted: 02/24/2010] [Indexed: 01/05/2023]
Abstract
The quantitative analysis of biochemical reactions and metabolites is at frontier of biological sciences. The recent availability of high-throughput technology data sets in biology has paved the way for new modelling approaches at various levels of complexity including the metabolome of a cell or an organism. Understanding the metabolism of a single cell and multi-cell organism will provide the knowledge for the rational design of growth conditions to produce commercially valuable reagents in biotechnology. Here, we demonstrate how equations representing steady state mass conservation, energy conservation, the second law of thermodynamics, and reversible enzyme kinetics can be formulated as a single system of linear equalities and inequalities, in addition to linear equalities on exponential variables. Even though the feasible set is non-convex, the reformulation is exact and amenable to large-scale numerical analysis, a prerequisite for computationally feasible genome scale modelling. Integrating flux, concentration and kinetic variables in a unified constraint-based formulation is aimed at increasing the quantitative predictive capacity of flux balance analysis. Incorporation of experimental and theoretical bounds on thermodynamic and kinetic variables ensures that the predicted steady state fluxes are both thermodynamically and biochemically feasible. The resulting in silico predictions are tested against fluxomic data for central metabolism in Escherichia coli and compare favourably with in silico prediction by flux balance analysis.
Collapse
|
28
|
Dauner M. From fluxes and isotope labeling patterns towards in silico cells. Curr Opin Biotechnol 2010; 21:55-62. [DOI: 10.1016/j.copbio.2010.01.014] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2009] [Revised: 01/23/2010] [Accepted: 01/31/2010] [Indexed: 10/19/2022]
|
29
|
Kleijn RJ, Buescher JM, Le Chat L, Jules M, Aymerich S, Sauer U. Metabolic fluxes during strong carbon catabolite repression by malate in Bacillus subtilis. J Biol Chem 2009; 285:1587-96. [PMID: 19917605 DOI: 10.1074/jbc.m109.061747] [Citation(s) in RCA: 86] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
Commonly glucose is considered to be the only preferred substrate in Bacillus subtilis whose presence represses utilization of other alternative substrates. Because recent data indicate that malate might be an exception, we quantify here the carbon source utilization hierarchy. Based on physiology and transcriptional data during co-utilization experiments with eight carbon substrates, we demonstrate that malate is a second preferred carbon source for B. subtilis, which is rapidly co-utilized with glucose and strongly represses the uptake of alternative substrates. From the different hierarchy and degree of catabolite repression exerted by glucose and malate, we conclude that both substrates might act through different molecular mechanisms. To obtain a quantitative and functional network view of how malate is (co)metabolized, we developed a novel approach to metabolic flux analysis that avoids error-prone, intuitive, and ad hoc decisions on (13)C rearrangements. In particular, we developed a rigorous approach for deriving reaction reversibilities by combining in vivo intracellular metabolite concentrations with a thermodynamic feasibility analysis. The thus-obtained analytical model of metabolism was then used for network-wide isotopologue balancing to estimate the intracellular fluxes. These (13)C-flux data revealed an extraordinarily high malate influx that is primarily catabolized via the gluconeogenic reactions and toward overflow metabolism. Furthermore, a considerable NADPH-producing malic enzyme flux is required to supply the biosynthetically required NADPH in the presence of malate. Co-utilization of glucose and malate resulted in a synergistic decrease of the respiratory tricarboxylic acid cycle flux.
Collapse
Affiliation(s)
- Roelco J Kleijn
- Institute of Molecular System Biology, ETH Zürich, CH-8093 Zürich, Switzerland
| | | | | | | | | | | |
Collapse
|
30
|
Allen DK, Libourel IGL, Shachar-Hill Y. Metabolic flux analysis in plants: coping with complexity. PLANT, CELL & ENVIRONMENT 2009; 32:1241-57. [PMID: 19422611 DOI: 10.1111/j.1365-3040.2009.01992.x] [Citation(s) in RCA: 95] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Theory and experience in metabolic engineering both show that metabolism operates at the network level. In plants, this complexity is compounded by a high degree of compartmentation and the synthesis of a very wide array of secondary metabolic products. A further challenge to understanding and predicting plant metabolic function is posed by our ignorance about the structure of metabolic networks even in well-studied systems. Metabolic flux analysis (MFA) provides tools to measure and model the functioning of metabolism, and is making significant contributions to coping with their complexity. This review gives an overview of different MFA approaches, the measurements required to implement them and the information they yield. The application of MFA methods to plant systems is then illustrated by several examples from the recent literature. Next, the challenges that plant metabolism poses for MFA are discussed together with ways that these can be addressed. Lastly, new developments in MFA are described that can be expected to improve the range and reliability of plant MFA in the coming years.
Collapse
Affiliation(s)
- Doug K Allen
- Michigan State University, Plant Biology Department, East Lansing, MI 48824, USA.
| | | | | |
Collapse
|
31
|
|