Acidophilic adaptation of family 11 endo-beta-1,4-xylanases: modeling and mutational analysis

Clustered surface amino acid residues modulate the acid stability of GH10 xylanase in fungi

Acidic xylanases are widely used in industries such as biofuels, animal feeding, and fruit juice clarification due to their tolerance to acidic environments. However, the factors controlling their acid stability, especially in GH10 xylanases, are only partially understood. In this study, we identified a series of thermostable GH10 xylanases with optimal temperatures ranging from 70 to 90 °C, and among these, five enzymes (Xyn10C, Xyn10RE, Xyn10TC, Xyn10BS, and Xyn10PC) exhibited remarkable stability at pH 2.0. Our statistical analysis highlighted several factors contributing to the acid stability of GH10 xylanases, including electrostatic repulsion, π-π stacking, ionic bonds, hydrogen bonds, and Van der Waals interactions. Furthermore, through mutagenesis studies, we uncovered that acid stability is influenced by a complex interplay of amino acid residues. The key amino acid sites determining the acid stability of GH10 xylanases were thus elucidated, mainly concentrated in two surface regions behind the enzyme active center. Notably, the critical residues associated with acid stability markedly enhanced Xyn10RE's thermostability by more than sixfold, indicating a potential acid-thermal interplay in GH10 xylanases. This study not only reported a series of valuable genes but also provided a range of modification targets for enhancing the acid stability of GH10 xylanases. KEY POINTS: • Five acid stable and thermostable GH10 xylanases were reported. • The key amino acid sites, mainly forming two enriched surface regions behind the enzyme active center, were identified responsible for acid stability of GH10 xylanases. • The finding revealed interactive amino acid sites, offering a pathway for synergistic enhancement of both acid stability and thermostability in GH10 xylanase modifications.

Strategies for tailoring pH performances of glycoside hydrolases

Glycoside hydrolases (GHs) exhibit high activity and stability under harsh conditions, such as high temperatures and extreme pHs, given their wide use in industrial biotechnology. However, strategies for improving the acidophilic and alkalophilic adaptations of GHs are poorly summarized due to the complexity of the mechanisms of these adaptations. This review not only highlights the adaptation mechanisms of acidophilic and alkalophilic GHs under extreme pH conditions, but also summarizes the recent advances in engineering the pH performances of GHs with a focus on four strategies of protein engineering, enzyme immobilization, chemical modification, and medium engineering (additives). The examples described here summarize the methods used in modulating the pH performances of GHs and indicate that methods integrated in different protein engineering techniques or methods are efficient to generate industrial biocatalysts with the desired pH performance and other adapted enzyme properties.

Improvement of alkalophilicity of an alkaline xylanase Xyn11A-LC from Bacillus sp. SN5 by random mutation and Glu135 saturation mutagenesis

BACKGROUND

Family 11 alkaline xylanases have great potential economic applications in the pulp and paper industry. In this study, we would improve the alkalophilicity of family 11 alkaline xylanase Xyn11A-LC from Bacillus sp. SN5, for the better application in this field.

RESULTS

A random mutation library of Xyn11A-LC with about 10,000 clones was constructed by error-prone PCR. One mutant, M52-C10 (V116A and E135V), with improved alkalophilicity was obtained from the library. Site-directed mutation showed that the mutation E135V was responsible for the alkalophilicity of the mutant. The variant E135V shifted the optimum pH of the wild-type enzyme from 7.5 to 8.0. Compared to the relative activities of the wild type enzyme, those of the mutant E135V increased by 17.5, 18.9, 14.3 and 9.5 % at pH 8.5, 9.0, 9.5 and 10.0, respectively. Furthermore, Glu135 saturation mutagenesis showed that the only mutant to have better alkalophilicity than E135V was E135R. The optimal pH of the mutant E135R was 8.5, 1.0 pH units higher than that of the wild-type. In addition, compared to the wild-type enzyme, the mutations E135V and E135R increased the catalytic efficiency (k _cat/K _m) by 57 and 37 %, respectively. Structural analysis showed that the residue at position 135, located in the eight-residue loop on the protein surface, might improve the alkalophilicity and catalytic activity by the elimination of the negative charge and the formation of salt-bridge.

CONCLUSIONS

Mutants E135V and E135R with improved alkalophilicity were obtained by directed evolution and site saturation mutagenesis. The residue at position 135 in the eight-residue loop on the protein surface was found to play an important role in the pH activity profile of family 11 xylanases. This study provided alkalophilic mutants for application in bleaching process, and it was also helpful to understand the alkaline adaptation mechanism of family 11 xylanases.

Sequence homolog-based molecular engineering for shifting the enzymatic pH optimum

Cell-free synthetic biology system organizes multiple enzymes (parts) from different sources to implement unnatural catalytic functions. Highly adaption between the catalytic parts is crucial for building up efficient artificial biosynthetic systems. Protein engineering is a powerful technology to tailor various enzymatic properties including catalytic efficiency, substrate specificity, temperature adaptation and even achieve new catalytic functions. However, altering enzymatic pH optimum still remains a challenging task. In this study, we proposed a novel sequence homolog-based protein engineering strategy for shifting the enzymatic pH optimum based on statistical analyses of sequence-function relationship data of enzyme family. By two statistical procedures, artificial neural networks (ANNs) and least absolute shrinkage and selection operator (Lasso), five amino acids in GH11 xylanase family were identified to be related to the evolution of enzymatic pH optimum. Site-directed mutagenesis of a thermophilic xylanase from Caldicellulosiruptor bescii revealed that four out of five mutations could alter the enzymatic pH optima toward acidic condition without compromising the catalytic activity and thermostability. Combination of the positive mutants resulted in the best mutant M31 that decreased its pH optimum for 1.5 units and showed increased catalytic activity at pH < 5.0 compared to the wild-type enzyme. Structure analysis revealed that all the mutations are distant from the active center, which may be difficult to be identified by conventional rational design strategy. Interestingly, the four mutation sites are clustered at a certain region of the enzyme, suggesting a potential "hot zone" for regulating the pH optima of xylanases. This study provides an efficient method of modulating enzymatic pH optima based on statistical sequence analyses, which can facilitate the design and optimization of suitable catalytic parts for the construction of complicated cell-free synthetic biology systems.

Structural Insight into and Mutational Analysis of Family 11 Xylanases: Implications for Mechanisms of Higher pH Catalytic Adaptation

To understand the molecular basis of higher pH catalytic adaptation of family 11 xylanases, we compared the structures of alkaline, neutral, and acidic active xylanases and analyzed mutants of xylanase Xyn11A-LC from alkalophilic Bacillus sp. SN5. It was revealed that alkaline active xylanases have increased charged residue content, an increased ratio of negatively to positively charged residues, and decreased Ser, Thr, and Tyr residue content relative to non-alkaline active counterparts. Between strands β6 and β7, alkaline xylanases substitute an α-helix for a coil or turn found in their non-alkaline counterparts. Compared with non-alkaline xylanases, alkaline active enzymes have an inserted stretch of seven amino acids rich in charged residues, which may be beneficial for xylanase function in alkaline conditions. Positively charged residues on the molecular surface and ionic bonds may play important roles in higher pH catalytic adaptation of family 11 xylanases. By structure comparison, sequence alignment and mutational analysis, six amino acids (Glu16, Trp18, Asn44, Leu46, Arg48, and Ser187, numbering based on Xyn11A-LC) adjacent to the acid/base catalyst were found to be responsible for xylanase function in higher pH conditions. Our results will contribute to understanding the molecular mechanisms of higher pH catalytic adaptation in family 11 xylanases and engineering xylanases to suit industrial applications.

Microbial heterotrophic production in an oligotrophic acidic geothermal lake: responses to organic amendments and terrestrial plant litter

Boiling Springs Lake (BSL) is an oligotrophic, acidic geothermal feature where even very low levels of microbial heterotrophic production still exceed autotrophy. To test whether allochthonous leaf litter (LL) inputs fuel this excess, we quantified leaf litterfall, leaching and decomposition kinetics, and measured the impact of organic amendments on production, germination and cell growth, using pyrosequencing to track changes in microbial community composition. Coniferous leaves in BSL exhibited high mass loss rates during leaching and decomposition, likely due to a combination of chemical hydrolysis and contributions of both introduced and endemic microbes. We measured very low in situ (3)H-thymidine incorporation over hours by the dominant chemolithotroph Acidimicrobium (13-65 μg C L(-1) day(-1)), which was inhibited by simple C sources (acetate, glucose). Longer term incubations with additions of 0.01-0.02% complex C/N sources induced germination of the Firmicute Alicyclobacillus within 1-2 days, as well as growth of Acetobacteraceae after 3-4 days. LL additions yielded the opposite successional patterns of these r-selected heterotrophs, boosting production to 30-150 μg C L(-1) day(-1). Growth and germination studies suggest both prokaryotes and fungi likely consume allochthonous organics, and might be novel sources of lignocellulose-degrading enzymes. A model of BSL's C budget supports the hypothesis that allochthonous inputs fuel seasonal microbial heterotrophy, but that dissolved organic C sources greatly exceed direct LL inputs.

Alkalophilic adaptation of XynB endoxylanase from Aspergillus niger via rational design of pKa of catalytic residues

Based on the strategy of changing pH-stability profiles by altering pKa values of catalytic residues, rational protein engineering was applied to improve alkalophilic adaptation of Aspergillus niger endoxylanase XynB. Computational predictions and molecular modeling were carried out using PROPKA server and SWISS-MODEL server, respectively. Three endoxylanase mutant of S108V, N151E, and Q178R, in which the pKa values of either catalytic glutamate residues shifted, were generated. In agreement with expectation, the variant of Q178R improved alkalophilic performances. The mutant Q178R raised the optimum pH of XynB from 5.5 to 6.0 and retained 37% of the maximum activity at pH 8.0. Interestingly, the pKa values of Glu84 and Glu175 in Q178R are 7.91 and 6.32, respectively. The pKa of Glu175 is lower than that of Glu84, as opposed to the fact that the pKa of Glu84 is lower than that of Glu175 in other GH11 xylanases. It indicated that Glu175 may convert into a nucleophile residue and Glu84 into an acid/base residue.

Molecular approaches for ameliorating microbial xylanases

In industrial processes, chemical catalysis is being replaced by enzyme catalysis, since the latter is environmentally benign, non-persistent and cost effective. Microbial xylanases have significant applications in textile, baking, food and feed industries, and in paper and pulp industries for reducing the chlorine requirement. The hazardous chlorine required for bleaching can be reduced up to 25-30% by including an enzymatic step in the pulp bleaching process. The paper pulp bleaching requires xylanases that are active at alkaline pH and elevated temperatures. The enzymes from the cultured microbes do not perform optimally in the paper industry due to their inadequate stability under the process conditions of high temperature and alkaline pH. This review, therefore, deals with the rationale of molecular approaches such as protein engineering for designing xylanases with improved characteristics to suit the process conditions in industries, and prospects and problems.

GH11 xylanases: Structure/function/properties relationships and applications

For technical, environmental and economical reasons, industrial demands for process-fitted enzymes have evolved drastically in the last decade. Therefore, continuous efforts are made in order to get insights into enzyme structure/function relationships to create improved biocatalysts. Xylanases are hemicellulolytic enzymes, which are responsible for the degradation of the heteroxylans constituting the lignocellulosic plant cell wall. Due to their variety, xylanases have been classified in glycoside hydrolase families GH5, GH8, GH10, GH11, GH30 and GH43 in the CAZy database. In this review, we focus on GH11 family, which is one of the best characterized GH families with bacterial and fungal members considered as true xylanases compared to the other families because of their high substrate specificity. Based on an exhaustive analysis of the sequences and 3D structures available so far, in relation with biochemical properties, we assess biochemical aspects of GH11 xylanases: structure, catalytic machinery, focus on their "thumb" loop of major importance in catalytic efficiency and substrate selectivity, inhibition, stability to pH and temperature. GH11 xylanases have for a long time been used as biotechnological tools in various industrial applications and represent in addition promising candidates for future other uses.

Molecular cloning and characterization of the novel acidic xylanase XYLD from Bispora sp. MEY-1 that is homologous to family 30 glycosyl hydrolases

We cloned and sequenced a xylanase gene named xylD from the acidophilic fungus Bispora sp. MEY-1 and expressed the gene in Pichia pastoris. The 1,422-bp full-length complementary DNA fragment encoded a 457-amino acid xylanase with a calculated molecular mass of 49.8 kDa. The mature protein of XYLD showed high sequence similarity to both glycosyl hydrolase (GH) families 5 and 30 but was more homologous to members of GH 30 based on phylogenetic analysis. XYLD shared the highest identity (49.9%) with a putative endo-1,6-beta-D-glucanase from Talaromyces stipitatus and exhibited 21.1% identity and 34.3% similarity to the well-characterized GH family 5 xylanase from Erwinia chrysanthemi. Purified recombinant XYLD showed maximal activity at pH 3.0 and 60 degrees C, maintained more than 60% of maximal activity when assayed at pH 1.5-4.0, and had good thermal stability at 60 degrees C and remained stable at pH 1.0-6.0. The enzyme activity was enhanced in the presence of Ni(2+) and beta-mercaptoethanol and inhibited by some metal irons (Hg(2+), Cu(2+), Pb(2+), Mn(2+), Li(+), and Fe(3+)) and sodium dodecyl sulfate. The specific activity of XYLD for beechwood xylan, birchwood xylan, 4-O-methyl-D-glucuronoxylan, and oat spelt xylan was 2,463, 2,144, 2,020, and 1,429 U mg(-1), respectively. The apparent K (m) and V (max) values for beechwood xylan were 5.6 mg ml(-1) and 3,622 micromol min(-1) mg(-1), respectively. The hydrolysis products of different xylans were mainly xylose and xylobiose.

Xylanase XYL1p from Scytalidium acidophilum: site-directed mutagenesis and acidophilic adaptation

The role of residues Asp60, Tyr35 and Glu141 in the pH-dependent activity of xylanase XYL1p from Scytalidium acidophilum was investigated by site-directed mutagenesis. These amino acids are highly conserved among the acidophilic family 11 xylanases and located near the catalytic site. XYL1p and its single mutants D60N, Y35W and E141A and three combined mutants DN/YW, DN/EA and YW/EA were over-expressed in Pichia pastoris and purified. Xylanase activities at different pH's and temperatures were determined. All mutations increased the pH optimum by 0.5-1.5 pH units. All mutants have lower specific activities except the E141A mutant that exhibited a 50% increase in specific activity at pH 4.0 and had an overall catalytic efficiency higher than the wild-type enzyme. Thermal unfolding experiments show that both the wild-type and E141A mutant proteins have a T(m) maximum at pH 3.5, the E141A mutant being slightly less stable than the wild-type enzyme. These mutations confirm the importance of these amino acids in the pH adaptation. Mutant E141A with its enhanced specific activity at pH 4.0 and improved overall catalytic efficiency is of possible interest for biotechnological applications.

Crystallographic analysis shows substrate binding at the -3 to +1 active-site subsites and at the surface of glycoside hydrolase family 11 endo-1,4-beta-xylanases

GH 11 (glycoside hydrolase family 11) xylanases are predominant enzymes in the hydrolysis of heteroxylan, an abundant structural polysaccharide in the plant cell wall. To gain more insight into the protein-ligand interactions of the glycone as well as the aglycone subsites of these enzymes, catalytically incompetent mutants of the Bacillus subtilis and Aspergillus niger xylanases were crystallized, soaked with xylo-oligosaccharides and subjected to X-ray analysis. For both xylanases, there was clear density for xylose residues in the -1 and -2 subsites. In addition, for the B. subtilis xylanase, there was also density for xylose residues in the -3 and +1 subsite showing the spanning of the -1/+1 subsites. These results, together with the observation that some residues in the aglycone subsites clearly adopt a different conformation upon substrate binding, allowed us to identify the residues important for substrate binding in the aglycone subsites. In addition to substrate binding in the active site of the enzymes, the existence of an unproductive second ligand-binding site located on the surface of both the B. subtilis and A. niger xylanases was observed. This extra binding site may have a function similar to the separate carbohydrate-binding modules of other glycoside hydrolase families.

Structure of a NADH-insensitive hexameric citrate synthase that resists acid inactivation

Acetobacter aceti converts ethanol to acetic acid, and strains highly resistant to both are used to make vinegar. A. aceti survives acetic acid exposure by tolerating cytoplasmic acidification, which implies an unusual adaptation of cytoplasmic components to acidic conditions. A. aceti citrate synthase (AaCS), a hexameric type II citrate synthase, is required for acetic acid resistance and, therefore, would be expected to function at low pH. Recombinant AaCS has intrinsic acid stability that may be a consequence of strong selective pressure to function at low pH, and unexpectedly high thermal stability for a protein that has evolved to function at approximately 30 degrees C. The crystal structure of AaCS, complexed with oxaloacetate (OAA) and the inhibitor carboxymethyldethia-coenzyme A (CMX), was determined to 1.85 A resolution using protein purified by a tandem affinity purification procedure. This is the first crystal structure of a "closed" type II CS, and its active site residues interact with OAA and CMX in the same manner observed in the corresponding type I chicken CS.OAA.CMX complex. While AaCS is not regulated by NADH, it retains many of the residues used by Escherichia coli CS (EcCS) for NADH binding. The surface of AaCS is abundantly decorated with basic side chains and has many fewer uncompensated acidic charges than EcCS; this constellation of charged residues is stable in varied pH environments and may be advantageous in the A. aceti cytoplasm.

Mapping of residues involved in the interaction between the Bacillus subtilis xylanase A and proteinaceous wheat xylanase inhibitors

The Bacillus subtilis xylanase A was subjected to site-directed mutagenesis, aimed at changing the interaction with Triticum aestivum xylanase inhibitor, the only wheat endogenous proteinaceous xylanase inhibitor interacting with this xylanase. The published structure of Bacillus circulans XynA was used to target amino acids surrounding the active site cleft of B.subtilis XynA for mutation. Twenty-two residues were mutated, resulting in 62 different variants. The catalytic activity of active mutants ranged from 563 to 5635 XU/mg and the interaction with T.aestivum xylanase inhibitor showed a similar variation. The results indicate that T.aestivum xylanase inhibitor interacts with several amino acid residues surrounding the active site of the enzyme. Three different amino acid substitutions in one particular residue (D11) completely abolished the interaction between T.aestivum xylanase inhibitor and B.subtilis xylanase A.

Synchrotron radiation small angle scattering studies of thermal stability of xylanase XYNII fromTrichoderma longibrachiatum

Xylanase XYNII from Trichoderma longibrachiatum is a small protein of the molecular weight 21 kDa, belonging to the family 11 of glycosyl hydrolases, which catalyses hydrolysis of xylan. This article reports thermal stability study of xylanase XYN II conformation in the temperature range 15-65 degrees C by the small angle synchrotron radiation scattering. The study has been performed at different pH conditions: at pH 4.0 (below the physiological optimum of the enzyme activity) at pH 5.8 close to the optimum for enzymatic activity and at pH 8.0. The radius of gyration and the pair distance distribution function p(r) have been analyzed to characterize the changes of the enzyme conformation on heating. In the environment of the pH close to that of the optimum for the enzymatic activity, xylanase shows the greatest thermal stability and undergoes denaturation only above 55 degrees C. In the acidic and basic environments, the enzyme stability is much lower and denaturation begins at 45 degrees C. On the basis of the SAXS data, the shape of the xylanase molecule in solution in different temperatures has been reconstructed using ab initio method and program DAMMIN. The shape of the xylanase molecule at room temperature is similar to the right hand, which is typically observed for xylanase crystal structure. In higher temperatures (close to the enzyme activity optimum), the conformation of the right hand is loosened and half opened.

Increased alkali stability in Trichoderma reesei endo-1,4-β-xylanase II by site directed mutagenesis

A number of engineered Trichoderma reesei endo-beta-1,4-xylanase (Xyn II) mutants were created and activity tests were performed for increased stability. The stability of the earlier characterized mutant Y5 (T2C, T28C, K58R, +191D) was further increased by the mutations creating the constructs P9 (N97R+F93W+H144K), P12 (H144C+N92C), P15 (F180Q+H144C+N92C) and P21 (H22K+F180Q+H144C+N92C). The resistance towards thermal inactivation at alkaline pH was increased in all of the mutants. Residual activity T(50%) was increased 4-5 degrees C for P9 at pH 9. The performance of the P9 mutant in sulphate pulp bleaching was also tested and was shown to increase brightness markedly compared to the reference. The bleaching results showed the industrial potential of the obtained mutant.

Improving the alkalophilic performances of the Xyl1 xylanase from Streptomyces sp. S38: structural comparison and mutational analysis

Endo-beta-1,4-xylanases of the family 11 glycosyl-hydrolases are catalytically active over a wide range of pH. Xyl1 from Streptomyces sp. S38 belongs to this family, and its optimum pH for enzymatic activity is 6. Xyn11 from Bacillus agaradhaerens and XylJ from Bacillus sp. 41M-1 share 85% sequence identity and have been described as highly alkalophilic enzymes. In an attempt to better understand the alkalophilic adaptation of xylanases, the three-dimensional structures of Xyn11 and Xyl1 were compared. This comparison highlighted an increased number of salt-bridges and the presence of more charged residues in the catalytic cleft as well as an eight-residue-longer loop in the alkalophilic xylanase Xyn11. Some of these charges were introduced in the structure of Xyl1 by site-directed mutagenesis with substitutions Y16D, S18E, G50R, N92D, A135Q, E139K, and Y186E. Furthermore, the eight additional loop residues of Xyn11 were introduced in the homologous loop of Xyl1. In addition, the coding sequence of the XylJ catalytic domain was synthesized by recursive PCR, expressed in a Streptomyces host, purified, and characterized together with the Xyl1 mutants. The Y186E substitution inactivated Xyl1, but the activity was restored when this mutation was combined with the G50R or S18E substitutions. Interestingly, the E139K mutation raised the optimum pH of Xyl1 from 6 to 7.5 but had no effect when combined with the N92D substitution. Modeling studies identified the possible formation of an interaction between the introduced lysine and the substrate, which could be eliminated by the formation of a putative salt-bridge in the N92D/E139K mutant.

Enzyme adaptation to alkaline pH: atomic resolution (1.08 A) structure of phosphoserine aminotransferase from Bacillus alcalophilus

The crystal structure of the vitamin B(6)-dependent enzyme phosphoserine aminotransferase from the obligatory alkaliphile Bacillus alcalophilus has been determined at 1.08 A resolution. The model was refined to an R-factor of 11.7% (R(free) = 13.9%). The enzyme displays a narrow pH optimum of enzymatic activity at pH 9.0. The final structure was compared to the previously reported structure of the mesophilic phosphoserine aminotransferase from Escherichia coli and to that of phosphoserine aminotransferase from a facultative alkaliphile, Bacillus circulans subsp. alkalophilus. All three enzymes are homodimers with each monomer comprising a two-domain architecture. Despite the high structural similarity, the alkaliphilic representatives possess a set of distinctive structural features. Two residues directly interacting with pyridoxal-5'-phosphate are replaced, and an additional hydrogen bond to the O3' atom of the cofactor is present in alkaliphilic phosphoserine aminotransferases. The number of hydrogen bonds and hydrophobic interactions at the dimer interface is increased. Hydrophobic interactions between the two domains in the monomers are enhanced. Moreover, the number of negatively charged amino acid residues increases on the solvent-accessible molecular surface and fewer hydrophobic residues are exposed to the solvent. Further, the total amount of ion pairs and ion networks is significantly reduced in the Bacillus enzymes, while the total number of hydrogen bonds is increased. The mesophilic enzyme from Escherichia coli contains two additional beta-strands in a surface loop with a third beta-strand being shorter in the structure. The identified structural features are proposed to be possible factors implicated in the alkaline adaptation of phosphoserine aminotransferase.

Functional importance of Asp37 from a family 11 xylanase in the binding to two proteinaceous xylanase inhibitors from wheat

Aspergillus niger xylanase is a target enzyme of the two wheat proteinaceous inhibitors, XIP-I and TAXI-I. We previously suggested that the xylanase "thumb" region was XIP-I binding site. Here, we expressed the Asp37Ala mutant in Pichia pastoris and showed that the mutation abolished the enzyme capacity to interact with both inhibitors, suggesting a direct contact at the active site. The mutant pH profile was altered, confirming the key role of Asp37 in determining the pH optima of glycoside hydrolase family 11. The results are consistent with a competitive inhibition mode and underline the strategic importance of Asp37 in the inhibition mechanism.