Skip to main content

Fine tuning the glycolytic flux ratio of EP-bifido pathway for mevalonate production by enhancing glucose-6-phosphate dehydrogenase (Zwf) and CRISPRi suppressing 6-phosphofructose kinase (PfkA) in Escherichia coli



Natural glycolysis encounters the decarboxylation of glucose partial oxidation product pyruvate into acetyl-CoA, where one-third of the carbon is lost at CO2. We previously constructed a carbon saving pathway, EP-bifido pathway by combining Embden-Meyerhof-Parnas Pathway, Pentose Phosphate Pathway and “bifid shunt”, to generate high yield acetyl-CoA from glucose. However, the carbon conversion rate and reducing power of this pathway was not optimal, the flux ratio of EMP pathway and pentose phosphate pathway (PPP) needs to be precisely and dynamically adjusted to improve the production of mevalonate (MVA).


Here, we finely tuned the glycolytic flux ratio in two ways. First, we enhanced PPP flux for NADPH supply by replacing the promoter of zwf on the genome with a set of different strength promoters. Compared with the previous EP-bifido strains, the zwf-modified strains showed obvious differences in NADPH, NADH, and ATP synthesis levels. Among them, strain BP10BF accumulated 11.2 g/L of MVA after 72 h of fermentation and the molar conversion rate from glucose reached 62.2%. Second, pfkA was finely down-regulated by the clustered regularly interspaced short palindromic repeats interference (CRISPRi) system. The MVA yield of the regulated strain BiB1F was 8.53 g/L, and the conversion rate from glucose reached 68.7%.


This is the highest MVA conversion rate reported in shaken flask fermentation. The CRISPRi and promoter fine-tuning provided an effective strategy for metabolic flux redistribution in many metabolic pathways and promotes the chemicals production.


It is an important challenge in metabolic engineering to reasonably allocate metabolic flux to achieve high yields of target products [1]. Traditional metabolic engineering methods modify and optimize an organism for production of chemicals by decreasing flow through competing pathways and introducing heterogeneous production pathways. As such, metabolic rewiring designs are necessary to increase flux towards essential metabolites, for example, overexpressing native pathways [2, 3], inhibition of competing pathways [4], increasing Coenzyme A (CoA) availability [5], and construction of pyruvate dehydrogenase bypass. Many strategies have been applied to improve production of chemicals [6].

A new strategy is to decrease the generation of harmful byproducts such as CO2 or increase the reuse of byproducts by constructing artificial synthetic pathways. With the rapid development of synthetic biology and molecular biotechnology, scientists have made great efforts to maximize microbial chemical yields focusing on enhancing the efficiency of CO2 fixation and decreasing CO2 emission. Many unnatural pathways have been constructed, such as CETCH [7], MCG [8], NOG [9], MOG [10], and so on. These pathways provide a variety of new ideas to use CO2 or one-carbon chemicals as carbon sources, and rewire metabolic pathways [11, 12].

In natural glycolysis, a variety of carbon sources are metabolized through the Embden-Meyerhoff-Parnas (EMP) pathway, which synthesizes C3 (pyruvate) and C2 (acetyl-CoA) metabolites. Acetyl-CoA is a precursor of almost all biosynthesis and energy metabolism pathways. It is normally produced via pyruvate decarboxylation, in which one-third of the carbon is lost as CO2. To improve carbon conversion efficiency, Xu et al. constructed the “EP-bifido pathway” in E. coli by introducing fxpk (encoding bifunctional phosphoketolase) and fbp (encoding fructose-1,6-bisphosphatase) to break down fructose 6-phosphate (F6P) into the theoretical maximum amount of acetyl-CoA from glucose. The edd gene in the Entner-Doudoroff (ED) pathway and the key enzyme-encoding gene pfkA (phosphofructose kinase A) in the EMP pathway were knocked out to save more F6P to break down into C2 metabolites; attenuate the flow of pyruvate to acetyl-CoA; and shift carbon flux from the EMP pathway to the pentose phosphate pathway (PPP) which supplies more NADPH. The EP-bifido pathway achieved high carbon yield of mevalonate (MVA) with 64.3 mol% [13], which is the highest MVA yield that have been reported.

PPP is an important energy metabolism pathway in all organisms. Increase of the dehydrogenase reactions of the PPP is effective in increasing the yield of NADPH-dependent products [14,15,16,17,18]. 13C metabolic flux analysis revealed that a MVA producing strain with EP-bifido pathway provided insufficient NADPH from oxidative PPP. In the EP-bifido pathway, the theoretical optimal carbon split between the EMP pathway and PPP in the EP-bifido route for MVA production is 1:6 [13]; thus the resulting maximum carbon MVA theoretical conversion rate of the EP-bifido pathway is 86% (mol/mol). In our previously constructed EP-bifido pathway, the detected carbon split ratio was 0.43:0.57, in which the PPP split ratio was too low. This carbon flux distribution needs to be further optimized.

In addition, the disruption of pfkA in the EP-bifido strain caused severe growth defect of the cells. This significantly limited the application of this synthetic pathway. Many reports have demonstrated this limitation of pfkA deficiency [19]. Metabolic engineering is usually static by blocking a competing pathway or introducing heterogeneous pathways permanently and continuously [20,21,22]. Sometimes this has a detrimental effect, for example in the early growth period carbon sources would ideally be dedicated to building biomass. Implementing flexible regulation would be valuable in engineering projects by rebalancing synthetic pathways to respond to the growth phase or the buildup of precursor metabolites [23,24,25].

Therefore, the artificial pathway must be optimized to be robust and the reducing power should be balanced. To overcome these challenges, here, we designed and constructed fine-tuned EP-bifido Escherichia coli strains and enhanced PPP by replacing promoter of zwf on the genome to a series of different strength promoters, and suppressing the expression of pfkA using the clustered regularly interspaced short palindromic repeats interference (CRISPRi) system. After CRISPRi was applied, we obtained a titer of 8.53 g/L MVA and a yield of 68.7% (mol/mol). This is the highest MVA yield reported in shaken flask fermentation.


Enhancement of PPP flux by increasing the expression of zwf in different degrees

First, to further optimize the flux ratio of EP-bifido pathway, we aimed to enhance the expression of the first key gene of the PPP, zwf, by replacing its promoter with relative stronger promoters. We selected five constitutive promoters with different strengths from the Anderson promoter library. The theoretical strength of each promoter is shown in Table 2. We compared the actual expression strength of these five synthetic promoters with the original zwf promoter by placing a gfp gene downstream of the promoters and cultivated the engineered strains for fluorescence intensity detection. The fluorescence/OD600 of detected at 16 h is shown in Fig. 1a. The strength of the promoters was relatively consistent with that stated by the Anderson promoter library. The strength of promoters BBa-J23100 and BBa-J23104 was relatively strong, and BBa-J23100 was the strongest. The strength of the original (native) zwf promoter is between that of BBa-J23108 and BBa-J23114, and is relatively weak.

Fig. 1
figure 1

The fermentation results of EP-bifido strains with modified zwf promoter replaced by different strength synthetic promoters. a The characterization of the strength of a series of synthetic promoters and zwf promoter. The fluorescence/OD600 of strains was detected at 16 h. MVA titre and molar conversion rate (b), cell growth (c) and glucose consumption (d) during 72 h fermentation

To detect the effect of PPP enhancement on MVA production, plasmids pBSA (expressing three enzymes catalysis acetyl-CoA to mevalonate) and pFF (carrying fbp and fxpk gene) were transformed into the five zwf-enhanced strains and cultivated with the control strain BW-P/pBSA pFF (abbreviated to BW-P BF). Strain BW-P10 BF showed almost the same growth and glucose consumption as the others, while the conversion rate of MVA was far higher than that in the control strain due to less byproducts generation. Promoters BBa-J23100 and BBa-J23108 resulted in the highest yield of MVA, 64.3% (mol/mol) and 62.3% (mol/mol) respectively, although the strength of the promoters did not show a perfectly positive correlation with the MVA yield. This proved that enhancing expression of gene zwf was effective for increasing the PPP flux.

13C-Metabolic flux analysis of changes in central carbon metabolism flux and energy metabolism

Strains BW-P10 BF and BW-P08 BF and control strain BW-P BF were chosen for metabolic flux analysis. With the zwf promoter replaced, the normalized data showed that the carbon flux through the oxidative part of the PPP was significantly increased, and the carbon flux through the TCA cycle was decreased, which was consistent with our expectations (Fig. 2). More carbon flux moved towards the EP-bifido pathway (see Additional file 2). The two zwf-expression-enhanced strains showed a large difference in TCA cycle flux, which may explain the growth difference between these strains (Fig. 1c).

In addition, the ATP, NADPH, and NADH synthesis capacity and glucose consumption of the three strains were compared based on the 13C-MFA data. After the EP-bifido pathway and the MVA synthetic pathway were introduced, the NADPH content and yield of the strain were significantly improved. pfkA deficiency shunted carbon flux to the PPP and the expression level of zwf was increased. Comparison of the zwf-expression-enhanced strains showed that overexpression of zwf enhanced NADPH synthesis, and the NADPH level was positively correlated with the promoter strength. Taking strains BW25113, BW-P BF, and BW-P10 BF as examples, the introduction of the EP-bifido pathway and overexpression of zwf changed the main source of NADPH: The main NADPH generating pathway shifted from isocitrate dehydrogenase in the TCA cycle to glucose-6-phosphate dehydrogenase in the PPP. This further proved that we have redirected part of the carbon flux of the EMP pathway to the PPP.

In addition, the production of NADH also changed significantly, as shown in Fig. 2c, d. Through zwf enhancement, the total amount and the yield of NADH decreased significantly. The NADH was produced distinctly in strain BW-P10 BF compared with wild-type strain BW25113: In strain BW25113, five dehydrogenases were the main source of NADH [glyceraldehyde-3-phosphate dehydrogenase (GAPDH), pyruvate dehydrogenase (PDH), malate dehydrogenase, α-ketoglutarate dehydrogenase, and succinate dehydrogenase]; in strain BW-P10 BF, NADH was mainly formed via GAPDH and PDH.

Fig. 2
figure 2

NADH, NADPH and ATP level and its derivations. a NADPH generation. b NADPH derivations of BW25113, BW-P BF, BW-P08 BF. c NADH generation. d NADH derivations of BW25113 and BW-P10 BF. d ATP generation. e ATP derivations. G6PDH, Glucose-6-phoshpate dehydrogenase; ICDH, Isolate dehydrogenase; GAPDH, Glyceraldehyde-3-phosphate dehydrogenase; ME, Malate enzyme; MDH, Malate dehydrogenase; PDH, Pyruate dehydrogenase; KDH, α-Ketoglutarate dehydrogenase; SDH, Succinate dehydrogenase; IDH, Isolate dehydrogenase

The EP-bifido pathway and MVA synthetic pathway expression increased the ATP yield from glucose, but the total amount of ATP decreased (Fig. 2e). The pfkA deficiency greatly impaired the EMP pathway and thereby blocked the three steps of substrate level phosphorylation absorption (glycerate-1, 3-diphosphate, phosphoenolpyruvate kinase, and succinyl CoA synthetase). This can also be verified from the origin ratio of ATP (Fig. 2f). Enhancing zwf expression resulted in increased ATP production and yield in strain BW-P08 BF compared with BW-P BF.

Down regulation of EMP pathway flux by targeting pfkA using CRISPRi system

To further rationally use the carbon source, we tried to suppress pfkA in a time-controlled way through exogenous induction and inhibition using CRISPRi system (Fig. 3). The CRISPRi gene regulation system requires only two components, dCas9 protein and a gRNA, to achieve regulation of the transcription level of any gene in the genome. The degree of suppression of gene expression can be controlled by adjusting the binding position and expression amount of the gRNA. Thus CRISPRi is widely used in the field of metabolic engineering and had shown a relatively good inhibition effect [26, 27].

Fig. 3
figure 3

Schematic diagram of CRISPRi on pfkA suppression. a Schematic diagram of pfkA transform. b Abridged general view of metabolic flux change after CRISPRi regulation. The red arrow indicated decreased metabolic flux, the green arrow indicated increased metabolic flux

To avoid the growth inhibition that may be caused by dCas9 from the CRISPRi gene regulation system, we selected a relatively low strength promoter, BBa-J23134, to promote dcas9. In order to obtain a different repression range, three different sgRNAs targeting the promoter or coding region of pfkA were designed. sgRNA1 were designed on the promoter region of pfkA, sgRNA2 and sgRNA3 targeted the coding chain of pfkA, at the region of 100 bp and 200 bp downstream of the initial codon, which may cause different repression effect [28]. After dcas9 and the sgRNAs were incorporated into pFF and pBSA respectively, six CRISPRi-regulated strains were generated. The fermentation results showed that the introduction of CRISPRi significantly inhibited the growth of cells and the glucose consumption was also reduced compared with that of the control strain BW25113 zwf-23100 pFF pBSA (abbreviated to BBF). This may be caused by the toxicity or leaky expression of dCas9. The CRISPRi-regulated strains were induced at 12 h by adding IPTG. The fermentation results showed that the three inhibitory sites we selected had different inhibitory effects (Fig. 4). sgRNA1 showed a better inhibition effect. Although strain BW25113 pFF-dCas9 pBSA-sgRNA1 produced only 8.53 g/L MVA, its MVA conversion rate reached 68.7%, which exceeded the previous best MVA conversion rate. Four control strains were also constructed to confirm the effect of the CRISPRi system on cell growth. Strains BW25113 zwf-23100 pBSA-sgRNA1 pFF, BW25113 zwf-23100 pBSA pFF-dCas9, BW25113 zwf-23100 pBSA-sgRNA1-dCas9 pFF and BW25113 zwf-23100 pBSA pFF, were abbreviated to BB1F, BBFd, BB1dF, and BBF, respectively. The results proved that CRISPRi can restrain cell growth effectively in these engineered strains (Additional file 1: Figure S1). The inhibitory effect of sgRNA1 on_pfkA is suitable for enhancing MVA fermentation in the EP-bifido system (Fig. 4). In the CRISPRi-regulated strains, we hardly detected any acetic acid, ethanol, or other byproducts during the fermentation process, which was in line with our expectations. The timely inhibition of pfkA reduces the flux of glycolysis, so there was no excessive carbon source overflow.

Fig. 4
figure 4

Fermentation of CRISPRi-regulated EP-bifido strains. a Glucose consumption of CRISPRi-regulated EP-bifido strains. b OD600 of CRISPRi-regulated EP-bifido strains. c MVA accumulation of CRISPRi-regulated EP-bifido strains. d MVA titre and yield of CRISPRi-regulated EP-bifido strains


Global warming is mainly due to excess CO2 emission; it is urgently necessary to find sustainable solutions to address this issue. Moreover, this wasted carbon may have a major impact on the overall economy of biobased products derived from fermentable carbon sources. Scientists have explored the possibility of using microbial systems to optimize carbon conservation during metabolic processes. The pyruvate decarboxylation step of glycolysis releases CO2 into the environment, resulting in 33% loss of carbon yield; this carbon loss has now been challenged by many scientists. We previously constructed an EP-bifido pathway in E. coli to reduce CO2 emissions and successfully applied it to produce a series of acetyl-CoA-derived compounds such as PHB, MVA and fatty acids. However, in the EP-bifido pathway, metabolic flux distribution between EMP and PPP has great optimization potential for maximum generation of NADPH; in addition, the pfkA knockout severely blocked EMP pathway and limited the growth of the engineered strains.

Here, our fine-tuning strategy to improve NADPH availability was to enhance the expression of zwf by replacing its promoter; we used five promoters with different strengths. The engineered strains BW-P08 BF and BW-P10 BF produced higher MVA titers than the control strain BW-P BF, 9.12 and 11.2 g/L, respectively. The MVA production by these strains did not show an obvious positive correlation with the zwf promoter strength. Since the MVA yield did not represent the flux distribution between the EMP pathway and the PPP, 13C-MFA was performed to detect the metabolic flux distribution in strains BW-P08 BF and BW-P10 BF (which had high MVA titers) and the control strain BW-P BF. The flux ratio between the PPP and the EMP pathway in strains BW-P08 BF and BW-P10 BF was much higher than that in strain BW-P BF (Fig. 5), indicating an improved shunt to the PPP. Also, enhanced zwf expression increased the total amount and molar yield of NADPH (Fig. 2). The NADPH-generating pathway shifted from isocitrate dehydrogenase in the TCA cycle to glucose-6-phosphate dehydrogenase in the PPP, further proving that carbon flux was redirected from the EMP pathway to the PPP. The NADH level also showed a decreased TCA cycle activity. In terms of the ATP level, enhancement of zwf expression increased the ratio of substrate level phosphorylation (Fig. 2) and reduced the energy supply ratio of pyruvate kinase. Thus, we identified the metabolic flux distribution following the fine-tuning of central metabolic nodes, which helps us to understand the impact on metabolism.

Fig. 5
figure 5

Metabolic flux diagram of zwf-strengthened EP-bifido strains. The metabolic flux shown for strains from top to bottom are BW-P10 BF, BW-P08 BF, BW-P BF

In dynamic regulation of metabolic pathways, the CRISPRi system has recently been used to improve flux through different pathways [29, 30]. One benefit of using the CRISPRi system over promoter replacement methods is that it does not require genome editing of the target gene, which remains a challenge. The introduction of the CRISPRi system is achieved by adding an inducer at a certain time to start the CRISPRi system. Here, we used it to adjust the inhibition level of pfkA, so as to achieve timely adjustment of EMP pathway/PPP flux. To reduce the growth inhibition caused by pfkA knockout, considering that the strain itself already harbors two plasmids and the CRISPRi system, we integrated dcas9 and the sgRNA onto the two plasmids respectively. After introduction of the CRISPRi system, cell growth and sugar consumption of the engineered strains were significantly decreased. The introduction of dCas9 may also have an inhibitory effect on bacterial growth. The relevant limitation of the CRISPRi system is therefore the toxicity of dCas9 expression in certain hosts [31, 32] that would affect the growth of pathway-expressing cells that typically already suffer from growth defects. The three targeting sites played a role in fine-tuning of pfkA and sgRNA1 showed the best inhibition effect. The MVA yield of strain BW25113 pFF-dCas9 pBSA-sgRNA1 was only 8.53 g/L, its yield reached 68.7%, exceeding the previous best conversion rate. sgRNA2 and sgRNA3 targeted the coding chain of pfkA at the region of 100 bp and 200 bp downstream of the initial codon, while sgRNA1 were designed on the 100 bp upstream of the promoter region of pfkA. The design of sgRNA1 can block the binding of RNA polymerase due to the steric hindrance of dcas9; the design of sgRNA2 and sgRNA3 can prevent the elongation of RNA polymerase. The actual inhibitory effects of three sgRNAs on the transcription of pfkA were not revealed, however, the strategy of CRISPRi modulation can achieve a broadly gene regulating range through designing various CRISPRi target sites.

The more relevant limitation of the CRISPRi system is the toxicity of dCas9 expression that would affect cell growth. Except for the CRISPRi system, gene repression at the transcriptional level usually combines multiple known promoter systems to construct the inverter to down-regulate genes. Moreover, at post-transcriptional and post-translational level, tools to down-regulate mRNA levels include anti-sense RNA (asRNA) and RNA interference (RNAi). These control methods that function at the post-transcriptional and post-translational levels are useful alternatives. In the future, there must be more convenient methods to down-regulate genes autonomously and dynamically.


This study showed glycolytic flux ratio fine-tuning strategies applying in an artificial carbon saving pathway for efficient MVA production. The strategies presented in this work serve as a guide to metabolic engineering projects requiring acetyl-CoA as a metabolic precursor.

Materials and methods

Media and culture

For plasmid preparation, E. coli strains were cultured at 37 °C on a rotary shaker (220 rpm) in test tubes containing 5 mL Luria–Bertani (LB) medium. For MVA production, 50-mL shake flask cultures were started by 2% inoculation from the 5-mL LB culture. The 50-mL cultures contained M9 minimal medium with 0.2% yeast extract containing 20 g/Lglucose and shaken at 37 °C in a rotary shaker (120 rpm) for 48 h. Overnight cultures were shaken at 37 °C in a rotary shaker (220 rpm). Antibiotics were added as follows: ampicillin (Amp) 100 μg/mL, spectinomycin (Spc) 50 μg/mL and chloromycetin (Cm) 25 μg/mL. For promoter integration and replacement procedure, strains were cultivated in SOB medium.

LB medium contains (g/L): tryptone(10), yeast extract(5) and NaCl(10). M9 medium contains (g/L): Na2HPO4·12H2O (15.138), KH2PO4 (3), NaCl (0.5) and NH4Cl (1). SOB medium contains (g/L): tryptone (20), yeast extract (5) and NaCl (5).

Strains and plasmids

All E. coli strains and plasmids used are listed in Table 1. Strain BW-P was used as the starting strain for further genetic manipulation [13]. All primers used for molecular manipulations are listed in Additional file 1: Table S1. All promoter used for genetic manipulation are listed in Table 2.

Table 1 Strains and plasmids
Table 2 Relative strength of Promoters

Plasmid construction for MVA production

To replace the original tac promoter of pBSA plasmid, five promoters BBa-J 23119, BBa-J 23100, BBa-J 23102, BBa-J 23104, BBa-J 23118 was designed into primers to construct plasmids pBSA-23119, pBSA-23100, pBSA-23102, pBSA-23104, pBSA-23118. Two primers were designed in the opposite direction. Five PCR amplicons were obtained using the original plasmid as the template with primier pcr-23119-F/R, pcr-23100-F/R, pcr-23102-F/R, pcr-23104-F/R, pcr-23118-F/R. DpnI was added to the PCR system in 37 °C for 1 h to remove the methylated template. A mixture of 50 μL was transformed into competent cells using chemical transformation. Colony PCR was performed by picking monoclonal from resistance plate to eliminate false positives and template interference. Finally, the five plasmids were transformed into the BW-P 23100 strain with the plasmid pFF.

Construction of CRISPRi suppression system

To select the CRISPRi inhibition site, three different sgRNAs was designed by targeting pfkA promoter sequence, 100 bp downstream of the promoter sequence, and 200 bp downstream of the initiation codon. dcas9 and sgRNAs were assembled into pFF and pBSA plasmids respectively downstream of an IPTG-induced promoter. Primer dCas9FF-F/R, dCas9-F/R were used to amplify dCas9 sequence. Primer sgRNA-F1, sgRNA1-R and sgRNA-F2, sgRNA1-R were used to amplify sgRNA1 sequence by PCR. The two amplified sequences were overlapped from homology arms. Sequences sgRNA2 and sgRNA3 were amplified as above. All constructed plasmids were electro-transformed into E. coli strains. To accomplished the timely control of pfkA by CRISPRi system, cells was induced by 200 μM IPTG after 12 h of fermentation.

Promoter replacement on E. coli genome

Promoter for replacing the promoter of zwf gene were selected from the Anderson promoter library ( Promoter replacement primers homoarm-F and homoarm-cm-R were designed using homology arms at about 300–500 bp at both ends of the target gene promoter, and plasmid pKD3 or pKD4 was used as a template to obtain recombinant fragments with kan or Cm resistance by using PKD3-cm-F/R or PKD4-cm-F/R. Three amplified sequences were overlapped from homology arms and resistance tag. All five replacing sequences were amplified as above.

The Red homologous recombination method was employed for gene integration. The pTKRED complementary plasmid was transformed into the target strain. The electrotransfection were performed by growing BW-P in 5 mL LB medium at 30 °C and shaking at 220 rpm for 12 h. 5-mL shake flask cultures using SOB broth were started with a 1% inoculation from the overnight culture. Isopropyl-β-D-thiogalactopyranoside (IPTG) was added to a final concentration of 0.5 mM to induce λ-prophage (bet, gam, and exo) gene expression. Cells were then incubated at 30 °C and shaking at 220 rpm until reaching an OD600 of 0.5 to 0.6. Cell were collected (2 mL), pelleted, and washed three times with cold sterile water to make them electrocompetent. ssDNA mixture (1 μM) was added to electrocompetent cells and electroporated at 2.5 kV. Add 1 mL LB liquid medium and cultivate for 1 h. After centrifuged, the collected bacteria were plated onto plates containing 25ug/ml kan or 18ug/ml spc for overnight incubation at 37 °C. Transformed strains were selected by their kanR phenotype and were verified by PCR.

Measurement of extracellular metabolites

A spectrophotometer was used to measure the optical density at 600 nm (OD600) of the bacterial culture. For extracellular metabolite analysis, 1 mL of culture was centrifuged at 12,000g for 2 min. The supernatant was filtered through a 0.22-μm syringe filter for high-performance liquid chromatography analysis. Glucose, MVA, acetate, and pyruvate were measured on an ion exchange column (HPX-87H; Bio-Rad Labs) with a differential refractive index detector (Shimadzu RID-10A). A 0.5-mL/min mobile phase using a 5 mM H2SO4 solution was applied to the column. The column was operated at 65 °C.


To investigate if the carbon flux was really redistributed to the newly constructed EP-bifido pathway, 13C-MFA was performed using 100% 1-13C1 glucose as the feeding substrate was added to a concentration of 10 g/L. Cells at the exponential growth phase were harvested by centrifugation at 7000 g for 5 min at 4 °C. The cell pellet was then washed twice with chemical defined medium and hydrolyzed in 6 M HCl for 24 h at 120 °C (Schwender et al. 2006). The resulting proteinogenic acids were derivatized with N-(tert-butyldimethylsilyl)-N-methyl-trifluoroacetamide containing tert-butyldimethylchlorosilane in acetonitrile at 105 °C for 1 h, and then analyzed by a GC–MS [Agilent 7890 A GC and 5975 C Mass Selective Detector (Agilent Technologies, Santa Clara, USA)] equipped with a DB-1column (Agilent Technologies). The data obtained from GC–MS were corrected by reduction of the natural abundance ratio of C, H, O, N, and Si isotopes [30]. Metabolic fluxes were estimated by minimizing the residual sum of squares between experimentally measured and model predicted 13C-enrichment using 13C-Flux software obtained from Dr. Wiechert [33].


  1. Keasling JD. Manufacturing molecules through metabolic engineering. Science. 2010;330:1355–8.

    Article  CAS  Google Scholar 

  2. Xu P, Ranganathan S, Fowler ZL, Maranas CD, Koffas MA. Genome-scale metabolic network modeling results in minimal interventions that cooperatively force carbon flux towards malonyl-CoA. Metab Eng. 2011;13:578–87.

    Article  CAS  Google Scholar 

  3. Hirokawa Y, Kubo T, Soma Y, Saruta F, Hanai T. Enhancement of acetyl-CoA flux for photosynthetic chemical production by pyruvate dehydrogenase complex overexpression in Synechococcus elongatus PCC 7942. Metab Eng. 2020;57:23–30.

    Article  CAS  Google Scholar 

  4. Centeno-Leija S, Huerta-Beristain G, Giles-Gomez M, Bolivar F, Gosset G, Martinez A. Improving poly-3-hydroxybutyrate production in Escherichia coli by combining the increase in the NADPH pool and acetyl-CoA availability. Antonie Van Leeuwenhoek. 2014;105:687–96.

    Article  CAS  Google Scholar 

  5. Ogata Y, Chohnan S. Prokaryotic type III pantothenate kinase enhances coenzyme A biosynthesis in Escherichia coli. J Gen Appl Microbiol. 2015;61:266–9.

    Article  CAS  Google Scholar 

  6. Ku JT, Chen AY, Lan EI. Metabolic engineering design strategies for increasing Acetyl-CoA Flux. Metabolites. 2020;10:166.

    Article  CAS  Google Scholar 

  7. Schwander T, Schada von Borzyskowski L, Burgener S, Cortina NS, Erb TJ. A synthetic pathway for the fixation of carbon dioxide in vitro. Science. 2016;354:900–4.

    Article  CAS  Google Scholar 

  8. Yu H, Li X, Duchoud F, Chuang DS, Liao JC. Augmenting the Calvin-Benson-Bassham cycle by a synthetic malyl-CoA-glycerate carbon fixation pathway. Nat Commun. 2008;2018:9.

    Google Scholar 

  9. Bogorad IW, Lin TS, Liao JC. Synthetic non-oxidative glycolysis enables complete carbon conservation. Nature. 2013;502:693–7.

    Article  CAS  Google Scholar 

  10. Bar-Even A, Noor E, Lewis NE, Milo R. Design and analysis of synthetic carbon fixation pathways. Proc Natl Acad Sci USA. 2010;107:8889–94.

    Article  CAS  Google Scholar 

  11. Hu G, Li Y, Ye C, Liu L, Chen X. Engineering microorganisms for enhanced CO2 sequestration. Trends Biotechnol. 2019;37:532–47.

    Article  CAS  Google Scholar 

  12. Francois JM, Lachaux C, Morin N. Synthetic biology applied to carbon conservative and carbon dioxide recycling pathways. Front Bioeng Biotechnol. 2019;7:446.

    Article  Google Scholar 

  13. Wang Q, Xu J, Sun Z, Luan Y, Li Y, Wang J, Liang Q, Qi Q. Engineering an in vivo EP-bifido pathway in Escherichia coli for high-yield acetyl-CoA generation with low CO2 emission. Metab Eng. 2019;51:79–87.

    Article  CAS  Google Scholar 

  14. Siedler S, Lindner SN, Bringer S, Wendisch VF, Bott M. Reductive whole-cell biotransformation with Corynebacterium glutamicum: improvement of NADPH generation from glucose by a cyclized pentose phosphate pathway using pfkA and gapA deletion mutants. Appl Microbiol Biotechnol. 2013;97:143–52.

    Article  CAS  Google Scholar 

  15. Chin JW, Cirino PC. Improved NADPH supply for xylitol production by engineered Escherichia coli with glycolytic mutations. Biotechnol Prog. 2011;27:333–41.

    Article  CAS  Google Scholar 

  16. Chemler JA, Fowler ZL, McHugh KP, Koffas MA. Improving NADPH availability for natural product biosynthesis in Escherichia coli by metabolic engineering. Metab Eng. 2010;12:96–104.

    Article  CAS  Google Scholar 

  17. Sundara Sekar B, Seol E, Park S. Co-production of hydrogen and ethanol from glucose in Escherichia coli by activation of pentose-phosphate pathway through deletion of phosphoglucose isomerase (pgi) and overexpression of glucose-6-phosphate dehydrogenase (zwf) and 6-phosphogluconate dehydrogenase (gnd). Biotechnol Biofuels. 2017;10:85.

    Article  Google Scholar 

  18. Kwon D-H, Kim M-D, Lee T-H, Oh Y-J, Ryu Y-W, Seo J-H. Elevation of glucose 6-phosphate dehydrogenase activity increases xylitol production in recombinant Saccharomyces cerevisiae. J Mol Catal B Enzym. 2006;43:86–9.

    Article  CAS  Google Scholar 

  19. Hollinshead WD, Rodriguez S, Martin HG, Wang G, Baidoo EE, Sale KL, Keasling JD, Mukhopadhyay A, Tang YJ. Examining Escherichia coli glycolytic pathways, catabolite repression, and metabolite channeling using Deltapfk mutants. Biotechnol Biofuels. 2016;9:212.

    Article  Google Scholar 

  20. Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, et al. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011;6:1290–307.

    Article  CAS  Google Scholar 

  21. Zhang YX, Perry K, Vinci VA, Powell K, Stemmer WP, del Cardayre SB. Genome shuffling leads to rapid phenotypic improvement in bacteria. Nature. 2002;415:644–6.

    Article  CAS  Google Scholar 

  22. Price ND, Reed JL, Palsson BO. Genome-scale models of microbial cells: evaluating the consequences of constraints. Nat Rev Microbiol. 2004;2:886–97.

    Article  CAS  Google Scholar 

  23. Farmer WR, Liao JC. Improving lycopene production in Escherichia coli by engineering metabolic control. Nat Biotechnol. 2000;18:533–7.

    Article  CAS  Google Scholar 

  24. Zhang F, Carothers JM, Keasling JD. Design of a dynamic sensor-regulator system for production of chemicals and fuels derived from fatty acids. Nat Biotechnol. 2012;30:354–9.

    Article  CAS  Google Scholar 

  25. Xu P, Li L, Zhang F, Stephanopoulos G, Koffas M. Improving fatty acids production by engineering dynamic pathway regulation and metabolic control. Proc Natl Acad Sci USA. 2014;111:11299–304.

    Article  CAS  Google Scholar 

  26. Cho S, Shin J, Cho BK. Applications of CRISPR/Cas system to bacterial metabolic engineering. Int J Mol Sci. 2018;19:1089.

    Article  Google Scholar 

  27. Larson MH, Gilbert LA, Wang X, Lim WA, Weissman JS, Qi LS. CRISPR interference (CRISPRi) for sequence-specific control of gene expression. Nat Protoc. 2013;8:2180–96.

    Article  CAS  Google Scholar 

  28. Bikard D, Jiang W, Samai P, Hochschild A, Zhang F, Marraffini LA. Programmable repression and activation of bacterial gene expression using an engineered CRISPR-Cas system. Nucleic Acids Res. 2013;41:7429–37.

    Article  CAS  Google Scholar 

  29. Kim S, Hahn J-S. Efficient production of 2,3-butanediol in Saccharomyces cerevisiae by eliminating ethanol and glycerol production and redox rebalancing. Metab Eng. 2015;31:94–101.

    Article  CAS  Google Scholar 

  30. Gordon GC, Korosh TC, Cameron JC, Markley AL, Begemann MB, Pfleger BF. CRISPR interference as a titratable, trans-acting regulatory tool for metabolic engineering in the cyanobacterium Synechococcus sp. strain PCC 7002. Metab Eng. 2016;38:170–9.

    Article  CAS  Google Scholar 

  31. Rock JM, Hopkins FF, Chavez A, Diallo M, Chase MR, Gerrick ER, Pritchard JR, Church GM, Rubin EJ, Sassetti CM, et al. Programmable transcriptional repression in mycobacteria using an orthogonal CRISPR interference platform. Nat Microbiol. 2017;2:1–9.

    Article  Google Scholar 

  32. Nielsen AAK, Voigt CA. Multi-input CRISPR/Cas genetic circuits that interface host regulatory networks. Mol Syst Biol. 2014;10:763.

    Article  Google Scholar 

  33. Weitzel M, Noh K, Dalman T, Niedenfuhr S, Stute B, Wiechert W. 13CFLUX2–high-performance software suite for (13)C-metabolic flux analysis. Bioinformatics. 2013;29:143–5.

    Article  CAS  Google Scholar 

  34. Kuhlman TE, Cox EC. Site-specific chromosomal integration of large synthetic constructs. Nucleic Acids Res. 2010;38:e92.

    Article  Google Scholar 

  35. Datsenko KA, Wanner BL. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci USA. 2000;97:6640–5.

    Article  CAS  Google Scholar 

Download references


We thank James Allen, PhD, from Liwen Bianji, Edanz Editing China (, for editing the English text of a draft of this manuscript.


This work was supported by grants from the National Key R&D Program of China (2019YFA0904900), the National Natural Science Foundation of China (31730003, 31670047 and 31670077), Young Scholars Program of Shandong University, the State Key Laboratory of Microbial Technology Open Projects Fund (M2019-11) and Guangdong Provincial Natural Science Foundation of China (no. 2018A030307013).

Author information

Authors and Affiliations



WQ and QQS designed the work. LY, XH and XY performed the experiments. SZJ and LY analyzed the 13C-MFA data. QQS encouraged this project. WQ and LY wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Zhijie Sun or Qian Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Fermentation of a series of CRISPRi-control strains. Table S1 Primers used in this study.

Additional file 2.

MFA simulated result.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, Y., Xian, H., Xu, Y. et al. Fine tuning the glycolytic flux ratio of EP-bifido pathway for mevalonate production by enhancing glucose-6-phosphate dehydrogenase (Zwf) and CRISPRi suppressing 6-phosphofructose kinase (PfkA) in Escherichia coli. Microb Cell Fact 20, 32 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: