1
|
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
|
2
|
Moiz B, Li A, Padmanabhan S, Sriram G, Clyne AM. Isotope-Assisted Metabolic Flux Analysis: A Powerful Technique to Gain New Insights into the Human Metabolome in Health and Disease. Metabolites 2022; 12:1066. [PMID: 36355149 PMCID: PMC9694183 DOI: 10.3390/metabo12111066] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Revised: 10/25/2022] [Accepted: 10/27/2022] [Indexed: 04/28/2024] Open
Abstract
Cell metabolism represents the coordinated changes in genes, proteins, and metabolites that occur in health and disease. The metabolic fluxome, which includes both intracellular and extracellular metabolic reaction rates (fluxes), therefore provides a powerful, integrated description of cellular phenotype. However, intracellular fluxes cannot be directly measured. Instead, flux quantification requires sophisticated mathematical and computational analysis of data from isotope labeling experiments. In this review, we describe isotope-assisted metabolic flux analysis (iMFA), a rigorous computational approach to fluxome quantification that integrates metabolic network models and experimental data to generate quantitative metabolic flux maps. We highlight practical considerations for implementing iMFA in mammalian models, as well as iMFA applications in in vitro and in vivo studies of physiology and disease. Finally, we identify promising new frontiers in iMFA which may enable us to fully unlock the potential of iMFA in biomedical research.
Collapse
Affiliation(s)
- Bilal Moiz
- Department of Bioengineering, University of Maryland, College Park, MD 20742, USA
| | - Andrew Li
- Department of Bioengineering, University of Maryland, College Park, MD 20742, USA
| | - Surya Padmanabhan
- 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
|
3
|
Bednarski TK, Rahim M, Young JD. In vivo 2H/ 13C flux analysis in metabolism research. Curr Opin Biotechnol 2021; 71:1-8. [PMID: 34048994 DOI: 10.1016/j.copbio.2021.04.005] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 03/29/2021] [Accepted: 04/26/2021] [Indexed: 12/12/2022]
Abstract
Identifying the factors and mechanisms that regulate metabolism under normal and diseased states requires methods to quantify metabolic fluxes of live tissues within their physiological milieu. A number of recent developments have expanded the reach and depth of isotope-based in vivo flux analysis, which have in turn challenged existing dogmas in metabolism research. First, minimally invasive techniques of intravenous isotope infusion and sampling have advanced in vivo metabolic tracer studies in animal models and human subjects. Second, recent breakthroughs in analytical instrumentation have expanded the scope of isotope labeling measurements and reduced sample volume requirements. Third, innovative modeling approaches and publicly available software tools have facilitated rigorous analysis of sophisticated experimental designs involving multiple tracers and expansive metabolomics datasets. These developments have enabled comprehensive in vivo quantification of metabolic fluxes in specific tissues and have set the stage for integrated multi-tissue flux assays.
Collapse
Affiliation(s)
- Tomasz K Bednarski
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN, USA
| | - Mohsin Rahim
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN, USA
| | - Jamey D Young
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, TN, USA; Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN, USA.
| |
Collapse
|
4
|
Volkova S, Matos MRA, Mattanovich M, Marín de Mas I. Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis. Metabolites 2020; 10:E303. [PMID: 32722118 PMCID: PMC7465778 DOI: 10.3390/metabo10080303] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 07/08/2020] [Accepted: 07/22/2020] [Indexed: 01/05/2023] Open
Abstract
Metabolic networks are regulated to ensure the dynamic adaptation of biochemical reaction fluxes to maintain cell homeostasis and optimal metabolic fitness in response to endogenous and exogenous perturbations. To this end, metabolism is tightly controlled by dynamic and intricate regulatory mechanisms involving allostery, enzyme abundance and post-translational modifications. The study of the molecular entities involved in these complex mechanisms has been boosted by the advent of high-throughput technologies. The so-called omics enable the quantification of the different molecular entities at different system layers, connecting the genotype with the phenotype. Therefore, the study of the overall behavior of a metabolic network and the omics data integration and analysis must be approached from a holistic perspective. Due to the close relationship between metabolism and cellular phenotype, metabolic modelling has emerged as a valuable tool to decipher the underlying mechanisms governing cell phenotype. Constraint-based modelling and kinetic modelling are among the most widely used methods to study cell metabolism at different scales, ranging from cells to tissues and organisms. These approaches enable integrating metabolomic data, among others, to enhance model predictive capabilities. In this review, we describe the current state of the art in metabolic modelling and discuss future perspectives and current challenges in the field.
Collapse
Affiliation(s)
| | | | | | - Igor Marín de Mas
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark; (S.V.); (M.R.A.M.); (M.M.)
| |
Collapse
|
5
|
Mairinger T, Hann S. Determination of Isotopologue and Tandem Mass Isotopologue Ratios Using Gas Chromatography Chemical Ionization Time of Flight Mass Spectrometry - Methodology and Uncertainty of Measurement. Methods Mol Biol 2020; 2088:1-16. [PMID: 31893367 DOI: 10.1007/978-1-0716-0159-4_1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
The accurate and precise analysis of isotopologue and tandem mass isotopologue ratios in heavy stable isotope labeling experiments is a critical part of assessing absolute intracellular metabolic fluxes. Resulting from feeding the organism of interest with a specifically isotope-labeled substrate, the principal characteristics of these labeling experiments are the metabolites' non-naturally distributed isotopologue patterns. For the purpose of inferring metabolic rates by maximizing the fit between a priori simulated and experimentally obtained labeling patterns, 13C is the preferred stable isotope of use.The analysis of the obtained labeling patterns can be approached by different mass spectrometric approaches. Gas chromatography (GC) features broad metabolite coverage and excellent separation efficiency of biologically relevant isomers. These advantages compensate for laborious derivatization steps and the resulting need for interference correction for natural abundant isotopes.Here, we describe a workflow based on GC-high resolution mass spectrometry with chemical ionization for the analysis of carbon-isotopologue distributions and some positional labeling information of primary metabolites. To study the associated measurement uncertainty of the resulting 13C labeling patterns, guidance to uncertainty estimation according to the EURACHEM guidelines with Monte-Carlo simulation is provided.
Collapse
Affiliation(s)
- Teresa Mairinger
- Department of Chemistry, University of Natural Resources and Life Sciences-BOKU Vienna, Vienna, Austria.
- Department of Environmental Chemistry, Eawag: Swiss Federal Institute of Aquatic Science and Technology, Duebendorf, Switzerland.
| | - Stephen Hann
- Department of Chemistry, University of Natural Resources and Life Sciences-BOKU Vienna, Vienna, Austria
- Austrian Center for Industrial Biotechnology, Vienna, Austria
| |
Collapse
|
6
|
Kappelmann J, Beyß M, Nöh K, Noack S. Separation of 13C- and 15N-Isotopologues of Amino Acids with a Primary Amine without Mass Resolution by Means of O-Phthalaldehyde Derivatization and Collision Induced Dissociation. Anal Chem 2019; 91:13407-13417. [PMID: 31577133 DOI: 10.1021/acs.analchem.9b01788] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Computational and experimental advances of recent years have culminated in establishing 13C-Metabolic Flux Analysis (13C-MFA) as a routine methodology to unravel the fluxome. As the acronym suggests, 13C-MFA has relied on the relative abundance of 13C-isotopes in metabolites for flux inference, most commonly measured by mass spectrometry. In this manuscript we expand the scope of labeling measurements to the case of simultaneous 13C- and 15N-labeling of amino acids. Analytically, the separation of isotopologues of this metabolite class can only be achieved at resolving power beyond 65,000. In this manuscript we harvest an overlooked property of the collision induced dissociation of amino acid adducts to discern 13C- and 15N- isotopologues of amino acids with a primary amine without separating them in the m/z domain.
Collapse
Affiliation(s)
- Jannick Kappelmann
- Institute of Bio- and Geosciences I, IBG-1: Biotechnology , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany
| | - Martin Beyß
- Institute of Bio- and Geosciences I, IBG-1: Biotechnology , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany
| | - Katharina Nöh
- Institute of Bio- and Geosciences I, IBG-1: Biotechnology , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany
| | - Stephan Noack
- Institute of Bio- and Geosciences I, IBG-1: Biotechnology , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany.,Bioeconomy Science Center (BioSC) , Forschungszentrum Jülich GmbH , 52425 Jülich , Germany
| |
Collapse
|
7
|
Lularevic M, Racher AJ, Jaques C, Kiparissides A. Improving the accuracy of flux balance analysis through the implementation of carbon availability constraints for intracellular reactions. Biotechnol Bioeng 2019; 116:2339-2352. [PMID: 31112296 DOI: 10.1002/bit.27025] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2019] [Revised: 03/29/2019] [Accepted: 04/25/2019] [Indexed: 11/07/2022]
Abstract
Constraint-based modeling methods, such as Flux Balance Analysis (FBA), have been extensively used to decipher complex, information rich -omics datasets to elicit system-wide behavioral patterns of cellular metabolism. FBA has been successfully used to gain insight in a wide range of applications, such as range of substrate utilization, product yields and to design metabolic engineering strategies to improve bioprocess performance. A well-known challenge associated with large genome-scale metabolic networks is that they result in underdetermined problem formulations. Consequently, rather than unique solutions, FBA and related methods examine ranges of reaction flux values that are consistent with the studied physiological conditions. The wider the reported flux ranges, the higher the uncertainty in the determination of basic reaction properties, limiting interpretability of and confidence in the results. Herein, we propose a new, computationally efficient approach that refines flux range predictions by constraining reaction fluxes on the basis of the elemental balance of carbon. We compared carbon constraint FBA (ccFBA) against experimentally-measured intracellular fluxes using the latest CHO GEM (iCHO1766) and were able to substantially improve the accuracy of predicted flux values compared with FBA. ccFBA can be used as a stand-alone method but is also compatible with and complimentary to other constraint-based approaches.
Collapse
Affiliation(s)
- Maximilian Lularevic
- Department of Biochemical Engineering, University College London, London, UK
- Research and Technology, Lonza Biologics PLC, Slough, UK
| | | | - Colin Jaques
- Research and Technology, Lonza Biologics PLC, Slough, UK
| | | |
Collapse
|
8
|
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
|
9
|
Schwaiger-Haber M, Hermann G, El Abiead Y, Rampler E, Wernisch S, Sas K, Pennathur S, Koellensperger G. Proposing a validation scheme for 13C metabolite tracer studies in high-resolution mass spectrometry. Anal Bioanal Chem 2019; 411:3103-3113. [PMID: 30972471 PMCID: PMC6526147 DOI: 10.1007/s00216-019-01773-7] [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: 10/12/2018] [Revised: 02/25/2019] [Accepted: 03/08/2019] [Indexed: 12/31/2022]
Abstract
13C metabolite tracer and metabolic flux analyses require upfront experimental planning and validation tools. Here, we present a validation scheme including a comparison of different LC methods that allow for customization of analytical strategies for tracer studies with regard to the targeted metabolites. As the measurement of significant changes in labeling patterns depends on the spectral accuracy, we investigate this aspect comprehensively for high-resolution orbitrap mass spectrometry combined with reversed-phase chromatography, hydrophilic interaction liquid chromatography, or anion-exchange chromatography. Moreover, we propose a quality control protocol based on (1) a metabolite containing selenium to assess the instrument performance and on (2) in vivo synthesized isotopically enriched Pichia pastoris to validate the accuracy of carbon isotopologue distributions (CIDs), in this case considering each isotopologue of a targeted metabolite panel. Finally, validation involved a thorough assessment of procedural blanks and matrix interferences. We compared the analytical figures of merit regarding CID determination for over 40 metabolites between the three methods. Excellent precisions of less than 1% and trueness bias as small as 0.01-1% were found for the majority of compounds, whereas the CID determination of a small fraction was affected by contaminants. For most compounds, changes of labeling pattern as low as 1% could be measured. Graphical abstract.
Collapse
Affiliation(s)
- Michaela Schwaiger-Haber
- Department of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090, Vienna, Austria.,Division of Nephrology, Department of Internal Medicine, University of Michigan, 1000 Wall St., Ann Arbor, MI, 48105, USA
| | - Gerrit Hermann
- Department of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090, Vienna, Austria.,ISOtopic solutions, Waehringer Str. 38, 1090, Vienna, Austria
| | - Yasin El Abiead
- Department of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090, Vienna, Austria.,Vienna Metabolomics Center (VIME), University of Vienna, Althanstraße 14, 1090, Vienna, Austria.,Chemistry Meets Microbiology, Althanstraße 14, 1090, Vienna, Austria
| | - Evelyn Rampler
- Department of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090, Vienna, Austria.,Vienna Metabolomics Center (VIME), University of Vienna, Althanstraße 14, 1090, Vienna, Austria.,Chemistry Meets Microbiology, Althanstraße 14, 1090, Vienna, Austria
| | - Stefanie Wernisch
- Division of Nephrology, Department of Internal Medicine, University of Michigan, 1000 Wall St., Ann Arbor, MI, 48105, USA
| | - Kelli Sas
- Division of Nephrology, Department of Internal Medicine, University of Michigan, 1000 Wall St., Ann Arbor, MI, 48105, USA
| | - Subramaniam Pennathur
- Division of Nephrology, Department of Internal Medicine, University of Michigan, 1000 Wall St., Ann Arbor, MI, 48105, USA.,Department of Molecular and Integrative Physiology, University of Michigan, 1000 Wall St, Ann Arbor, MI, 48105, USA
| | - Gunda Koellensperger
- Department of Analytical Chemistry, Faculty of Chemistry, University of Vienna, Waehringer Str. 38, 1090, Vienna, Austria. .,Vienna Metabolomics Center (VIME), University of Vienna, Althanstraße 14, 1090, Vienna, Austria. .,Chemistry Meets Microbiology, Althanstraße 14, 1090, Vienna, Austria.
| |
Collapse
|
10
|
Millard P, Delépine B, Guionnet M, Heuillet M, Bellvert F, Létisse F. IsoCor: isotope correction for high-resolution MS labeling experiments. Bioinformatics 2019; 35:4484-4487. [DOI: 10.1093/bioinformatics/btz209] [Citation(s) in RCA: 60] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2018] [Revised: 02/22/2019] [Accepted: 03/20/2019] [Indexed: 11/14/2022] Open
Abstract
Abstract
Summary
Mass spectrometry (MS) is widely used for isotopic studies of metabolism and other (bio)chemical processes. Quantitative applications in systems and synthetic biology require to correct the raw MS data for the contribution of naturally occurring isotopes. Several tools are available to correct low-resolution MS data, and recent developments made substantial improvements by introducing resolution-dependent correction methods, hence opening the way to the correction of high-resolution MS (HRMS) data. Nevertheless, current HRMS correction methods partly fail to determine which isotopic species are resolved from the tracer isotopologues and should thus be corrected. We present an updated version of our isotope correction software (IsoCor) with a novel correction algorithm which ensures to accurately exploit any chemical species with any isotopic tracer, at any MS resolution. IsoCor v2 also includes a novel graphical user interface for intuitive use by end-users and a command-line interface to streamline integration into existing pipelines.
Availability and implementation
IsoCor v2 is implemented in Python 3 and was tested on Windows, Unix and MacOS platforms. The source code and the documentation are freely distributed under GPL3 license at https://github.com/MetaSys-LISBP/IsoCor/ and https://isocor.readthedocs.io/.
Collapse
Affiliation(s)
- Pierre Millard
- LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse 31077, France
| | - Baudoin Delépine
- LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse 31077, France
| | - Matthieu Guionnet
- LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse 31077, France
- MetaToul-MetaboHUB, National Infrastructure of Metabolomics and Fluxomics, Toulouse, France
| | - Maud Heuillet
- LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse 31077, France
- MetaToul-MetaboHUB, National Infrastructure of Metabolomics and Fluxomics, Toulouse, France
| | - Floriant Bellvert
- LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse 31077, France
- MetaToul-MetaboHUB, National Infrastructure of Metabolomics and Fluxomics, Toulouse, France
| | - Fabien Létisse
- LISBP, Université de Toulouse, CNRS, INRA, INSA, Toulouse 31077, France
| |
Collapse
|
11
|
Correia DM, Sargo CR, Silva AJ, Santos ST, Giordano RC, Ferreira EC, Zangirolami TC, Ribeiro MPA, Rocha I. Mapping Salmonella typhimurium pathways using 13C metabolic flux analysis. Metab Eng 2019; 52:303-314. [PMID: 30529284 DOI: 10.1016/j.ymben.2018.11.011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2018] [Revised: 11/26/2018] [Accepted: 11/28/2018] [Indexed: 12/20/2022]
Abstract
In the last years, Salmonella has been extensively studied not only due to its importance as a pathogen, but also as a host to produce pharmaceutical compounds. However, the full exploitation of Salmonella as a platform for bioproduct delivery has been hampered by the lack of information about its metabolism. Genome-scale metabolic models can be valuable tools to delineate metabolic engineering strategies as long as they closely represent the actual metabolism of the target organism. In the present study, a 13C-MFA approach was applied to map the fluxes at the central carbon pathways of S. typhimurium LT2 growing at glucose-limited chemostat cultures. The experiments were carried out in a 2L bioreactor, using defined medium enriched with 20% 13C-labeled glucose. Metabolic flux distributions in central carbon pathways of S. typhimurium LT2 were estimated using OpenFLUX2 based on the labeling pattern of biomass protein hydrolysates together with biomass composition. The results suggested that pentose phosphate is used to catabolize glucose, with minor fluxes through glycolysis. In silico simulations, using Optflux and pFBA as simulation method, allowed to study the performance of the genome-scale metabolic model. In general, the accuracy of in silico simulations was improved by the superimposition of estimated intracellular fluxes to the existing genome-scale metabolic model, showing a better fitting to the experimental extracellular fluxes, whereas the intracellular fluxes of pentose phosphate and anaplerotic reactions were poorly described.
Collapse
Affiliation(s)
- Daniela M Correia
- Graduate Program of Chemical Engineering, Federal University of São Carlos, Rodovia Washington Luís, Km 235, São Carlos, SP 13565-905, Brazil
| | - Cintia R Sargo
- Graduate Program of Chemical Engineering, Federal University of São Carlos, Rodovia Washington Luís, Km 235, São Carlos, SP 13565-905, Brazil
| | - Adilson J Silva
- Graduate Program of Chemical Engineering, Federal University of São Carlos, Rodovia Washington Luís, Km 235, São Carlos, SP 13565-905, Brazil
| | - Sophia T Santos
- CEB-Centre of Biological Engineering, University of Minho, Campus De Gualtar, Braga 4710-057, Portugal
| | - Roberto C Giordano
- Graduate Program of Chemical Engineering, Federal University of São Carlos, Rodovia Washington Luís, Km 235, São Carlos, SP 13565-905, Brazil
| | - Eugénio C Ferreira
- CEB-Centre of Biological Engineering, University of Minho, Campus De Gualtar, Braga 4710-057, Portugal
| | - Teresa C Zangirolami
- Graduate Program of Chemical Engineering, Federal University of São Carlos, Rodovia Washington Luís, Km 235, São Carlos, SP 13565-905, Brazil
| | - Marcelo P A Ribeiro
- Graduate Program of Chemical Engineering, Federal University of São Carlos, Rodovia Washington Luís, Km 235, São Carlos, SP 13565-905, Brazil
| | - Isabel Rocha
- CEB-Centre of Biological Engineering, University of Minho, Campus De Gualtar, Braga 4710-057, Portugal; Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB-NOVA), Oeiras, Portugal.
| |
Collapse
|
12
|
Mairinger T, Sanderson J, Hann S. GC-QTOFMS with a low-energy electron ionization source for advancing isotopologue analysis in 13C-based metabolic flux analysis. Anal Bioanal Chem 2019; 411:1495-1502. [PMID: 30796486 DOI: 10.1007/s00216-019-01590-y] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2018] [Revised: 12/17/2018] [Accepted: 01/07/2019] [Indexed: 12/30/2022]
Abstract
For the study of different levels of (intra)cellular regulation and condition-dependent insight into metabolic activities, fluxomics experiments based on stable isotope tracer experiments using 13C have become a well-established approach. The experimentally obtained non-naturally distributed 13C labeling patterns of metabolite pools can be measured by mass spectrometric detection with front-end separation and can be consequently incorporated into biochemical network models. Here, despite a tedious derivatization step, gas chromatographic separation of polar metabolites is favorable because of the wide coverage range and high isomer separation efficiency. However, the typically employed electron ionization energy of 70 eV leads to significant fragmentation and consequently only low-abundant ions with an intact carbon backbone. Since these ions are considered a prerequisite for the analysis of the non-naturally distributed labeling patterns and further integration into modeling strategies, a softer ionization technique is needed. In the present work, a novel low energy electron ionization source is optimized for the analysis of primary metabolites and compared with a chemical ionization approach in terms of trueness, precision, and sensitivity.
Collapse
Affiliation(s)
- Teresa Mairinger
- Department of Chemistry, University of Natural Resources and Life Sciences-BOKU Vienna, Muthgasse 18, 1190, Vienna, Austria.
- EAWAG: Swiss Federal Institute of Aquatic Science and Technology, Ueberlandstrasse 133, 8600, Dübendorf, Switzerland.
| | - Jennifer Sanderson
- Agilent Technologies Inc, 5301 Stevens Creek Boulevard, Santa Clara, CA, 95051, USA
| | - Stephan Hann
- Department of Chemistry, University of Natural Resources and Life Sciences-BOKU Vienna, Muthgasse 18, 1190, Vienna, Austria
| |
Collapse
|
13
|
Metabolic engineering of Pichia pastoris. Metab Eng 2018; 50:2-15. [PMID: 29704654 DOI: 10.1016/j.ymben.2018.04.017] [Citation(s) in RCA: 132] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2018] [Revised: 04/16/2018] [Accepted: 04/23/2018] [Indexed: 12/11/2022]
Abstract
Besides its use for efficient production of recombinant proteins the methylotrophic yeast Pichia pastoris (syn. Komagataella spp.) has been increasingly employed as a platform to produce metabolites of varying origin. We summarize here the impressive methodological developments of the last years to model and analyze the metabolism of P. pastoris, and to engineer its genome and metabolic pathways. Efficient methods to insert, modify or delete genes via homologous recombination and CRISPR/Cas9, supported by modular cloning techniques, have been reported. An outstanding early example of metabolic engineering in P. pastoris was the humanization of protein glycosylation. More recently the cell metabolism was engineered also to enhance the productivity of heterologous proteins. The last few years have seen an increased number of metabolic pathway design and engineering in P. pastoris, mainly towards the production of complex (secondary) metabolites. In this review, we discuss the potential role of P. pastoris as a platform for metabolic engineering, its strengths, and major requirements for future developments of chassis strains based on synthetic biology principles.
Collapse
|