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 . 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 [ 24 ] . several recent studies also tried to utilize glycerol as a carbon source for the transformation of other valued products such as ethanol   and amino acids . 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 . 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 [ 810 ] . the dissimilation of glycerol in e. 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 ] . although the aerobic utilization of glycerol by e. 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 . since e. 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 [ 1618 ] . 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 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 , the production of polyhydroxybutanoate in yeast , 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 of e. 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 . the glycerol metabolic network of e. 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 a p / o ratio of 2 for nadh was assumed .the precursor demand for biomass formation was calculated according to the literature . the cell physiology of e. coli is strongly affected by oxygen from different levels such as the transcriptional regulatory which can not be represented simply in the metabolic network . for example , one component fumarate and nitrate reduction ( fnr ) protein is aerobic / anaerobic response regulator . 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 ( [ 4fe4s ] ) clusters in the dimeric fnr protein are converted to two [ 2fe2s ] clusters upon exposure to stoichiometric levels of oxygen . 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 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 . the tca cycle is broken at the alpha - ketoglutarate dehydrogenase step and the respiratory chain is assumed to be inactive . 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 detailed description of the model is listed in appendices a.1 and a.3 . in the present work , the elementary mode analysis was carried out for studying the aerobic and anaerobic metabolism of glycerol in e. coli by using metatool 5.1 . 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 ) . 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 ) . 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 . 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 ] . 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 . thus the pep - dependent dha kinase is the bottleneck of succinate production under anaerobic condition . 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 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 . 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 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 . r44 , r50 , and r53 which are catalyzed by citrate synthase , aldehyde dehydrogenase , malate synthase , and phosphate acetyltransferase , respectively . the influence of flux distribution at acetyl - coa node was shown in figure 7(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 . 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 figure 8 . 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 . 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 ( figure 8(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 . 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 . 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 . 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 . 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 .   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 , knocking out the pyruvate oxidase gene ( poxb ) and phosphate acetyltransferase gene ( pta ) is also expected to increase the succinate yield . 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 in e. 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 . enzrev ( reversible reactions ) 
							 r10r    r11r    r14r    r15r    r16r    r23r    r26r    r27r    r28r    r29r    r31r    r32r    r33r    r34r     r35r    r71r    r86r . enzirrev ( irreversible reactions ) 
							 r1    r2    r3     r6    r7    r13    r17    r21    r22    r24    r30    r40    r41    r42    r43    r44    r50    r51    r53    r60    r70    r80    r81    r82    r83    r84    r85 . metint ( internal metabolite declaration ) 
							 glycerol    dha    dhap    ga-3-p    propanediol    pg    pep    pyruvate    acetyl - coa    coash    oxaloacetate     citrate    isocitrate    a - ketoglutarate    succinate    fumarate    malate    glyoxylate    glucose-6-p    fructrose-6-p    fructrose-16-p    ribulose-5-p    xylulose-5-p    ribose-5-p    sed-7-p    erythrose-4-p    lactate    formate    acetate    ethanol    nad    nadh    atp    adp    nadp    nadph    co2 . metext ( external metabolite declaration ) 
							 glycerol_ext    ethanol_ext    acetate_ext    co2_ext     lactate_ext    succinate_ext     formate_ext biomass    propanediol_ext . glycerol specific metabolisms 
 r1:glycerol_ext = glycerol.r2:glycerol + nad = dha + nadh.r3:pep + dha = dhap + pyruvate.r6:dhap  +  2nadh = propanediol  +  nad.r7:dha  +  atp = dhap  +  adp . r10r : dhap = ga-3-p.r11r : dhap + ga-3-p = fructrose-16-p.r13:fructrose-16-p = fructrose-6-p.r14r : fructrose-6-p = glucose-6-p.r15r : ga-3-p + adp + nad = pg + atp + nadh.r16r : pg = pep.r17:pep + adp = pyr + atp . r21:pyr + coash = acetyl - coa + formate.r22:oxaloacetate + acetyl - coa = citrate + coash.r23r : citrate = isocitrate.r24:isocitrate + nadp = a - ketoglutarate + nadph + co2.r27r : succinate + nad = fumarate + nadh.r28r : fumarate = malate.r29r : malate + nad = oxaloacetate + nadh . pentose phosphate pathway 
 r30:glucose-6-p + 2nadp = ribulose-5-p + 2nadph + co2.r31r : ribulose-5-p = xylulose-5-p.r32r : ribulose-5-p = ribose_5_p.r33r : ribose-5-p + xylulose-5-p = sed-7-p + ga-3-p.r34r : ga-3-p + sed-7-p = erythrose-4-p + fructrose-6-p.r35r : erythrose-4-p + xylulose-5-p = ga-3-p + fructrose-6-p . anapleurotic reactions 
 r40:pep + co2 = oxaloacetate.r41:oxaloacetate + atp = pep + adp + co2.r42:malate + nadp = pyruvate + nadph + co2.r43:isocitrate = glyoxylate + succinate.r44:glyoxylate + acetyl - coa = malate + coash . 
  pep + co2 = oxaloacetate . oxaloacetate + atp = pep + adp + co2 . redox - associated reactions 
 r50:acetyl - coa + 2nadh = ethanol + 2nad + coash.r51:pyruvate + nadh = lactate + nad.r53:acetyl - coa + adp = acetate + coash + atp . r60:0.0206glucose-6-p + 0.0072fructrose-6-p + 0.0627ribose-5-p + 0.0361 erythrose-4-p + 0.0129ga-3-p + 0.1338pg + 0.0720pep + 0.2861pyruvate + 0.2930acetyl - coa + 0.1481 oxaloacetate + 0.1078 a - ketoglutarate + 1.6548 nadph + 1.7821atp + 0.3548 nad = 2.87 biomass + 1.6548 nadp + 0.2930 coash + 0.1678 co2 + 1.7821 adp + 0.3548 nadh . 0.0206glucose-6-p + 0.0072fructrose-6-p + 0.0627ribose-5-p + 0.0361 erythrose-4-p + 0.0129ga-3-p + 0.1338pg + 0.0720pep + 0.2861pyruvate + 0.2930acetyl - coa + 0.1481 oxaloacetate + 0.1078 a - ketoglutarate + 1.6548 nadph + 1.7821atp + 0.3548 nad = 2.87 biomass + 1.6548 nadp + 0.2930 coash + 0.1678 co2 + 1.7821 adp + 0.3548 nadh . oxidative phosphorylation / maintenance energy 
 r70:atp = adp.r71r : nadph + nad = nadh + nadp . membrane transport reactions 
 r80:lactate = lactate_ext.r81:formate = formate_ext.r82:acetate = acetate_ext.r83:ethanol = ethanol_ext.r84:succinate = succinate_ext.r85:propanediol = propanediol_ext.r86r : 
												 co2 = co2_ext . 
  lactate = lactate_ext . enzrev ( reversible reactions ) 
							 r10r    r11r    r14r    r15r    r16r    r23r    r26r    r27r    r28r    r29r    r31r    r32r    r33r    r34r     r35r    r71r    r86r . enzirrev ( irreversible reactions ) 
							 r1    r4    r5    r13    r17    r20    r22    r24    r25    r30    r40    r41    r42    r43    r44    r50    r51    r53    r60    r70    r72    r80    r81    r82    r83    r84    r87 . metint ( internal metabolite declaration ) 
							 glycerol    glycerol-3-p    dhap    ga-3-p    pg    pep    pyruvate    acetyl - coa    coash    oxaloacetate     citrate    isocitrate    a - ketoglutarate    succinyl - coa    succinate    fumarate    malate    glyoxylate    glucose-6-p    fructrose-6-p    fructrose-16-p    ribulose-5-p    xylulose-5-p    ribose-5-p    sed-7-p    erythrose-4-p    lactate    acetate    ethanol    nad    nadh    atp    adp    nadp    nadph    co2 o2 . metext ( external metabolite declaration ) 
							 glycerol_ext    ethanol_ext    acetate_ext    co2_ext     lactate_ext    succinate_ext     biomass    o2_ext . glycerol specific metabolisms 
 r1:glycerol_ext = glycerol.r4:glycerol + atp = glycerol-3-p + adp.r5:glycerol-3-p + nad = dhap + nadh . r10r : dhap = ga-3-p.r11r : dhap + ga-3-p = fructrose-16-p.r13:fructrose-16-p = fructrose-6-p.r14r : fructrose-6-p = glucose-6-p.r15r : ga-3-p + adp + nad = pg + atp + nadh.r16r : pg = pep.r17:pep + adp = pyruvate + atp . r20 : pyruvate + coash + nad = acetyl - coa + co2 + nadh.r22:oxaloacetate + acetyl - coa = citrate + coash.r23r : citrate = isocitrate.r24:isocitrate + nadp = a - ketoglutarate + nadph + co2.r25:a - ketoglutarate + nad + coash = succinyl - coa + nadh + co2.r26r : succinyl - coa + adp = succinate + atp + coash.r27r : succinate + nad = fumarate + nadh.r28r : fumarate = malate.r29r : malate + nad = oxaloacetate + nadh . a - ketoglutarate + nad + coash = succinyl - coa + nadh + co2 . pentose phosphate pathway 
 r30:glucose-6-p + 2nadp = ribulose-5-p + 2nadph + co2.r31r : ribulose-5-p = xylulose-5-p.r32r : ribulose-5-p = ribose_5_p.r33r : ribose-5-p + xylulose-5-p = sed-7-p + ga-3-p.r34r : ga-3-p + sed-7-p = erythrose-4-p + fructrose-6-p.r35r : erythrose-4-p + xylulose-5-p = ga-3-p + fructrose-6-p . anapleurotic reactions 
 r40:pep + co2 = oxaloacetate.r41:oxaloacetate + atp = pep + adp + co2.r42:malate + nadp = pyruvate + nadph + co2.r43:isocitrate = glyoxylate + succinate.r44:glyoxylate + acetyl - coa = malate + coash . 
  pep + co2 = oxaloacetate . oxaloacetate + atp = pep + adp + co2 . redox - associated reactions 
 r50:acetyl - coa + 2nadh = ethanol + 2nad + coash.r51:pyruvate + nadh = lactate + nad.r52:pyruvate = co2 + acetate.r53:acetyl - coa + adp = acetate + coash + atp . acetyl - coa + 2nadh = ethanol + 2nad + coash . pyruvate + nadh = lactate + nad . r60:0.0206glucose-6-p + 0.0072fructrose-6-p + 0.0627ribose-5-p + 0.0361 erythrose-4-p +    + 0.0129ga-3-p + 0.1338pg + 0.0720pep + 0.2861pyruvate + 0.2930acetyl - coa + 0.1481 oxaloacetate + 0.1078 a - ketoglutarate   + 1.6548 nadph + 1.7821atp + 0.3548 nad = 2.87 biomass + 1.6548 nadp + 0.2930 coash + 0.1678 co2 + 1.7821 adp + 0.3548 nadh . 0.0206glucose-6-p + 0.0072fructrose-6-p + 0.0627ribose-5-p + 0.0361 erythrose-4-p +    + 0.0129ga-3-p + 0.1338pg + 0.0720pep + 0.2861pyruvate + 0.2930acetyl - coa + 0.1481 oxaloacetate + 0.1078 a - ketoglutarate   + 1.6548 nadph + 1.7821atp + 0.3548 nad = 2.87 biomass + 1.6548 nadp + 0.2930 coash + 0.1678 co2 + 1.7821 adp + 0.3548 nadh . oxidative phosphorylation / maintenance energy : 
 r70:atp = adp.r71r : nadph + nad = nadh + nadp.r72:nadh + 2adp + 1/2o2 = nad + 2atp . membrane transport reactions 
 r80:lactate = lactate_ext.r82:acetate = acetate_ext.r83:ethanol = ethanol_ext.r84:succinate = succinate_ext.r86r : 
												 co2 = co2_ext.r87 : 
												 o2_ext = o2 .
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 by e. 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 .