Skip to main content
  • Research article
  • Open access
  • Published:

Analysis of l-glutamic acid fermentation by using a dynamic metabolic simulation model of Escherichia coli



Understanding the process of amino acid fermentation as a comprehensive system is a challenging task. Previously, we developed a literature-based dynamic simulation model, which included transcriptional regulation, transcription, translation, and enzymatic reactions related to glycolysis, the pentose phosphate pathway, the tricarboxylic acid (TCA) cycle, and the anaplerotic pathway of Escherichia coli. During simulation, cell growth was defined such as to reproduce the experimental cell growth profile of fed-batch cultivation in jar fermenters. However, to confirm the biological appropriateness of our model, sensitivity analysis and experimental validation were required.


We constructed an l-glutamic acid fermentation simulation model by removing sucAB, a gene encoding α-ketoglutarate dehydrogenase. We then performed systematic sensitivity analysis for l-glutamic acid production; the results of this process corresponded with previous experimental data regarding l-glutamic acid fermentation. Furthermore, it allowed us to predicted the possibility that accumulation of 3-phosphoglycerate in the cell would regulate the carbon flux into the TCA cycle and lead to an increase in the yield of l-glutamic acid via fermentation. We validated this hypothesis through a fermentation experiment involving a model l-glutamic acid-production strain, E. coli MG1655 ΔsucA in which the phosphoglycerate kinase gene had been amplified to cause accumulation of 3-phosphoglycerate. The observed increase in l-glutamic acid production verified the biologically meaningful predictive power of our dynamic metabolic simulation model.


In this study, dynamic simulation using a literature-based model was shown to be useful for elucidating the precise mechanisms involved in fermentation processes inside the cell. Further exhaustive sensitivity analysis will facilitate identification of novel factors involved in the metabolic regulation of amino acid fermentation.


Understanding metabolic behavior in terms of a system is important for the design and fabrication of useful substances by biotechnological approaches [1]. Due to the progress in -omics studies and bioinformatics, including computer simulation, we can now integrate knowledge from different levels such as that pertaining to gene expression, protein expression, and metabolite concentrations. This information is useful for designing a strain for producing a particular substance [2]. Dynamic metabolic modeling is a useful approach for studying the regulation of metabolism by integration of biological knowledge from different biochemical systems. For example, mathematical modeling of the tryptophan operon in Escherichia coli clearly showed that TrpR and TrpL were responsible for high and low intracellular tryptophan concentrations, respectively [3]. This knowledge formed the theoretical basis for engineering a tryptophan-production strain [4].

l-glutamic acid, a flavor enhancer that is produced worldwide in quantities of over 2 million metric tons per year, is a very important fermentation product and is typically produced by fermentation using Corynebacterium glutamicum[5, 6]. Studies on the mechanism underlying l-glutamic acid overproduction in C. glutamicum are in progress and have shown the importance of α-ketoglutarate dehydrogenase activity [7]. Additionally, E. coli has been used as a model microorganism for research and also as an industrial producer of useful substances, including amino acids and organic acids [8]. Specifically, the E. coli MG1655 ΔsucA strain, which lacks α-ketoglutarate dehydrogenase activity, has been used as a model l-glutamic acid-production strain [9, 10].

To optimize carbon flux through gene deletion or gene amplification is one of the key technologies for production of substances based on fermentation. Several algorithms and methodologies have been proposed for the identification of a target for molecular breeding in metabolic engineering. Flux balance analysis is commonly used to obtain the theoretical maximum yield and the optimal biosynthesis pathway [11, 12]. Additionally, OptKnock is used as a strategy for knocking out genes or pathways in order to find optimal biosynthesis pathways [13]. Elementary mode analysis is also used to analyze metabolic network analysis to obtain all possible combinations of reaction networks [14, 15]. These approaches are based on a stoichiometric matrix of chemical reactions and static analyses. In contrast, dynamic simulation of cell metabolism is expected to become a useful method for analyzing and elucidating not only a metabolic state but also all transient in vivo cellular bioprocesses; such dynamic modeling requires knowledge of kinetic parameters.

There are mainly 2 approaches to obtaining kinetic parameters. One is parameter estimation, used for determining the coefficients in power-law equations utilizing the S-system [16]. The other involves obtaining the Michaelis constants and catalytic constants from biochemical experiments, in which Michaelis-Menten-type equations are used for modeling [17, 18]. In a previous study, we constructed our own model for the glucose phosphotransferase system (PTS) in E. coli. It included the transcription of genes and their regulatory mechanisms, protein translation, and Michaelis-Menten-based approximations of enzymatic reactions, using parameters adopted from scientific literature. In a simulation study, this model indicated that amplification of ptsI would increase the specific glucose consumption rate; we subsequently validated this prediction experimentally [19]. We expanded this model and constructed a large-scale metabolic and regulatory model of E. coli central metabolism, which included the metabolic enzymes related to PTS, glycolysis, the pentose-phosphate pathway, the tricarboxylic acid (TCA) cycle, anaplerotic enzymes, and the glyoxylate shunt, as well as transcriptional regulation by the cyclic AMP receptor protein (CRP), making large colonies protein (Mlc), catabolite repressor/activator (Cra), pyruvate dehydrogenase complex repressor (PdhR), and acetate operon repressor (IclR) [10]. We also proposed a modeling approach to describe the successive transient phases in batch or fed-batch cultivation and successfully simulated the l-glutamic acid fermentation process using E. coli MG1655 ΔsucA as a model strain [10]. In this study, we improved our previous model, performed sensitivity analyses and simulation, and validated our model mathematically. This model led us to discover a new factor influencing l-glutamate production, which we then verified experimentally.


Sensitivity analysis for l-glutamic acid production

Using a dynamic simulation model, we calculated the difference in l-glutamic acid production yields in response to variation in gene copy number. We defined the ratio of this difference to the variation in copy number as sensitivity. Because the sensitivity value changes depending on copy number variation, we obtained a maximum sensitivity value by changing the copy number variation from 0.001 to 1000. Moreover, we defined a scaling factor in Table 1 as the copy number that provides the maximum sensitivity value. The maximum sensitivity of l-glutamic acid production to the increase or decrease in copy number of 44 individual gene species is shown in Table 1; in this table, the top 15 genes (or operons) are arranged in descending order of maximum sensitivities. Based on this result, we classified the 15 factors that affected l-glutamic acid production in our model into 3 groups. The first group consisted of factors for which the enhancements directly involved the biosynthetic pathway of l-glutamic acid, including the sugar uptake system. The genes icdA, gdhA, gltA, fba, tpiA, gapA, pgk, and pykF, corresponding to the 1st, 4th, 5th, 6th, 8th, 11th, 12th, and 14th items, respectively, in Table 1, encode the biosynthetic enzymes of l-glutamic acid. The genes ptsG and crr, corresponding to the 9th and 13th items in Table 1, respectively, encode components of the sugar-uptake system (PtsG and IIAGlc of PTS). Enhancement of these factors is estimated to result in positive changes in l-glutamic acid production. In previous studies, an increase in citrate synthetase or glutamate dehydrogenase activity has been shown to be a factor leading to an increase in l-glutamic acid production in C. glutamicum[2022]. Through sensitivity analysis, we confirmed a similar effect in E. coli. The maximum sensitivities for an increase in the copy number of gltA, which encodes citrate synthetase, and of gdhA, which encodes glutamate dehydrogenase, are 0.17565 and 0.21119, respectively (Table 1). Furthermore, in E. coli, we observed a 4.0 g/L increase of l-glutamic acid production in the gltA amplification strain and a similar 1.6 g/L increase in the icdA amplification strain in the fermentation experiment after 48 h cultivation (Table 2). These results suggested that there are similar metabolic principles for l-glutamic acid production in E. coli and in C. glutamicum.

Table 1 Maximum sensitivity for l -glutamic acid production
Table 2 L-Glutamic acid fermentation results of gltA and icdA gene amplification strain in E. coli

The second group of factors comprised those that decreased the carbon flux in the glyoxylate shunt, that is, aceBAK and iclR, which corresponded to the 2nd and 3rd items, respectively, in Table 1. The increase in iclR expression levels and the decrease in the aceBAK operon expression level led to a decrease in the carbon flux of the glyoxylate shunt. Furthermore, icdA, encoding isocitrate dehydrogenase, which caused maximum sensitivity, was also related to carbon flux distribution at the branch point of the TCA cycle and the glyoxylate shunt, indicating that the balance of the metabolic flux between the TCA cycle and the glyoxylate shunt is of critical importance for l-glutamic acid production.

For the third group of factors, including the pdhR (7th), gpmA (10th), and eno (15th) genes, it was difficult to determine the reason why decreased expression would lead to an increase in l-glutamic acid production. The reason why a 3-fold increase in the expression of pdhR resulted in maximum sensitivity for l-glutamic acid production is discussed later in the subsection “Perturbation analysis by pykF amplification.” Attenuation of the activity of phosphoglyceride mutase (GpmA) and enolase (Eno) led to a decrease in l-glutamic acid production, although both enzymes catalyzed reactions in the l-glutamate biosynthetic pathway. We assumed that a decrease in the expression of gpmA may lead to accumulation of 3-phosphoglycerate. A decrease in the expression of eno would also result in accumulation of 3-phosphoglycerate, since reactions catalyzed by GpmA and Eno are reversible. The detailed mechanism of the effect of 3-phosphoglycerate accumulation on glutamate production has been described later.

Perturbation analysis by pgk amplification

The top 5 entries in Table 1 are factors known to influence l-glutamic acid production. To identify novel factors using our simulation model, we analyzed the effects of enzyme concentration and gene copy number on l-glutamic acid production. In this analysis, we perturbed gene copy number by multiplying the copy number with factors from 0.001 to 1000 in a logarithmic manner. The results were plotted on graphs in which the horizontal axis corresponded to gene copy number and the vertical axis corresponded to l-glutamic acid and other metabolite concentrations. Through this analysis, we found that an increase in pgk copy number correlates with an increase in l-glutamic acid production (Figure 1).

Figure 1
figure 1

Simulated results of perturbation of the epd–pgk operon. The vertical axis in each graph shows the ratio of the concentration obtained for metabolites to the initial concentration of these metabolites; the horizontal axis shows the magnitude of perturbation. (A) Intermediates of the glycolysis pathway. (B) Intermediates of the pentose phosphate pathway. (C) Intermediates of the TCA cycle.

Amplification of pgk was a novel important factor for l-glutamic acid production, and we attempted to interpret the mechanism underlying this observation. In E. coli, pgk encodes phosphoglycerate kinase. Phosphoglycerate kinase is known to catalyze the reaction converting 1,3-bisphosphoglycerate and ADP to 3-phosphoglycerate and ATP [23]. An increase in the pgk copy number would cause accumulation of 3-phosphoglycerate, which seemed to be closely correlated with the relationship between the decreased expression of gpmA and eno and accumulation of 3-phosphoglycerate. It is known that 3-phosphoglycerate acts as an inhibitor of isocitrate dehydrogenase kinase/phosphatase, which is encoded by aceK in E. coli[24, 25]. Since unphosphorylated isocitrate dehydrogenase retains its activity, an increase in unphosphorylated isocitrate dehydrogenase is associated with an increase in l-glutamic acid production (Figure 2). In Figure 3, we have shown the accumulation of 3-phosphoglycerate and the consequential increase in l-glutamic acid production observed when simulating a100-fold amplification of pgk. If our model, simulation, and interpretation are appropriate, then the data indicate that more than 100 copies of pgk would increase l-glutamic acid production experimentally.

Figure 2
figure 2

Predicted mechanism underlying changes in l -glutamic acid production by pgk gene amplification. A schematic representation of the simplified metabolic pathway from GA3P to l-glutamic acid. The effect of increased 3-phosphoglycerate concentration on phosphorylation of isocitrate dehydrogenase is indicated.

Figure 3
figure 3

Simulation of 100-fold amplification of the pgk gene. Each color shows the relative change in concentration of each molecule. The circles next to the compound names indicate the relative change in concentration of the compound. The circles next to gene product names show the relative change in the concentration of the enzyme and mRNA encoded by the corresponding genes. The relative changes in mRNA and enzyme levels are presented in inner circles and outer circumferences, respectively.

l-glutamic acid production in a strain with amplification of pgk

To confirm the effect of pgk gene amplification in l-glutamic acid fermentation in E. coli, we constructed a pgk-expression plasmid by using the high copy-number plasmid pUC118. The gene pgk was cloned into this vector with its native regulatory sequence, promoter, and Shine-Dalgarno sequence [26]. The pgk expression vector was introduced into a model l-glutamate-producing strain, E. coli MG1655 ΔsucA. The results for l-glutamic acid fermentation are shown in Table 3. After 72 h of cultivation, no residual glucose was observed in all the strains. A 3.2 g/L increase in l-glutamic acid accumulation and a decrease in optical density at 600 nm by 2.0 were observed for cultures in which pUC118-pgk had been introduced, in comparison with those into which pUC118 had been introduced as a control. This result suggested that our hypothesis regarding the increase in l-glutamic acid production by pgk gene amplification was validated qualitatively.

Table 3 l -Glutamic acid fermentation results

During sensitivity analysis, an increase in gltA and gdhA expression levels was predicted to lead to an increase in l-glutamic acid production (Table 1). In a previous study, the RSF-PPG plasmid, containing the genes ppc, prpC, and gdhA was used in l-glutamic acid fermentation by Pantoea ananatis, a member of the Enterobacteriaceae family [2729]. The genes ppc and gdhA encode PEP carboxylase and glutamate dehydrogenase, respectively, and prpC encodes citrate synthase activity in addition to methylcitrate synthase activity [30]. By introducing the RSF-PPG plasmid into E. coli MG1655 ΔsucA/pUC118, a 4.6 g/L increase in l-glutamic acid accumulation was observed, in comparison with a strain into which pUC118 had been introduced as a control (Table 3). This result qualitatively supported our sensitivity analysis results. After72-h cultivation, there was almost no residual glucose in the medium. By introducing the pUC118-pgk plasmid into E. coli MG1655 ΔsucA/RSF-PPG, a 1.3 g/L increase in l-glutamic acid accumulation and a decrease of 1.0 in the optical density at 600 nm were observed, when compared with E. coli MG1655 ΔsucA/RSF-PPG (Table 3). This experimental result suggested that our hypothesis regarding pgk amplification was also valid in a higher l-glutamic acid-accumulating strain and that the increase in l-glutamic acid production was at least partially due to a decrease in biomass production.

Perturbation analysis by pykF amplification

An interesting phenomenon was observed when we perturbed the copy number of pykF, which encodes pyruvate kinase. A maximum increase in l-glutamic acid production was observed at 5-fold amplification of the pykF copy number (Table 1). Further enhancement of the pykF copy number led to a decrease in l-glutamic acid production. We propose that the mechanism underlying these findings involves the following. On increase in the pykF copy number, the conversion of pyruvate to PEP increases, resulting in a decrease in the acetyl-CoA supply from pyruvate and an increase in glucose incorporation via glucose PTS (Figure 4A). The decrease in acetyl-CoA supply was closely related to the decrease in isocitrate supply (Figure 4A). As a result of pykF amplification, which raised PEP consumption, the concentration of 3-phosphoglycerate would be decreased, which would lead to increased inactivation of isocitrate dehydrogenase by the activation of isocitrate dehydrogenase kinase (Figure 4B, C). Our simulation study suggested that, even though the amount of the activated form of isocitrate dehydrogenase was reduced, the carbon flux toward α-ketoglutarate could be slightly enhanced by reduction of the carbon flux into the glyoxylate shunt when the isocitrate concentration is low. As a result, the l-glutamic acid production would increase. When the copy number of pykF is amplified more than 5 times, it leads to a decrease in 3-phosphoglycerate concentration, while most of isocitrate dehydrogenase would be inactivated (Figure 4B). As a result, l-glutamic acid production would be decreased. We suggest that this is why moderate enhancement of pdhR expression could have the same effect on pyruvate concentration as pykF amplification.

Figure 4
figure 4

Simulation of pykF perturbation. The vertical axis in each graph shows the ratio of the concentration obtained for metabolites to the initial concentration of these metabolites or the ratio of enzymatic activity obtained to the initial enzymatic activity, while the horizontal axis shows the magnitude of perturbation. (A) Metabolites with a change rate between 0.5 and 2 fold. (B) Metabolites with a change rate between 0.1 and 10 fold; data for 3PG and PEP overlap. (C) Enzymatic activity change.


In our previous study, we constructed a dynamic simulation model of E. coli based on biological knowledge and reproduced the experimental cultivation results by parameter fitting [10]. In this study, we attempted to elucidate novel factors that affect l-glutamic acid fermentation by using dynamic simulation based on a computer-aided rational design of biochemical networks. First, we refined the model with respect to biomass production through model validation. Then, a precise sensitivity analysis was performed and revealed many factors that would be important for l-glutamic acid fermentation. For example, an increase in the expression of gltA, which encodes citrate synthase, icdA, which encodes isocitrate dehydrogenase, and a combined increase in the expression of both these genes were predicted to have a high impact on l-glutamic acid production. In fact, an increase in the expression of gltA or icdA enhanced l-glutamic acid production in E. coli (Table 2). These genes have already been utilized to optimize l-glutamic acid fermentation in an industrial strain, C. glutamicum[20, 21], thus supporting the accuracy of our dynamic simulation model for understanding E. coli metabolism. In a previous study, we proposed that the putative transcriptional regulator YdcI controls carbon flux into the TCA cycle in E. coli[31]. We observed that an increase of citrate synthase activity by deletion of ydcI led to an increase in l-glutamic acid production, and a decrease of citrate synthase activity due to ydcI amplification led to a decrease in l-glutamic acid production. Our sensitivity analysis results support a clear relationship between YdcI and l-glutamic acid production, because both sensitivity analysis and experimental results clearly showed that a change in citrate synthase expression levels exert a significant effect on l-glutamic acid production in E. coli. Integration of the transcriptional regulator YdcI into our dynamic simulation model is a subject for future studies.

In theoretical flux analysis, the metabolic flux distribution, which facilitates the maximal theoretical yield of l-glutamic acid in E. coli, indicated that flux through α-ketoglutarate dehydrogenase (encoded by sucAB) and the glyoxylate shunt pathway (encoded by aceBA) should be 0. Thus, to achieve maximal l-glutamic acid production from the MG1655 ΔsucA strain, deletion of the aceBAK operon, which encodes enzymes in the glyoxylate shunt pathway and isocitrate dehydrogenase kinase/phosphatase, or amplification of iclR, which encodes a negative regulator of the glyoxylate shunt, would be a plausible approach. In contrast, using the dynamic simulation model in addition to sensitivity analysis, we identified unexpected factors. This approach indicated that amplification of pgk and attenuation of gpmA and/or eno would induce accumulation of 3-phosphoglycerate, which inhibits the phosphorylation of isocitrate dehydrogenase, and consequently results in increased glutamate production. Furthermore, the effects of amplification of pykF or pdhR on l-glutamate production could be explained according to a different mechanism. This type of working hypothesis cannot be generated using conventional theoretical flux analysis because it requires that the modification systems that control enzyme activity be taken into account. Thus, a critical advantage in using the dynamic simulation model is to be able to take various modes of regulation into consideration simultaneously and comprehensively.

In biotechnology, both production yield and productivity are important. In general, it is not simple to maintain productivity while improving production strains because the achievable yield and productivity can vary, depending on the strains and production conditions. From an engineering point of view, improving yield or productivity from a current production strain is a practical issue. One of the key features that affect productivity is biomass formation. In this study, we improved our simulation model to describe biomass production as precisely as possible; however, it was not sufficiently accurate to predict productivity. In E. coli, the carbon flux through the TCA cycle is known to affect biomass production directly [32]. Our dynamic simulation predicted the changes in concentrations of each molecule with response to gene amplification or deletion, and based on simulation results, we could speculate whether these perturbations would affect cell growth through carbon flux into the TCA cycle. In this study, we predicted that the concentration of 3-phosphoglycerate could play an important role in controlling the carbon flux through the TCA cycle; we subsequently validated this hypothesis experimentally. We also estimated that metabolic regulation through 3-phosphoglycerate would contribute to the changes in biomass production and fermentation productivity. In future, we plan to describe these changes more precisely in our model.

Our model requires further refinement. When we perturbed the copy number of ppc, which encodes PEP carboxylase, the sensitivity was 0.02138 and the ranking of this gene was 17th. In a previous study on C. glutamicum, it had been experimentally shown that an increase in PEP carboxylase activity led to an increase in l-glutamic acid yield, with a reduction in organic acid byproducts [33]. We speculated that one of the reasons for the low sensitivity of ppc was related to the process of parameter tuning because the catalytic constant of PEP carboxylase was modified to 100 times greater than the reported value. Thus, an increase in PEP carboxylase was not sensitively related to l-glutamate concentration in our model [10]. However, there is a difference in anaplerotic pathway enzymes and their regulation in E. coli and C. glutamicum[34]. Together with these facts, we speculate that unknown mechanisms, related to PEP carboxylase or an anaplerotic pathway in E. coli, could exist and result in inadequate modeling. The citrate synthase reaction requires both oxaloacetate and acetyl-CoA as substrates. An increase in PEP carboxylase expression will intensify the carbon flux toward oxaloacetate. However, to enhance citrate synthase activity, acetyl-CoA would be required. In our simulation model, the current acetate production model may be too simplistic to describe the dynamic changes in acetyl-CoA concentration inside the cell because the regulation of acetate metabolism, including formation, excretion, uptake, and utilization of this substance, is quite complex in reality [10]. In the future, improvement in simulating the metabolism of acetate, including its excretion, is a priority for refining our model.

We recognize that our model is currently limited in terms of quantitative predictive power. According to our sensitivity analysis, amplification of pgk would increase glutamate production yield by 106%, whereas, experimentally, we found a 120% increase. One of the factors that affect the prediction seems to be the difference between the experimental conditions considered for simulation and those used for validation. In the simulation model, data for the parameters were obtained using jar fermenters; however, experimental validation involved the use of shake flasks [10]. We assume that the most important issue related to the predictive power is the modeling of cell growth. In l-glutamic acid production using E. coli MG1655 ΔsucA, the succinyl-CoA used for biomass formation should be supplied through the glyoxylate shunt. Amplification of pgk in E. coli MG1655 ΔsucA decreases the carbon flux toward the glyoxylate shunt pathway. Consequently, we observed 2 phenomena, viz., increased l-glutamic acid production and decreased biomass production. At present, this trade-off has not been accommodated in our simulation because the cell growth profile is fixed in accordance with experimental results obtained using jar fermenters [10]. To further refine our model, the cell growth profile should be allowed to be variable, based on the concentration of biomass precursor molecules. Experimentally verified biological evidence will continue to be appropriately incorporated into our model during further refinements.


In this study, we evaluated a literature-based dynamic metabolic pathway model of E. coli by computational analysis and verified it experimentally. Our kinetic metabolism model was particularly useful for analysis of feedback regulation systems, and it was stable and robust against perturbation. In future, we would need to improve cell growth modelling to improve the flexibility and quantitative prediction capability of the current model. This resource will contribute to metabolic engineering that predicts the key factors for substance production.


Model and modification of biomass production

Model construction is described in Usuda et al. [10]. The enzymatic reactions and transport processes were modelled based on Michaelis–Menten-type velocity equations. Gene expression was described by transcription and translation rate equations. The systems parameters and metabolite concentrations used as constants in the simulation were obtained largely from the literature. The initial values for metabolites, messenger RNA (mRNA), and proteins were set as described [10]. The parameters used for the enzymatic reactions, gene expression, and proteins levels were basically adopted from the literature. Promoter concentrations, rate constants of transcription, and mRNA degradation rates were estimated. The most important feature was that the RNA polymerase (RNAP) and ribosome concentrations were expressed as a function of the specific growth rate, μ, which varied during batch cultivation and was calculated from experiments.

We prepared a summation model for biomass formation, which comprised the consumption and production steps of key precursor substances required for biomass production. The key precursor substances and their quantitative composition reflected the actual biomass composition, and the biomass per gram was expressed as a summation of key precursor substances with stoichiometry coefficients [35]. This approach is well established in theoretical flux analysis [36], and we applied this methodology in our dynamic simulation model. Biomass formation was taken into account, based on the cell growth and the stoichiometric matrix, by subtracting the required amount of biomass formed from the precursor metabolite. In other words, if the stoichiometric coefficient of a key precursor of biomass composition is a positive value, the key precursor is used for biomass production, whereas if the stoichiometric coefficient of a key precursor is a negative value, the production of the key precursor is associated with biomass production. However, the choice of key precursor substances in this case would not be the same as in theoretical flux analysis [10]. In our previous model, the amounts of OAA and FUM molecules required for producing 1 g of biomass achieved negative values. This indicated that cell growth causes intracellular accumulation of OAA and FUM, which would not be in accordance with the metabolic process that we expect in E. coli and may have resulted in deviations from the experimental observations. Therefore, we revised the model so that the coefficients of OAA and FUM for producing biomass attained positive values.

We observed that lipid composition influenced the coefficients of OAA in our model. Odd-numbered fatty acids, as a part of biomass production, were set to 0 because, typically, fatty acids are synthesized as multiples of 2, and production of odd-numbered fatty acids is negligibly low. Saturated and unsaturated fatty acids with carbon numbers of 13, 15, and 17 were subjects for modification of our model. In a previous model [10], the FUM required for biomass production was itself produced not only from the TCA cycle, but also through the arginine biosynthesis pathway and the nucleotide biosynthesis pathway, as follows:

eq . A ) ARGSUC > ARG + FUM

Because energy supply through the TCA cycle is primarily important for biomass formation and the stoichiometry coefficient of FUM in the previous model had a negative value, we assumed that the FUM produced via the steps shown in eq. A and eq. B was immediately and completely oxidized through the TCA cycle. We therefore modified these equations as follows:

eq . A ) ARGSUC > ARG + 4 CO 2 6 NAD + + 6 NADH 2 FAD + + 2 FADH 2 ADP 2 Pi + 2 ATP ;
eq . B ) SAICAR > AICAR + 4 CO 2 6 NAD + + 6 NADH 2 FAD + + 2 FADH 2 ADP 2 Pi + 2 ATP

These modifications resolved the biological conundrum caused in our previous model by the negative stoichiometry coefficients of OAA and FUM in the biomass composition. For glycogen production, we modified the stoichiometry as follows:

1 glycogen = 8 G 6 P + 8 ATP


In our previous model, stoichiometric coefficients indicating the relation between key precursors for biomass production and other general metabolites were not integers. One of the reasons for this phenomenon is the use of a pseudoinverse matrix in the MATLAB script. For example, the balance equation for S7P was expressed as follows:

eq . C ) 1 S 7 P = 0.5 E 4 P 1 GAP + 0.5 R 5 P + 0.5 F 6 P + 0.5 X 5 P

In our model, S7P appeared in 2 reactions:

1 E 4 P + 1 F 6 P = 1 GAP + 1 S 7 P


1 X 5 P + 1 R 5 P = 1 GAP + 1 S 7 P

When we individually transformed these reactions, the 2 reactions were expressed as follows:

1 S 7 P = 1 GAP + 1 E 4 P + 1 F 6 P


1 S 7 P = 1 GAP + 1 R 5 P + 1 X 5 P ,

in which the coefficients were integers. However, when we combined these 2 equations into a single equation, we calculated the weighted average; thus the stoichiometry of the mass balance equation for S7P is not an integer, as shown in eq. C. In our MATLAB script, we transformed these reactions using a pseudoinverse matrix, thereby obtaining equivalent average coefficients. These non-integer coefficients are caused by metabolites that are at the branch-points of a metabolic reaction. To remove these apparent defects, we added metabolites at branch-points, like S7P, to the key substances required for biomass production. In this regard, however, these modifications were made merely for purposes of calculation and do not influence biomass production.

Calculation, sensitivity, and scale factor

We used the ode15s function in MATLAB for simulations over a calculation period of 840 minutes with 5 min intervals. The sensitivity analysis (results shown in Table 1) was performed as follows. Maximum sensitivity was defined as the maximum value for the ratio of l-glutamic acid yield to the variation in gene copy number or enzyme concentration. The initial l-glutamic acid yield was defined as Y0, and the l-glutamic acid yield obtained from the parameter-modified model was defined as Y. By changing each parameter and by multiplying with 0.001 to 1000 at an exponential rate, a maximum value of (Y – Y0)/Y0 was explored. Thus, maximum sensitivity was defined as maximum value of (Y – Y0)/Y0. The scale factor was defined as the model parameter that resulted in maximum sensitivity (Table 1). The initial model parameter was defined as X0, and the modified parameter, giving the maximum sensitivity, was defined as X; the scale factor was defined as X/X0.

Strains, plasmids, and medium

The strains and plasmids used in this study are shown in Table 4. For pgk cloning, primers with the sequences 5′-gcggatccctgtaaaagccaatgaatgtc-3′ and 5′-gcaagcttattacgccaggttttacgaa-3′ were used to amplify pgk from E. coli genomic DNA. After Bam HI and Hin dIII digestion of the PCR fragment, it was cloned into pUC118 (Takara Bio.). E. coli DH5α was used as the cloning host and E. coli MG1655 ΔsucA was used as the l-glutamic acid fermentation strain. To construct the pTWV228-P tac -T trp vector, the synthesized tac promoter region and tryptophan operon terminator region were cloned into Bam HI- and Kpn I-digested pTWV228 (Takara Bio) vector. For gltA cloning, primers with sequences 5′-cacaaggagactcccatggctgatacaaaagcaaaactc-3′ and 5′-gaactggcggctcccttaacgcttgatatcgcttttaaa-3′ were used to amplify gltA from E. coli genomic DNA. For icdA cloning, primers with sequences 5′-cacaaggagactcccatggaaagtaaagtagttgttccg-3′ and 5′-gaactggcggctcccttacatgttttcgatgatcgcgtc-3′ were used to amplify icdA from E. coli genomic DNA. By using an in-fusion cloning kit (Clontech), these fragments were cloned into a Sma I-digested pTWV228-P tac -T trp vector to construct the pTWV228-P tac -gltA-T trp and pTWV228-P tac -icdA-T trp expression vectors, respectively. Shake-flask culture was used for fermentation. The composition of the medium for l-glutamic acid fermentation was as follows: 40 g/L of glucose, 1 g/L of MgSO4·7H2O, 20 g/L of (NH4)2SO4, 1 g/L of KH2PO4, 10 mg/L of FeSO4·7H2O, 10 mg/L of MnSO4·7H2O, 2 g/L of yeast extract, and 0.6 g/flask of CaCO3. Chloramphenicol (30 mg/L), Tetracycline (25 mg/L) or ampicillin (100 mg/L) was added to the media as required to select for the corresponding markers in the bacterial chromosome or plasmid. Residual glucose and l-glutamic acid concentrations were measured by using an enzymatic electrode, using BF-5 (Oji Keisokuki).

Table 4 Strains and plasmids









6-phosphogluconate dehydrogenase


Aspartate aminotransferase


Acetyl coenzyme A


Gene encoding isocitrate lyase


Gene encoding malate synthase


Gene encoding a subunit of the E1p component of pyruvate dehydrogenase complex


Gene encoding a subunit of the E2 component of pyruvate dehydrogenase complex


Gene encoding isocitrate dehydrogenase kinase/phosphatase




Gene encoding aconitase




Gene encoding aconitase




External acetate


Adenosine 5′-diphosphate


5-aminoimidazole-4-carboxamide ribotide




Adenosine monophosphate








Gene encoding aspartate aminotransferase


Adenosine 5′-triphosphate


Gene encoding malate dehydrogenase


Cyclic adenosine monophosphate


External cyclic adenosine monophosphate


cAMP efflux protein




Carbon dioxide


Coenzyme A


cAMP phosphodiesterase


Gene encoding cAMP phosphodiesterase


Catabolite repressor/activator


cAMP receptor protein


Gene encoding cAMP receptor protein


Gene encoding Crr a subunit of the glucose phosphotransferase system


Citrate synthase


Adenylate cyclase


Gene encoding adenylate cyclase




Dihydroxyacetone phosphate


2,3-bisphosphoglycerate-dependent phosphoglycerate mutase


Erythrose 4-phosphate


Phosphotransferase system enzyme I




Gene encoding enolase


Gene encoding erythrose 4-phosphate dehydrogenase


Fructose 6-phosphate


Flavin adenine dinucleotide


Reduced flavin adenine dinucleotide


Fructose bisphosphate aldolase


Gene encoding fructose bisphosphate aldolase


Fructose 1,6-bisphosphatase


Gene encoding fructose 1,6-bisphosphatase


Fructose 1,6-bisphosphate




External formate


Fumarate reductase


Gene encoding fumarate reductase flavoprotein


Gene encoding fumarate reductase iron-sulfur protein


Gene encoding fumarate reductase membrane protein


Gene encoding fumarate reductase membrane protein


Gene encoding fructose phosphotransferase system permease FruA subunit


Gene encoding fructose phosphotransferase system permease FruB subunit




Gene encoding 1-phosphofructokinase


Gene encoding catabolite repressor activator






Gene encoding fumarase


Glucose 6-phosphate




Glyceraldehyde 3-phosphate


Glyceraldehyde 3-phosphate dehydrogenase


Gene encoding glycerate 3-phosphate dehydrogenase


Glutamate dehydrogenase


Gene encoding glutamate dehydrogenase


External glucose


Gene encoding citrate synthase




External glutamate




Gene encoding 6-phosphogluconate dehydrogenase


Gene encoding 2,3-bisphosphoglycerate-dependet phosphoglycerate mutase


Histidine protein


Isocitrate dehydrogenase


Gene encoding isocitrate dehydrogenase


Isocitrate dehydrogenase kinase/phosphatase


Phosphorylated isocitrate dehydrogenase




Isocitrate lyase


Acetate operon regulator


Gene encoding isocitrate lyase regulator


IIA subunit of glucose phosphotransferase system


IICB subunit of glucose phosphotransferase system


2,3-bisphosphoglycerate-independent phosphoglycerate mutase


α-ketoglutarate dehydrogenase




Malate dehydrogenase


Gene encoding malate dehydrogenase


Making large colony protein


Gene encoding making large colony protein


Malate synthase


Nicotinamide adenine dinucleotide


Reduced nicotinamide adenine dinucleotide


Malic enzyme NAD+−linked


Nicotinamide adenine dinucleotide phosphate


Reduced nicotinamide adenine dinucleotide phosphate


Malic enzyme NADP+−linked




Phosphoenolpyruvate carboxykinase


Gene encoding phosphoenolpyruvate carboxykinase


Pyruvate dehydrogenase


Pyruvate dehydrogenase complex regulator


Gene encoding pyruvate dehydrogenase complex regulator




Phosphoenolpyruvate carboxylase




Gene encoding 6-phosphofructokinase


Phosphoglycerate mutase


Gene encoding phosphoglycerate mutase


Phosphoglycerate kinase


Gene encoding phosphoglycerate kinase


Phosphoglycerate mutase




Gene encoding phosphoenolpyruvate carboxylase


Pentose phosphate pathway


Phosphoenolpyruvate synthetase


Gene encoding phosphoenolpyruvate synthetase


Gene encoding IICBGlc


Gene encoding HPr


Gene encoding EI


Pyruvate kinase I


Gene encoding pyruvate kinase I




Ribose 5-phosphate


Ribulose 5-phosphate


Ribulose-5-phosphate 3-epimerase


Gene encoding ribulose-5-phosphate 3-epimerase


Ribose-5-phosphate isomerase A


Gene encoding ribose-5-phosphate isomerase A


Sedoheptulose 7-phosphate




Succinyl coenzyme A synthetase


Succinate dehydrogenase


Gene encoding succinate dehydrogenase flavoprotein


Gene encoding succinate dehydrogenase iron-sulfur protein


Gene encoding succinate dehydrogenase membrane protein


Gene encoding succinate dehydrogenase membrane protein


Gene encoding malic enzyme NAD+−linked


Gene encoding a subunit of the E1 component of α-ketoglutarate dehydrogenase complex


Gene encoding the SucB subunit of the α-ketoglutarate dehydrogenase complex


Gene encoding succinyl coenzyme A synthetase β subunit


Succinyl coenzyme A


Gene encoding succinyl coenzyme A synthetase α subunit


Transaldolase B


Gene encoding transaldolase B


Gene encoding transketolase


Transketolase I


Transketolase II


Triose phosphate isomerase


Gene encoding triose phosphate isomerase


Xylulose 5-phosphate


Gene encoding 2,3-bisphosphoglycerate-independent phosphoglycerate mutase


glucose Gene encoding 6-phosphate-1-dehydrogenase.


  1. Kurata H, Masaki K, Sumida Y, Iwasaki R: CADLIVE dynamic simulator: direct link of biochemical networks to dynamic models. Genome Res. 2005, 15: 590-600. 10.1101/gr.3463705.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Toyoda T, Wada A: Omic space: coordinate-based integration and analysis of genomic phenomic interactions. Bioinformatics. 2004, 20: 1759-1765. 10.1093/bioinformatics/bth165.

    Article  PubMed  Google Scholar 

  3. Santillan M, Mackey MC: Dynamic regulation of the tryptophan operon: a modeling study and comparison with experimental data. Proc Natl Acad Sci USA. 2001, 98: 1364-1369. 10.1073/pnas.98.4.1364.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Ikeda M: Towards bacterial strains overproducing L-tryptophan and other aromatics by metabolic engineering. Appl Microbiol Biotechnol. 2006, 69: 615-626. 10.1007/s00253-005-0252-y.

    Article  PubMed  Google Scholar 

  5. Kimura E: Metabolic engineering of glutamate production. Adv Biochem Eng Biotechnol. 2003, 79: 37-57.

    PubMed  Google Scholar 

  6. Sano C: History of glutamate production. Am J Clin Nutr. 2009, 90: 728S-732S. 10.3945/ajcn.2009.27462F.

    Article  PubMed  Google Scholar 

  7. Asakura Y, Kimura E, Usuda Y, Kawahara Y, Matsui K, Osumi T, Nakamatsu T: Altered metabolic flux due to deletion of odhA causes L-glutamate overproduction in Corynebacterium glutamicum. Appl Environ Microbiol. 2007, 73: 1308-1319. 10.1128/AEM.01867-06.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Wendisch VF, Bott M, Eikmanns BJ: Metabolic engineering of Escherichia coli and Corynebacterium glutamicum for biotechnological production of organic acids and amino acids. Curr Opin Microbiol. 2006, 9: 268-274. 10.1016/j.mib.2006.03.001.

    Article  PubMed  Google Scholar 

  9. Imaizumi A, Kojima H, Matsui K: The effect of intracellular ppGpp levels on glutamate and lysine overproduction in Escherichia coli. J Biotechnol. 2006, 125: 328-337. 10.1016/j.jbiotec.2006.03.015.

    Article  PubMed  Google Scholar 

  10. 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: 17-30. 10.1016/j.jbiotec.2010.02.018.

    Article  PubMed  Google Scholar 

  11. Varma A, Boesch BW, Palsson BO: Biochemical production capabilities of Escherichia coli. Biotechnol Bioeng. 1993, 42: 59-73. 10.1002/bit.260420109.

    Article  PubMed  Google Scholar 

  12. Feist AM, Herrgård MJ, Thiele I, Reed JL, Palsson BØ: Reconstruction of biochemical networks in microorganisms. Nat Rev Microbiol. 2009, 7: 129-143.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Pharkya P, Burgard AP, Maranas CD: Exploring the overproduction of amino acids using the bilevel optimization framework OptKnock. Biotechnol Bioeng. 2003, 84: 887-899. 10.1002/bit.10857.

    Article  PubMed  Google Scholar 

  14. Stelling J, Klamt S, Bettenbrock K, Schuster S, Gilles ED: Metabolic network structure determines key aspects of functionality and regulation. Nature. 2002, 420: 190-193. 10.1038/nature01166.

    Article  PubMed  Google Scholar 

  15. Klamt S, Stelling J, Ginkel M, Gilles ED: FluxAnalyzer: exploring structure, pathways, and flux distributions in metabolic networks on interactive flux maps. Bioinformatics. 2003, 19: 261-269. 10.1093/bioinformatics/19.2.261.

    Article  PubMed  Google Scholar 

  16. Kitayama T, Kinoshita A, Sugimoto M, Nakayama Y, Tomita M: A simplified method for power-law modelling of metabolic pathways from time-course data and steady-state flux profiles. Theor Biol Med Model. 2006, 3: 24-10.1186/1742-4682-3-24.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Chassagnole C, Noisommit-Rizzi N, Schmid JW, Mauch K, Reuss M: Dynamic modeling of the central carbon metabolism of Escherichia coli. Biotechnol Bioeng. 2002, 79: 53-73. 10.1002/bit.10288.

    Article  PubMed  Google Scholar 

  18. Mogilevskaya E, Bagrova N, Plyusnina T, Gizzatkulov N, Metelkin E, Goryacheva E, Smirnov S, Kosinsky Y, Dorodnov A, Peskov K, Karelina T, Goryanin I, Demin O: Kinetic modeling as a tool to integrate multilevel dynamic experimental data. Methods Mol Biol. 2009, 563: 197-218. 10.1007/978-1-60761-175-2_11.

    Article  PubMed  Google Scholar 

  19. 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-

    Article  PubMed  PubMed Central  Google Scholar 

  20. Shiio I, Ozaki H: Regulation of nicotinamide adenine dinucleotide phosphate-specific glutamate dehydrogenase from Brevibacterium flavum, a glutamate-producing bacterium. J Biochem. 1970, 68: 633-647.

    PubMed  Google Scholar 

  21. Shiio I, Ozaki H, Ujigawa K: Regulation of citrate synthase in Brevibacterium flavum, a glutamate-producing bacterium. J Biochem. 1977, 82: 395-405.

    PubMed  Google Scholar 

  22. Eikmanns BJ, Thum-Schmitz N, Eggeling L, Lüdtke KU, Sahm H: Nucleotide sequence, expression and transcriptional analysis of the Corynebacterium glutamicum gltA gene encoding citrate synthase. Microbiology. 1994, 140: 1817-1828. 10.1099/13500872-140-8-1817.

    Article  PubMed  Google Scholar 

  23. Nellemann LJ, Holm F, Atlung T, Hansen FG: Cloning and characterization of the Escherichia coli phosphoglycerate kinase (pgk) gene. Gene. 1989, 77: 185-191. 10.1016/0378-1119(89)90373-9.

    Article  PubMed  Google Scholar 

  24. LaPorte DC, Koshland DE: A protein with kinase and phosphatase activities involved in regulation of tricarboxylic acid cycle. Nature. 1982, 300: 458-460. 10.1038/300458a0.

    Article  PubMed  Google Scholar 

  25. Miller SP, Chen R, Karschnia EJ, Romfo C, Dean A, LaPorte DC: Locations of the regulatory sites for isocitrate dehydrogenase kinase/phosphatase. J Biol Chem. 2000, 275: 833-839. 10.1074/jbc.275.2.833.

    Article  PubMed  Google Scholar 

  26. Bardey V, Vallet C, Robas N, Charpentier B, Thouvenot B, Mougin A, Hajnsdorf E, Régnier P, Springer M, Branlant C: Characterization of the molecular mechanisms involved in the differential production of erythrose-4-phosphate dehydrogenase, 3-phosphoglycerate kinase and class II fructose-1,6-bisphosphate aldolase in Escherichia coli. Mol Microbiol. 2005, 57: 1265-1287. 10.1111/j.1365-2958.2005.04762.x.

    Article  PubMed  Google Scholar 

  27. Katashkina JI, Hara Y, Golubeva LI, Andreeva IG, Kuvaeva TM, Mashko SV: Use of the lambda Red-recombineering method for genetic engineering of Pantoea ananatis. BMC Mol Biol. 2009, 10: 34-10.1186/1471-2199-10-34.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Takikawa R, Hara Y: L-amino acid producing microorganism and method for producing an L-amino acid. US patent #7919284

  29. Hara Y, Kadotani N, Izui H, Katashkina JI, Kuvaeva TM, Andreeva IG, Golubeva LI, Malko DB, Makeev VJ, Mashko SV, Kozlov YI: The complete genome sequence of Pantoea ananatis AJ13355, an organism with great biotechnological potential. Appl Microbiol Biotechnol. 2012, 93: 331-341. 10.1007/s00253-011-3713-5.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Gerike U, Hough DW, Russell NJ, Dyall-Smith ML, Danson MJ: Citrate synthase and 2-methylcitrate synthase: structural, functional and evolutionary relationships. Microbiology. 1998, 144: 929-935. 10.1099/00221287-144-4-929.

    Article  PubMed  Google Scholar 

  31. Nishio Y, Suzuki T, Matsui K, Usuda Y: Metabolic control of the TCA cycle by the YdcI transcriptional regulator in Escherichia coli. J Microb Biochem Technol. 2013, 5: 59-67.

    Google Scholar 

  32. Waegeman H, Beauprez J, Moens H, Maertens J, De Mey M, Foulquié-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-10.1186/1471-2180-11-70.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Sato H, Orishimo K, Shirai T, Hirasawa T, Nagahisa K, Shimizu H, Wachi M: Distinct roles of two anaplerotic pathways in glutamate production induced by biotin limitation in Corynebacterium glutamicum. J. Biosci Bioeng. 2008, 106: 51-58. 10.1263/jbb.106.51.

    Article  PubMed  Google Scholar 

  34. Sauer U, Eikmanns BJ: The PEP-pyruvate-oxaloacetate node as the switch point for carbon flux distribution in bacteria. FEMS Microbiol Rev. 2005, 29: 765-794. 10.1016/j.femsre.2004.11.002.

    Article  PubMed  Google Scholar 

  35. Neidhardt FC, Curtiss R: Escherichia coli and Salmonella cellular and molecular biology. 1999, Washington D.C.: ASM Press

    Google Scholar 

  36. Van Dien SJ, Iwatani S, Usuda Y, Matsui K: Theoretical analysis of amino acid-producing Escherichia coli using a stoichiometric model and multivariate linear regression. J Biosci Bioeng. 2006, 102: 34-40. 10.1263/jbb.102.34.

    Article  PubMed  Google Scholar 

  37. Neidhardt FC, Ingraham JL, Schaechter M: Physiology of the bacterial cell: a molecular approach. 1990, Sunderland, Mass: Sinauer Associates

    Google Scholar 

Download references


We thank Yoshihiko Hara and Rie Takikawa for providing experimental materials and thank Yoshinori Tajima for useful discussion. We also thank Yukiko Iwata and Takahiro Sonoda for providing excellent technical support. This study was funded by Ajinomoto Co. Inc., and no external funds were used.

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Yousuke Nishio or Soichi Ogishima.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YN carried out the fermentation studies, and SO, MI, TM performed the simulation study. YY, YU, HT participated in the design of the study and helped to draft manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Nishio, Y., Ogishima, S., Ichikawa, M. et al. Analysis of l-glutamic acid fermentation by using a dynamic metabolic simulation model of Escherichia coli. BMC Syst Biol 7, 92 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: