Structural Analysis of a Glycoside Hydrolase Family 11 Xylanase from Neocallimastix patriciarum

Background: Thermophilic xylanases are valuable in many industrial applications. Results: The structures of a xylanase XynCDBFV and its complex with xylooligosaccharides were determined, and its N-terminal region (NTR) contributes to thermostability. Conclusion: NTR may stabilize the overall protein folding of XynCDBFV. Significance: The structural and functional investigation of unprecedented NTR of XynCDBFV provides a new insight into the molecular basis of thermophilic xylanases. The catalytic domain of XynCDBFV, a glycoside hydrolase family 11 (GH11) xylanase from ruminal fungus Neocallimastix patriciarum previously engineered to exhibit higher specific activity and broader pH adaptability, holds great potential in commercial applications. Here, the crystal structures of XynCDBFV and its complex with substrate were determined to 1.27–1.43 Å resolution. These structures revealed a typical GH11 β-jelly-roll fold and detailed interaction networks between the enzyme and ligands. Notably, an extended N-terminal region (NTR) consisting of 11 amino acids was identified in the XynCDBFV structure, which is found unique among GH11 xylanases. The NTR is attached to the catalytic core by hydrogen bonds and stacking forces along with a disulfide bond between Cys-4 and Cys-172. Interestingly, the NTR deletion mutant retained 61.5% and 19.5% enzymatic activity at 55 °C and 75 °C, respectively, compared with the wild-type enzyme, whereas the C4A/C172A mutant showed 86.8% and 23.3% activity. These results suggest that NTR plays a role in XynCDBFV thermostability, and the Cys-4/Cys-172 disulfide bond is critical to the NTR-mediated interactions. Furthermore, we also demonstrated that Pichia pastoris produces XynCDBFV with higher catalytic activity at higher temperature than Escherichia coli, in which incorrect NTR folding and inefficient disulfide bond formation might have occurred. In conclusion, these structural and functional analyses of the industrially favored XynCDBFV provide a molecular basis of NTR contribution to its thermostability.

The catalytic domain of XynCDBFV, a glycoside hydrolase family 11 (GH11) xylanase from ruminal fungus Neocallimastix patriciarum previously engineered to exhibit higher specific activity and broader pH adaptability, holds great potential in commercial applications. Here, the crystal structures of XynCD-BFV and its complex with substrate were determined to 1.27-1.43 Å resolution. These structures revealed a typical GH11 ␤-jelly-roll fold and detailed interaction networks between the enzyme and ligands. Notably, an extended N-terminal region (NTR) consisting of 11 amino acids was identified in the XynCDBFV structure, which is found unique among GH11 xylanases. The NTR is attached to the catalytic core by hydrogen bonds and stacking forces along with a disulfide bond between Cys-4 and Cys-172. Interestingly, the NTR deletion mutant retained 61.5% and 19.5% enzymatic activity at 55°C and 75°C, respectively, compared with the wild-type enzyme, whereas the C4A/C172A mutant showed 86.8% and 23.3% activity. These results suggest that NTR plays a role in XynCDBFV thermostability, and the Cys-4/Cys-172 disulfide bond is critical to the NTR-mediated interactions. Furthermore, we also demonstrated that Pichia pastoris produces XynCDBFV with higher catalytic activity at higher temperature than Escherichia coli, in which incorrect NTR folding and inefficient disulfide bond formation might have occurred. In conclusion, these structural and functional analyses of the industrially favored XynCDBFV provide a molecular basis of NTR contribution to its thermostability.
Xylans are the major hemicellulose components in the plant cell wall and account for nearly one-third of all renewable organic carbon source on earth (1,2). Xylans are heteropolysaccharides composed of ␤-1,4-glycosidic bond-linked xylose units as a backbone chain that is usually decorated by different side groups such as methyl group, acetyl group, and other sugar molecules (1,3). Due to its structural complexity, a set of enzymes is required for complete xylan decomposition, including endo-1,4-␤-xylanase, ␤-xylosidase, acetylxylan esterase, arabinase, and ␣-glucuronidase (3,4). Among them, the glycoside hydrolase endo-1,4-xylanase (xylanase, EC 3.2.1.8) is the key enzyme that catalyzes random hydrolysis of the xylan backbone to small fragments by cleaving the ␤-1,4-glycosidic bonds (1). Xylanases have been widely applied in many industries such as feed manufacture, paper and pulp processing, and food industry (5)(6)(7). These biotechnological treatments usually involve harsh conditions and demand enzymes with good thermostability, broad pH adaptability, and high specific activity. Therefore, numerous research projects have been carried out in search of novel xylanases with favorable properties and to improve performance of available enzymes via direct evolution or rational design approaches.
Previously, a GH11 xylanase from an anaerobic ruminal fungus, Neocallimastix patriciarum, was isolated and characterized (14,15). The catalytic domain of the enzyme (Xyn-CD) exerts optimal activity at pH 6.0, and an alkalophilic mutant XynCDBFV with seven amino acid substitutions was later created via direct evolution using error-prone PCR (16). The recombinant XynCDBFV protein is among the highest active xylanases and possesses an optimal temperature of 65°C, broad pH adaptability, and remarkable tolerance at pH 10.0 when expressed in E. coli. Therefore, the XynCDBFV is an attractive candidate to be developed as an industrial product. In the present study the recombinant XynCDBFV was expressed in the industrial strain P. pastoris and crystallized. The overall protein fold and ligand complex structure are analyzed in detail. Based on these data, potential factors contributing to the enzyme thermostability are proposed.

EXPERIMENTAL PROCEDURES
Gene Cloning and Mutagenesis-The synthesized gene encoding XynCDBFV, an engineered mutant of Xyn-CD from N. patriciarum (GenBank TM accession number AF123252) (16), was amplified by using PCR with a forward primer of 5Ј-CCCGAATTCCAAAGTTTCTGTAGTTCAGCTTCT-3Ј and a reverse primer of 5Ј-CCCGCGGCCGCTTAATCACCAA-TGTAAACCTTTGCGTA-3Ј. The gene was then cloned into the vector pPICZ␣A for the P. pastoris system by EcoRI and NotI to yield pPICZ␣A/xynCDBFV. The substituted mutants including E109A, C4A, C172A, and C4A/C172A were prepared by using the QuikChange site-directed mutagenesis kit (Agilent) with pPICZ␣A/xynCDBFV as the template. The genes encoding the deleted mutants of ⌬6 (deletion of Gln-1-Ser-6) and ⌬11 (deletion of Gln-1-Gly-11) were generated by PCR with full-length xynCDBFV gene as the template. These truncated genes were then cloned into the vector pPICZ␣A by using EcoRI and NotI to yield pPICZ␣A/xynCDBFV-⌬6 and pPICZ␣A/xynCDBFV-⌬11. The sequences of the mutated primers are listed in the supplemental Table S1.
Alternatively, the xynCDBFV gene was amplified by using PCR and cloned into the vector pET32 Xa/LIC for E. coli expression system. This vector has designed a His tag before the N terminus of targeted gene for purification purpose. The specific primers used here were 5Ј-GGTATTGAGGGTCGCCA-AAGTTTCTGTAGTTCAGCT-3Ј (forward) and 5Ј-AGAG-GAGAGTTAGAGCCTTAATCACCAATGTAAACCTT-TGC-3Ј (reverse).
Protein Expression and Purification-These above plasmids, except pET32 Xa/LIC-xynCDBFV, were linearized by PmeI and then individually transformed into X33 strain of P. pastoris by electroporation. The transformants were selected on the YPD (yeast extract peptone dextrose) plates containing 100 g/ml zeocin (Invitrogen). The selected clones were inoculated and amplified in 50 ml of buffered glycerol-complex medium (BMGY; 1% yeast extract, 2% peptone, 100 mM potassium phosphate, pH 6.0, 1.34% yeast nitrogen base (YNB) with ammonium sulfate without amino acids, 4ϫ 10 -5% biotin, and 1% glycerol) at 30°C for 1 day. Then the culture medium of cultured cells was replaced by 20 ml of buffered methanol-complex medium (BMMY; 1% yeast extract, 2% peptone, 100 mM potassium phosphate, pH 6.0, 1.34% yeast nitrogen base (YNB) with ammonium sulfate without amino acids, 4ϫ 10 -5% biotin, and 0.5% methanol) to induce protein expression. For protein purification, the supernatants were collected and dialyzed twice against the buffer containing 25 mM Tris, pH 7.5. In addition, the proteins were treated by endoglycosidase H (New England Biolabs) during the dialysis procedure. The proteins were then purified by FPLC system using diethylaminoethyl (DEAE) column (GE Healthcare) and eluted using a linear gradient of 0 -250 mM NaCl in the buffer containing 25 mM Tris, pH 7.5. The purified proteins were finally concentrated to 10 mg/ml in 25 mM Tris, pH 7.5, 150 mM NaCl by using Amicon Ultra-15 Centrifugal Filter Units (Millipore), and the purity (Ͼ95%) was checked by SDS-PAGE.
On the other hand the pET32 Xa/LIC-xynCDBFV plasmid was transformed into BL21 (DE3) strain of E. coli. Then the transformed cells were propagated and induced by adding isopropyl 1-thio-␤-D-galactopyranoside for protein expression. The protein was purified by FPLC system using a nickel nitriloacetic acid column. The His-tagged protein was eluted using a gradient of 0 -250 mM imidazole in the buffer containing 25 mM Tris, pH 7.5, 150 mM NaCl. The eluted protein was then dialyzed against the buffer containing 25 mM Tris, pH 7.5, 150 mM NaCl, with Factor Xa enzyme to remove the His tag. The mixture was loaded onto another nickel nitriloacetic acid column, and the untagged protein was eluted in the buffer without imidazole for purification of untagged protein. The eluted protein was then dialyzed against the buffer containing 25 mM Tris, pH 7.5, for the DEAE column. Next, the purified protein was eluted using a gradient of 0 -250 mM NaCl in the 25 mM Tris, pH 7.5, for the DEAE column. The purified proteins were finally concentrated to 10 mg/ml in 25 mM Tris, pH 7.5, 150 mM NaCl, and the purity was checked by SDS-PAGE.
Crystallization and Data Collection-The XynCDBFV and E109A proteins were crystallized by using sitting-drop vapor diffusion method and Crystal Screen kit (Hampton Research). The XynCDBFV protein crystals were obtained from the reservoir solution containing 0.1 M Tris, pH 8.5, and 2 M ammonium sulfate at room temperature for 1 day. The E109A protein crystals were obtained from the reservoir solution containing 0.1 M sodium cacodylate, pH 6.5, 0.2 M ammonium sulfate, 26% PEG8000, and 5% glycerol at room temperature for 2 days. The hexagonal E109A-xylotriose-bound crystal was prepared by soaking with 10 mM xylotriose (Megazyme) in the reservoir solution for 1 h. The cryoprotectants for XynCDBFV and E109A crystals contained 0.1 M Tris, pH 8.5, 2 M ammonium sulfate and 10% glycerol and 0.12 M sodium cacodylate, pH 6.5, 0.24 M ammonium sulfate, 31% PEG8000, and 5% glycerol, respectively. All of the x-ray diffraction data were collected at beam line BL13B1 of National Synchrotron Radiation Research Center in Hsinchu, Taiwan. The diffraction images were processed by using HKL2000 (17).
Structure Determination and Refinement-The crystal structure of XynCDBFV was solved by the molecular replacement method with the program Phaser (18) from the CCP4 suite (19) using the hypothetical XynCDBFV model generated from the structure of Bacillus subtilis B230 Xyn11X (PDB code 1IGO; 47% sequence identity with XynCDBFV) by the SWISS-MODEL website (20, 21) as a search model. Subsequent model building and structural refinement were carried out by using the programs COOT (22) and REFMAC5 (23), respectively. The complex structure of E109A-xylotriose was determined by the molecular replacement method with Phaser using refined XynCDBFV structure as a search model. The structural refinements were finished by the programs COOT (22) and . Some data collection and statistics are summarized in Table 1. All of the structural diagrams were drawn by using PyMOL (24).
Xylanase Activity Assay-The detection of xylanase activity was based on the determination of reducing sugar by using dinitrosalicylic acid method (25,26). In general, the protein solution in 50 mM sodium acetate buffer, pH 5.3, with a proper concentration was mixed with 1% xylan substrate (beechwood, Sigma) in the proportion of 1 to 9 and then incubated at 55°C for 10 min. The reaction was stopped by adding 1% dinitrosalicylic acid solution and incubated at 100°C in boiling water for 10 min. The absorbance of A 540 nm was measured for calculation of the enzyme activity. The standard curve for calibrating the enzyme activity was determined by 0 -0.6 mg/ml xylose solutions. One unit of activity is defined as the amount of enzyme that releases 1 mol of reducing sugar equivalent to xylose per minute.

RESULTS
Overall Structure-The crystal structure of XynCDBFV was solved to a resolution of 1.27 Å by molecular replacement using the Xyn11X from B. subtilis B230 as a search model. The data collection and refinement statistics of the wild type, E109A mutant (Ͻ0.2% activity, data not shown), and its complex with xylotriose (XTI, E109A-XTI) structures are listed in Table 1. The root mean square deviations (r.m.s.d.) of C␣ atoms between XynCDBFV and the E109A mutant is 0.374 Å, indicating that the mutation did not cause a significant alteration in the protein structure. To investigate whether the ligand binding could cause significant conformational change, the ligand-free (XynCDBFV and E109A mutant) and ligand-bound (E109A-XTI) structures were superimposed. The r.m.s.d. between these structures ranged from 0.192 to 0.436 Å for all C␣ atoms, suggesting that the ligand binding did not result in significant conformational change.
As shown in Fig. 1A, the protein structure of XynCDBFV displays a ␤-jelly-roll fold, which is typical of GH11 family enzymes. It comprises two anti-parallel ␤-sheets, linked by interconnecting loops and ␣-helices. The topology of GH11 enzymes is also described as the shape of a right hand (27). The two ␤-sheets mimic the palm and fingers, whereas the loop ␤12-␤13 is like a thumb. As shown in Fig. 1B, there are two disulfide bridges in the XynCDBFV structure, including DS1 (Cys-4 -Cys-172), which connects the N-terminal ␣1 helix to the strand ␤14, and DS2 (Cys-50 -Cys-60), which joins the strand ␤5 to the strand ␤6. The tunnel-like active site cleft (Fig.  1C) is formed by the curved inner ␤-sheet consisting of ␤2, ␤3, ␤6, ␤15, ␤8, ␤9, ␤10, ␤12, and ␤13 strands (Fig. 1, A and C). The catalytic residues Glu-109 and Glu-202 predicted from sequence alignment of GH11 xylanases were found embedded in this region and located on the strand ␤10 and ␤15, respectively (Fig. 1, A and D). Substrate Binding Site-The electron density maps in the tunnel-like cleft of the E109A-XTI structure clearly indicated the presence of three ␤-1,4-linked xylosyl moieties in the ϩ1 to ϩ3 subsites and two in the Ϫ2 to Ϫ3 subsites ( Fig. 2A). The sugar molecules were modeled into the wild-type XynCDBFV to display the detailed interactions between the substrates and the surrounding residues (Fig. 2B). Trp-32 and Trp-125 provide stacking forces to sugar units in the subsites Ϫ2 and ϩ3, respectively. Trp-32 further forms a hydrogen bond to the hydroxyl group of the Ϫ3 sugar. Glu-30, Arg-61, and Tyr-100 form hydrogen bonds to the Ϫ2 sugar unit. Asp-57, Tyr-94, Tyr-111, Arg-148, Gln-161, and Glu-202 (one of catalytic residues) form hydrogen bonds to sugar moiety in the ϩ1 subsite. Arg-148 also forms a hydrogen bond to the ϩ2 sugar unit. The Ϫ1 sugar moiety was not observed from our crystal structure, but the subsite can be predicted from the superimposition of the XynCDBFV-XTI model and the complex structure of XynII from Trichoderma reesei (PDB code 4HK9) (28). The superimposition indicates that several residues might participate in the Ϫ1 subsite formation including Asp-57, Arg-148, Pro-151, and Gln-161 along with two catalytic residues, Glu-109 and Glu-202 (Fig. 2C). Among these ligand-interacting residues, Trp-32, Tyr-100, Tyr-111, Gln-161, and Pro-151 are strictly conserved in GH11 xylanase, indicating their vital roles in catalytic reactions.
The Unique N-terminal Region of XynCDBFV-From protein sequence alignment and structure superimposition, we found that XynCDBFV carries an extended N-terminal region (NTR) that is unique in the GH11 family (Fig. 3, A and B). The NTR is composed of 11 amino acids and spans the convex side of the palm ␤-sheet, whereas other GH11 enzymes start from the strand ␤1 or ␤2 (Fig. 3B). There are several interactions formed between the NTR and the nearby ␤-strands including ␤4, ␤5, ␤14, ␤16, and ␤17 (Fig. 3C). Two stacking interactions are formed between NTR-Phe-3 and ␤17-Tyr-218 and between NTR-His-9 and ␤4-Tyr-43 (Fig. 3C). NTR-His-9 is also at a hydrogen bonding distance to ␤5-Ser-47 and ␤16-Asp-215. Notably, NTR-Cys-4 forms the DS1 disulfide bond with ␤14-Cys-172. Through these interactions, the NTR is stably attached to the catalytic core of XynCDBFV instead of hanging freely as a flexible segment.
Functional Analysis of the NTR-To investigate the functional significance of NTR, we constructed the NTR-truncated mutants, including ⌬6 (Gln-1-Ser-6 deletion) and ⌬11 (Gln-1- Gly-11 deletion) (Fig. 4A), and analyzed their performance. As shown in Fig. 4B, the enzymatic activities of mutants ⌬6 and ⌬11 were reduced to 80.9 and 61.5% at 55°C, respectively, compared with the wild-type protein. Remarkably, the variations between wild-type and mutant enzymes were even larger at higher temperatures. As the wild-type XynCDBFV showed enhanced activity at 65°C (140%) and 75°C (151.8%), the relative activities of ⌬6 and ⌬11 were reduced to 33.7 and 20.7% at 65°C and 21 and 19.6% at 75°C (Fig. 4B). These results clearly demonstrated that the NTR plays a role in the catalytic activity of XynCDBFV, and its presence is especially important for the enzyme to exert thermophilic functions. As mentioned, there is a disulfide bond DS1 formed between the NTR and strand ␤14. To further investigate the role of the DS1 linkage, mutants C4A, C172A, and C4A/C172A were constructed. Similar to the results of the NTR-truncated mutants, the DS1-removed mutants showed lower activities comparing to the wild-type enzyme at 55°C, and their activities further declined when the FIGURE 2. The electron density maps of ligands and the detailed interaction networks in the active site of XynCDBFV. A, the refined oligosaccharide models of E109A-XTI are superimposed on the 2F o Ϫ F c electron density maps, which are contoured at 1 and 2.5 levels and shown in cyan and red, respectively. B, a stereo view of detailed interactions between the active-site residues of XynCDBFV and ligands is shown with side chains of amino acids as thin green sticks and sugar moieties as gray thick sticks. The dashed lines denote potential hydrogen bonds between the enzyme and ligands. C, a stereo view of detailed interactions of active site residues and ligands in the glycone subsites. XynCDBFV-XTI model and T. reesei xylanase complex structure (PDB code 4HK9) are superimposed and the xylobiose (XBI) from the former and XTI from the latter are shown as stick models and colored in gray and salmon, respectively.  (47). The secondary structural elements (␣-helices and ␤-strands) of XynCD-BFV are shown above the sequences. Strictly conserved residues are highlighted by a red background, and conservatively residues are boxed. The figure was produced by ESPript (48). B, the structure of XynCDBFV is superimposed with other crystal structures of GH11 xylanases. The colors used here are red for XynCDBFV, cyan for 1IGO, magenta for 2DCJ, yellow for 1H4H, light pink for 1F5J, green for 2C1F, blue for 1M4W, purple for 1PVX, orange for 1H1A, dark green for 2NQY, pink for 1HIX, and light yellow for 3ZSE. C, the side chains of residues participating in the detailed interactions of NTR and catalytic core of XynCDBFV are shown as thin sticks and colored in blue and green, respectively. The disulfide bond, labeled as DS1, is shown as thick sticks in yellow. The dashed lines denote potential hydrogen bonds.
tested temperature was increased to 65°C (39.3-44.8%) and 75°C (20.8 -23.3%) (Fig. 4C). Protein expression levels of all mutants were similar to that of wild type (data not shown). These results suggested that the NTR is important for the catalytic activity and is an essential element for thermophilicity of XynCDBFV, and its function is highly dependent on the presence of disulfide bond DS1.
In the beginning we tried to use E. coli-expressed protein (E-XynCDBFV) for structural study because the E. coli expression system is easier to handle than P. pastoris expression system. The E-XynCDBFV was easily expressed and purified but failed to be crystallized despite the protein being catalytically active (Table 2). When the purified proteins were subjected to SDS-PAGE analysis, the E-XynCDBFV formed a significant amount of dimer and trimer in the absence of reducing agent, suggesting the presence of one or more free Cys residues on the protein surface (Fig. 5). Reducing agent was supplemented during E-XynCDBFV crystallization to prevent formation of nonspecific disulfide bonds, but no crystal was obtained. However, the crystal structure solved in this study indicated that the four Cys residues in XynCDBFV form two pairs of disulfide bonds, and there is no free thiol group exposed to the bulk solvent. Indeed, the recombinant P-XynCDBFV, which was successfully crystallized, showed much less tendency of dimer formation (Fig. 5).
Intriguingly, the specific activity of E-XynCDBFV was lower than that of P-XynCDBFV by 20% at 55°C and only slightly increased at higher experimental temperature (up to 18% increase), whereas the P-XynCDBFV showed a 55.7 and 63.5% increase in specific activity at 65 and 75°C, respectively ( Table  2). These observations were reminiscent of the results of analyzing NTR-truncated and DS1-removed mutants. Therefore, we suspect that the XynCDBFV was misfolded in E. coli to a certain extent so that the Cys residues were exposed. The disruption of the DS2 linkage is rather unlikely because it is embedded in the protein interior and its exposure should cause deleterious effects to the enzyme. It would be more likely that the NTR was incorrectly folded and disrupted DS1 formation. The observation that DS1-removed mutant C4A/C172A of FIGURE 4. Enzymatic activity analysis of wild-type and mutant XynCDBFV. A, schematic and schematic diagrams of the NTR-truncated mutants. The disulfide bond (Cys-4 -Cys-172) is shown as yellow sticks. B and C, the culture supernatant of P. pastoris expressing wild-type XynCDBFV was tested for activity and compared with those of NTR-truncated mutants (B) and DS1-removed mutants (C). The enzymatic activity was determined at the indicated temperatures, and the relative activity of each sample was presented as a percentage of the value of wild-type XynCDBFV at 55°C.  P-XynCDBFV hardly formed a dimer may support our hypothesis (Fig. 5). In conclusion, the failure of obtaining the E-XynCDBFV crystals and the lower activity at higher temperatures might be attributed to the inefficient formation of the DS1 linkage in the E. coli system, and P. pastoris should be a more ideal system to express XynCDBFV protein with better performance.

DISCUSSION
Thermostability is one of the most valued properties for industrial enzymes; thus, investigating the strategies utilized by thermophilic enzymes is an intense area of research. XynCD-BFV exhibits remarkably high specific activity, which is appreciated for a thermophilic enzyme and is worthy of further investigations. Here the crystal structure of XynCDBFV was solved at high resolution, revealing an unprecedented NTR element that is functionally related to enzymatic activity and thermophilicity. The underlying cause of NTR-mediated thermophilicity is thus an intriguing issue to address. Previous studies have found correlations between the thermostability/thermophilicity of GH11 enzymes and the increased intramolecular interactions including hydrogen bonding, salt bridges, disulfide bonds, and hydrophobic forces in them (10 -13, 29, 30). Although various kinds of interactions were proposed to contribute to enzyme thermostability, it seems that making the protein fold into a more compact state is a fundamental principle. Notably, heat-mediated unfolding of a GH11 enzyme is proposed to initiate from the N-terminal region (31). Therefore, introducing an N-terminal disulfide bond to stabilize the first ␤-strand has been reported as a useful strategy to enhance enzyme thermostability (32)(33)(34)(35). In XynCDBFV, the NTR adheres to several ␤-strands and spans the back of the outer ␤-sheet. The Cys-4 -Cys-172 DS1 bond on the most N-terminal helix ␣1 serves as a molecular rivet to fix the NTR to strand ␤14. Through these interactions, the NTR clamps the catalytic core into a more restrained architecture, which may prevent its destabilization at higher temperatures. This is the first report to reveal an extended NTR feature and its function in the GH11 family, and it shall be interesting to examine whether this element could also provide benefits when incorporated into other GH11 enzymes.
Among the 26 crystal structures of GH11 xylanases, only 6 possess one single disulfide bond (8). Three of them are in an equivalent position that connects the cord and strand ␤12, including the xylanases from Aspergillus kawachii (36), Aspergillus niger (37), and Scytalidium acidophilum (38). For xylanases from Paecilomyces varioti Bainier (30) and Thermomyces lanuginosus (12), the disulfide bond is located between helix ␣4 and strand ␤11. Another xylanase from an uncultured bacterium forms the disulfide bond between the N-terminal strands ␤1 and ␤4 (39). The XynCDBFV structure in this study, revealing two disulfide bridges DS1 (Cys-4 -Cys-172) and DS2 (Cys-50 -Cys-60), is the first GH11 structure that includes multiple disulfide bonds in its catalytic domain. A recent study inferred that XynCDBFV has only one disulfide bond formed between Cys-4 and Cys-172 based on the thiol titration analysis (40). The discrepancy might be a result of different experimental approaches. In the present study the clear electron density map has demonstrated the linkage between Cys-50 and Cys-60 and thus validates the existence of DS2 bond in the XynCDBFV structure. Notably, removing the DS2 bond from XynCDBFV did not significantly affect its thermostability (40). These findings suggest different roles played by DS1 and DS2 in XynCDBFV.
The N-terminal deletion experiments showed the ⌬11 mutant has lower activity than the ⌬6 mutant at 55 and 65°C, raising the possibility that one or more residues between Ala-7 and Gln-12 in NTR may also be involved in interactions for XynCDBFV thermostability aside from the DS1 disulfide bond. The most possible candidate is His-9, which provides hydrogen bonding to Ser-47 and Asp-216 as well as stacking force to Tyr-43. The forces provided by His-9 are weaker than the covalent disulfide linkages and are more readily dissociated at higher temperatures, which may explain the similar enzymatic activities exhibited by ⌬6 and ⌬11 at 75°C (Fig. 4B). This assumption needs to be further validated in the future.
Our results clearly demonstrated that the P. pastoris system is superior to the E. coli system in obtaining XynCDBFV proteins with better performance. Since it was developed, P. pastoris has become a very useful heterologous expression system for both industrial and academic purposes that numerous soluble enzymes with higher activities and better thermal profiles were produced with this system (41)(42)(43). The P. pastoris is particularly attractive in making eukaryotic proteins because yeast cells provide correct post-translational modifications such as phosphorylation, glycosylation, and disulfide bond formation (44,45). XynCDBFV possesses two putative N-glycosylation sites, and P-XynCDBFV showed a glycosylated pattern (data not shown). But the possibility of glycosylation-mediated thermostability is excluded by the fact that the enzymatic performance was not influenced by deglycosylation treatments (the parental proteins in Fig. 4 versus the deglycosylated protein in Table 2). Instead, the correct folding of NTR and the DS1 linkage formation are likely responsible for the different thermal profiles of E-XynCDBFV and P-XynCDBFV.
Previous and present studies show that utilizing the E. coli expression system yields E-XynCDBFV protein that exhibits its highest activity at no more than 65°C and is less active than the P-XynCDBFV at every temperatures tested (16). Similar phenomenon was observed in our previous study that analyzed the structure of rhodostomin, which contains six pairs of disulfide bonds. NMR examinations indicated that the E. coli-produced rhodostomin was misfolded despite being capable of inhibiting platelet aggregation (46). In contrast, P. pastoris produced correctly folded rhodostomin with all of six intramolecular disulfide pairings formed, which also exhibited higher activity than that from E. coli. Consequently, choosing a suitable production host is critical for XynCDBFV to exhibit its optimal performance, and the utilization of an industrial strain shall also benefit its further commercial applications. Nevertheless, the enhanced thermostability of XynCDBFV requires the presence of NTR in the first place.