Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2008, p. 3634-3643, Vol. 74, No. 12
0099-2240/08/$08.00+0 doi:10.1128/AEM.02708-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Department of Chemical Engineering and Materials Science,1 BioTechnology Institute, University of Minnesota, 240 Gortner Laboratory, 1479 Gortner Avenue, St. Paul, Minnesota 551082
Received 30 November 2007/ Accepted 8 April 2008
|
|
|---|
|
|
|---|
A metabolic network can function according to many different pathway options. Such redundancy of pathways enables cells to compete efficiently and to survive under changing environmental conditions (39). Elementary mode (EM) analysis has emerged as a powerful systems biological tool that rigorously dissects a metabolic network into its basic building blocks (35, 36). The metabolism of a functioning cell has to be viewed as a weighted average of the fluxes through all fundamental pathways (EMs) that its metabolic network supports (40, 43). These EMs therefore represent the inherent building blocks of the metabolic structure. The set of EMs is in fact the parts list for cell function encoded at a higher level in the hierarchy of biological complexity. The quest for the minimal cell can be conducted, therefore, at this functional level. DNA replication and protein synthesis appear to be only supporting functions that ensure the self renewal of cell functionality.
Ethanol has emerged as an important renewable and sustainable energy source that can reduce our reliance on fossil resources. It can be produced from inexpensive, abundant, and renewable feedstocks, including cellulosic and lignocellulosic biomass, by fermentation using microorganisms, and significant advances have already been made (6, 17, 20, 30, 44, 45). The key to this technology is developing efficient and robust microorganisms to convert the biomass-derived hexoses and pentoses to ethanol at the best possible yields and high productivities (33).
To construct a minimal cell that is dedicated to produce ethanol in the most efficient way, we took a top-down approach and started with a complex, functioning Escherichia coli cell that has at its disposition all the numerous metabolic pathway options characteristic of a wild-type cell. We chose E. coli as a model organism to demonstrate this approach because it can degrade a variety of pentoses and hexoses and its genetics can be easily engineered according to rational strain design with well-established molecular techniques.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. List of strains and plasmids
|
Growth in batch bioreactors.
Batch bioreactor experiments were conducted in 10-liter Braun bioreactors (Biostat MD; B. Braun Biotech International, Melsungen, Germany) with a working volume of 6 liters under anaerobic conditions. The temperature and agitation rate were set at 37°C and 200 rpm, respectively. Single colonies were picked from freshly streaked plates and grown overnight in 15-ml tubes containing 5 ml of rich medium. The cultures were then transferred to 250-ml capped shake flasks containing 100 ml of rich medium. Exponential cultures grown in shake flasks (37°C, 225 rpm) were then used for inoculation. The media used for inoculation and for the bioreactors were identical. The initial optical density measured at a 600-nm wavelength (OD600) after inoculation in all batch bioreactors was 0.05. To maintain anaerobic growth conditions, nitrogen was sparged into bioreactors through a 0.2-µm filter at a volumetric flow rate of 100 ml/min at least 4 h before inoculation and throughout the fermentation. The exhaust gas was first passed through an exhaust gas condenser and then a 0.2-µm filter and a pressure regulator and finally into a Prima
-B mass spectrometer (ThermoOnix, Houston, TX) to analyze gas composition. The reactor gauge pressure was set at 1 lb/in2 to minimize air diffusion into the bioreactors and hence maintain anaerobic growth conditions. pH was controlled at 6.5 by using 6 M NaOH and 40% H3PO4. The anaerobic growth conditions could be confirmed by the absence of oxygen signals from mass spectroscopy. Fermentation was completed when H3PO4 started being added to the bioreactors.
For growth conducted in baffled shake flasks, the procedures used were as previously reported (41). For growth in anaerobic shake tubes, anaerobic 30-ml glass tubes were filled with 20 ml of defined medium, sparged with N2, and sealed with aluminum crimples on a rubber stopper. The initial OD600 after inoculation was set at 0.01. Measurement of the optical density of the anaerobic tubes was carried out directly by using the Spectronic 20D+ (model no. 333182000; Houston, TX) that has a tube adapter.
Analytical techniques.
The optical densities of the cultures were measured at a wavelength of 600 nm in 1-cm cuvettes using a Hewlett Packard 8452A diode array spectrophotometer (Palo Alto, CA). Ten milliliters of a culture was withdrawn periodically from a bioreactor and immediately processed to determine the cell dry weight and the levels of secreted metabolites in the supernatant. First, the sample was spun at 7,000 x g at 4°C for 10 min. Then its supernatant was stored at –20°C for later analysis, and the cell pellet was washed once with deionized water, vacuum filtered, and weighed in a weighing dish after being dried in a 65°C oven for at least 1 day. The formula for weight conversion of optical density is 1 OD600 = 0.259 g/liter (R2 = 0.942). Metabolite concentrations were determined by using a high-pressure liquid chromatography system (Shimadz10A; Shimadzu, Columbia, MD) equipped with an autosampler (SIL-10AF), a cation exchange column (HPX-87H; Bio-Rad, Hercules, CA), and two detectors in series consisting of a UV-visual spectroscopy detector (SPD-10A) and a refractive index detector (RID-10A). Samples from cell supernatants were first filtered through a 0.22-µm filter unit. Then the samples were loaded into the column operated at 65°C. A 5-mM H2SO4 solution was used as the mobile phase and run isocratically at a flow rate of 0.5 ml/min.
Yield calculation and carbon balance.
Ethanol yields on sugars (YETOH/sugars) were determined by the formula YETOH/sugars = rETOH/rsugars (g of ethanol/g of sugars), where rETOH (g of ethanol/liter/h) and rsugars (g of sugars/liter/h) represent the ethanol production rate and sugar consumption rate, respectively. In all experiments, the YETOH/sugars appeared to be constant, since the linear regression of ethanol produced (g/liter) and sugars consumed (g/liter) yielded a perfect fit with R2 > 0.99. Yields of other by-products on sugars were computed in the same way.
The percent of carbon recovery was calculated by the following equation:
), where qX, qSuc, qLac, qAce, qETOH, qFor, and
(Cmol of product/Cmol of sugars/h) are the specific rates of biomass, succinic acid, lactic acid, acetic acid, ethanol, formic acid, and carbon dioxide, respectively. The specific rates of CO2 were not directly measured but were estimated from the formation of ethanol, succinic acid, and acetic acid. Specifically, the formation of ethanol and acetic acid leads to gains of equimolar amounts of CO2, but the formation of succinic acid under anaerobic growth conditions leads to losses of equimolar amounts of CO2. The carbon recovery rate was close to 100%, as shown in Table S2 in the supplemental material.
E. coli metabolic network.
A metabolic network was constructed for E. coli that can grow on pentoses and hexoses, including D-(+)-xylose, L-(+)-arabinose, D-(+)-glucose, D-(+)-galactose, and D-(+)-mannose, by using available public databases (22) and published literature (8, 28), as shown in Fig. 1 and Table S3 in the supplemental material. We considered that the transport of glucose and mannose in the model was facilitated by the phosphoenolpyruvate sugar transferase system (PTS), while the uptake of galactose, xylose, and arabinose was mediated by ABC transporters with high affinity (26). The constructed model represents the intermediary metabolism of E. coli. The validity of the model was experimentally tested in previous studies (41, 43). The pyruvate decarboxylase reaction that converts pyruvate to acetaldehyde was also considered in constructing the model. The pyruvate decarboxylase does not exist in E. coli but was introduced into E. coli through the plasmid pLOI297 (ATCC 68239) (3). The introduction of pyruvate decarboxylase into the strain contributed to the use of the precursor pyruvate for the ethanol production pathway and mimicked the ethanol-producing pathway of native ethanologenic strains such as Zymomonas mobilis and Saccharomyces cerevisiae. It has been reported that pyruvate carboxylase has a higher affinity to pyruvate than pyruvate formate lyase, which converts pyruvate to acetyl coenzyme A (CoA) (19). Furthermore, the conversion of pyruvate into acetaldehyde is not dependent on the cofactor CoA, which can be highly regulated and become limiting. Some reactions in sugar degradation pathways that occur in series without branches have been lumped to simplify the model without affecting the analysis. We calculated all EMs of the E. coli metabolic network using METATOOL 5.0, the currently available, fast, and flexible Matlab-based software package designed to handle complex metabolic networks (42).
![]() View larger version (23K): [in a new window] |
FIG. 1. Metabolic map of E. coli central metabolic network. Deleted reactions in TCS083 are shown next to the symbol X. Listed are glucose-6-phosphate-1-dehydrogenase (PPP1, zwf), NADH dehydrogenase II (OPM4r, ndh), NAD/NADP-dependent malate enzyme (ANA2, sfcA/maeB), D-lactate dehydrogenase (FEM3, ldhA), fumarate reductase (TCA10, frdA), pyruvate oxidase (FEM2, poxB), and phosphate acetyltransferase (FEM7, pta). ETOH, ethanol.
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 2. EMs that utilize different pentoses and hexoses as carbon sources and ranges of ethanol and biomass yields
|
![]() View larger version (45K): [in a new window] |
FIG. 2. Effect of reaction deletions on the numbers of anaerobic EMs for growth on xylose (A and B) and glucose (C and D). The bars in panels A and C specify the numbers of EMs for strains with deletions of the indicated genes. In each group of bars, the numbers of (i) total modes, (ii) modes that make ethanol, (iii) modes that produce biomass, and (iv) modes that make both biomass and ethanol are listed. The possible maximal and minimal ethanol and biomass yields for xylose (B) and glucose (D) are shown. Note that the minimal yields are pushed toward the upper theoretical limit with increasing numbers of deleted genes. ETOH, ethanol.
|
|
View this table: [in a new window] |
TABLE 3. Stoichiometric equations for the efficient ethanol-producing EMs that use pentoses and hexosesa
|
Coutilization of pentose and hexose.
The coutilization of xylose and glucose poses an interesting situation, since the metabolism of individual sugars exhibits opposite flux distributions in parts of the metabolic network. The EM analysis identified 92,593 EMs, which are considerably more than the sum of the total number of EMs that support utilization of each sugar alone (Table 2). The new EMs appear due to the coutilization of both glucose and xylose. Deletion of the same set of 7 reactions resulted in 18 remaining anaerobic EMs; 12 utilize either glucose or xylose individually as presented above, and the other 6 coutilize glucose and xylose. Deletion of this set of reactions effectively eliminated all lower yielding pathways and pushed the range of ethanol yields toward the upper limit (0.36 to 0.51 g ethanol/g sugars) (Table 2).
The set of EMs possible on a mixture of sugars consisting of glucose and pentoses can be further reduced to six by creating a strain unable to use glucose. Such a pentose-specific strain can be realized by removing, in addition to the previously described genes, the glucose phosphotransferase system (ptsG), glucose kinase (glk), and the mannose phosphotransferase system (manX), thus preventing glucose transport into the cell or glucose phosphorylation in the cell.
Strain comparisons.
Over the past 15 to 20 years, after the expression of the foreign ethanol-producing pathway from Zymomonas mobilis in E coli was successful (3), several ethanologenic E. coli strains have been developed to improve ethanol production through numerous rounds of modification. The modifications typically rely on intuitive understanding of cell metabolism and cell behavior (10, 11, 24, 30, 44). Unlike these approaches, our approach was rational, based on EM analysis to design cells with a minimal functionality dedicated to ethanol production. To demonstrate differences between the rationally designed strain TCS083 and other strains developed previously via intuition, we applied EM analysis to describe the effect of the gene mutations in these strains on the range of ethanol yields and on the reduction of inefficient ethanol-producing pathways (Table 4). The results show that the developed strains FBR3 (11), FBR5 (10), KO11 (LYO1) (30, 44), and LY168 (21) still contain a large portion of inefficient ethanol-producing EMs that support a large range of ethanol yields. The existence of inefficient ethanol-producing pathways in these strains can potentially reduce the ethanol yield. For instance, the characterization of KO11 in chemostat studies showed that during growth on xylose, KO11 loses its hyperethanologenicity at the expense of cell growth and acetate synthesis (12). Compared to each of these engineered strains, TCS083 contains a larger number and a different set of deleted reactions. It should be mentioned that unlike TCS083, all of these engineered strains with the exception of LY168 have pflB deleted. This gene belongs to the operon of pflABCD, which encodes pyruvate formate lyase that is active under anaerobic growth conditions (1, 46). According to EM analysis, deletion of this enzyme prevents growth under anaerobic conditions because the synthesis of acetyl-CoA, required for biomass synthesis, is blocked. For this reason, TCS083 does not contain a pflB deletion. However, the normal growth of FBR3, FBR5, and KO11 (LY01) likely still occurs, presumably due to the activity of other subunits of pyruvate formate lyase and the subsequent metabolic evolution of engineered strains (11, 30).
|
View this table: [in a new window] |
TABLE 4. Elementary mode analysis of various ethanologenic E. coli strains using glucose as a carbon source under anaerobic conditionsa
|
![]() View larger version (21K): [in a new window] |
FIG. 3. PCR of deleted genes for TCS083 and CT1101, with the wild type used as a positive control. Panels A and C show deleted genes tested for TCS083, including zwf, ndh, sfcA, maeB, ldhA, frdA, poxB, and pta. Panels B and D display deleted genes tested for CT1101, including zwf, ndh, sfcA, maeB, ldhA, frdA, poxB, pta, glk, manX, and ptsG. For each gene tested, the left lane shows the location of an amplified gene for the wild type and the right lane for the mutant. A shift to a smaller band size occurring in a lane of a mutant indicates that the tested gene is deleted. The arrow in each lane points to the location of the expected band size.
|
We tested the performance of strain TCS083/pLOI297 together with MG1655/pLOI297 as a control, first on individual xylose and glucose sugars and then on sugar mixtures in controlled batch bioreactors. Figure 4A and B shows the batch reactor time profiles of the wild-type MG1655/pLOI297 and TCS083/pLOI297 for growth on 80 g/liter glucose. MG1655/pLOI297 achieved an ethanol yield of 0.46 ± 0.01 (g ethanol/g glucose) and an ethanol titer of 36.53 ± 0.31 g/liter. Under identical growth conditions, TCS083/pLOI297 reached a cell dry weight of 3.11 ± 0.47 g/liter. Because the stoichiometries of the two biomass-producing pathways and the four remaining pathways that produce only ethanol are known (Table 3), we could compute on the basis of biomass formed that 16% ± 2% of the total glucose was consumed in the pathway supporting growth and 84% ± 2% in the pathways producing only ethanol. Furthermore, with these weighting factors together with the pathway stoichiometries, we could predict that 80 g/liter of consumed glucose should result in a final ethanol titer of 39.20 ± 0.34 g/liter. This value is in excellent agreement with the measured ethanol titer of 38.77 ± 0.63 (g/liter) (Table 5).
![]() View larger version (30K): [in a new window] |
FIG. 4. Time profiles for glucose, xylose, cell dry weight (cdw), and ethanol for the wild type (upper three panels) and mutant TCS083/pLOI297 (lower three panels). The strains were anaerobically cultivated in controlled 10-liter bioreactors where they were sparged with nitrogen. The initial sugar concentration was 80 g/liter. In an experiment with mixed sugars, 40 g/liter of each sugar was provided. Note that the figure shown is representative of a single-batch bioreactor run of duplicate experiments and that each data point represents the mean ± standard deviation of three independent measurements.
|
|
View this table: [in a new window] |
TABLE 5. Prediction of ethanol titers of TCS083/pLOI297 and CT1101/pLOI297a
|
Determination of the ethanol production kinetics on a mixture of 40 g/liter glucose and 40 g/liter xylose showed that MG1655/pLOI297 consumed first glucose and then xylose in a sequential manner. As shown in Fig. 4E, it took MG1655/pLOI297 only about 12 h to completely consume glucose but 72 h to consume just 20% of 40 g/liter xylose. The slow xylose consumption rate of the wild type after the glucose was consumed was probably due to a combination of cellular inhibitory effects resulting from the formation of by-products such as succinic acid, lactic acid, acetic acid, ethanol, and formic acid. MG1655/pLOI297 achieved an ethanol yield of 0.40 ± 0.02 (g ethanol/g sugar) and an ethanol titer of 18.66 ± 1.01 (g/liter). In contrast, TCS083/pLOI297 did not show pronounced diauxic growth behavior. It could simultaneously consume both glucose and xylose and completely utilize all sugars within 48 h. TCS083/pLOI297 reached a biomass concentration of 3.54 ± 0.16 g/liter (Fig. 4F). Assuming that both sugars equally contributed to the biomass formation reaction, this translates to the use of 32% ± 1% of the sugars for growth and the conversion of the rest into ethanol at theoretical values to reach a final titer of 39.46 ± 0.57 g/liter. The experimentally determined final ethanol titer was 38.81 ± 0.91 (g/liter) (Table 5).
Catabolite repression.
Fig. 5C shows the consumption of only pentoses by CT1101/pLOI297 in a mixture of xylose and glucose. This strain evidently has all catabolite repression removed and can grow on pentoses even in the presence of glucose. The strain reached a biomass concentration of 1.48 ± 0.02 g/liter. According to pathway stoichiometries, 39 ± 1% of the total xylose consumed was used for growth. The computed final ethanol titer was 19.32 ± 0.73 g/liter and compared well with the experimental value of 20.03 ± 0.18 (g/liter) (Table 5).
![]() View larger version (34K): [in a new window] |
FIG. 5. Time profiles for glucose, xylose, cell dry weight (cdw), and ethanol for mutant CT1101/pLOI297 growing on a sugar mixture (lower panel) (see the legend to Fig. 4 for growth conditions). The relationship between consumed xylose and consumed glucose is shown for anaerobic growth (upper left panel) and aerobic growth (upper right panel) for the wild-type MG1655 and mutant strains TCS083 and CT1101 and a 1:4 initial mixture of strains TCS083 and CT1101.
|
|
|
|---|
The theoretical framework for designing a minimal cell is based on EM analysis. It can dissect a metabolic network into unique, nondecomposable pathways that present all possible physiological states of cells at steady-state conditions (35, 36). A different approach based on flux balance analysis using optimization frameworks such as OptKnock and MOMA has also been applied by other research groups to identify gene knockouts to optimize both growth and production of a metabolite (2, 7, 13, 37). This approach does not identify the complete set of alternative optimal solutions nor other subsets of suboptimal solutions. Therefore, even though the set of deleted genes identified by flux balance analysis can identify an optimal pathway, there is no guarantee that the engineered cells can function according to this optimal pathway.
A remarkable phenotype of TCS083/pLOI297 is the ability to utilize both xylose and glucose simultaneously (Fig. 4F and 5A). This phenotype clearly does not exist in the wild type (Fig. 4E, 5A), which shows pronounced diauxic growth behavior presumably mediated through glucose catabolite repression (25). It was previously shown that in E. coli, deletion of the ptsG gene can remove catabolite repression at the expense of a lower glucose uptake rate (29). In our minimal strains, ptsG is still present, but the specific growth rate on glucose is reduced, indicating that catabolite repression is closely linked to the growth dynamics, and likely to the glycolytic flux, that determine the concentrations and the effects of specific regulation factors in the strains. This finding is consistent with experiments that have shown that blockage of the glycolytic flux triggers the degradation of ptsG mRNA (23). Moreover, when it grows at a high specific growth rate, our minimal strain does show pronounced catabolite repression and preferential glucose consumption under aerobic growth conditions comparable to those of the wild-type strain (Fig. 5B).
TCS083/pLOI297 was rationally designed for efficient ethanol production based on EM analysis. The performance of TCS083/pLOI297 was tested experimentally. The mutant outperformed the wild type in the efficient conversion of sugars into ethanol and closely matched the theoretical prediction. Based on published data, we compared the mutant TCS083/pLOI297 with other engineered ethanologenic E. coli strains such as KO11 (44) and FBR5/pLOI297 (29) under similar growth conditions. For a similar period of xylose fermentation, both KO11 and FBR5/pLOI297 achieved a similar ethanol yield of 0.46 ± 0.01 (g/g), for a range of 75 to 90 g/liter xylose consumed, while TCS083/pLOI297 reached a higher ethanol yield of 0.49 ± 0.01 (g/g). In addition, FBR5/pLOI297 obtained an ethanol yield of 0.46 ± 0.03 (g/g) for a similar period of glucose fermentation, while TCS083/pLOI297 achieved a higher ethanol yield of 0.49 ± 0.01 (g/g). Furthermore, TCS083/pLOI297 can simultaneously consume pentoses and hexoses, while neither KO11 nor FBR5 can. Overall, TCS083/pLOI297 appears to perform more favorably than KO11 and FBR5/pLOI297 under similar growth conditions. The performance of TCS083/pLOI297 will be further characterized and compared with those of KO11 and FBR5/pLOI297 under different growth conditions.
TCS083 is expected to be useful for other biotechnological applications. Under completely aerobic growth conditions, the biomass yields of TCS083 are maximized, similar to those of TCS062, since CO2 is the only by-product and its formation is minimized. In large-scale reactor operations, oxygen-limiting growth conditions are typically created due to mixing inhomogeneities, resulting in the secretion of acetic acid. In such cases, TCS083 should perform better than TCS062, since acetate-producing pathways are disrupted. Indeed, shake flask experiments, in which completely aerobic conditions cannot be achieved, have already confirmed that TCS083 does indeed outperform not only TCS062 and the wild-type strain but also the industrial strain BL21 (38) and the reduced-genome strain MDS42 (32) typically used for protein production (Fig. 6).
![]() View larger version (37K): [in a new window] |
FIG. 6. Growth characteristics of different E. coli strains, including BL21 MDS42, MG1655, TCS062, and TCS083. (A) Specific growth rates. (B) Biomass yields on glucose. (C) Acetate yields on glucose. The experiments were conducted in baffled shake flasks containing defined medium supplied with various pentoses and hexoses under aerobic conditions. Each value represents the mean ± standard deviation of the results of triplicate experiments.
|
The design of a strain with minimal metabolic functionality must always be coordinated with the purpose of the functionality. Here, the purpose was the most efficient production of ethanol. However, the methodology applied is general and should prove useful for the design and construction of many different minimal strains tailored for applications in biotechnology and for fundamental studies. Such strains are attractive due to the efficiency and simplicity of their metabolic pathways.
We thank the Minnesota Supercomputing Institute for use of their supercomputing facilities.
Published ahead of print on 18 April 2008. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
DE3) and Escherichia coli JM109. Biotechnol. Bioeng. 49:421-428.[Medline]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»