ORIGINAL RESEARCH article
Synergistic function of four novel thermostable glycoside hydrolases from a long-term enriched thermophilic methanogenic digester
- 1School of Biotechnology, Jiangnan University, Wuxi, China
- 2Department of Energy and Resources Engineering, College of Engineering, Peking University, Beijing, China
- 3Institute of Engineering (Baotou), College of Engineering, Peking University, Baotou, China
- 4Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
In biofuel production from lignocellulose, low thermostability and product inhibition strongly restrict the enzyme activities and production process. Application of multiple thermostable glycoside hydrolases, forming an enzyme “cocktail”, can result in a synergistic action and therefore improve production efficiency and reduce operational costs. Therefore, increasing enzyme thermostabilities and compatibility are important for the biofuel industry. In this study, we reported the screening, cloning and biochemical characterization of four novel thermostable lignocellulose hydrolases from a metagenomic library of a long-term dry thermophilic methanogenic digester community, which were highly compatible with optimal conditions and specific activities. The optimal temperatures of the four enzymes, β-xylosidase, xylanase, β-glucosidase, and cellulase ranged from 60 to 75°C, and over 80% residual activities were observed after 2 h incubation at 50°C. Mixtures of these hydrolases retained high residual synergistic activities after incubation with cellulose, xylan, and steam-exploded corncob at 50°C for 72 h. In addition, about 55% dry weight of steam-exploded corncob was hydrolyzed to glucose and xylose by the synergistic action of the four enzymes at 50°C for 48 h. This work suggested that since different enzymes from a same ecosystem could be more compatible, screening enzymes from a long-term enriching community could be a favorable strategy.
Enzymes catalyzing the hydrolysis of cellulose and hemicellulose are a source of worldwide interest for alternative bio-energy development and other environmental applications, such as treating wastewater from the dying, textile and pulping industries. Many enzymes are used in these processes, including exoglucanase (EC 22.214.171.124), endoglucanase (EC 126.96.36.199), and β-glucosidase (EC 188.8.131.52) for cellulose hydrolysis; xylanase (EC 184.108.40.206) and β-xylosidase (EC 220.127.116.11) to hydrolyze xylan; and α-L-arabinofuranosidase (EC 18.104.22.168) and α-glucuronidase (EC 22.214.171.124) to hydrolyze side chains in hemicellulose (Beg et al., 2001). To utilize these enzymes in biofuel production, two processes are generally used: separate hydrolysis and fermentation (SHF), and simultaneous saccharification and fermentation (SSF). Temperatures of hydrolysis in SHF and SSF are around 50 and 32°C, respectively (Tomas-Pejo et al., 2008). SHF has been reported to be more efficient and cost-effective than SSF (Magnus et al., 2012) because it uses thermostable enzymes, which improves the solubility of substrates and enhances the mass transfer rate, while reducing cooling requirements and the probability of bacterial contamination (Krahe et al., 1996). However, product inhibition can be a major source of inefficiency in SHF, especially the inhibition on cellulase and xylanase (Zhang and Lynd, 2004). Therefore, combinations of relevant enzymes, forming a so-called enzyme cocktail, have been intensively investigated in recent decades (Zhou et al., 2009; Banerjee et al., 2010). For example, the combination of exoglucanase, endoglucanase, and β-glucosidase is used to increase cellulose degradation efficiency (Gusakov et al., 2007), and the combination of xylanase and β-xylosidase has been shown to enhance hemicellulose hydrolysis (Fan et al., 2009). In general, these cocktails are formed from widely used commercial enzymes, such as the cellobiohydrolases I and II and endoglucanases I and II from the fungus Trichoderma reesei (Henrissat et al., 1985). However, since the host microorganisms from which these enzymes were originally isolated are from different habitats and have different optimal growth conditions, the isolated enzymes generally also have different optimal conditions, and may not be compatible in the same cocktail.
Metagenomic approaches for screening novel enzymes have the potential to increase the compatibility of the candidate enzymes in a cocktail, because these metagenome-oriented enzymes can be from the same microbial community, living in the same habitat, and even in the same niche. For example, hundreds of lignocellulose hydrolases were screened from metagenomic libraries of the rumen (Hess et al., 2011), the gut system of termites (Warnecke et al., 2007), soil (Mcandrew et al., 2013), and biogas digesters (Yan et al., 2013). Nevertheless, more hydrolases are believed to be awaiting discovery in the “microbial dark matter” (Rinke et al., 2013). So far, most of this research has focused on the desirable properties of the discovered enzymes, such as thermostability, ionic strength, alkaline tolerance, activity at cold temperatures, and wide substrate specificity. Comparatively few efforts have been made to investigate their compatibility and synergistic actions (Del Pozo et al., 2012; Geng et al., 2012; Gruninger et al., 2014).
For obtaining thermostable enzymes, we constructed a metagenomic library from a dry thermophilic methanogenic digester community enriched mainly with paper for over 10 years (Tang et al., 2011). We isolated and identified four novel thermostable hydrolases and found them to have high compatibility as a synergistic cocktail for hydrolyzing cellulose, xylan, and steam-exploded corncob. This work further proved that metagenomic screening is an effective way for finding compatible enzymes.
Materials and Methods
Construction, Functional Screening and In Silico Analysis of Fosmid Metagenomic Library
The genomic DNA of the thermophilic methanogenic digester community, as described previously (Tang et al., 2011), was extracted using FastDNA® SPIN Kit for soil (MP, Santa Ana, America). The DNA purification, ligation, and transformation were performed according to the manufacturer's instructions of Copy Control™ Fosmid Library Production Kit (Epicentre, Madison, America). A fosmid library was constructed, and the analysis of insert DNA fragments was carried out as described previously (Hu et al., 2010). Activities of cellulase and xylanase were screened from the library with the Congo red method (Wood et al., 1988). To screen β-glucosidase activity, the ammonium ferric citrate-esculin hydrate method (Eberhart et al., 1964) was used, while pNPX (p-nitrophenyl-β-D-xyloside) was used to identify β-xylosidase (Poutanen and Puls, 1988). Fosmid clones with positive activities were sequenced using the GS FLX system (Roche, Basel, Switzerland). The sequences were then assembled with the putative open reading frames (ORFs) predicted using SoftBerry fgenesb (Warnecke et al., 2007). The deduced proteins were then annotated against the NCBI non-redundant database. The Pfam (Finn et al., 2014) and CAZY (Lombard et al., 2014) databases were used to identify conserved domains and glycoside hydrolase families, respectively. Signal peptides were identified using SignalP 4.1 (Petersen et al., 2011). Homology models of the proteins were predicted from the SWISS-MODEL server (Arnold et al., 2006) and visualized using PyMOL (Schrodinger, 2010). Active sites of the proteins were identified by multiple sequence alignment (Mcwilliam et al., 2013). The composition of amino acid residues was analyzed online (http://www.bio-soft.net/sms/). In total, four genes from three fosmid clones were selected for further study: Xyl522, coding for β-xylosidase; Xyn526, coding for xylanase; Bgl8520, coding for β-glucosidase; and Cel1753, coding for cellulase. Phylogenetic analysis was performed by MEGA5.1 with the best hits in BLASTP and characterized protein sequences from the same protein family, using the neighbor-joining algorithm.
Cloning, Expression, Purification, and Characterization of Target Enzymes
The four target genes were PCR amplified from the corresponding fosmids with the primers listed in Supplementary Table 1 using PrimeSTAR HS DNA Polymerase (Takara, Dalian, China) according to the manufacturer's instructions. After purification and double digestion, the PCR amplicons were cloned into pET-28a (+) (Qiagen, Hilden, Germany). Recombinants of Xyl522, Xyn526, Bgl8520, and Cel1753 were further transformed into E. coli BL21 (DE3) (CoWin Biotech, Beijing, China). Recombinant cells were cultured in Lysogeny Broth medium with 50 μg/mL kanamycin and induced by 1 mM IPTG at 37°C for 4 h. Cells were harvested and subjected to freeze-thaw lysis. Next, the cell lysate was centrifuged at 4°C, 12,000 × g for 30 min, and the supernatant was collected for protein purification though Ni-NTA resin (GE Healthcare, Uppsala, Sweden) according to the manufacturer's instructions. The protein concentration was determined using the BCA Protein Assay Kit (Tiangen, Beijing, China).
Enzyme activities were detected as follows. pNPX (25 mM final concentration) and pNPG (p-nitrophenyl-β-D-glucoside, 25 mM final concentration) were mixed with 1 μg Xyl522 and 1 μg Bgl8520, respectively, each to a final volume of 40 μL (pH 5.0). After incubation at 50°C for 5 min, 50 μL 1 M Na2CO3 was added to terminate the reaction. The concentrations of p-nitrophenol product were determined by detecting the absorbance at 410 nm using spectrophotometer UV1700 (Shimadzu, Kyoto, Japan). The specific activities (U) of β-xylosidase and β-glucosidase were shown in micromole p-nitrophenol produced from pNPX and pNPG per minute per milligram of protein. Similarly, beech wood xylan (1%, w/v final concentration) and CMC (carboxymethylcellulose, 1%, w/v final concentration) were mixed with 5 μg Xyn526 and 5 μg Cel1753, respectively, to a final volume of 200 μL (pH 5.0). After incubation at 50°C for 15 min, the reaction mixtures were added to 500 μL dinitrosalicylic acid and boiled for 5 min to terminate the reactions (Miller, 1959). Concentrations of released xylose or glucose were determined by analyzing the absorbance of the reaction mixtures at 540 nm. Activities of cellulase and xylanase were defined as micromole xylose and glucose produced per minute per milligram of protein. The determination of filter paper activity (FPA) was performed as previously described (Zhang et al., 2010); FPA was defined as micromole glucose produced per minute per milligram of protein.
The pH and temperature optima of the four enzymes were determined in the range of pH 3.0–12.0 (Xyl522, pH 3.5–7.0; Xyn526, pH 3.5–11.0; Bgl8520, pH 4.0–9.0; and Cel1753, pH 3.5–9.0) and 20–90°C (Xyl522 and Bgl8520, 20–85°C; Xyn526, 20–90°C; and Cel1753, 20–80°C). To maintain suitable pH values, citric acid-Na2HPO4 buffer (pH 3.0–8.0), Tris-HCl (pH 8.0–10.0) and Na2HPO4-NaOH (pH 10.0–12.0) were used. For determination of optimal temperatures, pH was maintained at 5.0, for consistency with conditions in the methanogenic digester community from which the enzymes originate. To test for pH stability, the four enzymes were incubated in the buffer solutions mentioned above with corresponding pH values at 4°C for 24 h, after which their residual enzyme activities were measured as described above. For example, Xyn526 was incubated at 4°C in buffer solutions ranging from pH 3.5 to 12.0 for 24 h, after which mixtures containing 5 μg xylanase protein were used for activity determination as described above. To test the thermostability of the four enzymes, Xyl522 and Bgl8520 were incubated at 50, 55, 60, and 75°C; Xyn526 incubated at 50, 65, 70, and 75°C; and Cel1753 incubated at 50, 60, 65, and 70°C for 2 h in buffer solutions of pH 4.5, 6.0, 5.0, and 5.0, respectively. The residual activities of all enzymes were tested at 50°C and pH 5.0. Briefly, mixtures containing 1 μg Xyl522 and Bgl8520 each, or 5 μg Xyn526 and Cel1753 each were sampled every 20 min until 120 min, and enzyme activity was determined for samples taken at each time point. In addition, CMC, microcrystalline cellulose (MIC), beech wood xylan, pNPG, pNPC (pNP-β-D-cellobioside), pNPX, and pNPA (pNP-α-L-arabinofuranoside) were used to investigate the substrate specificity of the enzymes. Km and Vmax were calculated by Lineweaver–Burk plot based on enzyme activities at their optimal conditions with different substrate concentrations (1–50 mM for pNPG, pNPX, pNPA, and pNPC; 1–50 g/L for xylan beech wood, xylan oat spelt, CMC, and MIC).
Synergistic Action of Enzymes on CMC, Xylan, and Steam-exploded Corncob
The activities of Cel1753 and Bgl8520 cocktail against CMC (1%, w/v) and Xyn526 and Xyl522 cocktail against xylan (1%, w/v) were assayed by detecting the products using thin-layer chromatography (TLC). After the mixtures were incubated at 50°C, pH 5.0 for 4 h, respectively, 5 μL of each culture was spotted on TLC Silica gel 60 (Merck, Darmstadt, Germany) (Murashima et al., 2003). TLC plates were developed with ethyl acetate/methanol/acetic acid/water (12:3:3:1, v/v) as the mobile phase, and dyed with chromogenic agent containing 4 g diphenylamine, 4 mL phenylamine, 20 mL 85% phosphoric acid, and 200 mL acetone for 15 s followed by heating at 80°C for 10 min.
Steam-exploded corncob was prepared according to a previously described protocol (Teng et al., 2010) with a dry weight of 32%, which was composed of cellulose (~35%), hemicellulose (~38%) and lignin (~18%). The corncob biomass (2%, w/v) were incubated with Cel1753 (40 FPA/g, dry weight, as follows), Xyn526 (200 U/g), a cocktail of Cel1753 (40 FPA/g), and Bgl8520 (50 U/g), a cocktail of Xyn526 (200 U/g) and Xyl522 (50 U/g), and a cocktail of all four enzymes at 50°C, pH 5.0 for 72 h. Enzyme dosages used above were the minimum ones to achieve the saturated hydrolysis of corncobs by itself. The concentrations of glucose, cellobiose, xylose, and xylobiose were measured using a Prominence Liquid Chromatograph-20A (Shimadzu) according to the manufacturer's instructions. Briefly, a fluorometric detector and a shim-pack ISA-07 column were used. The system was operated at 65°C and a flow rate of 0.6 mL/min using H3BO3-KOH as the mobile phase with 1% w/v L-arginine for product derivatization. Calculation of saccharification efficiency was performed as described previously (Del Pozo et al., 2012).
Nucleotide Sequence Accession Numbers
The GenBank accession numbers of Xyl522, Xyn526, Bgl8520, and Cel1753 genes were KM982177, KM982178, KM982179, and KM982180, respectively.
Screening of the Metagenomic Library and Sequence Analysis
The fosmid metagenomic library contained about 92,000 clones. The average length of the insert DNA fragments was 36.0 kb. Altogether, we identified 49, 37, and 92 clones harboring cellulase, xylanase and β-glucosidase activities from 9,700 clones, respectively. According to the average insert size of fosmid clones, the hit rates (gene/Mbp) of cellulase, xylanase and β-glucosidase were 1/8.1, 1/10.8, and 1/4.4, respectively, which were higher than those metagenomes from a 1.5 year enriched biogas digester (Yan et al., 2013), from buffalo rumen (Duan et al., 2009) and from soil (Kim et al., 2008), suggesting that long-term enrichment could increase the possibility of identifying the target genes. Three positive fosmids, designated F52, F85, and F175, were found positive for xylanase, β-glucosidase, and cellulase activities, respectively. Complete sequence analysis revealed that fosmids F52, F85, and F175 had inserts of 38,176, 34,045, and 37,655 bp, and contained 26, 29, and 33 ORFs, respectively (Supplementary Figures 1–3). From the inserted DNA fragments in fosmids F52, F85, and F175, putative genes coding for β-xylosidase, xylanase, β-glucosidase, and cellulase and their glycoside hydrolase (GH) family were annotated (Supplementary Table 2) and designated as Xyl522, Xyn526, Bgl8520, and Cel1753, respectively. The four putative enzymes shared the highest amino acid sequence similarities with β-xylosidase from Clostridium stercorarium (63%, Genbank accession WP_015357972.1), 1,4-β-xylanase from Caldicoprobacter oshimai (56%, Genbank accession WP_025746830.1), β-glucosidase from bacterium UASB270 (60%, Genbank accession WP_021653982.1), and cellulase from an uncultured bacterium (59%, Genbank accession AEV59736.1). Xyl522 and Bgl8520 did not have signal peptides, while Xyn526 and Cel1753 had, suggesting that Xyl522 and Bgl8520 were intracellular, while Xyn526 and Cel1753 were functional only when secreted. It is notable that the alanine, proline, arginine, and glutamic acid residue contents in our proteins were generally around 1–4% higher than those of the mesothermal enzymes, whereas the glycine, serine, lysine and aspartic acid content was around 1–4% lower (Supplementary Table 4) (Morana et al., 2007; Alvarez et al., 2013; Han et al., 2013; Mcandrew et al., 2013). These differences in relative content of specific amino acid residues could be related to the temperature adaptabilities of these enzymes.
The predicted models of enzyme structure (Supplementary Figure 4) suggested that the four enzymes all had a classical (α/β)8 TIM barrel fold structure. Xyl522, a member of the GH 3 family, also contained an immunoglobulin-like region, which has been proven to contribute to both substrate binding and dimerization in this family (Mcandrew et al., 2013). Additionally, the distance between the two active sites in Xyl522 was calculated as 5.4 Å, In Bgxa1, a multifunctional β-glucosidase/β-xylosidase/α-arabinosidase, the same distance is 6.1 Å (Gruninger et al., 2014). Although Bgxa1 and Xyl522 are from the same GH 3 family, this distinction might suggest different substrate specificities.
Cloning, Expression, Purification, and Characterization of Xyl522, Xyn526, Bgl8520, and Cel1753
Xyl522, Xyn526, Bgl8520, and Cel1753, without their signal peptides, were heterologously expressed and purified. Plate-based function screening confirmed each protein was responsible for the assumed activities, further supporting the results of sequence alignment. Xyn526 and Cel1753 showed 79.6 U/mg and 174.0 U/mg specific activities for the hydrolysis of beech wood xylan and CMC, respectively. However, Xyn526 had relatively weak hydrolyzing activity against pNPX, pNPG, and CMC (Table 1). Cel1753 also exhibited hydrolyzing activity for MIC and filter paper, suggesting that it has potential for hydrolyzing cellulose alone. These results revealed that Cel1753 harbored a high proportion of exoglucanase activity relative to its endoglucanase activity (21.8%). This is in contrast to most cellulases, which do not possess both activities (Kim et al., 2008; Ko et al., 2011). Xyl522 could hydrolyze pNPX more effectively than pNPA or pNPG. Bgl8520 could hydrolyze pNPG 50 times more efficiently than pNPC (Table 2), which combined with its relatively low Km for pNPG indicated a higher affinity and efficiency for glycosidic bonds in cellobiose. In addition, Bgl8520 has a higher affinity toward pNPC than other β-glucosidases whose Km are normally between 4.8 and 17.6 mM (Del Pozo et al., 2012), indicating its strong capacity on cellotriose hydrolysis. It is notable that no obvious activities of Xyl522 and Bgl8520 were detected against CMC, MIC, or xylan. Our results therefore suggest that pNPG, xylan from beech wood, pNPX, and CMC were the optimal substrates for Xyl522, Xyn526, Bgl8520, and Cel1753, respectively.
The pH value influenced hydrolysis activity and enzyme stability differently among the four enzymes. The optimal pH values were 4.5, 5.0, 6.0, and 5.0 for Xyl522, Xyn526, Bgl8520, and Cel1753, respectively (Figure 1A). Xyl522 had a narrow pH range, with over 80% hydrolysis activity between pH 4.0 and 5.0, and dramatically lower activity at pH values outside this range. In contrast, Xyn526, Cel1753, and Bgl8520 kept over 80% hydrolysis activity within the pH ranges 4.0–11.0, 3.5–8.0, and 4.5–7.5, respectively. Xyn526 and Cel1753 had high pH stabilities, retaining greater than 80% enzyme activity after incubation at 4°C for 24 h in pH 4.0–12.0 and pH 3.5–8.0, respectively (Figure 1B).
Figure 1. pH optima (A) and stability (B) of Xyl522, Xyn526, Bgl8520, and Cel1753. Activities of Xyl522, Xyn526, Bgl8520, and Cel1753 are shown as solid square, open square, open circle, and solid circle, respectively. Enzymes activities were determined with pNPX, xylan, pNPG, and CMC as substrates for Xyl522 (pH 3.5–7.0), Xyn526 (pH 3.5–11.0), Bgl8520 (pH 4.0–9.0), and Cel1753 (pH 3.5–9.0), respectively. Standard deviations are shown from triplicate measurements.
Xyn526 had the highest optimal temperature of 75°C (Figure 2A); the corresponding values for Xyl522, Bgl8520, and Cel1753 were 60, 65, and 60°C, respectively. These optimal temperatures were generally higher than lignocellulose hydrolases isolated from other metagenomic libraries so far, the majority of which were below 60°C (Duan and Feng, 2010; Yan et al., 2013). In addition, Xyl522, Xyn526, Bgl8520, and Cel1753 had notably high temperature stabilities, retaining over 80% residual activity after incubation at 50°C for 2 h (Figure 2), with thermal half-life being 14, 48, 10, and 32 h, respectively (Supplementary Table 3). Even after 2 h incubation at 70°C, Xyn526 still had 82.3% residual activity toward xylan, whose thermal half-life was impressively 8 h under the same condition. In comparison, the thermostabilities of many other metagenomic enzymes are relatively weak. For example, cellulase En1 from a 40°C biogas digester metagenome (Yan et al., 2013) and glucoside hydrolase Bgxa1 from a dairy cow rumen metagenome (Gruninger et al., 2014) could only retain 80 and 10% activities respectively after incubation at 50°C for 1 h. Xylanase Xyn10N18 from a dairy cow rumen metagenome (Gong et al., 2013) and β-xylosidase RuBG3B from a yak rumen metagenome (Bao et al., 2012) lost all of their activity under the same conditions. Similarly, the alanine, proline, arginine, and glutamic acid content in these proteins were also generally less than those of our four enzymes, while the glycine, serine, lysine, and aspartic acid content was higher. For instance, the percentages of alanine, proline, and arginine residues in cellulase En1 were 4, 2, and 2% less respectively than those in Cel1753; glutamic acid content was the same between the two enzymes. In contrast, the percentages of aspartic acid and lysine residuals in cellulase En1 were 1 and 3% higher respectively than those in Cel1753, although they shared the same percentage of glycine and serine content.
Figure 2. Temperature optima (A) and stability (B–E) of Xyl522, Xyn526, Bgl8520, and Cel1753. (A) Activities of Xyl522, Xyn526, Bgl8520, and Cel1753 are shown as solid square, open square, open circle, and solid circle, respectively. (B–E) Activities of enzymes at 50, 55, 60, 65, 70, and 75°C are shown as solid square, open square, solid circle, open circle, solid diamond, and open diamond, respectively. Enzymes activities were determined with pNPX, xylan, pNPG, and CMC as substrates for Xyl522, Xyn526, Bgl8520, and Cel1753, respectively. Standard deviations are shown from triplicate measurements.
Synergistic Hydrolysis of CMC, Xylan, and Steam-exploded Corncob
Cel1753 could hydrolyze CMC to produce mainly cello-oligosaccharides along with small amounts of glucose and cellobiose, whereas the Bgl8520 could hardly hydrolyze CMC (Figure 3A). When a cocktail of Bgl8520 and Cel1753 was added to CMC, production of glucose increased markedly compared with Cel1753 working alone, whereas the production of cellobiose and other cello-oligosaccharides decreased, suggesting the efficient degradation of CMC. Similarly, the cocktail of Xyn526 and Xyl522 resulted in higher xylose production from beech wood xylan, with xylobiose completely hydrolyzed (Figure 3B).
Figure 3. TLC spots of the four enzymes against CMC (A) and beech wood xylan (B). The enzymes were incubated with substrate at 50°C, pH 5.0 for 4 h. The products were developed by TLC and stained. G1, G2, G3, and G4 are glucose, cellobiose, cellotriose, and cellotetraose, respectively. X1, X2, and X3 are xylose, xylobiose, and xylotriose, respectively. NC: negative control without addition of enzymes.
Synergistic action of these enzymes was also detected when steam-exploded corncob was used as substrate. Compared with Cel1753 alone, addition of Bgl8520 significantly increased the yield of glucose from steam-exploded corncob from 11.8 to 52.3%, while the yield of cellobiose in this cocktail decreased dramatically (Figure 4A). Moreover, glucose production was much quicker in the first 12 h, revealing the synergistic action of the two enzymes and its benefit. Similarly, the supplementation of Xyn526 with Xyl522 resulted in a 56.2% increase in the yield of xylose. This supplementation also simultaneously accelerated the hydrolysis of xylan by consuming xylobiose, which might be attributed to a reduction of Xyn526 inhibition (Figure 4B). When all four enzymes were pooled in a cocktail, xylose and glucose production increased by around 8.3 and 20.5%, respectively, resulting in total hydrolysis efficiency for cellulose and hemicellulose of 71.5 and 87.1% in 48 h at 50°C, respectively (Figure 4). This incremental effect on catalytic efficiency was higher than seen in previous studies of metagenomic enzymes. For example, supplying a metagenomic β-glucosidase into commercial Trichoderma reesei cellulase resulted in a relative increase of 20% glucose yield from corn stover at 50°C (Del Pozo et al., 2012). The mixture of two metagenomic enzymes, β-xylosidase and xylanase, resulted in a relative increase of 37% yield at 40°C from beech wood xylan compared with xylanase alone (Gruninger et al., 2014).
Figure 4. Activities of different enzyme cocktails against steam-exploded corncob. (A) Hydrolysis capacities of Cel1753 (open square: cellubiose production; solid square: glucose production), combination of Bgl8520 and Cel1753 (open circle: cellubiose production; solid circle: glucose production), and the four-enzyme cocktail on steam-exploded corncob (solid diamond: glucose production). (B) Hydrolysis capacities of Xyn526 (solid square: xylobiose production; open square: xylose production), combination of Xyl522 and Xyn526 (solid circle: xylobiose production; open circle: xylose production), and the four-enzyme cocktail on steam-exploded corncob (solid diamond: xylose production). Standard deviations are shown from triplicate measurements.
In the biofuel industry, increasing the thermostability of enzymes is important for reducing operational cost. Many attempts have been made to isolate thermostable enzymes, from organisms such as Clostridium absonum (Rani and Nand, 2000) and Thermoascus aurantiacus (Hong et al., 2007). Because of product inhibition, using of a single enzyme usually leads to low efficiency in the lignocellulose hydrolysis process (Gusakov and Sinitsyn, 1992). However, enzyme cocktails composed of different enzymes can form a “working team”, which is similar as the different enzymes in a single cell (Nie et al., 2014), to sequentially degrade lignocellulose and its intermediates. These enzyme cocktails have proven to be a suitable solution to product inhibition and have therefore drawn worldwide research interest (Chundawat et al., 2011; Gao et al., 2011). To form an efficient cocktail, candidate enzymes need to have compatible pH and temperature optima as well as compatible specific activities. Cocktails have generally been based on commercial enzymes isolated from different microorganisms, which may have different optimal growth conditions. Recently, the isolation of novel enzymes, including thermostable enzymes from metagenomic libraries, has been extensively studied, leading to the isolation of many enzymes with special characteristics such as being cold- (Wierzbicka-Wos et al., 2013) or heat-active (Schroder et al., 2014), and ionic liquid- (Ilmberger et al., 2012) or product-tolerant (Lu et al., 2013). Microorganisms in the same metagenome generally live together for a relatively long time, and should have similar tolerance for environmental conditions, such as temperature, pH, and oxidation-reduction potentials, which may result in isolation of enzymes with similar optimal conditions. Metagenomic screening may therefore offer a way of increasing the compatibility of different candidate enzymes in a cocktail. However, to date there have been few reports of successful trials using a cocktail of enzymes from the same metagenome. Therefore, we identified four novel thermostable enzymes and reported the first trial of a combination of compatible enzymes from the same metagenome.
The four enzymes we isolated from a long-term enriched thermophilic microbial community had higher thermostability (Figure 2) than those from metagenomes sampled from mesophilic soil, sea sediment, and animal gastrointestinal tracts, whose optimal temperatures were generally under 60°C (Duan and Feng, 2010), and whose thermostabilities were relatively poor (Feng et al., 2007; Gong et al., 2012). Although a β-glucosidase has been previously found (from screening a hydrothermal spring metagenome) with an optimal temperature of 90°C (Schroder et al., 2014), its observed relative activity at 50°C was only 20%. Similarly, although one β-glucosidase (LAB25g2 from cow rumen) retains as much as 82% activity after 5 days of incubation at 50°C, it required several days to consume cellobiose completely, suggesting its high thermostability with relatively low hydrolysis efficiency (Del Pozo et al., 2012). We attribute the high thermostabilities of our four enzymes to the fact that they were isolated from a thermophilic methanogenic digester community. This community has been enriched with paper at 53°C for longer than 10 years (Tang et al., 2011), which should have provided a stronger high-temperature selective pressure than that found in the gastrointestinal tracts of ruminants (Yan et al., 2013). Furthermore, the unusual amino acid residue composition of these enzymes may also contribute to their thermostabilities. Lu and colleagues, 1998 reported that there are remarkable differences in amino acids composition between mesophilic and thermophilic proteins. In this research, although the thermostable enzymes we isolated share high amino acid sequence similarities with the mesothermal stable enzymes (Supplementary Figures 5–8), the thermostable enzymes have higher relative ratios of alanine, proline, arginine, and glutamic acid residues and lower ratios of glycine, serine, lysine, and aspartic acid residues (Supplementary Table 4). Alanine tends to form a stable helix, while the reduced number of glycine residues could lessen flexibility of loops and therefore provide heat resistance (Yip et al., 1995). The structural entropy of proline is less than other amino acids as it requires more energy to unfold; its presence can therefore enhance the thermostability of a protein (Vogt et al., 1997). The presence of arginine, which has a large side chain with a delocalized positive charge, can promote the formation of salt bridges and hydrogen bonds to enhance protein stability (Kumar et al., 2000). Lower levels of hydrophilic amino acids such as serine and lysine have also been linked with thermostability (Argos et al., 1979; Lu et al., 1998). Therefore, we propose that the unusual amino acid composition of the four enzymes we isolated did not change their molecular conformation, but did increase their thermostability (Figure 2).
Using a cocktail composed of the four lignocellulose hydrolases, we observed enhanced hydrolysis of steam-exploded corncob after incubation at 50°C for 72 h, suggesting the compatible thermostabilities and synergistic cooperation of these enzymes. It is reasonable that Xyl522 and Bgl8520 could prevent product inhibition of Xyn526 and Cel1753, respectively. As cellulose is often embedded in hemicellulose and lignin (Jeoh et al., 2007), the hydrolytic activity of Xyn526 on hemicellulose could expose more cellulose to Cel1753, along with Xyl522 and Bgl8520 in the four-enzyme cocktail, which could further improve the yield of monosaccharide and stimulate lignocellulose hydrolysis in this community (Figure 4). Indeed, Xyl522 and Xyn526 were cloned from the same fosmid clone (F52), suggesting they might be from the same host organism and potentially hydrolyze xylan synergistically in situ. The presence of a signal peptide in Xyn526 indicates that it is secreted out of cells to hydrolyze xylan. The resulting oligosaccharides, such as xylo-oligosaccharides and arabino-oligosaccharides, could be then transported into the cytoplasm by carbohydrate ABC transporters, which were found to be enriched in F52. Glucoside hydrolases including Xyl522 could then degrade these oligosaccharides into monosaccharides like xylose and arabinose for further metabolism (Supplementary Figure 1). In addition, the extracellular Cel1753 could hydrolyze cellulose outside cells, which might also be transported into the cells with enriched ABC transporters, such as those found in the F85 fosmid clone (Supplementary Figure 2). These could then be hydrolyzed by non-signal peptide β-glucosidases, such as Bgl8520. It is reasonable to assume that such a complex cooperation can only be achieved in a cocktail of enzymes from the same well-enriched metagenome, because all the derived enzymes had similar optimal temperatures, thermostabilities, and pH values, as well as compatible specific activities, and this favorable compatibility was the basis of their synergism and high efficiency. It is notable that our enzyme cocktail was relatively weak in hydrolysis capacity compared with cocktail Accellerase®Trio™, which is composed of cellulase, xylanase, and β-glucosidase, and can hydrolyze 75–90% of glucan and xylan at 50°C in 72 h. In order to improve the catalytic efficiency of our cocktail to reach industrial standards, we also need to use the strategy of developing industrial enzymes, such as molecular modification (Banerjee et al., 2010). In addition, since the original thermophilic methanogenic digester community could nearly completely convert paper to methane within 72 h (Tang et al., 2011), there should be more powerful enzymes, for more powerful enzyme cocktail, still in the dark matter waiting to be explored.
In summary, four novel thermostable glycoside hydrolases were isolated and identified from a metagenomic library constructed from a long-term dry thermophilic methanogenic digester community. The isolated enzymes and their mixed cocktail were highly effective in the hydrolysis of lignocellulose such as steam-exploded corncob. Since enzymes from the same metagenome may have high compatibilities and act synergistically, these experiments offer a way of screening compatible enzymes to form highly efficient enzyme cocktails.
MW, YN, ZS, and XW designed and coordinated the study and write the manuscript. MW, GS, GL, and BZ carried out the experiments and analyzed the results. LL and XW contributed preparing the manuscript. All authors read and approved the final manuscript. ZS, YN, and XW contributed equally to this work.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This study was supported by National Natural Science Foundation of China (31225001), and National High Technology Research and Development Program (“863”Program: 2012AA02A703).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2015.00509/abstract
SHF, separate hydrolysis and fermentation; SSF, simultaneous saccharification and fermentation; GH, glycoside hydrolase; CMC, carboxymethylcellulose; pNPX, p-nitrophenyl-β-D-xyloside; pNPG, p-nitrophenyl-β-D-glucoside; MIC, microcrystalline cellulose; pNPA, pNP-α-L-arabinofuranoside; pNPC, pNP-β-D-cellobioside; ORF, open reading frame; FPA, filter paper activity; TLC, thin-layer chromatography.
Alvarez, T. M., Paiva, J. H., Ruiz, D. M., Cairo, J. P., Pereira, I. O., Paixao, D. A., et al. (2013). Structure and function of a novel cellulase 5 from sugarcane soil metagenome. PLoS ONE 8:e83635. doi: 10.1371/journal.pone.0083635
Arnold, K., Bordoli, L., Kopp, J., and Schwede, T. (2006). The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics 22, 195–201. doi: 10.1093/bioinformatics/bti770
Banerjee, G., Car, S., Scott-Craig, J. S., Borrusch, M. S., Aslam, N., and Walton, J. D. (2010). Synthetic enzyme mixtures for biomass deconstruction: production and optimization of a core set. Biotechnol. Bioeng. 106, 707–720. doi: 10.1002/bit.22741
Bao, L., Huang, Q., Chang, L., Sun, Q., Zhou, J., and Lu, H. (2012). Cloning and characterization of two β-glucosidase/xylosidase enzymes from yak rumen metagenome. Appl. Biochem. Biotechnol. 166, 72–86. doi: 10.1007/s12010-011-9405-x
Chundawat, S. P., Beckham, G. T., Himmel, M. E., and Dale, B. E. (2011). Deconstruction of lignocellulosic biomass to fuels and chemicals. Annu. Rev. Chem. Biomol. 2, 121–145. doi: 10.1146/annurev-chembioeng-061010-114205
Del Pozo, M. V., Fernandez-Arrojo, L., Gil-Martinez, J., Montesinos, A., Chernikova, T. N., Nechitaylo, T. Y., et al. (2012). Microbial β-glucosidases from cow rumen metagenome enhance the saccharification of lignocellulose in combination with commercial cellulase cocktail. Biotechnol. Biofuels. 5:73. doi: 10.1186/1754-6834-5-73
Duan, C. J., Xian, L., Zhao, G. C., Feng, Y., Pang, H., Bai, X. L., et al. (2009). Isolation and partial characterization of novel genes encoding acidic cellulases from metagenomes of buffalo rumens. J. Appl. Microbiol. 107, 245–256. doi: 10.1111/j.1365-2672.2009.04202.x
Fan, Z., Wagschal, K., Lee, C. C., Kong, Q., Shen, K. A., Maiti, I. B., et al. (2009). The construction and characterization of two xylan-degrading chimeric enzymes. Biotechnol. Bioeng. 102, 684–692. doi: 10.1002/bit.22112
Feng, Y., Duan, C. J., Pang, H., Mo, X. C., Wu, C. F., Yu, Y., et al. (2007). Cloning and identification of novel cellulase genes from uncultured microorganisms in rabbit cecum and characterization of the expressed cellulases. Appl. Microbiol. Biotechnol. 75, 319–328. doi: 10.1007/s00253-006-0820-9
Gao, D., Uppugundla, N., Chundawat, S. P., Yu, X., Hermanson, S., Gowda, K., et al. (2011). Hemicellulases and auxiliary enzymes for improved conversion of lignocellulosic biomass to monosaccharides. Biotechnol. Biofuels. 4:5. doi: 10.1186/1754-6834-4-5
Geng, A., Zou, G., Yan, X., Wang, Q., Zhang, J., Liu, F., et al. (2012). Expression and characterization of a novel metagenome-derived cellulase Exo2b and its application to improve cellulase activity in Trichoderma reesei. Appl. Microbiol. Biotechnol. 96, 951–962. doi: 10.1007/s00253-012-3873-y
Gong, X., Gruninger, R. J., Qi, M., Paterson, L., Forster, R. J., Teather, R. M., et al. (2012). Cloning and identification of novel hydrolase genes from a dairy cow rumen metagenomic library and characterization of a cellulase gene. BMC Res. Notes 5:566. doi: 10.1186/1756-0500-5-566
Gong, X., Gruniniger, R. J., Forster, R. J., Teather, R. M., and Mcallister, T. A. (2013). Biochemical analysis of a highly specific, pH stable xylanase gene identified from a bovine rumen-derived metagenomic library. Appl. Microbiol. Biotechnol. 97, 2423–2431. doi: 10.1007/s00253-012-4088-y
Gruninger, R. J., Gong, X., Forster, R. J., and Mcallister, T. A. (2014). Biochemical and kinetic characterization of the multifunctional β-glucosidase/β-xylosidase/α-arabinosidase, Bgxa1. Appl. Microbiol. Biotechnol. 98, 3003–3012. doi: 10.1007/s00253-013-5191-4
Gusakov, A. V., Salanovich, T. N., Antonov, A. I., Ustinov, B. B., Okunev, O. N., Burlingame, R., et al. (2007). Design of highly efficient cellulase mixtures for enzymatic hydrolysis of cellulose. Biotechnol. Bioeng. 97, 1028–1038. doi: 10.1002/bit.21329
Gusakov, A. V., and Sinitsyn, A. P. (1992). A theoretical analysis of cellulase product inhibition: effect of cellulase binding constant, enzyme/substrate ratio, and β-glucosidase activity on the inhibition pattern. Biotechnol. Bioeng. 40, 663–671. doi: 10.1002/bit.260400604
Han, Q., Robinson, H., Cao, L., Qian, C., Wang, Q., Xie, L., et al. (2013). Biochemical characterization and crystal structure of a GH10 xylanase from termite gut bacteria reveal a novel structural feature and significance of its bacterial Ig-like domain. Biotechnol. Bioeng. 110, 3093–3103. doi: 10.1002/bit.24982
Hess, M., Sczyrba, A., Egan, R., Kim, T. W., Chokhawala, H., Schroth, G., et al. (2011). Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science 331, 463–467. doi: 10.1126/science.1200387
Hong, J., Tamaki, H., and Kumagai, H. (2007). Cloning and functional expression of thermostable β-glucosidase gene from Thermoascus aurantiacus. Appl. Microbiol. Biotechnol. 73, 1331–1339. doi: 10.1007/s00253-006-0618-9
Hu, Y., Fu, C., Yin, Y., Cheng, G., Lei, F., Yang, X., et al. (2010). Construction and preliminary analysis of a deep-sea sediment metagenomic fosmid library from Qiongdongnan Basin, South China Sea. Mar. Biotechnol. 12, 719–727. doi: 10.1007/s10126-010-9259-1
Ilmberger, N., Meske, D., Juergensen, J., Schulte, M., Barthen, P., Rabausch, U., et al. (2012). Metagenomic cellulases highly tolerant towards the presence of ionic liquids-linking thermostability and halotolerance. Appl. Microbiol. Biotechnol. 95, 135–146. doi: 10.1007/s00253-011-3732-2
Jeoh, T., Ishizawa, C. I., Davis, M. F., Himmel, M. E., Adney, W. S., and Johnson, D. K. (2007). Cellulase digestibility of pretreated biomass is limited by cellulose accessibility. Biotechnol. Bioeng. 98, 112–122. doi: 10.1002/bit.21408
Kim, S. J., Lee, C. M., Han, B. R., Kim, M. Y., Yeo, Y. S., Yoon, S. H., et al. (2008). Characterization of a gene encoding cellulase from uncultured soil bacteria. FEMS Microbiol. Lett. 282, 44–51. doi: 10.1111/j.1574-6968.2008.01097.x
Ko, K.-C., Han, Y., Choi, J. H., Kim, G.-J., Lee, S.-G., and Song, J. J. (2011). A novel bifunctional endo-/exo-type cellulase from an anaerobic ruminal bacterium. Appl. Microbiol. Biotechnol. 89, 1453–1462. doi: 10.1007/s00253-010-2949-9
Lu, J., Du, L., Wei, Y., Hu, Y., and Huang, R. (2013). Expression and characterization of a novel highly glucose-tolerant β-glucosidase from a soil metagenome. Acta Biochim. Biophys. Sin. 45, 664–673. doi: 10.1093/abbs/gmt061
Magnus, A., Kim, O., Tommaso, D. F., Laura, R., Merja, P., Gunnar, L., et al. (2012). Challenges in enzymatic hydrolysis and fermentation of pretreated Arundo donax revealed by a comparison between SHF and SSF. Process Biochem. 47, 1452–1459. doi: 10.1016/j.procbio.2012.05.016
Mcandrew, R. P., Park, J. I., Heins, R. A., Reindl, W., Friedland, G. D., D'haeseleer, P., et al. (2013). From soil to structure, a novel dimeric β-glucosidase belonging to glycoside hydrolase family 3 isolated from compost using metagenomic analysis. J. Biol. Chem. 288, 14985–14992. doi: 10.1074/jbc.M113.458356
Morana, A., Paris, O., Maurelli, L., Rossi, M., and Cannio, R. (2007). Gene cloning and expression in Escherichia coli of a bi-functional β-D-xylosidase/α-L-arabinosidase from Sulfolobus solfataricus involved in xylan degradation. Extremophiles 11, 123–132. doi: 10.1007/s00792-006-0020-7
Murashima, K., Kosugi, A., and Doi, R. H. (2003). Synergistic effects of cellulosomal xylanase and cellulases from Clostridium cellulovorans on plant cell wall degradation. J. Bacteriol. 185, 1518–1524. doi: 10.1128/JB.185.5.1518-1524.2003
Nie, Y., Liang, J. L., Fang, H., Tang, Y. Q., and Wu, X. L. (2014). Characterization of a CYP153 alkane hydroxylase gene in a Gram-positive Dietzia sp. DQ12-45-1b and its “team role” with alkW1 in alkane degradation. Appl. Microbiol. Biotechnol. 98, 163–173. doi: 10.1007/s00253-013-4821-1
Poutanen, K., and Puls, J. (1988). Characteristics of Trichoderma reesei β-xylosidase and its use in the hydrolysis of solubilized xylans. Appl. Microbiol. Biotechnol. 28, 425–432. doi: 10.1007/BF00268208
Rinke, C., Schwientek, P., Sczyrba, A., Ivanova, N. N., Anderson, I. J., Cheng, J. F., et al. (2013). Insights into the phylogeny and coding potential of microbial dark matter. Nature 499, 431–437. doi: 10.1038/nature12352
Schroder, C., Elleuche, S., Blank, S., and Antranikian, G. (2014). Characterization of a heat-active archaeal β-glucosidase from a hydrothermal spring metagenome. Enzyme Microb. Technol. 57, 48–54. doi: 10.1016/j.enzmictec.2014.01.010
Tang, Y. Q., Ji, P., Hayashi, J., Koike, Y., Wu, X. L., and Kida, K. (2011). Characteristic microbial community of a dry thermophilic methanogenic digester: its long-term stability and change with feeding. Appl. Microbiol. Biotechnol. 91, 1447–1461. doi: 10.1007/s00253-011-3479-9
Teng, C., Yan, Q., Jiang, Z., Fan, G., and Shi, B. (2010). Production of xylooligosaccharides from the steam explosion liquor of corncobs coupled with enzymatic hydrolysis using a thermostable xylanase. Bioresour. Technol. 101, 7679–7682. doi: 10.1016/j.biortech.2010.05.004
Tomas-Pejo, E., Oliva, J. M., Ballesteros, M., and Olsson, L. (2008). Comparison of SHF and SSF processes from steam-exploded wheat straw for ethanol production by xylose-fermenting and robust glucose-fermenting Saccharomyces cerevisiae strains. Biotechnol. Bioeng. 100, 1122–1131. doi: 10.1002/bit.21849
Warnecke, F., Luginbuhl, P., Ivanova, N., Ghassemian, M., Richardson, T. H., Stege, J. T., et al. (2007). Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite. Nature 450, 560–565. doi: 10.1038/nature06269
Wierzbicka-Wos, A., Bartasun, P., Cieslinski, H., and Kur, J. (2013). Cloning and characterization of a novel cold-active glycoside hydrolase family 1 enzyme with β-glucosidase, β-fucosidase and β-galactosidase activities. BMC Biotechnol. 13:22. doi: 10.1186/1472-6750-13-22
Wood, P. J., Erfle, J. D., and Teather, R. M. (1988). Use of complex formation between Congo Red and polysaccharides in detection and assay of polysaccharide hydrolases. Method. Enzymol. 160, 59–74. doi: 10.1016/0076-6879(88)60107-8
Yan, X., Geng, A., Zhang, J., Wei, Y., Zhang, L., Qian, C., et al. (2013). Discovery of (hemi-) cellulase genes in a metagenomic library from a biogas digester using 454 pyrosequencing. Appl. Microbiol. Biotechnol. 97, 8173–8182. doi: 10.1007/s00253-013-4927-5
Yip, K., Stillman, T., Britton, K., Artymiuk, P., Baker, P., Sedelnikova, S., et al. (1995). The structure of Pyrococcus furiosus glutamate dehydrogenase reveals a key role for ion-pair networks in maintaining enzyme stability at extreme temperatures. Structure 3, 1147–1158. doi: 10.1016/S0969-2126(01)00251-9
Zhang, J., Zhong, Y., Zhao, X., and Wang, T. (2010). Development of the cellulolytic fungus Trichoderma reesei strain with enhanced β-glucosidase and filter paper activity using strong artifical cellobiohydrolase 1 promoter. Bioresour. Technol. 101, 9815–9818. doi: 10.1016/j.biortech.2010.07.078
Zhang, Y. H. P., and Lynd, L. R. (2004). Toward an aggregated understanding of enzymatic hydrolysis of cellulose: noncomplexed cellulase systems. Biotechnol. Bioeng. 88, 797–824. doi: 10.1002/bit.20282
Zhou, J., Wang, Y. H., Chu, J., Luo, L. Z., Zhuang, Y. P., and Zhang, S. L. (2009). Optimization of cellulase mixture for efficient hydrolysis of steam-exploded corn stover by statistically designed experiments. Bioresour. Technol. 100, 819–825. doi: 10.1016/j.biortech.2008.06.068
Keywords: cellulase, xylanase, β-xylosidase, β-glucosidase, metagenome, enzyme compatibility, enzyme cocktail, biofuels production
Citation: Wang M, Lai G-L, Nie Y, Geng S, Liu L, Zhu B, Shi Z and Wu X-L (2015) Synergistic function of four novel thermostable glycoside hydrolases from a long-term enriched thermophilic methanogenic digester. Front. Microbiol. 6:509. doi: 10.3389/fmicb.2015.00509
Received: 03 April 2015; Accepted: 08 May 2015;
Published: 22 May 2015.
Edited by:Cheng Zhong, Tianjin University of Science and Technology, China
Reviewed by:Amy Michele Grunden, North Carolina State University, USA
Pravin Malla Shrestha, University of California, Berkeley, USA
Copyright © 2015 Wang, Lai, Nie, Geng, Liu, Zhu, Shi and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yong Nie and Xiao-Lei Wu, College of Engineering, Peking University, Yiheyuan Road, Beijing 100871, China, email@example.com; firstname.lastname@example.org;
Zhongping Shi, School of Biotechnology, Jiangnan University, Lihu Road, Wuxi 214122, China, email@example.com