Directed evolution of a cellobiose utilization pathway in Saccharomyces cerevisiae by simultaneously engineering multiple proteins
© Eriksen et al.; licensee BioMed Central Ltd. 2013
Received: 27 January 2013
Accepted: 3 June 2013
Published: 26 June 2013
The optimization of metabolic pathways is critical for efficient and economical production of biofuels and specialty chemicals. One such significant pathway is the cellobiose utilization pathway, identified as a promising route in biomass utilization. Here we describe the optimization of cellobiose consumption and ethanol productivity by simultaneously engineering both proteins of the pathway, the β-glucosidase (gh1-1) and the cellodextrin transporter (cdt-1), in an example of pathway engineering through directed evolution.
The improved pathway was assessed based on the strain specific growth rate on cellobiose, with the final mutant exhibiting a 47% increase over the wild-type pathway. Metabolite analysis of the engineered pathway identified a 49% increase in cellobiose consumption (1.78 to 2.65 g cellobiose/(L · h)) and a 64% increase in ethanol productivity (0.611 to 1.00 g ethanol/(L · h)).
By simultaneously engineering multiple proteins in the pathway, cellobiose utilization in S. cerevisiae was improved. This optimization can be generally applied to other metabolic pathways, provided a selection/screening method is available for the desired phenotype. The improved in vivo cellobiose utilization demonstrated here could help to decrease the in vitro enzyme load in biomass pretreatment, ultimately contributing to a reduction in the high cost of biofuel production.
KeywordsCellobiose utilization β-glucosidase Cellodextrin transporter Directed evolution Protein engineering Pathway engineering Pathway optimization Pathway libraries
The saccharification of lignocellulosic biomass into fermentable sugars is recognized as one of the most expensive operations in biofuel process economics [1, 2]. Hydrolysis of the polysaccharide constituents of biomass into simple sugars requires a diverse set of enzymes, namely endoglucanases (EC 184.108.40.206), exoglucanases (EC220.127.116.11), and β-glucosidases (EC 18.104.22.168) . In vitro endo- and exoglucanases collectively break down cellulose into cellobiose, a β-1,4-glucose disaccharide. Then, β-glucosidases hydrolyze cellobiose into two glucose monomers, which are subsequently fermented by the cell. To lessen the need of these in vitro hydrolytic enzymes, researchers have investigated the use of cellobiose itself as a fermentable sugar . Reducing the price of this in vitro enzyme pretreatment would lower the cost of biofuel production and make biofuels more economically competitive with petroleum-derived fuels [1, 2].
Three main strategies for heterologous cellobiose utilization have recently been developed. The first strategy involves cell-surface display of extracellular β-glucosidases [4–8]. In this strategy, the cellobiose is hydrolyzed extracellularly and then the glucose is transported into the cell and metabolized. A second strategy, the phosphorolytic pathway, relies on heterologous expression of a cellodextrin transporter and a cellobiose phosphorylase [9–11]. The cellobiose is transported into the cell via the cellodextrin transporter and the intracellular phosphorylase cleaves the disaccharide with an inorganic phosphate, producing a glucose molecule and an α-glucose-1-phosphate, which is quickly metabolized. The hydrolytic pathway is a third strategy, which involves heterologous expression of a cellodextrin transporter and an intracellular β-glucosidase [12–14]. After being transported into the cell by the cellodextrin transporter, the intracellular cellobiose is hydrolyzed by the β-glucosidase into two glucose molecules, which are then metabolized by the cell. Though each strategy has its own advantages, no method is as efficient as glucose assimilation, which is one of the most rapid sugar utilization rates. Therefore, further engineering is required to improve the cellobiose utilization rates to rival glucose consumption. The hydrolytic cellobiose utilization pathway has recently been investigated by our laboratory for optimization through combinatorial transcriptional engineering . Though this method and other transcriptional engineering techniques have been successful in optimizing metabolic pathways [15–18], optimal gene expression cannot overcome the inherent inefficiencies of the proteins within the pathway . Therefore, we sought to apply protein engineering strategies to improve the pathway performance.
Engineering proteins within a pathway has shown to be very successful in improving titer of the desired product . Often these methods tend to isolate the protein from the pathway and engineer it independently from other enzymes within the pathway. In this study, we have investigated a complementary strategy for pathway engineering, wherein we simultaneously modify multiple proteins within the context of the pathway and cellular metabolism. One of the advantages of this strategy is it can identify a mutant pathway with an increased and balanced flux. Engineering pathways with a balanced flux remains a large challenge in metabolic and pathway engineering [17, 18, 21–25]. By simultaneously screening mutants of all pathway proteins, mutants are not just chosen for a high protein activity, but instead are selected based on an optimized flux and balanced protein activity through the pathway. This method is also advantageous for its simple and efficient library creation which can generate mutations in each protein, allowing for a comprehensive exploration of the potential diversity of the pathway, possibly identifying mutations which could synergistically increase the pathway flux. This approach can be widely applicable to any system where a high-throughput screening/selection method for the desired phenotype is available. The implementation of this directed evolution strategy could further enhance researchers’ abilities to optimize pathways, offering a new avenue for metabolic engineering strategies to include pathway-scale approaches for protein engineering.
Library construction and screening
Characterization of the strains expressing the mutant cellobiose utilizing pathways
Parameters for the strains harboring the improved pathways identified through directed evolution
Specific Growth Rate
0.0694 ± 0.003
0.0919 ± 0.006
0.102 ± 0.002
1.78 ± 0.06
2.33 ± 0.02
2.65 ± 0.02
(g cellobiose/(L · h))
0.611 ± 0.02
0.815 ± 0.03
1.00 ± 0.03
(g ethanol/(L · h))
(g ethanol / g cellobiose)
0.4073 ± 0.02
0.4236 ± 0.0008
0.4363 ± 0.004
Specific activity measured from the crude lysate of the engineered β-glucosidase
0.00729 ± 0.0001
0.0021 ± 0.002
0.193 ± 0.01
0.0668 ± 0.001
0.130 ± 0.01
0.124 ± 0.004
0.219 ± 0.004
0.0738 ± 0.006
L173H + H23L
0.156 ± 0.007
0.170 ± 0.01
Specific activity measurements from the cellodextrin transporter as determined through radioactive-labeled cellobiose uptake rate assay
Cellobiose uptake rate
0.006 ± 0.001
1.77 ± 0.07
2.09 ± 0.05
2.63 ± 0.2
D433G + C82S
2.82 ± 0.2
Efficient production of biofuels and specialty chemicals hinges on optimization of the metabolic pathways associated with the desired product. In this endeavor, a variety of successful pathway optimization methods, mainly based on transcriptional engineering, have been developed [15–18, 41–44]. Though successful, these methods are unable to overcome inefficiencies based on innate enzyme characteristics such as Vmax and Km, therefore protein engineering strategies must also be applied . To improve these properties, often the enzymes are removed from the pathway context and engineered in vitro or in vivo for the highest activity. Engineering each enzyme independently can be arduous and there is no guarantee the engineered enzyme will reduce the bottleneck when it is reinstated into the pathway. We are introducing a new approach to combine pathway and protein engineering, involving the evolution of all proteins within the pathway simultaneously. Through this strategy, no a priori information about pathway bottlenecks or enzymes is required for the directed evolution. Mutant proteins of the pathway are not chosen for high activity, but instead interrogated as a whole for a balanced, increased activity and flux which is demonstrated by an improved pathway phenotype. This combinatorial mutagenesis library strategy allows for a comprehensive exploration of the pathway potential diversity, which could identify unexpected synergistic effects [42, 45]. Recent advances in large-scale library creation with high efficiency have made the multi-protein library a possibility [15, 36, 46–48]. Our study illustrates an in vivo library creation for simultaneous engineering of multiple proteins in a pathway through directed evolution in S. cerevisiae. To the best of our knowledge, there has been only one example of mutagenesis of multiple proteins through directed evolution, which was illustrated in Escherichia coli for arsenate resistance . This strategy has not been widely implemented and has not been applied to other systems. Our work optimizes a sugar utilization pathway for increased metabolite production with a simple one-step assembly method, allowing for quick and efficient rounds of evolution. This represents the first example of pathway optimization through multiple protein engineering in S. cerevisiae.
The final optimized pathway (R2) increased cellobiose consumption by 49% and ethanol productivity by 64%. Characterization of the strains expressing the mutations which conferred these improved phenotypes illustrated that the highest specific growth rates on cellobiose were achieved only when the β-glucosidase L173H mutation is present. In the cellodextrin transporter, the most significant mutation was the C82S mutation. It was shown that the optimized pathway involved an overall increase of protein activity: the β-glucosidase had an increased specificity towards cellobiose and the cellodextrin transporter exhibited a higher overall activity.
The β-glucosidase from glycoside hydrolase family 1 (GH1) contains the standard (α/β)-barrel structure of the enzyme family. The active site is located in a deep cleft formed by the connection loops at the C-terminus of the β-sheets of the TIM barrel with negatively charged residues surrounding the bottom of the active site . The conserved catalytic residues of the GH1 family involve the catalytic acid/base Glu 166 and the nucleophile Glu 377, which cleave the glycosidic bond via the common double displacement mechanism . This mechanism acts on β-1,4-linked glucose derivatives; the enzyme is promiscuous, being active towards a range of substrates with varying affinities . In this study, the wild-type GH1-1 β-glucosidase enzyme exhibited an innate substrate preference for p-NPG over cellobiose. The mutation L173H shifted the substrate specificity towards cellobiose. Homology modeling suggested the H173 residue could have a direct hydrogen bond with the hydroxyl group of the C1 atom of the cellobiose (Figure 4c), which was not predicted with the wild-type β-glucosidase. It is possible that this new hydrogen bond could increase the affinity of the enzyme for cellobiose by firmly keeping the substrate in place during hydrolysis, a hypothesis previously established by studies in substrate recognition for the β-glucosidase family [26, 29]. There has been significant interest in substrate specificity of the GH1 family, which is suggested to occur via an extensive network of polar interactions with well-ordered water molecules and enzyme-substrate bonding [26–29]. A specific study of interest involved a β-glucosidase from Trichoderma reesei [GenBank Accession number AB003110.1] which was recently investigated to identify critical residues within the substrate entrance region . Through rational design, substitutions at residue 172 were shown to improve substrate recognition, thermostability, and enzyme activity. This β-glucosidase has a 73% sequence identify to the GH1-1 enzyme studied here. Sequence alignment and structural modeling predicted this residue to be in the same position as the L173 residue from this work, corroborating the importance of this residue position. The H23L mutation is hypothesized to be an overall activity enhancer, with only slightly improved enzyme activity towards both substrates p-NPG and cellobiose. H23 is located on the periphery of the enzyme, distant from the active site (Figure 4a).
Engineering sugar transporters has recently been shown to greatly enhance pathway productivity [11, 30], a result also illustrated here. In our directed evolution effort, the combined D433G and C82S mutations of the cellodextrin transporter increased the specific activity by 60% compared to the wild-type. The specific activity is directly proportional to the Vmax of the transporters, which is the protein property that has been most often enhanced by sugar transporter engineering [11, 30]. It is probable that this property is the most likely to be improved due to the screening conditions, which were performed at saturating conditions. The cellodextrin transporter model predicted that the D433G mutation is located on a large outside loop. It is possible that this loop could have been causing a steric hindrance for the mass transfer of cellodextrin to the transporter, which could have been alleviated by the altered loop structure associated with this mutation. The C82S mutation, which conferred the highest increase in activity, is predicted to be located in the first transmembrane helix (Figure 5). The location in the helix and the increased activity suggest that this mutation could be directly associated with cellobiose interactions or the protein complex stability. Young et al. investigated directed evolution of xylose transporters and identified mutations which tended to cluster around the first transmembrane helix , a similar finding to this study.
Cellobiose utilization has recently become a significant consideration in economical biofuel production. Efficient utilization could decrease the expensive in vitro enzyme load in biomass pretreatment. We successfully optimized the hydrolytic cellobiose utilization pathway through a new directed evolution strategy: simultaneously engineering multiple proteins within the context of the pathway. The optimized pathway involved a β-glucosidase mutant with an increased specificity towards cellobiose and a cellodextrin transporter mutant with an increased overall activity. By applying directed evolution to the entire pathway, engineered proteins can be found which synergistically improve the phenotype either through a balanced flux of the pathway or through overall protein activity improvement. Simultaneous multi-protein engineering is expected to expand researchers’ abilities to optimize biosynthetic metabolic pathways.
Materials and methods
Strains, media, and culture conditions
The industrial Saccharomyces cerevisiae strain Still Spirits (Classic) Turbo Distiller’s Yeast was purchased from Homebrew Heaven (Everett, WA). Yeast strains were cultivated in YP media (1% yeast extract, 2% peptone) with 2% glucose (YPD) or 2-8% cellobiose (YPC). YPC with 8% cellobiose was used in fermentation analysis while 2% cellobiose was used for plate screening. S. cerevisiae strains were cultured at 30°C with orbital shaking at 250 rpm for aerobic growth or 100 rpm for oxygen-limited conditions in non-baffled flasks. As needed, 200 μg/mL G418 (KSE Scientific, Durham, NC) supplemented the YPD for pRS-KanMX plasmid selection. Escherichia coli DH5α (Cell Media Facility, University of Illinois at Urbana-Champaign, Urbana, IL) was used for recombinant DNA manipulations. E. coli strains were cultured in Luria broth (LB) (Fischer Scientific, Pittsburgh, PA) at 37°C and 250 rpm, supplemented with 100 μg/mL ampicillin. Yeast and bacterial strains were stored in 15% glycerol at −80°C. All chemicals were purchased from Sigma Aldrich or Fisher Scientific.
Plasmid and strain construction
Restriction enzymes were purchased from New England Biolabs (Ipswich, MA). All cloning work was performed via DNA Assembler, which relies on the in vivo homologous recombination mechanism in yeast . Yeast plasmids were isolated using Zymoprep Yeast Plasmid Miniprep II kit (Zymo Research, Irvine, CA), then transformed into E. coli for isolation of high purity DNA. E. coli plasmids were isolated using Qiagen Spin Plasmid Mini-Prep Kit (Qiagen, Valencia, CA). PCR fragments were purified by Qiagen QIAQuick Gel Purification Kit.
The β-glucosidase gene gh1-1 [GenBank Accession number XM_951090] from Neurospora crassa was expressed using a PYK1 promoter and an ADH1 terminator. The cellobiose transporter gene cdt-1 [GenBank Accession number XM_958708] from N. crassa was expressed using a TEF1 promoter and a PGK1 terminator as previously constructed . To transfer the pathway to the pRS-KanMX plasmid , primers kanMX-PYKp-F and PGKt-kanMX-R (Additional file 1: Table S1) were designed to amplify the full cellobiose utilization pathway. To facilitate the creation of a library of cellobiose utilization pathways, a pRS-KanMX helper plasmid was constructed. The helper plasmid contained the PYK1 promoter and the PGK1 terminator, separated by a unique restriction enzyme recognition site Bam HI for plasmid linearization. The helper plasmid was later linearized and used as a backbone for the library assembly. The PYK1 promoter was amplified using primers kanMX-PYKp-F and PYKp-BamHI-PGKt-R (Additional file 1: Table S1). The PGK1 terminator was amplified using primers PYKp-BamHI-PGKt-F and PGKt-kanMX-R (Additional file 1: Table S1).
The ADH1 terminator and TEF1 promoter were not subjected to random mutagenesis, thus this cassette was amplified independently using primers internal-BGL-ADHt-F and internal-TEFp-CDT-R (Additional file 1: Table S1) by standard PCR amplification. The gh1-1 and cdt-1 genes were subjected to error prone PCR with error rates of 2–3 basepair mutations per gene, resulting in one amino acid mutations per protein on average. Each fragment was amplified with 60–80 basepair homologies to upstream and downstream DNA sequences. The three fragments were co-transformed with the linearized helper plasmid into S. cerevisiae. The library size obtained was 104, consistent with libraries previously generated by this method [15, 36]. To confirm diversity of the library, ten colonies were randomly picked from the YPD + G418 plate and their plasmids were isolated and sequenced.
The library was screened on YPC plates for large colonies [15, 36] (Additional file 1: Figure S1). The growth rate of the largest colonies was quantified by small-scale fermentation. The colonies were seeded overnight in YPD + G418, washed twice with sterile water, and then inoculated into 10 mL YPC unbaffled flasks at 100 rpm and 30°C (Additional file 1: Figure S2). After confirmation of the top 5 clones with the fastest growth rate, the plasmids were isolated and retransformed. After retransformation, the strains were seeded overnight in YPD + G418, washed twice with sterile water, and then inoculated into 50 mL YPC in 250 mL unbaffled flasks at 100 rpm and 30°C. The cultures were sampled and analyzed for growth and metabolite production as described previously [15, 36]. After confirmation of the strains with improved phenotype, the plasmids were isolated and sequenced to identify the mutations within the genes of the pathway. The genes were used as a template for a second round of directed evolution.
Seed cultures were inoculated from plates of freshly streaked frozen culture stocks. YPD + G418 seed cultures were used for the individual and combined mutant pathway fermentation tests. YPC seed cultures were used for final mutant pathway analysis. The seed cultures were grown at 30°C and 250 rpm overnight, then washed with sterile water twice before inoculating 50 mL YPC in 250 mL un-baffled flasks to an initial OD of 0.2. The cultures were incubated at 30°C and 100 rpm orbital shaking for oxygen-limited growth . The cultures were sampled and analyzed for growth and metabolite production as previously described [15, 36]. For the metabolite analysis, end-point values at the time required to consume more than 95% of total sugar was used to determine rate, yield, and productivity.
Construction of single mutants
The single mutant genes were constructed through site-directed mutagenesis and the mega primer PCR method with primers listed in Additional file 1: Table S1. After PCR amplification, the genes were transformed into the Classic strain along with the ADHt1/TEFp1 cassette into a linearized pRS-KanMX helper plasmid. The final plasmid was purified and sequenced for confirmation.
β-Glucosidase enzyme activity assays
Classic strains harboring an empty vector, the wild-type, and mutant pathways were grown to mid-exponential phase and washed three times with potassium phosphate buffer (pH 7). A final cell mass equivalent to an OD600 of 20 was harvested. Cell-free extracts were prepared through YPER Extraction Reagent (Thermo Scientific, Rockport, IL), 125 μL of YPER was used to lyse the cells for 20 minutes at 25°C with vigorous shaking at 700 rpm in a thermomixer (Eppendorf, Germany). After lysing, the cell debris was pelleted for 10 minutes with 15,000 rpm at 4°C. The total protein concentration was determined via the BCA protein assay kit (Pierce, Rockford, IL), the standard manufacturer protocol was followed.
The lysate was tested for p-NPG activity with 1 mM p-NPG in 100 mM potassium phosphate buffer at pH 7. The colorimetric change was monitored at 405 nm in a 96-well Biotech Synergy 2 plate-reader (Winooski, VT) for 30 minutes at 30°C. The amount of p-nitrophenol (p-NP) released was quantified by a p-NP standard in100 mM potassium phosphate buffer pH 7. One unit of activity is equivalent to 1 μmol of p-NP released per min. For cellobiose-based enzymatic assay, the linear range of the β-glucosidase was determined. The reactions were carried out at 60 mM cellobiose in 100 mM potassium phosphate buffer pH 7. After addition of lysate, the reaction was allowed to react for 15, 30, 45, and 60 minutes at 30°C. The reaction was stopped by boiling at 100°C for 10 minutes. The samples were then centrifuged at 15,000 rpm for 5 minutes before being stored on ice. The amount of glucose which had been produced in the allotted time frame was then measured using the D-glucose kit (R-Biopharm, Germany). Standard manufacturer’s instructions were followed. To determine the β-glucosidase protein abundance in each strain, the wild-type and mutant proteins were fused to GFP. There was no statistically significant difference in β-glucosidase abundance of each strain (Additional file 1: Figure S4). The total protein concentration, as determined through BCA assay, was comparable in each lysate. Therefore, the ratio of β-glucosidase to total protein was assumed constant, and hence the activity was normalized to total protein concentration. One unit (U) is defined as one micromole of glucose produced per minute.
CDT activity assay
The cellodextrin transporter was assayed using the oil-stop protocol previously reported [12, 38]. Cultures of the wild-type and mutant strains were grown in YPC to an OD of 15–20, washed three times with ice-cold assay buffer (30 mM MES-NaOH + 50 mM ethanol), and then normalized to an OD of 20. 50 μL of cells were added to 50 μL of [3H]-cellobiose (Movarek Biochemicals, Brea, CA) at 30°C and layered over 100 μL of silicone oil (Sigma 85419), incubated for 10, 20, 40, and 80 seconds. The cells were then centrifuged through the oil at 15,000 rpm for 30 seconds. After being placed in ethanol/dry-ice bath, the cell pellets were solubilized in NaOH overnight. The amount of [3H]-cellobiose present in the cells was then quantified by a Beckman Coulter LS 6500 liquid scintillation counter (Brea, CA). The amount of labeled cellobiose that was transported into the cell was plotted against the time of reaction. One unit (U) is defined as one micromole of cellobiose taken up by the cell per min. The rate of the cellobiose uptake was normalized by the transporter abundance determined through GFP fluorescence measurements (Additional file 1: Figure S4) and the gram cell dry weight (gcdw).
A homology model of the β-glucosidase enzyme was constructed to identify the structure-function relationships of the mutations discovered. The gene encoding the β-glucosidase from Trichoderma reesei [PDB accession code 3AHY] afforded the highest homology to the gh1-1 gene from N. crassa with 73% sequence identity and was used as a template for homology modeling. The structure model of the β-glucosidase from N. crassa was constructed using the modeling program Molecular Operating Environment (Chemical Computing Group, Montreal, Canada). After constructing the homology model, the substrate in the co-crystal structure of Neotermes koshunensis β-glucosidase [PDB accession code 3VIK] was docked into the model . The model was energy minimized before ligand interactions were investigated. To identify the effects of the mutations, the mutations were introduced to the model and the energy was minimized again before investigating the ligand interactions.
Gram cell dry weight.
This work was supported by the Energy Biosciences Institute. DTE acknowledges a National Science Foundation Graduate Research Fellowship. Ryan Cobb is acknowledged for assistance with the molecular modeling. Dr. George Schmitz is appreciated for helpful discussions. The Roy J. Carver Biotechnology Center (CBC) is acknowledged for use of the flow cytometer. Jiazhang Lian is also acknowledged for his help in the cellodextrin transporter assays.
- Wyman CE: What is (and is not) vital to advancing cellulosic ethanol. Trends Biotechnol. 2007, 25: 153-157. 10.1016/j.tibtech.2007.02.009.View Article
- Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS: Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev. 2002, 66: 506-577. 10.1128/MMBR.66.3.506-577.2002.View Article
- Elkins JG, Raman B, Keller M: Engineered microbial systems for enhanced conversion of lignocellulosic biomass. Curr Opin Biotechnol. 2010, 21: 657-662. 10.1016/j.copbio.2010.05.008.View Article
- van Rooyen R, Hahn-Hagerdal B, La Grange DC, van Zyl WH: Construction of cellobiose-growing and fermenting Saccharomyces cerevisiae strains. J Biotechnol. 2005, 120: 284-295. 10.1016/j.jbiotec.2005.06.013.View Article
- Nakamura N, Yamada R, Katahira S, Tanaka T, Fukuda H, Kondo A: Effective xylose/cellobiose co-fermentation and ethanol production by xylose-assimilating S. cerevisiae via expression of β-glucosidase on its cell surface. Enzyme and Microbial Technology. 2008, 43: 233-236. 10.1016/j.enzmictec.2008.04.003.View Article
- Guo ZP, Zhang L, Ding ZY, Gu ZH, Shi GY: Development of an industrial ethanol-producing yeast strain for efficient utilization of cellobiose. Enzyme Microb Technol. 2011, 49: 105-112. 10.1016/j.enzmictec.2011.02.008.View Article
- Saitoh S, Hasunuma T, Tanaka T, Kondo A: Co-fermentation of cellobiose and xylose using β-glucosidase displaying diploid industrial yeast strain OC-2. Appl Microbiol Biotechnol. 2010, 87: 1975-1982. 10.1007/s00253-010-2714-0.View Article
- Machida M, Ohtsuki I, Fukui S, Yamashita I: Nucleotide sequences of Saccharomycopsis fibuligera genes for extracellular β-glucosidases as expressed in Saccharomyces cerevisiae. Appl Environ Microbiol. 1988, 54: 3147-3155.
- Sadie CJ, Rose SH, den Haan R, van Zyl WH: Co-expression of a cellobiose phosphorylase and lactose permease enables intracellular cellobiose utilisation by Saccharomyces cerevisiae. Appl Microbiol Biotechnol. 2011, 90: 1373-1380. 10.1007/s00253-011-3164-z.View Article
- Aeling KA, Salmon KA, Laplaza JM, Li L, Headman JR, Hutagalung A, Picataggio S: Co-fermentation of xylose and cellobiose by an engineered Saccharomyces cerevisiae. J Ind Microbiol Biotechnol. 2012, 39: 1597-1604. 10.1007/s10295-012-1169-y.View Article
- Ha SJ, Galazka JM, Joong Oh E, Kordic V, Kim H, Jin YS, Cate JH: Energetic benefits and rapid cellobiose fermentation by Saccharomyces cerevisiae expressing cellobiose phosphorylase and mutant cellodextrin transporters. Metab Eng. 2012, 15: 134-143.View Article
- Galazka JM, Tian C, Beeson WT, Martinez B, Glass NL, Cate JHD: Cellodextrin transport in yeast for improved biofuel production. Science. 2010, 330: 84-86. 10.1126/science.1192838.View Article
- Ha SJ, Galazka JM, Kim SR, Choi JH, Yang X, Seo JH, Glass NL, Cate JH, Jin YS: Engineered Saccharomyces cerevisiae capable of simultaneous cellobiose and xylose fermentation. Proc Natl Acad Sci U S A. 2011, 108: 504-509. 10.1073/pnas.1010456108.View Article
- Li S, Du J, Sun J, Galazka JM, Glass NL, Cate JHD, Yang X, Zhao H: Overcoming glucose repression in mixed sugar fermentation by co-expressing a cellobiose transporter and a β-glucosidase in Saccharomyces cerevisiae. Mol Biosyst. 2010, 6: 2129-2132. 10.1039/c0mb00063a.View Article
- Du J, Yuan Y, Si T, Lian J, Zhao H: Customized optimization of metabolic pathways by combinatorial transcriptional engineering. Nucleic Acids Res. 2012, 40: e142- 10.1093/nar/gks549.View Article
- Alper H, Fischer C, Nevoigt E, Stephanopoulos G: Tuning genetic control through promoter engineering. Proc Natl Acad Sci USA. 2005, 102: 12678-12683. 10.1073/pnas.0504604102.View Article
- Salis HM, Mirsky EA, Voigt CA: Automated design of synthetic ribosome binding sites to control protein expression. Nat Biotechnol. 2009, 27: 946-950. 10.1038/nbt.1568.View Article
- Pfleger BF, Pitera DJ, Smolke CD, Keasling JD: Combinatorial engineering of intergenic regions in operons tunes expression of multiple genes. Nat Biotechnol. 2006, 24: 1027-1032. 10.1038/nbt1226.View Article
- Yoshikuni Y, Dietrich JA, Nowroozi FF, Babbitt PC, Keasling JD: Redesigning enzymes based on adaptive evolution for optimal function in synthetic metabolic pathways. Chem Biol. 2008, 15: 607-618. 10.1016/j.chembiol.2008.05.006.View Article
- Eriksen DT, Lian J, Zhao H: Protein design for pathway engineering. J Struct Biol. 2013, 10.1016/j.jsb.2013.03.011.
- Dueber JE, Wu GC, Malmirchegini GR, Moon TS, Petzold CJ, Ullal AV, Prather KL, Keasling JD: Synthetic protein scaffolds provide modular control over metabolic flux. Nat Biotechnol. 2009, 27: 753-759. 10.1038/nbt.1557.View Article
- Anthony JR, Anthony LC, Nowroozi F, Kwon G, Newman JD, Keasling JD: Optimization of the mevalonate-based isoprenoid biosynthetic pathway in Escherichia coli for production of the anti-malarial drug precursor amorpha-4, 11-diene. Metab Eng. 2009, 11: 13-19. 10.1016/j.ymben.2008.07.007.View Article
- Lutke-Eversloh T, Stephanopoulos G: Combinatorial pathway analysis for improved L-tyrosine production in Escherichia coli: identification of enzymatic bottlenecks by systematic gene overexpression. Metab Eng. 2008, 10: 69-77. 10.1016/j.ymben.2007.12.001.View Article
- Bond-Watts BB, Bellerose RJ, Chang MC: Enzyme mechanism as a kinetic control element for designing synthetic biofuel pathways. Nat Chem Biol. 2011, 7: 222-227. 10.1038/nchembio.537.View Article
- Shen CR, Lan EI, Dekishima Y, Baez A, Cho KM, Liao JC: Driving forces enable high-titer anaerobic 1-butanol synthesis in Escherichia coli. Appl Environ Microbiol. 2011, 77: 2905-2915. 10.1128/AEM.03034-10.View Article
- Isorna P, Polaina J, Latorre-Garcia L, Canada FJ, Gonzalez B, Sanz-Aparicio J: Crystal structures of Paenibacillus polymyxa β-glucosidase B complexes reveal the molecular basis of substrate specificity and give new insights into the catalytic machinery of family I glycosidases. J Mol Biol. 2007, 371: 1204-1218. 10.1016/j.jmb.2007.05.082.View Article
- Lee HL, Chang CK, Jeng WY, Wang AH, Liang PH: Mutations in the substrate entrance region of β-glucosidase from Trichoderma reesei improve enzyme activity and thermostability. Protein Eng Des Sel. 2012, 25: 733-740. 10.1093/protein/gzs073.View Article
- Jeng WY, Wang NC, Lin MH, Lin CT, Liaw YC, Chang WJ, Liu C, Liang PH, Wang AHJ: Structural and functional analysis of three β-glucosidases from bacterium Clostridium cellulovorans, fungus Trichoderma reesei and termite Neotermes koshunensis. J Struct Biol. 2011, 173: 46-56. 10.1016/j.jsb.2010.07.008.View Article
- Nam KH, Sung MW, Hwang KY: Structural insights into the substrate recognition properties of β-glucosidase. Biochem Biophys Res Commun. 2010, 391: 1131-1135. 10.1016/j.bbrc.2009.12.038.View Article
- Young EM, Comer AD, Huang H, Alper HS: A molecular transporter engineering approach to improving xylose catabolism in Saccharomyces cerevisiae. Metab Eng. 2012, 14: 401-411. 10.1016/j.ymben.2012.03.004.View Article
- Bokma E, Koronakis E, Lobedanz S, Hughes C, Koronakis V: Directed evolution of a bacterial efflux pump: adaptation of the E. coli TolC exit duct to the Pseudomonas MexAB translocase. FEBS Lett. 2006, 580: 5339-5343. 10.1016/j.febslet.2006.09.005.View Article
- Dunlop MJ, Dossani ZY, Szmidt HL, Chu HC, Lee TS, Keasling JD, Hadi MZ, Mukhopadhyay A: Engineering microbial biofuel tolerance and export using efflux pumps. Mol Syst Biol. 2011, 7: 487-View Article
- Ren C, Chen T, Zhang J, Liang L, Lin Z: An evolved xylose transporter from Zymomonas mobilis enhances sugar transport in Escherichia coli. Microb Cell Fact. 2009, 8: 66-75. 10.1186/1475-2859-8-66.View Article
- Shao ZY, Zhao H, Zhao HM: DNA assembler, an in vivo genetic method for rapid construction of biochemical pathways. Nucleic Acids Res. 2009, 37: e61- 10.1093/nar/gkp168.View Article
- Firth AE, Patrick WM: Statistics of protein library construction. Bioinformatics. 2005, 21: 3314-3315. 10.1093/bioinformatics/bti516.View Article
- Kim B, Du J, Eriksen DT, Zhao H: Combinatorial design of a highly efficient xylose utilizing pathway for cellulosic biofuels production in Saccharomyces cerevisiae. Appl Environ Microbiol. 2013, 79: 931-941. 10.1128/AEM.02736-12.View Article
- Kim SR, Kwee NR, Kim H, Jin YS: Feasibility of xylose fermentation by engineered Saccharomyces cerevisiae overexpressing endogenous aldose reductase (GRE3), xylitol dehydrogenase (XYL2), and xylulokinase (XYL3) from Scheffersomyces stipitis. FEMS Yeast Res. 2013, 13: 312-321. 10.1111/1567-1364.12036.View Article
- Arendt CS, Ri K, Yates PA, Ullman B: Genetic selection for a highly functional cysteine-less membrane protein using site saturation mutagenesis. Anal Biochem. 2007, 365: 185-193. 10.1016/j.ab.2007.03.030.View Article
- Tusnady GE, Simon I: Principles governing amino acid composition of integral membrane proteins: application to topology prediction. J Mol Biol. 1998, 283: 489-506. 10.1006/jmbi.1998.2107.View Article
- Tusnady GE, Simon I: The HMMTOP transmembrane topology prediction server. Bioinformatics. 2001, 17: 849-850. 10.1093/bioinformatics/17.9.849.View Article
- Alper H, Stephanopoulos G: Global transcription machinery engineering: a new approach for improving cellular phenotype. Metab Eng. 2007, 9: 258-267. 10.1016/j.ymben.2006.12.002.View Article
- Wang HH, Isaacs FJ, Carr PA, Sun ZZ, Xu G, Forest CR, Church GM: Programming cells by multiplex genome engineering and accelerated evolution. Nature. 2009, 460: 894-898. 10.1038/nature08187.View Article
- Warnecke TE, Lynch MD, Karimpour-Fard A, Lipscomb ML, Handke P, Mills T, Ramey CJ, Hoang T, Gill RT: Rapid dissection of a complex phenotype through genomic-scale mapping of fitness altering genes. Metab Eng. 2010, 12: 241-250. 10.1016/j.ymben.2009.12.002.View Article
- Warner JR, Reeder PJ, Karimpour-Fard A, Woodruff LBA, Gill RT: Rapid profiling of a microbial genome using mixtures of barcoded oligonucleotides. Nat Biotechnol. 2010, 28: 856-862. 10.1038/nbt.1653.View Article
- Alper H, Miyaoku K, Stephanopoulos G: Construction of lycopene-overproducing E. coli strains by combining systematic and combinatorial gene knockout targets. Nat Biotechnol. 2005, 23: 612-616. 10.1038/nbt1083.View Article
- Wingler LM, Cornish VW: Reiterative recombination for the in vivo assembly of libraries of multigene pathways. Proc Natl Acad Sci USA. 2011, 108: 15135-15140. 10.1073/pnas.1100507108.View Article
- Ramon A, Smith HO: Single-step linker-based combinatorial assembly of promoter and gene cassettes for pathway engineering. Biotechnol Lett. 2011, 33: 549-555. 10.1007/s10529-010-0455-x.View Article
- Quan J, Tian J: Circular polymerase extension cloning of complex gene libraries and pathways. PLoS One. 2009, 4: e6441- 10.1371/journal.pone.0006441.View Article
- Crameri A, Dawes G, Rodriguez E, Silver S, Stemmer WP: Molecular evolution of an arsenate detoxification pathway by DNA shuffling. Nat Biotechnol. 1997, 15: 436-438. 10.1038/nbt0597-436.View Article
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 cited.