Metabolic Adaptation Processes That Converge to Optimal Biomass Flux Distributions

In simple organisms like E.coli, the metabolic response to an external perturbation passes through a transient phase in which the activation of a number of latent pathways can guarantee survival at the expenses of growth. Growth is gradually recovered as the organism adapts to the new condition. This adaptation can be modeled as a process of repeated metabolic adjustments obtained through the resilencings of the non-essential metabolic reactions, using growth rate as selection probability for the phenotypes obtained. The resulting metabolic adaptation process tends naturally to steer the metabolic fluxes towards high growth phenotypes. Quite remarkably, when applied to the central carbon metabolism of E.coli, it follows that nearly all flux distributions converge to the flux vector representing optimal growth, i.e., the solution of the biomass optimization problem turns out to be the dominant attractor of the metabolic adaptation process.     

不同浓度外源乙酸钠对耐乙酸大肠埃希菌DA19代谢和关键酶酶活的影响

为研究外源乙酸钠对大肠埃希菌DA19生长代谢的影响,将该菌株在氮源限制基本培养基及添加不同浓度乙酸钠的氮源限制基本培养基中连续培养,测定稳态时生长代谢参数和胞内关键酶酶活。与MN培养基相比,葡萄糖比消耗速率和延胡索酸比生成速率随外源乙酸钠质量浓度增加而逐渐下降,丙酮酸比生成速率则随外源乙酸钠质量浓度增加而明显增加,而乙酸比生成速率则明显降低(除9 g/L乙酸钠外)。磷酸果糖激酶、异柠檬酸脱氢酶、异柠檬酸裂解酶、苹果酸脱氢酶、磷酸烯醇式丙酮酸羧化酶和乙酸激酶酶活随外源乙酸钠质量浓度增加而呈先下降后上升的趋势,而6-磷酸葡萄糖脱氢酶则随着外源乙酸钠质量浓度增加而逐渐降低。为了应对外源乙酸钠压力,大肠埃希菌DA19的生长代谢和中心代谢途径酶活都发生了明显改变。     

Effects of Growth Mode and Pyruvate Carboxylase on Succinic Acid Production by Metabolically Engineered Strains of Escherichia coli

Escherichia coli NZN111, which lacks activities for pyruvate-formate lyase and lactate dehydrogenase, and AFP111, a derivative which contains an additional mutation in ptsG (a gene encoding an enzyme of the glucose phophotransferase system), accumulate significant levels of succinic acid (succinate) under anaerobic conditions. Plasmid pTrc99A-pyc, which expresses the Rhizobium etli pyruvate carboxylase enzyme, was introduced into both strains. We compared growth, substrate consumption, product formation, and activities of seven key enzymes (acetate kinase, fumarate reductase, glucokinase, isocitrate dehydrogenase, isocitrate lyase, phosphoenolpyruvate carboxylase, and pyruvate carboxylase) from glucose for NZN111, NZN111/pTrc99A-pyc, AFP111, and AFP111/pTrc99A-pyc under both exclusively anaerobic and dual-phase conditions (an aerobic growth phase followed by an anaerobic production phase). The highest succinate mass yield was attained with AFP111/pTrc99A-pyc under dual-phase conditions with low pyruvate carboxylase activity. Dual-phase conditions led to significant isocitrate lyase activity in both NZN111 and AFP111, while under exclusively anaerobic conditions, an absence of isocitrate lyase activity resulted in significant pyruvate accumulation. Enzyme assays indicated that under dual-phase conditions, carbon flows not only through the reductive arm of the tricarboxylic acid cycle for succinate generation but also through the glyoxylate shunt and thus provides the cells with metabolic flexibility in the formation of succinate. Significant glucokinase activity in AFP111 compared to NZN111 similarly permits increased metabolic flexibility of AFP111. The differences between the strains and the benefit of pyruvate carboxylase under both exclusively anaerobic and dual-phase conditions are discussed in light of the cellular constraint for a redox balance.     

Citrate production from n-alkane by Candida lipolytica in reference to carbon fluxes in vivo

Summary Citrate production by Candida lipolytica from n-alkane as sole carbon source was studied with a NH ₄ ⁺ -limited chemostat culture. This organism produced citrate and isocitrate simultaneously, the ratio of these acids remaining nearly unchanged within the range of dilution rates used. However, the specific production rates of these acids were maximal around a particular dilution rate. For each steady-state culture, specific activities of some key enzymes in tricarboxylate and glyoxylate cycles were measured in vitro. In order to correlate the enzyme activites with the change in specific production rate of each acid, the intracellular carbon fluxes were assessed by solving algebraic equations derived from carbon balance with respect to a simplified metabolic map of n-alkane metabolism. By comparing the fluxes assessed in vivo with the enzyme activities measured in vitro, the metabolic pathway(s) most likely to control the production of acids and cell material were investigated.     

Metabolic Analysis of Escherichia coli in the Presence and Absence of the Carboxylating Enzymes Phosphoenolpyruvate Carboxylase and Pyruvate Carboxylase

Fermentation patterns of Escherichia coli with and without the phosphoenolpyruvate carboxylase (PPC) and pyruvate carboxylase (PYC) enzymes were compared under anaerobic conditions with glucose as a carbon source. Time profiles of glucose and fermentation product concentrations were determined and used to calculate metabolic fluxes through central carbon pathways during exponential cell growth. The presence of the Rhizobium etli pyc gene in E. coli (JCL1242/pTrc99A-pyc) restored the succinate producing ability of E. coli ppc null mutants (JCL1242), with PYC competing favorably with both pyruvate formate lyase and lactate dehydrogenase. Succinate formation was slightly greater by JCL1242/pTrc99A-pyc than by cells which overproduced PPC (JCL1242/pPC201, ppc⁺), even though PPC activity in cell extracts of JCL1242/pPC201 (ppc⁺) was 40-fold greater than PYC activity in extracts of JCL1242/pTrc99a-pyc. Flux calculations indicate that during anaerobic metabolism the pyc⁺ strain had a 34% greater specific glucose consumption rate, a 37% greater specific rate of ATP formation, and a 6% greater specific growth rate compared to the ppc⁺ strain. In light of the important position of pyruvate at the juncture of NADH-generating pathways and NADH-dissimilating branches, the results show that when PPC or PYC is expressed, the metabolic network adapts by altering the flux to lactate and the molar ratio of ethanol to acetate formation.     

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.     

Nodule carbohydrate metabolism and polyols involvement in the response of Medicago sativa to salt stress

Alterations of plant growth, chlorophyll fluorescence parameters, nodule carbon metabolism and polyols concentration as result of salt stress were examined in alfalfa (Medicago sativa). Plants, in symbiosis with Sinorhizobium meliloti GR4 strain, were grown under controlled conditions for 35 days (DAS) and subjected to 150 mM of NaCl stress. Plant biomass (PDW) and nitrogen fixation rate (NFR) were markedly affected by salt stress conditions; the highest reductions of PDW (50%) and NFR (40%) were registered at 84 DAS and 56 DAS, respectively. In addition, salinity affected the chlorophyll fluorescence parameters, decreased initial chlorophyll fluorescence (F₀) and increased the optimum quantum yield of PSII (F_v/F_m ratio). The enzyme activities sucrose synthase activity and phosphoenolpyruvate carboxylase, responsible for the carbon supply to the bacteroids by the formation of dicarboxylates, were drastically inhibited by salinity, mainly at 56 DAS with the beginning of flowering. The content of total soluble sugars and proline increased under salt stress, and these concentrations were higher in nodule than in leaf. This last result suggests that the nodule is an organ specially protected in order to maintain its functioning, even under stress conditions. Besides, the content of myoinositol and pinitol in leaves and nodules changed with the plant growth stage and the saline treatment. Under salinity stress, the concentrations of pinitol in nodule were higher than in leaf, which supports the central function of this molecule in the adaptive response of nodules to salt stress. The increase of pinitol synthesis in nodule of M. sativa under salt stress could be one of the adaptive features used by the plant.     

Global metabolic response of Escherichia coli to gnd or zwf gene-knockout, based on 13C-labeling experiments and the measurement of enzyme activities

An integrated study on cell growth, enzyme activities and carbon flux redistribution was made to investigate how the central metabolism of Escherichia coli changes with the knockout of genes in the oxidative pentose phosphate pathway (PPP). Mutants deficient in glucose-6-phosphate dehydrogenase and 6-phosphogluconate dehydrogenase were constructed by disrupting the zwf and gnd genes and were grown in minimal media with two different carbon sources, such as glucose or pyruvate. It was shown that the knockout of either gnd or zwf gene did not affect the cell growth rate significantly, but the cellular metabolism was changed. While the specific substrate uptake rate and the specific carbon dioxide evolution rate for either mutant grown on glucose were higher than those obtained for the parent strain, these two rates were markedly decreased in mutants grown on pyruvate. The measurement of enzyme activities implied a significant change in metabolism, when alternative pathways such as the Entner–Doudoroff pathway (EDP) and the malic enzyme pathway were activated in the gnd mutant grown on glucose. As compared with the parent strain, the activities of phosphoglucose isomerase were increased in mutants grown on glucose but decreased in mutants grown on pyruvate. The metabolic flux redistribution obtained based on ¹³C-labeling experiments further indicated that the direction of the flux through the non-oxidative PPP was reversed in response to the gene knockout. Moreover, the knockout of genes caused an increased flux through the tricarboxlic acid cycle in mutants grown on glucose but caused a decrease in the case of using pyruvate. There was also a negative correlation between the fluxes through malic enzyme and isocitrate dehydrogenase in the mutants; and a positive correlation was found between the fluxes through malic enzyme and phosphoenolpyruvate carboxylase.Electronic Supplementary Material Supplementary material is available in the online version of this article at     

A systematic investigation of <Emphasis Type="Italic">Escherichia coli</Emphasis> central carbon metabolism in response to superoxide stress

Background  

The cellular responses of bacteria to superoxide stress can be used to model adaptation to severe environmental changes. Superoxide stress promotes the excessive production of reactive oxygen species (ROS) that have detrimental effects on cell metabolic and other physiological activities. To antagonize such effects, the cell needs to regulate a range of metabolic reactions in a coordinated way, so that coherent metabolic responses are generated by the cellular metabolic reaction network as a whole. In the present study, we have used a quantitative metabolic flux analysis approach, together with measurement of gene expression and activity of key enzymes, to investigate changes in central carbon metabolism that occur in Escherichia coli in response to paraquat-induced superoxide stress. The cellular regulatory mechanisms involved in the observed global flux changes are discussed.     

Structural Control of Metabolic Flux

Organisms have to continuously adapt to changing environmental conditions or undergo developmental transitions. To meet the accompanying change in metabolic demands, the molecular mechanisms of adaptation involve concerted interactions which ultimately induce a modification of the metabolic state, which is characterized by reaction fluxes and metabolite concentrations. These state transitions are the effect of simultaneously manipulating fluxes through several reactions. While metabolic control analysis has provided a powerful framework for elucidating the principles governing this orchestrated action to understand metabolic control, its applications are restricted by the limited availability of kinetic information. Here, we introduce structural metabolic control as a framework to examine individual reactions'' potential to control metabolic functions, such as biomass production, based on structural modeling. The capability to carry out a metabolic function is determined using flux balance analysis (FBA). We examine structural metabolic control on the example of the central carbon metabolism of Escherichia coli by the recently introduced framework of functional centrality (FC). This framework is based on the Shapley value from cooperative game theory and FBA, and we demonstrate its superior ability to assign “share of control” to individual reactions with respect to metabolic functions and environmental conditions. A comparative analysis of various scenarios illustrates the usefulness of FC and its relations to other structural approaches pertaining to metabolic control. We propose a Monte Carlo algorithm to estimate FCs for large networks, based on the enumeration of elementary flux modes. We further give detailed biological interpretation of FCs for production of lactate and ATP under various respiratory conditions.     

Parameter estimation and dynamic control analysis of central carbon metabolism in <Emphasis Type="Italic">Escherichia coli</Emphasis>

The central carbon metabolic pathway is the most important among metabolic pathways in all microorganisms since it produces energy and precursors for biosynthesis. In this study, a dynamic model for central carbon metabolism in Escherichia coli (E. coli) consisting of the phosphotransferase (PTS) system, glycolysis, pentose-phosphate pathway (PPP), and storage materials was obtained by ameliorating the model proposed by Chassagnole et al. (2002). In order to improve the performance of the model, principal parameters were estimated through the experimental measurements of intracellular concentrations of metabolites under transient conditions. Through dynamic metabolic control analysis (MCA), the tendencies of the metabolic fluxes at branch points were investigated, and the key parameters and enzyme kinetics that most dominantly affected the productivity of the desired metabolites were determined.     

Metabolic responses to pyruvate kinase deletion in lysine producing Corynebacterium glutamicum

Background

Pyruvate kinase is an important element in flux control of the intermediate metabolism. It catalyzes the irreversible conversion of phosphoenolpyruvate into pyruvate and is under allosteric control. In Corynebacterium glutamicum, this enzyme was regarded as promising target for improved production of lysine, one of the major amino acids in animal nutrition. In pyruvate kinase deficient strains the required equimolar ratio of the two lysine precursors oxaloacetate and pyruvate can be achieved through concerted action of the phosphotransferase system (PTS) and phosphoenolpyruvate carboxylase (PEPC), whereby a reduced amount of carbon may be lost as CO₂ due to reduced flux into the tricarboxylic acid (TCA) cycle. In previous studies, deletion of pyruvate kinase in lysine-producing C. glutamicum, however, did not yield a clear picture and the exact metabolic consequences are not fully understood.

Results

In this work, deletion of the pyk gene, encoding pyruvate kinase, was carried out in the lysine-producing strain C. glutamicum lysC^fbr, expressing a feedback resistant aspartokinase, to investigate the cellular response to deletion of this central glycolytic enzyme. Pyk deletion was achieved by allelic replacement, verified by PCR analysis and the lack of in vitro enzyme activity. The deletion mutant showed an overall growth behavior (specific growth rate, glucose uptake rate, biomass yield) which was very similar to that of the parent strain, but differed in slightly reduced lysine formation, increased formation of the overflow metabolites dihydroxyacetone and glycerol and in metabolic fluxes around the pyruvate node. The latter involved a flux shift from pyruvate carboxylase (PC) to PEPC, by which the cell maintained anaplerotic supply of the TCA cycle. This created a metabolic by-pass from PEP to pyruvate via malic enzyme demonstrating its contribution to metabolic flexibility of C. glutamicum on glucose.

Conclusion

The metabolic flux analysis performed illustrates the high flexibility of the metabolic network of C. glutamicum to compensate for external perturbation. The organism could almost maintain its growth and production performance through a local redirection of the metabolic flux, thereby fulfilling all anabolic and catabolic needs. The formation of the undesired overflow metabolites dihydroxyacetone and glycerol, in the deletion mutant, however, indicates a limiting capacity of the metabolism down-stream of their common precursor glyceraldehyde 3-phosphate and opens possibilities for further strain engineering.     

Optimizing Metabolite Production Using Periodic Oscillations

Methods for improving microbial strains for metabolite production remain the subject of constant research. Traditionally, metabolic tuning has been mostly limited to knockouts or overexpression of pathway genes and regulators. In this paper, we establish a new method to control metabolism by inducing optimally tuned time-oscillations in the levels of selected clusters of enzymes, as an alternative strategy to increase the production of a desired metabolite. Using an established kinetic model of the central carbon metabolism of Escherichia coli, we formulate this concept as a dynamic optimization problem over an extended, but finite time horizon. Total production of a metabolite of interest (in this case, phosphoenolpyruvate, PEP) is established as the objective function and time-varying concentrations of the cellular enzymes are used as decision variables. We observe that by varying, in an optimal fashion, levels of key enzymes in time, PEP production increases significantly compared to the unoptimized system. We demonstrate that oscillations can improve metabolic output in experimentally feasible synthetic circuits.     

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.     

Induction of Crassulacean Acid Metabolism in the Facultative Halophyte Mesembryanthemum crystallinum by Abscisic Acid

The facultative halophyte, Mesembryanthemum crystallinum, shifts its mode of carbon assimilation from the C₃ pathway to Crassulacean acid metabolism (CAM) in response to water stress. In this study, exogenously applied abscisic acid (ABA), at micromolar concentrations, could partially substitute for water stress in induction of CAM in this species. ABA at concentrations of 5 to 10 micromolar, when applied to leaves or to the roots in hydroponic culture or in soil, induced the expression of CAM within days (as indicated by the nocturnal accumulation of total titratable acidity and malate). After applying ABA there was also an increase in phosphoenolpyruvate carboxylase and NADP-malic enzyme activities. The degree and time course of induction by ABA were comparable to those induced by salt and water stress. Electrophoretic analyses of leaf soluble protein indicate that the increases in phosphoenolpyruvate carboxylase activity during the induction by ABA, salt, and water stress are due to an increase in the quantity of the enzyme protein. ABA may be a factor in the stress-induced expression of CAM in M. crystallinum, serving as a functional link between stress and biochemical adaptation.     

Characterization of physiological responses to 22 gene knockouts in Escherichia coli central carbon metabolism

Understanding the impact of gene knockouts on cellular physiology, and metabolism in particular, is centrally important to quantitative systems biology and metabolic engineering. Here, we present a comprehensive physiological characterization of wild-type Escherichia coli and 22 knockouts of enzymes in the upper part of central carbon metabolism, including the PTS system, glycolysis, pentose phosphate pathway and Entner–Doudoroff pathway. Our results reveal significant metabolic changes that are affected by specific gene knockouts. Analysis of collective trends and correlations in the data using principal component analysis (PCA) provide new, and sometimes surprising, insights into E. coli physiology. Additionally, by comparing the data-to-model predictions from constraint-based approaches such as FBA, MOMA and RELATCH we demonstrate the important role of less well-understood kinetic and regulatory effects in central carbon metabolism.