- Open Access
Metabolic flux responses to genetic modification for shikimic acid production by Bacillus subtilis strains
© Liu et al.; licensee BioMed Central Ltd. 2014
- Received: 30 November 2013
- Accepted: 21 February 2014
- Published: 14 March 2014
Shikimic acid (SA) is a key chiral starting molecule for the synthesis of the neuramidase inhibitor GS4104 against viral influenza. Microbial production of SA has been extensively investigated in Escherichia coli, and to a less extent in Bacillus subtilis. However, metabolic flux of the high SA-producing strains has not been explored. In this study, we constructed with genetic manipulation and further determined metabolic flux with 13C-labeling test of high SA-producing B. subtilis strains.
B. subtilis 1A474 had a mutation in SA kinase gene (aroI) and accumulated 1.5 g/L of SA. Overexpression of plasmid-encoded aroA, aroB, aroC or aroD in B. subtilis revealed that aroD had the most significantly positive effects on SA production. Simultaneous overexpression of genes for 3-deoxy-D-arabinoheptulosonate-7-phosphate synthase (aroA) and SA dehydrogenase (aroD) in B. subtilis BSSA/pSAAroA/pDGSAAroD resulted in SA production of 3.2 g/L. 13C-Metabolic flux assay (MFA) on the two strains BSSA/pHCMC04/pDG148-stu and BSSA/pSAAroA/pDGSAAroD indicated the carbon flux from glucose to SA increased to 4.6% in BSSA/pSAAroA/pDGSAAroD from 1.9% in strain BSSA/pHCMC04/pDG148-stu. The carbon flux through tricarboxylic acid cycle significantly reduced, while responses of the pentose phosphate pathway and the glycolysis to high SA production were rather weak, in the strain BSSA/pSAAroA/pDGSAAroD. Based on the results from MFA, two potential targets for further optimization of SA production were identified. Experiments on genetic deletion of phosphoenoylpyruvate kinase gene confirmed its positive influence on SA production, while the overexpression of the transketolase gene did not lead to increase in SA production.
Of the genes involved in shikimate pathway in B. subtilis, aroD exerted most significant influence on SA accumulation. Overexpression of plasmid-encoded aroA and aroD doubled SA production than its parent strain. MFA revealed metabolic flux redistribution among phosphate pentose pathway, glycolysis, TCA cycle in the low and high SA-producing B. subtilis strains. The high SA producing strain BSSA/pSAAroA/pDGSAAroD had increased carbon flux into shikimate pathway and reduced flux into TCA cycle.
- Shikimic acid production
- Shikimate pathway
- Bacillus subtilis
- Metabolic flux assay (MFA)
Shikimic acid (3,4,5-trihydroxy-1-cyclohexene-1-carboxylic acid, SA), is a key chiral starting material for the synthesis of the antiviral neuraminidase inhibitor GS4104 [1, 2]. It is a key metabolic intermediate of the shikimate pathway for biosynthesis of aromatic amino acids (L-Phe, L-Trp, and L-Tyr) and many alkaloids in plants and microorganisms [3–5]. As a commercial product, SA has been extracted from the fruits of the Illicium plant. However, microbial fermentation as an alternative process for SA production has attracted more and more interests.
In this study, we investigated the metabolic fluxes with 13C-labeled glucose in two genetically modified, SA-producing strains of B. subtilis, namely BSSA/pHCMC04/pDG148-stu and BSSA/pSAAroA/pDGSAAroD. Our results indicated that significant metabolic flux changes occurred for the shikimate pathway, tricarboxylic acid (TCA) cycle, and reactions involving PEP and E4P generation and consumption. Based on the 13C-labeling metabolic flux assay (MFA), we further created mutants for SA production. The highest SA-producing B. subtilis strain produced 3.46 g/L of SA during batch cultivation in flasks.
Strain construction and SA productivity
Production of shikimic acid by various recombinant strains of Bacillus subtilis
1.50 ± 0.22
1.81 ± 0.18
Parent strain carrying pHCMC04
1.42 ± 0.20
Overexpression of aroA
1.74 ± 0.23
Overexpression of aroB
1.77 ± 0.15
Overexpression of aroC
2.30 ± 0.17
Overexpression of aroD
3.20 ± 0.07
Overexpression of aroD plus aroA
2.99 ± 0.05
Overexpression of aroD plus aroB
2.90 ± 0.03
Overexpression of aroD plus aroC
2.91 ± 0.04
Overexpression of aro D on two plasmids
3.11 ± 0.03
Overexpression of tkt in BSSA/pSAAroA/pDGSAAroD
3.46 ± 0.04
Knockout of pyk in BSSA/pSAAroA/pDGSAAroD
Strain 1A474 was used to construct high SA producers. For this purpose, we adopted a strategy to overexpress aro genes involving shikimate pathway in strain 1A474. The gene aroA, aroB, aroC, and aroD were cloned from the genome of strain 168 and were individually overexpressed in strain 1A474. SA production by the resulting recombinant strains was determined and was compared to their parent strains. Results showed that overexpression of plasmid-encoded aroA (strain BSSA/pSAAroA), aroB (strain BSSA/pSAAroB), or aroC (strain BSSA/pSAAroC), exhibited no significant increase in SA production, however, the overexpression of plasmid-encoded aroD in strain BSSA/pSAAroD resulted in higher production of SA (2.3 g/L, an increase of 53%), compared to its parent strain (1.5 g/L) (Table 1).
In order to improve SA production further, we co-overexpressed aroA, aroB, or aroC with aroD in B. subtilis. Results showed that co-overexpression of aroA (strain BSSA/pSAAroA/pDGSAAroD), aroB (strain BSSA/pSAAroB/pDGSAAroD), or aroC (strain BSSA/pSAAroC/pDGSAAroD) with aroD indeed increased SA production compared to strain BSSA/pHCMC04/ pDGSAAroD and parent strain 1A474 (Table 1). We also found that double overexpression of aroD on two plasmids further increased SA production compared to single overexpression of aroD (Table 1). Strain BSSA/pSAAroA/pDGSAAroD that co-overexpressed aroA and aroD was the best SA producer, with SA production of 3.2 g/L during batch cultivation fermentation. Compared to the parent strain 1A474, SA production was doubled in strain BSSA/pSAAroA/pDGSAAroD (Table 1).
Phenotypic characterization of growth, glucose consumption and metabolite production with strains BSSA474a and BSSA47407
Determination of growth (OD 600 ), specific cell growth rate (μ), specific glucose consumption rate (q glu ), specific shikimic acid production rate (q sa ), and specific 3-dehydroshikimic acid production rate (q dhs ) of B. subtilis strains BSSA/pHCMC04/ pDG148-stu and BSSA/pSAAroA/pDGSAAroD
qglu(mmol g-1 h-1)
qsa(mmol g-1 h-1)
qdhs(mmol g-1 h-1)
0.434 ± 0.023
5.199 ± 0.120
0.101 ± 0.011
0.132 ± 0.013
0.332 ± 0.015
4.666 ± 0.132
0.159 ± 0.014
0.048 ± 0.007
Metabolic flux responses to high production of SA in strain BSSA47407
Metabolic flux analysis by GC/LC-MS
Overall metabolic fluxes and responses of shikimate pathway
We observed a significant flux redistribution of several metabolic reactions in responses to the overexpression of plasmid-encoded aro genes in the high SA-producing strain BSSA/pSAAroA/pDGSAAroD. The overall conversion of glucose into shikimate pathway increased to 6.8% in strain BSSA/pSAAroA/pDGSAAroD from 4.4% in strain BSSA/pHCMC04/pDG148-stu. Consequently, the metabolic flux to SA increased from 1.9% in strain BSSA/pHCMC04/pDG148-stu to 4.6% in strain BSSA/pSAAroA/pDGSAAroD. We observed that the production rate of 3-dehydroshikimic acid decreased to 0.048 mmol g-1 h-1 in strain BSSA/pSAAroA/pDGSAAroD from 0.132 mmol g-1 h-1 in the parent strain BSSA/pHCMC04/pDG148-stu (Table 2). Correspondingly, the metabolic flux from E4P and PEP into DHS decreased from 2.5% in the parent strain BSSA/pSAAroA/pDGSAAroD to 2.2% in strain BSSA/pSAAroA/pDGSAAroD. The reversible reaction converting DHS to SA and being catalyzed by SA dehydrogenase (AroD) increased its flux by 2.4 folds (from 1.9% in BSSA/pHCMC04/pDG148-stu to 4.6% in BSSA/pSAAroA/pDGSAAroD) in strain BSSA47407 (Figure 3).
Flux shifts of TCA cycle, EMP and PPP pathways
As shown in Figure 3, the flux through the TCA cycle was 81.6% relative to the glucose uptake rate in BSSA/pHCMC04/pDG148-stu, while this flux decreased to 69.7% in BSSA/pSAAroA/pDGSAAroD. The overall flux of TCA cycle reduced significantly in the high SA-producing strain BSSA/pSAAroA/pDGSAAroD, suggesting a weak TCA cycle during high production of SA with B. subtilis. In contrast, fluxes of the EMP and pentose phosphate pathways did not shift significantly, except the reactions involving E4P generation and consumption. In the PPP pathway, transketolase, encoded by tkt, catalyzes reversible reaction (see Additional file 1: Appendix S1) and is that to E4P flux. We did not observe significant changes of flux driven by the transaldolase (v 6). However, the net flux of reaction encoded by tkt decreased from 4.6% in the parent strain BSSA/pHCMC04/pDG148-stu to 1.3% in the high SA-producing strain BSSA/pSAAroA/pDGSAAroD, suggesting that this reaction was severely disturbed and might be interesting to see how the SA synthesis would be influenced if tkt were overexpressed.
Supply of PEP for SA synthesis
PEP is one of the initial substrates for shikimate synthesis. In almost all known microorganisms, the DAHP synthase catalyzes the condensation of PEP and E4P, which is the first reaction of the shikimate pathway. It is reasonable to consider any changes of PEP generation and consumption in the high SA-producing BSSA/pSAAroA/pDGSAAroD. Two reactions were considered in this study for generation of PEP: 1) The enolase (NP_391270.1, EC:188.8.131.52) catalyzes glycerate-2-phosphate into PEP (v 8, see Additional file 1: Appendix S1). 2) The PEP-forming ATP:oxaloacetate carboxylyase (also named PEP carboxykinase, PckA, NP_390934.2, EC:184.108.40.206) catalyzes the conversion of oxaloacetate into PEP (v 10, see Additional file 1: Appendix S1). As shown in Figure 3, flux into reaction v 8 decreased from 160.3% in BSSA/pHCMC04/pDG148-stu to 150.2% in BSSA/pSAAroA/pDGSAAroD, and flux into reaction v 10 reduced from 44.3% in BSSA/pHCMC04/pDG148-stu to 21.4% in BSSA/pSAAroA/pDGSAAroD. Meanwhile, we also observed that the conversion of PEP into pyruvate, which is encoded by pyk, also decreased, from 194.2% in BSSA/pHCMC04/pDG148-stu to 159.0% in BSSA/pSAAroA/pDGSAAroD. Taking all these flux information, we considered that the supply of PEP might in shortage in the high SA-producing BSSA/pSAAroA/pDGSAAroD. In addition, we observed that the flux of reaction (v 9) in the SA-producing BSSA/pSAAroA/pDGSAAroD was still high (159%). For the purpose to increase supply of PEP to shikimate synthesis, knockout of this pyk might be a choice to increase SA production.
Validation of MFA results: effects of tkt overexpression and pyk deletion on SA production
Based on the previous results, we selected two targets for further genetic manipulation and to improve the SA production in B. subtilis. For the purpose to increase carbon flux to PEP, efforts were made to overexpress tkt in BSSA/pSAAroA/pDGSAAroD. However, the resulting strain BSSA/Ωtkt::pSAtkt did not show significant changes in SA production (Table 1). Further efforts were make to knock out the pyruvate kinase thus would diminish the flux of PEP into pyruvate. The resulting strain BSSA/Ωpyk::pSApyk produced 3.46 g/L of SA, an increase of 8% compared to its parent strain BSSA/pSAAroA/pDGSAAroD.
In this study, we investigated the influence of overexpression of aro genes on the accumulation of SA, and we constructed high SA producers by overexpression of aroA and aroD in B. subtilis. Our results demonstrated that overexpression of plasmid-encoded aroD increased the accumulation of SA, suggesting that AroD-driven reduction of 3-dehydroshikimate into SA was rate-limiting step for SA accumulation in B. subtilis. This is different from that of E. coli, of which the condensation of PEP and E4P into DAHP is the rate-limiting reaction and this reaction is subjected to extensive feedback regulations [19–21].
13C-MFA revealed carbon flux redistributions into shikimate pathway and TCA cycle, and reactions involving generation and consumption of PEP and E4P in genetically modified B. subtilis for SA production. Based on the MFA results in this study and previous investigation on E. coli[7, 8, 13], tkt and pyk were selected to be targets for genetically modification and for improvement of SA production with B. subtilis. Previous investigations showed that transformation with plasmid carrying tktA that encodes transketolase I increased the availability of E4P, and therefore increased SA and aromatic amino acid productions in E. coli. However, this study showed that overexpression of tkt did not increased SA production in B. subtilis, a remarkable difference from that in E. coli. Since transketolase catalyzes a reversible reaction (see Additional file 1: Appendix S1), i.e., it generates E4P in one direction but consumes E4P in the other direction, we deduced that the tkt-encoded transketolase in B. subtilis catalyzed reaction more favored to direction of fructose 6-phosphate and glyceraldehyde 3-phosphate, and that the tktI-encoded transketolase in E. coli more favored the formation of E4P and xylulose 5-phosphate. Confirmation of this deduction needs further biochemical evidences on catalytic properties of the transketolases from B. subtilis and E. coli. These results suggest that MFA data need to be carefully interpreted and conclusions from MFA need to be experimentally confirmed, and application of MFA data from one bacterial species may be different from that of another bacterial species.
Based on MFA results as well as previous investigations in E. coli, a second target was spotted for knockout of pyk, so that consumption of PEP by this Pyk-driven reaction would be vanished. Pyruvate kinase, encoded by pyk, is an enzyme involved in glycolysis and catalyzes the consumption of one molecule of PEP, yielding one molecule of pyruvate and one molecule of ATP. Previous investigation revealed that PEP increased its intracellular pools when pyk was deleted . We observed that SA production increased by 8% in B. subtilis, from 3.2 g/L in the control strain BSSA/pSAAroA/pDGSAAroD to 3.46 g/L in strain BSSA/Ωpyk::pSApyk. This result is again different from that observed in E. coli pyk mutant, which did not showed obvious increase of intracellular PEP level . Careful exploration of metabolic networks in B. subtilis and E. coli suggested that they have very different reactions at the PEP-pyruvate-oxaloacetate node, particularly the reactions for PEP and oxaloacetate generation/consumption . As demonstrated in previous studies, the pyk mutant of B. subtilis maintained high intracellular PEP levels  and produced less acetic acid than wild type did. This invoked the enthusiasm on exploitation of B. subtilis pyk mutants for production of various aromatic compounds. For example, B. subtilis strain BSTZ0402, an inducible pyk mutant, had been successfully constructed and yielded three folds higher folic acids compared to the parent strain B. subtilis 168 . Additional works on mining the 13C-MFA data and exploitation of metabolic network from B. subtilis for improving SA productivity are in progress.
Based on the data obtained in this study, it was concluded that overexpression of plasmid-encoded aroD exerted most significant influence on SA accumulation in B. subtilis. Co-everexpression of plasmid-encoded aroA and aroD doubled SA production than its parent strain BSSA/pHCMC04. So far as we know, this is the first time to quantify the flux flow in SA producing B. subtilis strains. MFA revealed metabolic flux redistribution among phosphate pentose pathway, glycolysis, TCA cycle in the low and high SA-producing B. subtilis strains. The high SA producing strain BSSA/pSAAroA/pDGSAAroD had increased carbon flux into shikimate pathway and reduced flux into TCA cycle. Based on the MFA results, tkt and pyk were targeted for improving SA production by B. subtilis. Overexpression of plasmid-encoded tkt did not increase SA production, and the deletion of pyk increased SA production by 8% in B. subtilis. The results showed that B. subtilis and E. coli had very different responses to genetic manipulations of the aro genes involving shikimate pathway and the tkt and pyk involving PEP/E4P generation and consumption.
Microorganisms, plasmids, cultivation, and determination of cell growth
Bacillus subtilis strains, plasmids, and primers used in this study
1A474 carrying pHCMC04
1A474 carrying pSAAroA
1A474 carrying pSAAroB
1A474 carrying pSAAroC
1A474 carrying pSAAroD
1A474 carrying pHCMC04 and pDG148-stu
1A474 carrying pHCMC04 and pDGSAAroD
1A474 carrying pSAAroA and pDGSAAroD
1A474 carrying pSAAroB and pDGSAAroD
1A474 carrying pSAAroC and pDGSAAroD
1A474 carrying pSAAroD and pDGSAAroD
BSSA/pSAAroA/pDGSAAroD harboring Ωtkt::pSATKT
BSSA/pSAAroA/pDGSAAroD harboring Ωpyk::pSAPYK
bla P xyl A xylR cat
neo ble Pspac lacI bla
bla Pspac lacZ lacI ermC
pHCMC04 derivative, aroA expressed from P xyl A
pHCMC04 derivative, aroB expressed from P xyl A
pHCMC04 derivative, aroC expressed from P xyl A
pHCMC04 derivative, aroD expressed from P xyl A
pDG148-stu derivative, aroD expressed from Pspac
bla Pspac-tkt-lacZ lacI ermC
bla Pspac-pyk’-lacZ lacI ermC
Cloning of aroA into pHCMC04
Cloning of aroB into pHCMC04
Cloning of aroC into pHCMC04
Cloningof aroD into pHCMC04
Cloning of aroD into pDG148-stu
Cloning of tkt into pMUTIN4
Confirming the integration of pSATKT
Cloning 5’ end of pyk into pMUTIN4
Confirming the integration of pSAPYK
Cell growth was monitored by optical density measurements at 600 nm (OD600) using a UV/visible spectrophotometer. Cellular dry weights were determined with duplicate samples of culture broth, by precipitation of cells and washed with distilled water and lyophilization.
Determination of SA and 3-dehydroshikimic acid by high performance liquid chromatography (HPLC)
B. subtilis was cultivated in LB broth at 37°C and at 200 rpm, overnight, and cells from culture were used as inoculum. SA production was determined in triplicates after 90 h cultivation at 37°C and 200 rpm in 250-mL flasks that contained 50 mL of fermentation broth. The fermentation broth was modified from Iomantas et al.  by deletion of maltose, and had the following components (g/L): K2HPO4 · 3H2O (18.3); KH2PO4 (6); Urea (6); sucrose (80); MgSO4 · 7H2O (1); FeSO4 (0.01); MnSO4 (0.01); Yeast Extract (4). For induction of the target gene overexpression in recombinant B. subtilis strains, IPTG at final concentration of 0.05 mM and xylose at final concentration of 5 g/L were added into the fermentation broth.
SA and 3-dehydroskimic acid accumulation in culture supernatant was determined with an HPLC system (Agilent 1200 series) equipped with a SB C18 column (4.6 mm × 250 mm × 5 μm). Five ml samples of fermentation broth were centrifuged at 13,000 rpm for 5 min on Microfuge® 18 centrifuge (BECKMAN COULTERTM) to remove cells, and the supernatant was applied for HPLC analysis. The HPLC was run with mixture of solution A (phosphoric acid in water, pH 2.5) and solution B (methanol) as eluant and was operated at a flow rate of 0.35 mL/min. SA and 3-dehydroshikimic acid were detected at 215 nm with adiode array detector. The following gradient was used: at 0–7.5 min, 95% of solution A and 5% of solution B; at 7.5-8.0 min, 100% of solution B; at 8.0-15 min, 100% of solution B; 15.0-15.5 min, 95% of solution A and 5% of solution B; 15.5-22.5 min, 95% of solution A and 5% of solution B. Standard SA and 3-hydroskimic acid were eluted at 5.411 and 6.241 min, respectively, under these conditions.
Determination of glucose concentrations
Glucose concentrations were enzymatically determined with glucose oxidase kit by using a bioanalyzer, according to producer’s description (SBA-40C, Shandong Academy of Sciences).
DNA extraction, amplification, plasmid construction and genetic transformation
Plasmid and chromosomal DNAs were isolated using the Plasmid Minspin HP Kit and the SiMax™ Genome DNA Kit, respectively. Purification of DNA fragments from agarose gels was done with the TIANgel Midi Purification Kit. Restriction enzymes, ligases and other DNA-manipulating enzymes were used according to their manufacturer’s instructions.
The aro genes involving shikimate pathway, i.e., aroA (Genbank accession number, NP_390853.1), aroB (NP_390151.1), aroC (NP_390189.1), and aroD (NP_390444.2) were PCR amplified from genomic DNA of B. subtilis strain 168 using primers listed in Table 3. The PCR products were digested with Spe I and Xma I and then cloned into the same digested pHCMC04 . All plasmids (Table 3) constructed in this study were verified by DNA sequencing. Transformation of chemically competent cells of E. coli DH5 and of B. subtilis were carried out according to the manufacturer’s manual or by the method of Spizizen . In order to simultaneous overexpression of aroA and aroD in B. subtilis, aroD was cloned into the pDG148-Stu vector , and the resulted plasmid was named pDGSAAroD (Table 3). Expression of the cloned aro genes in B. subtilis was induced by either 5 g/L of xylose for pHCMC04 derivatives or 0.05 mM IPTG for pDG148-Stu derivative.
Cloning and overexpression of tkt (NP_389672.10) in B. subtilis were performed with the plasmid pMUTIN4 . The tkt with its ribosome binding site (RBS) was amplified by using primers listed in Table 3 and was cloned into the Bam HI-Xma I treated backbone of pMUTIN4. Transformation of E. coli DH5α and B. subtilis BSSA/pSAAroA/pDGSAAroD and selection of transformants were performed as previously described. The B. subtilis transformants were selected on LB plates containing 15 ug/mL erythromycin, 5 ug/mL chloramphenicol, 30 ug/mL neomycin and then analyzed by PCR to ensure the integration of a single copy of the plasmid into the target gene on the chromosome.
Construction of an inducible pyk mutant
The plasmid pMUTIN4 was also used to produce knockout mutations . The fragment of 5’ end of pyk (NP_390796.1) with its RBS was amplified from Bacillus subtilis 168 chromosomal DNA using primers listed in Table 3, and was cloned into the Bam HI-Xma I treated backbone of pMUTIN4. Transformation of E. coli DH5α and B. subtilis BSSA/pSAAroA/pDGSAAroD were performed as previously described. The B. subtilis transformants were selected on LB plates containing 15 μg/mL erythromycin, 5 μg/mL chloramphenicol, 30 μg/mL neomycin and 1 mM IPTG and then analyzed by PCR to ensure that the pyk gene had been knocked out.
Growth and cultivation of B. subtilis with 13C-labeled glucose
The chemically defined M9 minimal medium  was applied for the 13C labeling experiment. Glucose was added to M9 minimal medium as either entirely the [1-13C]-labeled isotope isomer (99%; Cambridge Isotope Laboratories Inc., Andover, MA) or as a mixture of 20% (w/w) [U-13C] (99%; Cambridge Isotope Laboratories Inc.) and 80% (w/w) unlabeled glucose. B. subtilis strains pre-cultivated overnight in 5 mL of M9 minimal medium supplemented with 5 g/L unlabeled glucose, and the precultures (2 mL) were used as inoculum for 13C-labeling experiments. The 13C-labeling experiments were performed in a 250-mL flask that contained 50 mL of M9 minimal medium, which was supplemented with either 20:80 (w/w) mixture of [U-13C] glucose:unlabeled glucose or 100% [1-13C] glucose. Temperature and the shaking rate of the incubator were held constant at 37°C and 200 rpm, respectively. When cells reached exponential phase, aliquots of 0.5 mL was aseptically collected from the culture at intervals of 0.5 h. Measure the OD600 on a spectrometer to monitor the growth of the cells. For extracellular metabolite (shikimic acid and 3-dehydroshikimic acid) analysis by HPLC, culture samples were centrifuged for 5 min at 4°C and 20,000 g to remove the cells. The specific growth rate (μ) was defined as the increase in cell biomass concentration per unit time. The specific glucose consumption rate (qglu) was calculated during the exponential growth phase as the differential change in glucose with time (t) normalized to the biomass concentration. The specific SA production rate (qsa), and the specific 3-dehydroshikicmic acid production rate (qdhs) were calculated during the exponential growth phase as the differential change in SA, and 3-dehydroshikicmic acid with time (t) normalized to the biomass concentration, respectively.
Preparation of cellular hydrolysate from 13C-labeled biomass and measurements of 13C-labeled amino acids with gas chromatography and mass spectroscopy (GC-MS)
Preparation of cellular hydrolysate from 13C-labeled biomass for the GC-MS measurement followed standard protocols . For GC-MS measurement, the amino acids were derivatized with TBDMS and 1% tert-butyldimethylchlorosilane at 85°C for 60 min. The amounts of biomass used for GC-MS measurement was approximately 10 mg of cellular dry weight. GC-MS was carried out using an Agilent GC-7890A gas chromatograph equipped with an Agilent HP-5MS column (30 m × 0.25 mm × 0.25 mm) that was directly connected to an MS-5975C mass spectrometer (Agilent). The following amino acids were detected: Ala, Asp, Glu, Gly, His, Ile, Leu, Lys, Met, Phe, Pro, Ser, Thr, Tyr, and Val.
Determination of 13C-labeled SA and 3-dehydroshikimic acid with LC-MS
Agilent 1290 UHPLC system equipped with ZORBAX 300SB-C18 Column (2.1 mm × 100 mm × 1.8 μm) and Agilent 6530 Q-TOF mass spectrometer fitted with an electrospray ionization (ESI) source were used to measure 13C-labeling states of SA and 3-dehydroshikimic acid. SA and 3-dehydroshikimic acid was eluted with mixture of 0.1% formic acid in water (solution C) and acetonitrile (solution D) at ratio of 98:2 (solution C: solution D) and at a flow rate of 0.15 mL/min. MS data was collected in negative ion mode. The capillary voltage was set to 3500 kV, and the fragmentor voltage at 90 V. The drying gas temperature was maintained at 350°C with a flow rate of 11 L/min and a nebulizer pressure of 40 psi.
Metabolic net-flux analysis
13C MFA was a well-established method, and details concerning the general modeling framework of classical stationary 13C labeling experiments could be found elsewhere . In the case of U-13C labeling, several flux ratios, including (1) Pyr originating from malate, (2) PEP from OAA, (3) OAA from PEP, (4) PEP from transketolase, (5) PEP from the PP pathway, (6) pentose-5-phosphate from Glu-6-P, (7) OAA from the glyoxylate shunt, (8) erythrose-4-phosphate from pentose-5-phosphate, (9) Ser from Gly, (10) Gly from Ser, and (11) PEP and E4P from shikimic acid were calculated from the MDV M of precursors (Additional file 1: Table S2) with probabilistic equations. The calculating procedures were based on the general principle of flux ratio analysis [30–32] and were performed using a home-made Matlab program. Positional labeling pattern was also analyzed from cells grown exclusively on [1-13C] glucose to obtain precise information about the in vivo activities of the PPP and EMP pathways.
For quantification of carbon fluxes in B. subtilis, a bioreaction network was constructed. This network includes the reactions of the EMP, PP, and shikimate pathways, as well as the tricarboxylic acid (TCA) cycle and the glyoxylate shunt. The reactions catalyzed by PEP carboxykinase (v 10), pyruvate carboxylase (v 11), and malic enzyme (v 12) were also included. The networks of active pathways identified by flux ratio analysis were used for flux quantification. The following enzyme reactions were considered reversible: phosphoglucose isomerase (v 2), transketolase (v 5 and v 7), transaldolase (v 6), the sequence of glycolytic reactions leading from triose 3-phosphate to PEP (v 8) and malate dehydrogenase (v 17). The metabolic reactions considered as shown in the Additional file 1: Appendix S1.
All the above data, including (1) the stoichiometric reaction matrix, (2) the flux ratios derived from metabolic flux ratio analysis (Additional file 1: Table S3), (3) physiological data, and (4) precursor requirements for biomass synthesis, were utilized to calculate absolute fluxes in millimoles per gram of biomass per hour. Net fluxes were then estimated with Matlab by solving the stoichiometric matrix. The resulting flux distribution was the best fit to available data from both metabolite balances and tracer experiment-based flux ratios within the specified metabolic model.
Monte Carlo approach  was applied for statistical analysis.
This work was supported by 973 Project from Ministry of Science and Technology (No. 2012CB7211-04).
- Sheu TG, Deyde VM, Okomo-Adhiambo M, Garten RJ, Xu X, Bright RA, Butler EN, Wallis TR, Klimov AI, Gubareva LV: Surveillance for neuraminidase inhibitor resistance among human influenza A and B viruses circulating worldwide from 2004 to 2008. Antimicrob Agents Chemother. 2008, 52: 3284-3292. 10.1128/AAC.00555-08.View ArticleGoogle Scholar
- Dharan NJ, Gubareva LV, Meyer JJ, Okomo-Adhiambo M, McClinton RC, Marshall SA, George KS, Epperson S, Brammer L, Klimov AI, Bresee JS, Fry AM: Infections with oseltamivir-resistant influenza A(H1N1) virus in the United States. Jama-j Am Med Assoc. 2009, 301: 1034-1041. 10.1001/jama.2009.294.View ArticleGoogle Scholar
- Bradley D: Star role for bacteria in controlling flu pandemic?. Nat Rev Drug Discov. 2005, 4: 945-946. 10.1038/nrd1917.View ArticleGoogle Scholar
- Enrich LB, Scheuermann ML, Mohadjer A, Matthias KR, Eller CF, Newman MS, Fujinaka M, Poon T: Liquidambar styraciflua: a renewable source of shikimic acid. Tetrahedron Lett. 2008, 49: 2503-2505. 10.1016/j.tetlet.2008.02.140.View ArticleGoogle Scholar
- Lingens F: Biosynthesis of aromatic amino acids and its regulation. Angewandte Chemie-International Edition. 1968, 7: 350-10.1002/anie.196803501. 10.1002/anie.196803501.View ArticleGoogle Scholar
- Kramer M, Bongaerts J, Bovenberg R, Kremer S, Muller U, Orf S, Wubbolts M, Raeven L: Metabolic engineering for microbial production of shikimic acid. Metab Eng. 2003, 5: 277-283. 10.1016/j.ymben.2003.09.001.View ArticleGoogle Scholar
- Escalante A, Calderon R, Valdivia A, de Anda R, Hernandez G, Ramirez OT, Gosset G, Bolivar F: Metabolic engineering for the production of shikimic acid in an evolved Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system. Microb Cell Fact. 2010, 9: 21- 10.1186/1475-2859-9-21.View ArticleGoogle Scholar
- Dell KA, Frost JW: Identification and removal of impediments to biocatalytic synthesis of aromatics from D-glucose - rate-limiting enzymes in the common pathway of aromatic amino-acid biosynthesis. J Am Chem Soc. 1993, 115: 11581-11589. 10.1021/ja00077a065.View ArticleGoogle Scholar
- Knop DR, Draths KM, Chandran SS, Barker JL, von Daeniken R, Weber W, Frost JW: Hydroaromatic equilibration during biosynthesis of shikimic acid. J Am Chem Soc. 2001, 123: 10173-10182. 10.1021/ja0109444.View ArticleGoogle Scholar
- Chandran SS, Yi J, Draths KM, von Daeniken R, Weber W, Frost JW: Phosphoenolpyruvate availability and the biosynthesis of shikimic acid. Biotechnol Prog. 2003, 19: 808-814. 10.1021/bp025769p.View ArticleGoogle Scholar
- Yi J, Draths KM, Li K, Frost JW: Altered glucose transport and shikimate pathway product yields in E-coli. Biotechnol Prog. 2003, 19: 1450-1459. 10.1021/bp0340584.View ArticleGoogle Scholar
- Patnaik R, Liao JC: Engineering of Escherichia-Coli central metabolism for aromatic metabolite production with near theoretical yield. Appl Environ Microb. 1994, 60: 3903-3908.Google Scholar
- Patnaik R, Spitzer RG, Liao JC: Pathway engineering for production of aromatics in Escherichia-Coli - confirmation of stoichiometric analysis by independent modulation of Arog, Tkta, and Pps activities. Biotechnol Bioeng. 1995, 46: 361-370. 10.1002/bit.260460409.View ArticleGoogle Scholar
- Lu JL, Liao JC: Metabolic engineering and control analysis for production of aromatics: role of transaldolase. Biotechnol Bioeng. 1997, 53: 132-138. 10.1002/(SICI)1097-0290(19970120)53:2<132::AID-BIT2>3.0.CO;2-P.View ArticleGoogle Scholar
- Draths KM, Knop DR, Frost JW: Shikimic acid and quinic acid: replacing isolation from plant sources with recombinant microbial biocatalysis. J Am Chem Soc. 1999, 121: 1603-1604. 10.1021/ja9830243.View ArticleGoogle Scholar
- Rodriguez A, Martinez JA, Baez-Viveros JL, Flores N, Hernandez-Chavez G, Ramirez OT, Gosset G, Bolivar F: Constitutive expression of selected genes from the pentose phosphate and aromatic pathways increases the shikimic acid yield in high-glucose batch cultures of an Escherichia coli strain lacking PTS and pykF. Microb Cell Fact. 2013, 12: 86- 10.1186/1475-2859-12-86.View ArticleGoogle Scholar
- Iomantas YAV, Abalakina EG, Polanuer BM, Yampolskaya TA, Bachina TA, Kozlov YI: Method for producing shikimic acid. 2002, US patent 6436664 B1Google Scholar
- Sauer U, Hatzimanikatis V, Hohmann HP, Manneberg M, van Loon APGM, Bailey JE: Physiology and metabolic fluxes of wild-type and riboflavin-producing Bacillus subtilis. Appl Environ Microb. 1996, 62: 3687-3696.Google Scholar
- Davies WD, Davidson BE: The nucleotide sequence of aroG, the gene for 3-deoxy-D-arabinoheptulosonate-7-phosphate synthetase (phe) in Escherichia coli K12. Nucleic Acids Res. 1982, 10: 4045-4058. 10.1093/nar/10.13.4045.View ArticleGoogle Scholar
- Zurawski G, Gunsalus RP, Brown KD, Yanofsky C: Structure and regulation of aroH, the structural gene for the tryptophan-repressible 3-deoxy-D-arabino-heptulosonic acid-7-phosphate synthetase of Escherichia coli. J Mol Biol. 1981, 145: 47-73. 10.1016/0022-2836(81)90334-X.View ArticleGoogle Scholar
- Umbarger HE: Amino-acid biosynthesis and its regulation. Annu Rev Biochem. 1978, 47: 533-606. 10.1146/annurev.bi.47.070178.002533.View ArticleGoogle Scholar
- Fry B, Zhu T, Domach MM, Koepsel RR, Phalakornkule C, Ataai MM: Characterization of growth and acid formation in a Bacillus subtilis pyruvate kinase mutant. Appl Environ Microb. 2000, 66: 4045-4049. 10.1128/AEM.66.9.4045-4049.2000.View ArticleGoogle Scholar
- Sauer U, Eikmanns BJ: The PEP-pyruvate-oxaloacetate node as the switch point for carbon flux distribution in bacteria. FEMS Microbiol Rev. 2005, 29: 765-794. 10.1016/j.femsre.2004.11.002.View ArticleGoogle Scholar
- Zhu T, Pan Z, Domagalski N, Koepsel R, Ataai MM, Domach MM: Engineering of Bacillus subtilis for enhanced total synthesis of folic acid. Appl Environ Microb. 2005, 71: 7122-7129. 10.1128/AEM.71.11.7122-7129.2005.View ArticleGoogle Scholar
- Nguyen HD, Nguyen QA, Ferreira RC, Ferreira LCS, Tran LT, Schumann W: Construction of plasmid-based expression vectors for Bacillus subtilis exhibiting full structural stability. Plasmid. 2005, 54: 241-248. 10.1016/j.plasmid.2005.05.001.View ArticleGoogle Scholar
- Anagnostopoulos C, Spizizen J: Requirements for Transformation in Bacillus subtilis. J Bacteriol. 1961, 81: 741-746.Google Scholar
- Joseph P, Fantino JR, Herbaud ML, Denizot F: Rapid orientated cloning in a shuttle vector allowing modulated gene expression in Bacillus subtilis. FEMS Microbiol Lett. 2001, 205: 91-97. 10.1111/j.1574-6968.2001.tb10930.x.View ArticleGoogle Scholar
- Vagner V, Dervyn E, Ehrlich SD: A vector for systematic gene inactivation in Bacillus subtilis. Microbiology-Uk. 1998, 144: 3097-3104. 10.1099/00221287-144-11-3097.View ArticleGoogle Scholar
- Fuhrer T, Fischer E, Sauer U: Experimental identification and quantification of glucose metabolism in seven bacterial species. J Bacteriol. 2005, 187: 1581-1590. 10.1128/JB.187.5.1581-1590.2005.View ArticleGoogle Scholar
- Nanchen A, Fuhrer T, Sauer U: Determination of metabolic flux ratios from 13C-experiments and gas chromatography–mass spectrometry data: protocol and principles. Methods Mol Biol. 2007, 358: 177-197. 10.1007/978-1-59745-244-1_11.View ArticleGoogle Scholar
- Fischer E, Sauer U: Metabolic flux profiling of Escherichia coli mutants in central carbon metabolism using GC-MS. Eur J Biochem. 2003, 270: 880-891. 10.1046/j.1432-1033.2003.03448.x.View ArticleGoogle Scholar
- Yang C, Hua Q, Baba T, Mori H, Shimizu K: Analysis of Escherichia coli anaplerotic metabolism and its regulation mechanisms from the metabolic responses to altered dilution rates and phosphoenolpyruvate carboxykinase knockout. Biotechnol Bioeng. 2003, 84: 129-144. 10.1002/bit.10692.View ArticleGoogle Scholar
- Mollney M, Wiechert W, Kownatzki D, de Graaf AA: Bidirectional reaction steps in metabolic networks: IV. Optimal design of isotopomer labeling experiments. Biotechnol Bioeng. 1999, 66: 86-103. 10.1002/(SICI)1097-0290(1999)66:2<86::AID-BIT2>3.0.CO;2-A.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.