3.1 Strain construction to eliminate and replace PPi-dependent glycolytic reactions
In our previous work, we created strain LL1660 (Fig. 1), where the PPi-pfk gene was replaced with an ATP-pfk gene [17], eliminating a major sink for PPi. We also observed the spontaneous occurrence of a large partial genome duplication. This partial genome duplication region encodes proteins that may serve as alternative sinks for PPi, including a membrane-bound pyrophosphatase gene (HppA - Clo1313_0823) and pyruvate:phosphate dikinase gene (PpdK - Clo1313_0949).
We therefore hypothesized that the partial genome duplication event might have occurred to address an imbalance in the supply and demand of PPi created by the mutations in strain LL1660. To address this potential problem, we set out to construct a new strain where PPi-consuming reactions in glycolysis were eliminated, and where the excess PPi generated by metabolism was hydrolyzed using a soluble PPase [25]. Initially, we attempted to express soluble PPase in wild-type C. thermocellum, but we did not get any colonies in at least two transformations. This was not surprising, given our prior work showing the importance of PPi as a cofactor for the PFK reaction in glycolysis [11, 15].
Next, we tried to express the soluble PPase in a strain with ATP-linked PFK activity (strain LL1649), which we hypothesized would be more compatible with the metabolic effects of the PPase. In this strain, we occasionally obtained a small number of colonies with our PPase expression plasmid. However, in those colonies, the PPase gene was invariably inactivated either by a transposon insertion or inactivating nonsense mutation. In this phase, we tested PPases from Geobacillus stearothermophilus, Geobacillus thermoglucosidasius, Thermoplasma acidophilum, Thermus aquaticus, and Thermus thermophilus (plasmids pLL1501 to pLL1505, Table 2). Out of these, only PPase from Geobacillus thermoglucosidasius and Thermus aquaticus showed PPase activity (data not shown).
We further hypothesized that in addition to eliminating the need for PPi in the PFK reaction, we also needed to eliminate the need for PPi in the PPDK reaction (which converts PEP to pyruvate) [15, 34]. To do this, we designed an integration construct to simultaneously introduce pyruvate kinase and a soluble PPase onto the chromosome, in this case using the PPases from G. thermoglucosidasius and T. aquaticus since we had in some instances observed PPase activity from cell extracts of C thermocellum that expressed those PPases [25]. Successful strain construction was obtained only with T. aquaticus PPase containing construct, pLL1507 (Table 2), which was validated by the gain of both pyruvate kinase and cytosolic PPase activity (Fig. 2). Our ability to observe functional PPase expression only in strains with both ATP-PFK and PYK activity provides further evidence of progress toward eliminating PPi as a key energy carrying cofactor in C. thermocellum.
Strain LL1592, which uses the native C. thermocellum glycolytic pathway, with T. saccharolyticum pyruvate-to-ethanol pathway [17] was selected as reference strain for engineering. Strain LL1649 was derived from strain LL1592, by introduction of the T. saccharolyticum transaldolase, followed by simultaneous deletion of the ADP-glucose synthase operon and expression of T. saccharolyticum ATP-pfk. Strain LL1689 was derived from strain LL1649 by the simultaneous integration of the T. saccharolyticum pyruvate kinase and the cytosolic pyrophosphatase from Thermus aquaticus. The purpose of introducing a heterologous pyruvate kinase was to allow subsequent deletion of the ppdk gene. The purpose of introducing a soluble PPase was to serve as an alternative sink for PPi, once both the PPi-pfk and ppdk genes were deleted. With these modification in place, we were then able to successfully delete glycolytic sinks for PPi, including the ppdk gene (resulting in strain LL1710), and then, subsequently, the PPi-pfk gene (resulting in strain LL1711) (Fig. 1).
Biochemical assays were performed to confirm the effect of the genetic modifications of the strains used in this study. As expected, strain LL1592 and strain LL1711 exclusively have PPi-PFK and ATP-PFK activity, respectively. Strains LL1649, LL1689, and LL1710 have both ATP- and PPi-linked PFK activity, as expected (PFK reaction, Fig. 2). Functional expression of the pyk gene was confirmed by the presence of PYK activity in strains LL1689, LL1710 and LL1711 (PYK reaction, Fig. 2). The effect of the ppdk gene deletion was confirmed in strains LL1710 and LL1711 based on elimination of PPDK activity (PPDK reaction, Fig. 2). Low levels of PPase activity were still detected in strains LL1592 and LL1649 (Paired t-test, p > 0.5); this could be attributed to the native membrane bound PPase enzyme being present in the cell extracts. Previous studies have reported the presence of membrane-bound pyrophosphatase in C. thermocellum cells [25]. Nonetheless, about 5- to 13-fold higher PPase activities (Paired t-test, p < 0.005) were observed in strains LL1689, LL1710 and LL1711, which we attribute to expression of a soluble PPase. As such, it was determined that strain LL1711 had been engineered to not require PPi as a cofactor in glycolysis, and that it now possessed a more canonical glycolytic pathway [12].
3.2 Use of ATP-dependent pathway alters thermodynamics of glycolysis
To determine the effect the genetic modifications had on the reversibility of the PFK reaction, we performed 13C labeling on C. thermocellum strains LL1590, LL1592 and LL1711. Initially, we were planning to compare only the parent strain (LL1592) and the PPi-free glycolysis strain (LL1711); however, whole genome sequencing of strain LL1711 revealed a frameshift mutation in the wild type C. thermocellum adhE gene corresponding to amino acid position E784 that eliminated its function (Figure S1). We therefore included a control strain (LL1590) derived from the parent strain with a targeted deletion of the wild type C. thermocellum adhE gene [18], to better understand the effects of our changes to the glycolytic pathway without the confounding effects of the modification to ethanol production pathway. Note that all three strains (LL1590, LL1592, and LL1711) still produce ethanol via the heterologous T. saccharolyticum adhEG544D (Fig. 1).
C. thermocellum does not readily consume glucose [35], and thus cellobiose is the most commonly-used soluble sugar for laboratory growth experiments. Strains were cultured on a 1:1 ratio of naturally labeled to uniformly 13C-labeled cellobiose. We confirmed that 13C cellobiose and naturally labeled cellobiose are taken up equally by measuring the ratio of M + 0 to M + 6 isotopomers for glucose 6-phosphate, which ranged from 0.8–1.2 (Additional File 2). Samples were collected at three different phases (early-log, late-log and stationary) during growth (Figure S2). As previously described [12], the forward glycolysis reaction should generate a 1:1 mixture of fully labeled (M + 6) and fully unlabeled (M + 0) glucose-6-phosphate (G6P), fructose-6-phosphate (F6P) and fructose-1,6-bisphosphate (FBP). The forward flux of fructose-bisphosphate aldolase reaction (FBA) will generate 50% unlabeled (M + 0) and 50% fully labeled (M + 3) dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (G3P). However, the reverse FBA reaction flux generates a 1:2:1 mixture of M + 0, M + 3 and M + 6 FBP. The reverse flux from phosphofructokinase (PFK) and phosphoglycerate isomerase (PGI) reactions transfer these M + 3 species into the F6P and G6P pools. Appearance of the M + 3 isotopomer in the F6P and G6P pools allows us to uniquely identify reverse flux through the PFK reaction. If the PFK reaction were completely irreversible, we would expect to see no M + 3 species in either the F6P or G6P pools.
As expected, based on our prior work [17], we observed that the M + 3 pool in the glycolytic intermediates G6P and F6P decreased substantially in the PPi-free glycolysis strain (LL1711) vs. the parent and control strains (LL1590 and LL1592), indicating that the genetic modifications introduced into LL1711 reduced the reversibility of the PFK reaction (G6P, F6P; Fig. 3). We saw an unexpected decrease in the M + 3 fraction of FBP in mid-log phase in the PPi-free glycolysis strain (LL1711) (Fig. 3, Figure S3). This suggests reduced reversibility of Fba reaction, although the reason for this change is not known.
In addition to comparing the isotope ratios, we also compared the ratio of the FBP and F6P metabolite pool sizes. This gives us an additional indication of the thermodynamic driving force of the reaction. A high thermodynamic driving force at the PFK reaction would result in a high FBP/F6P ratio. Comparisons of FBP/F6P ratio showed a lower ratio for the parent and control strains that use the PPi-PFK reaction (LL1590 and LL1592) and a much higher ratio for the PPi-free glycolysis strain (LL1711), which uses the ATP-PFK reaction (Fig. 3). This ratio is 5- to 17-fold and 74- to 120-fold higher in ATP-pfk strain compared to PPi-pfk strains during log phases and stationary phases, respectively. The higher FBP/F6P ratio is similar to what is observed in organisms with canonical glycolytic pathways, such as E. coli and T. saccharolyticum [12, 36]. Taken together, these two lines of evidence (M + 3 abundance changes and FBP/F6P ratio changes) support the conclusion that we have increased the thermodynamic driving force of the PFK reaction in our PPi-free glycolysis strain.
3.3 Expression of PPase results in lower PPi levels
To determine the effect of PPase expression in PPi levels of C. thermocellum strains with wild type glycolysis (strains LL1590 and LL1592) and PPi-free glycolysis (strain LL1711), we quantified the concentration of PPi in cells.
We observed a significant decrease in PPi concentration in the PPi-free glycolysis strain (strain LL1711) (Fig. 4) during the early-log and late-log growth phases. During the stationary phase, the difference between strain LL1592 and strain LL1711 isn’t as significant, which might be due to the depletion of the PPi pool during this phase. We see consistently higher levels of PPi during the early-log and late-log phases in strains LL1590 and LL1592, which then drop during the stationary phase in both strains. The constant low level of PPi in strain LL1711 during all the growth phases might be the basal level of PPi in cells, even with the expression of PPase.
3.4 Increased driving force leads to higher level of fermentation products in ATP-pfk strain
Having demonstrated increased thermodynamic driving force in glycolysis in the PPi-free glycolysis strain, we explored the impact of fermentation behavior and ethanol production using batch fermentations of all three strains with high substrate concentrations (100 g/L cellobiose). The fermentations had three distinct phases: 1. growth coupled fermentation, where both the cell growth and ethanol production are at the highest rate, 2. growth uncoupled fermentation, where cell growth stops and the rate of ethanol production decreases, and 3. cessation of both growth and ethanol production. During this final phase, cellobiose is hydrolyzed to glucose at a slow rate, but metabolism appears to be inactive.
The final ethanol titer was 329 ± 8 mM (15.1 ± 0.4 g/L) for strain LL1590, 403 ± 11 mM (18.6 ± 0.5 g/L) for strain LL1592, and 455 ± 12 mM (21.0 ± 0.6 g/L) for strain LL1711. Titer of another fermentation product, acetate, was 65 ± 0 mM (3.8 ± 0 g/L) for strain LL1590, 51 ± 9 mM (3.0 ± 0.5 g/L) for strain LL1592 and 88 ± 7 (5.2 ± 0.4 g/L) for strain LL1711. (Fig. 5, Figure S4, Additional File 3). Comparison of the fermentation products between LL1592 and LL1590 suggests that deletion of C. thermocellum adhE impairs the production of ethanol. Taking the effect of the native adhE inactivation into account (i.e. comparing LL1711 to LL1590), allows us to observe the effect of PPi-free glycolysis on ethanol production without the confounding effects of changes in the ethanol production pathway. In this comparison (LL1711 vs LL1590), we observe a 38% increase in ethanol titer, and a 35.4% increase in acetate titer. Together, these suggest that our PPi-free glycolysis strain exhibits increased glycolytic flux, resulting in an increased abundance of acetyl-CoA, which is then converted to more ethanol and acetate.
3.5 Changes in glycolysis metabolite levels at elevated substrate concentrations in pH-controlled fermentations
To better understand the changes in central metabolism that allowed for increased ethanol titer in the PPi-free glycolysis strain, we measured intracellular glycolytic metabolites over the course of the fermentation using the same fermentation conditions from Fig. 5 (100 g/L cellobiose, pH-controlled batch fermentations). Since preliminary data showed similar pattern of glycolytic metabolites (Figure S5) and 13C labeling patterns were similar for both the PPi-pfk strains (Fig. 3), and strain LL1590 produced lower levels of ethanol than strain LL1592 (Fig. 5), we compared glycolytic metabolite levels in strains LL1592 and LL1711 .
Three major changes were observed in the glycolytic metabolite profiles between these two strains. First, in the parent strain (LL1592, with wild type glycolysis) we see gradual accumulation of hexose phosphates (G6P, F6P) as the fermentation progresses. In the PPi-free glycolysis strain (LL1711), hexose phosphates were still observed to accumulate, but to a lesser extent as the fermentation progressed. Second, during the growth-uncoupled fermentation phase, we observed a dramatic increase in the levels of lower glycolysis metabolites (DHAP, 3PG, and PEP) in the PPi-free glycolysis strain. Finally, in the phase where fermentation and ethanol production stop, the parent strain shows a large accumulation of upper glycolysis metabolites (glucose, G6P, and F6P), while the PPi-free glycolysis strain shows low levels of almost all metabolites (Fig. 6, Figure S6).
These intracellular metabolites data over the course of fermentation also underscores the importance of conducting metabolomics studies throughout the fermentation process. Most previous metabolomics studies in C. thermocellum and other microbes have focused on metabolite levels at a single-time point, typically during the mid-exponential growth phase [12, 36, 37]. However, analyzing metabolites at multiple time points during fermentation reveals important dynamic changes in the metabolite pools.