Quantification of central metabolic fluxes in the facultative methylotroph methylobacterium extorquens AM1 using 13C-label tracing and mass spectrometry

The metabolic fluxes of central carbon metabolism were measured in chemostat-grown cultures of Methylobacterium extorquens AM1 with methanol as the sole organic carbon and energy source and growth-limiting substrate. Label tracing experiments were carried out using 70% (13)C-methanol in the feed, and the steady-state mass isotopomer distributions of amino acids derived from total cell protein were measured by gas chromatography coupled to mass spectrometry. Fluxes were calculated from the isotopomer distribution data using an isotopomer balance model and evolutionary error minimization algorithm. The combination of labeled methanol with unlabeled CO(2), which enters central metabolism in two different reactions, provided the discriminatory power necessary to allow quantification of the unknown fluxes within a reasonably small confidence interval. In wild-type M. extorquens AM1, no measurable flux was detected through pyruvate dehydrogenase or malic enzyme, and very little flux through alpha-ketoglutarate dehydrogenase (1.4% of total carbon). In contrast, the alpha-ketoglutarate dehydrogenase flux was 25.5% of total carbon in the regulatory mutant strain phaR, while the pyruvate dehydrogenase and malic enzyme fluxes remained insignificant. The success of this technique with growth on C(1) compounds suggests that it can be applied to help characterize the effects of other regulatory mutations, and serve as a diagnostic tool in the metabolic engineering of methylotrophic bacteria.     

Flux analysis of central metabolic pathways in Geobacter metallireducens during reduction of soluble Fe(III)-nitrilotriacetic acid

We analyzed the carbon fluxes in the central metabolism of Geobacter metallireducens strain GS-15 using 13C isotopomer modeling. Acetate labeled in the first or second position was the sole carbon source, and Fe-nitrilotriacetic acid was the sole terminal electron acceptor. The measured labeled acetate uptake rate was 21 mmol/g (dry weight)/h in the exponential growth phase. The resulting isotope labeling pattern of amino acids allowed an accurate determination of the in vivo global metabolic reaction rates (fluxes) through the central metabolic pathways using a computational isotopomer model. The tracer experiments showed that G. metallireducens contained complete biosynthesis pathways for essential metabolism, and this strain might also have an unusual isoleucine biosynthesis route (using acetyl coenzyme A and pyruvate as the precursors). The model indicated that over 90% of the acetate was completely oxidized to CO2 via a complete tricarboxylic acid cycle while reducing iron. Pyruvate carboxylase and phosphoenolpyruvate (PEP) carboxykinase were present under these conditions, but enzymes in the glyoxylate shunt and malic enzyme were absent. Gluconeogenesis and the pentose phosphate pathway were mainly employed for biosynthesis and accounted for less than 3% of total carbon consumption. The model also indicated surprisingly high reversibility in the reaction between oxoglutarate and succinate. This step operates close to the thermodynamic equilibrium, possibly because succinate is synthesized via a transferase reaction, and the conversion of oxoglutarate to succinate is a rate-limiting step for carbon metabolism. These findings enable a better understanding of the relationship between genome annotation and extant metabolic pathways in G. metallireducens.     

Flux Analysis of Central Metabolic Pathways in Geobacter metallireducens during Reduction of Soluble Fe(III)-Nitrilotriacetic Acid

We analyzed the carbon fluxes in the central metabolism of Geobacter metallireducens strain GS-15 using ¹³C isotopomer modeling. Acetate labeled in the first or second position was the sole carbon source, and Fe-nitrilotriacetic acid was the sole terminal electron acceptor. The measured labeled acetate uptake rate was 21 mmol/g (dry weight)/h in the exponential growth phase. The resulting isotope labeling pattern of amino acids allowed an accurate determination of the in vivo global metabolic reaction rates (fluxes) through the central metabolic pathways using a computational isotopomer model. The tracer experiments showed that G. metallireducens contained complete biosynthesis pathways for essential metabolism, and this strain might also have an unusual isoleucine biosynthesis route (using acetyl coenzyme A and pyruvate as the precursors). The model indicated that over 90% of the acetate was completely oxidized to CO₂ via a complete tricarboxylic acid cycle while reducing iron. Pyruvate carboxylase and phosphoenolpyruvate (PEP) carboxykinase were present under these conditions, but enzymes in the glyoxylate shunt and malic enzyme were absent. Gluconeogenesis and the pentose phosphate pathway were mainly employed for biosynthesis and accounted for less than 3% of total carbon consumption. The model also indicated surprisingly high reversibility in the reaction between oxoglutarate and succinate. This step operates close to the thermodynamic equilibrium, possibly because succinate is synthesized via a transferase reaction, and the conversion of oxoglutarate to succinate is a rate-limiting step for carbon metabolism. These findings enable a better understanding of the relationship between genome annotation and extant metabolic pathways in G. metallireducens.     

GC-MS analysis of amino acids rapidly provides rich information for isotopomer balancing

Gas chromatography-mass spectrometry (GC-MS) is a rapid method that provides rich information on isotopomer distributions for metabolic flux analysis. First, we established a fast and reliable experimental protocol for GC-MS analysis of amino acids from total biomass hydrolyzates, and common experimental pitfalls are discussed. Second, a suitable interface for the use of mass isotopomer data is presented. For this purpose, a general, matrix-based correction procedure that accounts for naturally occurring isotopes is introduced. Simulated and experimentally determined mass distributions of unlabeled amino acids showed a deviation of less than 3% for 89% of the analyzed fragments. Third, to investigate general properties of GC-MS-based isotopomer balancing, altered flux ratios through glycolysis and pentose phosphate pathway and/or exchange fluxes were simulated. Different fluxes were found to exert specific and significant influence on the mass isotopomer distributions, thus indicating that GC-MS data contain valuable information for metabolic flux analysis. Fourth, comparison of different methods revealed that GC-MS analysis provides the largest number of independent constraints on amino acid isotopomer abundance, followed by correlation spectroscopy and fractional enrichment analysis by nuclear magnetic resonance.     

Flux quantification in central carbon metabolism of Catharanthus roseus hairy roots by 13C labeling and comprehensive bondomer balancing

Methods for accurate and efficient quantification of metabolic fluxes are desirable in plant metabolic engineering and systems biology. Toward this objective, we introduce the application of "bondomers", a computationally efficient and intuitively appealing alternative to the commonly used isotopomer concept, to flux evaluation in plants, by using Catharanthus roseus hairy roots as a model system. We cultured the hairy roots on (5% w/w U-(13)C, 95% w/w naturally abundant) sucrose, and acquired two-dimensional [(13)C, (1)H] and [(1)H, (1)H] NMR spectra of hydrolyzed aqueous extract from the hairy roots. Analysis of these spectra yielded a data set of 116 bondomers of beta-glucans and proteinogenic amino acids from the hairy roots. Fluxes were evaluated from the bondomer data by using comprehensive bondomer balancing. We identified most fluxes in a three-compartmental model of central carbon metabolism with good precision. We observed parallel pentose phosphate pathways in the cytosol and the plastid with significantly different fluxes. The anaplerotic fluxes between phosphoenolpyruvate and oxaloacetate in the cytosol and between malate and pyruvate in the mitochondrion were relatively high (60.1+/-2.5 mol per 100 mol sucrose uptake, or 22.5+/-0.5 mol per 100 mol mitochondrial pyruvate dehydrogenase flux). The development of a comprehensive flux analysis tool for this plant hairy root system is expected to be valuable in assessing the metabolic impact of genetic or environmental changes, and this methodology can be extended to other plant systems.     

Modeling isotopomer distributions in biochemical networks using isotopomer mapping matrices

Within the last decades NMR spectroscopy has undergone tremendous development and has become a powerful analytical tool for the investigation of intracellular flux distributions in biochemical networks using (13)C-labeled substrates. Not only are the experiments much easier to conduct than experiments employing radioactive tracer elements, but NMR spectroscopy also provides additional information on the labeling pattern of the metabolites. Whereas the maximum amount of information obtainable with (14)C-labeled substrates is the fractional enrichment in the individual carbon atom positions, NMR spectroscopy can also provide information on the degree of labeling at neighboring carbon atom positions by analyzing multiplet patterns in NMR spectra or using 2-dimensional NMR spectra. It is possible to quantify the mole fractions of molecules that show a specific labeling pattern, i.e., information of the isotopomer distribution in metabolite pools can be obtained. The isotopomer distribution is the maximum amount of information that in theory can be obtained from (13)C-tracer studies. The wealth of information contained in NMR spectra frequently leads to overdetermined algebraic systems. Consequently, fluxes must be estimated by nonlinear least squares analysis, in which experimental labeling data is compared with simulated steady state isotopomer distributions. Hence, mathematical models are required to compute the steady state isotopomer distribution as a function of a given set of steady state fluxes. Because 2(n) possible labeling patterns exist in a molecule of n carbon atoms, and each pattern corresponds to a separate state in the isotopomer model, these models are inherently complex. Model complexity, so far, has restricted usage of isotopomer information to relatively small metabolic networks. A general methodology for the formulation of isotopomer models is described. The model complexity of isotopomer models is reduced to that of classical metabolic models by expressing the 2(n) isotopomer mass balances of a metabolite pool in a single matrix equation. Using this approach an isotopomer model has been implemented that describes label distribution in primary carbon metabolism, i.e., in a metabolic network including the Embden-Meyerhof-Parnas and pentose phosphate pathway, the tricarboxylic acid cycle, and selected anaplerotic reaction sequences. The model calculates the steady state label distribution in all metabolite pools as a function of the steady state fluxes and is applied to demonstrate the effect of selected anaplerotic fluxes on the labeling pattern of the pathway intermediates. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55:831-840, 1997.     

Quantitative analysis of metabolic fluxes in Escherichia coli, using two-dimensional NMR spectroscopy and complete isotopomer models.

Complete isotopomer models that simulate distribution of label in 13C tracer experiments are applied to the quantification of metabolic fluxes in the primary carbon metabolism of E. coli under aerobic and anaerobic conditions. The concept of isotopomer mapping matrices (IMMs) is used to simplify the formulation of isotopomer mass balances by expressing all isotopomer mass balances of a metabolite pool in a single matrix equation. A numerically stable method to calculate the steady-state isotopomer distribution in metabolic networks in introduced. Net values of intracellular fluxes and the degree of reversibility of enzymatic steps are estimated by minimization of the deviations between experimental and simulated measurements. The metabolic model applied includes the Embden-Meyerhof-Parnas and the pentose phosphate pathway, the tricarboxylic acid cycle, anaplerotic reaction sequences and pathways involved in amino acid synthesis. The study clearly demonstrates the value of complete isotopomer models for maximizing the information obtainable from 13C tracer experiments. The approach applied here offers a completely general and comprehensive analysis of carbon tracer experiments where any set of experimental data on the labeling state and extracellular fluxes can be used for the quantification of metabolic fluxes in complex metabolic networks.     

Probing in vivo metabolism by stable isotope labeling of storage lipids and proteins in developing Brassica napus embryos

Developing embryos of Brassica napus accumulate both triacylglycerols and proteins as major storage reserves. To evaluate metabolic fluxes during embryo development, we have established conditions for stable isotope labeling of cultured embryos under steady-state conditions. Sucrose supplied via the endosperm is considered to be the main carbon and energy source for seed metabolism. However, in addition to 220 to 270 mM carbohydrates (sucrose, glucose, and fructose), analysis of endosperm liquid revealed up to 70 mM amino acids as well as 6 to 15 mM malic acid. Therefore, a labeling approach with multiple carbon sources is a precondition to quantitatively reflect fluxes of central carbon metabolism in developing embryos. Mid-cotyledon stage B. napus embryos were dissected from plants and cultured for 15 d on a complex liquid medium containing (13)C-labeled carbohydrates. The (13)C enrichment of fatty acids and amino acids (after hydrolysis of the seed proteins) was determined by gas chromatography/mass spectrometry. Analysis of (13)C isotope isomers of labeled fatty acids and plastid-derived amino acids indicated that direct glycolysis provides at least 90% of precursors of plastid acetyl-coenzyme A (CoA). Unlabeled amino acids, when added to the growth medium, did not reduce incorporation of (13)C label into plastid-formed fatty acids, but substantially diluted (13)C label in seed protein. Approximately 30% of carbon in seed protein was derived from exogenous amino acids and as a consequence, the use of amino acids as a carbon source may have significant influence on the total carbon and energy balance in seed metabolism. (13)C label in the terminal acetate units of C(20) and C(22) fatty acids that derive from cytosolic acetyl-CoA was also significantly diluted by unlabeled amino acids. We conclude that cytosolic acetyl-CoA has a more complex biogenetic origin than plastidic acetyl-CoA. Malic acid in the growth medium did not dilute (13)C label incorporation into fatty acids or proteins and can be ruled out as a source of carbon for the major storage components of B. napus embryos.     

Profiling of central metabolism in human cancer cells by two-dimensional NMR, GC-MS analysis, and isotopomer modeling

Tracking metabolic profiles has the potential to reveal crucial enzymatic steps that could be targeted in the drug discovery process. It is of special importance for various types of cancer known to be associated with substantial rewiring of metabolic networks. Here we introduce an integrated approach for the analysis of metabolome that allows us to simultaneously assess pathway activities (fluxes) and concentrations of a large number of the key components involved in central metabolism of human cells. This is accomplished by in vivo labeling with [U-¹³C]glucose followed by two-dimensional nuclear magnetic resonance (NMR) spectroscopy and gas chromatography-mass spectrometry (GC-MS) analysis. A comprehensive isotopomer model was developed, which enabled us to compare fluxes through the key central metabolic pathways including glycolysis, pentose phosphate pathway, tricarboxylic acid cycle, anaplerotic reactions, and biosynthetic pathways of fatty acids and amino acids. The validity and strength of this approach is illustrated by its application to a number of perturbations to breast cancer cells, including exposure to hypoxia, drug treatment, and tumor progression. We observed significant differences in the activities of specific metabolic pathways resulting from these perturbations and providing new mechanistic insights. Based on these findings we conclude that the developed metabolomic approach constitutes a promising analytical tool for revealing specific metabolic phenotypes in a variety of cell types and pathological conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Chen Yang and Adam D. Richardson contributed equally to this work.     

Identification of metabolic fluxes in hepatic cells from transient 13C-labeling experiments: Part I. Experimental observations

An experimental set-up for acquiring metabolite and transient (13)C-labeling data in mammalian cells is presented. An efficient sampling procedure was established for hepatic cells cultured in six-well plates as a monolayer attached to collagen, which allowed simultaneous quenching of metabolism and extraction of the intracellular intermediates of interest. Extracellular concentrations of glucose, amino acids, lactate, pyruvate, and urea were determined by GC-MS procedures and were used for estimation of metabolic uptake and excretion rates. Sensitive LC-MS and GC-MS methods were used to quantify the intracellular intermediates of tricarboxylic acid cycle, glycolysis, and pentose phosphate pathway and for the determination of isotopomer fractions of the respective metabolites. Mass isotopomer fractions were determined in a transient (13)C-labeling experiment using (13)C-labeled glucose as substrate. The absolute amounts of intracellular metabolites were obtained from a non-labeled experiment carried out in exactly the same way as the (13)C-labeling experiment, except that the media contained naturally labeled glucose only. Estimation of intracellular metabolic fluxes from the presented data is addressed in part II of this contribution.     

Central carbon metabolism of Saccharomyces cerevisiae explored by biosynthetic fractional (13)C labeling of common amino acids.

Aerobic and anaerobic central metabolism of Saccharomyces cerevisiae cells was explored in batch cultures on a minimal medium containing glucose as the sole carbon source, using biosynthetic fractional (13)C labeling of proteinogenic amino acids. This allowed, firstly, unravelling of the network of active central pathways in cytosol and mitochondria, secondly, determination of flux ratios characterizing glycolysis, pentose phosphate cycle, tricarboxylic acid cycle and C1-metabolism, and thirdly, assessment of intercompartmental transport fluxes of pyruvate, acetyl-CoA, oxaloacetate and glycine. The data also revealed that alanine aminotransferase is located in the mitochondria, and that amino acids are synthesized according to documented pathways. In both the aerobic and the anaerobic regime: (a) the mitochondrial glycine cleavage pathway is active, and efflux of glycine into the cytosol is observed; (b) the pentose phosphate pathways serve for biosynthesis only, i.e. phosphoenolpyruvate is entirely generated via glycolysis; (c) the majority of the cytosolic oxaloacetate is synthesized via anaplerotic carboxylation of pyruvate; (d) the malic enzyme plays a key role for mitochondrial pyruvate metabolism; (e) the transfer of oxaloacetate from the cytosol to the mitochondria is largely unidirectional, and the activity of the malate-aspartate shuttle and the succinate-fumarate carrier is low; (e) a large fraction of the mitochondrial pyruvate is imported from the cytosol; and (f) the glyoxylate cycle is inactive. In the aerobic regime, 75% of mitochondrial oxaloacetate arises from anaplerotic carboxylation of pyruvate, while in the anaerobic regime, the tricarboxylic acid cycle is operating in a branched fashion to fulfill biosynthetic demands only. The present study shows that fractional (13)C labeling of amino acids represents a powerful approach to study compartmented eukaryotic systems.     

Comparative study on central metabolic fluxes of <Emphasis Type="Italic">Bacillus megaterium</Emphasis> strains in continuous culture using <Superscript>13</Superscript>C labelled substrates

Fluxes of central carbon metabolism [glycolysis, pentose phosphate pathway (PPP), tricarboxylic acid cycle (TCA cycle), biomass formation] were determined for several Bacillus megaterium strains (DSM319, WH320, WH323, MS941) in C- and N-limited chemostat cultures by ¹³C labelling experiments. The labelling patterns of proteinogenic amino acids were analysed by GC/MS and therefrom flux ratios at important nodes within the metabolic network could be calculated. On the basis of a stoichiometric metabolic model flux distributions were estimated for the different B. megaterium strains used at various cultivation conditions. Generally all strains exhibited similar metabolic flux distributions, however, several significant changes were found in (1) the glucose flux entering the PPP via the oxidative branch, (2) the reversibilities within the PPP, (3) the relative fluxes of pyruvate and acetyl-CoA fed to the TCA cycle, (4) the fluxes around the pyruvate node involving a futile cycle.     

Determination of full 13C isotopomer distributions for metabolic flux analysis using heteronuclear spin echo difference NMR spectroscopy

13C-isotopomer labeling experiments play an increasingly important role in the analysis of intracellular metabolic fluxes for genetic engineering purposes. 13C NMR spectroscopy is a key technique in the experimental determination of isotopomer distributions. However, only subsets of isotopomers can be quantitated using this technique due to redundancies in the scalar coupling patterns and due to invisibility of the 12C isotope in NMR. Therefore, we developed and describe in this paper a 1H NMR spectroscopy method that allows to determine the complete isotopomer distribution in metabolites having a backbone consisting of up to at least four carbons. The proposed pulse sequences employ up to three alternately applied frequency-selective inversion pulses in the 13C channel. In a first application study, the complete isotopomer distribution of aspartate isolated from [1-13C]ethanol-grown Ashbya gossypii was determined. A tentative model of the central metabolism of this organism was constructed and used for metabolic flux analysis. The aspartate isotopomer NMR data played a key role in the successful determination of the flux through the peroxisomal glyoxylate pathway. The new NMR method can be highly instrumental in generating the data upon which isotopomer labeling experiments for flux analysis, that are becoming increasingly important, are based.     

Dynamic flux cartography of hairy roots primary metabolism

A dynamic model for plant cell and hairy root primary metabolism is presented. The model includes nutrient uptake (Pi, sugars, nitrogen sources), the glycolysis and pentose phosphate pathways, the TCA cycle, amino acid biosynthesis, respiratory chain, biosynthesis of cell building blocks (structural hexoses, organic acids, lipids, and organic phosphated molecules). The energy shuttles (ATP, ADP) and cofactors (NAD/H, NADP/H) are also included. The model describes the kinetics of 44 biochemical reactions (fluxes) of the primary metabolism of plant cells and includes 41 biochemical species (metabolites, nutrients, biomass components). Multiple Michaelis-Menten type kinetics are used to describe biochemical reaction rates. Known regulatory phenomena on metabolic pathways are included using sigmoid switch functions. A visualization framework showing fluxes and metabolite concentrations over time is presented. The visualization of fluxes and metabolites is used to analyze simulation results from Catharanthus roseus hairy root 50 d batch cultures. The visualization of the metabolic system allows analyzing split ratios between pathways and flux time-variations. For carbon metabolism, the cells were observed to have relatively high and stable fluxes for the central carbon metabolism and low and variable fluxes for anabolic pathways. For phosphate metabolism, a very high free intracellular Pi turnover rate was observed with higher flux variations than for the carbon metabolism. Nitrogen metabolism also exhibited large flux variations. The potential uses of the model are also discussed.     

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.     

Determination of the fluxes in the central metabolism of Corynebacterium glutamicum by nuclear magnetic resonance spectroscopy combined with metabolite balancing

To determine the in vivo fluxes of the central metabolism we have developed a comprehensive approach exclusively based on the fundamental enzyme reactions known to be present, the fate of the carbon atoms of individual reactions, and the metabolite balance of the culture. No information on the energy balance is required, nor information on enzyme activities, or the directionalities of reactions. Our approach combines the power of (1)H-detected (13)C nuclear magnetic resonance spectroscopy to follow individual carbons with the simplicity of establishing carbon balances of bacterial cultures. We grew a lysine-producing strain of Corynebacterium glutamicum to the metabolic and isotopic steady state with [1-(13)C]glucose and determined the fractional enrichments in 27 carbon atoms of 11 amino acids isolated from the cell. Since precursor metabolites of the central metabolism are incorporated in an exactly defined manner in the carbon skeleton of amino acids, the fractional enrichments in carbons of precursor metabolites (oxaloacetate, glyceraldehyde 3-phosphate, erythrose 4-phosphate, etc.) became directly accessible. A concise and generally applicable mathematical model was established using matrix calculus to express all metabolite mass and carbon labeling balances. An appropriate all-purpose software for the iterative solution of the equations is supplied. Applying this comprehensive methodology to C. glutamicum, all major fluxes within the central metabolism were determined. The result is that the flux through the pentose phosphate pathway is 66.4% (relative to the glucose input flux of 1.49 mmol/g dry weight h), that of entry into the tricarboxylic acid cycle 62.2%, and the contribution of the succinylase pathway of lysine synthesis 13.7%. Due to the large amount and high quality of measured data in vivo exchange reactions could also be quantitated with particularly high exchange rates within the pentose phosphate pathway for the ribose 5-phosphate transketolase reaction. Moreover, the total net flux of the anaplerotic reactions was quantitated as 38.0%. Most importantly, we found that in vivo one component within these anaplerotic reactions is a back flux from the carbon 4 units of the tricarboxylic acid cycle to the carbon 3 units of glycolysis of 30.6%. (c) 1996 John Wiley & Sons, Inc.     

Staphylococcus aureus physiological growth limitations: insights from flux calculations built on proteomics and external metabolite data

Comparing proteomics and metabolomics allows insights into Staphylococcus aureus physiological growth. We update genome and proteome information and deliver strain-specific metabolic models for three S. aureus strains (COL, N315, and Newman). We find a number of differences in metabolism and enzymes. Growth experiments (glucose or combined with oxygen limitation) were conducted to measure external metabolites. Fluxes of the central metabolism were calculated from these data with low error. In exponential phase, glycolysis is active and amino acids are used for growth. In later phases, dehydroquinate synthetase is suppressed and acetate metabolism starts. There are strain-specific differences for these phases. A time series of 2-D gel protein expression data on COL strain delivered a second data set (glucose limitation) on which fluxes were calculated. The comparison with the metabolite-predicted fluxes shows, in general, good correlation. Outliers point to different regulated enzymes for S. aureus COL under these limitations. In exponential growth, there is lower activity for some enzymes in upper glycolysis and pentose phosphate pathway and stronger activity for some in lower glycolysis. In transition phase, aspartate kinase is expressed to meet amino acid requirements and in later phases there is high expression of glyceraldehyde-3-phosphate dehydrogenase and lysine synthetase. Central metabolite fluxes and protein expression of their enzymes correlate in S. aureus.     

Dissection of central carbon metabolism of hemoglobin-expressing Escherichia coli by 13C nuclear magnetic resonance flux distribution analysis in microaerobic bioprocesses

Escherichia coli MG1655 cells expressing Vitreoscilla hemoglobin (VHb), Alcaligenes eutrophus flavohemoprotein (FHP), the N-terminal hemoglobin domain of FHP (FHPg), and a fusion protein which comprises VHb and the A. eutrophus C-terminal reductase domain (VHb-Red) were grown in a microaerobic bioreactor to study the effects of low oxygen concentrations on the central carbon metabolism, using fractional (13)C-labeling of the proteinogenic amino acids and two-dimensional [(13)C, (1)H]-correlation nuclear magnetic resonance (NMR) spectroscopy. The NMR data revealed differences in the intracellular carbon fluxes between E. coli cells expressing either VHb or VHb-Red and cells expressing A. eutrophus FHP or the truncated heme domain (FHPg). E. coli MG1655 cells expressing either VHb or VHb-Red were found to function with a branched tricarboxylic acid (TCA) cycle. Furthermore, cellular demands for ATP and reduction equivalents in VHb- and VHb-Red-expressing cells were met by an increased flux through glycolysis. In contrast, in E. coli cells expressing A. eutrophus hemeproteins, the TCA cycle is running cyclically, indicating a shift towards a more aerobic regulation. Consistently, E. coli cells displaying FHP and FHPg activity showed lower production of the typical anaerobic by-products formate, acetate, and D-lactate. The implications of these observations for biotechnological applications are discussed.     

A flux model of glycolysis and the oxidative pentosephosphate pathway in developing Brassica napus embryos

Developing oilseeds synthesize large quantities of triacylglycerol from sucrose and hexose. To understand the fluxes involved in this conversion, a quantitative metabolic flux model was developed and tested for the reaction network of glycolysis and the oxidative pentose phosphate pathway (OPPP). Developing Brassica napus embryos were cultured with [U-13C6]glucose, [1-13C]glucose, [6-13C]glucose, [U-13C12]sucrose, and/or [1,2-13C2]glucose and the labeling patterns in amino acids, lipids, sucrose, and starch were measured by gas chromatography/mass spectrometry and NMR. Data were used to verify a reaction network of central carbon metabolism distributed between the cytosol and plastid. Computer simulation of the steady state distribution of isotopomers in intermediates of the glycolysis/OPPP network was used to fit metabolic flux parameters to the experimental data. The observed distribution of label in cytosolic and plastidic metabolites indicated that key intermediates of glycolysis and OPPP have similar labeling in these two compartments, suggesting rapid exchange of metabolites between these compartments compared with net fluxes into end products. Cycling between hexose phosphate and triose phosphate and reversible transketolase velocity were similar to net glycolytic flux, whereas reversible transaldolase velocity was minimal. Flux parameters were overdetermined by analyzing labeling in different metabolites and by using data from different labeling experiments, which increased the reliability of the findings. Net flux of glucose through the OPPP accounts for close to 10% of the total hexose influx into the embryo. Therefore, the reductant produced by the OPPP accounts for at most 44% of the NADPH and 22% of total reductant needed for fatty acid synthesis.