Sequence and structure analysis of HWxyn11
The sequence of HWxyn11 gene from Hortaea werneckii EXF-12619 has one open reading frame of 993 bp encoding 331 amino acids with a putative signal peptide of 20 amino acids (Fig. 1). The signal peptide was removed in cloning and expression subsequently. Based on the BLAST analysis, xylanase HWxyn11 shows the highest similarity to GH 11 family and exhibits the highest amino acid identity of 58.08% with xylanase rXyn11A from Thermobifida fusca (PDB codes: 3zse.1) (Zhao et al. 2015). The sequence alignment and conserved domain prediction show that HWxyn11 possesses only one catalytic domain of GH 11 and two glutamate residues are important in the catalytic reaction as same as other GH 11 xylanases.
The structure of HWxyn11 is a classical right-handed palm consisting of α-helixs and ß-sheets, involving two glutamates acting as catalytic residues as other xylanases belonging to GH 11 (Fig. 2). The hydrophobic faces are packed against each other to form the hydrophobic core of protein. Exceptionally, HWxyn11 contains a propeptide in front of protein catalytic domain at N-terminus and an unfolded random curly with a length of 109 amino acids at C-terminus.
Expression and purification of HWxyn11
After inducing with 0.5% (v/v) methanol at 29 ℃ for 108 h, the mature HWxyn11 without signal peptide was successfully expressed in Pichia pastoris X-33. The protein was purified to apparent homogeneity by Ni-NTA affinity chromatography. The specific activity of HWxyn11 was 99.94 U/mg. Results of sodium dodecyl sulfate-polyacrylamide gel electrophoresis indicated that recombinant HWxyn11 was approximately 60 kDa and significantly larger than the predicted value of 33.78 kDa. However, the molecular weight of HWxyn11 did not change after deglycosylation treatment with endoglycosidase F1 to remove N-glycosylation (Fig. 3). The results indicated that increased molecular weight of HWxyn11 was due to extensive O-glycosylation modifications.
Biochemical characteristics of HWxyn11
As shown in Fig. 4a, HWxyn11 revealed its peak activity at 45°C. HWxyn11 is a mesophilic xylanase, maintaining more than 80% relative enzyme activity with a wide range from 35 ℃ to 60 ℃. Results of thermostability tests indicated that HWxyn11 showed high stability at lower temperature, preserving nearly 88% of its initial activity at 45 ℃ after 60 min. However, the enzyme activity was only 19.34% after incubating at 55 ℃ for 60 min (Fig. 4b). The half-life of HWxyn11 at 55 ℃ is 22.95 min.
Consistent with most reported xylanases of GH 11, HWxyn11 exhibited highest activity under weakly acidic conditions at pH 5.0(Luong et al. 2023; Kim et al. 2023; Zhang et al. 2023). HWxyn11 maintained more than 75% of its maximal enzyme activity in a range of pH 4.0–7.0 (Fig. 4c). Furthermore, HWxyn11 showed exceptional stability between pH 3.0–11.0, retaining over 87% of its residual enzyme activity after 24 h incubation at 4 ℃ (Fig. 4d).
Figure 4e shows the effects of metal ions on the enzyme activity of HWxyn11. Fe2+, Al3+ had inhibitory effect on recombinant HWxyn11. However, Na+ improved the enzyme activity to 1.16 times at a concentration of 10 mM. Other metal ions such as K+, Ca2+, Mg2+, Fe3+, Cu2+, Zn2+, Co+, Ni2+, Mn2+ and Li+ all showed little influence on the enzyme activity of HWxyn11, and residual enzyme activity is all above 78% after incubation at 4 ℃ for 24 h.
As shown in Fig. 4f, HWxyn11 is a halophilic xylanase. The stability of HWxyn11 was barely affected in 0–5.0 M concentration of NaCl after incubating at 4 ℃ for 24 h. The enzyme activity retained about 83% of when hydrolyzing beechwood xylan at a final concentration of 5.0 M NaCl.
Construction of mutants based on predicted 3-dimensional structure and molecular dynamic simulation
As shown in Fig. 5, the RMSF values and B-factor for N-terminal and C-terminal regions of HWxyn11 are significantly higher than those in other areas. Generally, the higher the RMSF values and B-factor, the more flexible the amino acid residues are(Sun et al. 2019). Both regions mentioned above exhibit greater structural flexibility and lower thermostability.
By sequence comparison and structural domain analysis, there is a relatively low similarity of N-terminal and C-terminal sequences between HWxyn11 and thermophilic xylanase rXyn11A. Moreover, rXyn11A is composed of a GH 11 xylanase structural domain and a CBM 2. However, HWxyn11 contains a GH 11 xylanase structural domain and a long flexible loop region at C-terminus instead of CBM (Fig. 6). It is inferred that significant differences in N-terminal and C-terminal sequences and structural domain composition of HWxyn11 may be important factors affecting its thermostability.
The structure of wild-type HWxyn11 was superimposed with the structure of rXyn11A with results shown in Fig. 7. HWxyn11 contains a short peptide that spatially distances itself from adjacent ß-sheet. Additionally, its length and rotation angle of ß-sheet positioned similarly in N-terminus of rXyn11A are different. Molecular dynamics simulations analysis also showed greater fluctuations in the N-terminus of HWxyn11, suggesting that the presence of this peptide could impact dynamic stability of protein conformation. The substantial N-terminal structural differences caused by this peptide may be the reason for the accelerated unfolding and inferior thermostability of HWxyn11 under high temperatures.
Therefore, we chose to modify N-terminal and C-terminal regions of HWxyn11 for constructing mutants. Two types of mutants were constructed. HWxAs were composed of HWxyn11 amino acid sequence with different lengths of N-terminal sequence of rXyn11A (Fig. 6a). HWxyn11-CBMs were constructed by fusing different CBMs from thermophilic bacteria at C-terminus.
Temperature characteristics of mutants
All mutants were successfully expressed and purified to apparent homogeneity by Ni-NTA affinity chromatography (Fig. 8). The effect of temperature on enzyme activity of mutants are presented in Fig. 9.
The optimal temperature for both HWxyn11 and HWxAs was 45 ℃ (Fig. 9a). HWxA31 showed over 60% relative enzyme activity and HWxA40 had 30% after incubation at 55 ℃ for 60 min. The wild-type HWxyn11 had 19% residual enzyme activity. In contrast, HWxA6, HWxA15 and HWxA26 showed lower thermostability compared with wild-type HWxyn11 (Fig. 9b). The half-life of HWxA31 (77.9 min) was approximately 3.4-fold higher than that of HWxyn11 (23.0 min), indicating that HWxA31 was more stable than HWxyn11 at 55 ℃.
The optimal temperature for both HWxyn11 and HWxyn11-CBMs was also 45 ℃ (Fig. 9c). And the residual enzyme activity of HWxyn11-CBM 6 exhibited over 50% at 55 ℃ for 60 min, other CBMs contributed to thermostability of HWxyn11 with varying degrees (Fig. 9d). The half-life of HWxyn11-CBM 6 was 65.4 min.
The combined mutant HWxA31-CBM 6, constructed by both N-terminal replacement and C-terminal fusion strategies, showed an optimal reaction temperature of 50 ℃ (Fig. 9e). Its residual enzyme activity was 65% after incubation at same condition as above xylanases, proving the thermostability of combined mutant is better than that of other mutants (Fig. 9f). Additionally, the half-life of HWxA31-CBM 6 (111.8 min) represented an 4.8-fold higher than that of HWxyn11 (23.0 min), showing that HWxA31-CBM 6 has better stability than other mutants at 55 ℃.
Molecular dynamics simulations analysis of mutants
To obtain deep insight into the influence of mutants on structure, molecular dynamics simulations of HWxA31, HWxyn11-CBM 6 and HWxA31-CBM 6 were done at 328 K (55 ℃) (Fig. 10). RMSF represents the influence of fluctuation of single amino acid residue on the whole kinetics of protein and RMSD value refers to root mean square deviation of each atom between protein conformation. Initial conformation under simulated high temperature can be used to evaluate thermostability of protein(Malau and Sianturi 2017). Both RMSF and RMSD are the parameters to characterize structural stability of protein. The lower values of RMSF and RMSD, the stronger structural rigidity and the higher heat resistance of protein. As shown in Fig. 10, RMSD and RMSF of mutants were lower than HWxyn11 at 55 ℃, especially combined mutant of HWxA31-CBM 6, indicating that declined flexibility of these areas. Molecular dynamics simulations analysis showed that the thermostability of mutants by N-terminal replacement and C-terminal CBM fusion were better than wild-type HWxyn11. Our findings were consistent with previous thermostability analysis.