# Elementary Mode Analysis for the Rational Design of Efficient Succinate Conversion from Glycerol byEscherichia coli

Authors: Zhen Chen; Hongjuan Liu; Jianan Zhang; Dehua Liu
Journal: Journal of Biomedicine and Biotechnology
 (2010)

---

## Abstract

By integrating the restriction of oxygen and redox sensing/regulatory system, elementary mode analysis was used to predict the metabolic potential of glycerol for succinate production byE. coli under either anaerobic or aerobic conditions. It was found that although the theoretical maximum succinate yields under both anaerobic and aerobic conditions are 1.0 mol/mol glycerol, the aerobic condition was considered to be more favorable for succinate production. Although increase of the oxygen concentration would reduce the succinate yield, the calculation suggests that controlling the molar fraction of oxygen to be under 0.65 mol/mol would be beneficial for increasing the succinate productivity. Based on the elementary mode analysis, the rational genetic modification strategies for efficient succinate production under aerobic and anaerobic conditions were obtained, respectively. Overexpressing the phosphoenolpyruvate carboxylase or heterogonous pyruvate carboxylase is considered to be the most efficient strategy to increase the succinate yield.


---

## Body

## 1. Introduction

Glycerol has become an abundant and inexpensive carbon source due to its generation as an inevitable byproduct of biodiesel production. Over the past few years, the price of crude glycerol has decreased 10-fold due to the tremendous growth of the biodiesel industry [1]. Much effort has been paid for the development of processes to convert crude glycerol into higher-value products to maximize the full economic potential of biodiesel process. For example, the transformation of glycerol into 1,3-propanediol has been extensively studied in the past few years [2–4].Several recent studies also tried to utilize glycerol as a carbon source for the transformation of other valued products such as ethanol [5] and amino acids [6]. Succinate is traditionally produced from sugars and suffers the limitation due to the availability of reducing equivalents. Compared with glucose, glycerol has a higher reduced state and also several microorganisms such as E. coli can transform glycerol into succinate [7]. So, the byproduct glycerol is a potential substrate for the succinate production. Despite few attempts in the past, no industrially competitive organisms can effectively produce succinate from glycerol so far. In the light of the new powerful tools of metabolic engineering, the quest for targeted development of strains that can effectively utilize glycerol for succinate production is strongly revived. E. coli is one of the most promising organisms since it can directly utilize glycerol and it has been traditionally developed for succinate production [8–10].The dissimilation of glycerol inE. coli is catalyzed by proteins encoded by glp regulon under aerobic conditions. Glycerol is first phosphorylated into glycerol 3-phophate (G3P) by ATP-dependent glycerol kinase encoded by glpK gene, and then glycerol 3-phophate is converted into dihydroxyacetone phosphate (DHAP) by aerobic G3P dehydrogenase encoded by glpD gene (Figure 1) [11, 12].Figure 1
Central metabolic network of glycerol in wildtypeE. coli. The dashed arrows represent the particular pathways in anaerobic conditions. Reversible reactions are represented by a double-headed arrow. Key genes associated with the pathway are included.Although the aerobic utilization of glycerol byE. coli has been known for a long time, the fermentative pathway of glycerol has just been clarified recently [13, 14]. It has been suggested that the feasibility of fermenting glycerol into fuels and other reduced chemicals is through the inducing of its native 1,2-propanediol fermentative pathway without using external electron acceptors. In this pathway, glycerol is converted to dihydroxyacetone (DHA) by NAD+-linked glycerol dehydrogenase (GDH), and the DHA is phosphorylated to DHAP via the ATP-dependent or phosphoenolpyruvate (PEP)-dependent DHA kinase (DHAK). DHAP is then reduced into 1,2-propanediol or enter glycolysis [15].SinceE. coli can utilize glycerol in both aerobic and anaerobic conditions, it is necessary to analyze the potential and the feasibility of engineering E. coli for the succinate production in either condition. A careful metabolic pathway analysis is very helpful in such kind of estimation and rational strain development. Elementary mode analysis is one of the most powerful tools for metabolic pathway analysis using for the metabolic properties study of cellular systems [16–18]. Elementary mode analysis allows the calculation of a solution space that contains all possible steady-state flux distributions of a network. The stoichiometry of the metabolic network, including carbon as well as cofactor requirements, is fully considered in elementary mode analysis. On the other hand, it also allows determining the overall capacity, that is, theoretical maximum yield, of a cellular system and studying the effects of any genetic modification. Based on such studies, rational design can be obtained for the efficient production and genetic modification. Moreover, knowledge of the theoretical maximum yield allows estimating the potential economic efficiency of a process. Recently, elementary mode analysis has been used for genome scale metabolic studies dealing with, for example, the rational design of methionine production in E. coli and C. glutamicum [19], the production of polyhydroxybutanoate in yeast [20], and growth-related aspects in Saccharomyces cerevisiae [21, 22] and E. coli [23, 24].In this work, the elementary mode analysis was carried out for succinate production by comparing the metabolic networks ofE. coli in anaerobic and aerobic conditions. The pathways involved in the flux scenario representing optimal succinate production were investigated and the effect of oxygen level on succinate production and biomass was also discussed. Furthermore, the rational design for the genetic modification of the E. coli to enhance the succinate production was developed. This work is considered useful for the further strain improvement and metabolic regulation in the succinate production by E. coli with glycerol as the substrate.


## 2. Materials and Methods

### 2.1. Metabolic Reaction Network

The glycerol metabolic network ofE. coli was constructed (Figure 1) based on KEGG database (http://www.genome.jp/kegg/metabolism.html) as well as biochemical and physiological literatures [13, 14, 25]. It includes glycerol dissimilation pathways, glycolysis pathway (EMP), pentose phosphate pathway (PPP), tricarboxylic acid (TCA) cycle, biosynthesis pathway, anaplerosis, and respiratory chain. The metabolic network was depicted in Figure 1. For the interconversion of NADH and NADPH, a cytosolic transhydrogenase transferring protons from NADPH to NAD+ and a membrane-bound transhydrogenase reducing NADP+ by oxidation of NADH were implemented [26]. For ATP production in the respiratory chain, a P/O ratio of 2 for NADH was assumed [27].The precursor demand for biomass formation was calculated according to the literature [28]. The biomass term is represents as Cmol basis (CN0.24S0.008). The cell physiology of E. coli is strongly affected by oxygen from different levels such as the transcriptional regulatory which cannot be represented simply in the metabolic network. For example, one component fumarate and nitrate reduction (FNR) protein is aerobic/anaerobic response regulator [25]. FNR appears to sense oxygen directly through a redox-sensitive iron-sulphur cluster in the protein and is active only during anaerobic growth. The two iron-sulphur ([4Fe–4S]2+) clusters in the dimeric FNR protein are converted to two [2Fe–2S]2+ clusters upon exposure to stoichiometric levels of oxygen. Active FNR protein activates and represses target genes in response to anaerobiosis. It acts as a positive regulator of genes expressed under anaerobic fermentative conditions such as aspartase, formate dehydroganases, fumarate reductase, and pyruvate formate lyase. To account for the effects of FNR regulator and other experimental discovery, the following constraints are used to discriminate the metabolic networks under aerobic and anaerobic network.For the anaerobic model, glycerol is assumed to be dissimilated into DHA by NAD+-dependent glycerol dehydrogenase and then DHA is phosphorylated into DHAP by ATP-dependent or PEP-dependent DHA kinase [13, 14]. The pathway from DHAP to 1,2-propanediol is considered to be active under anaerobic condition [13]. The pyruvate dehydrogenase complex is inactive under anaerobic condition and thus pyruvate-formate lyase was the only active enzyme that catalyzes the transformation of pyruvate into acetyl-CoA [29]. The TCA cycle is broken at the alpha-ketoglutarate dehydrogenase step and the respiratory chain is assumed to be inactive [30]. The detailed description of the model is listed in Appendices A.1 and A.2.For the aerobic model, glycerol is firstly phosphorylated into G3P by ATP-dependent glycerol kinase and then G3P is transferred into DHAP by NAD+-dependent G3P dehydrogenase [11, 12]. The 1, 2-propanediol pathway is assumed to be inactive. Pyruvate oxidase is active under aerobic condition which will transfer pyruvate into acetate. The detailed description of the model is listed in Appendices A.1 and A.3.


### 2.2. Computational Methods

In the present work, the elementary mode analysis was carried out for studying the aerobic and anaerobic metabolism of glycerol inE. coli by using METATOOL 5.1 [31]. The script files and compiled shared library of METATOOL 5.1 can be downloaded from the METATOOL website (http://pinguin.biologie.uni-jena.de/bioinformatik/networks/metatool/metatool5.1/metatool5.1.html). The mathematical details of the algorithm were described elsewhere [32]. Metabolic pathway analysis resulted in tens to hundreds of elementary flux modes for each situation investigated. For each of these flux modes, the fluxes were calculated as relative molar values normalized to the glycerol uptake rate and were expressed as mol/mol (glycerol).


## 2.1. Metabolic Reaction Network

The glycerol metabolic network ofE. coli was constructed (Figure 1) based on KEGG database (http://www.genome.jp/kegg/metabolism.html) as well as biochemical and physiological literatures [13, 14, 25]. It includes glycerol dissimilation pathways, glycolysis pathway (EMP), pentose phosphate pathway (PPP), tricarboxylic acid (TCA) cycle, biosynthesis pathway, anaplerosis, and respiratory chain. The metabolic network was depicted in Figure 1. For the interconversion of NADH and NADPH, a cytosolic transhydrogenase transferring protons from NADPH to NAD+ and a membrane-bound transhydrogenase reducing NADP+ by oxidation of NADH were implemented [26]. For ATP production in the respiratory chain, a P/O ratio of 2 for NADH was assumed [27].The precursor demand for biomass formation was calculated according to the literature [28]. The biomass term is represents as Cmol basis (CN0.24S0.008). The cell physiology of E. coli is strongly affected by oxygen from different levels such as the transcriptional regulatory which cannot be represented simply in the metabolic network. For example, one component fumarate and nitrate reduction (FNR) protein is aerobic/anaerobic response regulator [25]. FNR appears to sense oxygen directly through a redox-sensitive iron-sulphur cluster in the protein and is active only during anaerobic growth. The two iron-sulphur ([4Fe–4S]2+) clusters in the dimeric FNR protein are converted to two [2Fe–2S]2+ clusters upon exposure to stoichiometric levels of oxygen. Active FNR protein activates and represses target genes in response to anaerobiosis. It acts as a positive regulator of genes expressed under anaerobic fermentative conditions such as aspartase, formate dehydroganases, fumarate reductase, and pyruvate formate lyase. To account for the effects of FNR regulator and other experimental discovery, the following constraints are used to discriminate the metabolic networks under aerobic and anaerobic network.For the anaerobic model, glycerol is assumed to be dissimilated into DHA by NAD+-dependent glycerol dehydrogenase and then DHA is phosphorylated into DHAP by ATP-dependent or PEP-dependent DHA kinase [13, 14]. The pathway from DHAP to 1,2-propanediol is considered to be active under anaerobic condition [13]. The pyruvate dehydrogenase complex is inactive under anaerobic condition and thus pyruvate-formate lyase was the only active enzyme that catalyzes the transformation of pyruvate into acetyl-CoA [29]. The TCA cycle is broken at the alpha-ketoglutarate dehydrogenase step and the respiratory chain is assumed to be inactive [30]. The detailed description of the model is listed in Appendices A.1 and A.2.For the aerobic model, glycerol is firstly phosphorylated into G3P by ATP-dependent glycerol kinase and then G3P is transferred into DHAP by NAD+-dependent G3P dehydrogenase [11, 12]. The 1, 2-propanediol pathway is assumed to be inactive. Pyruvate oxidase is active under aerobic condition which will transfer pyruvate into acetate. The detailed description of the model is listed in Appendices A.1 and A.3.


## 2.2. Computational Methods

In the present work, the elementary mode analysis was carried out for studying the aerobic and anaerobic metabolism of glycerol inE. coli by using METATOOL 5.1 [31]. The script files and compiled shared library of METATOOL 5.1 can be downloaded from the METATOOL website (http://pinguin.biologie.uni-jena.de/bioinformatik/networks/metatool/metatool5.1/metatool5.1.html). The mathematical details of the algorithm were described elsewhere [32]. Metabolic pathway analysis resulted in tens to hundreds of elementary flux modes for each situation investigated. For each of these flux modes, the fluxes were calculated as relative molar values normalized to the glycerol uptake rate and were expressed as mol/mol (glycerol).


## 3. Results and Discussion

### 3.1. Elementary Mode Analysis of Glycerol Metabolism under Anaerobic Condition

Under anaerobic condition, the metabolic network model got 55 elementary flux modes. The relationship between the yields of products and biomass was shown in Figure2. The maximum molar yield of biomass under anaerobic condition is 0.187 mol/mol in which the respective yields of 1,2-propanediol, ethanol, and formate were 0.248 mol/mol, 0.495 mol/mol, and 0.57 mol/mol and no succinate, acetate, and lactate were produced. It was found that the cell growth was always associated with the production of 1,2-propandediol, ethanol, and formate. Therefore, the production of 1,2-propanediol, ethanol, and formate was necessary for the biomass synthesis during the glycerol metabolism. The biomass synthesis process consumes ATP and produces reducing equivalents (NADH). Both ATP and NAD need to be regenerated through the production of other byproducts (the biomass synthesis equation in Appendix A.2). For the glycerol metabolism in anaerobic condition, only the glycerol to 1,2-propanediol pathway can consume extra NADH (see (1)) and thus provide the mean to consume the reducing equivalents generated during the synthesis of biomass. The conversion of glycerol to ethanol and formate (see (2)), a redox-balanced pathway, fulfills energy requirements by generating ATP via substrate-level phosphorylation (Appendix A.4). The calculation results were consistent with the experiment observation that 1,2-propanediol and ethanol were growth-associated products [13, 14].
(1)glycerol+NADH→1,2-propanediol+NAD,(2)glycerol+ADP→ethanol+formate+ATP.Relationship between the yields of biomass and byproducts for the obtained elementary modes ofE. coli under anaerobic conditions. (a) Succinate, (b) Lactate, (c) Acetate, (d) 1,2-Propanediol, (e) Formate, and (f) Ethanol. The enclosed regions represent the possible solution space. The fluxes were normalized by glycerol uptake rate and expressed as mol/mol (glycerol).
(a)(b)(c)(d)(e)(f)The maximum succinate yield under anaerobic condition is 1.0 mol/mol when CO2 or carbonate salts are added as cosubstrates and the optimal flux distribution for succinate production was shown in Figure 3. In this case, there was no production of biomass and other byproducts. The key points for this mode were that the phosphorylation of DHA was only catalyzed by ATP-dependent DHA kinase and PEP was totally carboxylated into oxaloacetate by PEP carboxylase and the latter was further transferred into succinate. This required a very high activity of ATP-dependent DHA kinase and PEP carboxylase. However, it was reported that the PEP-dependent DHA kinase plays the main role in E. coli which dramatically reduced the yield of succinate. With single PEP-dependent DHA kinase function, there will be no succinate production [14]. Thus the PEP-dependent DHA kinase is the bottleneck of succinate production under anaerobic condition.Figure 3
The optimum flux distribution of glycerol metabolism for succinate production inE. coli under anaerobic conditions when only ATP-dependent DHA kinase plays a function.


### 3.2. Element Ary Mode Analysis of Glycerol Metabolism under Aerobic Condition

The aerobic metabolic network model got 259 elementary flux modes. The relationship between the yields of products and biomass was shown in Figure5. It was indicated that the flux distribution modes under aerobic condition were completely different from that of anaerobic condition. The maximum molar yield of biomass under aerobic condition was 0.725 mol/mol in which only CO2 was produced. The results suggested that aerobic condition was more favorable for biomass formation and the most effective way for biomass formation was through the TCA cycle. The cell growth process was no longer associated with the production of 1,2-propanedio, formate, and ethanol. Comparing with the anaerobic condition, the succinate production modes were more diversely distributed throughout the feasible solution space, and the corresponding yield of succinate at the same biomass yield was higher which indicated that the high potential of succinate production is associated with the cell growth by metabolic modification under aerobic condition.The maximum succinate yield under aerobic condition (an aerobic mode suggests oxygen consumption) is also 1.0 mol/mol which requires the CO2 or carbonate salts to be added as cosubstrates, and the optimal flux distribution for succinate production was shown in Figure 6. The optimal flux distribution modes were quite similar for the aerobic and anaerobic conditions. The key point for obtaining high succinate yield was considered as that PEP was totally carboxylated into oxaloacetate by PEP carboxylase and the latter was further transferred into succinate. This required a very high activity of PEP carboxylase.The network robustness and its sensitivity to perturbation were critical to the optimal metabolic pathway. The sensitivity of succinate yield to the flux ratios at the key branch nodes PEP and acetyl-CoA was considered in this work. PEP was consumed in reactions R40 and R17 which are catalyzed by PEP carboxylase and Pyruvate kinase, respectively; the flux ratio was denoted as(3)R(40,17)=R40(R40+R17).The influence of flux distribution at PEP node was shown in Figure7(a). The succinate yield increased from 0.5 mol/mol to 1.0 mol/mol when R(40,17) increased from 0 to 1.0. The results showed that increasing the flux distribution from PEP to oxaloacetate was beneficial for succinate production. For acetyl-CoA node, it was consumed in reactions R22, R44, R50, and R53 which are catalyzed by citrate synthase, aldehyde dehydrogenase, malate synthase, and phosphate acetyltransferase, respectively. The flux ratio was denoted as
(4)R(22,44,50,53)=R22(R22+R44+R50+R53).The influence of flux distribution at acetyl-CoA node was shown in Figure7(b). The succinate yield decreased from 1.0 mol/mol to 0.75 mol/mol when R(22,44,50,53) increased from 0 to 1.0. The results indicated that decreasing the flux distribution from acetyl to TCA cycle was beneficial for succinate production since the carbon would be lost during the TCA cycle.


### 3.3. Effect of Oxygen Consumption on the Production of Succinate and Biomass

The effect of oxygen consumption on succinate and biomass production was further investigated according to the elementary mode analysis; the relationship between the molar fractions of oxygen consumption and succinate or biomass under aerobic conditions was calculated and shown in Figure8. The theoretical succinate yield decreased as the oxygen consumption ratio increased (Figure 8(a)). The elementary modes of succinate yield distributed only on the left part of the solution space, which indicated that the higher consumption of oxygen was unfavorable for the succinate production. This is reasonable from the redox consideration. The higher consumption of oxygen results in more NAD which is critical for succinate production, to be utilized for the ATP production through oxidative phosphorylation. Especially as the molar oxygen consumption fraction of oxygen was more than 1.75 mol/mol, the elementary mode that produced succinate did not exist.The theoretical biomass yield increased firstly and then decreased when the oxygen consumption ratio increased (Figure8(b)). This is reasonable and consistent with the results of electron conservation because the rational increase of oxygen flux would be favorable for ATP synthesis through oxidative phosphorylation which is essential for biomass synthesis. However, high oxygen flux also results in less carbon source and reducing equivalents available for biomass because more carbon source would be oxidized to CO2. The maximum biomass was achieved when the molar fraction of oxygen was 0.65 mol/mol. Since higher cell concentration is beneficial for increase of the productivity, the optimal regulation strategy is controlling the molar fraction of oxygen consumption less than 0.65 mol/mol.


### 3.4. Rational Design to Improve the Succinate Production by Genetic Modifications

Comparing the results of elementary flux mode analysis above, although the maximum succinate yields were 1.0 mol/mol under both anaerobic and aerobic conditions, the aerobic condition seemed to be more favorable for succinate production in fact. The slow cell growth under anaerobic condition hindered the practical application of glycerol fermentation for succinate production [33]. What is more, the association of cell growth with the production of 1,2-propanediol, ethanol, and formate reduced the total yield of succinate.As discussed above, to improve the yield of succinate yield under anaerobic condition, the substitution of PEP-dependent DHA kinase into ATP-dependent DHA kinase and overexpressing the PEP carboxylase would be a prior consideration. An alternative choice is to express the heterogeneous pyruvate carboxylase inE. coli. The overexpression of pyruvate carboxylase could redistribute the flux of pyruvate into oxaloacetate for succinate production. The optimal flux distribution for succinate production in such case could also reach 1.0 mol/mol which was shown in Figure 4. It has been reported that expressing the ATP-dependent DHA kinase and pyruvate carboxylase could both increase the yield of succinate [14]. Another consideration to increase the succinate yield is reducing the byproducts production. Since the production of 1,2-propanediol, ethanol, and formate is necessary for biomass synthesis, the proper strategy is reducing the acetate and lactate production. Knockout of the phosphate acetyltransferase gene (pta) and lactate dehydrogenase gene (ldh) would be expected to increase the succinate production.Figure 4
The alternative optimum flux distribution of glycerol metabolism for succinate production inE. coli under anaerobic conditions when only PEP-dependent DHA kinase plays a function and heterogeneous PEP carboxylase (pyc) are introduced and overexpressed in E. coli.Relationship between the yields of biomass and byproducts for the obtained elementary modes ofE. coli under aerobic conditions. (a) Succinate, (b) Lactate, (c) Acetate, and (d) Ethanol. The enclosed regions represent the possible solution space. The fluxes were normalized by glycerol uptake rate and expressed as mol/mol (glycerol).
(a)(b)(c)(d)Figure 6
The optimum flux distribution of glycerol metabolism for succinate production inE. coli under aerobic conditions.Sensitivity of succinate yield to the relative fluxes at the (a) PEP node and (b) AcCoA node under aerobic conditions. The PEP node involves the catabolic reactions of R40 and R17 which are catalyzed by PEP carboxylase and pyruvate kinase, respectively. The AcCoA node involves the catabolic reactions of R22, R50, R44, and R53 which are catalyzed by citrate synthase, aldehyde dehydrogenase, malate synthase, and phosphate acetyltransferase, respectively.
(a)(b)Effect of oxygen consumption on the succinate production and biomass formation under aerobic conditions. The enclosed regions represent the possible solution space. The fluxes were normalized by glycerol uptake rate and expressed as mol/mol (glycerol).
(a)(b)For the case of aerobic condition, overexpressing the PEP carboxylase or expressing pyruvate carboxylase would be a prior consideration as discussed above. Knockdown of the isocitrate dehydrogenase gene (icd) would enhance the succinate production since the flux flowed from isocitrate to alpha-ketoglutarate and succinyl-CoA would result in the carbon lost. Since acetate is the main byproduct under aerobic condition [34], knocking out the pyruvate oxidase gene (poxB) and phosphate acetyltransferase gene (pta) is also expected to increase the succinate yield.


## 3.1. Elementary Mode Analysis of Glycerol Metabolism under Anaerobic Condition

Under anaerobic condition, the metabolic network model got 55 elementary flux modes. The relationship between the yields of products and biomass was shown in Figure2. The maximum molar yield of biomass under anaerobic condition is 0.187 mol/mol in which the respective yields of 1,2-propanediol, ethanol, and formate were 0.248 mol/mol, 0.495 mol/mol, and 0.57 mol/mol and no succinate, acetate, and lactate were produced. It was found that the cell growth was always associated with the production of 1,2-propandediol, ethanol, and formate. Therefore, the production of 1,2-propanediol, ethanol, and formate was necessary for the biomass synthesis during the glycerol metabolism. The biomass synthesis process consumes ATP and produces reducing equivalents (NADH). Both ATP and NAD need to be regenerated through the production of other byproducts (the biomass synthesis equation in Appendix A.2). For the glycerol metabolism in anaerobic condition, only the glycerol to 1,2-propanediol pathway can consume extra NADH (see (1)) and thus provide the mean to consume the reducing equivalents generated during the synthesis of biomass. The conversion of glycerol to ethanol and formate (see (2)), a redox-balanced pathway, fulfills energy requirements by generating ATP via substrate-level phosphorylation (Appendix A.4). The calculation results were consistent with the experiment observation that 1,2-propanediol and ethanol were growth-associated products [13, 14].
(1)glycerol+NADH→1,2-propanediol+NAD,(2)glycerol+ADP→ethanol+formate+ATP.Relationship between the yields of biomass and byproducts for the obtained elementary modes ofE. coli under anaerobic conditions. (a) Succinate, (b) Lactate, (c) Acetate, (d) 1,2-Propanediol, (e) Formate, and (f) Ethanol. The enclosed regions represent the possible solution space. The fluxes were normalized by glycerol uptake rate and expressed as mol/mol (glycerol).
(a)(b)(c)(d)(e)(f)The maximum succinate yield under anaerobic condition is 1.0 mol/mol when CO2 or carbonate salts are added as cosubstrates and the optimal flux distribution for succinate production was shown in Figure 3. In this case, there was no production of biomass and other byproducts. The key points for this mode were that the phosphorylation of DHA was only catalyzed by ATP-dependent DHA kinase and PEP was totally carboxylated into oxaloacetate by PEP carboxylase and the latter was further transferred into succinate. This required a very high activity of ATP-dependent DHA kinase and PEP carboxylase. However, it was reported that the PEP-dependent DHA kinase plays the main role in E. coli which dramatically reduced the yield of succinate. With single PEP-dependent DHA kinase function, there will be no succinate production [14]. Thus the PEP-dependent DHA kinase is the bottleneck of succinate production under anaerobic condition.Figure 3
The optimum flux distribution of glycerol metabolism for succinate production inE. coli under anaerobic conditions when only ATP-dependent DHA kinase plays a function.


## 3.2. Element Ary Mode Analysis of Glycerol Metabolism under Aerobic Condition

The aerobic metabolic network model got 259 elementary flux modes. The relationship between the yields of products and biomass was shown in Figure5. It was indicated that the flux distribution modes under aerobic condition were completely different from that of anaerobic condition. The maximum molar yield of biomass under aerobic condition was 0.725 mol/mol in which only CO2 was produced. The results suggested that aerobic condition was more favorable for biomass formation and the most effective way for biomass formation was through the TCA cycle. The cell growth process was no longer associated with the production of 1,2-propanedio, formate, and ethanol. Comparing with the anaerobic condition, the succinate production modes were more diversely distributed throughout the feasible solution space, and the corresponding yield of succinate at the same biomass yield was higher which indicated that the high potential of succinate production is associated with the cell growth by metabolic modification under aerobic condition.The maximum succinate yield under aerobic condition (an aerobic mode suggests oxygen consumption) is also 1.0 mol/mol which requires the CO2 or carbonate salts to be added as cosubstrates, and the optimal flux distribution for succinate production was shown in Figure 6. The optimal flux distribution modes were quite similar for the aerobic and anaerobic conditions. The key point for obtaining high succinate yield was considered as that PEP was totally carboxylated into oxaloacetate by PEP carboxylase and the latter was further transferred into succinate. This required a very high activity of PEP carboxylase.The network robustness and its sensitivity to perturbation were critical to the optimal metabolic pathway. The sensitivity of succinate yield to the flux ratios at the key branch nodes PEP and acetyl-CoA was considered in this work. PEP was consumed in reactions R40 and R17 which are catalyzed by PEP carboxylase and Pyruvate kinase, respectively; the flux ratio was denoted as(3)R(40,17)=R40(R40+R17).The influence of flux distribution at PEP node was shown in Figure7(a). The succinate yield increased from 0.5 mol/mol to 1.0 mol/mol when R(40,17) increased from 0 to 1.0. The results showed that increasing the flux distribution from PEP to oxaloacetate was beneficial for succinate production. For acetyl-CoA node, it was consumed in reactions R22, R44, R50, and R53 which are catalyzed by citrate synthase, aldehyde dehydrogenase, malate synthase, and phosphate acetyltransferase, respectively. The flux ratio was denoted as
(4)R(22,44,50,53)=R22(R22+R44+R50+R53).The influence of flux distribution at acetyl-CoA node was shown in Figure7(b). The succinate yield decreased from 1.0 mol/mol to 0.75 mol/mol when R(22,44,50,53) increased from 0 to 1.0. The results indicated that decreasing the flux distribution from acetyl to TCA cycle was beneficial for succinate production since the carbon would be lost during the TCA cycle.


## 3.3. Effect of Oxygen Consumption on the Production of Succinate and Biomass

The effect of oxygen consumption on succinate and biomass production was further investigated according to the elementary mode analysis; the relationship between the molar fractions of oxygen consumption and succinate or biomass under aerobic conditions was calculated and shown in Figure8. The theoretical succinate yield decreased as the oxygen consumption ratio increased (Figure 8(a)). The elementary modes of succinate yield distributed only on the left part of the solution space, which indicated that the higher consumption of oxygen was unfavorable for the succinate production. This is reasonable from the redox consideration. The higher consumption of oxygen results in more NAD which is critical for succinate production, to be utilized for the ATP production through oxidative phosphorylation. Especially as the molar oxygen consumption fraction of oxygen was more than 1.75 mol/mol, the elementary mode that produced succinate did not exist.The theoretical biomass yield increased firstly and then decreased when the oxygen consumption ratio increased (Figure8(b)). This is reasonable and consistent with the results of electron conservation because the rational increase of oxygen flux would be favorable for ATP synthesis through oxidative phosphorylation which is essential for biomass synthesis. However, high oxygen flux also results in less carbon source and reducing equivalents available for biomass because more carbon source would be oxidized to CO2. The maximum biomass was achieved when the molar fraction of oxygen was 0.65 mol/mol. Since higher cell concentration is beneficial for increase of the productivity, the optimal regulation strategy is controlling the molar fraction of oxygen consumption less than 0.65 mol/mol.


## 3.4. Rational Design to Improve the Succinate Production by Genetic Modifications

Comparing the results of elementary flux mode analysis above, although the maximum succinate yields were 1.0 mol/mol under both anaerobic and aerobic conditions, the aerobic condition seemed to be more favorable for succinate production in fact. The slow cell growth under anaerobic condition hindered the practical application of glycerol fermentation for succinate production [33]. What is more, the association of cell growth with the production of 1,2-propanediol, ethanol, and formate reduced the total yield of succinate.As discussed above, to improve the yield of succinate yield under anaerobic condition, the substitution of PEP-dependent DHA kinase into ATP-dependent DHA kinase and overexpressing the PEP carboxylase would be a prior consideration. An alternative choice is to express the heterogeneous pyruvate carboxylase inE. coli. The overexpression of pyruvate carboxylase could redistribute the flux of pyruvate into oxaloacetate for succinate production. The optimal flux distribution for succinate production in such case could also reach 1.0 mol/mol which was shown in Figure 4. It has been reported that expressing the ATP-dependent DHA kinase and pyruvate carboxylase could both increase the yield of succinate [14]. Another consideration to increase the succinate yield is reducing the byproducts production. Since the production of 1,2-propanediol, ethanol, and formate is necessary for biomass synthesis, the proper strategy is reducing the acetate and lactate production. Knockout of the phosphate acetyltransferase gene (pta) and lactate dehydrogenase gene (ldh) would be expected to increase the succinate production.Figure 4
The alternative optimum flux distribution of glycerol metabolism for succinate production inE. coli under anaerobic conditions when only PEP-dependent DHA kinase plays a function and heterogeneous PEP carboxylase (pyc) are introduced and overexpressed in E. coli.Relationship between the yields of biomass and byproducts for the obtained elementary modes ofE. coli under aerobic conditions. (a) Succinate, (b) Lactate, (c) Acetate, and (d) Ethanol. The enclosed regions represent the possible solution space. The fluxes were normalized by glycerol uptake rate and expressed as mol/mol (glycerol).
(a)(b)(c)(d)Figure 6
The optimum flux distribution of glycerol metabolism for succinate production inE. coli under aerobic conditions.Sensitivity of succinate yield to the relative fluxes at the (a) PEP node and (b) AcCoA node under aerobic conditions. The PEP node involves the catabolic reactions of R40 and R17 which are catalyzed by PEP carboxylase and pyruvate kinase, respectively. The AcCoA node involves the catabolic reactions of R22, R50, R44, and R53 which are catalyzed by citrate synthase, aldehyde dehydrogenase, malate synthase, and phosphate acetyltransferase, respectively.
(a)(b)Effect of oxygen consumption on the succinate production and biomass formation under aerobic conditions. The enclosed regions represent the possible solution space. The fluxes were normalized by glycerol uptake rate and expressed as mol/mol (glycerol).
(a)(b)For the case of aerobic condition, overexpressing the PEP carboxylase or expressing pyruvate carboxylase would be a prior consideration as discussed above. Knockdown of the isocitrate dehydrogenase gene (icd) would enhance the succinate production since the flux flowed from isocitrate to alpha-ketoglutarate and succinyl-CoA would result in the carbon lost. Since acetate is the main byproduct under aerobic condition [34], knocking out the pyruvate oxidase gene (poxB) and phosphate acetyltransferase gene (pta) is also expected to increase the succinate yield.


## 4. Conclusions

Nowadays, the conversion of byproduct glycerol has attracted more and more attention with the development of biodesiel industrial. The potentials of using glycerol for succinate production inE. coli under the anaerobic and aerobic conditions were compared by using elementary mode analysis in this work. The aerobic conditions seem to be more favorable for succinate production and the maximum succinate yield was 1.0 mol/mol. Although increase of the oxygen concentration would reduce the succinate yield, controlling the molar fraction of oxygen under 0.65 mol/mol would be beneficial for increasing the succinate productivity. According to the elementary mode analysis, the rational design was obtained for improving the succinate production by genetic modification under aerobic and anaerobic conditions, respectively. The results are considered useful for further investigation on the succinate metabolism of E. coli. The information also is beneficial for the efficient production of succinate from glycerol by E. coli.Table 1
Reactions and enzymes involved in Figure1.
ReactionsGenesEnzymesAerobic/anaerobic specificityReferencesR1glpFglycerol facilitator[1]R2gldAGlycerol dehydrogenase IIanaerobic[2]R3dhaKLMPTS-dependent Dihydroxyacetone kinaseanaerobic[2]R4glpKglycerol kinaseaerobic[1]R5glpDGlycerol-3-phosphate dehydrogenaseaerobic[1]R6sets of reactionsanaerobic[2]R7dhaKATP-dependent Dihydroxyacetone kinase[2]R10tpi12Triose-phosphate isomeraseKEGGR11fbaFructose-bisphosphate aldolaseKEGGR13fbpFructose-bisphosphataseKEGGR14gpiGlucose-6-phosphate isomeraseKEGGR15sets of reactionsKEGGR16sets of reactionsKEGGR17pykAFPyruvate kinaseKEGGR20pdhPyruvate dehydrogenaseaerobic[3]R21pflPyruvate formate-lyaseanaerobic[4]R22gltACitrate synthaseKEGGR23acnBAconitate hydrataseKEGGR24icdIsocitrate dehydrogenaseKEGGR25sucABOxoglutarate dehydrogenaseaerobic[5]R26sucCDSuccinate-CoA ligaseaerobic[5]R27sdhABCDSuccinate dehydrogenaseKEGGR28fumABCFumarate hydrataseKEGGR29mdhMalate dehydrogenaseKEGGR30sets of reactionsKEGGR31rpeRibulose-phosphate 3-epimeraseKEGGR32rpiARibose-5-phosphate isomeraseKEGGR33tktTransketolaseKEGGR34tktTransaldolaseKEGGR35tktTransketolaseKEGGR40ppcPEP carboxylaseKEGGR41pckPEP carboxykinaseKEGGR42malEmalic enzymeKEGGR43aceAIsocitrate lyaseKEGGR44aceBMalate synthaseKEGGR50adhEAldehyde dehydrogenaseKEGGR51ldhALactate dehydrogenaseKEGGR52poxBpyruvate oxidaseaerobic[6]R53pta,ackAPhosphate acetyltransferase, Acetate kinaseKEGGR60Biomass formation[7]R70ATP-hydrolysis[7]R71Transhydrogenase[7]R72respiratory chain 1aerobic[7]R80Membrane transport reactionR81Membrane transport reactionR82Membrane transport reactionR83Membrane transport reactionR84Membrane transport reactionR85Membrane transport reactionTable 2
Degree of reduction of considered substrates and products.
Glycerol14O2-4CO20Succinate14Ethanol12Lactate12Formate21,2-propanediol16Acetate8



---

*Source: 518743-2010-09-19.xml*