The intracellular carbon flux distribution in wild-type and pyruvate kinase-deficient Escherichia coli was estimated using biosynthetically directed fractional 13 C labeling experiments with [U-13 C 6 ]glucose in glucoseor ammonia-limited chemostats, two-dimensional nuclear magnetic resonance (NMR) spectroscopy of cellular amino acids, and a comprehensive isotopomer model. The general response to disruption of both pyruvate kinase isoenzymes in E. coli was a local flux rerouting via the combined reactions of phosphoenolpyruvate (PEP) carboxylase and malic enzyme. Responses in the pentose phosphate pathway and the tricarboxylic acid cycle were strongly dependent on the environmental conditions. In addition, high futile cycling activity via the gluconeogenic PEP carboxykinase was identified at a low dilution rate in glucose-limited chemostat culture of pyruvate kinase-deficient E. coli, with a turnover that is comparable to the specific glucose uptake rate. Furthermore, flux analysis in mutant cultures indicates that glucose uptake in E. coli is not catalyzed exclusively by the phosphotransferase system in glucose-limited cultures at a low dilution rate. Reliability of the flux estimates thus obtained was verified by statistical error analysis and by comparison to intracellular carbon flux ratios that were independently calculated from the same NMR data by metabolic flux ratio analysis.The central carbon pathways fulfill anabolic and catabolic functions by providing cofactors and building blocks for macromolecular synthesis as well as energy. While some singlegene knockout mutations in central metabolism preclude growth on glucose, a majority of such variations can potentially be compensated for either by isoenzymes (39) or by a rerouting of carbon fluxes through alternative pathways (7,15,20,42). Usually, the phenotypes of knockout mutants is characterized by quantitative physiological analysis, and conclusions on intracellular metabolism are then based on qualitative interpretation of these results. To reveal cause and effect relationships, however, it is important to gain deeper insight into the complex metabolic responses at the level of intracellular metabolite concentrations and fluxes. These intracellular carbon fluxes, or in vivo reaction rates, are per se nonmeasurable quantities that cannot usually be inferred directly from in vitro enzyme activities because not all in vivo effector concentrations are known. Hence, the intracellular reaction rates are commonly estimated by methods of metabolic flux analysis, which provide a holistic perspective on metabolism (52, 54).The most common approach is based on flux balancing analysis within a stoichiometric model of cellular metabolism (26,41, 50, 52). In this approach, uptake and secretion rates and biosynthetic requirements are balanced in a stoichiometric model, assuming quasi-steady-state material balances on the intermediates. In many cases, however, limited extracellular data and stoichiometric constraints lead to underdetermined systems, so that addition...

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.

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.

The energetic efficiency of microbial growth is significantly reduced in cultures growing under glucose excess compared to cultures growing under glucose limitation, but the magnitude to which different energy-dissipating processes contribute to the reduced efficiency is currently not well understood. We introduce here a new concept for balancing the total cellular energy flux that is based on the conversion of energy and carbon fluxes into energy equivalents, and we apply this concept to glucose-, ammonia-, and phosphate-limited chemostat cultures of riboflavin-producing Bacillus subtilis. Based on [U-13 C 6 ]glucose-labeling experiments and metabolic flux analysis, the total energy flux in slow-growing, glucose-limited B. subtilis is almost exclusively partitioned in maintenance metabolism and biomass formation. In excess-glucose cultures, in contrast, uncoupling of anabolism and catabolism is primarily achieved by overflow metabolism, while two quantified futile enzyme cycles and metabolic shifts to energetically less efficient pathways are negligible. In most cultures, about 20% of the total energy flux could not be assigned to a particular energy-consuming process and thus are probably dissipated by processes such as ion leakage that are not being considered at present. In contrast to glucoseor ammonia-limited cultures, metabolic flux analysis revealed low tricarboxylic acid (TCA) cycle fluxes in phosphate-limited B. subtilis, which is consistent with CcpA-dependent catabolite repression of the cycle and/or transcriptional activation of genes involved in overflow metabolism in the presence of excess glucose. ATPdependent control of in vivo enzyme activity appears to be irrelevant for the observed differences in TCA cycle fluxes.The very basis of microbial growth resides in balanced fluxes through anabolic and catabolic reactions. These metabolic fluxes are highly variable and change with the environmental conditions and the rate of growth, since faster-growing cells demand a higher rate of metabolism. To delineate these influences, metabolic flux responses are typically studied in chemostat cultures that are maintained under different nutrient limitations. When microorganisms are limited for their energy source (usually the carbon source), catabolism is tightly coupled to anabolism and high biomass yields on the carbon source are achieved (40). Compared to those with carbon (C) limitation, excess-C cultures exhibit generally high rates of carbon consumption and low yields of biomass and thus have a low energetic growth efficiency (11,31,32,57). Most frequently excess-C cultures in chemostats are limited by other cellular macroelements such as nitrogen (N) and phosphorus (P) but also potassium, the predominant intracellular cation (11). Potassium or P limitation invokes usually a strong uncoupling of catabolic and anabolic processes that lead consequently to low biomass yields on the energy source, while N limitation has more moderate effects on cellular physiology (31, 68). Although mostly studied in the gr...

13C metabolic flux analysis (MFA) is based on carbon-labeling experiments where a specifically (13)C labeled substrate is fed. The labeled carbon atoms distribute over the metabolic network and the label enrichment of certain metabolic pools is measured by using different methods. Recently, MS methods have been dramatically improved-large and precise datasets are now available. MS data has to be preprocessed and corrected for natural stable mass isotopes. In this article we present (1). a new elegant method to correct MS measurement data for natural stable mass isotopes by infinite dimensional matrix calculus and (2). we statistically analyze and discuss a reconstruction of labeling pattern in metabolic precursors from biosynthesis molecules. Moreover, we establish a new method for consistency checking of MS spectra that can be applied for automatic error recognition in high-throughput flux analysis procedures. Preprocessing the measurement data changes their statistical properties which have to be considered in the subsequent parameter fitting process for (13)C MFA. We show that correcting for stable mass isotopes leads to rather small correlations. On the other hand, a direct reconstruction of a precursor labeling pattern from an aromatic amino acid measurement turns out to be critical. Reasonable results are only obtained if additional, independent information about the labeling of at least one precursor is available. A versatile MatLab tool for the rapid correction and consistency checking of MS spectra is presented. Practical examples for the described methods are also given.

Rate equations for measured extracellular rates and macromolecular composition data were combined with a stoichiometric model to describe riboflavin production with an industrial Bacillus subtilis strain using errors in variables regression analysis. On the basis of this combined stoichiometric growth model, we explored the topological features of the B. subtilis metabolic reaction network that was assembled from a large amount of literature. More specifically, we simulated maximum theoretical yields of biomass and riboflavin, including the associated flux regimes. Based on the developed model, the importance of experimental data on building block requirements for maximum yield and flux calculations were investigated. These analyses clearly show that verification of macromolecular composition data is important for optimum flux calculations.

Metabolome analysis, the analysis of large sets of intracellular metabolites, has become an important systems analysis method in biotechnological and pharmaceutical research. In metabolic engineering, the integration of metabolome data with fluxome and proteome data into large-scale mathematical models promises to foster rational strategies for strain and cell line improvement. However, the development of reproducible sampling procedures for quantitative analysis of intracellular metabolite concentrations represents a major challenge, accomplishing (i) fast transfer of sample, (ii) efficient quenching of metabolism, (iii) quantitative metabolite extraction, and (iv) optimum sample conditioning for subsequent quantitative analysis. In addressing these requirements, we propose an integrated sampling procedure. Simultaneous quenching and quantitative extraction of intracellular metabolites were realized by short-time exposure of cells to temperatures < or =95 degrees C, where intracellular metabolites are released quantitatively. Based on these findings, we combined principles of heat transfer with knowledge on physiology, for example, turnover rates of energy metabolites, to develop an optimized sampling procedure based on a coiled single tube heat exchanger. As a result, this sampling procedure enables reliable and reproducible measurements through (i) the integration of three unit operations into a one unit operation, (ii) the avoidance of any alteration of the sample due to chemical reagents in quenching and extraction, and (iii) automation. A sampling frequency of 5 s(-)(1) and an overall individual sample processing time faster than 30 s allow observing responses of intracellular metabolite concentrations to extracellular stimuli on a subsecond time scale. Recovery and reliability of the unit operations were analyzed. Impact of sample conditioning on subsequent IC-MS analysis of metabolites was examined as well. The integrated sampling procedure was validated through consistent results from steady-state metabolite analysis of Escherichia coli cultivated in a chemostat at D = 0.1 h(-)(1).

In mammalian cell cultures, ammonia that is released into the medium as a result of glutamine metabolism and lactate that is excreted due to incomplete glucose oxidation are both known to essentially inhibit the growth of cells. For some cell lines, for example, hybridoma cells, excreted ammonia also has an effect on product formation. Although glutamine has been generally considered as the major energy source for mammalian cells, it was recently found that various adherent cell lines (MDCK, CHO-K1, and BHK21) can grow as well in glutamine-free medium, provided glutamine is substituted with pyruvate. In such a medium the level of both ammonia and lactate released was significantly reduced. In this study, metabolic flux analysis (MFA) was applied to Madin Darby Canine Kidney (MDCK) cells cultivated in glutamine-containing and glutamine-free medium. The results of the MFA allowed further investigation of the influence of glutamine substitution with pyruvate on the metabolism of MDCK cells during different growth stages of adherent cells, e.g., early exponential and late contact-inhibited phase. Pyruvate seemed to directly enter the TCA cycle, whereas most of the glucose consumed was excreted as lactate. Although the exact mechanisms are not clear so far, this resulted in a reduction of the glucose uptake necessary for cellular metabolism in glutamine-free medium. Furthermore, consumption of ATP by futile cycles seemed to be significantly reduced when substituting glutamine with pyruvate. These findings imply that glutamine-free medium favors a more efficient use of nutrients by cells. However, a number of metabolic fluxes were similar in the two cultivations considered, e.g., most of the amino acid uptake and degradation rates or fluxes through the branch of the TCA cycle converting R-ketoglutarate to malate, which is responsible for the mitochondrial ATP synthesis. Besides, the specific rate of cell growth was approximately the same in both cultivations. Thus, the switch from glutamine-containing to glutamine-free medium with pyruvate provided a series of benefits without dramatic changes of cellular metabolism.