13C‐metabolic flux analysis for batch culture of Escherichia coli and its pyk and pgi gene knockout mutants based on mass isotopomer distribution of intracellular metabolites

Since most bio‐production processes are conducted in a batch or fed‐batch manner, the evaluation of metabolism with respect to time is highly desirable. Toward this aim, we applied ¹³C‐metabolic flux analysis to nonstationary conditions by measuring the mass isotopomer distribution of intracellular metabolites. We performed our analysis on batch cultures of wild‐type Escherichia coli, as well as on Pyk and Pgi mutants, obtained the fluxes and metabolite concentrations as a function of time. Our results for the wild‐type indicated that the TCA cycle flux tended to increase during growth on glucose. Following glucose exhaustion, cells controlled the branch ratio between the glyoxylate pathway and the TCA cycle, depending on the availability of acetate. In the Pyk mutant, the concentrations of glycolytic intermediates changed drastically over time due to the dumping and feedback inhibition caused by PEP accumulation. Nevertheless, the flux distribution and free amino acid concentrations changed little. The growth rate and the fluxes remained constant in the Pgi mutant and the glucose‐6‐phosphate dehydrogenase reaction was the rate‐limiting step. The measured fluxes were compared with those predicted by flux balance analysis using maximization of biomass yield or ATP production. Our findings indicate that the objective function of biosynthesis became less important as time proceeds on glucose in the wild‐type, while it remained highly important in the Pyk mutant. Furthermore, ATP production was the primary objective function in the Pgi mutant. This study demonstrates how cells adjust their metabolism in response to environmental changes and/or genetic perturbations in the batch cultivation. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Genome-scale metabolic flux analysis of Streptomyces lividans growing on a complex medium

Constraint-based metabolic modeling comprises various excellent tools to assess experimentally observed phenotypic behavior of micro-organisms in terms of intracellular metabolic fluxes. In combination with genome-scale metabolic networks, micro-organisms can be investigated in much more detail and under more complex environmental conditions. Although complex media are ubiquitously applied in industrial fermentations and are often a prerequisite for high protein secretion yields, such multi-component conditions are seldom investigated using genome-scale flux analysis. In this paper, a systematic and integrative approach is presented to determine metabolic fluxes in Streptomyces lividans TK24 grown on a nutritious and complex medium. Genome-scale flux balance analysis and randomized sampling of the solution space are combined to extract maximum information from exometabolome profiles. It is shown that biomass maximization cannot predict the observed metabolite production pattern as such. Although this cellular objective commonly applies to batch fermentation data, both input and output constraints are required to reproduce the measured biomass production rate. Rich media hence not necessarily lead to maximum biomass growth. To eventually identify a unique intracellular flux vector, a hierarchical optimization of cellular objectives is adopted. Out of various tested secondary objectives, maximization of the ATP yield per flux unit returns the closest agreement with the maximum frequency in flux histograms. This unique flux estimation is hence considered as a reasonable approximation for the biological fluxes. Flux maps for different growth phases show no active oxidative part of the pentose phosphate pathway, but NADPH generation in the TCA cycle and NADPH transdehydrogenase activity are most important in fulfilling the NADPH balance. Amino acids contribute to biomass growth by augmenting the pool of available amino acids and by boosting the TCA cycle, particularly when using glutamate and aspartate. Depletion of glutamate and aspartate causes a distinct shift in fluxes of the central carbon and nitrogen metabolism. In the current work, hurdles encountered in flux analysis at a genome-scale level are addressed using hierarchical flux balance analysis and uniform sampling of the constrained solution space. This general framework can now be adopted in further studies of S. lividans, e.g., as a host for heterologous protein production.     

Parametric analysis of metabolic fluxes of α-amylase and protease-producing Bacillus subtilis

Parametric analysis was applied for a metabolic flux model for the fed-batch culture of Bacillus subtilis producing recombinant α-amylase and protease. The metabolic flux model was formulated as a linear programming problem consisting of 49 reactions (decision variables) and 50 metabolites (equality constraints). This study was aimed to determine the response of the metabolic fluxes and objective function value of minimizing the difference between ATP consumption and ATP production (ATP balance). With regard to intracellular metabolite accumulation, the objective function value was least sensitive to variation in succinate and most sensitive to variation in malate. Amongst the variations in the accumulation rates of extracellular metabolites, the objective function value was least sensitive to variation in glutamate and most sensitive to variation in starch hydrolysis and triglyceride synthesis. A 10% variation in metabolite accumulation rates caused a maximum of 13.8% variation (standard error = 3.8%) in the objective function value.     

Quantitative metabolic flux analysis revealed uneconomical utilization of ATP and NADPH in Acremonium chrysogenum fed with soybean oil

A metabolic network was constructed for the Acremonium chrysogenum cultivation fed with soybean oil. Metabolic flux analysis indicated that the shift from exponential growth to rapid cephalosporin C (CPC) formation was accompanied by 1.63- and 5-fold carbon flux enlargement in TCA cycle and glyoxylate by-pass, respectively. The flux via pentose phosphate pathway branch was little affected during the rapid CPC formation period; the contributory explanation was that 35.6% of NADPH was consumed in the dissimilation of fatty acids. Estimation of NADPH, ATP generation, and consumption demonstrated that, with soybean oil as carbon source in rapid CPC formation phase, the NADPH consumed in fatty acid catabolism was fourfold greater than that used in the CPC biosynthesis-relevant part; simultaneously, more than 90% energy spent was not directly related to the CPC formation. Therefore, the improvement of CPC production yield through optimization of the NADPH, ATP generation, and consumption was put forward.     

Optimisation of medium composition for clavulanic acid production by Streptomyces clavuligerus

Among four different commercially available nitrogen sources containing soybean derivatives, a protein extract of soybean gave the highest yield for clavulanic acid production by Streptomyces clavuligerus. A statistical method based on factorial design of experiments was applied to optimise the medium. An empirical model was obtained by applying response surface statistical analysis. The analysis of variance showed that concentrations of protein extract of soybean and glycerol and the interaction between these two variables were significant at 95% level of confidence. The maximum clavulanic acid concentration obtained in 72 h was 1.2 g l^–1.     

Stoichiometric model and flux balance analysis for a mixed culture of Leptospirillum ferriphilum and Ferroplasma acidiphilum

The oxidation process of sulfide minerals in natural environments is achieved by microbial communities from the Archaea and Bacteria domains. A metabolic reconstruction of two dominant species, Leptospirillum ferriphilum and Ferroplasma acidiphilum, which are always found together as a mixed culture in this natural environments, was made. The metabolic model, composed of 152 internal reactions and 29 transport reactions, describes the main interactions between these species, assuming that both use ferrous iron as energy source, and F. acidiphilum takes advantage of the organic compounds secreted by L. ferriphilum for chemomixotrophic growth. A first metabolic model for a mixed culture used in bacterial leaching is proposed in this article, which pretends to represent the characteristics of the mixed culture in a simplified manner. It was evaluated with experimental data through flux balance analysis (FBA) using as objective function the maximization of biomass. The growth yields on ferrous iron obtained for each microorganism are consistent with experimental data, and the flux distribution obtained allows understanding of the metabolic capabilities of both microorganisms growing together in a bioleaching process. The model was used to simulate the growth of F. acidiphilum on different substrates, to determine in silico which compounds maximize cell growth, and which are essential. Knockout simulations were carried out for L. ferriphilum and F. acidiphilum metabolic models, predicting key enzymes of central metabolism. The results of this analysis are consistent with experimental data from literature, showing a robust behavior of the metabolic model. © 2014 American Institute of Chemical Engineers Biotechnol. Prog., 31:307–315, 2015     

Quantification of proteomic and metabolic burdens predicts growth retardation and overflow metabolism in recombinant Escherichia coli

Escherichia coli has been the host organism most frequently investigated for efficient recombinant protein production. However, the production of a foreign protein in recombinant E. coli often leads to growth deterioration and elevated secretion of acetic acid. Such observed phenomena have been widely linked with cell stress responses and metabolic burdens originated particularly from the increased energy demand. In this study, flux balance analysis and dynamic flux balance analysis were applied to investigate the observed growth physiology of recombinant E. coli, incorporating the proteome allocation theory and an adjustable maintenance energy level (ATPM) to capture the proteomic and energetic burdens introduced by recombinant protein synthesis. Model predictions of biomass growth, substrate consumption, acetate excretion, and protein production with two different strains were in good agreement with the experimental data, indicating that the constraint on the available proteomic resource and the change in ATPM might be important contributors governing the growth physiology of recombinant strains. The modeling framework developed in this work, currently with several limitations to overcome, offers a starting point for the development of a practical, model-based tool to guide metabolic engineering decisions for boosting recombinant protein production.     

调控大肠杆菌胞内ATP和NADH水平促进琥珀酸生产

琥珀酸作为一种重要的C4平台化合物,广泛应用于食品、化学、医药等领域。利用大肠杆菌(Escherichia coli)发酵生产琥珀酸受胞内辅因子不平衡的影响,存在产率低、生产强度低、副产物多等问题。为此,对不同氧气条件下琥珀酸产量和化学计量学分析发现,微厌氧条件下E.coli FMME-N-26高效积累琥珀酸需要借助三羧酸循环(tricarboxylic acid cycle,TCA)为还原性三羧酸途径(reductive tricarboxylic acid pathway,r-TCA)提供足够的ATP和NADH。通过减少ATP消耗、强化ATP合成、阻断NADH竞争途径和构建NADH回补路径等代谢工程策略,组合调控胞内ATP与NADH含量,获得工程菌株E.coli FW-17。通过发酵条件优化,菌株E.coli FW-17在5 L发酵罐能积累139.52 g/L琥珀酸,比出发菌株提高了17.81%,乙酸浓度为1.40 g/L,降低了67.59%。进一步在1000 L发酵罐中进行放大实验,琥珀酸产量和乙酸浓度分别为140.2 g/L和1.38 g/L。     

Insights into pH‐induced metabolic switch by flux balance analysis

Lactate accumulation in mammalian cell culture is known to impede cellular growth and productivity. The control of lactate formation and consumption in a hybridoma cell line was achieved by pH alteration during the early exponential growth phase. In particular, lactate consumption was induced even at high glucose concentrations at pH 6.8, whereas highly increased production of lactate was obtained at pH 7.8. Consequently, constraint‐based metabolic flux analysis was used to examine pH‐induced metabolic states in the same growth state. We demonstrated that lactate influx at pH 6.8 led cells to maintain high fluxes in the TCA cycle and malate‐aspartate shuttle resulting in a high ATP production rate. In contrast, under increased pH conditions, less ATP was generated and different ATP sources were utilized. Gene expression analysis led to the conclusion that lactate formation at high pH was enabled by gluconeogenic pathways in addition to facilitated glucose uptake. The obtained results provide new insights into the influence of pH on cellular metabolism, and are of importance when considering pH heterogeneities typically present in large scale industrial bioreactors. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:347–357, 2015     

Metabolic engineering of Escherichia coli for the production of fumaric acid

Fumaric acid is a naturally occurring organic acid that is an intermediate of the tricarboxylic acid cycle. Fungal species belonging to Rhizopus have traditionally been employed for the production of fumaric acid. In this study, Escherichia coli was metabolically engineered for the production of fumaric acid under aerobic condition. For the aerobic production of fumaric acid, the iclR gene was deleted to redirect the carbon flux through the glyoxylate shunt. In addition, the fumA, fumB, and fumC genes were also deleted to enhance fumaric acid formation. The resulting strain was able to produce 1.45 g/L of fumaric acid from 15 g/L of glucose in flask culture. Based on in silico flux response analysis, this base strain was further engineered by plasmid‐based overexpression of the native ppc gene, encoding phosphoenolpyruvate carboxylase (PPC), from the strong tac promoter, which resulted in the production of 4.09 g/L of fumaric acid. Additionally, the arcA and ptsG genes were deleted to reinforce the oxidative TCA cycle flux, and the aspA gene was deleted to block the conversion of fumaric acid into L ‐aspartic acid. Since it is desirable to avoid the use of inducer, the lacI gene was also deleted. To increase glucose uptake rate and fumaric acid productivity, the native promoter of the galP gene was replaced with the strong trc promoter. Fed‐batch culture of the final strain CWF812 allowed production of 28.2 g/L fumaric acid in 63 h with the overall yield and productivity of 0.389 g fumaric acid/g glucose and 0.448 g/L/h, respectively. This study demonstrates the possibility for the efficient production of fumaric acid by metabolically engineered E. coli. Biotechnol. Bioeng. 2013; 110: 2025–2034. © 2013 Wiley Periodicals, Inc.     

Poly-3-hydroxybutyrate production by washed cells of Alcaligenes eutrophus; purification, characterisation and potential regulatory role of citrate synthase

Washed cells prepared from carbon-limited continuous cultures of Alcaligenes eutrophus synthesised poly-3-hydroxybutyrate (PHB) rapidly when supplied with glucose, dl-lactate or l-lactate. Unlike growing cultures, washed cells excreted significant amounts of pyruvate. The combined rates of PHB production (qPHB) and pyruvate excretion (qPyr) were linearly related to the rate of carbon substrate utilisation (qS), showing that washed cells behaved similarly to growing cultures when corrected for the absence of non-PHB biomass production. The addition of formate (as a potential source of NADH and/or ATP) significantly stimulated both qPHB and qPyr, but slightly decreased qS and substantially decreased the flux of carbon through the tricarboxylic acid cycle (qTCA). Citrate synthase activity of broken cells was inhibited by physiological concentrations of NADH, but not of ATP, in a manner that was not reversible by AMP. Citrate synthase was purified and shown to be a “large” form of the enzyme (M _r 227,000), comprising a single type of subunit (M _r 47,000) as found in several other gram-negative aerobes. The potential role of citrate synthase in the regulation of PHB production via its ability to control carbon flux into the tricarboxylic acid cycle is discussed. Received: 14 March 1997 / Accepted: 9 July 1997     

Predictive Potential of Flux Balance Analysis of Saccharomyces cerevisiae Using as Optimization Function Combinations of Cell Compartmental Objectives

Background

The main objective of flux balance analysis (FBA) is to obtain quantitative predictions of metabolic fluxes of an organism, and it is necessary to use an appropriate objective function to guarantee a good estimation of those fluxes.

Methodology

In this study, the predictive performance of FBA was evaluated, using objective functions arising from the linear combination of different cellular objectives. This approach is most suitable for eukaryotic cells, owing to their multiplicity of cellular compartments. For this reason, Saccharomyces cerevisiae was used as model organism, and its metabolic network was represented using the genome-scale metabolic model iMM904. As the objective was to evaluate the predictive performance from the FBA using the kind of objective function previously described, substrate uptake and oxygen consumption were the only input data used for the FBA. Experimental information about microbial growth and exchange of metabolites with the environment was used to assess the quality of the predictions.

Conclusions

The quality of the predictions obtained with the FBA depends greatly on the knowledge of the oxygen uptake rate. For the most of studied classifications, the best predictions were obtained with “maximization of growth”, and with some combinations that include this objective. However, in the case of exponential growth with unknown oxygen exchange flux, the objective function “maximization of growth, plus minimization of NADH production in cytosol, plus minimization of NAD(P)H consumption in mitochondrion” gave much more accurate estimations of fluxes than the obtained with any other objective function explored in this study.     

Flux control-based design of furfural-resistance strains of Saccharomyces cerevisiae for lignocellulosic biorefinery

We have previously developed a dynamic flux balance analysis of Saccharomyces cerevisiae for elucidation of genome-wide flux response to furfural perturbation (Unrean and Franzen, Biotechnol J 10(8):1248–1258, 2015). Herein, the dynamic flux distributions were analyzed by flux control analysis to identify target overexpressed genes for improved yeast robustness against furfural. The flux control coefficient (FCC) identified overexpressing isocitrate dehydrogenase (IDH1), a rate-controlling flux for ethanol fermentation, and dicarboxylate carrier (DIC1), a limiting flux for cell growth, as keys of furfural-resistance phenotype. Consistent with the model prediction, strain characterization showed 1.2- and 2.0-fold improvement in ethanol synthesis and furfural detoxification rates, respectively, by IDH1 overexpressed mutant compared to the control. DIC1 overexpressed mutant grew at 1.3-fold faster and reduced furfural at 1.4-fold faster than the control under the furfural challenge. This study hence demonstrated the FCC-based approach as an effective tool for guiding the design of robust yeast strains.

     

Optimization of glycerol feeding for clavulanic acid production by Streptomyces clavuligerus with glycerol feeding

Glycerol at 10–20 g l^–1 increased clavulanic acid production by Streptomyces clavuligerus in shake-flask culture. The biosynthesis of clavulanic acid continued for longer by feeding glycerol and production increased to 250 mg l^–1 compared with 115 mg l^–1 without feeding. In fermenter batch culture, degradation of clavulanic acid began after 72 h. With glycerol feeding in fed-batch culture, clavulanic acid production was not only increased further to about 280 mg l^–1 but also remained stable up to 130 h.     

Analysis of metabolic pathways and fluxes in a newly discovered thermophilic and ethanol‐tolerant Geobacillus strain

A recently discovered thermophilic bacterium, Geobacillus thermoglucosidasius M10EXG, ferments a range of C5 (e.g., xylose) and C6 sugars (e.g., glucose) and is tolerant to high ethanol concentrations (10%, v/v). We have investigated the central metabolism of this bacterium using both in vitro enzyme assays and ¹³C‐based flux analysis to provide insights into the physiological properties of this extremophile and explore its metabolism for bio‐ethanol or other bioprocess applications. Our findings show that glucose metabolism in G. thermoglucosidasius M10EXG proceeds via glycolysis, the pentose phosphate pathway, and the TCA cycle; the Entner–Doudoroff pathway and transhydrogenase activity were not detected. Anaplerotic reactions (including the glyoxylate shunt, pyruvate carboxylase, and phosphoenolpyruvate carboxykinase) were active, but fluxes through those pathways could not be accurately determined using amino acid labeling. When growth conditions were switched from aerobic to micro‐aerobic conditions, fluxes (based on a normalized glucose uptake rate of 100 units (g DCW)^?1 h^?1) through the TCA cycle and oxidative pentose phosphate pathway were reduced from 64 ± 3 to 25 ± 2 and from 30 ± 2 to 19 ± 2, respectively. The carbon flux under micro‐aerobic growth was directed to ethanol, L ‐lactate (>99% optical purity), acetate, and formate. Under fully anerobic conditions, G. thermoglucosidasius M10EXG used a mixed acid fermentation process and exhibited a maximum ethanol yield of 0.38 ± 0.07 mol mol^?1 glucose. In silico flux balance modeling demonstrates that lactate and acetate production from G. thermoglucosidasius M10EXG reduces the maximum ethanol yield by approximately threefold, thus indicating that both pathways should be modified to maximize ethanol production. Biotechnol. Bioeng. 2009;102: 1377–1386. © 2008 Wiley Periodicals, Inc.     

Identification of metabolic model: Citrate production from glucose by Candida lipolytica

The concept of mass balance was used to analyze the metabolic pathways of citrate production by Candida lipolytica from glucose. Specific rates of glucose consumption, citrate and isocitrate productions, carbon dioxide evolution, and cellular syntheses of protein and carbohydrate were observed in an NH₄⁺-limited chemostat culture. These data permitted one to assess the carbon flux in vivo by solving simultaneous carbon balance equations with respect to intermediary metabolite pools in the steady State. Among the three models considered here, model I (which coordinates the pyruvate carboxylation with the tricarboxylic acid cycle, but disregards the glyoxylate cycle) was considered plausible because the carbon flux calculated so far was acceptable. On the other hand, models II and III (which overlook the pyruvate carboxylation and the 2-oxoglutarate dehydrogenation, respectively) were found to be most unlikely because of the unusual flux assessed from these models.     

Fed-batch culture of Escherichia coli for L-valine production based on in silico flux response analysis

We have previously reported the development of a 100% genetically defined engineered Escherichia coli strain capable of producing L ‐valine from glucose with a high yield of 0.38 g L ‐valine per gram glucose (0.58 mol L ‐valine per mol glucose) by batch culture. Here we report a systems biological strategy of employing flux response analysis in bioprocess development using L ‐valine production by fed‐batch culture as an example. Through the systems‐level analysis, the source of ATP was found to be important for efficient L ‐valine production. There existed a trade‐off between L ‐valine production and biomass formation, which was optimized for the most efficient L ‐valine production. Furthermore, acetic acid feeding strategy was optimized based on flux response analysis. The final fed‐batch cultivation strategy allowed production of 32.3 g/L L ‐valine, the highest concentration reported for E. coli. This approach of employing systems‐level analysis of metabolic fluxes in developing fed‐batch cultivation strategy would also be applicable in developing strategies for the efficient production of other bioproducts. Biotechnol. Bioeng. 2011; 108:934–946. © 2010 Wiley Periodicals, Inc.     

ATP limitation in a pyruvate formate lyase mutant of <Emphasis Type="Italic">Escherichia coli</Emphasis> MG1655 increases glycolytic flux to <Emphasis Type="SmallCaps">d</Emphasis>-lactate

A derivative strain of Escherichia coli MG1655 for d-lactate production was constructed by deleting the pflB, adhE and frdA genes; this strain was designated “CL3.” Results show that the CL3 strain grew 44% slower than its parental strain under nonaerated (fermentative) conditions due to the inactivation of the main acetyl-CoA production pathway. In contrast to E. coli B and W3110 pflB derivatives, we found that the MG1655 pflB derivative is able to grow in mineral media with glucose as the sole carbon source under fermentative conditions. The glycolytic flux was 2.8-fold higher in CL3 when compared to the wild-type strain, and lactate yield on glucose was 95%. Although a low cell mass formed under fermentative conditions with this strain (1.2 g/L), the volumetric productivity of CL3 was 1.31 g/L h. In comparison with the parental strain, CL3 has a 22% lower ATP/ADP ratio. In contrast to wild-type E. coli, the ATP yield from glucose to lactate is 2 ATP/glucose, so CL3 has to improve its glycolytic flux in order to fulfill its ATP needs in order to grow. The aceF deletion in strains MG1655 and CL3 indicates that the pyruvate dehydrogenase (PDH) complex is functional under glucose-fermentative conditions. These results suggest that the pyruvate to acetyl-CoA flux in CL3 is dependent on PDH activity and that the decrease in the ATP/ADP ratio causes an increase in the flux of glucose to lactate.