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.     

Genetic characterization of a single bifunctional enzyme for fumarate reduction and succinate oxidation in Geobacter sulfurreducens and engineering of fumarate reduction in Geobacter metallireducens

The mechanism of fumarate reduction in Geobacter sulfurreducens was investigated. The genome contained genes encoding a heterotrimeric fumarate reductase, FrdCAB, with homology to the fumarate reductase of Wolinella succinogenes and the succinate dehydrogenase of Bacillus subtilis. Mutation of the putative catalytic subunit of the enzyme resulted in a strain that lacked fumarate reductase activity and was unable to grow with fumarate as the terminal electron acceptor. The mutant strain also lacked succinate dehydrogenase activity and did not grow with acetate as the electron donor and Fe(III) as the electron acceptor. The mutant strain could grow with acetate as the electron donor and Fe(III) as the electron acceptor if fumarate was provided to alleviate the need for succinate dehydrogenase activity in the tricarboxylic acid cycle. The growth rate of the mutant strain under these conditions was faster and the cell yields were higher than for wild type grown under conditions requiring succinate dehydrogenase activity, suggesting that the succinate dehydrogenase reaction consumes energy. An orthologous frdCAB operon was present in Geobacter metallireducens, which cannot grow with fumarate as the terminal electron acceptor. When a putative dicarboxylic acid transporter from G. sulfurreducens was expressed in G. metallireducens, growth with fumarate as the sole electron acceptor was possible. These results demonstrate that, unlike previously described organisms, G. sulfurreducens and possibly G. metallireducens use the same enzyme for both fumarate reduction and succinate oxidation in vivo.     

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.     

Pathway confirmation and flux analysis of central metabolic pathways in Desulfovibrio vulgaris hildenborough using gas chromatography-mass spectrometry and Fourier transform-ion cyclotron resonance mass spectrometry

Flux distribution in central metabolic pathways of Desulfovibrio vulgaris Hildenborough was examined using 13C tracer experiments. Consistent with the current genome annotation and independent evidence from enzyme activity assays, the isotopomer results from both gas chromatography-mass spectrometry (GC-MS) and Fourier transform-ion cyclotron resonance mass spectrometry (FT-ICR MS) indicate the lack of an oxidatively functional tricarboxylic acid (TCA) cycle and an incomplete pentose phosphate pathway. Results from this study suggest that fluxes through both pathways are limited to biosynthesis. The data also indicate that >80% of the lactate was converted to acetate and that the reactions involved are the primary route of energy production [NAD(P)H and ATP production]. Independently of the TCA cycle, direct cleavage of acetyl coenzyme A to CO and 5,10-methyl tetrahydrofuran also leads to production of NADH and ATP. Although the genome annotation implicates a ferredoxin-dependent oxoglutarate synthase, isotopic evidence does not support flux through this reaction in either the oxidative or the reductive mode; therefore, the TCA cycle is incomplete. FT-ICR MS was used to locate the labeled carbon distribution in aspartate and glutamate and confirmed the presence of an atypical enzyme for citrate formation suggested in previous reports [the citrate synthesized by this enzyme is the isotopic antipode of the citrate synthesized by the (S)-citrate synthase]. These findings enable a better understanding of the relation between genome annotation and actual metabolic pathways in D. vulgaris and also demonstrate that FT-ICR MS is a powerful tool for isotopomer analysis, overcoming the problems with both GC-MS and nuclear magnetic resonance spectroscopy.     

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.     

Metabolic flux analysis with a comprehensive isotopomer model in Bacillus subtilis

Fluxes in central carbon metabolism of a genetically engineered, riboflavin-producing Bacillus subtilis strain were investigated in glucose-limited chemostat cultures at low (0.11 h(-1)) and high (0.44 h(-1)) dilution rates. Using a mixture of 10% [U-(13)C] and 90% glucose labeled at natural abundance, (13)C-labeling experiments were carried out to provide additional information for metabolic flux balancing. The resulting labeling pattern in the proteinogenic amino acids were analyzed by two-dimensional [(13)C, (1)H] nuclear magnetic resonance (NMR) spectroscopy. To account rigorously for all available data from these experiments, we developed a comprehensive isotopomer model of B. subtilis central metabolism. Using this model, intracellular carbon net and exchange fluxes were estimated on the basis of validated physiological data and biomass composition in combination with 2D NMR data from 45 individual carbon atom spectra in the amino acids. Glucose catabolism proceeded primarily via glycolysis but pentose phosphate pathway fluxes increased with increasing growth rate. Moreover, significant back fluxes from the TCA cycle to the lower part of glycolysis via the gluconeogenic PEP carboxykinase were detected. The malic enzyme reaction, in contrast, was found to be inactive. A thorough statistical analysis was performed to prove the reliability of the isotopomer balance model and the obtained results. Specifically, a chi(2) test was applied to validate the model and the chi-square criterion was used to explore the sensitivity of model predictions to the experimental data.     

Unraveling the metabolism of HEK-293 cells using lactate isotopomer analysis

HEK-293 is the most extensively used human cell line for the production of viral vectors and is gaining increasing attention for the production of recombinant proteins by transient transfection. To further improve the metabolic characterization of this cell line, we have performed cultures using ¹³C-labeled substrates and measured the resulting mass isotopomer distributions in lactate by LC/MS. Simultaneous metabolite and isotopomer balancing allowed improvement and validation of the metabolic model and quantification of key intracellular pathways. We have determined the amounts of glucose carbon channeled through the PPP, incorporated into the TCA cycle for energy production and lipids biosynthesis, as well as the cytosolic and mitochondrial malic enzyme fluxes. Our analysis also revealed that glutamine did not significantly contribute to lactate formation. An improved and quantitative understanding of the central carbon metabolism is greatly needed to pursue the rational development of engineering approaches at both the cellular and process levels.     

13C-NMR, 1H-NMR and gas-chromatography mass-spectrometry studies of the biosynthesis of 13C-enriched L-lysine by Brevibacterium flavum

The basic metabolic pathways of lysine biosynthesis in Brevibacterium flavum, a strain which excretes excessive amounts of L-lysine, have been followed by using two 13C-labeled precursors. 13C- and 1H-NMR spectroscopies in conjunction with gas chromatography mass spectrometry (GC-MS) have revealed the various metabolic pathways leading to L-[13C]lysine. Discrete metabolic pathways give rise to distinct labeling patterns. L-Lysine resulting from [1-13C]glucose fermentation is relatively specifically labeled: L-[3,5-13C]lysine is the main product. Experimental and theoretical approaches based on the 13C-enrichment values of intracellular glutamate, a major intermediate metabolite, allowed us to assess the relative contribution of the major metabolic pathways forming lysine. The labeling pattern of glutamate reflects the isotope distribution in 2-oxoglutarate. When [2-13C]acetate is used as the sole carbon source in the culture, the energy-producing steps of the Krebs cycle are essential. The higher activity of the Krebs cycle, when endogenous carbohydrates are exhausted from the culture, is indicated by the increased 13C enrichment in C-1 of lysine and reveal a high content of isotopomers of four, five and six 13C atoms in the lysine molecule, pointing out that the four-carbon intermediates of the cycle are being derived from the glyoxylate shunt pathway. Such a phenomenon does not occur in glucose fermentation. GC-MS analyses of 13C enrichments and isotopomer distributions in metabolites and end products are in good agreement with the predicted contribution of each metabolic pathway. This new methodological approach of combined NMR and GC-MS has been demonstrated to be applicable to various other metabolic studies.     

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.     

Application of metabolic flux analysis for the identification of metabolic bottlenecks in the biosynthesis of penicillin-G

A detailed stoichiometric model was developed for growth and penicillin-G production in Penicillium chrysogenum. From an a priori metabolic flux analysis using this model it appeared that penicillin production requires significant changes in fluxes through the primary metabolic pathways. This is brought about by the biosynthesis of carbon precursors for the beta-lactan nucleus and an increased demand for NADPH, mainly for sulfate reduction. As a result, significant changes in flux partitioning occur around four principal nodes in primary metabolism. These are located at: (1) glucose-6-phosphate; (2) 3-phosphoglycerate; (3) mitochondrial pyruvate; and (4) mitochondrial isocitrate. These nodes should be regarded as potential bottlenecks for increased productivity. The flexibility of these principal nodes was investigated by experimental manipulation of the fluxes through the central metabolic pathways using a high-producing strain of P. chrysogenum. Metabolic fluxes were manipulated through growth of the cells on different substrates in carbon-limited chemostat culture. Metabolic flux analysis, based on measured input and output fluxes, was used to calculate the fluxes around the principal nodes. It was found that, for growth on glucose, ethanol, and acetate, the flux partitioning around these nodes differed significantly. However, this had hardly any effect on penicillin productivity, showing that primary carbon metabolism is not likely to contain potential bottlenecks. Further experiments were performed to manipulate the total metabolic demand for the cofactor nicotinamide adenine dinucleotide phosphate (NADPH). NADPH demand was increased stepwise by cultivating the cells on glucose or xylose as the carbon source combined with either ammonia or nitrate as the nitrogen source, which resulted in a stepwise decrease of penicillin production. This clearly shows that, in penicillin fermentation, possible limitations in primary metabolism reside in the supply/regeneration of cofactors (NADPH) rather than in the supply of carbon precursors.     

Metabolic flux analysis of Escherichia coli K12 grown on 13C-labeled acetate and glucose using GC-MS and powerful flux calculation method

A new algorithm was developed for the estimation of the metabolic flux distribution based on GC-MS data of proteinogenic amino acids. By using a sensitive GC-MS protocol as well as by combining the global search algorithm such as the genetic algorithm with the local search algorithm such as the Levenberg-Marquardt algorithm, not only the distribution of the net fluxes in the entire network, but also certain exchange fluxes which contribute significantly to the isotopomer distribution could be quantified. This mass isotopomer analysis could identify the biochemical changes involved in the regulation where acetate or glucose was used as a main carbon source. The metabolic flux analysis clearly revealed that when the specific growth rate increased, only a slight change in flux distribution was observed for acetate metabolism, indicating that subtle regulation mechanism exists in certain key junctions of this network system. Different from acetate metabolism, when glucose was used as a carbon source, as the growth rate increased, a significant increase in relative pentose phosphate pathway (PPP) flux was observed for Escherichia coli K12 at the expense of the citric acid cycle, suggesting that when growing on glucose, the flux catalyzed by isocitrate dehydrogenase could not fully fulfill the NADPH demand for cell growth, causing the oxidative PPP to be utilized to a larger extent so as to complement the NADPH demand. The GC-MS protocol as well as the new algorithm demonstrated here proved to be a powerful tool for characterizing metabolic regulation and can be utilized for strain improvement and bioprocess optimization.     

<Superscript>13</Superscript>C metabolite profiling to compare the central metabolic flux in two yeast strains

¹³C metabolite profiling to quantify the dynamic changes of central carbon metabolites was attempted using mass isotopomer distribution analysis in two yeast strains, Saccharomyces cerevisiae and Kluyveromyces marxianus. Mass and isotopomer balances of the intermediates were examined and calculated in both yeast species and central carbon metabolic fluxes were successfully determined. Metabolic fluxes of pentose phosphate pathway in K. marxianus were 1.66 times higher than S. cerevisiae. The flux difference was also supported by relatively high abundance of partially labeled fructose 6-phosphate and 3-phosphoglycerate as well as an increased concentration of labeled L-valine in K. marxianus. Metabolic flux analysis combined with dynamic metabolite profiling has provided better understanding in the central carbon metabolic pathways of two model organisms and can be applied as a method to analyze more complicated metabolic networks in other organisms.     

The relationships between the metabolic fluxes and 13C-labeled isotopomer distribution for the flux analysis of the main metabolic pathways

It has been known that ¹³C-labeling technique is quite useful in estimating the metabolic fluxes. Although the program-based flux analysis is powerful, it is not easy to be confident with the result obtained without experiences and exhaustive trial and errors based on statistical analysis for the confidence intervals in practice. It is, therefore, quite important to grasp the relationship between the fluxes and the ¹³C-labeled isotopomer distribution to get deeper insight into the metabolic flux analysis. In the present research, it was shown explicitly how the isotopomer distribution changes with respect to the fluxes in relation to the labeling patterns of the substrate, where either labeled glucose, acetate, or pyruvate was used as a carbon source. Some of the analytical expressions were derived based on the matrix representation, and they were utilized for analysis. It was shown that the isotopomer pattern does not necessarily change uniformly with respect to fluxes, but changes in a complicated way in particular for the case of using pyruvate as a carbon source where some isotopomers do not necessarily change monotonically. It was shown to be quite important to grasp how the isotopomer pattern changes with respect to fluxes and the labeling pattern of the substrate for flux determination and the experimental design. It was also shown that the mixture of [1-¹³C] acetate and [2-¹³C] acetate should not be used from the information index point of view. Some of the experimental data were evaluated from the present approach. It was also shown that the isotopomer distribution is less sensitive to the bidirectional fluxes in the reversible pathway.     

13C and 1H nuclear magnetic resonance study of glycogen futile cycling in strains of the genus Fibrobacter

We investigated the carbon metabolism of three strains of Fibrobacter succinogenes and one strain of Fibrobacter intestinalis. The four strains produced the same amounts of the metabolites succinate, acetate, and formate in approximately the same ratio (3.7/1/0.3). The four strains similarly stored glycogen during all growth phases, and the glycogen-to-protein ratio was close to 0.6 during the exponential growth phase. 13C nuclear magnetic resonance (NMR) analysis of [1-13C]glucose utilization by resting cells of the four strains revealed a reversal of glycolysis at the triose phosphate level and the same metabolic pathways. Glycogen futile cycling was demonstrated by 13C NMR by following the simultaneous metabolism of labeled [13C]glycogen and exogenous unlabeled glucose. The isotopic dilutions of the CH2 of succinate and the CH3 of acetate when the resting cells were metabolizing [1-13C]glucose and unlabeled glycogen were precisely quantified by using 13C-filtered spin-echo difference 1H NMR spectroscopy. The measured isotopic dilutions were not the same for succinate and acetate; in the case of succinate, the dilutions reflected only the contribution of glycogen futile cycling, while in the case of acetate, another mechanism was also involved. Results obtained in complementary experiments are consistent with reversal of the succinate synthesis pathway. Our results indicated that for all of the strains, from 12 to 16% of the glucose entering the metabolic pathway originated from prestored glycogen. Although genetically diverse, the four Fibrobacter strains studied had very similar carbon metabolism characteristics.     

Determination of metabolic fluxes in a non-steady-state system

Estimation of fluxes through metabolic networks from redistribution patterns of (13)C has become a well developed technique in recent years. However, the approach is currently limited to systems at metabolic steady-state; dynamic changes in metabolic fluxes cannot be assessed. This is a major impediment to understanding the behaviour of metabolic networks, because steady-state is not always experimentally achievable and a great deal of information about the control hierarchy of the network can be derived from the analysis of flux dynamics. To address this issue, we have developed a method for estimating non-steady-state fluxes based on the mass-balance of mass isotopomers. This approach allows multiple mass-balance equations to be written for the change in labelling of a given metabolite pool and thereby permits over-determination of fluxes. We demonstrate how linear regression methods can be used to estimate non-steady-state fluxes from these mass balance equations. The approach can be used to calculate fluxes from both mass isotopomer and positional isotopomer labelling information and thus has general applicability to data generated from common spectrometry- or NMR-based analytical platforms. The approach is applied to a GC-MS time-series dataset of (13)C-labelling of metabolites in a heterotrophic Arabidopsis cell suspension culture. Threonine biosynthesis is used to demonstrate that non-steady-state fluxes can be successfully estimated from such data while organic acid metabolism is used to highlight some common issues that can complicate flux estimation. These include multiple pools of the same metabolite that label at different rates and carbon skeleton rearrangements.     

d-Pentose metabolism byBacteroides vulgatus strain 8482

Bacteroides vulgatus strain 8482 metabolizedd-arabinose by a mechanism involving a 32 (top to bottom) cleavage of the arabinose carbon skeleton. During growth in the presence of 1-¹⁴C-d-arabinose, acetate, propionate, and succinate were labeled, but during growth in the presence of 5-labeledd-arabinose, only labeled acetate and succinate were formed. The metabolism ofd-ribose by strain 8482 differed from that ford-arabinose. Strain 8482 converted glycolic acid and glycine to acetate and succinate, but not propionate, by a mechanism involving cleavage of the glycine and glycolic acid carbon skeletons and equilibration of carbons 1 and 2 of glycolic acid and glycine with nonequivalent metabolic pools. The metabolism ofd-arabinose,d-ribose,d-glycine, andd-glycolic acid by strain 8482 was similar, in some respects, to that ofBacteroides fragilis strain 2044, but differed substantially from the metabolism of the same substances byBacteroides ruminicola strain B₁4.     

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.     

Experimental identification and quantification of glucose metabolism in seven bacterial species

The structurally conserved and ubiquitous pathways of central carbon metabolism provide building blocks and cofactors for the biosynthesis of cellular macromolecules. The relative uses of pathways and reactions, however, vary widely among species and depend upon conditions, and some are not used at all. Here we identify the network topology of glucose metabolism and its in vivo operation by quantification of intracellular carbon fluxes from 13C tracer experiments. Specifically, we investigated Agrobacterium tumefaciens, two pseudomonads, Sinorhizobium meliloti, Rhodobacter sphaeroides, Zymomonas mobilis, and Paracoccus versutus, which grow on glucose as the sole carbon source, represent fundamentally different metabolic lifestyles (aerobic, anaerobic, photoheterotrophic, and chemoheterotrophic), and are phylogenetically distinct (firmicutes, gamma-proteobacteria, and alpha-proteobacteria). Compared to those of the model bacteria Escherichia coli and Bacillus subtilis, metabolisms of the investigated species differed significantly in several respects: (i) the Entner-Doudoroff pathway was the almost exclusive catabolic route; (ii) the pentose phosphate pathway exhibited exclusively biosynthetic functions, in many cases also requiring flux through the nonoxidative branch; (iii) all aerobes exhibited fully respiratory metabolism without significant overflow metabolism; and (iv) all aerobes used the pyruvate bypass of the malate dehydrogenase reaction to a significant extent. Exclusively, Pseudomonas fluorescens converted most glucose extracellularly to gluconate and 2-ketogluconate. Overall, the results suggest that metabolic data from model species with extensive industrial and laboratory history are not representative of microbial metabolism, at least not quantitatively.