Modeling and simulation of the main metabolism in Escherichia coli and its several single-gene knockout mutants with experimental verification
© Kadir et al; licensee BioMed Central Ltd. 2010
Received: 16 July 2010
Accepted: 19 November 2010
Published: 19 November 2010
It is quite important to simulate the metabolic changes of a cell in response to the change in culture environment and/or specific gene knockouts particularly for the purpose of application in industry. If this could be done, the cell design can be made without conducting exhaustive experiments, and one can screen out the promising candidates, proceeded by experimental verification of a select few of particular interest. Although several models have so far been proposed, most of them focus on the specific metabolic pathways. It is preferred to model the whole of the main metabolic pathways in Escherichia coli, allowing for the estimation of energy generation and cell synthesis, based on intracellular fluxes and that may be used to characterize phenotypic growth.
In the present study, we considered the simulation of the main metabolic pathways such as glycolysis, TCA cycle, pentose phosphate (PP) pathway, and the anapleorotic pathways using enzymatic reaction models of E. coli. Once intracellular fluxes were computed by this model, the specific ATP production rate, the specific CO2 production rate, and the specific NADPH production rate could be estimated. The specific ATP production rate thus computed was used for the estimation of the specific growth rate. The CO2 production rate could be used to estimate cell yield, and the specific NADPH production rate could be used to determine the flux of the oxidative PP pathway. The batch and continuous cultivations were simulated where the changing patterns of extracellular and intra-cellular metabolite concentrations were compared with experimental data. Moreover, the effects of the knockout of such pathways as Ppc, Pck and Pyk on the metabolism were simulated. It was shown to be difficult for the cell to grow in Ppc mutant due to low concentration of OAA, while Pck mutant does not necessarily show this phenomenon. The slower growth rate of the Ppc mutant was properly estimated by taking into account the lower specific ATP production rate. In the case of Pyk mutant, the enzyme level regulation was made clear such that Pyk knockout caused PEP concentration to be up-regulated and activated Ppc, which caused the increase in MAL concentration and backed up reduced PYR through Mez, resulting in the phenotypic growth characteristics similar to the wild type.
It was shown to be useful to simulate the main metabolism of E. coli for understanding metabolic changes inside the cell in response to specific pathway gene knockouts, considering the whole main metabolic pathways. The comparison of the simulation result with the experimental data indicates that the present model could simulate the effect of the specific gene knockouts to the changes in the metabolisms to some extent.
One of the most challenging goals of metabolic engineering is to design the cell metabolism based on the analysis of metabolic regulation. For this, it is strongly desired to develop a mathematical model which can describe the dynamic behaviour of the cell in response to the changes in the culture environment and/or specific genetic modifications. Although an attempt has been made to develop a platform for the whole cell model , the total cell model has not yet been developed. If such a model could be developed, it becomes possible to check the metabolism of a specific gene knockout on the metabolism and fermentation characteristics without conducting many exhaustive experiments, and allow for the screening out of the preferred candidates for performance improvement, followed by experimental verification only for the selected candidates.
Some of the mathematical models which can describe the dynamic behaviour of the intracellular metabolite concentrations of the central metabolic pathways have been developed for Saccharomyces cerevisiae [2–4]. The measurement of the intracellular metabolite concentrations for the pulse addition of glucose during continuous culture has been made, and the time profile was compared to the predicted dynamic simulation together with model parameter identification [3, 5–7]. The kinetic equations for the glycolysis and the pentose phosphate (PP) pathway have also been developed for E. coli to simulate the transient data obtained by the fast sampling system . These models do not contain TCA cycle and the fermentative pathways, and thus cannot simulate the typical aerobic batch culture.
In the present research, therefore, we considered several kinetic models for the TCA cycle, anapleorotic pathways as well as the glycolysis and the PP pathway to simulate the time profiles of the batch and continuous cultures. Moreover, most of the kinetic models developed so far can express only enzyme level regulation due to the change in the concentrations of substrate and product as well as various effectors. Thus, the application of the conventional model is limited in practice to some extent. Recently, several mathematical models which describe the effects of global regulators on the metabolic pathway reactions for catabolite repression for substrate uptake  and for suc mutant for glutamate production  have been proposed, which pay attention to particular pathways. Recently, we estimated the flux changes during batch culture of E. coli based on 13C-labeling experiment using CE-TOF/MS . It is quite important to estimate the flux changes of the main metabolic pathways, allowing opportunity for the proper analysis of the energy metabolism and cell synthesis. Although 13C-metabolic flux analysis has proven to be quite useful [12, 13], it is a method of the analysis of a static physiological state of the organism and does not have predictability characteristics. It is highly important and indeed useful to be able to predict cell growth characteristics. In the present study, therefore, we attempted to develop a new model for the cell growth rate with the advantages of considering the metabolic fluxes and an enzymatic model. Furthermore, incorporating the relationship between ATP production rate obtained by the intracellular fluxes of the main metabolic pathways and the cell growth rate together with some rule-based approach for gene-level regulation. Once we could simulate the whole main metabolic pathways, we may be able to compute CO2 production rate and NAD(P)H production rate as well as ATP production rate. In particular, we attempted to simulate several single-gene knockout mutants to show the utility of the model and its limitations. Some of the experiments were also conducted to verify the simulation result, which seems to be quite important for the practical use of modeling and simulation.
Kinetic rate equations
Kinetic parameter values
μ m = 0.6, K s = 0.1, μ mA = 0.9, K sA = 0.01, X m = 2.3
KG 6p= 0.5mM
K PEP = 3.26mM
K NADH = 0.02mM, K eq = 0.63
K ATP = 22.5mM, K ADP = 0.26mM, L Pyk = 1000,
k3 = 0.0542mM, k4 = 0.8139mM, k5 = 0.0939mM, k6 = 0.2693mM
K ATPinh = 3.0mM
During the fitting procedure, it was found that some of the initial concentrations of the intracellular metabolites such as OAA, PEP and PYR affected the growth cultivation characteristics. Additional file 1a shows the effect of initial OAA concentration, while Additional file 1b shows the effect of initial PEP concentration on the fermentation characteristics. Since those will change depending on the pre-cultivation, we assumed those to be also adjustable parameters for the simulation of batch culture. Note that initial OAA concentration was low around 0.01 mmol/gDCW, while initial PEP concentration was around 0.1 mmol/gDCW. The low value of intracellular OAA concentration coincides with the fact that anapleorotic pathway is required for OAA for cell synthesis.
The changes in the intracellular metabolite concentrations for batch culture of wild type as shown in Figure 3c are given in Additional file 2, indicating that as glucose is consumed, the concentrations of metabolites G6P/F6P, FDP and GAP/DHAP increase. After glucose depletes these metabolite concentrations decreased sharply, however for PEP concentration, PTS flux became zero, and thus PEP concentration tended to increase for short a period. Unlike other glycolysis intermediates as mentioned above, the PEP concentration slowly decreased. This is due to the exchange with OAA via Ppc and Pck. The change in PYR concentration affected AcCoA concentration, where the latter increased after glucose was depleted. The concentrations of the TCA cycle intermediates changed in a similar fashion. The glyoxylate pathway was activated after the glucose concentration became low level due to rule 6 of Additional file 3: Table S1.
In the present study, we considered to simulate the main metabolism of E. coli. By considering the whole main metabolic pathways such as glycolysis, PP pathway and TCA cycle, we could estimate the specific ATP production rate, the specific NADPH production rate, and specific CO2 production rate. Those could be used to estimate the specific growth rate and the cell yield etc., and those were affected by the specific gene knockout. Thus, the simulation result is quite useful in identifying the mechanism of metabolic changes in terms of intracellular metabolite concentrations and the fluxes caused by the specific gene knockout. It should be noted that the flux balance analysis (FBA) has been known to be useful to predict the steady-state fluxes. However, since FBA is based on the mass balances with stoichiomatric matrix, it cannot treat the dynamic behavior. The effects of metabolite concentrations on the fluxes cannot be investigated by FBA, while it can be done by the kinetic equation-based approach. Moreover, although the metabolic flux analysis based on 13C-labeling experiments is powerful [12, 13], it is method of analysis and cannot predict the flux changes caused by the change in culture environment and/or specific gene knockout.
It should be noted that the modeling and simulation of the whole central metabolic pathways are quite important, since the fluxes obtained allow us to estimate how ATP production rate and CO2 production rate etc. change with respect to time and/or the specific gene knockout, whereas part of the metabolic pathways such as only EMP pathway and/or the PP pathway cannot do this. In the present approach, since we considered the whole main metabolic pathways, the specific ATP production rate could be estimated for substrate level phosphorylation and oxidative phosphorylation with respect to time, and this could be utilized to estimate the cell growth rate. Moreover, the specific CO2 production rate could also be computed, and this can be used to estimate the yield.
Figure 2b indicates that the specific NADPH production rate, v NADPH is growth rate-dependent. It is usually difficult to predict the PP pathway fluxes using only enzymatic equations. Once we could estimate the cell growth rate based on ATP production rate (catabolism) as stated above, the PP flux may be estimated based on the relationship between the cell growth rate and NADPH production rate (anabolism) as given in Figure 2b. This reflected in the flux computations. Here we assumed that Entner Doudoroff pathway was inactive and the fluxes of G6PDH and 6PGDH were equal at steady state without loss of generality. Note that the computational algorithm for computing the flux through oxidative PP pathway requires repetitive evaluation as follows: Once the specific glucose consumption rate was given, all the intracellular fluxes can be computed using kinetic models for the main metabolic pathways. Then the specific ATP production rate was estimated by Eq.(3) with (4), and the specific growth rate was estimated by Eq.(5), using the relationship found from Figure 2a. The v NADPH was computed by Eq.(8) with the relationship of Figure 2b. The fluxes through G6PDH and 6PGDH were estimated from Eq.(7). Since the mass balance at G6P changes by the newly determined flux through G6PDH, the intracellular fluxes from G6P to Pgi were recalculated again and so on. Note that the fluxes of the oxidative PP pathways such as G6PDH, 6PGDH and/or ICDH were estimated, NADPH/NADP+ concentrations can also be estimated if NADPH consumption rates at GDH (glutamate dehydrogenase) and other amino acids were taken into account.
As stated before, the predictability is one of the primary concerns for the simulation from the practical application point of view. Figure 5a and Additional file 4 indicate that Ppc knockout caused Pyk flux to be up-regulated and PEP concentration to be up-regulated. The increase of PEP concentration also caused GAP/DHAP concentration to be up-regulated. The Ppc knockout also caused OAA concentration to be significantly down regulated. Among the data investigated, only AcCoA concentration was different between the simulation and experimental results, where it became lower for the experiment, while it became higher for the simulation. This is due to the fact that AcCoA concentration increased by the increase in PDH flux, and the reduced flux of CS caused by low OAA concentration.
Different from Ppc knockout, Pck knockout resulted in significant reduction in the Ppc flux, and PEP concentration decreased in the simulation, as also shown by the experiment . Moreover, the PEP, PYR, and AcCoA concentrations together with ICit concentration were down-regulated. This may be due to the reduction of AcCoA concentration which acts as an activator of Ppc . Moreover, the glyoxylate pathway was activated in order to backup OAA.
As shown in Additional file 7, most of the glycolytic metabolites such as G6P/F6P, GAP/DHAP and PEP were higher for Pyk mutant, and the simulation result is consistent with the experimental data. The inhibition of Pfk by PEP is not so high, and the glycolysis fluxes became rather higher due to the increased PTS flux caused by the increase in PEP concentration. The increased flux of Ppc caused an increase in OAA and MAL concentrations. The concentration of PYR is supposed to be lower as observed in the experiment, but it shows the reverse result in the simulation. This may be due to the increase in PTS flux caused by the increase in PEP concentration. If we consider the mass balance for PYR, VPyk is zero, VPTS increased, while VMez and VPDH increased, and PYR concentration may decrease [11, 17–19], while the simulation result shows reverse. The concentrations of all the TCA cycle metabolites such as ICIT, αKG, FUM, and SUC increased in the simulation, whereas the experimental data indicates the reverse trends. This may be due to the increase in PYR for the simulation. In consistency with the experimental data [11, 17, 18], the simulation result indicates that Pyk knockout has little effect on the cell growth characteristics (Figure 7a.). However, significant changes were observed in the metabolite concentrations and fluxes inside the cell, where the increased PEP concentration increased the flux of Ppc, which caused the increase of OAA and MAL concentrations and Mez flux to backup PYR. Those are consistent with enzyme activity data  and the flux result obtained by the 13C-labeling experiment.
In the present investigation, we focused on the regulation at the important branch point of PEP in view of energy generation and anapleorotic requirement. This may be extended to others mutants. However, the present model cannot simulate single gene knockouts for enzymes encoded by multiple genes, such as the operons acnAB encoding Acn, sucAB, lpdA encoding αKGDH, sdhCDAB encoding SDH, and fumABC encoding Fum in the TCA cycle. This is because, the present model does not include such gene level regulation.
The simulation results indicate the usefulness in understanding the metabolic changes based on intracellular metabolite concentrations and fluxes, and how those affected the fermentation characteristics in response to specific pathway gene knockouts for both batch and continuous cultures. In conclusion, the present model could predict the metabolic changes due to these specific pathway knockout to some extent, and thus it is useful in understanding the metabolic regulation of the specific gene knockout mutant in practice.
Strains used and culture conditions
The strains used were Escherichia coli BW25113 (lacIq rrn BT14 ΔlacZWJ16 hsdR514 ΔaraBADAH22 ΔrhaBADLD78) and its ppc (JWK3928), pck (JWK3366) and pykF (JWK1666) mutants. The M9 synthetic medium was used for both batch and continuous cultivations, where it consists of 48 mM Na2HPO4, 22 mM KH2PO4, 10 mM NaCl, 40 mM (NH4)zSO4 and 4 g/l of glucose. The following components were filter-sterilized separatory and then added (per liter of final volume): 1 ml 1 M MgSO4, 1 ml vitamin B1 (1 mgl stock), 1 ml 0.1 mM CaCl2, and 10 ml trace element solution containing (per liter): 0.55g CaCl2 1g FeCl3, 0.1 mg/l MnCl2.4H2O, 0.17 g ZnCl2, 0.043 g CuCl2.2H2O, 0.06 g CoCl2.2H2O, and 0.06 g Na2MoO4.2H2O). Batch and continuous cultures were conducted at 37°C in a working volume of 1 l in a 2l fermentor (M-100, Tokyo Rikakikai, Japan) equipped with pH and temperature sensors. The air flow rate was maintained at 1 l min-1 and 300-400 RPM, which ensured the aerobic condition where the dissolved oxygen concentration was kept at 3-4 ppm. The pH of the culture was controlled at 7.0 ± 0.1 by automatic addition of either 2.0 M HCl or 2.0 M NaOH with a pH controller. The CO2 and O2 concentrations in the off-gas were measured by an off-gas analyzer (DZX-2562, Able Co. Japan).
Measurement for extracellular metabolites concentrations
The 5 ml of samples were taken from the culture broth and rapidly quenched into 15 ml of 60% (v/v) aqueous methanol containing 70 mM HEPES, where quenching solution was kept at -80°C before used. The cells were separated from the culture by centrifugation at 10,000 ×g for 15 min at 0°C. To extract intracellular metabolites from the cell pellet, 500 ml of 50% methanol was added, and the cells were re-suspended by vortexing the mixture. Then 2 ml of 35% of perchloric acid was added, which was pre-cooled on ice, and vortexed again for 10 s. After one freeze-thaw cycle, proteins and cell fragments were removed by centrifugation at 12,000 ×g for 30 min at 0°C. The clear suparnatant was neutralized by adding collected supernatant, and then 895 μm of 5 M K2KO3 was added. Afterward, precipitated perchlorates were removed by another centrifugation step (12,000 ×g at 0°C for 10 min), the clear supernatant was collected, and they were stored as 200 μm aliquots at -20°C for further analysis [20, 21].
The intracellular metabolite concentrations were determined by enzymatic methods, where those were performed mainly as described by [22, 23] and others [3, 20, 21, 24] after some modifications. Enzymatic determinations of metabolites were performed in 50 mM tri ethanolamine buffer (TEA) at pH 7.6. For example, G6P concentration was determined in the presence of 5 mM MgSO4 and 0.48 mM NADP+ after addition of 0.7 U ml-1 of G6PDH. The measurements were made using Luminescence Spectrophotometer (LS55, Perkin Elmer, UK). The chemical determination was made at 37°C .
13C-labeling experiments and sample preparation
The biomass sample was kept on ice for 2-3 minutes, and the sample was centrifuged at 6,000 rpm at 2°C for 15 minutes . The cell pellets were washed three times with 20 mM Tris-HCl at pH 7.6 and suspended in 10 ml of 6 M HCl. The mixture was then hydrolyzed at 105°C for 15 hours in a sealed glass tube. During acid hydrolysis, tryptophan and cysteine were oxidized, and glutamine and asparagine were deaminated. The hydrolysate was evaporated to dryness. The dried material was dissolved in milli-Q water and filtered through a 0.22 μm pore-size filter and evaporated to dryness. About 1.5 ml acetonitrile was added in the dried hydrolyte and incubated at room temperature overnight. After the color of liquid turned a color of pale yellow, it was filtrated through 0.22 μm pore-size filter. The filtrate was then derivatized by N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide (MTBST-FA) (Aldrich, USA) at 110C for 30 minutes and was transferred to GC-MS sample tube for analysis .
13C-labeling experiments were initiated after the culture reached steady state, which was inferred from the stable oxygen and carbon dioxide concentrations in the off-gas and stable OD in the effluent medium for at least twice as long as the residence time. The feed medium with 4 g/l of unlabeled glucose was then replaced by an identical medium containing 0.4 g of [U-13C] glucose, 0.4 g of [1-13C] glucose, and 3.2 g of naturally labeled glucose per liter. Biomass samples for GC-MS analysis were taken after one residence time. Sample preparation, analytical procedures for GC-MS analysis, and flux computation are given elsewhere [16–18, 25].
Metabolic flux analysis
Preparation of biomass hydrolysates and recording of the GC-MS spectra (PerklinElmer, Germany) were made as described previously [25–27]. The program Turbomass Gold (Perklin Elmer,Germany) was used for peak assignment and MS data processing. The main idea is to perform isotopomer balance on carbon atoms in order to track the fate of the labelled carbon atoms from the substrate . Isotopomer balance enables us to determine the isotopomer distributions of the intracellular metabolites in the central metabolic network. Since the isotopomer distribution of amino acids can be inferred from the isotopomer distributions of their corresponding precursors, the GC-MS signals for the amino acids can then be simulated. A set of flux distributions is then determined by minimizing the differences between the experimental and simulated data. For the estimation of GC-MS signals using isotopomer balance, three types of corrections were made to take into account the effects of natural abundance, non-steady state condition, and skewing effect. First, the isotopomer distributions of the input substrate were corrected for natural abundance in the unlabeled glucose and impurity in the labeled glucose. Second, the isotopic steady state condition is normally not attained at the time of harvesting biomass. Thus, the simulated data have to be corrected based on the actual harvesting time by assuming first order washout dynamics . Finally, the simulated GC-MS data were corrected for the natural isotope abundance of O, N, H, Si, S, and C atoms in the derivatizing agent using the correction matrices .
where [·] denotes the concentration, and the abbreviations for the metabolites are explained in Nomenclature. μ is the specific growth rate and v i are the intracellular fluxes. The superscript "ex" means extra cellular. Note that the last term with μ in each equation denotes the dilution effect due to the increase in cell volume as the cell grows . As shown in Figure 1, GAP and DHAP were lumped together. A sequence of enzymatic reactions by GAPDH, Pgk, Pgm and Eno can be considered to be in equilibrium, and thus it was assumed to be one lumped reaction from GAP to PEP, for simplicity. The overall simulation procedure is given in Additional file 8.
Modeling for the cell growth rate
where X is the cell concentration and X m is the final value of X in the batch culture, and n is the model parameter. X m may be set as S(0)Y S/X where S (0) is the initial substrate concentration, and Y S/X is the yield coefficient. Although this equation can be used to simulate the dynamic behavior in batch culture, Eq. (2b) may not be able express the cell growth rate for the mutants. We, therefore, introduce another idea by taking into account the effect of ATP generation on the cell growth rate based on experimental observation.
where: (P/O) and (P/O)' are the P/O ratios for NADH and FADH2, respectively.
Where v ATP (•) is the specific ATP production rate computed from the fluxes, and v ATP is the adjustable model parameter.
Kinetic equations for the enzyme kinetics
The kinetic equations for the various enzyme-catalyzed reactions considered in the present study are given in Table 1. The specific growth rate was assumed to be expressed by Eq.(1), where Eq.(1a) in Table 1 was used for the cell growth phase as state above, while Eq.(1b) was used for the late growth phase where acetate was used as a carbon source.
The glucose transport via phosphotransferase system (PTS) has been extensively investigated and is catalyzed by a sequence of enzymes such as EI, HPr, EIIAGlc and EIICBGlc. These were encoded by such genes as ptsHI, crr and ptsG. It may be assumed that EIICBGlc can be either phosphorylated first or bind with glucose first, although the phosphorylation of EIIB may be facilitated by binding EIIC to glucose. Based on these assumptions, the kinetic rate equation may be derived for the glucose transport through PTS as given by Eq.(2) in Table 1 , where Glcex and PEP are the substrates for the reaction, and PYR and G6P are the products, which inhibit the reaction rate as those are accumulated.
The phophoglucose isomerase (Pgi) is the first enzyme in glycolysis, and the two Pgi isoenzymes have been lumped into one rate equation . The phosphofructokinase (Pfk) is encoded by pfkA and pfkB in E. coli, where Pfk-1 encoded by pfkA is dominant, accounting for 90% of the total activity . Here, we, therefore, considered only Pfk-1, where it is inhibited by PEP , and its reaction rate may be expressed as given by Eq.(4) of Table 1, where F6P is the substrate for this reaction, and PEP inhibits vPfk as it accumulates . For FDP aldolase reaction, Eq.(5) of Table 1 was considered , where FDP is the substrate and GAP is the product. The minus term in the numerator is for the backward reaction. For GAPDH reaction, Eq.(6) of Table 1  was considered, where GAP is the substrate, and PEP is the product for this reaction. Note that the reactions from GAP to PEP were lumped into GAPDH reaction for simplicity. This reaction rate was assumed to be inhibited by NADH as [NADH]/[NAD] ratio increases. The pyruvate kinase (Pyk) reaction is catalyzed by two isoenzymes, PykI and PykII, encoded by pykF and pykA, respectively. PykI is activated by FDP and inhibited by ATP, whereas PykII is activated by AMP. Here we lumped these together, and the rate equations may be expressed as given by Eq.(7) of Table 1, where PEP and ADP are the substrates for this reaction and inhibited by ATP . Moreover, Eq.(7) indicates that FDP and AMP are the activators of the reaction.
It is quite important to correctly simulate the fluxes at the branch point of PEP, where the enzyme activities of Pyk and Ppc determine the fluxes. Based on the experimental data , it is known that the reaction catalyzed by Ppc exhibits a hyperbolic function with respect to PEP concentration, and the reaction rate is usually low without any activator, where AcCoA is a very potent activator, and FDP alone exhibits no activation, but it produces a strong synergistic activation with AcCoA. Based on these observations, Eq.(8) of Table 1 may be considered . Consider now the reaction in the reverse direction, namely Pck. The kinetics of the Pck in E. coli may be assumed to follow the rapid equilibrium mechanism . Thus, we may express the rate equation as given by Eq.(9) of Table 1, where OAA and ATP are the substrates for this reaction, and the reaction rate is inhibited by its products such as PEP and ADP .
For the rate equation of PDHc, the hyperbolic function (Michaelis-Menten type) was considered with respect to PYR, and it was assumed to be inhibited by NADH as [NADH]/[NAD] ratio increases as given by Eq.(10) of Table 1 . The acetate formation is of important concern in practice for the cell growth and the specific metabolite production in E. coli. The cells undergo a metabolism switch associated with the production and utilization of acetate. When enough glucose is present around the cell, the cell produces and excretes acetate. Once glucose is consumed or becomes low level, the cells utilize acetate instead of excreting it. This acetate-associated metabolic switch occurs just as the cells begin to decelerate growth, or fall into stationary phase in the batch culture. The first pathway for acetate production is Pta-Ack pathway, where the reaction catalyzed by Pta and Ack proceeds through an unstable, high energy, acetyl phosphate (AceP) intermediate. For the Pta reaction, Eq.(11) of Table 1 was considered [40, 41]. The rate equation for Ack was considered as expressed by Eq.(12) of Table 1 . During the stationary phase when glucose was depleted, the cells scavenge for the acetate by using the 2nd pathway of Acs. This reaction catalyzed by Acs proceeds through an enzyme-bound acetyladenylate (Acetyl-AMP) intermediate. Acs has high affinity to acetate while Ack-Pta has low affinity. The rate equation for Acs was assumed to be expressed by the simple Michaelis-Menten equation as given by Eq.(13) of Table 1 .
For citrate synthase reaction, Eq.(14) in Table 1  was considered, where AcCoA and OAA are the substrates, and the reaction is inhibited by NADH. Another important branch point is at isocitrate (ICIT) branch point in the TCA cycle, where ICDH (high affinity to ICIT, Km = 8 mM) and Icl (low affinity to ICIT, Km = 406 mM) compete for ICIT. This branch point is under both gene level regulation by aceBAK operon and the enzyme level regulation by the phosphorylation/dephosphorylation of ICDH. Here, we considered the rate equation for ICDH as given by Eq.(15) of Table 1 , where ICIT is the substrate, and the reaction is inhibited by NAD(P)H as [NAD(P)H]/[NAD(P)] ratio increases and as [αKG] increases. For KGDH reaction, Eq.(18) of Table 1 was considered, where αKG and CoA are the substrates, and the reaction is assumed to be inhibited by NADH as [NADH]/[NAD] ratio increases . For SDH reaction, Eq.(19) of Table 1 was considered, where SUC is the substrate and FUM is the product . For fumarase reaction, Eq.(20) of Table 1 was considered, where FUM is the substrate, and the reaction is inhibited by its product MAL . For MDH reaction, Eq.(21) of Table 1 was considered , where MAL is the substrate, and the reaction is inhibited by its products OAA and NADH as these concentrations increase. For malic enzyme reaction, Eq.(22) of Table 1 was considered .
As for the glyoxylate pathway, Eq. (16) of Table 1 was considered for Icl, ICIT being the substrate of reaction, inhibited by high accumulation of its products SUC and GOX . For malate synthase reaction, Eq.(17) of Table 1 was considered, where GOX and AcCoA are the substrates, and the reaction is inhibited by its product MAL .
For G6PDH and 6PGDH reactions in the PP pathway, Eqs.(23) and (24) of Table 1 were considered, where G6P is the substrate of the former reaction inhibited by NADPH, and 6PG is the substrate of the latter reaction inhibited by NADPH . Other equations for PP pathway such as Rpe, Rpi, TktA, TktB, and Tal were considered in the form of mass-action kinetic, given in Eqs. (25)-(29) in Table 1 .
and the total CO2 production rate was computed by the product of net production of CO2 and biomass concentration, i.e. [X].
Once we identified k NADPH value from Figure 2b, we may compute k NADPH after μ was computed. Then, the flux of G6PDH may be estimated.
As stated in introduction, it is not easy, even if it is not impossible, to express gene-level regulation by mathematical equations. Here, we alternatively introduced several rules for the simulation based on the experimental observations and the roles of global regulators such as given in Additional file 3: Table S1, where the rules 1-3 in Additional file 3: Table S1 came from the experimental observations, while Rule 4 and Rule 5 are known by gene level regulations [47, 48]. Rule 6 came from our previous analysis .
2KG: 2-Keto-D-gluconate; 6PG: 6-Phosphogluconolactone; ACE: Acetate; AcP: Acetyl phosphate; AcCoA: Acetyl-CoA; ASP: Aspartate; ADP: Adenosine diphosphate; ATP: Adenosine-5'-triphosphate; AMP: Adenosine monophosphate; DHAP: Dihydroxyacetone phosphate; E4P: Erythrose 4-phosphate; F6P: Fructose 6-phosphate; FDP: Fructose 1,6-bisphosphate; FUM: Fumarate; G6P: Glucose-6-phosphate; GAP: Glyceraldehyde 3-phosphate; ICIT: Isocitrate; NAD/NADH: Nicotinamide adenine dinucleotide; NADP/NADPH: Nicotinamide adenine dinucleotide phosphate; MAL: Malate; OAA: Oxaloacetate; P: Phosphate; PEP: Phosphoenolpyruvate; PYR: Pyruvate; R5P: Ribulose 5-phosphate; Ru5P: Ribose 5-phosphate; S7P: Sedoheptulose 7-phosphate; SUC: Succinate; X5P: Xylulose 5-phosphate
2KGDH: 2-Keto-D-gluconate Dehydrogenase; Ack: Acetate kinase; Acs: Acetyl coenzyme A synthetase; AcP: Acetyl phosphate; Aldo: Aldolase; CS: Citrate synthase; EI: Cytoplasmatic protein (enzyme1); EII, EIIB, EIIC,: Carbohydrate specific (enzyme II); EIIAGlc,/EIICBGlc: Enzyme for glucose; Eno: Enolase; Fum: Fumarase; G6PDH: Glucose-6-phosphate dehydrogenase; GAPDH: Glyceraldehyde 3-phosphate dehydrogenase; ICDH: Isocitrate dehydrogenase; HPr: Histidine containing protien; Icl: Isocitrate lyase; MDH: Malate dehydrogenase; Mez: Malic enzyme; Ms: Malate synthase; Pck: Phosphoenolpyruvate carboxykinase; PDH: Pyruvate dehydrogenase; Pfk/Pfk-1: Phosphofructokinase-1; Pgi: Phosphoglucose isomerase/Glucosephosphate isomerase; Pgk: Phosphoglycerate kinase; Pgm: Phosphoglycerate mutase; Ppc: PEP carboxylase; Pps: Phosphoenolpyruvate synthase; Pta: Phosphotransacetylase; PTS: Phosphotransferase system; Pyk: Pyruvate kinase; Rpe: Ribulose phosphate 3-epimerase; Rpi: Ribulose 5-phosphate 3-isomerase; SDH: Succinate dehydrogenase; Tal: Transaldolase; TktA: TransketolaseI; TktB: TransketolaseII
aceBAK: operon which encodes the metabolic and regulatory enzymes of the glyoxylate bypass; cra: Catabolite repressor/activator; crr: Catabolite repressor; fadR: Fatty acid metabolism regulator; icdA: Isocitrate dehydrogenase A; iclR: Transcriptional repressor IclR; pckA: Phosphoenolpyruvate carboxykinase gene; pfkA, pfkB: Phosphofructokinase gene; ppsA: Phosphoenolpyruvate synthase gene; ptsG, ptsHI: Pts gene; pykA, pykF: Pyruvate kinase gene
This research was supported in part by Strategic International Cooperative program, Japan Science and Technology Agency (JST) and Biotechnology and Biological Sciences Research Council (BBSRC).
- Tomita M, Hashimoto K, Takahashi K, Shimizu TS, Matsuzaki Y, Miyoshi F, Saito K, Tanida S, Yugi K, Venter JC, Hutchison CA: E-CELL: Software environment for whole-cell simulation. Bioinformatics. 1999, 15: 72-84.View ArticleGoogle Scholar
- Rizzi M, Baltes M, Theobald U, Reuss M: In vivo analysis of metabolic dynamic in Saccharomyces cerevisiae: II. Mathematical model. Biotechnol Bioeng. 1997, 55: 592-608.View ArticleGoogle Scholar
- Theobald U, Mailinger W, Baltes M, Rizzi M, Reuss M: In vivo analysis of metabolic dynamic in Saccharomyces cerevisiae: I. Experimental observations. Biotechnol Bioeng. 1997, 55: 305-316.View ArticleGoogle Scholar
- Vaseghi S, Baumeister A, Rizzi M, Reuss M: In vivo dynamics of the pentose phosphate pathway in Saccharomyces cerevisiae. Metab Eng. 1999, 1: 128-140.View ArticleGoogle Scholar
- Buziol S, Bashir I, Baumeister A, Claassen W, Rizzi M, Mailinger W, Reuss M: A new bioreactor coupled rapid stopped-flow-sampling-technique for measurement of transient metabolites in time windows of milliseconds. Proceeding of the Fourth International Congress on Biochemical Engineering. Stuggart: Fraunhofer IRB. 2000, 79-83.Google Scholar
- Buziol S, Bashir I, Baumeister A, Claassen W, Rizzi M, Mailinger W, Reuss M: A new bioreactor coupled rapid stopped-flow sampling technique for measurement of metabolite dynamics on a subsecond time scale. Biotecnol Bioeng. 2001, 80: 632-636.View ArticleGoogle Scholar
- Theobald U, Mailinger W, Reuss M, Rizzi M: In vivo analysis of glucose induced fast changes in yeast adenine nucleotide pool applying a rapid sampling technique. Anal Biochem. 1993, 214: 31-37.View ArticleGoogle Scholar
- Chassagnole C, Noisommit-Rizzi N, Schmid JW, Mauch K, Reuss M: Dynamic modeling of the central carbon metabolism of Escherichia coli. Biotech and Bioeng. 2002, 79: 53-73.View ArticleGoogle Scholar
- Kotte O, Zaugg JB, Heinemann : Bacterial adaption through distributed senseing of metabolic fluxes. Molecular Systems Biology. 2010, 6: 355-View ArticleGoogle Scholar
- Usuda Y, Nishio Y, Iwatani S, Van Dien SJ, Imaizumi A, Shimbo K, Kageyama N, Iwahata D, Miyano H, Matsui K, Usuda Y, Nishio : Dynamic modeling of Escherichia coli metabolic and regulatory systems for amino-acid production. J Biotechnol. 2010, 147 (1): 17-30.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. 2010Google Scholar
- Shimizu K: Toward systematic metabolic engineering based on the analysis of metabolic regulation by the integration of different levels of information. Biochem Eng. 2010, 46: 235-251.View ArticleGoogle Scholar
- Wittmann C: Fluxome analysis using GC-MS. Microb Cell Factories. 2007, 6: 6-View ArticleGoogle Scholar
- Ishii N, Nakahigashi K, Baba T, Robert M, Soga T, Kanai A, Hirasawa T, Naba M, Hirai K, Hoque A, Ho PY, Kakazu Y, Sugawara K, Igarashi S, Harada S, Masuda T, Sugiyama N, Togashi T, Hasegawa M, Takai Y, Yugi K, Arakawa K, Iwata N, Toya Y, Nakayama Y, Nishioka T, Shimizu K, Mori H, Tomita M: Multiple high-throughput analyses monitor the response of E. coli to perturbations. Science. 2007, 316: 593-597.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. Biotech and Bioeng. 2003, 84: 129-144.View ArticleGoogle Scholar
- Peng L, Arauzo-Bravo MJ, Shimizu K: Metabolic flux analysis for a ppc mutant Escherichia coli based on 13C-labelling experiments together with enzyme activity assays and intracellular metabolite measurements. FEM Microbiology Letters. 2004, 235: 17-23.View ArticleGoogle Scholar
- Siddiquee KAZ, Arauzo-Bravo MJ, Shimizu K: Metabolic flux analysis of pykF gene knockout Escherichia coli based on 13C labeled experiment together with measurements of enzyme activities and intracellular metabolite concentrations. Appl Microbiol Biotechnol. 2004, 63: 407-417.View ArticleGoogle Scholar
- Siddique KAZ, Arauzo-Bravo MJ, Shimizu K: Effect of a pyruvate kinase (pykF-gene) knockout mutation on the control of gene expression and metabolic fluxes in Escherichia coli. FEMS Microbiology Letters. 2004, 235: 25-33.View ArticleGoogle Scholar
- Hoque MA, Ushiyama H, Tomita M, Shimizu K: Dynamic responses of the intracellular metabolite concentrations of the wild type and pykA mutant Escherichia coli against pulse addition of glucose or NH3 under those limiting continuous cultures. Biochemical Eng Journal. 2005, 26: 38-49.View ArticleGoogle Scholar
- de Koning W, van Dam K: A method for the determination of changes in glycolytic metabolites in yeast on a subsecond time scale using extraction at neutral pH. Anal Biochem. 1992, 204: 118-123.View ArticleGoogle Scholar
- Schaefer U, Boos W, Takors R, Weuster-Botz D: Automated sampling device for monitoring intracellular metabolite dynamics. Anal Biochem. 1999, 270: 88-96.View ArticleGoogle Scholar
- Bergmeyer HU: Methods of Enzymatic Analysis. 1984, 6: Verlag Chemie, Weinheim, Germany, thirdGoogle Scholar
- Bergmeyer HU: Methods of Enzymatic Analysis. 1984, 7: Verlag Chemie, Weinheim, Germany, thirdGoogle Scholar
- Sidney PC, Kaplan NO: Methods in Enzymology, Citric Acid Cycle. Edited by: Lowenstein JM. 1969, XIII: Academic PressGoogle Scholar
- Zhao J, Shimizu K: Metabolic flux analysis of E. coli K12 grown on 13C-labeled acetate and glucose using GC-MS and powerful flux calculation method. J Biotechnol. 2003, 101: 101-117.View ArticleGoogle Scholar
- Zhao J, Baba T, Mori H, Shimizu K: Effect of zwf Gene Knockout on the Metabolism of Escherichia coli Grown on Glucose or Acetate. Metabolic Eng. 2004, 6: 164-174.View ArticleGoogle Scholar
- Zhao J, Baba T, Mori H, Shimizu K: Analysis of metabolic and physiological responses to gnd knockout in E. coli by using 13C tracer experiment and enzyme activity measurement. FEMS Microb Lett. 2004, 220: 295-301.Google Scholar
- Szyperski T: Biosynthetically directed fractional 13C-labeling of proteinogenic amino acids. An efficient analytical tool to investigate intermediary metabolism. Eur J Biochem. 1995, 232: 433-448.View ArticleGoogle Scholar
- Van Winden WA, Heijnen JJ, Verheijen PJT: Cumulative bondomers: a new concept in flux analysis from 2D[13C, 1H] COSY NMR data. Biotechnol And Bioeng. 2002, 80: 731-745.View ArticleGoogle Scholar
- Wittmann C, Heinzle E: Mass spectrometry for metabolic flux analysis. Biotechnol And Bioeng. 1999, 62: 739-750.View ArticleGoogle Scholar
- Gennis RB, Stewart V: Respiration. In Escherichia coli and Salmonella. Edited by: Neidhardt FC. 1996, 217-261. ASM Press: Washington, DCGoogle Scholar
- Harold FM, Maloney PC: Escherichia coli and Salmonella. Edited by: Neidhardt FC. 1996, 283-306. ASM Press, Washington, D. CGoogle Scholar
- Liao JC, Hou SH, Chao YP: Pathway analysis, engineering and physiological considerations for redirecting central metabolism. Biotech Bioeng. 1996, 52: 129-140.View ArticleGoogle Scholar
- Schreyer R, Bock A: Phosphoglucose isomerase from Escherichia coli K 10:purification, properties and formation under aerobic and anaerobic. Arch Microbiol. 1980, 127 (3): 289-298.View ArticleGoogle Scholar
- Karp P, Riley M, Paley S, Pellegrini-Toole A, Krummenacker M: EcoCyc: Encyclopedia of Escherichia coli genes and metabolism. Nucleic Acids Res. 1997, 25: 43-50.View ArticleGoogle Scholar
- Kotlarz D, Garreau H, Buc H: Regulation of the amount and of the activity of phosphofructokinase and pyruvate kinase in Escherichia coli. Biochem Biophys Acta. 1975, 381 (2): 257-268.View ArticleGoogle Scholar
- Izui K, Taguchi M, Morikawa M, Katsuki H: Regulation of Escherichia coli phosphoenolpyruvate carboxylase by multiple effectors in vivo. Kinetic studies with a reaction system containing physiological concentrations of ligands. J Biochem. 1981, 90: 1321-1331.Google Scholar
- Lee B, Yen J, Yang L, Liao JC: Incorporating qualitative knowledge in enzyme kinetic models using fuzzy logic. Biotech and Bioeng. 1999, 62: 722-729.View ArticleGoogle Scholar
- Krebs A, Bridger WA: The kinetic properties of phosphoenolpyruvate carboxykinase of Escherichia coli. J Biochem. 1980, 58: 309-318.Google Scholar
- Hoefnagel MHN, Starrenburg MJC, Martens DE, Hugenholtz J, Kleerebezem M, Van Swam II, Bongers R, Westerhoff HV, Soep JL: Metabolic engineering of lactic acid bacteria, the combined approach: kinetic modeling, metabolic control and experimental analysis. Microbiology. 2002, 148: 1003-1013.View ArticleGoogle Scholar
- Henkin J, Abales RH: Evidence against an acyl-enzyme intermediate in the reaction catalyzed by clostridial phosphotransacetylase. Biochemistry. 1976, 15: 3475-3479.View ArticleGoogle Scholar
- Fung E, Wong WW, Suen JK, Bulter T, Lee Sun-gu, Liao JC: A synthetic gene-metabolism oscillator. Nature. 2005, 435: 118-122.View ArticleGoogle Scholar
- Mogilevskaya EA, Lebedeva GV, Demin OV: Kinetic model of E. coli citrate syntase functioning. Proceeding of the 12th Int Conf Mathematics Computer Education. 2006, V3: 934-937.Google Scholar
- Mogilevskaya EA, Lebedeva GV, Goryanin II, Demin OV: Kinetic model of functioning and regulation of Escherichia coli isocitrate dehydrogenase. Biophysics. 2007, 52: 30-39.View ArticleGoogle Scholar
- Wright BE, Butler MH, Albe KR: System analysis of the tricarboxylic acid cycle in the Dictyostelium discoideum I. The basic for model construction. J of Bio Chem. 1992, 267: 3101-3105.Google Scholar
- Singh VK, Ghosh I: Kinetic modeling of tricarboxylic acid cycle and glyoxylate bypass in Mycobacterium tuberculosis, and its application to assessment of drug target. Theor Biol and Med Model. 2006, 3: 27-View ArticleGoogle Scholar
- Peng L, Shimizu K: Effect of fadR gene knockout on the metabolism of Escherichia coli based on analyses of protein expressions, enzyme activities and intracellular metabolite concentrations. Enzyme Microbiol Technol. 2006, 38: 512-520.View ArticleGoogle Scholar
- Sarkar D, Siddiquee KA, Arauzo-Bravo MJ, Oba T, Shimizu K: Effect of cra gene knockout together with edd and iclR genes knockout on the metabolism in Escherichia coli. Arch Microbiol. 2008, 190: 559-571.View ArticleGoogle Scholar
- Ishii N, Suga Y, Hagiya A, Watanabe H, Mori H, Yoshino M, Tomita M: Dynamic Simulation of an in vitro multi-enzyme system. FEBS Let. 2006, 581: 413-420.View ArticleGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.