Insight into the functional roles of Glu175 in the hyperthermostable xylanase XYL10C-ΔN through structural analysis and site-saturation mutagenesis
Biotechnology for Biofuels volume 11, Article number: 159 (2018)
Improving the hydrolytic performance of hemicellulases to degrade lignocellulosic biomass is of considerable importance for second-generation biorefinery. Xylanase, as the crucial hemicellulase, must be thermostable and have high activity for its potential use in the bioethanol industry. To obtain excellent xylanase candidates, it is necessary to understand the structure–function relationships to provide a meaningful reference to improve the enzyme properties. This study aimed to investigate the catalytic mechanism of a highly active hyperthermophilic xylanase variant, XYL10C-ΔN, for hemicellulose degradation.
By removing the N-terminal 66 amino acids, the variant XYL10C-ΔN showed a 1.8-fold improvement in catalytic efficiency and could hydrolyze corn stover more efficiently in hydrolysis of corn stover; however, it retained similar thermostability to the wild-type XYL10C. Based on the crystal structures of XYL10C-ΔN and its complex with xylobiose, Glu175 located on loop 3 was found to be specific to GH10 xylanases and probably accounts for the excellent enzyme properties by interacting with Lys135 and Met137 on loop 2. Site-saturation mutagenesis confirmed that XYL10C-ΔN with glutamate acid at position 175 had the highest catalytic efficiency, specific activity, and the broadest pH-activity profile. The functional roles of Glu175 were also verified in the mutants of another two GH10 xylanases, XylE and XynE2, which showed increased catalytic efficiencies and wider pH-activity profiles.
XYL10C-ΔN, with excellent thermostability, high catalytic efficiency, and great lignocellulose-degrading capability, is a valuable candidate xylanase for the biofuel industry. The mechanism underlying improved activity of XYN10C-ΔN was thus investigated through structural analysis and functional verification, and Glu175 was identified to play the key role in the improved catalytic efficiency. This study revealed the importance of a key residue (Glu175) in XYN10C-ΔN and provides a reference to modify GH10 xylanases for improved catalytic performance.
The rising global demand for renewable fuels and the climate change associated with the consumption of fossil fuels are considered as the biggest social challenges over the next few decades [1, 2]. Bioethanol produced from polysaccharides such as starch can substitute for fossil fuel and has been identified as the first-generation biofuel; however, its use so far has been restricted as it is not sustainable. For second-generation biofuels, renewable resources including wood, grass, agricultural and forest residues, and fiber sludge are widely used for bioethanol production [3,4,5]. Biomass mainly comprises cellulose, hemicellulose and lignin, in which the cellulose fibers are embedded in a matrix of hemicellulose and lignin. Therefore, to convert biomass into fermentable sugars efficiently, the synergistic actions of hemicellulases and cellulases are very important [6,7,8].
Hemicellulose is the second most abundant polysaccharide on earth. Efficient conversion of hemicellulose to other value-added fermentation products has both environmental and economic benefits . Xylan, the major component of hemicellulose, consists of a backbone of xylose units, linked via β-1,4-glycosidic bonds and some branched chains. Owing to the structural complexity, complete degradation of xylan requires a series of enzymes. Among them, xylanase (endo-1,4-d-xylanase, EC 188.8.131.52) is a key enzyme that randomly cleaves the internal β-1,4-d-xylosidic linkages of xylan to yield different chain lengths of xylo-oligosaccharides . Based on their amino acid sequence similarities, three-dimensional (3-D) structures, and hydrophobic cluster analysis, xylanases are usually classified into glycoside hydrolase (GH) families 10 and 11 . GH10 xylanases have a typical (β/α)8 barrel fold structure (also known as the typical triosephosphate isomerase barrel) with catalytic residues at the end of β-sheets 4 and 7 [10, 11], while GH11 xylanases display a β-jelly roll structure and employ the same retaining mechanism as GH10 xylanases [12, 13].
Currently, the commercial uses of xylanases are mainly in the paper, food, and animal feed industries [9, 14], however, xylanases are increasingly being recognized as important enzymes in biorefining of lignocellulosic biomass [15, 16]. In the production of bioethanol from lignocellulosic materials, xylanase can promote the hydrolysis of cellulose by degrading xylan and increasing the access of cellulase to the cellulose surface . GH10 xylanases have a broad substrate specificity and, in combination with cellulases, degrade wheat straw and corn fiber very efficiently [18,19,20]. The poor catalytic performance of most wild-type xylanases restricts their application in bioprocesses ; therefore, it is desirable to determine the catalytic mechanism of xylanase and improve its properties by protein engineering. Some structural elements, such as local loops and terminal regions, and key residues involved in the formation of protein–protein interactions (e.g., hydrogen bonds, salt bridges, and sigma π), fold on triggered trimerization, which affects the catalytic process of proteins, making them potential foci for protein engineering [22,23,24,25,26,27]. Among these factors, some key residues in the active loops and the N-terminus play certain roles in the catalytic efficiency of xylanases [25, 26]. Their functions have been verified by site-directed and saturation mutagenesis, based on the simple but powerful Darwinian evolution principles of mutation and selection [28, 29].
To meet the rigorous requirements of industrial processes, and to decrease the production-cost, highly active xylanases with excellent thermostability are of great interest. The GH10 xylanase, XYL10C, from acidophilic Bispora sp. MEY-1 is an attractive candidate industrial enzyme that has remarkable activity and favorable thermostability . The present study aimed to construct a truncated mutant, XYL10C-ΔN, and compare its enzymatic properties with those of the wild-type to reveal the underlying mechanisms of its high catalytic efficiency, and to form the basis for improving other GH10 xylanases.
Characterization of the N-terminus-truncated mutant XYL10C-ΔN
Based on the multiple sequence alignment of GH10 xylanases (Additional file 1: Figure S1), XYL10C was observed to harbor an extra N-terminal sequence of 66 amino acid residues (except for the 18-residue signal peptide). The mutant version of XYL10C with the N-terminal sequence deleted, XYL10C-ΔN, was expressed in Pichia pastoris and purified (Additional file 1: Figure S2). In comparison with the wild-type XYL10C (optimally active at 85 °C and pH 4.5), XYL10C-ΔN showed a slight downshift of 5 °C and 0.5 in the optimal temperature and pH (Table 1). The two enzymes had similar kinetic stabilities at 80 °C and 90 °C, but different thermodynamic stabilities (Additional file 1: Figure S3). Compared with XYL10C, XYL10C-ΔN had a decreased Tm value of 73.8 °C (vs. 82.0 °C). The enzymes had similar pH stability (Additional file 1: Figure S4) and employed an endo-mode of action to cleave xylotriose, xylotetraose, xylopentaose, xylohexaose, and beechwood xylan into mainly xylose and xylobiose (Additional file 1: Figure S5). Significant differences were observed for their substrate affinities and catalytic performances. XYL10C-ΔN had higher Km, kcat/Km (catalytic efficiency), and specific activity values (Table 1), which were 1.3-, 2.7-, and 1.8-fold higher than those of XYL10C, respectively.
Analysis of the crystal structures of XYL10C-ΔN and its complex
The crystal structure of XYL10C-ΔN was solved at a high resolution of 1.6 Å by molecular replacement using the structure of XylE (PDB code: 4F8X) from Penicillium scopiformis  as a search model. The structure determination procedures are listed in Additional file 1: Table S1. Based on the analysis of the X-ray diffraction pattern (Additional file 1: Figure S6), the XYL10C-ΔN crystals belonged to the monoclinic space group of C2. As shown in Additional file 1: Figure S7, there are two monomers, which have the topology of a single (β/α)8 barrel typical of GH10 xylanases, and are denoted as chain A and B in the asymmetric unit. Strong protein interactions were detected between chain A and B, including hydrogen bonds and van der Waals forces. These tightly fixed subunits probably have a better ability to bind the substrate at the catalytic pocket and maintain the stability of the monomer.
To fully understand the substrate binding and catalysis mechanisms, co-crystallization of XYL10C-ΔN and xylobiose was also completed at a resolution of 1.7 Å. As shown in Fig. 1, there is one xylobiose molecule in the open and extended active-site groove of XYL10C-ΔN, and two β-1,4-linked xylosyl moieties in subsites − 1 and − 2. Superimposition of the XYL10C-ΔN-xylobiose complex structure and the XYL10C-ΔN structure yielded a root-mean-square deviation (RMSD) value of 0.24 Å for all Cα atoms, suggesting that ligand binding did not alter the protein conformation. The xylobiose can form several hydrogen bonds, such as with Ala131, Asn132, His168, Asn218, Glu300, and Trp372 around the catalytic pocket, which provide strong binding forces to stabilize the enzyme–substrate complex and to maintain effective catalysis. Among them, Ala131, Lys135, and Met138 are located on loop 2, which is the active loop involved in substrate binding and catalysis. Thus, the amino acid residues that can affect the swing of loop 2 will affect the catalytic efficiency of the enzyme. By docking with long-chain xyloheptaose (Additional file 1: Figure S8), Asn269, Tyr272, His302, and Arg336 were also found to interact with the substrate. These residues are conserved throughout evolution and indispensable for substrate binding and catalysis.
Identification of the key residue
As shown in the structure of the XYL10C-ΔN-xylobiose complex (Fig. 1c), the side chain N atom of Lys135, located on the loop 2, forms three hydrogen bonds with the hydroxyl oxygen, glycosidic oxygen, and O5 of xylobiose. These interactions limit the swing of loop 2, and further affect the interactions of the substrate and other functional amino acids. The OE2 of Glu175 shows a distance of approximately 2.9 Å to the side chain N atom of Lys135, which permits the formation of salt and hydrogen bonds. Moreover, the carboxy oxygen of Glu175 can form a hydrogen bond with the nitrogen of Met137 located on loop 2. In addition, loop 3, located in the catalytic pocket, would be immobilized because of the counteracting force. The strictly conserved Trp172 in loop 3 can form a hydrogen bond with the catalytic Glu219 and exerts a hydrophobic force with xylopyranose. Glu175 is located in the center of this powerful interaction network and contributes to the stable conformation of loop 3, the loop responsible for the rapid catalytic reaction. Multiple sequence alignment of 50 GH10 xylanases indicated that four proteins have Glu at this site (Table 2). Thus, residue Glu175 of XYL10C-ΔN, which interacts strongly to Lys135 and Met137 of loop 2, and is in close proximity to loop 3, might have the roles of stabilizing the substrate conformation in the catalytic center and affecting the catalytic efficiency of XYL10C-ΔN. To verify the functional roles of Glu175 in the catalysis of GH10 xylanases, site-saturation mutagenesis was performed.
Characterization of XYL10C-ΔN mutants
XYL10C-ΔN mutants with Glu175 replacement were successfully constructed, expressed, and purified (Additional file 1: Figure S2). Compared with XYL10C-ΔN, the variant enzymes showed a temperature optimum change of − 5 to + 15 °C and a pH optimum change of − 0.5 to + 1.0 (Additional file 1: Table S2). Moreover, the mutants had narrower pH-activity profiles (Fig. 2a). Using beechwood xylan as the substrate, the mutants had increased Km values (1.1–4.8-fold) and decreased kcat/Km values (1.0–37.3 times), as well as decreased specific activities (0.8–32.5 times) (Fig. 3a and Additional file 1: Table S3). The catalytic efficiency of XYL10C-ΔN was significantly different from that of all other site-mutated variants (P < 0.05). The Km values of XYL10C-ΔN and XYL10C-ΔN-E175N were very similar, but were significantly different from those of other variants. The results indicated that residue substitution at position 175 probably altered the intramolecular interaction, the protein structure, or the substrate conformation, which consequently changed the enzyme’s properties.
Functional validation of Glu175 in other GH10 xylanases
Two xylanases, XylE and XynE2, were selected to verify the roles of Glu175 in the catalytic performance and pH adaptability of GH10 xylanases. Compared with the wild-type proteins, residue substitution caused no changes in the temperature and pH optima of most mutants (65–70 °C and 5.0–5.5; Additional file 1: Table S4); however, mutants XylE-Q116E and XynE2-Q85E showed broader pH-activity profiles (Fig. 2b, c), increased catalytic efficiencies (by 1.2 and 9.2-fold, respectively), and specific activities (by 1.1 and 7.7-fold, respectively) (Table 1, Additional file 1: Table S5 and Fig. 3b). The catalytic efficiency and Km value of XylE-Q16E were different from those of all other variants (P < 0.05). These results confirmed the functional roles of Glu175 in the catalytic performance and pH-activity profile of GH10 xylanases. However, the importance of Glu175 in the catalytic efficiency of other GH10 xylanases remains to be determined.
Enzymatic hydrolysis of corn stover
Steam-exploded corn stover mainly comprises cellulose and hemicellulose. When using steam-exploded corn stover as the substrate, the hydrolytic capabilities of commercial cellulase from T. reesei and its combination with XYL10C or XYL10C-ΔN were assessed at pH 5.0 and 50 °C with agitation for various times. As shown in Fig. 4, 100 U of T. reesei cellulase alone released 8.0 µmol of reducing sugars during 24 h of incubation. When combined with XYL10C or XYL10C-ΔN (30 U of xylanase and 70 U of cellulase), the reducing sugar amounts released were 9.1 and 9.4 µmol, and during the first 2 h of reaction, the sugar production rates increased by 2 and 1.6 times, respectively. The results indicated that T. reesei cellulase and XYL10C or XYL10C-ΔN have a synergistic activity to degrade steam pretreated corn stover, and that XYL10C-ΔN is more efficient than the wild-type for the hydrolysis of corn stover.
The study of the relationship between enzyme structure and function provides an important reference for efficient protein engineering . Artificial enzyme engineering is a powerful strategy to pinpoint functional determinants and to rapidly improve enzyme fitness with regard to its physical or biochemical properties [33,34,35,36]. Thermostable xylanases with high catalytic efficiency are highly desirable for many industrial processes. In our previous study, a hyperthermophilic xylanase of GH10, XYL10C, from Bispora sp. MEY-1 showed great application potential in the bioethanol industry . Multiple sequence alignment indicated that XYL10C contains an extra N-terminal sequence of 66 amino acid residues. N-terminal motifs play an important role in protein function, such as binding, thermostability, and specific activity [37,38,39]. In the present study, after removal of this N-terminal sequence, the variant XYL10C-ΔN retained its hyperthermophilic characteristic but showed a 1.8- and 2.7-fold improvement in catalytic efficiency and specific activity, respectively. XYL10C-ΔN shares high-sequence identity with the Talaromyces leycettanus xylanase TlXyn10A (54%)  and Penicillium canescens xylanase XylE (53%) , but has a higher specific activity (8700 vs. 2240 and 50 U/mg) and catalytic efficiency (8800 vs. 1626 and 92 mL/s/mg). Moreover, XYL10C-ΔN exhibits similar thermostability to TlXyn10A (no enzyme activity loss after 1 h-incubation at 80 °C), which is better than XylE (completely inactive after 30 min incubation at 70 °C). Compared with other GH10 xylanases, such as reAuXyn10A from Aspergillus usamii (5448 mg/mL) , the H179F mutant of XynA (1030 mL/s/mg) from Geobacillus stearothermophilus , MpXyn10A from Malbranchea pulchella , and GtXyn10 from Gloeophyllum trabeum , XYL10C-ΔN also has a superior catalytic performance and thermostability. Therefore, it is important to reveal the underlying mechanism of highly active XYL10C-ΔN for industrial purposes.
GH10 xylanases have conserved amino acids located in the shallow slot around the top face of the protein molecule, including two catalytic glutamates and five to seven substrate binding subsites [45,46,47,48]. Among them, the highly conserved subsites − 1, − 2, and + 1  have several diverse interactions with xylose moieties, and thus play a critical role in substrate recognition and binding. Structural analysis of the XYL10C-ΔN-xylobiose complex indicated five key residues (Asn132, Lys135, His168, Asn218, and Trp372) that are involved in substrate binding via hydrogen bonds. However, ligand binding in the aglycon region of the substrate binding cleft is determined by the hydrophobic interactions between aromatic amino acids and xylobiose rings. Despite the conserved sequence and structure near the glycon subsites, the proximal glycon subsites may display variations. For example, Cellvibrio japonicus xylanase 10C, having an insertion of Tyr and a substitution of Glu/Gly at subsites − 3 and − 2, showed very low activity on less polymerized substrates . XYL10C-ΔN and XylE, with 53% sequence similarity, share similar three-dimensional structures; however, their enzyme properties differ considerably [30, 31]. As shown in the structural alignment of XYL10C-ΔN and XylE (Fig. 1d), there are several important but different residues in the catalytic pockets. Among them, Glu175 of XYL10C-ΔN on loop 3, which is in close proximity to many important substrate binding sites and the substrate, might play a part in the formation of the catalytic pocket. In particular, it can form a variety of force interactions with Lys135 on loop 2, the loop that interacts with the substrate to stabilize the conformation. The loop structures, especially the flexible ones, are involved in substrate binding, enzyme catalysis, and product release of triose-phosphate isomerase (TIM)-(β/α)8 enzymes [49, 50]. The residue at position 175 is highly conserved in evolution. Most known GH10 xylanases have glutamine at this site, while XYL10C-ΔN has glutamate instead. Site-saturation mutagenesis of XYL10C-ΔN indicated that this site is closely associated with the catalytic efficiency, and Glu, Met, Ser, and Gln were the dominant residues for the catalytic efficiency and specific activity.
Xylanases with high specific activity, good catalytic efficiency, and excellent stability are required for biorefinery uses. Recent studies have shown that the longer the xylo-oligosaccharides, the greater the extent of cellulase inhibition . In addition, cellulases stuck on cellulose microfibrils during enzymatic hydrolysis show a decreased degradation rate . Xylanase can assist in the release of cellulases stuck on the substrate  and relieve xylo-oligosaccharide-induced cellulase inhibition [54, 55]. To date, most industrial processes have used GH11 xylanases for biorefining purposes . The XYL10C-ΔN produced in this study has the potential to be applied for biorefining. It not only shows a synergistic activity with cellulase to degrade corn stover, but also enhances the reducing sugar-producing rate and increases the total amount of reducing sugar. We inferred that XYL10C-ΔN could remove xylan obstacles that impede the processive movement of cellulases, involving cleaving xylan into smaller chain products, which would ultimately promote corn stover degradation [17, 57]. Therefore, XYL10C-ΔN has great potential for application in biomass degradation and biofuels.
In the present study, we obtained a xylanase variant XYL10C-ΔN with excellent thermostability, high catalytic efficiency, and synergistic action with cellulase in the degradation of straw. A key residue on loop 3, Glu175, was found to interact with Lys135 and Met137 of loop 2, and could account for the improved catalytic efficiency and broadened pH-activity profile. This study not only provides an excellent candidate xylanase for potential industrial use, but also provides the theoretical and practical basis for the modification of GH10 xylanases.
The gene fragment coding for the truncated XYL10C (GenBank Accession Number: FJ492963) without the N-terminal 66 amino acid residues, designated xyl10c-ΔN, was amplified by PCR with a primer set (5′-GAATTCTGGGGTCTTAATAATGCAGCTCGAGCCG-3′ and 5′-GCGGCCGCTCATGGACTTTCCGCCTTATGTTGCAAAGCCTG-3′, containing the EcoRI and NotI restriction sites). The gene fragment was then cloned into the expression vector pPIC9 at the EcoRI and NotI restriction sites using the T4 ligase. The recombinant plasmid pPIC9-xyl10c-ΔN was linearized by BglII and transformed into P. pastoris GS115 competent cells using a Gene Pulser X cell Electroporation System (Bio-Rad). Recombinant expression and fermentation were conducted following the method previously described .
The culture supernatants were collected by centrifugation (12,000×g, 4 °C, and 10 min), followed by concentration through a Vivaflow ultrafiltration membrane (Vivascience) with a molecular weight cut-off of 10 kDa. The crude enzymes were loaded onto a FPLC HiTrap Q Sepharose XL 5 mL column (GE Healthcare) that was equilibrated with 20 mM McIlvaine buffer (pH 6.5). Enzymes were eluted utilizing a linear gradient of NaCl (1.0 M) in the same buffer at a flow rate of 4.0 mL/min. Fractions exhibiting xylanase activities were collected and subjected to sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). The protein concentration was determined using the Bradford method with albumin from bovine serum as a standard. To remove N-glycans, approximately 600 mg of purified recombinant XYL10C-ΔN was treated by endo-β-N-acetylglucosaminidase H (Endo Hf, New England Biolabs) in 50 mM sodium citrate (pH 5.6) at 300 K for 1 day. The reaction mixture was then dialyzed twice against 50 mM sodium citrate (pH 6.0) containing 20 mM MES, followed by elution through a DEAE column in 200 mM NaCl. The eluted XYL10C-ΔN was re-dialyzed in 50 mM sodium citrate (pH 7.5) containing 25 mM Tris–HCl, 150 mM NaCl, and concentrated to 66 mg/mL using Amicon Ultra-15 Centrifugal Filter Units (Millipore). The protein purity was higher than 95% as shown by SDS-PAGE.
Crystallization and data collection
The initial crystallization screen was performed manually using the 768 different reservoir conditions from Hampton Research kits (Laguna Niguel) and sitting-drop vapor-diffusion method. Briefly, 2 μL of XYL10C-ΔN solution (25 mM Tris–HCl [pH 7.5] and 150 mM NaCl; 66 mg/mL) was mixed with 2 μL of reservoir solution in 24-well Cryschem Plates (Hampton Research), and equilibrated against 300 μL of reservoir solution at 298.15 K. The initial crystals of XYL10C-ΔN were obtained within 5 days using the Index condition No. 83 (0.2 M MgCl2·6H2O, 0.1 M Bis–Tris [pH 6.5], and 25% [w/v] PEG 3350). The crystallization conditions were then optimized to 0.4 M MgCl2·6H2O, 0.1 M Bis–Tris (pH 6.5), and 19% (w/v) PEG 3350. At day 5, the crystals reached to a size of approximately 1.0 mm × 0.4 mm × 0.3 mm. The crystals that diffracted to 1.5 Å resolution were cryo-protected with a solution containing 0.5 M MgCl2·6H2O, 0.15 M Bis–Tris (pH 5.5), 25% (w/v) polyethylene glycol 3350 and 10% (w/v) glycerol. Complete X-ray intensity data were collected from a single crystal at beam line BL15A1 of the National Synchrotron Radiation Research Center of China. A 0.5° oscillation angle was used, the exposure time was 3 s, and the distance between crystal and detector was 200 mm. The diffraction images were processed using the program HKL-2000 . The XYL10C-ΔN structure was solved using the molecular replacement (MR) method of Phaser program  in the CCP4 suite with the endo-1,4-β-xylanase XylE from Phialocephala scopiformis (53% sequence identity) as a search model (PDB code: 4F8X). Initial structure refinement using REFMAC5 resulted in a model with the Rwork and Rfree values of 0.125 and 0.175, respectively.
Site-saturation mutagenesis and enzyme production
Based on the structure analysis and sequence evolution analysis of 50 GH10 xylanases (Table 2), Glu175 of XYL10C-ΔN was found to be special and probably play a key role in the efficient catalysis. Overlap PCR was then performed to substitute Glu175 of XYL10C-ΔN with other amino acid residues. The specific primers are shown in Additional file 1: Table S6. Construction of recombinant plasmids and heterologous expression and purification followed the same procedures as described above.
Biochemical characterization of the XYL10C-ΔN and its mutants
Beechwood xylan from Sigma was used as the substrate. The pH-activity profiles of the wild-type and variant enzymes were determined by measuring the xylanase activity after incubation at 80°C in McIlvaine buffer (pH 2.5–8.5) and glycine–NaOH (9.0–11.0) containing 1.0% beechwood xylan for 10 min. To determine the optimal temperature, the enzyme activities were measured at temperatures ranging from 30 to 95 °C in McIlvaine buffer of each optimal pH for 10 min. The amounts of reducing sugar released were determined using the 3,5-dinitrosalicylic acid (DNS) method . Thermal stability was determined by measuring the residual activities under standard conditions (pH 4.0, 80 °C, and 10 min) after incubation of the enzymes at 80 and 90 °C for 0, 2, 5, 10, 20, 30, and 60 min. The melting temperature (Tm, thermodynamic stabilities) of XYL10C and XYL10C-ΔN were determined on a MicroCal™ VP-Capillary DSC (GE Healthcare). Each protein sample, 200 μg, were dissolved in 0.5 mL of 20-mM McIlvaine buffer (pH 6.5). The degassed protein samples and controls were treated with a heating rate of 120 °C/h over the temperature range of 20–110 °C. Each experiment was repeated at least three times.
Analysis of the hydrolysis products
The hydrolysis products of xylo-oligosaccharides and beechwood xylan by XYL10C and XYL10C-ΔN were determined using the high-performance anion-exchange chromatography (HPAEC; Thermo Fisher Scientific, Sunnyvale, CA) equipped with a CarboPac PA200 guard column (3 mm × 250 mm). Each enzyme (0.2 U) was added into the 0.4 mL 50 mM Na2HPO4-citric acid (pH 4.5) containing 2 μg/mL of xylobiose to xylohexaose or 1% (w/v) beechwood xylan, and incubated at 70 °C for 16 h. The hydrolysis products were resolved in a mobile phase of 100 mM NaOH using xylose to xylohexaose as standards.
Determination of the kinetic values
Beechwood xylan of 1–10 mg/mL in 0.1 M Na2HPO4-citric acid buffer was used as the substrate for determination of the Km, Vmax and kcat values. The enzymatic activities of purified recombinant XYL10C-ΔN and its mutants were assayed under standard conditions (80 °C and pH 4.0) for 5 min. The data were plotted and calculated according to the Lineweaver–Burk method. The catalytic efficiency (kcat/Km) of each enzyme was calculated for comparison. Significant differences of the Km and kcat/Km values were tested using the ANOVA of SPSS 15.0 at the P < 0.05 level.
Verification of the functional roles of Glu175 in other xylanases
To verify the functional roles of Glu175 in other GH10 xylanases, the corresponding residue of XylE from P. scopiformis and XynE2 from Anoxybacillus sp. E2 were replaced with glutamic acid or other amino acid residues. The gene fragments of xylE and xynE2 (GenBank Accession Number: FJ860894 and GQ240232) without the signal peptide-coding sequence were synthesized by Biomed (Beijing, China). The gene fragments harboring mutations were obtained by overlap PCR with specific primers (Additional file 1: Table S6). For the XylE and its mutants, construction of recombinant vector, expression in P. pastoris, purification and characterization followed the same procedures as described above. The standard assay conditions for XylE and its mutants were 70 °C and pH 5.0. For XynE2 and its mutant, the gene fragment and vector pET-22b(+) (Novagen) were digested by NcoI and HindШ, ligated by T4 ligase, and transformed into E. coli BL21 (DE3) competent cells for expression. The standard assay conditions for XynE2 and its mutant were 65 °C and pH 7.5.
Enzymatic hydrolysis of corn stover
Corn stover was treated with 15% (w/w) ammonia solution for 24 h as previously reported . Pretreated corn stover was solved in 50 mM sodium citrate (pH 5.0) at the concentration of 2% (w/v), followed by the addition of 100 U of T. reesei cellulase or 30 U of XYL10C or XYL10C-ΔN and 70 U of T. reesei cellulase. The reactions were incubated at 50 °C with an agitation rate of 165 rpm for 24 h. The amounts of reducing sugars released were determined at 1, 2, 4, 8, 12, 18, and 24 h using the DNS method. The experiment was performed in triplicate.
sodium dodecyl sulfate-polyacrylamide gel electrophoresis
open reading frame
minimal dextrose medium
buffered glycerol-complex medium
buffered methanol-complex medium
albumin from bovine serum
- Endo H f :
Victor DG, Leape JP. Global climate agreement: after the talks. Nature. 2015;527:439–41.
Wahlstrom RM, Suurnakki A. Enzymatic hydrolysis of lignocellulosic polysaccharides in the presence of ionic liquids. Green Chem. 2015;17:694–714.
Goldbeck R, Damásio ARL, Gonçalves TA, Machado CB, Paixão DAA, Wolf LD. Development of hemicellulolytic enzyme mixtures for plant biomass deconstruction on target biotechnological applications. Appl Microbiol Biotechnol. 2014;98:8513–25.
Sims RE, Mabee W, Saddler JN, Taylor M. An overview of second generation biofuel technologies. Bioresour Technol. 2010;101:1570–80.
Braga CM, Delabona PS, Lima DJ, Paixã£O DA, Pradella JG, Farinas CS. Addition of feruloyl esterase and xylanase produced on-site improves sugarcane bagasse hydrolysis. Bioresour Technol. 2014;170:316–24.
Hu J, Arantes V, Saddler JN. The enhancement of enzymatic hydrolysis of lignocellulosic substrates by the addition of accessory enzymes such as xylanase: is it an additive or synergistic effect? Biotechnol Biofuels. 2011;4:36.
Harris PV, Welner D, McFarland KC, Re E, Navarro Poulsen JC, Brown K, Salbo R, Ding H, Vlasenko E, Merino S, Xu F, Cherry J, Larsen S, Lo Leggio L. Stimulation of lignocellulosic biomass hydrolysis by proteins of glycoside hydrolase family 61: structure and function of a large, enigmatic family. Biochemistry. 2010;49:3305–16.
Banerjee G, Car S, Scott-Craig JS, Borrusch MS, Bongers M, Walton JD. Synthetic multi-component enzyme mixtures for deconstruction of lignocellulosic biomass. Bioresour Technol. 2010;101:9097–105.
Collins T, Gerday C, Feller G. Xylanase, xylanase families and extremophilic xylanase. FEMS Microbiol Rev FEMS Microbiol Lett. 2005;29:3–23.
Taibi Z, Saoudi B, Boudelaa M, Trigui H, Belghith H, Gargouri A, Ladjama A. Purification and biochemical characterization of a highly thermostable xylanase from Actinomadura sp. Strain Cpt20 isolated from poultry compost. Appl Biochem Biotechnol. 2012;166:663–79.
Henrissat B, Callebaut I, Fabrega S, Lehn P, Mornon JP, Davies G. Conserved catalytic machinery and the prediction of a common fold for several families of glycosyl hydrolases. Proc Natl Acad Sci USA. 1996;93:7090–4.
Wakarchuk WW, Campbell RL, Sung WL, Davoodi J, Yaguchi M. Mutational and crystallographic analyses of the active site residues of the Bacillus circulans xylanase. Protein Sci. 1994;3:467–75.
Gary S, Stephen GW, Nham TN, Lawrence PM, Lothar Z, Gary DB. Sugar ring distortion in the glycosyl-enzyme intermediate of a family G/11 xylanase. Biochemistry. 1999;38:5346–54.
Kumar V, Marã-N-Navarro J, Shukla P. Thermostable microbial xylanases for pulp and paper industries: trends, applications and further perspectives. World J Microbiol Biotechnol. 2016;32:1–10.
Juturu V, Jin CW. Insight into microbial hemicellulases other than xylanases: a review. J Chem Technol Biotechnol. 2013;88:353–63.
Bhattacharya AS, Bhattacharya A, Pletschke BI. Synergism of fungal and bacterial cellulases and hemicellulases: a novel perspective for enhanced bio-ethanol production. Biotechnol Lett. 2015;37:1–13.
Zhang J, Siikaaho M, Puranen T, Ming T, Tenkanen M, Viikari L. Thermostable recombinant xylanases from Nonomuraea flexuosa and Thermoascus aurantiacus show distinct properties in the hydrolysis of xylans and pretreated wheat straw. Biotechnol Biofuels. 2011;4:12.
Hu J, Arantes V, Pribowo A, Saddler JN. The synergistic action of accessory enzymes enhances the hydrolytic potential of a “cellulase mixture” but is highly substrate specific. Biotechnol Biofuels. 2013;6:112.
Kim JH, Irwin D, Wilson DB. Purification and characterization of Thermobifida fusca xylanase 10B. Can J Microbiol. 2004;50:835–43.
Beaugrand J, Chambat G, Wong VW, Goubet F, Rémond C, Paës G, Benamrouche S, Debeire P, O’Donohue M, Chabbert B. Impact and efficiency of GH10 and GH11 thermostable endoxylanases on wheat bran and alkali-extractable arabinoxylans. Carbohydr Res. 2004;339:2529–40.
Cai H, Shi P, Bai Y, Huang H, Yuan T, Yang P, Luo H, Meng K, Yao B. A novel thermoacidophilic family 10 xylanase from Penicillium pinophilum C1. Process Biochem. 2011;46:2341–6.
Singh RK, Tiwari MK, Kim IW, Chen Z, Lee JK. Probing the role of sigma π interaction and energetics in the catalytic efficiency of endo-1,4-β-xylanase. Appl Environ Microbiol. 2012;78:8817–21.
Khajeh K, Ranjbar B, Naderimanesh H, Ebrahim HA, Nematgorgani M. Chemical modification of bacterial alpha-amylases: changes in tertiary structures and the effect of additional calcium. Biochim Biophys Acta. 2001;1548:229–37.
Moraïs S, Barak Y, Hadar Y, Wilson DB, Shoham Y, Lamed R, Bayer EA. Assembly of xylanases into designer cellulosomes promotes efficient hydrolysis of the xylan component of a natural recalcitrant cellulosic substrate. MBio. 2011;2:549–77.
Song L, Dumon C, Siguier B, André I, Eneyskaya E, Kulminskaya A, Bozonnet S, O’Donohue MJ. Impact of an N-terminal extension on the stability and activity of the GH11 xylanase from Thermobacillus xylanilyticus. J Biotechnol. 2014;174:64–72.
Hoque MA, Zhang Y, Chen L, Yang GY, Khatun MA, Chen HF, Hao L, Feng Y. Stepwise loop insertion strategy for active site remodeling to generate novel enzyme functions. ACS Chem Biol. 2017;12:1188–93.
Wang X, Ge H, Zhang D, Wu S, Zhang G. Oligomerization triggered by foldon: a simple method to enhance the catalytic efficiency of lichenase and xylanase. BMC Biotechnol. 2017;17:57.
Stemmer WPC. DNA shuffling by random fragmentation and reassembly: in vitro recombination for molecular evolution. Proc Natl Acad Sci USA. 1994;91:10747–51.
Leisola M, Turunen O. Protein engineering: opportunities and challenges. Appl Microbiol Biotechnol. 2007;75:1225–32.
Luo H, Li J, Yang J, Wang H, Yang Y, Huang H, Shi P, Yuan T, Fan Y, Yao B. A thermophilic and acid stable family-10 xylanase from the acidophilic fungus Bispora sp. MEY-1. Extremophiles. 2009;13:849–57.
Fedorova TV, Chulkin AM, Vavilova EA, Maisuradze IG, Trofimov AA, Zorov IN, Khotchenkov VP, Polyakov KM, Benevolensky SV, Koroleva OV. Purification, biochemical characterization, and structure of recombinant endo-1,4-β-xylanase XylE. Biochemistry (Mosc). 2012;77:1190–8.
Widmann M, Radloff R, Pleiss J. The thiamine diphosphate dependent enzyme engineering database: a tool for the systematic analysis of sequence and structure relations. BMC Biochem. 2010;11:9.
Ventorim RZ, Ta DO, Trevizano LM, Am DS, Guimarães VM. Impact of the removal of N-terminal non-structured amino acids on activity and stability of xylanases from Opinomyces sp. PC-2. Int J Biol Macromol. 2018;106:312–9.
Dias FM, Goyal A, Gilbert HJ, Prates JAM, Ferreira LM, Fontes CM. The N-terminal family 22 carbohydrate-binding module of xylanase 10b of Clostridium themocellum is not a thermostabilizing domain. FEMS Microbiol Lett. 2004;238:71–8.
Najmudin S, Pinheiro BA, Prates JA, Gilbert HJ, Romão MJ, Fontes CM. Putting an N-terminal end to the Clostridium thermocellum xylanase Xyn10B story: crystal structure of the CBM22-1-GH10 modules complexed with xylohexaose. J Struct Biol. 2010;172:353–62.
Li Z, Xue X, Zhao H, Yang P, Luo H, Zhao J, Huang H, Yao B. A C-terminal proline-rich sequence simultaneously broadens the optimal temperature and pH ranges and improves the catalytic efficiency of glycosyl hydrolase family 10 ruminal xylanases. Appl Environ Microbiol. 2014;80:3426–32.
Hantke K, Braun V. Covalent binding of lipid to protein. Diglyceride and amide-linked fatty acid at the N-terminal end of the murein-lipoprotein of the Escherichia coli outer membrane. Eur J Biochem. 1973;34:284–96.
Yin X, Li JF, Wang JQ, Tang CD, Wu MC. Enhanced thermostability of a mesophilic xylanase by N-terminal replacement designed by molecular dynamics simulation. J Sci Food Agric. 2013;93:3016–23.
Rudolph K, Bauer P, Schmid B, Mueller-Uri F, Kreis W. Truncation of N-terminal regions of Digitalis lanata progesterone 5β-reductase alters catalytic efficiency and substrate preference. Biochimie. 2014;101:31–8.
Wang X, Huang H, Xie X, Ma R, Bai Y, Zheng F, You S, Zhang B, Xie H, Yao B, Luo H. Improvement of the catalytic performance of a hyperthermostable GH10 xylanase from Talaromyces leycettanus JCM12802. Bioresour Technol. 2016;222:277–84.
Wang J, Tan Z, Wu M, Li J, Wu J. Improving the thermostability of a mesophilic family 10 xylanase, AuXyn10A, from Aspergillus usamii by in silico design. J Ind Microbiol Biotechnol. 2014;41:1217–25.
Wang Y, Feng S, Zhan T, Huang Z, Wu G, Liu Z. Improving catalytic efficiency of endo-β-1,4-xylanase from Geobacillus stearothermophilus by directed evolution and H179 saturation mutagenesis. J Biotechnol. 2013;168:341–7.
Ribeiro LF, Lucas RCD, Vitcosque GL, Ribeiro LF, Ward RJ, Rubio MV, Damásio AR, Squina FM, Gregory RC, Walton PH, Jorge JA, Prade RA, Buckeridge MS, Polizeli Mde L. A novel thermostable xylanase GH10 from Malbranchea pulchella expressed in Aspergillus nidulans with potential applications in biotechnology. Biotechnol Biofuels. 2014;7:115.
Wang X, Luo H, Yu W, Ma R, You S, Liu W, Hou L, Zheng F, Xie X, Yao B. A thermostable Gloeophyllum trabeum xylanase with potential for the brewing industry. Food Chem. 2016;199:516–23.
Lo LL, Kalogiannis S, Bhat MK, Pickersgill RW. High resolution structure and sequence of T. aurantiacus xylanase I: implications for the evolution of thermostability in family 10 xylanases and enzymes with (β)α-barrel architecture. Proteins. 1999;36:295–306.
Pell G, Taylor EJ, Gloster TM, Turkenburg JP, Fontes CM, Ferreira LM, Nagy T, Clark SJ, Davies GJ, Gilbert HJ. The mechanisms by which family 10 glycoside hydrolases bind decorated substrates. J Biol Chem. 2004;279:9597–605.
Chu Y, Tu T, Penttinen L, Xue X, Wang X, Yi Z, Gong L, Rouvinen J, Luo H, Hakulinen N. Insights into the roles of non-catalytic residues in the active site of a GH10 xylanase with activity on cellulose. J Biol Chem. 2017;292:19315–27.
White A, Tull D, Johns K, Withers SG, Rose DR. Crystallographic observation of a covalent catalytic intermediate in a β-glycosidase. Nat Struct Mol Biol. 1996;3:149–54.
Borchert TV, Abagyan R, Kishan KV, Zeelen JP, Wierenga RK. The crystal structure of an engineered monomeric triosephosphate isomerase, monoTIM: the correct modelling of an eight-residue loop. Structure. 1993;1:205–13.
Berlow RB, Igumenova TI, Loria JP. Value of a hydrogen bond in triosephosphate isomerase loop motion. Biochemistry. 2007;46:6001–10.
Kumar R, Wyman CE. Effect of xylanase supplementation of cellulase on digestion of corn stover solids prepared by leading pretreatment technologies. Bioresour Technol. 2009;100:4203–13.
Igarashi K, Uchihashi T, Koivula A, Wada M, Kimura S, Okamoto T, Penttilä M, Ando T, Samejima M. Traffic jams reduce hydrolytic efficiency of cellulase on cellulose surface. Science. 2011;333:1279–82.
Eriksson T, Börjesson J, Tjerneld F. Mechanism of surfactant effect in enzymatic hydrolysis of lignocellulose. Enzyme Microb Technol. 2002;31:353–64.
Qing Q, Yang B, Wyman CE. Xylooligomers are strong inhibitors of cellulose hydrolysis by enzymes. Bioresour Technol. 2010;101:9624–30.
Qing Q, Wyman CE. Supplementation with xylanase and β-xylosidase to reduce xylo-oligomer and xylan inhibition of enzymatic hydrolysis of cellulose and pretreated corn stover. Biotechnol Biofuels. 2011;4:1–12.
Henrissat B, Davies G. Structural and sequence-based classification of glycoside hydrolases. Curr Opin Struct Biol. 1997;7:637–44.
Kolenov, Aacute K, Biely P. Mode of action of endo-beta-1,4-xylanases of families 10 and 11 on acidic xylooligosaccharides. J Biotechnol. 2006;121:338–45.
Otwinowski Z, Minor W. Processing of X-ray diffraction data collected in oscillation mode. Methods Enzymol. 1997;276:307–26.
Mccoy AJ, Grossekunstleve RW, Adams PD, Winn MD, Storoni LC, Read RJ. Phaser crystallographic software. J Appl Crystallogr. 2007;40:658–74.
Miller GL. Use of dinitrosalicylic acid reagent for determination of reducing sugar. Anal Biochem. 1959;31:426–8.
Li X, Kim T, Nghiem N. Bioethanol production from corn stover using aqueous ammonia pretreatment and two-phase simultaneous saccharification and fermentation (TPSSF). Bioresour Technol. 2010;101:5910–6.
SY performed the major experiments including the site-directed mutagenesis and enzyme characterization, and drafted the manuscript; C-CC completed the crystal structure; TT, XW, and RM participated in the discussion and revised the manuscript; YW, HH, YB, and HC performed data analysis and guided the experiment. R-TG, HL, and BY were the corresponding authors. They designed the study, supervised the work and analyzed the experimental data. All authors reviewed the results of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The dataset supporting the conclusions of this article are included within the article (and its Additional file 1).
Consent for publication
All authors consent the manuscript for publication in Biotechnology for Biofuels.
Ethical approval and consent to participate
This research was supported by the National Natural Science Foundation of China (31472127), the Collaborative Innovation Project of Chinese Academy of Agricultural Sciences (CAAS-XTCX2016011-04-5), and the National Chicken Industry Technology System of China (CARS-41).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
X-Ray data collection and structure refinement statistics. Table S2. Temperature and pH optima of XYL10C-ΔN and its mutants with beechwood xylan as the substrate. Table S3. Specific activity of XYL10C-ΔN and its mutants with beechwood xylan as the substrate under standard conditions (80 °C and pH 4.0). Table S4. Temperature and pH optima of XylE and its mutants with beechwood xylan as the substrate. Table S5. Specific activity of XylE and its mutants with beechwood xylan as the substrate under standard conditions (70 °C and pH 5.0). Table S6. Primers used in this study. Figure S1. Multiple sequence alignments of XYL10C-ΔN and other GH10 xylanases using the ClustalW and ESPript. Strictly conserved residues are marked in red. Similar residues are shown in black bold characters and marked in yellow. The catalytic residues are indicated by red diamonds. Some important amino acids are indicated by red triangles. α-Helices and β-strands are represented by black coils and arrows, respectively. Figure S2. SDS-PAGE analysis of the purified recombinant XYL10C, XYL10C-ΔN, XylE, XynE2, and their mutants. Lanes: M, the standard protein molecular weight markers; 1, 3, 10, and 45, the crude enzyme of wild-type XYL10C, XYL10C-ΔN, XylE and XynE2; 2, 4, and 11, the deglycosylated XYL10C, XYL10C-ΔN, and XylE; 5-9 and 13-26, the saturated mutants of XYL10C-ΔN; 12 and 27-44, the saturated mutants of XylE; and 46, the mutant XynE2-Q85E. Figure S3. Comparison of thermostability of XYL10C and XYL10C-ΔN. A, the kinetic stability assayed at 80 °C and 90 °C; B, the thermodynamic stability (Tm values). Figure S4. pH stability of XYL10C and XYL10C-ΔN. Figure S5. HPLC analysis of the hydrolysis products by XYL10C (A) and XYL10C-ΔN (B). X1, xylose; X2, xylobiose; X3, xylotriose; X4, xylotetraose; X5, xylopentaose; X6, xylohexaose; beechwood, beechwood xylan. Figure S6. The X-ray diffraction pattern of the XYL10C-ΔN crystal. The resolution edges are shown by different concentric circles. Figure S7. The dimeric structure (Chains A and B) of the XYL10C-ΔN. A, The hydrogen bonds between monomers. The residues involved in forming hydrogen bonds are indicated, and the red dotted line represents hydrogen bond. B, The van der Waals force network between monomers. The residues in blue are from chain A, and those red are from chain B. The residues involved are represented by yellow dashed lines. Figure S8. The XYL10C-ΔN docked with xyloheptaose. The enzyme and substrate interactions are indicated by yellow dotted lines. Figure S9. The active-site architecture of mutant XylE-Q116E.
About this article
Cite this article
You, S., Chen, CC., Tu, T. et al. Insight into the functional roles of Glu175 in the hyperthermostable xylanase XYL10C-ΔN through structural analysis and site-saturation mutagenesis. Biotechnol Biofuels 11, 159 (2018). https://doi.org/10.1186/s13068-018-1150-8