- Open Access
Examining Escherichia coli glycolytic pathways, catabolite repression, and metabolite channeling using Δpfk mutants
Biotechnology for Biofuelsvolume 9, Article number: 212 (2016)
Glycolysis breakdowns glucose into essential building blocks and ATP/NAD(P)H for the cell, occupying a central role in its growth and bio-production. Among glycolytic pathways, the Entner Doudoroff pathway (EDP) is a more thermodynamically favorable pathway with fewer enzymatic steps than either the Embden–Meyerhof–Parnas pathway (EMPP) or the oxidative pentose phosphate pathway (OPPP). However, Escherichia coli do not use their native EDP for glucose metabolism.
Overexpression of edd and eda in E. coli to enhance EDP activity resulted in only a small shift in the flux directed through the EDP (~20 % of glycolysis flux). Disrupting the EMPP by phosphofructokinase I (pfkA) knockout increased flux through OPPP (~60 % of glycolysis flux) and the native EDP (~14 % of glycolysis flux), while overexpressing edd and eda in this ΔpfkA mutant directed ~70 % of glycolytic flux through the EDP. The downregulation of EMPP via the pfkA deletion significantly decreased the growth rate, while EDP overexpression in the ΔpfkA mutant failed to improve its growth rates due to metabolic burden. However, the reorganization of E. coli glycolytic strategies did reduce glucose catabolite repression. The ΔpfkA mutant in glucose medium was able to cometabolize acetate via the citric acid cycle and gluconeogenesis, while EDP overexpression in the ΔpfkA mutant repressed acetate flux toward gluconeogenesis. Moreover, 13C-pulse experiments in the ΔpfkA mutants showed unsequential labeling dynamics in glycolysis intermediates, possibly suggesting metabolite channeling (metabolites in glycolysis are pass from enzyme to enzyme without fully equilibrating within the cytosol medium).
We engineered E. coli to redistribute its native glycolytic flux. The replacement of EMPP by EDP did not improve E. coli glucose utilization or biomass growth, but alleviated catabolite repression. More importantly, our results supported the hypothesis of channeling in the glycolytic pathways, a potentially overlooked mechanism for regulating glucose catabolism and coutilization of other substrates. The presence of channeling in native pathways, if proven true, would affect synthetic biology applications and metabolic modeling.
Escherichia coli have three native glycolytic pathways: EMPP, EDP, and OPPP. The EMPP employs ten enzymatic steps to yield two pyruvates, two ATP, and two NADH per glucose molecule , while OPPP serves as an oxidation route for NADPH synthesis. In E. coli, glucose metabolism mainly relies on the EMPP and the OPPP, while the EDP primarily remains inactive except during growth with gluconate . The EDP utilizes only five enzymes to produce one pyruvate, one glyceraldehyde-3-phosphate, and one NADPH per glucose molecule (Fig. 1). With further conversion of glyceraldehyde-3-phosphate via the lower EMPP, the two pathways (EDP and EMPP) result in the same net production of pyruvate. However, the EMPP contains thermodynamic bottlenecks, comprising of fructose 1,6-bisphosphate aldolase and triose-phosphate isomerase. The EDP avoids both these unfavorable reactions (at the expense of ATP yield), and requires substantially less enzymatic protein than the EMPP . Its end-products (glyceraldehyde 3-phosphate and pyruvate) are precursors of the non-mevalonate pathway, and thus, EDP upregulation has been proven to improve the yields of isoprenoids [4–6]. In addition, the EDP may also alleviate oxidative stress [7, 8] and improve NADPH generation without the loss of carbon (i.e., CO2) .
The EDP relies on two unique enzymes phosphogluconate dehydratase (edd) and 2-dehydro-3-deoxyphosphogluconate aldolase (eda), which are separate from the enzymes shared with the EMPP and OPPP. In wild-type (WT) E. coli, EDP flux is negligible. Δpgi (G6P→F6P, encoding glucose-6-phosphate isomerase), Δpfk (F6P→FBP, encoding phosphofructokinase), or overexpression of OPPP may redirect glycolytic fluxes and increase EDP activity [10, 11]. Δpgi mutant has measureable EDP flux (up to ~15 % of glycolysis), but grew and consumed glucose more slowly than the WT strain. Δpgi mutant may recover its growth rate by either increasing acetate overflow or activating the glyoxylate shunt after adaptive evolution . Moreover, some bacterial species with exclusive EDP activity (e.g., Rhodococcus opacus) can coutilize glucose with other substrates . Thus, we were curious if the presence of a highly active EDP could alleviate E. coli glucose catabolite repression and promote cell growth. To reorganize the glycolytic pathways, we engineered and characterized ΔpfkA mutants (deletion of phosphofructokinase I, important in the regulation of glycolysis-gluconeogenesis, Fig. 1). Through growth experiments and 13C-labeling, we elucidated E. coli physiological changes and its capability for simultaneous utilization of carbon substrates. Finally, ΔpfkA mutants also allowed us to evaluate the hypothesis of the existence of metabolite channeling in glycolytic pathways .
Overexpression of EDP in E. coli
Overexpressing edd and eda in the wild-type (WT) strain reduced the flux (~20 %) through EMPP and OPPP. Despite the increased EDP activity, the growth rate of the EDP overexpressing strain decreased by ~30 % compared to WT (Table 1). To confirm the tradeoff from metabolic burdens imposed by antibiotics and plasmid/protein expression, we replaced EDP genes with the gene encoding for yellow fluorescent protein (YFP). The YFP-expressing strain had a similar reduced growth. A recent study found that the insertion of EDP (pgi, zwf, pgl, edd, eda from Z. mobilis) into E. coli chromosome also repressed the growth of the engineered strain .
To further improve EDP flux, the Keio collection ΔpfkA (phosphofructokinase I, F6P→FBP) mutant was used . The phosphofructokinase has two isoenzymes (pfkA and pfkB) with pfkA as the primary enzyme responsible for the conversion of F6P to FBP . 13C-fingerprinting revealed that the ΔpfkA mutant (JW3887) distributed glucose flux through the OPPP (~62 %), the EDP (~14 %) and the EMPP (~24 %) (Fig. 1). The pfkA knockout increased cell lag phase, and reduced both cell growth rate and acetate overflow (Additional file 1: Figure S1). This is because the glycolytic flux reorganization can cause metabolic burdens, cofactor imbalances, and reduced carbon yield due to CO2 loss from the high OPPP flux . Overexpressing EDP in the ΔpfkA mutant (resulting in strain WH04) raised the EDP flux to ~72 % and reduced the EMPP flux to ~18 % (Fig. 1). WH04 grew much faster than WH02 (the ∆pfkA mutant with YFP overexpression) (Table 1), demonstrating a beneficial impact of EDP in the absence of EMPP.
Removal of glucose, carbon catabolite repression
Lignocellulose hydrolysate contains glucose, xylose, and acetate. However, carbon catabolite repression (CCR) inhibits hosts from consuming diverse substrates simultaneously (i.e., E. coli mainly consumes glucose first in a glucose/xylose medium). The CCR is often explained by glucose inhibition on the synthesis of enzymes involved in catabolism of other carbon sources, while recent studies suggest that the presence of subpopulations within E. coli cultures can have different carbon utilization hierarchies [18, 19]. To investigate glucose catabolite repression, we grew the strains with a mixed carbon source of 10 g/L of glucose and 6 g/L of xylose, and measured the consumption of the two sugars (Fig. 2a, b). The ΔpfkA mutant (JW3887) simultaneously consumed glucose and xylose, and xylose coutilization nearly doubled the growth rate compared to its glucose-only culture (Fig. 2a). WH04 can also uptake xylose with glucose, and xylose addition reduced the strain’s lag phase (Fig. 2b). Moreover, the assimilation of xylose into biomass was determined through 13C-experiments (cultures fed with 13C6-glucose and unlabeled xylose). The fraction of xylose incorporated into proteinogenic amino acids during the exponential phase was measured (Fig. 2c). For both JW3887 and WH04, the xylose incorporation was ~50 % (Fig. 2c), much higher utilization than in the WT strain. Therefore, the removal of one pivotal EMP enzyme, phosphofructokinase I, alleviated glucose catabolite repression.
Acetate, another common component in lignocellulosic hydrolysate, has to be metabolized through gluconeogenesis . ΔpfkA is expected to alleviate the EMPP repression on gluconeogenesis, and improve the coutilization of glucose and acetate. To test this hypothesis, the mutants (JW3887 and WH04) were grown with 10 g/L of glucose and supplemented with 6 g/L of sodium acetate during the early exponential phase (Fig. 3a–c). Although acetate may interfere with intracellular pH and hinder E. coli growth [21, 22], JW3887 was able to coutilize acetate with little effect on its growth. Acetate incorporation into biomass measured via cultures with unlabeled glucose and fully labeled acetate revealed that the mutant utilized acetate for biomass synthesis with enrichment of 13C into TCA and glycolysis-derived proteinogenic amino acids (Fig. 3d). In contrast, the glycolytic amino acids: Ala, Phe, and Ser, were unlabeled in WT under similar mixed carbon cultures, these amino acids being derived from PYR, PEP, and 3PG, respectively. Therefore, ΔpfkA (JW3887) could coutilize acetate to generate TCA cycle metabolites and glycolytic intermediates via upregulation of gluconeogenesis. However, introduction of EDP into the ΔpfkA mutant inhibited gluconeogenesis, as revealed with a lack of 13C Ala/Phe/Ser in WH04 biomass (Fig. 3d).
Dynamic labeling pattern for sugar phosphates in ΔpfkA mutants
We grew WT and the ΔpfkA mutants (WH04 and JW3887) with unlabeled glucose into the mid-exponential growth phase, and then, we added uniformly 13C6-labeled glucose into the cultures. The resulting kinetics of 13C-labeling incorporation in key metabolites was examined. The WT demonstrated very fast metabolite turnover, such that most metabolites’ labeling reached isotopic steady state in 15 s (Fig. 4a and Additional file 2: Fig. S2). This observation is consistent with the rapid EMPP turnover rates previously reported . In contrast, WH04 had lower rates of labeling incorporation. Interestingly, G6P and 6PG appeared to be labeled slower than their downstream metabolites (3PG and PEP) (Fig. 4b, Additional file 3: Fig. S3). Similarly, the 13C-pulse experiment for JW3887 also revealed that 3PG and R5P reached high isotopic ratios more quickly than their higher pathway metabolites (G6P and 6PG), revealing unusual labeling patterns in glycolytic intermediates (Additional file 4: Fig. S4).
The unsequential 13C labeling patterns for glycolysis intermediates in ΔpfkA mutants could be explained by metabolite channeling (Fig. 5a): intermediates in metabolic pathways are passed from enzyme to enzyme without equilibration within the cellular medium . For example, a large unlabeled metabolite pool formed outside of the ‘channel’ could dilute the labeled metabolite concentration, resulting in a slow labeling incorporation of the bulk metabolite, as measured by LC–MS methods. The downstream metabolites, however, retain their fast rates with their production primarily from the metabolite pool within the ‘channel’. In a previous study, in vivo evidence discovered that EMPP intermediates are mainly concentrated within a ‘channel’ with minimal mixing with the cytoplasmic pool in E. coli . In vivo studies in eukaryotes and in vitro studies also suggest that EMPP and OPPP metabolites (e.g., 6PG) are highly channeled [25, 26]. In this study, disruption of EMPP channel by the knockout of pfkA caused a bottleneck in the conversion of F6P to FBP, which likely allowed hexose6P metabolites to accumulate within the cytosol. LC–MS peak abundances and isotopomer ratio analysis (Fig. 4c, d) further confirmed that ∆pfkA (JW3887) and WH04 had larger G6P and 6PG pool sizes than those in the WT. During 13C-pulse experiments, the unlabeled G6P and 6PG amassed in the cytosol could significantly slow the measured labeling of bulk hexose6P, causing 13C to appear more gradually than in their downstream metabolites (3PG and PEP). There could be other reasons for the observed labeling patterns in ∆pfkA mutants, such as activation of unknown pathways associated with glycolysis (e.g., Δedd surprisingly increases glycogen accumulation ), experimental artifacts during 13C-metabolite sampling/extraction or heterogenity due to different cell sub-populations within the mutant cultures. However, channeling is a logical explanation given the labeling pattern difference between the WT and Δpfk mutants for only certain intermediates.
Fluxome response to changes in glycolytic pathway utilization
Flux balance analysis (FBA) was used to predict growth and flux distribution under different glycolytic pathways and carbon substrates . Using the measured glycolytic flux ratios, substrate utilization, and acetate production as constraints, we performed FBA simulations on the mutant strains using maximal biomass growth as the objective function (Fig. 6a, b). FBA-predicted growth rate of WH04 agreed with the experimentally measured the glucose-based growth rate, while xylose and acetate co-utilization models show certain discrepanies. The measured growth rates of JW3887 were generally below the growth rates predicted by FBA (Fig. 6c). The knockout of pfkA may have unknown effects that cause a suboptimal metabolism for biomass synthesis. Unlike the WT strain, FBA predicted that ΔpfkA mutant and WH04 upregulated the glyoxylate shunt and maintained low activity of anaplerotic pathways for glucose metabolism (Fig. 6a). Compared to WT, the absolute fluxes of energy production (ATP, NADH, and NADPH) rate in mutants were reduced (Fig. 6b). FBA improved our overview of the central metabolic responses to our glycolysis flux reorganization.
Glycolytic pathways hold considerable control over carbon utilization and biosynthetic efficiency. The knowledge of their regulation is crucial for designing optimal microbial production strains. For example, the desire to use cheap cellulosic feedstock has prompted metabolic engineers to develop microbial hosts that utilize both 6- and 5-carbon sugars, simultaneously [29, 30]. However, carbon catabolite repression (CCR) prevents the simultaneous consumption of other carbon sources. CCR is believed to act in E. coli by preferentially directing resources to consume glucose more efficiently by regulating its sugar transporters and carbon degradation pathways . Approaches such as adaptive evolution and overexpression of xylose transport and catabolic genes have been tested to minimize CCR . The knockout of phosphoenolpyruvate-dependent glucose phosphotransferase system (PTS, the primary transport system for glucose) has also been shown to be an effective way to enable glucose and xylose coutilization [32–35], but has a negative effect on glucose uptake and growth rates. In this study, ΔpfkA mutants also showed no CCR and employed xylose co-consumption to increase its growth rate. ΔpfkA caused an increased pool size of G6P that may have inhibited PTS and decreased biomass accumulation [36, 37].
Gluconeogenesis is usually inactive under high glucose concentrations as it is regulated by phosphofructokinase and pyruvate kinase. The lack of EMPP flux could facilitate gluconeogenesis for utilization of non-sugar based substrates. This is seen in R. opacus, which lacks pfk activity and shows simultaneous gluconeogenesis and EDP for sugar and phenol co-utilization . In this study, ΔpfkA mutants actively utilized acetate from the glucose medium. This observation is encouraging since acetate is a notorious inhibitor in E. coli cultivation, because acetate can freely diffuse through the cell membrane and disrupt the intracellular pH [21, 38, 39]. 13C-labeling analysis of Asp and Glu revealed that both JW3887 and WH04 actively employed the TCA cycle to degrade 13C-acetate in the presence of glucose. However, increased EDP flux in WH04 repressed the gluconeogenesis activity, since EDP overexpression drives high flux toward the lower portion of the EMPP.
Despite being more thermodynamically favorable and requiring ~5-fold less enzymatic protein , the E. coli mutants grew slower than WT after overexpression of EDP genes. One reason could be the lower ATP yield of EDP as compared to EMPP (this explanation is questionable since respiration remains the main ATP production route in all aerobic cultures via our FBA simulations). Another explanation is metabolite channeling. The intracellular environment contains high amounts of macromolecules (30–60 % of cell volume) creating a crowding effect within the cell that can reduce metabolite diffusion, and thus, channeling is an important mechanism for in vivo enzyme reactions . Channeling can be naturally accomplished in the cell via electrostatic interactions, intramolecular tunnels, and small spatial distances between enzymes . In eukaryotic cells, channeling is present for the TCA cycle and other pathways, as eukaryotes have organized their pathways within organelles [25, 26]. Pathways in E. coli have also been found to demonstrate forms of channeling, such as for fatty acid type II synthesis pathway . As shown in this study and previous studies , E. coli EMPP channeling would significantly improve reaction thermodynamics and overcome kinetic obstacles. Thereby, it could be an important factor in WT and WH03 preference for EMPP rather than the EDP route.
The channeling process is difficult to study non-invasively, particularly in the central metabolism. In an early study by Shearer et al., E. coli was engineered so that the mutants could uptake both 14C-glucose and unlabeled sugar phosphates . The authors found that the presence of unlabeled FBP in the culture medium did not dilute 14C EMPP intermediates or 14CO2, indicating a high degree of EMPP channeling. In this study, we performed 13C-pulse experiments to measure the labeling dynamics in WT and the ΔpfkA mutants (Additional files 2, 3, 4: Figure S2–S4). Our results supported the hypothesis that metabolites are not equilibrated in the cytoplasm (Fig. 5a, b). Previous 13C-pulse studies also examined the labeling dynamics for wild-type E. coli and B. subtilis, but the metabolites reached isotopic steady state within seconds of 13C-pulse [44, 45]. In this study, ΔpfkA mutants had slower metabolite transfer rates, and accumulated the upstream hexose6Ps outside of the channel (Fig. 5b). The resulting presence of two pools of the same metabolite allowed 13C-pulse experiments to detect unsequential metabolite labeling through glycolysis.
Channeling may have several implications on our understanding of cell physiology. First, metabolite channeling avoids diffusion limitations and thus significantly improves bioconversion efficiency . Channeling would explain the robustness of natively evolved pathways for biosynthesis, encouraging synthetic biologists to replicate channeling as a strategy to improve heterologous pathway efficiency (e.g., colocalization/compartmentalization of overexpressed enzymes) . Second, it necessitates that we rethink the in vivo Gibbs free energy or kinetics typically reported in studies, since the global metabolite pool measurement does not reflect the local substrate concentrations (The notion of ‘one perfectly mixed solution’ may be a poor reflection of the cell ‘in vivo’). This presents complications with metabolic modeling. Third, metabolite channeling could affect steady-state 13C-metabolic flux analysis, because channeling may prevent the carbon randomization during conversion of a symmetric metabolite or introduce bypass routes for 13C labeling . Finally, we hypothesize that channeling could be another factor behind CCR. At high growth rates, the EMPP is capable of reaching fast glucose catabolic rates, such that metabolites derived from other substrates (e.g., xylose) cannot sufficiently diffuse into the EMPP channel (i.e., glucose catabolite repression). During slow growth conditions, glucose catabolic flux decreases. This favors diffusion of metabolites from cytosol into the EMPP channel and encourage the coutilization of the secondary substrate (Fig. 5c, d). In reality, inhibition of EMPP flux (e.g., removal of glucose transporter ΔptsG or glucose-limiting chemostat cultures) allows E. coli to cometabolize glycerol with glucose .
We have rewired the E. coli central metabolic network through a pfkA knockout and EDP overexpression. The altered strains alleviated glucose catabolite repression, and could be beneficial in biosynthesis from renewable cellulosic hydrolysates. Physiological analyses of different glycolytic strategy clearly revealed that despite the theoretical prediction that EDP could be a preferred glycolysis pathway based on both thermodynamics and the cost of enzymatic protein, the cell prefers EMPP for optimal growth. The results of 13C-pulse experiment revealed a form of glycolysis channeling that, if proven, could significantly affect our understanding of reaction thermodynamics, flux analysis, and synthetic biology applications.
Glucose, [1-13C] glucose, [U-13C] glucose, [U-13C] acetate, and all other chemicals unless otherwise stated were purchased from Sigma Aldrich (St. Louis, MO). Q5 High Fidelity 2X Master Mix was obtained from New England Biolabs (Ipswich, MA), and all other enzymes were purchased through Thermo Scientific (Waltham, MA).
Construction of strains and plasmids
All plasmids were developed from a BglBrick expression vector, pBbE5c-YFP, which contains a ColE1 origin of replication, chloramphenicol resistance and LacUV5 promoter . The genes, phosphogluconate dehydratase (edd) and 2-dehydro-3-deoxyphosphogluconate aldolase (eda) were PCR-amplified from genomic DNA of E. coli K-12 MG1655. The plasmid pEDW was assembled from the PCR-amplified genomic region from edd to eda into pBbE5c backbone via Gibson assembly. Primers used for plasmid construction are the following:
Table 1 lists the resulting mutant strains, their JBEI ICE registry numbers, and growth rates. All strains generated at JBEI are available thought the JBEI Registry (https://acs-registry.jbei.org; ). The plasmid was first transformed into DH10B for plasmid propagation, and then verified by sequencing completed by Quintara (San Francisco, CA). After verification for no mutations, the plasmid was transformed into Keio collection strains, BW25113 (WT) and the ∆pfkA mutant (JW3887) , obtained from the Yale E. coli Genetic Stock Center (New Haven, CT).
Media composition and growth conditions
Seed cultures were grown in LB for ~12 h, then inoculated (2 % v/v) in M9 minimal medium overnight. Strains were then inoculated in fresh M9 minimal medium (~2 % v/v). The M9 medium (pH = 7–7.5) consisted of 11 g/L of 5× M9 salts (Sigma Aldrich, St. Louis, MO), 2 mM MgSO4·7H2O, 0.1 mM CaCl2·2H2O, 1 mg/L of Thiamine HCl, 3 mg/L of FeSO4·7H2O, and trace minerals. Initial sugar concentrations were 10 g/L of glucose and 6 g/L of xylose. For mutant cultures, 30 mg/L chloramphenicol and 10 mg/L kanamycin were supplied as needed. 50 µM Isopropyl β-d-1-thiogalactopyranoside (IPTG) was used to induce plasmid borne genes at early exponential growth phase. To study acetate coutilization, 6 g/L sodium acetate (or 6 g/L NaCl served as a control) was added at early phase of mutant cultures (OD600 ~ 0.2, pH ~ 7.5) with the IPTG induction. All cultivations were aerobic and conducted in 5 mL cultures in test tubes on a rotary shaker at 250 rpm at 37 °C. For labeling experiments, unlabeled glucose was replaced with 100 % [1-13C] glucose or 100 % [U-13C] glucose, and unlabeled acetate was replaced with [U-13C] acetate.
Strains were grown in approximately 70 mL cultures in 250-mL Erlenmeyer flasks on a rotary shaker at 250 rpm with 2 g/L unlabeled glucose M9 media. In all 13C-pulse experiments, cultures in exponential growth phase (pH = 6–7) were pulsed with 5 mL–60 g/L 13C6-glucose stock solution. To measure the 13C-incorporation into glycolysis metabolites over time, cultures were harvested at different time points (from 15 s to 3 min). The harvested samples were quenched using the procedure of Fast-Cooling . Specifically, ~10 mL culture was poured into a 50-mL falcon tube containing 2 mL ice-cold M9 medium (no carbon source) and the tube was immediately immersed in liquid N2. To facilitate heat transfer and avoid ice formation in the sample solution during the liquid N2 bath, the sample solutions were manually agitated using a digital thermometer. The sample solution could be cooled to 0 °C in 10 s. The ice-cold samples were centrifuged at 0 °C for 3 min, and the pellets were stored at −80 °C until metabolite extraction.
Supernatant samples were taken in parallel experiments separate from cultures used for optical density measurement to reduce the loss of volume. At each time point, at least 150 µL of culture was extracted, centrifuged to remove cell biomass, and stored at −20 °C. Samples were diluted by 2× and then filtered through 0.22 μm centrifugal filters. Samples were run on an Agilent Technologies 1200 series HPLC equipped with an Aminex H column, and the concentrations were estimated based on standard curves .
Amino acid extraction and GC–MS analysis
Amino acid extraction and GC–MS analysis were performed as described previously . Briefly, cell pellets from 5 mL cultures taken during the exponential phase were washed with 0.9 % (w/v) of NaCl solution, and then hydrolyzed in 6 M HCl at 100 °C. The resulting amino acids were derivatized by N-tert-butyldimethylsilyl-N-methyltrifluoroacetamide prior to GC–MS analysis. For isotopic tracing, we used the fragments [M-15]+ or [M-57]+ (containing the entire amino acid backbone) and [M-159]+ or [M-85]+ (containing the amino acid backbone after loss of its first carbon). The natural isotopic abundance of derivatized amino acids was corrected using a reported algorithm . The mass isotopomer distributions of alanine and serine for the strains used are given in Additional file 5: Table S1.
1-13C glucose was used to measure in vivo glucose catabolism: the OPPP cleaves the labeled carbon as CO2; the EDP produces the first carbon labeled pyruvate, while EMPP results in the labeled carbon present in the third position of pyruvate . GC–MS analysis of alanine (synthesized from pyruvate) can reveal this positional labeling as further described previously [8, 13, 54]. Equations (1–3) were solved for estimation of the flux ratios between the three pathways (without considering metabolite channeling effect):
where v EDP, v EMPP, and v OPPP are the estimated fluxes for the EDP, EMPP, and OPPP, respectively. Ala xxx are the fractions of the isotopomer determined by fragmentation data, with 1 denoting a 13Carbon isotope. For example, Ala 100 is the fraction of alanine labeled at carbon position 1, a unique isotopomer produced by the EDP. The calculated flux ratios are estimations, which assume the reactions are irreversible.
Metabolite extraction and LC–MS/MS analysis
After quenching, the cell pellet was suspended in 1 mL of cooled 7:3 methanol:chloroform mixture. These samples were then placed on a rotary shaker at 250 rpm, overnight at 4 °C. The methanol layer was separated with the addition of 500 µL of water, then filtered through an Amicon Ultra centrifuge filter (3000 Da; EMD Millipore, Billerica, MA), lyophilized (−50 °C) and reconstituted in 100 µL of acetonitrile–water (6:4, v/v). The samples were then analyzed via liquid chromatography-mass spectrometry (LC–MS) using a SeQuant Zic-pHILIC column (EMD Millipore, Billerica, MA, USA) in an Agilent Technologies 1200 Series HPLC system. The mobile phase was composed of 15 mM ammonium carbonate (Sigma-Aldrich, St. Louis, MO, USA) in water (solvent A) and 15 mM ammonium carbonate in 75 % acetonitrile and 25 % water (solvent B). A flow rate of 0.2 mL/min was used, unless stated otherwise. Metabolites were separated via gradient elution under the following conditions: 100 % B (0 min), 82 % B (4.4 min), 72 % B (7.7 min), 60 % B (9.7 min), 100 % B (10.2 min), 100 % B (12.5 min), 100 % B (13 min, 0.4 mL/min), and 100 % B (18.5 min, 0.4 mL/min). The HPLC system was coupled to an Agilent Technologies 6210 series time-of-flight mass spectrometer (for LC–TOF MS) via a MassHunter workstation (Agilent Technologies, USA). A split ratio of 1:4 was used throughout. Drying and nebulizing gases were set to 10 L/min and 25 lb/in.2, respectively, and a drying-gas temperature of 300 °C was used throughout. Electrospray ionization (ESI) was conducted in the negative ion mode and a capillary voltage of 3500 V was utilized. The acquisition range was from 70 to 1000 m/z, and the acquisition rate was 0.86 spectra/s. Metabolite mass isotopomer distribution was determined based on the ratio of the integrated peak area of the chosen isotopomer to the sum of all the integrated peak areas of the possible isotopomers for the given metabolite.
Estimation of pool size via isotope-ratio approach
The relative pool size of key metabolites (G6P, 6PG, and 3PG) was measured via an isotope-based ratio approach, as modified from . Wild-type (BW25113) was cultivated in 5 g/L of 13C6-glucose and 1 g/L of fully labeled 13C-sodium bicarbonate, which generated 13C-metabolites as internal standards. The fully labeled WT biomass and the unlabeled mutant cultures were quenched using liquid N2 as described earlier, then mixed together at 1:1 volume ratio for metabolite extraction. The ratio of WT metabolite pool size (labeled peak area) to the mutant strains metabolite pool size (unlabeled peak area) was determined using LC–MS, which provided a qualitative comparison of metabolite pool sizes between WT and mutants.
Flux balance analysis constrained by measured flux ratios
The E. coli Genome Scale Model iJO1366 (2251 metabolic reactions) was adapted for metabolic modeling of mutants under the different growth conditions . Flux balance analysis (FBA) used the following objective function:
Subject to: S · v = 0
where S represents the stoichiometric matrix, v represents the matrix of reaction rates (fluxes), while matrices lb and ub are the lower and upper bounds, respectively. Glucose, xylose, and acetate uptake rates were fixed to measured average values with 5 % presumed variation. Flux ratios determined by 13C-analysis were used to constrain fluxes through EMPP, EDP, and OPPP (Fig. 1). As the activities of pfkA and pfkB were lumped in iJO1366, the PFK reaction was not constrained. All generated fluxes were normalized to total substrate uptake rate.
Romano AH, Conway T. Evolution of carbohydrate metabolic pathways. Res Microbiol. 1996;147:448–55.
Eisenberg RC, Dobrogosz WJ. Gluconate metabolism in Escherichia coli. J Bacteriol. 1967;93:941–9.
Flamholz A, Noor E, Bar-Even A, Liebermeister W, Milo R. Glycolytic strategy as a tradeoff between energy yield and protein cost. Proc Natl Acad Sci. 2013;110:10039–44.
Ng CY, Farasat I, Maranas CD, Salis HM. Rational design of a synthetic Entner–Doudoroff pathway for improved and controllable NADPH regeneration. Metab Eng. 2015;29:86–96.
Li C, Ying L-Q, Zhang S-S, Chen N, Liu W-F, Tao Y. Modification of targets related to the Entner–Doudoroff/pentose phosphate pathway route for methyl-d-erythritol 4-phosphate-dependent carotenoid biosynthesis in Escherichia coli. Microb Cell Fact. 2015;14:1–12.
Liu H, Wang Y, Tang Q, Kong W, Chung W-J, Lu T. MEP pathway-mediated isopentenol production in metabolically engineered Escherichia coli. Microb Cell Factor. 2014;13:1–8.
Chavarría M, Nikel PI, Pérez-Pantoja D, de Lorenzo V. The Entner–Doudoroff pathway empowers Pseudomonas putida KT2440 with a high tolerance to oxidative stress. Environ Microbiol. 2013;15:1772–85.
Klingner A, Bartsch A, Dogs M, Wagner-Dobler I, Jahn D, Simon M, Brinkhoff T, Becker J, Wittmann C. Large-scale 13C-flux profiling reveals conservation of the Entner-Doudoroff pathway as a glycolytic strategy among marine bacteria that use glucose. Appl Environ Microbiol. 2015;81:2408–22.
He L, Xiao Y, Gebreselassie N, Zhang F, Antoniewicz MR, Tang YJ, Peng L. Central metabolic responses to the overproduction of fatty acids in Escherichia coli based on 13C-metabolic flux analysis. Biotechnol Bioeng. 2014;111:575–85.
Toya Y, Ishii N, Nakahigashi K, Hirasawa T, Soga T, Tomita M, Shimizu K. 13C-metabolic flux analysis for batch culture of Escherichia coli and its pyk and pgi gene knockout mutants based on mass isotopomer distribution of intracellular metabolites. Biotechnol Prog. 2010;26:975–92.
Seol E, Sekar BS, Raj SM, Park S. Co-production of hydrogen and ethanol from glucose by modification of glycolytic pathways in Escherichia coli—from Embden-Meyerhof-Parnas pathway to pentose phosphate pathway. Biotechnol J. 2016;11:249–56.
Fong SS, Nanchen A, Palsson BO, Sauer U. Latent pathway activation and increased pathway capacity enable Escherichia coli adaptation to loss of key metabolic enzymes. J Biol Chem. 2006;281:8024–33.
Hollinshead WD, Henson WR, Abernathy M, Moon TS, Tang YJ. Rapid metabolic analysis of Rhodococcus opacus PD630 via parallel 13C-metabolite fingerprinting. Biotechnol Bioeng. 2016;113:91–100.
Shearer G, Lee JC, Koo JA, Kohl DH. Quantitative estimation of channeling from early glycolytic intermediates to CO2 in intact Escherichia coli. FEBS J. 2005;272:3260–9.
Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006.
Fraenkel DG. Mutants in glucose metabolism. Annu Rev Biochem. 1986;55:317–37.
Seol E, Ainala SK, Sekar BS, Park S. Metabolic engineering of Escherichia coli strains for co-production of hydrogen and ethanol from glucose. Int J Hydrogen Energy. 2014;39(33):19323–30.
Koirala S, Wang X, Rao CV. Reciprocal regulation of l-arabinose and d-xylose metabolism in Escherichia coli. J Bacteriol. 2016;198:386–93.
Beisel CL, Afroz T. Rethinking the hierarchy of sugar utilization in bacteria. J Bacteriol. 2016;198:374–6.
Zhao J, Shimizu K. Metabolic flux analysis of Escherichia coli K12 grown on 13C-labeled acetate and glucose using GC–MS and powerful flux calculation method. J Biotechnol. 2003;101:101–17.
Xiao Y, Ruan Z, Liu Z, Wu SG, Varman AM, Liu Y, Tang YJ. Engineering Escherichia coli to convert acetic acid to free fatty acids. Biochem Eng J. 2013;76:60–9.
Luli GW, Strohl WR. Comparison of growth, acetate production, and acetate inhibition of Escherichia coli strains in batch and fed-batch fermentations. Appl Environ Microbiol. 1990;56:1004–11.
Nöh K, Grönke K, Luo B, Takors R, Oldiges M, Wiechert W. Metabolic flux analysis at ultra short time scale: isotopically non-stationary 13C labeling experiments. J Biotechnol. 2007;129:249–67.
Williams TCR, Sweetlove LJ, Ratcliffe RG. Capturing metabolite channeling in metabolic flux phenotypes. Plant Physiol. 2011;157:981–4.
Malaisse WJ, Zhang Y, Sener A. Enzyme-to-enzyme channeling in the early steps of glycolysis in rat pancreatic islets. Endocrine. 2004;24:105–9.
Debnam PM, Shearer G, Blackwood L, Kohl DH. Evidence for channeling of intermediates in the oxidative pentose phosphate pathway by soybean and pea nodule extracts, yeast extracts, and purified yeast enzymes. Eur J Biochem. 1997;246:283–90.
Long CP, Gonzalez JE, Sandoval NR, Antoniewicz MR. Characterization of physiological responses to 22 gene knockouts in Escherichia coli central carbon metabolism. Metab Eng. 2016;37:102–13.
Burgard AP, Maranas CD. Probing the performance limits of the Escherichia coli metabolic network subject to gene additions or deletions. Biotechnol Bioeng. 2001;74:364–75.
Kim J-H, Block DE, Mills DA. Simultaneous consumption of pentose and hexose sugars: an optimal microbial phenotype for efficient fermentation of lignocellulosic biomass. Appl Microbiol Biotechnol. 2010;88:1077–85.
Nieves LM, Panyon LA, Wang X. Engineering sugar utilization and microbial tolerance toward lignocellulose conversion. Front Bioeng Biotechnol. 2015;3:17.
Gorke B, Stulke J. Carbon catabolite repression in bacteria: many ways to make the most out of nutrients. Nat Rev Microbiol. 2008;6:613–24.
Lu H, Zhao X, Wang Y, Ding X, Wang J, Garza E, Manow R, Iverson A, Zhou S. Enhancement of d-lactic acid production from a mixed glucose and xylose substrate by the Escherichia coli strain JH15 devoid of the glucose effect. BMC Biotechnol. 2016;16:1–10.
Su B, Wu M, Zhang Z, Lin J, Yang L. Efficient production of xylitol from hemicellulosic hydrolysate using engineered Escherichia coli. Metab Eng. 2015;31:112–22.
Jung I-Y, Lee J-W, Min W-K, Park Y-C, Seo J-H. Simultaneous conversion of glucose and xylose to 3-hydroxypropionic acid in engineered Escherichia coli by modulation of sugar transport and glycerol synthesis. Bioresour Technol. 2015;198:709–16.
Chiang C-J, Lee HM, Guo HJ, Wang ZW, Lin L-J, Chao Y-P. Systematic approach to engineer Escherichia coli pathways for co-utilization of a glucose–xylose mixture. J Agric Food Chem. 2013;61:7583–90.
Chassagnole C, Noisommit-Rizzi N, Schmid JW, Mauch K, Reuss M. Dynamic modeling of the central carbon metabolism of Escherichia coli. Biotechnol Bioeng. 2002;79:53–73.
Morita T, El-Kazzaz W, Tanaka Y, Inada T, Aiba H. Accumulation of glucose 6-phosphate or fructose 6-phosphate is responsible for destabilization of glucose transporter mRNA in Escherichia coli. J Biol Chem. 2003;278:15608–14.
Ding J, Holzwarth G, Penner MH, Patton-Vogt J, Bakalinsky AT. Overexpression of acetyl-CoA synthetase in Saccharomyces cerevisiae increases acetic acid tolerance. FEMS Microbiol Lett. 2015;362:1–7.
Repaske DR, Adler J. Change in intracellular pH of Escherichia coli mediates the chemotactic response to certain attractants and repellents. J Bacteriol. 1981;145:1196–208.
Zhang YHP. Substrate channeling and enzyme complexes for biotechnological applications. Biotechnol Adv. 2011;29:715–25.
Wheeldon I, Minteer SD, Banta S, Barton SC, Atanassov P, Sigman M. Substrate channelling as an approach to cascade reactions. Nat Chem. 2016;8:299–309.
Janßen HJ, Steinbüchel A. Fatty acid synthesis in Escherichia coli and its applications towards the production of fatty acid based biofuels. Biotechnol Biofuels. 2014;7:7.
Noor E, Bar-Even A, Flamholz A, Reznik E, Liebermeister W, Milo R. Pathway thermodynamics highlights kinetic obstacles in central metabolism. PLoS Comput Biol. 2014;10:e1003483.
Link H, Kochanowski K, Sauer U. Systematic identification of allosteric protein-metabolite interactions that control enzyme activity in vivo. Nat Biotechnol. 2013;31:357–61.
Millard P, Massou S, Wittmann C, Portais J-C, Létisse F. Sampling of intracellular metabolites for stationary and non-stationary 13C metabolic flux analysis in Escherichia coli. Anal Biochem. 2014;465:38–49.
Sarria S, Wong B, Martín HG, Keasling JD, Peralta-Yahya P. Microbial synthesis of pinene. ACS Synth Biol. 2014;3:466–75.
Yao R, Xiong D, Hu H, Wakayama M, Yu W, Zhang X, Shimizu K. Elucidation of the co-metabolism of glycerol and glucose in Escherichia coli by genetic engineering, transcription profiling, and 13C metabolic flux analysis. Biotechnol Biofuels. 2016;9:175.
Lee TS, Krupa RA, Zhang F, Hajimorad M, Holtz WJ, Prasad N, Lee SK, Keasling JD. BglBrick vectors and datasheets: a synthetic biology platform for gene expression. J Biol Eng. 2011;5:12–3.
Ham TS, Dmytriv Z, Plahar H, Chen J, Hillson NJ, Keasling JD. Design, implementation and practice of JBEI-ICE: an open source biological part registry platform and tools. Nucleic Acids Res. 2012;40:e141.
Fu Y, Yoon JM, Jarboe L, Shanks JV. Metabolic flux analysis of Escherichia coli MG1655 under octanoic acid (C8) stress. Appl Microbiol Biotechnol. 2015;99:4397–408.
Rodriguez S, Denby CM, Van Vu T, Baidoo EEK, Wang G, Keasling JD. ATP citrate lyase mediated cytosolic acetyl-CoA biosynthesis increases mevalonate production in Saccharomyces cerevisiae. Microb Cell Fact. 2016;15:48.
You L, Page L, Feng X, Berla B, Pakrasi HB, Tang YJ. Metabolic pathway confirmation and discovery through 13C-labeling of proteinogenic amino acids. J Vis Exp. 2012;59:e3583.
Wahl SA, Dauner M, Wiechert W. New tools for mass isotopomer data evaluation in 13C flux analysis: mass isotope correction, data consistency checking, and precursor relationships. Biotechnol Bioeng. 2004;85:259–68.
Borodina I, Schöller C, Eliasson A, Nielsen J. Metabolic network analysis of Streptomyces tenebrarius, a streptomyces species with an active entner–doudoroff pathway. Appl Environ Microbiol. 2005;71:2294–302.
Bennett BD, Yuan J, Kimball EH, Rabinowitz JD. Absolute quantitation of intracellular metabolite concentrations by an isotope ratio-based approach. Nat Protoc. 2008;3:1299–311.
Orth JD, Conrad TM, Na J, Lerman JA, Nam H, Feist AM, Palsson BØ. A comprehensive genome-scale reconstruction of Escherichia coli metabolism—2011. Mol Syst Biol. 2011;7(1):535.
AM, WDH, and YJT designed experiments. WDH, EEB, GW, and SR conducted experiments. All authors analyzed the data, wrote the manuscript. All authors read and approved the final manuscript.
We would like to thank Victor Chubukov, Mary Abernathy, and Tolutola Oyetunde for their assistance with paper. The authors are especially thankful for useful discussion about channeling with Prof. Daniel H. Kohl at Washington University.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files] and all strains are available through JBEI ICE Registry.
The authors declare that they have no competing interests. JDK has a financial interest in Amyris and Lygos.
This work was supported by NSF (DBI 1356669) and DOE-BER (DESC0012722) given to YJT and by the Joint BioEnergy Institute (http://www.jbei.org/), which is supported by the US Department of Energy, Office of Science, Office of Biological and Environmental Research, through contract DE-AC02-05CH11231 between Lawrence Berkeley National Laboratory and the US Department of Energy and DOE Office of Science Graduate Fellowship (SCGF2015) Program to WDH. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes.