13C-Labeled metabolic flux analysis of a fed-batch culture of elutriated Saccharomyces cerevisiae

This study addresses the question of whether observable changes in fluxes in the primary carbon metabolism of Saccharomyces cerevisiae occur between the different phases of the cell division cycle. To detect such changes by metabolic flux analysis, a 13C-labeling experiment was performed with a fed-batch culture inoculated with a partially synchronized cell population obtained through centrifugal elutriation. Such a culture exhibits dynamic changes in the fractions of cells in different cell cycle phases over time. The mass isotopomer distributions of free intracellular metabolites in central carbon metabolism were measured by liquid chromatography-mass spectrometry. For four time points during the culture, these distributions were used to obtain the best estimates for the metabolic fluxes. The obtained flux fits suggested that the optimally fitted split ratio for the pentose phosphate pathway changed by almost a factor of 2 up and down around a value of 0.27 during the experiment. Statistical analysis revealed that some of the fitted flux distributions for different time points were significantly different from each other, indicating that cell cycle-dependent variations in cytosolic metabolic fluxes indeed occurred.     

Shewanella oneidensis MR-1 fluxome under various oxygen conditions

The central metabolic fluxes of Shewanella oneidensis MR-1 were examined under carbon-limited (aerobic) and oxygen-limited (microaerobic) chemostat conditions, using 13C-labeled lactate as the sole carbon source. The carbon labeling patterns of key amino acids in biomass were probed using both gas chromatography-mass spectrometry (GC-MS) and 13C nuclear magnetic resonance (NMR). Based on the genome annotation, a metabolic pathway model was constructed to quantify the central metabolic flux distributions. The model showed that the tricarboxylic acid (TCA) cycle is the major carbon metabolism route under both conditions. The Entner-Doudoroff and pentose phosphate pathways were utilized primarily for biomass synthesis (with a flux below 5% of the lactate uptake rate). The anaplerotic reactions (pyruvate to malate and oxaloacetate to phosphoenolpyruvate) and the glyoxylate shunt were active. Under carbon-limited conditions, a substantial amount (9% of the lactate uptake rate) of carbon entered the highly reversible serine metabolic pathway. Under microaerobic conditions, fluxes through the TCA cycle decreased and acetate production increased compared to what was found for carbon-limited conditions, and the flux from glyoxylate to glycine (serine-glyoxylate aminotransferase) became measurable. Although the flux distributions under aerobic, microaerobic, and shake flask culture conditions were different, the relative flux ratios for some central metabolic reactions did not differ significantly (in particular, between the shake flask and aerobic-chemostat groups). Hence, the central metabolism of S. oneidensis appears to be robust to environmental changes. Our study also demonstrates the merit of coupling GC-MS with 13C NMR for metabolic flux analysis to reduce the use of 13C-labeled substrates and to obtain more-accurate flux values.     

Responses of the central metabolism in Escherichia coli to phosphoglucose isomerase and glucose-6-phosphate dehydrogenase knockouts

The responses of Escherichia coli central carbon metabolism to knockout mutations in phosphoglucose isomerase and glucose-6-phosphate (G6P) dehydrogenase genes were investigated by using glucose- and ammonia-limited chemostats. The metabolic network structures and intracellular carbon fluxes in the wild type and in the knockout mutants were characterized by using the complementary methods of flux ratio analysis and metabolic flux analysis based on [U-(13)C]glucose labeling and two-dimensional nuclear magnetic resonance (NMR) spectroscopy of cellular amino acids, glycerol, and glucose. Disruption of phosphoglucose isomerase resulted in use of the pentose phosphate pathway as the primary route of glucose catabolism, while flux rerouting via the Embden-Meyerhof-Parnas pathway and the nonoxidative branch of the pentose phosphate pathway compensated for the G6P dehydrogenase deficiency. Furthermore, additional, unexpected flux responses to the knockout mutations were observed. Most prominently, the glyoxylate shunt was found to be active in phosphoglucose isomerase-deficient E. coli. The Entner-Doudoroff pathway also contributed to a minor fraction of the glucose catabolism in this mutant strain. Moreover, although knockout of G6P dehydrogenase had no significant influence on the central metabolism under glucose-limited conditions, this mutation resulted in extensive overflow metabolism and extremely low tricarboxylic acid cycle fluxes under ammonia limitation conditions.     

Isotopically nonstationary 13C metabolic flux analysis in resting and activated human platelets

Platelet metabolism is linked to platelet hyper- and hypoactivity in numerous human diseases. Developing a detailed understanding of the link between metabolic shifts and platelet activation state is integral to improving human health. Here, we show the first application of isotopically nonstationary ¹³C metabolic flux analysis to quantitatively measure carbon fluxes in both resting and thrombin activated platelets. Metabolic flux analysis results show that resting platelets primarily metabolize glucose to lactate via glycolysis, while acetate is oxidized to fuel the tricarboxylic acid cycle. Upon activation with thrombin, a potent platelet agonist, platelets increase their uptake of glucose 3-fold. This results in an absolute increase in flux throughout central metabolism, but when compared to resting platelets they redistribute carbon dramatically. Activated platelets decrease relative flux to the oxidative pentose phosphate pathway and TCA cycle from glucose and increase relative flux to lactate. These results provide the first report of reaction-level carbon fluxes in platelets and allow us to distinguish metabolic fluxes with much higher resolution than previous studies.     

Shewanella oneidensis MR-1 Fluxome under Various Oxygen Conditions

The central metabolic fluxes of Shewanella oneidensis MR-1 were examined under carbon-limited (aerobic) and oxygen-limited (microaerobic) chemostat conditions, using ¹³C-labeled lactate as the sole carbon source. The carbon labeling patterns of key amino acids in biomass were probed using both gas chromatography-mass spectrometry (GC-MS) and ¹³C nuclear magnetic resonance (NMR). Based on the genome annotation, a metabolic pathway model was constructed to quantify the central metabolic flux distributions. The model showed that the tricarboxylic acid (TCA) cycle is the major carbon metabolism route under both conditions. The Entner-Doudoroff and pentose phosphate pathways were utilized primarily for biomass synthesis (with a flux below 5% of the lactate uptake rate). The anaplerotic reactions (pyruvate to malate and oxaloacetate to phosphoenolpyruvate) and the glyoxylate shunt were active. Under carbon-limited conditions, a substantial amount (9% of the lactate uptake rate) of carbon entered the highly reversible serine metabolic pathway. Under microaerobic conditions, fluxes through the TCA cycle decreased and acetate production increased compared to what was found for carbon-limited conditions, and the flux from glyoxylate to glycine (serine-glyoxylate aminotransferase) became measurable. Although the flux distributions under aerobic, microaerobic, and shake flask culture conditions were different, the relative flux ratios for some central metabolic reactions did not differ significantly (in particular, between the shake flask and aerobic-chemostat groups). Hence, the central metabolism of S. oneidensis appears to be robust to environmental changes. Our study also demonstrates the merit of coupling GC-MS with ¹³C NMR for metabolic flux analysis to reduce the use of ¹³C-labeled substrates and to obtain more-accurate flux values.     

Collisional fragmentation of central carbon metabolites in LC-MS/MS increases precision of ¹³C metabolic flux analysis

Experimental determination of fluxes by (13)C-tracers relies on detection of (13)C-patterns in metabolites or by-products. In the field of (13)C metabolic flux analysis, the most recent developments point toward recording labeling patterns by liquid chromatography (LC)-mass spectrometry (MS)/MS directly in intermediates in central carbon metabolism (CCM) to increase temporal resolution. Surprisingly, the flux studies published so far with LC-MS measurements were based on intact metabolic intermediates-thus neglected the potential benefits of using positional information to improve flux estimation. For the first time, we exploit collisional fragmentation to obtain more fine-grained (13)C-data on intermediates of CCM and investigate their impact in (13)C metabolic flux analysis. For the case study of Bacillus subtilis grown in mineral medium with (13)C-labeled glucose, we compare the flux estimates obtained by iterative isotopologue balancing of (13)C-data obtained either by LC-MS/MS for solely intact intermediates or LC-MS/MS for intact and fragmented intermediates of CCM. We show that with LC-MS/MS data, fragment information leads to more precise estimates of fluxes in pentose phosphate pathway, glycolysis, and to the tricarboxylic acid cycle. Additionally, we present an efficient analytical strategy to rapidly acquire large sets of (13)C-patterns by tandem MS, and an in-depth analysis of the collisional fragmentation of primary intermediates. In the future, this catalogue will enable comprehensive in silico calculability analyses to identify the most sensitive measurements and direct experimental design.     

Respirometric 13C flux analysis--Part II: in vivo flux estimation of lysine-producing Corynebacterium glutamicum

A novel method for metabolic flux studies of central metabolism which is based on respirometric (13)C flux analysis, i.e., parallel (13)C tracer studies with online CO(2) labeling measurements is applied to flux quantification of a lysine-producing mutant of Corynebacterium glutamicum. For this purpose, 3 respirometric (13)C labeling experiments with [1-(13)C(1)], [6-(13)C(1)] and [1,6-(13)C(2)] glucose were carried out in parallel. All fluxes comprising the reactions of glycolysis, of TCA cycle, of C3- and C4-metabolite interconversion and of lysine biosynthesis as well as the net reactions in the pentose phosphate pathway could be quantified solely using experimental data obtained from CO(2) labeling and extracellular rate measurements. At key branch points, 68+/-5% of glucose 6-phosphate were observed to be metabolized into pentose phosphate pathway and 48+/-1% of pyruvate into TCA cycle via pyruvate dehydrogenase. The results showed a good agreement with the previous studies using (13)C tracer cultivation and GC/MS analysis of proteinogenic amino acids. Also, respiratory quotient calculated from flux estimates using redox balance showed a high accordance with the value determined directly from the measured specific rates of O(2) consumption and CO(2) production. The results strongly support that the respirometric (13)C metabolic flux analysis is suited as an alternative to the conventional methods to study functional and regulatory activities of cells. The developed method is applicable to study growing or non-growing cells, primary and secondary metabolism and immobilized cells. Due to the non-accumulating nature of CO(2) labeling and instantaneous nature of the resulting fluxes, the method can also be used for dynamic profiling of metabolic activities. Therefore, it is complementary to conventional methods for metabolic flux analysis.     

The metabolic architecture of plant cells. Stability of central metabolism and flexibility of anabolic pathways during the growth cycle of tomato cells

The changes in the intermediary metabolism of plant cells were quantified according to growth conditions at three different stages of the growth cycle of tomato cell suspension. Eighteen fluxes of central metabolism were calculated from (13)C enrichments after near steady-state labeling by a metabolic model similar to that described in Dieuaide-Noubhani et al. (Dieuaide-Noubhani, M., Raffard, G., Canioni, P., Pradet, A., and Raymond, P. (1995) J. Biol. Chem. 270, 13147-13159), and 10 net fluxes were obtained directly from end-product accumulation rates. The absolute flux values of central metabolic pathways gradually slowed down with the decrease of glucose influx into the cells. However, the relative fluxes of glycolysis, the pentose-P pathway, and the tricarboxylic acid cycle remained unchanged during the culture cycle at 70, 28, and 40% of glucose influx, respectively, and the futile cycle of sucrose remained high at about 6-fold the glucose influx, independently from carbon nutritional conditions. This natural resistance to flux alterations is referred to as metabolic stability. The numerous anabolic pathways, including starch synthesis, hexose accumulation, biosynthesis of wall polysaccharides, and amino and organic acid biosynthesis were comparatively low and variable. The phosphoenolpyruvate carboxylase flux decreased 5-fold in absolute terms and 2-fold in relation to the glucose influx rate during the culture cycle. We conclude that anabolic fluxes constitute the flexible part of plant cell metabolism that can fluctuate in relation to cell demands for growth.     

Understanding flux in plant metabolic networks

The revolutionary growth in our ability to identify the 'parts list' of cellular infrastructure in plants in detail, and to alter it with precision, challenges us to develop methods to quantify how these parts function. For components of metabolism, this means mapping fluxes at the level of metabolic networks. Advances in experimental, analytical and software tools for metabolic flux analysis now allow maps of the fluxes through central metabolism to be obtained from the results of stable-isotope-labeling experiments. Such maps have led to notable successes in understanding and engineering metabolic function in microorganisms. Recent studies in plants are giving insight into particular fluxes, such as those of the pentose phosphate pathway, and into general phenomena, such as substrate- or futile-cycles and compartmentation. The importance of experimental design and statistical analysis have been illustrated, and analyses of fluxes in heterotrophic plant tissues have been carried out recently.     

Identification of Enzymes and Quantification of Metabolic Fluxes in the Wild Type and in a Recombinant Aspergillus oryzae Strain

Two α-amylase-producing strains of Aspergillus oryzae, a wild-type strain and a recombinant containing additional copies of the α-amylase gene, were characterized with respect to enzyme activities, localization of enzymes to the mitochondria or cytosol, macromolecular composition, and metabolic fluxes through the central metabolism during glucose-limited chemostat cultivations. Citrate synthase and isocitrate dehydrogenase (NAD) activities were found only in the mitochondria, glucose-6-phosphate dehydrogenase and glutamate dehydrogenase (NADP) activities were found only in the cytosol, and isocitrate dehydrogenase (NADP), glutamate oxaloacetate transaminase, malate dehydrogenase, and glutamate dehydrogenase (NAD) activities were found in both the mitochondria and the cytosol. The measured biomass components and ash could account for 95% (wt/wt) of the biomass. The protein and RNA contents increased linearly with increasing specific growth rate, but the carbohydrate and chitin contents decreased. A metabolic model consisting of 69 fluxes and 59 intracellular metabolites was used to calculate the metabolic fluxes through the central metabolism at several specific growth rates, with ammonia or nitrate as the nitrogen source. The flux through the pentose phosphate pathway increased with increasing specific growth rate. The fluxes through the pentose phosphate pathway were 15 to 26% higher for the recombinant strain than for the wild-type strain.     

Metabolism of peptide amino acids by Chinese hamster ovary cells grown in a complex medium

Metabolic flux analysis is a useful tool for unraveling relationships between metabolism and cell function. Material balancing can be used to provide estimates of major metabolic pathway fluxes, provided all significant metabolite uptake and production rates are measured. Potential sources of metabolizable material in many serum-free media formulations are low molecular weight digests of biological material such as yeast extracts and plant or animal tissue hydrolysates. These digests typically contain large amounts of peptides, which may be utilized as amino acids. This article demonstrates the need for accounting for amino acids liberated from peptides in order to accurately estimate pathway fluxes in Chinese hamster ovary cells grown in a complex (hydrolysate containing) medium. A simplified model of central carbon metabolism provides the framework for analyzing external metabolite measurements. Redundant measurements are included to ensure the consistency of data and assumed biochemistry by comparing redundant measurements with their predicted values from a minimum data set, and by expressing the degree of agreement using a statistical "consistency index." The consistency index tests whether redundancies are satisfied within expected experimental error. For chemostat steady states of CHO cultures grown in a hydrolysate-supplemented medium, consistent data were obtained only when amino acids liberated from peptides were taken into account.     

Metabolic flux analysis using stoichiometric models for Aspergillus niger: comparison under glucoamylase-producing and non-producing conditions

Aspergillus niger AB1.13 cultures with glucoamylase production (with D-glucose as substrate) and without glucoamylase production (with D-xylose as substrate) were characterized by metabolic flux analysis. Two comprehensive metabolic models for d-glucose- as well as for D-xylose-consumption were used to quantify and compare the metabolic fluxes through the central pathways of carbon metabolism at different pH-values. The models consist of the most relevant metabolic pathways for A. niger including glycolysis, pentose-phosphate pathway, citrate cycle, energy metabolism and anaplerotic reactions comprising two intracellular compartments, the cytoplasm and mitochondrion. When D-xylose was used as the sole carbon source, the relative flux of the substrate through the oxidative pentose-phosphate pathway (PPP) via G6P-dehydrogenase was unaffected by the pH-value of the culture medium. About 30% of D-xylose consumed was routed through the oxidative PPP. In contrast, the flux of D-glucose (i.e., under glucoamylase-producing conditions) through the oxidative PPP was remarkably higher and, in addition was significantly affected by the pH-value of the culture medium (40% at pH 5.5, 56% at pH 3.7, respectively). Summarizing, the flux through the PPP under glucoamylase producing conditions was 30-90% higher than for non-producing conditions.     

Metabolic flux profiling of Escherichia coli mutants in central carbon metabolism using GC-MS.

We describe here a novel methodology for rapid diagnosis of metabolic changes, which is based on probabilistic equations that relate GC-MS-derived mass distributions in proteinogenic amino acids to in vivo enzyme activities. This metabolic flux ratio analysis by GC-MS provides a comprehensive perspective on central metabolism by quantifying 14 ratios of fluxes through converging pathways and reactions from [1-13C] and [U-13C]glucose experiments. Reliability and accuracy of this method were experimentally verified by successfully capturing expected flux responses of Escherichia coli to environmental modifications and seven knockout mutations in all major pathways of central metabolism. Furthermore, several mutants exhibited additional, unexpected flux responses that provide new insights into the behavior of the metabolic network in its entirety. Most prominently, the low in vivo activity of the Entner-Doudoroff pathway in wild-type E. coli increased up to a contribution of 30% to glucose catabolism in mutants of glycolysis and TCA cycle. Moreover, glucose 6-phosphate dehydrogenase mutants catabolized glucose not exclusively via glycolysis, suggesting a yet unidentified bypass of this reaction. Although strongly affected by environmental conditions, a stable balance between anaplerotic and TCA cycle flux was maintained by all mutants in the upper part of metabolism. Overall, our results provide quantitative insight into flux changes that bring about the resilience of metabolic networks to disruption.     

Central carbon metabolism influences cellulase production in Bacillus licheniformis

Bacillus licheniformis that can produce cellulase including endo glucanase and glucosidase is an important industrial microbe for cellulose degradation. The purpose of this research was to assess the effect of endo glucanase gene bglC and glucosidase gene bglH on the central metabolic flux in B. licheniformis. bglC and bglH were knocked out using homologous recombination method, respectively, and the corresponding knockout strains were obtained for ¹³C metabolic flux analysis. A significant change was observed in metabolic fluxes after ¹³C metabolic flux ratio analysis. In both of the knockout strains, the increased fluxes of the pentose phosphate pathway and malic enzyme reaction enabled an elevated supply of NADPH which provided enough reducing power for the in vivo synthesis reactions. The fluxes through tricarboxylic acid cycle and anaplerotic reactions increased fast in the two knockout strains, which meant more energy generated. The changed fluxes in central carbon metabolism provided a holistic view of the physiological status in B. licheniformis and possible targets for further strain engineering.

Significance and Impact of the Study

Cellulase is very important in the field of agriculture and bioenergy because of its degrading effect on cellulosic biomass. This study presented the effect of central carbon metabolism on cellulase production in Bacillus licheniformis. The study also provided a holistic view of the physiological status in B. licheniformis. The shifted metabolism provided a quantitative evaluation of the biosynthesis of cellulase and a priority ranked target list for further strain engineering.     

Response of the central metabolism in Corynebacterium glutamicum to the use of an NADH-dependent glutamate dehydrogenase

The extensive use of 13C enrichments in precursor metabolites for flux quantification does not rely on NADPH stoichiometries and can therefore be used to quantify reducing power fluxes. As an application of this concept, the NADPH fluxes were quantified in an L-lysine producer of Corynebacterium glutamicum grown into metabolic and isotopic steady state with [1-13C]glucose. In this case, where the organism's NADPH-dependent glutamate dehydrogenase consumes reducing power, the NADPH flux generated is 210% (molar flux relative to glucose uptake rate) with its major part (72% of the total) generated via the pentose phosphate pathway activity. An isogenic strain in which the glutamate dehydrogenase of C. glutamicum was replaced by the NADH-dependent glutamate dehydrogenase of Peptostreptococcus asaccharolyticus was made and the metabolite fluxes were again estimated. The major response to this local perturbation is a drastically reduced NADPH generation of only 139%. Most of the NADPH (62% of the total) is now generated via the tricarboxylic acid cycle activity. This shows the extraordinary flexibility of the central metabolism and provides a picture of the global regulatory properties of the central metabolism. Furthermore, a detailed analysis of the fluxes and exchange fluxes within the anaplerotic reactions is given. It is hypothesized that these reactions might also serve to balance the total reducing power budget as well as the energy budget within the cell.     

Using Bioconductor Package BiGGR for Metabolic Flux Estimation Based on Gene Expression Changes in Brain

Predicting the distribution of metabolic fluxes in biochemical networks is of major interest in systems biology. Several databases provide metabolic reconstructions for different organisms. Software to analyze flux distributions exists, among others for the proprietary MATLAB environment. Given the large user community for the R computing environment, a simple implementation of flux analysis in R appears desirable and will facilitate easy interaction with computational tools to handle gene expression data. We extended the R software package BiGGR, an implementation of metabolic flux analysis in R. BiGGR makes use of public metabolic reconstruction databases, and contains the BiGG database and the reconstruction of human metabolism Recon2 as Systems Biology Markup Language (SBML) objects. Models can be assembled by querying the databases for pathways, genes or reactions of interest. Fluxes can then be estimated by maximization or minimization of an objective function using linear inverse modeling algorithms. Furthermore, BiGGR provides functionality to quantify the uncertainty in flux estimates by sampling the constrained multidimensional flux space. As a result, ensembles of possible flux configurations are constructed that agree with measured data within precision limits. BiGGR also features automatic visualization of selected parts of metabolic networks using hypergraphs, with hyperedge widths proportional to estimated flux values. BiGGR supports import and export of models encoded in SBML and is therefore interoperable with different modeling and analysis tools. As an application example, we calculated the flux distribution in healthy human brain using a model of central carbon metabolism. We introduce a new algorithm termed Least-squares with equalities and inequalities Flux Balance Analysis (Lsei-FBA) to predict flux changes from gene expression changes, for instance during disease. Our estimates of brain metabolic flux pattern with Lsei-FBA for Alzheimer’s disease agree with independent measurements of cerebral metabolism in patients. This second version of BiGGR is available from Bioconductor.