Modeling and simulation of the redox regulation of the metabolism in Escherichia coli at different oxygen concentrations
© The Author(s) 2017
Received: 13 March 2017
Accepted: 5 July 2017
Published: 14 July 2017
Microbial production of biofuels and biochemicals from renewable feedstocks has received considerable recent attention from environmental protection and energy production perspectives. Many biofuels and biochemicals are produced by fermentation under oxygen-limited conditions following initiation of aerobic cultivation to enhance the cell growth rate. Thus, it is of significant interest to investigate the effect of dissolved oxygen concentration on redox regulation in Escherichia coli, a particularly popular cellular factory due to its high growth rate and well-characterized physiology. For this, the systems biology approach such as modeling is powerful for the analysis of the metabolism and for the design of microbial cellular factories.
Here, we developed a kinetic model that describes the dynamics of fermentation by taking into account transcription factors such as ArcA/B and Fnr, respiratory chain reactions and fermentative pathways, and catabolite regulation. The hallmark of the kinetic model is its ability to predict the dynamics of metabolism at different dissolved oxygen levels and facilitate the rational design of cultivation methods. The kinetic model was verified based on the experimental data for a wild-type E. coli strain. The model reasonably predicted the metabolic characteristics and molecular mechanisms of fnr and arcA gene-knockout mutants. Moreover, an aerobic–microaerobic dual-phase cultivation method for lactate production in a pfl-knockout mutant exhibited promising yield and productivity.
It is quite important to understand metabolic regulation mechanisms from both scientific and engineering points of view. In particular, redox regulation in response to oxygen limitation is critically important in the practical production of biofuel and biochemical compounds. The developed model can thus be used as a platform for designing microbial factories to produce a variety of biofuels and biochemicals.
KeywordsKinetic modeling Fermentation Dissolved oxygen limitation Redox regulation ArcA Fnr Respiratory chain NADH/NAD+ ratio Escherichia coli
Microbial production of biofuels and biochemicals from renewable feedstocks has received considerable recent attention from environmental protection and energy production perspectives. A limited number of cell factory platforms have been employed for the industrial production of a wide range of fuels and chemicals. Escherichia coli is probably the most widely used cellular factory due to its high growth rate and well-characterized physiology . Many biofuels and biochemicals, such as ethanol and lactate, are produced by fermentation under oxygen-limited conditions. One method in particular, dual-phase cultivation method, combines the advantages afforded by aerobic and micro-aerobic (or anaerobic) conditions [2, 3]. In dual-phase processes, cultivation is initiated with an aerobic culture to increase the biomass (contributing to productivity), and it is followed by anaerobic or micro-aerobic cultivation to facilitate efficient production of the target product. It is, therefore, highly desirable to evaluate the metabolic characteristics at different dissolved oxygen (DO) concentrations. For this purpose, appropriate quantitative models that can simulate such cultivations are needed.
Of the various modeling approaches currently available to cellular metabolism, flux balance analysis (FBA) approach has been extensively employed, but restricts to stoichiometric equations at the steady state, and thus it is difficult to simulate the dynamic changes in metabolic fluxes. On the other hand, a kinetic modeling approach can reproduce the dynamics of metabolite concentrations and fluxes in response to changes in genetic and environmental conditions , because it takes into account the mechanism of complex reactions such as allosteric modulation [5, 6], enzyme modification , and gene expression regulation by transcription factors (TFs) [8–11].
Of the various types of metabolic regulation, carbon catabolite regulation has been extensively modeled by a number of researchers [7, 12–19] to elucidate the mechanism of carbon uptake and metabolism. A detailed kinetic model of central carbon metabolism in E. coli that incorporates a constrained optimization method for parameter estimation on a supercomputer was recently developed . As compared with other kinetic models, this model enabled more accurate prediction of the dynamics of wild-type (WT) cells and multiple-gene-knockout mutants in batch culture. However, from the perspective of practical applications to develop cellular factories for biofuel and biochemical production, the effect of oxygen limitation on redox regulation with carbon catabolite regulation is critically important. Considerable effort has been expended in this regard in the Systems Understanding of Microbial Oxygen Metabolism (SUMO) project [21–27].
For the proper modeling on the respiratory chain and the redox regulation, we have to consider the basic regulation mechanisms. Oxygen serves as the final electron acceptor of the respiratory chain . In E. coli, two major oxidases, cytochrome bo (Cyo) and cytochrome bd (Cyd), transfer electrons from quinol to oxygen [29, 30]. Cyo has a low affinity for oxygen but a high reaction rate and functions primarily under aerobic conditions. By contrast, Cyd has high oxygen affinity but a lower reaction rate and functions primarily under micro-aerobic conditions. On the other hand, the dehydrogenases NADH dehydrogenase-I (Nuo) and NADH dehydrogenase-II (Ndh) oxidize electron donors such as NADH and FADH2 by reducing quinone to quinol . The function of the respiratory chain is the successive transport of electrons from electron donors to electron acceptors with translocation of protons from the cytoplasm to the periplasmic space via the inner membrane. The resulting proton gradient (proton motive force) drives ATP synthesis. This series of reactions proceeds when oxygen is available, such as under aerobic or micro-aerobic conditions.
At limited oxygen concentrations, the transcription factors Fnr and ArcA/B play essential roles in metabolic regulation in E. coli . The direct oxygen sensor Fnr regulates the expression of metabolic pathway genes under anaerobic conditions , whereas ArcA/B regulates these genes under both micro-aerobic and anaerobic conditions [32, 33]. The ArcA/B system is a two-component system: ArcB is a membrane-bound sensor kinase and ArcA is the cognate response regulator. ArcB auto-phosphorylates, and then trans-phosphorylates ArcA when oxygen is limited . Phosphorylated ArcA in turn either activates or represses the expression of metabolic pathway genes. In addition, phosphorylated ArcA represses cyoABCD, which encodes Cyo, and activates cydAB, which encodes Cyd in the respiratory chain. Note that quinone inhibits the auto-phosphorylation of ArcB , which in turn represses the activity of ArcA.
The redox ratio (i.e., NADH/NAD+) increases as the activity of the respiratory chain decreases in response to oxygen limitation. The excretion rates of fermentation products such as lactate, ethanol, succinate, formate (CO2 and H2 also), and acetate are influenced by this redox ratio. NADH is reoxidized to generate NAD+ via these fermentative pathways to enable continuation of metabolism under micro-aerobic and anaerobic conditions. Lactate is formed by lactate dehydrogenase (LDH), whereas ethanol is formed by acetaldehyde dehydrogenase (ALDH) and alcohol dehydrogenase (ADH). Succinate is formed from phosphoenol pyruvate (PEP) via phosphoenolpyruvate carboxylase (Ppc) through the reverse pathway of the normal tricarboxylic acid (TCA) cycle from oxaloacetate to succinate, whereas the succinate dehydrogenase (SDH) pathway is reversed by fumarate reductase (Frd). Formate is formed by pyruvate formate-lyase (Pfl), and acetate is formed by phosphoacetyl transferase (Pta) and acetate kinase (Ack).
In the present study, we developed a kinetic model that describes the dynamics of the metabolism in response to different DO levels by taking into account the roles of transcription factors such as ArcA/B and Fnr, the respiratory chain reactions, the fermentative pathways as mentioned above, as well as catabolite regulation.
Modeling primary metabolism
To properly model primary metabolism, the metabolic regulation mechanisms must be incorporated. Enzyme-level regulation can be represented by incorporating the effectors (metabolites) into the corresponding kinetic models. For example, in E. coli, fructose-1,6-bisphosphate (FBP) is the feed-forward activator of pyruvate kinase (Pyk) and Ppc, whereas PEP is the feedback inhibitor of phosphofructokinase (Pfk). These effectors were incorporated in the corresponding kinetic models (Additional file 1).
Model parameters were adjusted so that the model can reproduce the experimental behavior of WT strain in the batch cultures under both micro-aerobic and aerobic conditions [41, 42], whereas other parameters, including the Michaelis–Menten and dissociation constants, were retained as those given in the references (Additional file 2). MATLAB (MathWorks) was used for all simulations. The ode15s was adopted as an ordinary differential equation solver.
Experimental verification of the kinetic model
Effect of DO level on the metabolic characteristics in the WT strain
Changes in the metabolism of the WT strain with respect to DO level
As shown in Fig. 3b, the changes in the concentrations of intracellular metabolites and fluxes were simulated with respect to DO level in the WT strain. The metabolic characteristics were evaluated by classifying the DO level into four categories: (I) anaerobic condition (DO = 0%) in which both Fnr and ArcA are active; (II) micro-aerobic conditions (0% < DO < 7%) in which both Fnr and ArcA are active; (III) micro-aerobic conditions (7% ≤ DO < 20%) in which ArcA is primarily active and Fnr is inactive; and (IV) aerobic conditions (DO ≥ 20%) in which neither Fnr nor ArcA is active. TFFnr and TFArcA in Fig. 3b were calculated by Eqs. 8 and 9, respectively. The change in the typical carbon metabolism is illustrated for these categories in Additional file 3: Figure S4.
The specific oxygen uptake rate (qOUR), which indicates the rate of oxygen consumption via Cyo and Cyd reactions, was simulated to be high under condition IV, whereas it decreased under conditions III, II, and I (Fig. 3b). As the DO level decreased, the Cyd flux increased and then decreased at <7% DO. This up and down behavior can be attributed to the activation of Cyd synthesis by ArcA under condition III, whereas Cyd synthesis was repressed by Fnr under conditions I and II (Fig. 1). The Cyo flux was simulated to be higher than the Cyd flux under condition IV, whereas the Cyd flux was more dominant than the Cyo flux under condition II. These simulation results are supported by the experimental fact that the affinity of Cyd to oxygen is higher than that of Cyo . Since quinone is produced by Cyo and Cyd, quinone decreases with a decrease in DO level. This phosphorylates ArcB and then ArcA, resulting in the increase in the ArcA activity under conditions I, II, and III.
Among the enzymes associated with consumption of pyruvate, Pfl, pyruvate dehydrogenase (PDH), and LDH, play critical roles in determining the metabolite formation pattern. As the DO level decreased, the Pfl flux was simulated to increase under condition III because ArcA activated the Pfl reaction (Fig. 1). The Pfl flux was further enhanced by Fnr and ArcA under conditions I and II. In contrast, as the DO level decreased, the PDH flux decreased under condition III because ArcA represses the aceE/F genes that encode PDH (Fig. 1). The Pfl and PDH fluxes were both active under condition III, which was consistent with the experimental data . The LDH flux exhibited an up and down behavior with respect to DO level. As the DO level decreased, the LDH flux increased under condition III, whereas it declined steeply under condition II.
The NADH/NAD+ ratio increased steeply with decreasing DO level under conditions I and II because NADH is hard to be consumed by the NADH dehydrogenases in the respiratory chain. This simulation result was consistent with the experimental data . A high NADH/NAD+ ratio promoted the ADH reaction, which resulted in enhanced ethanol production under conditions I and II. Since Fnr activates the Frd flux (Fig. 1), succinate production was enhanced under conditions I and II.
To obtain a better understanding of the mechanisms by which ATP and NADH are produced or consumed under the categorized DO conditions examined, the specific production/consumption rates of ATP and NADH were simulated, as illustrated in Fig. 3c, d. DO levels of 0, 3, 8, and 40% were selected as the representatives of conditions I, II, III, and IV, respectively. The specific ATP production rate decreased in the order of conditions IV, III, II, and I (Fig. 3c). Additional file 3: Figure S5 indicates the relationship between the specific ATP production rate and the specific growth rate. Once the specific ATP production rate was calculated by Eq. 1, the specific growth rate was estimated by Eq. 3. This linear relationship between the specific ATP production rate and the specific growth rate held not only under aerobic conditions but also under micro-aerobic and anaerobic conditions with a correlation coefficient of 0.92 (p < 0.05), as experimentally observed [25, 37, 38, 41, 42, 47].
The DO level affected the specific ATP production rate (Fig. 3c). ATP was primarily synthesized by respiration under condition IV. By contrast, substrate-level phosphorylation by glycolysis and acetate formation became dominant under conditions I and II. NADH was consumed by the NADH dehydrogenases (Nuo and Ndh) in the respiratory chain under condition IV, whereas NADH was primarily consumed by ethanol formation under conditions I and II (Fig. 3d). This simulation result demonstrates that the ADH flux increased under conditions I and II due to a high NADH/NAD+ ratio (Fig. 3b). In fact, it was experimentally shown that ethanol is produced under micro-aerobic and anaerobic conditions [42, 47, 50]. The reaction of reductive TCA arm via malate dehydrogenase (MDH) consumed NADH under conditions I and II (Fig. 3d). The resultant fumarate/malate were supplied as the substrates for the reaction of Fnr-activated Frd (Fig. 3b), producing succinate (Fig. 3a). This simulation result was consistent with the experimental observation [47, 50]. NADH was produced by glycolysis, the PDH reaction, and the TCA cycle under condition IV (Fig. 3d), whereas NADH production by the PDH flux and TCA cycle declined significantly under conditions I and II because ArcA represses the PDH flux and both ArcA and Fnr repress the TCA cycle.
Additional file 3: Figure S3A shows the carbon balances of the extracellular products, CO2, and biomass at different DO levels (conditions I, II, III, and IV) in the WT strain. The metabolic modes changed significantly depending on DO level. Most of glucose was converted to biomass, CO2, and acetate under condition IV. On the other hand, biomass and CO2 production were decreased under condition I.
Prediction of the metabolic characteristics of an fnr-knockout mutant
Figure 4b shows the effect of DO level on the intracellular metabolic fluxes, redox status, and transcriptional activities. The Pfl flux was predicted to be lower in the fnr-knockout mutant than in the WT strain (Fig. 3b) under conditions I and II. The LDH flux increased under condition III and slightly decreased under conditions I and II in the fnr-knockout mutant, but it was higher than that of the WT strain (Fig. 3b). As the DO decreased, the NADH/NAD+ ratio increased under conditions III, II, and I, which resulted in the increased ADH flux, while the Frd flux was zero due to a lack of Fnr.
The simulated specific production/consumption rates of ATP and NADH are shown in Fig. 4c, d. The profiles of the specific ATP production rates of the fnr-knockout mutant were almost the same as those of the WT strain (Fig. 3c), whereas the specific NADH consumption rate in the lactate and ethanol formation through LDH and ADH and the respiratory pathway somewhat differed from that of the WT strain under condition I (Figs. 3d, 4d). The NADH consumption rate through LDH in the fnr-knockout mutant was higher than that of the WT strain (as discussed later), whereas the NADH consumption rate through ADH was lower than that of the WT strain. As compared with the WT strain, the NADH consumption rate by the NADH dehydrogenases in the respiratory chain increased under condition II because the lack of Fnr de-repressed the NADH dehydrogenase reactions.
Additional file 3: Figure S3B shows the carbon balances of the metabolic products including CO2 and biomass in the fnr-knockout mutant at different DO levels. The carbon balances differed between the WT strain and the fnr-knockout mutant under conditions I and II (Additional file 3: Figure S3A, B). More glucose carbon was converted into lactate in the fnr-knockout mutant than in the WT strain.
Prediction of the metabolic characteristics of an arcA-knockout mutant
Figure 5b shows the effect of the DO level on the intracellular metabolic fluxes, redox status, and transcriptional activities. The Cyd flux was lower than that of the WT strain under conditions III and II (Figs. 3b, 5b), as experimentally observed . As the DO level decreased under conditions III and II, the PDH flux slightly decreased. The decrease in the PDH flux was small compared to that of the WT strain and the fnr-knockout mutant (Figs. 3b, 4b, 5b) because the PDH flux is not repressed in the arcA-knockout mutant. The NADH/NAD+ ratio in the arcA-knockout mutant was higher than that of the WT strain and the fnr-knockout mutant under condition II (Figs. 3b, 4b, 5b), as experimentally observed . As DO level decreased, the NADH/NAD+ ratio increased and then declined slightly as experimentally observed . The simulated NADH/NAD+ ratio of the arcA-knockout mutant was higher than its experimental ratio, although the simulated NADH/NAD+ ratios of the WT strain and the fnr-knockout mutant were relatively consistent with their experimental ratios. While the activities of the PDH, citrate synthase (CS), and isocitrate dehydrogenase (ICDH) enzymes are allosterically inhibited by NADH to suppress an excess production of NADH [56, 57], the present model did not implement such allosteric inhibitions. The neglect of the allosteric inhibitions relatively reproduced the NADH/NAD+ ratios of the WT strain and the fnr-knockout mutant because their NADH level was not so high as that of the arcA-knockout mutant, but would overestimate the NADH/NAD+ ratio of the arcA-knockout mutant. As the DO level decreased, the Frd flux steeply increased and then declined (Fig. 5b). The simulation result of the Frd flux showed the similar trend as that of the NADH/NAD+ ratio due to the fact that NADH is oxidized at MDH with Frd.
The specific production/consumption rates of ATP and NADH were simulated as shown in Fig. 5c, d. The specific ATP production rate increased in the arcA-knockout mutant under condition III as compared with the WT strain (Figs. 3c, 5c) because Cyo is activated in the arcA-knockout mutant (Fig. 1). For this, the qOUR for the arcA-knockout mutant was higher than that of the WT strain under condition III (Figs. 3b, 5b), as experimentally observed . The specific NADH production rate in the TCA cycle was slightly higher than that of the WT strain under condition III (Figs. 3d, 5d), as experimentally observed , because the TCA cycle is not repressed in this mutant. The specific NADH consumption rate through ethanol formation (by ALDH and ADH) was higher than that of the WT strain under conditions I and II (Figs. 3d, 5d).
Additional file 3: Figure S3C shows the carbon balances of the metabolic products including CO2 and biomass in the arcA-knockout mutant at different DO levels. More glucose carbon was converted to biomass, and more CO2 was produced in the arcA-knockout mutant than the WT strain due to de-repression of the PDH and TCA cycle under condition III (Additional file 3: Figure S3A, C).
Rational design of a method for lactate production by a pfl-knockout mutant
Dual-phase cultivation was designed to enhance the target metabolite production, starting with an aerobic cultivation to promote the cell growth, followed by an anaerobic or micro-aerobic condition to facilitate the target metabolite production. The switching time when the culture condition is changed from aerobic to micro-aerobic condition is generally a key parameter for enhanced productivity. The effect of the switching time on lactate yield (g of product/g of substrate consumed) and productivity (g/l of product concentration/h of cultivation time) was simulated for the pfl-knockout mutant when 10 g/l glucose was supplied as a carbon source (Fig. 6b). The DO levels of 40 and 1% were set to the aerobic and micro-aerobic conditions, respectively. The symbols in Fig. 6b represent the productivity of lactate obtained from the experiments [44, 58]. As expected, the yield was the highest when the cells were cultured consistently under the micro-aerobic condition, although the productivity was low. The productivity was improved to 0.81 g/l/h (at a switching time of 4.5 h) by the dual-phase cultivation, as compared to 0.38 g/l/h under the micro-aerobic condition throughout the cultivation (Fig. 6b).
Advantages of the proposed kinetic model
There are several modeling approaches to simulate the fermentation characteristics. Khodayari et al.  simulated the succinate overproduction by E. coli under both aerobic and anaerobic conditions using the kinetic model-based k-OptForce method with ensemble modeling approach and parameterization based on the data obtained from multiple mutant strains. Their model was able to predict the metabolism that improves the succinate yield under aerobic condition but failed to predict it under anaerobic condition. It is essential to predict the dynamics of the cell growth and metabolite production over a broad range of DO levels and to understand the metabolic regulation mechanisms for the rational design of useful metabolite production. To meet these requirements, we have developed a kinetic model that implements redox regulation by Fnr and ArcA into central carbon metabolism [15, 20]. An advantage in the proposed model is to accurately simulate metabolisms under anaerobic, micro-aerobic, and aerobic conditions.
The model was constructed and verified using available experimental data for the WT strain [24, 29, 41–44, 47–51, 54]. The model-predicted behaviors were validated by the experimental data of the fnr-knockout mutant , the arcA-knockout mutant [32, 52, 53, 55, 60], and the pfl-knockout mutant .
To achieve efficient production of a target metabolite, the dual-phase cultivation method was investigated to improve the lactate production using the pfl-knockout mutant. This investigation revealed the importance of the optimal switching time from aerobic to micro-aerobic conditions to maximize the productivity (Fig. 6b). In addition, the trade-off between yield and productivity must be considered in practice, because the yield decreases with increased duration of the aerobic period (Fig. 6b).
Regulation mechanisms underlying the metabolic changes in response to DO level
In the simulation of the WT strain, the LDH flux exhibited up and down changes with respect to DO level (Fig. 3b). The LDH flux increased more under condition III than under condition IV. Under condition III, ArcA repressed the PDH flux while increasing the Pfl flux. Although the total flux from pyruvate to acetyl-CoA (PDH flux + Pfl flux) was almost the same between under conditions III and IV, the increased NADH/NAD+ ratio increased the LDH flux under condition III. On the other hand, the LDH flux decreased under condition II because pyruvate, the substrate of the LDH reaction, was consumed by the Fnr-enhanced Pfl reaction.
Lactate production increased in the fnr-knockout mutant under conditions I and II (Figs. 3a, 4a) as compared with the WT strain, as experimentally observed . Since the Pfl flux was reduced in the fnr-knockout mutant under conditions I and II (Figs. 3b, 4b), the total flux from pyruvate to acetyl-CoA was also reduced as compared to the WT strain, resulting in the accumulation of pyruvate. Pyruvate was converted to lactate, accompanied by NADH consumption. On the other hand, the arcA-knockout mutant also exhibited higher lactate production than the WT strain around 4% DO under condition II (Figs. 3a, 5a) because the total flux from pyruvate to acetyl-CoA was reduced as compared to the WT strain as experimentally observed .
At lower oxygen levels under conditions I and II, the arcA-knockout mutant exhibited a marked increase in ethanol production (Figs. 3a, 5a). Although the total flux from pyruvate to acetyl-CoA was almost the same as that of the WT strain under conditions I and II, the NADH/NAD+ ratio in the arcA-knockout mutant was much higher than those in the WT strain and the fnr-knockout mutant due to the high flux of PDH (Figs. 3b, 4b, 5b). This resulted in enhanced ethanol production. These simulation results were consistent with the experimental observation except for condition I (anaerobic condition) . While the simulated ethanol production flux of the arcA-knockout mutant was higher than that of the WT strain, the experimental ethanol production flux in the arcA-knockout mutant was comparable to that in the WT strain under anaerobic condition. The discrepancy in the ethanol production under anaerobic condition would be due to the overestimation of the NADH/NAD+ ratio in the arcA-knockout mutant. This overestimation results from the fact that the simulated reductive pathway flux through Ppc-MDH/Fum-Frd in the arcA-knockout mutant was lower than the experimental flux under anaerobic condition. Such underestimation of the reductive flux may be caused by the neglect of some effectors [61, 62] on the Ppc reaction responsible for the MDH and Frd fluxes. The present model includes the effect of FBP on the Ppc activity, but did not include the effects of acetyl-CoA, malate, and aspartate.
Toward virtual metabolism
Synthetic biology aims to understand the mechanisms governing the dynamic behaviors of biochemical networks in response to environmental stresses or genetic variations and facilitate the rational design or engineering of cells at the gene-regulation level. Synthetic biology approaches consist of the construction of a rigorously defined biochemical network map, development of mathematical models, experimental validation of these models, and analysis and rational design of biological systems, ultimately leading to computer-aided design of cells [63–65]. The proposed kinetic model was constructed according to this synthetic approach to provide a platform for the rational design of biofuel and biochemical production by E. coli and for further modeling efforts, including extension to amino acid, nucleotide, lipid, and polysaccharide metabolisms, as well as cell physiology. A comprehensive dynamic model, called ‘virtual E. coli,’ is expected to reproduce the complex dynamics of a series of genetic mutants under different conditions, such as consumption of multiple sugars, nitrogen, and phosphate starvation, osmotic pressure, and changes in pH. In addition, the kinetic model of the E. coli central carbon metabolism would be a feasible reference model for constructing the kinetic models of a variety of microbes, because central carbon metabolisms are relatively conserved across them.
On the other hand, another characteristic of microbes is their metabolic variety due to evolution under various growth conditions on earth. For example, yeast produces ethanol via pyruvate decarboxylase (PDC) and ADH, Clostridia employs acetone–butanol–ethanol (ABE) pathway, and Zymomonas has the Entner–Doudoroff (ED) pathway. Since the detailed metabolic pathways depend on the microbes, it is essential to take into account their differences to construct their kinetic models, while using the E. coli kinetic model as a reference model.
It is quite important to understand metabolic regulation mechanisms from both scientific and engineering points of view. In particular, redox regulation in response to oxygen limitation is critically important in the practical production of biofuel and biochemical compounds. Therefore, we developed a kinetic model with enzymatic and transcriptional regulations to predict the dynamics of metabolism at different DO levels. Transcription factor activities, metabolite concentrations, and fluxes of the WT strain and fnr- and arcA-knockout mutants were simulated to validate the model. Using this kinetic model, a rational operation strategy for the pfl-knockout mutant was designed to enhance lactate production. A dual-phase strategy was considered that involves initial cultivation under aerobic condition to enhance the cell growth rate, with subsequent cultivation under anaerobic or micro-aerobic condition to enhance the lactate production.
Primary metabolic pathway and transport system
EI enzyme I, EIIA enzyme IIA, ED pathway, Entner–Doudoroff pathway, HPr histidine-phosphorylatable protein, PP pathway, pentose phosphate pathway, PTS phosphotransferase system, TCA cycle tricarboxylic acid cycle.
ACAL acetaldehyde, AcCoA acetyl-CoA, CIT citrate, DHAP dihydroxy acetone phosphate, E4P erythrose-4-phosphate, ETH ethanol, FBP fructose-1,6-bisphosphate, FOR formate, F6P fructose-6-phosphate, FUM fumarate, G6P glucose-6-phosphate, GAP glyceraldehyde-3-phosphate, GLC glucose, GOX glyoxylate, ICI isocitrate, αKG α-ketoglutarate, LAC lactate, MAL malate, OAA oxaloacetate, PEP phosphoenol pyruvate, 6PG 6-phosphogluconate, 6PGL 6-phosphogluconolactone, PYR pyruvate, Q quinone, QH 2 quinol, R5P ribose-5-phosphate, RU5P ribulose-5-phosphate, S7P sedoheptulose-7-phosphate, SUC succinate.
Ack acetate kinase, Acs acetyl coenzyme A synthetase, ADH alcohol dehydrogenase, ALDH acetaldehyde dehydrogenase, Cya adenylate cyclase, Cyd cytochrome bd, Cyo cytochrome bo, CS citrate synthase, Eno enolase, Fba fructose-1,6-bisphosphate aldolase, Fbp fructose bisphosphatase, Frd fumarate reductase, Fum fumarase, G6PDH glucose-6-phosphate dehydrogenase, GAPDH glyceraldehyde-3-phosphate dehydrogenase, Glk glucokinase, ICDH isocitrate dehydrogenase, Icl isocitrate lyase, αKGDH α-ketoglutarate dehydrogenase, LDH lactate dehydrogenase, MDH malate dehydrogenase, Mez malic enzyme, MS malate synthase, Ndh NADH dehydrogenase-II, Nuo NADH dehydrogenase-I, Pck phosphoenolpyruvate carboxykinase, PDH pyruvate dehydrogenase, Pfk phosphofructokinase, Pfl pyruvate formate-lyase, PGDH 6-phosphogluconate dehydrogenase, Pgk phosphoglucokinase, Pgm phosphoglucomutase, Ppc phosphoenolpyruvate carboxylase, Pps phosphoenolpyruvate synthase, Pta phosphotransacetylase, Pyk pyruvate kinase, Rpe ribulose phosphate epimerase, Rpi ribose phosphate isomerase, SDH succinate dehydrogenase, Tal transaldolase, TktA transketolase I, TktB transketolase II.
YM designed the research, developed the kinetic model, performed the simulations, analyzed the data, and wrote the manuscript. HK analyzed the data and wrote the manuscript. Both authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets supporting the results of this article are included within the article and its Additional files 1, 2, and 3. The proposed kinetic model (Matlab program) is freely available at our site: http://www.cadlive.jp/cadlive_main/Softwares/KineticModel/Ecolimetabolism.html.
This work was supported by a Grant-in-Aid for Scientific Research (B) (16H02898) from the Japan Society for the Promotion of Science (JSPS) and was partially supported by the developing key technologies for discovering and manufacturing pharmaceuticals used for next-generation treatments and diagnoses, both from the Ministry of Economy, Trade and Industry, Japan (METI), and from the Japan Agency for Medical Research and Development (AMED).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Huffer S, Roche CM, Blanch HW, Clark DS. Escherichia coli for biofuel production: bridging the gap from promise to practice. Trends Biotechnol. 2012;30(10):538–45.View ArticleGoogle Scholar
- Lange J, Takors R, Blombach B. Zero-growth bioprocesses—a challenge for microbial production strains and bioprocess engineering. Eng Life Sci. 2016. doi:10.1002/elsc.201600108.Google Scholar
- Vemuri GN, Eiteman MA, Altman E. Succinate production in dual-phase Escherichia coli fermentations depends on the time of transition from aerobic to anaerobic conditions. J Ind Microbiol Biotechnol. 2002;28(6):325–32.View ArticleGoogle Scholar
- Matsuoka Y, Shimizu K. Current status and future perspectives of kinetic modeling for the cell metabolism with incorporation of the metabolic regulation mechanism. Bioresour Bioprocess. 2015;2:4.View ArticleGoogle Scholar
- Link H, Kochanowski K, Sauer U. Systematic identification of allosteric protein-metabolite interactions that control enzyme activity in vivo. Nat Biotechnol. 2013;31(4):357–61.View ArticleGoogle Scholar
- Machado D, Herrgård MJ, Rocha I. Modeling the contribution of allosteric regulation for flux control in the central carbon metabolism of E. coli. Front Bioeng Biotechnol. 2015;3:154.Google Scholar
- Kremling A, Bettenbrock K, Laube B, Jahreis K, Lengeler JW, Gilles ED. The organization of metabolic reaction networks. III. Application for diauxic growth on glucose and lactose. Metab Eng. 2001;3(4):362–79.View ArticleGoogle Scholar
- Hardiman T, Lemuth K, Keller MA, Reuss M, Siemann-Herzberg M. Topology of the global regulatory network of carbon limitation in Escherichia coli. J Biotechnol. 2007;132(4):359–74.View ArticleGoogle Scholar
- Matsuoka Y, Shimizu K. Metabolic regulation in Escherichia coli in response to culture environments via global regulators. Biotechnol J. 2011;6(11):1330–41.View ArticleGoogle Scholar
- Perrenoud A, Sauer U. Impact of global transcriptional regulation by ArcA, ArcB, Cra, Crp, Cya, Fnr, and Mlc on glucose catabolism in Escherichia coli. J Bacteriol. 2005;187(9):3171–9.View ArticleGoogle Scholar
- Shimizu K. Toward systematic metabolic engineering based on the analysis of metabolic regulation by the integration of different levels of information. Biochem Eng J. 2009;46(3):235–51.View ArticleGoogle Scholar
- Bettenbrock K, Fischer S, Kremling A, Jahreis K, Sauter T, Gilles ED. A quantitative approach to catabolite repression in Escherichia coli. J Biol Chem. 2006;281(5):2578–84.View ArticleGoogle Scholar
- Kadir TAA, Mannan AA, Kierzek AM, McFadden J, Shimizu K. Modeling and simulation of the main metabolism in Escherichia coli and its several single-gene knockout mutants with experimental verification. Microb Cell Fact. 2010;9:88.View ArticleGoogle Scholar
- Kotte O, Zaugg JB, Heinemann M. Bacterial adaptation through distributed sensing of metabolic fluxes. Mol Syst Biol. 2010;6:355.View ArticleGoogle Scholar
- Matsuoka Y, Shimizu K. Catabolite regulation analysis of Escherichia coli for acetate overflow mechanism and co-consumption of multiple sugars based on systems biology approach using computer simulation. J Biotechnol. 2013;168(2):155–73.View ArticleGoogle Scholar
- Nishio Y, Usuda Y, Matsui K, Kurata H. Computer-aided rational design of the phosphotransferase system for enhanced glucose uptake in Escherichia coli. Mol Syst Biol. 2008;4:160.View ArticleGoogle Scholar
- Usuda Y, Nishio Y, Iwatani S, Van Dien SJ, Imaizumi A, Shimbo K, Kageyama N, Iwahata D, Miyano H, Matsui K. Dynamic modeling of Escherichia coli metabolic and regulatory systems for amino-acid production. J Biotechnol. 2010;147(1):17–30.View ArticleGoogle Scholar
- Kurata H, Maeda K, Matsuoka Y. Dynamic modeling of metabolic and gene regulatory systems toward developing virtual microbes. J Chem Eng Jpn. 2014;47(1):1–9.View ArticleGoogle Scholar
- Kremling A, Geiselmann J, Ropers D, de Jong H. Understanding carbon catabolite repression in Escherichia coli using quantitative models. Trends Microbiol. 2015;23(2):99–109.View ArticleGoogle Scholar
- Jahan N, Maeda K, Matsuoka Y, Sugimoto Y, Kurata H. Development of an accurate kinetic model for the central carbon metabolism of Escherichia coli. Microb Cell Fact. 2016;15(1):112.View ArticleGoogle Scholar
- Bettenbrock K, Bai H, Ederer M, Green J, Hellingwerf KJ, Holcombe M, Kunz S, Rolfe MD, Sanguinetti G, Sawodny O, et al. Towards a systems level understanding of the oxygen response of Escherichia coli. Adv Microb Physiol. 2014;64:65–114.View ArticleGoogle Scholar
- Ederer M, Steinsiek S, Stagge S, Rolfe MD, Ter Beek A, Knies D, de Mattos MJT, Sauter T, Green J, Poole RK, et al. A mathematical model of metabolism and regulation provides a systems-level view of how Escherichia coli responds to oxygen. Front Microbiol. 2014;5:124.View ArticleGoogle Scholar
- Rolfe MD, Ocone A, Stapleton MR, Hall S, Trotter EW, Poole RK, Sanguinetti G, Green J. Systems analysis of transcription factor activities in environments with stable and dynamic oxygen concentrations. Open Biol. 2012;2(7):120091.View ArticleGoogle Scholar
- Rolfe MD, Ter Beek A, Graham AI, Trotter EW, Asif HMS, Sanguinetti G, de Mattos JT, Poole RK, Green J. Transcript profiling and inference of Escherichia coli K-12 ArcA activity across the range of physiologically relevant oxygen concentrations. J Biol Chem. 2011;286(12):10147–54.View ArticleGoogle Scholar
- Steinsiek S, Frixel S, Stagge S, Bettenbrock K. Sumo: characterization of E. coli MG1655 and frdA and sdhC mutants at various aerobiosis levels. J Biotechnol. 2011;154(1):35–45.View ArticleGoogle Scholar
- Steinsiek S, Stagge S, Bettenbrock K. Analysis of Escherichia coli mutants with a linear respiratory chain. PLoS ONE. 2014;9(1):e87307.View ArticleGoogle Scholar
- Trotter EW, Rolfe MD, Hounslow AM, Craven CJ, Williamson MP, Sanguinetti G, Poole RK, Green J. Reprogramming of Escherichia coli K-12 metabolism during the initial phase of transition from an anaerobic to a micro-aerobic environment. PLoS ONE. 2011;6(9):e25501.View ArticleGoogle Scholar
- Borisov VB, Verkhovsky MI. Oxygen as acceptor. EcoSal Plus. 2009. doi:10.1128/ecosalplus.3.2.7.Google Scholar
- Alexeeva S, Hellingwerf KJ, Teixeira de Mattos MJ. Quantitative assessment of oxygen availability: perceived aerobiosis and its effect on flux distribution in the respiratory chain of Escherichia coli. J Bacteriol. 2002;184(5):1402–6.View ArticleGoogle Scholar
- Unden G, Bongaerts J. Alternative respiratory pathways of Escherichia coli: energetics and transcriptional regulation in response to electron acceptors. Biochim Biophys Acta. 1997;1320(3):217–34.View ArticleGoogle Scholar
- Kang Y, Weber KD, Qiu Y, Kiley PJ, Blattner FR. Genome-wide expression analysis indicates that FNR of Escherichia coli K-12 regulates a large number of genes of unknown function. J Bacteriol. 2005;187(3):1135–60.View ArticleGoogle Scholar
- Alexeeva S, Hellingwerf KJ, Teixeira de Mattos MJ. Requirement of ArcA for redox regulation in Escherichia coli under microaerobic but not anaerobic or aerobic conditions. J Bacteriol. 2003;185(1):204–9.View ArticleGoogle Scholar
- Gunsalus RP. Control of electron flow in Escherichia coli: coordinated transcription of respiratory pathway genes. J Bacteriol. 1992;174(22):7069–74.View ArticleGoogle Scholar
- Kwon O, Georgellis D, Lin ECC. Phosphorelay as the sole physiological route of signal transmission by the arc two-component system of Escherichia coli. J Bacteriol. 2000;182(13):3858–62.View ArticleGoogle Scholar
- Georgellis D, Kwon O, Lin EC. Quinones as the redox signal for the arc two-component system of bacteria. Science. 2001;292(5525):2314–6.View ArticleGoogle Scholar
- Tomashek JJ, Brusilow WSA. Stoichiometry of energy coupling by proton-translocating ATPases: a history of variability. J Bioenerg Biomembr. 2000;32(5):493–500.View ArticleGoogle Scholar
- Nanchen A, Schicker A, Sauer U. Nonlinear dependency of intracellular fluxes on growth rate in miniaturized continuous cultures of Escherichia coli. Appl Environ Microbiol. 2006;72(2):1164–72.View ArticleGoogle Scholar
- Yao R, Hirose Y, Sarkar D, Nakahigashi K, Ye Q, Shimizu K. Catabolic regulation analysis of Escherichia coli and its crp, mlc, mgsA, pgi and ptsG mutants. Microb Cell Fact. 2011;10:67.View ArticleGoogle Scholar
- Henkel SG, Ter Beek A, Steinsiek S, Stagge S, Bettenbrock K, de Mattos MJT, Sauter T, Sawodny O, Ederer M. Basic regulatory principles of Escherichia coli’s electron transport chain for varying oxygen conditions. PLoS ONE. 2014;9(9):e107640.View ArticleGoogle Scholar
- Potzkei J, Kunze M, Drepper T, Gensch T, Jaeger KE, Büchs J. Real-time determination of intracellular oxygen in bacteria using a genetically encoded FRET-based biosensor. BMC Biol. 2012;10:28.View ArticleGoogle Scholar
- Toya Y, Nakahigashi K, Tomita M, Shimizu K. Metabolic regulation analysis of wild-type and arcA mutant Escherichia coli under nitrate conditions using different levels of omics data. Mol BioSyst. 2012;8(10):2593–604.View ArticleGoogle Scholar
- Zhu J, Shimizu K. Effect of a single-gene knockout on the metabolic regulation in Escherichia coli for D-lactate production under microaerobic condition. Metab Eng. 2005;7(2):104–15.View ArticleGoogle Scholar
- Toya Y, Ishii N, Nakahigashi K, Hirasawa T, Soga T, Tomita M, Shimizu K. 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. Biotechnol Prog. 2010;26(4):975–92.Google Scholar
- Zhu J, Shimizu K. The effect of pfl gene knockout on the metabolism for optically pure D-lactate production by Escherichia coli. Appl Microbiol Biotechnol. 2004;64(3):367–75.View ArticleGoogle Scholar
- Valgepea K, Adamberg K, Nahku R, Lahtvee PJ, Arike L, Vilu R. Systems biology approach reveals that overflow metabolism of acetate in Escherichia coli is triggered by carbon catabolite repression of acetyl-CoA synthetase. BMC Syst Biol. 2010;4:166.View ArticleGoogle Scholar
- Bernal V, Castano-Cerezo S, Canovas M. Acetate metabolism regulation in Escherichia coli: carbon overflow, pathogenicity, and beyond. Appl Microbiol Biotechnol. 2016;100(21):8985–9001.View ArticleGoogle Scholar
- Gonzalez JE, Long CP, Antoniewicz MR. Comprehensive analysis of glucose and xylose metabolism in Escherichia coli under aerobic and anaerobic conditions by 13C metabolic flux analysis. Metab Eng. 2016;39:9–18.View ArticleGoogle Scholar
- Alexeeva S, de Kort B, Sawers G, Hellingwerf KJ, de Mattos MJ. Effects of limited aeration and of the ArcAB system on intermediary pyruvate catabolism in Escherichia coli. J Bacteriol. 2000;182(17):4934–40.View ArticleGoogle Scholar
- de Graef MR, Alexeeva S, Snoep JL, Teixeira de Mattos MJ. The steady-state internal redox state (NADH/NAD) reflects the external redox state and is correlated with catabolic adaptation in Escherichia coli. J Bacteriol. 1999;181(8):2351–7.Google Scholar
- Chen XW, Alonso AP, Allen DK, Reed JL, Shachar-Hill Y. Synergy between 13C-metabolic flux analysis and flux balance analysis for understanding metabolic adaption to anaerobiosis in E. coli. Metab Eng. 2011;13(1):38–48.View ArticleGoogle Scholar
- Kim HJ, Hou BK, Lee SG, Kim JS, Lee DW, Lee SJ. Genome-wide analysis of redox reactions reveals metabolic engineering targets for d-lactate overproduction in Escherichia coli. Metab Eng. 2013;18:44–52.View ArticleGoogle Scholar
- Waegeman H, Beauprez J, Moens H, Maertens J, De Mey M, Foulquie-Moreno MR, Heijnen JJ, Charlier D, Soetaert W. Effect of iclR and arcA knockouts on biomass formation and metabolic fluxes in Escherichia coli K12 and its implications on understanding the metabolism of Escherichia coli BL21 (DE3). BMC Microbiol. 2011;11:70.View ArticleGoogle Scholar
- Levanon SS, San KY, Bennett GN. Effect of oxygen on the Escherichia coli ArcA and FNR regulation systems and metabolic responses. Biotechnol Bioeng. 2005;89(5):556–64.View ArticleGoogle Scholar
- Zhu J, Shalel-Levanon S, Bennett G, San KY. Effect of the global redox sensing/regulation networks on Escherichia coli and metabolic flux distribution based on C-13 labeling experiments. Metab Eng. 2006;8(6):619–27.View ArticleGoogle Scholar
- Govantes F, Orjalo AV, Gunsalus RP. Interplay between three global regulatory proteins mediates oxygen regulation of the Escherichia coli cytochrome d oxidase (cydAB) operon. Mol Microbiol. 2000;38(5):1061–73.View ArticleGoogle Scholar
- Kim Y, Ingram LO, Shanmugam KT. Dihydrolipoamide dehydrogenase mutation alters the NADH sensitivity of pyruvate dehydrogenase complex of Escherichia coli K-12. J Bacteriol. 2008;190(11):3851–8.View ArticleGoogle Scholar
- Molgat GF, Donald LJ, Duckworth HW. Chimeric allosteric citrate synthases: construction and properties of citrate synthases containing domains from two different enzymes. Arch Biochem Biophys. 1992;298(1):238–46.View ArticleGoogle Scholar
- Liu HM, Kang JH, Qi QS, Chen GJ. Production of lactate in Escherichia coli by redox regulation genetically and physiologically. Appl Biochem Biotech. 2011;164(2):162–9.View ArticleGoogle Scholar
- Khodayari A, Chowdhury A, Maranas CD. Succinate overproduction: a case study of computational strain design using a comprehensive Escherichia coli kinetic model. Front Bioeng Biotechnol. 2015;2:76.View ArticleGoogle Scholar
- Nikel PI, Zhu J, San KY, Mendez BS, Bennett GN. Metabolic flux analysis of Escherichia coli creB and arcA mutants reveals shared control of carbon catabolism under microaerobic growth conditions. J Bacteriol. 2009;191(17):5538–48.View ArticleGoogle Scholar
- Yang C, Hua Q, Baba T, Mori H, Shimizu K. Analysis of Escherichia coli anaplerotic metabolism and its regulation mechanisms from the metabolic responses to altered dilution rates and phosphoenolpyruvate carboxykinase knockout. Biotechnol Bioeng. 2003;84(2):129–44.View ArticleGoogle Scholar
- Kai Y, Matsumura H, Inoue T, Terada K, Nagara Y, Yoshinaga T, Kihara A, Tsumura K, Izui K. Three-dimensional structure of phosphoenolpyruvate carboxylase: a proposed mechanism for allosteric inhibition. Proc Natl Acad Sci USA. 1999;96(3):823–8.View ArticleGoogle Scholar
- Kurata H, Maeda K, Onaka T, Takata T. BioFNet: biological functional network database for analysis and synthesis of biological systems. Brief Bioinform. 2013;15(5):699–709.View ArticleGoogle Scholar
- Kurata H, Masaki K, Sumida Y, Iwasaki R. CADLIVE dynamic simulator: direct link of biochemical networks to dynamic models. Genome Res. 2005;15(4):590–600.View ArticleGoogle Scholar
- Kurata H, Matoba N, Shimizu N. CADLIVE for constructing a large-scale biochemical network based on a simulation-directed notation and its application to yeast cell cycle. Nucleic Acids Res. 2003;31(14):4071–84.View ArticleGoogle Scholar