Gas chromatography-mass spectrometry analysis of 13C labeling in sugars for metabolic flux analysis

Metabolic flux analysis, using 13C labeled substrates, has become a powerful methodology for quantifying intracellular fluxes. Most often, analysis is restricted to nuclear magnetic resonance or mass spectrometry measurement of 13C label incorporation into protein amino acids. However, amino acid isotopomer distribution insufficiently covers the entire network of central metabolism, especially in plant cells with highly compartmented metabolism, and analysis of other metabolites is required. Analysis of label in saccharides provides complementary data to better define fluxes around hexose, pentose, and triose phosphate pools. Here, we propose a gas chromatography-mass spectrometry (GC-MS) method to analyze 13C labeling in glucose and fructose moieties of sucrose, free glucose, fructose, maltose, inositol, and starch. Our results show that saccharide labeling for isotopomer quantification is better analyzed by chemical ionization than by electron ionization. The structure of the generated fragments was simulated and validated using labeled standards. The method is illustrated by analysis of saccharides extracted from developing rapeseed (Brassica napus L.) embryos. It is shown that glucose 6-phosphate isomerase and plastidial glucose 6-phosphate transport reactions are not at equilibrium, and light is shed on the pathways leading to fructose, maltose, and inositol synthesis.     

13C metabolic flux analysis: optimal design of isotopic labeling experiments

1. Download : Download high-res image (292KB)
2. Download : Download full-size image

Highlights► Design of isotopic labeling experiments determines the precision of flux estimates. ► New approaches were developed to select optimal isotopic tracers and measurements. ► Tandem mass spectrometry provides more informative data for ¹³C flux analysis. ► Parallel labeling experiments can improve the quality of flux estimates. ► New methods are needed for optimal design of parallel labeling experiments.     

Parallel labeling experiments for pathway elucidation and 13C metabolic flux analysis

1. Download : Download high-res image (170KB)
2. Download : Download full-size image

     

13C metabolic flux analysis.

Metabolic flux analysis using 13C-labeled substrates has become an important tool in metabolic engineering. It allows the detailed quantification of all intracellular fluxes in the central metabolism of a microorganism. The method has strongly evolved in recent years by the introduction of new experimental procedures, measurement techniques, and mathematical data evaluation methods. Many of these improvements require advanced skills in the application of nuclear magnetic resonance and mass spectrometry techniques on the one hand and computational and statistical experience on the other hand. This minireview summarizes these recent developments and sketches the major practical problems. An outlook to possible future developments concludes the text.     

A transient isotopic labeling methodology for 13C metabolic flux analysis of photoautotrophic microorganisms

Metabolic flux analysis is increasingly recognized as an integral component of systems biology. However, techniques for experimental measurement of system-wide metabolic fluxes in purely photoautotrophic systems (growing on CO(2) as the sole carbon source) have not yet been developed due to the unique problems posed by such systems. In this paper, we demonstrate that an approach that balances positional isotopic distributions transiently is the only route to obtaining system-wide metabolic flux maps for purely autotrophic metabolism. The outlined transient (13)C-MFA methodology enables measurement of fluxes at a metabolic steady-state, while following changes in (13)C-labeling patterns of metabolic intermediates as a function of time, in response to a step-change in (13)C-label input. We use mathematical modeling of the transient isotopic labeling patterns of central intermediates to assess various experimental requirements for photoautotrophic MFA. This includes the need for intracellular metabolite concentration measurements and isotopic labeling measurements as a function of time. We also discuss photobioreactor design and operation in order to measure fluxes under precise environmental conditions. The transient MFA technique can be used to measure and compare fluxes under different conditions of light intensity, nitrogen sources or compare strains with various mutations or gene deletions and additions.     

A universal framework for 13C metabolic flux analysis.

A general methodology is presented for the modeling, simulation, design, evaluation, and statistical analysis of (13)C-labeling experiments for metabolic flux analysis. The universal software framework 13C-FLUX was implemented to support all steps of this process. Guided by the example of anaplerotic flux determination in Corynebacterium glutamicum, the technical details of the model setup, experimental design, and data evaluation are discussed. It is shown how the network structure, the input substrate composition, the assumptions about fluxes, and the measurement configuration are specified within 13C-FLUX. Based on the network model, different experimental designs are computed depending on the goal of the investigations. Finally, a specific experiment is evaluated and the various statistical methods used to analyze the results are briefly explained. The appendix gives some details about the software implementation and availability.     

Parallel labeling experiments with [U-(13)C]glucose validate E. coli metabolic network model for (13)C metabolic flux analysis

(13)C-metabolic flux analysis (MFA) is a widely used method for measuring intracellular metabolic fluxes in living cells. (13)C MFA relies on several key assumptions: (1) the assumed metabolic network model is complete, in that it accounts for all significant enzymatic and transport reactions; (2) (13)C-labeling measurements are accurate and precise; and (3) enzymes and transporters do not discriminate between (12)C- and (13)C-labeled metabolites. In this study, we tested these inherent assumptions of (13)C MFA for wild-type E. coli by parallel labeling experiments with [U-(13)C]glucose as tracer. Cells were grown in six parallel cultures in custom-constructed mini-bioreactors, starting from the same inoculum, on medium containing different mixtures of natural glucose and fully labeled [U-(13)C]glucose, ranging from 0% to 100% [U-(13)C]glucose. Macroscopic growth characteristics of E. coli showed no observable kinetic isotope effect. The cells grew equally well on natural glucose, 100% [U-(13)C]glucose, and mixtures thereof. (13)C MFA was then used to determine intracellular metabolic fluxes for several metabolic network models: an initial network model from literature; and extended network models that accounted for potential dilution effects of isotopic labeling. The initial network model did not give statistically acceptable fits and produced inconsistent flux results for the parallel labeling experiments. In contrast, an extended network model that accounted for dilution of intracellular CO(2) by exchange with extracellular CO(2) produced statistically acceptable fits, and the estimated metabolic fluxes were consistent for the parallel cultures. This study illustrates the importance of model validation for (13)C MFA. We show that an incomplete network model can produce statistically unacceptable fits, as determined by a chi-square test for goodness-of-fit, and return biased metabolic fluxes. The validated metabolic network model for E. coli from this study can be used in future investigations for unbiased metabolic flux measurements.     

Bioreaction network topology and metabolic flux ratio analysis by biosynthetic fractional 13C labeling and two-dimensional NMR spectroscopy

Biosynthetically directed fractional 13C labeling of the proteinogenic amino acids is achieved by feeding a mixture of uniformly 13C-labeled and unlabeled carbon source compounds into a bioreaction network. Analysis of the resulting labeling pattern enables both a comprehensive characterization of the network topology and the determination of metabolic flux ratios. Attractive features with regard to routine applications are (i) an inherently small demand for 13C-labeled source compounds and (ii) the high sensitivity of two-dimensional [13C,1H]-correlation nuclear magnetic resonance spectroscopy for analysis of 13C-labeling patterns. A user-friendly program, FCAL, is available to allow rapid data analysis. This novel approach, which recently also has been employed in conjunction with metabolic flux balancing to obtain reliable estimates of in vivo fluxes, enables efficient support of metabolic engineering and biotechnology process design.     

Towards high throughput metabolic flux analysis in plants

Research on plant metabolism is currently experiencing the common use of various omics methods creating valuable information on the concentrations of the cell's constituents. However, little is known about in vivo reaction rates, which can be determined by Metabolic Flux Analysis (MFA), a combination of isotope labeling experiments and computer modeling of the metabolic network. Large-scale applications of this method so far have been hampered by tedious procedures of tissue culture, analytics, modeling and simulation. By streamlining the workflow of MFA, the throughput of the method could be significantly increased. We propose strategies for these improvements on various sub-steps which will move flux analysis to the medium-throughput range and closer to established methods such as metabolite profiling. Furthermore, this may enable novel applications of MFA, for example screening plant populations for traits related to the flux phenotype.     

Optimization of 13C isotopic tracers for metabolic flux analysis in mammalian cells

Mammalian cells consume and metabolize various substrates from their surroundings for energy generation and biomass synthesis. Glucose and glutamine, in particular, are the primary carbon sources for proliferating cancer cells. While this combination of substrates generates static labeling patterns for use in (13)C metabolic flux analysis (MFA), the inability of single tracers to effectively label all pathways poses an obstacle for comprehensive flux determination within a given experiment. To address this issue we applied a genetic algorithm to optimize mixtures of (13)C-labeled glucose and glutamine for use in MFA. We identified tracer combinations that minimized confidence intervals in an experimentally determined flux network describing central carbon metabolism in tumor cells. Additional simulations were used to determine the robustness of the [1,2-(13)C(2)]glucose/[U-(13)C(5)]glutamine tracer combination with respect to perturbations in the network. Finally, we experimentally validated the improved performance of this tracer set relative to glucose tracers alone in a cancer cell line. This versatile method allows researchers to determine the optimal tracer combination to use for a specific metabolic network, and our findings applied to cancer cells significantly enhance the ability of MFA experiments to precisely quantify fluxes in higher organisms.     

High-throughput metabolic flux analysis based on gas chromatography-mass spectrometry derived 13C constraints

13C-constrained flux balancing analysis based on gas chromatography-mass spectrometry data is presented here as a simple and robust method for the estimation of intracellular carbon fluxes. In this approach, the underdetermined system of metabolite balances deduced from stoichiometric relations and measured extracellular rates is complemented with 13C constraints from metabolic flux ratio analysis. Fluxes in central carbon metabolism of exponentially growing Escherichia coli were estimated by 13C-constrained flux balancing from three different 13C-labeled glucose experiments. The best resolution of the network was achieved using 13C constraints derived from [U-13C]glucose and [1-13C]glucose experiments. The corresponding flux estimate was in excellent agreement with a solution that was independently obtained with a comprehensive isotopomer model. This new methodology was also demonstrated to faithfully capture the intracellular flux distribution in E. coli shake flasks and 1-ml deep-well microtiter plates. Due to its simplicity, speed, and robustness, 13C-constrained metabolic flux balancing is promising for routine and high-throughput analysis on a miniaturized scale.     

Metabolic flux analysis at ultra short time scale: isotopically non-stationary 13C labeling experiments

A novel approach to (13)C metabolic flux analysis (MFA) is presented using cytosolic metabolite pool sizes and their (13)C labeling data from an isotopically non-stationary (13)C labeling experiment (INST-CLE). The procedure is demonstrated with an E. coli wild type strain grown at fed batch conditions. The intra cellular labeling dynamics are excited by a sudden step increase of the (13)C portion in the substrate feed. Due to unchanged saturation of the substrate uptake system, the metabolic fluxes remain constant during the following sampling time period of only 16s, in which 20 samples are taken by an automated rapid sampling device immediately stopping metabolism by methanol quenching. Subsequent cell disruptive sample preparation and LC-MS/MS enabled simultaneous determination of pool sizes and mass isotopomers of intra cellular metabolites requiring detection limits in the nM range. Based on this data the new computational flux analysis tool 13CFLUX/INST is used to determine the intra cellular fluxes based on a complex carbon labeling network model. The measured data is in good agreement with the model predictions, thus proving the applicability of the new isotopically non-stationary (13)C metabolic flux analysis (INST-(13)C-MFA) concept. Moreover, it is shown that significant new information with respect to flux identifiability, non-measurable pool sizes, data consistency, or large storage pools can be taken from the novel kind of experimental data. This offers new insight into the biological operation of the metabolic network in vivo.     

Comparative 13C metabolic flux analysis of pyruvate dehydrogenase complex-deficient, L-valine-producing Corynebacterium glutamicum

L-Valine can be formed successfully using C. glutamicum strains missing an active pyruvate dehydrogenase enzyme complex (PDHC). Wild-type C. glutamicum and four PDHC-deficient strains were compared by (13)C metabolic flux analysis, especially focusing on the split ratio between glycolysis and the pentose phosphate pathway (PPP). Compared to the wild type, showing a carbon flux of 69% ± 14% through the PPP, a strong increase in the PPP flux was observed in PDHC-deficient strains with a maximum of 113% ± 22%. The shift in the split ratio can be explained by an increased demand of NADPH for l-valine formation. In accordance, the introduction of the Escherichia coli transhydrogenase PntAB, catalyzing the reversible conversion of NADH to NADPH, into an L-valine-producing C. glutamicum strain caused the PPP flux to decrease to 57% ± 6%, which is below the wild-type split ratio. Hence, transhydrogenase activity offers an alternative perspective for sufficient NADPH supply, which is relevant for most amino acid production systems. Moreover, as demonstrated for L-valine, this bypass leads to a significant increase of product yield due to a concurrent reduction in carbon dioxide formation via the PPP.     

Quantifying plant phenotypes with isotopic labeling & metabolic flux analysis

1. Download : Download high-res image (297KB)
2. Download : Download full-size image

     

Elucidating the role of copper in CHO cell energy metabolism using 13C metabolic flux analysis

¹³C‐metabolic flux analysis was used to understand copper deficiency‐related restructuring of energy metabolism, which leads to excessive lactate production in recombinant protein‐producing CHO cells. Stationary‐phase labeling experiments with U‐¹³C glucose were conducted on CHO cells grown under high and limiting copper in 3 L fed‐batch bioreactors. The resultant labeling patterns of soluble metabolites were measured by GC‐MS and used to estimate metabolic fluxes in the central carbon metabolism pathways using OpenFlux. Fluxes were evaluated 300 times from stoichiometrically feasible random guess values and their confidence intervals calculated by Monte Carlo simulations. Results from metabolic flux analysis exhibited significant carbon redistribution throughout the metabolic network in cells under Cu deficiency. Specifically, glycolytic fluxes increased (25%–79% relative to glucose uptake) whereas fluxes through the TCA and pentose phosphate pathway (PPP) were lower (15%–23% and 74%, respectively) compared with the Cu‐containing condition. Furthermore, under Cu deficiency, 33% of the flux entering TCA via the pyruvate node was redirected to lactate and malate production. Based on these results, we hypothesize that Cu deficiency disrupts the electron transport chain causing ATP deficiency, redox imbalance, and oxidative stress, which in turn drive copper‐deficient CHO cells to produce energy via aerobic glycolysis, which is associated with excessive lactate production, rather than the more efficient route of oxidative phosphorylation. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:1179–1186, 2015     

Metabolic flux analysis as a tool in metabolic engineering of plants

Methods of metabolic flux analysis (MFA) provide insights into the theoretical capabilities of metabolic networks and allow probing the in vivo performance of cellular metabolism. In recent years, an increasing awareness has developed that network analysis methods within the systems biology toolbox are serving to improve our understanding and ability to manipulate metabolism. In this minireview the potential of MFA to increase the chances of success in metabolic engineering of plants is presented, recent progress related to engineering and flux analysis in central metabolism of plants is discussed, and some recent advances in flux analysis methodology are highlighted.     

13C metabolic flux analysis for larger scale cultivation using gas chromatography-combustion-isotope ratio mass spectrometry

¹³C-based metabolic flux analysis (¹³CMFA) is limited to smaller scale experiments due to very high costs of labeled substrates. We measured ¹³C enrichment in proteinogenic amino acid hydrolyzates using gas chromatography-combustion-isotope ratio mass spectrometry (GC-C-IRMS) from a series of parallel batch cultivations of Corynebacterium glutamicum utilizing mixtures of natural glucose and [1-¹³C] glucose, containing 0%, 0.5%, 1%, 2%, and 10% [1-¹³C] glucose. Decreasing the [1-¹³C] glucose content, kinetic isotope effects played an increasing role but could be corrected. From the corrected ¹³C enrichments in vivo fluxes in the central metabolism were determined by numerical optimization. The obtained flux distribution was very similar to those obtained from parallel labeling experiments using conventional high labeling GC-MS method and to published results. The GC-C-IRMS-based method involving low labeling degree of expensive tracer substrate, e.g. 1%, is well suited for larger laboratory and industrial pilot scale fermentations.     

13C tracer experiments and metabolite balancing for metabolic flux analysis: comparing two approaches

Conventional metabolic flux analysis uses the information gained from determination of measurable fluxes and a steady-state assumption for intracellular metabolites to calculate the metabolic fluxes in a given metabolic network. The determination of intracellular fluxes depends heavily on the correctness of the assumed stoichiometry including the presence of all reactions with a noticeable impact on the model metabolite balances. Determination of fluxes in complex metabolic networks often requires the inclusion of NADH and NADPH balances, which are subject to controversial debate. Transhydrogenation reactions that transfer reduction equivalents from NADH to NADPH or vice versa can usually not be included in the stoichiometric model, because they result in singularities in the stoichiometric matrix. However, it is the NADPH balance that, to a large extent, determines the calculated flux through the pentose phosphate pathway. Hence, wrong assumptions on the presence or activity of transhydrogenation reactions will result in wrong estimations of the intracellular flux distribution. Using 13C tracer experiments and NMR analysis, flux analysis can be performed on the basis of only well established stoichiometric equations and measurements of the labeling state of intracellular metabolites. Neither NADH/NADPH balancing nor assumptions on energy yields need to be included to determine the intracellular fluxes. Because metabolite balancing methods and the use of 13C labeling measurements are two different approaches to the determination of intracellular fluxes, both methods can be used to verify each other or to discuss the origin and significance of deviations in the results. Flux analysis based entirely on metabolite balancing and flux analysis, including labeling information, have been performed independently for a wild-type strain of Aspergillus oryzae producing alpha-amylase. Two different nitrogen sources, NH4+ and NO3-, have been used to investigate the influence of the NADPH requirements on the intracellular flux distribution. The two different approaches to the calculation of fluxes are compared and deviations in the results are discussed. Copyright 1998 John Wiley & Sons, Inc.     

MetaFluxNet: the management of metabolic reaction information and quantitative metabolic flux analysis

SUMMARY: MetaFluxNet is a program package for managing information on the metabolic reaction network and for quantitatively analyzing metabolic fluxes in an interactive and customized way. It allows users to interpret and examine metabolic behavior in response to genetic and/or environmental modifications. As a result, quantitative in silico simulations of metabolic pathways can be carried out to understand the metabolic status and to design the metabolic engineering strategies. The main features of the program include a well-developed model construction environment, user-friendly interface for metabolic flux analysis (MFA), comparative MFA of strains having different genotypes under various environmental conditions, and automated pathway layout creation. AVAILABILITY: http://mbel.kaist.ac.kr/ Supplementary information: A manual for MetaFluxNet is available as PDF file.     

Combining metabolic flux analysis tools and 13C NMR to estimate intracellular fluxes of cultured astrocytes

In this work, brain cell metabolism was investigated by (13)C NMR spectroscopy and metabolic flux analysis (MFA). Monotypic cultures of astrocytes were incubated with labeled glucose for 38 h, and the distribution of the label was analyzed by (13)C NMR spectroscopy. The analysis of the spectra reveals two distinct physiological states characterized by different ratios of pyruvate carboxylase to pyruvate dehydrogenase activities (PC/PDH). Intracellular flux distributions for both metabolic states were estimated by MFA using the isotopic information and extracellular rate measurements as constraints. The model was subsequently checked with the consistency index method. From a biological point of view, the occurrence of the two physiological states appears to be correlated with the presence or absence of extracellular glutamate. Concerning the model, it can be stated that the metabolic network and the set of constraints adopted provide a consistent and robust characterization of the astrocytic metabolism, allowing for the calculation of central intracellular fluxes such as pyruvate recycling, the anaplerotic flux mediated by pyruvate carboxylase, and the glutamine formation through glutamine synthetase.