Hydrophobic interaction network analysis for thermostabilization of a mesophilic xylanase

Improvement of thermostability and catalytic efficiency of xylanase from Myceliophthora thermophilar by N-terminal and C-terminal truncation

Introduction

Extracting xylanase from thermophilic filamentous fungi is a feasible way to obtain xylanase with good thermal stability.

Methods

The transcriptomic data of Myceliophthora thermophilic destructive ATCC42464 were differentially expressed and enriched. By comparing the sequences of Mtxylan2 and more than 10 xylanases, the N-terminal and C-terminal of Mtxylan2 were truncated, and three mutants 28N, 28C and 28NC were constructed.

Results and discussion

GH11 xylan Mtxylan2 was identified by transcriptomic analysis, the specific enzyme activity of Mtxylan2 was 104.67 U/mg, and the optimal temperature was 65°C. Molecular modification of Mtxylan2 showed that the catalytic activity of the mutants was enhanced. Among them, the catalytic activity of 28C was increased by 9.3 times, the optimal temperature was increased by 5°C, and the residual enzyme activity remained above 80% after 30 min at 50-65°C, indicating that redundant C-terminal truncation can improve the thermal stability and catalytic performance of GH11 xylanase.

Shifting the pH profiles of Staphylococcus epidermidis lipase (SEL) and Staphylococcus hyicus lipase (SHL) through generating chimeric lipases by DNA shuffling strategy

Enzymes are often required to function in a particular reaction condition by the industrial procedure. In order to identify critical residues affecting the optimum pH of Staphylococcal lipases, chimeric lipases from homologous lipases were generated via a DNA shuffling strategy. Chimeric 1 included mutations of G166S, K212E, T243A, H271Y. Chimeric 2 consisted of substitutions of K212E, T243A, H271Y. Chimeric 3 contained substitutions of K212E, R359L. From the screening results, the pH profiles for chimeric 1 and 2 lipases were shifted from pH 7 to 6. While the pH of chimeric 3 was shifted to 8. It seems the mutation of K212E in chimeric 1 and 2 decreased the pH to 6 by changing the electrostatic potential surface. Furthermore, chimeric 3 showed 10 ˚C improvement in the optimum temperature due to the rigidification of the catalytic loop through the hydrophobic interaction network. Moreover, the substrate specificity of chimeric 1 and 2 was increased towards the longer carbon length chains due to the mutation of T243A adjacent to the lid region through increasing the flexibility of the lid. Current study illustrated that directed evolution successfully modified lipase properties including optimum pH, temperature and substrate specificity through mutations, especially near catalytic and lid regions.

Characterization, Structural Analysis, and Thermal Stability Mutation of a New Zearalenone-Degrading Enzyme Mined from Bacillus subtilis

Zearalenone (ZEN) is a widespread mycotoxin that causes serious damage to animal husbandry and poses a threat to human health. A screen of ZEN-degrading soil bacteria yielded Bacillus subtilis YT-4, which yielded 80% ZEN degradation after 6 h and 95% after 36 h. The gene sequence encoding the degradative enzyme ZENY was mined from the genome of YT-4 and expressed in yeast. ZENY is an α/β-hydrolase with an optimal enzyme activity at 37 °C and pH 8. By breaking the lactone ring of ZEN, it produces ZENY-C18H24O5 with a molecular weight of 320.16 g/mol. Sequence comparison and molecular docking analyses identified the catalytic ZENY triad 99S-245H-123E and the primary ZEN-binding mode within the hydrophobic pocket of the enzyme. To improve the thermal stability of the enzyme for industrial applications, we introduced a mutation at the N-terminus, specifically replacing the fifth residue N with V, and achieved a 25% improvement in stability at 45 °C. These findings aim to achieve ZEN biodegradation and provide insight into the structure and function of ZEN hydrolases.

Computer-Aided Reconstruction and Application of Bacillus halodurans S7 Xylanase with Heat and Alkali Resistance

β-1,4-Endoxylanase is the most critical hydrolase for xylan degradation during lignocellulosic biomass utilization. However, its poor stability and activity in hot and alkaline environments hinder its widespread application. In this study, BhS7Xyl from Bacillus halodurans S7 was improved using a computer-aided design through isothermal compressibility (βT) perturbation engineering and by combining three thermostability prediction algorithms (ICPE-TPA). The best variant with remarkable improvement in specific activity, heat resistance (70 °C), and alkaline resistance (both pH 9.0 and 70 °C), R69F/E137M/E145L, exhibited a 4.9-fold increase by wild-type in specific activity (1368.6 U/mg), a 39.4-fold increase in temperature half-life (458.1 min), and a 57.6-fold increase in pH half-life (383.1 min). Furthermore, R69F/E137M/E145L was applied to the hydrolysis of agricultural waste (corncob and hardwood pulp) to efficiently obtain a higher yield of high-value xylooligosaccharides. Overall, the ICPE-TPA strategy has the potential to improve the functional performance of enzymes under extreme conditions for the high-value utilization of lignocellulosic biomass.

Analyzing Current Trends and Possible Strategies to Improve Sucrose Isomerases' Thermostability

Due to their ability to produce isomaltulose, sucrose isomerases are enzymes that have caught the attention of researchers and entrepreneurs since the 1950s. However, their low activity and stability at temperatures above 40 °C have been a bottleneck for their industrial application. Specifically, the instability of these enzymes has been a challenge when it comes to their use for the synthesis and manufacturing of chemicals on a practical scale. This is because industrial processes often require biocatalysts that can withstand harsh reaction conditions, like high temperatures. Since the 1980s, there have been significant advancements in the thermal stabilization engineering of enzymes. Based on the literature from the past few decades and the latest achievements in protein engineering, this article systematically describes the strategies used to enhance the thermal stability of sucrose isomerases. Additionally, from a theoretical perspective, we discuss other potential mechanisms that could be used for this purpose.

Improved thermostability of type I pullulanase from Bacillus thermoliquefaciens by error-prone PCR

Pullulanase (PulB) is a starch-debranching enzyme. In order to improve its catalytic performance, random mutagenesis was performed on the pullulanase gene derived from Bacillus thermoliquefaciens. Two rounds of error-prone PCR were carried out. Mutant T252S was screened in the first round of error-prone library, which had the highest catalytic activity. During the second round of mutations, mutant enzyme G250P/T252S/G253T/N255K was screened, which had further improved catalytic activity and the best thermostability. Compared with the parent enzyme, the specific activity of mutant enzyme G250P/T252S/G253T/N255K increased by 1.9 times, K_m decreased by 22.7 %, k_cat increased by 28.7 %, and k_cat/K_m increased by 68.4 %. The thermostability of the mutant enzyme improved significantly, showing that the half-life at 60 °C was extended to 7.5 h, which was 87.5 % higher than that of the parent enzyme. The mutation sites in these two rounds were concentrated in the 250-255 regions, indicating that this region was an important region affecting the catalytic activity and Thermostability. The reasons for the change of enzymtic properties was also preliminarily analyzed through three-dimensional simulation.

Three Molecular Modification Strategies to Improve the Thermostability of Xylanase XynA from Streptomyces rameus L2001

Glycoside hydrolase family 11 (GH11) xylanases are the preferred candidates for the production of functional oligosaccharides. However, the low thermostability of natural GH11 xylanases limits their industrial applications. In this study, we investigated the following three strategies to modify the thermostability of xylanase XynA from Streptomyces rameus L2001 mutation to reduce surface entropy, intramolecular disulfide bond construction, and molecular cyclization. Changes in the thermostability of XynA mutants were analyzed using molecular simulations. All mutants showed improved thermostability and catalytic efficiency compared with XynA, except for molecular cyclization. The residual activities of high-entropy amino acid-replacement mutants Q24A and K104A increased from 18.70% to more than 41.23% when kept at 65 °C for 30 min. The catalytic efficiencies of Q24A and K143A increased to 129.99 and 92.26 mL/s/mg, respectively, compared with XynA (62.97 mL/s/mg) when using beechwood xylan as the substrate. The mutant enzyme with disulfide bonds formed between Val3 and Thr30 increased the t_1/2^{60 °C} by 13.33-fold and the catalytic efficiency by 1.80-fold compared with the wild-type XynA. The high thermostabilities and hydrolytic activities of XynA mutants will be useful for enzymatic production of functional xylo-oligosaccharides.

Engineering the kinetic stability of a β-trefoil protein by tuning its topological complexity

Kinetic stability, defined as the rate of protein unfolding, is central to determining the functional lifetime of proteins, both in nature and in wide-ranging medical and biotechnological applications. Further, high kinetic stability is generally correlated with high resistance against chemical and thermal denaturation, as well as proteolytic degradation. Despite its significance, specific mechanisms governing kinetic stability remain largely unknown, and few studies address the rational design of kinetic stability. Here, we describe a method for designing protein kinetic stability that uses protein long-range order, absolute contact order, and simulated free energy barriers of unfolding to quantitatively analyze and predict unfolding kinetics. We analyze two β-trefoil proteins: hisactophilin, a quasi-three-fold symmetric natural protein with moderate stability, and ThreeFoil, a designed three-fold symmetric protein with extremely high kinetic stability. The quantitative analysis identifies marked differences in long-range interactions across the protein hydrophobic cores that partially account for the differences in kinetic stability. Swapping the core interactions of ThreeFoil into hisactophilin increases kinetic stability with close agreement between predicted and experimentally measured unfolding rates. These results demonstrate the predictive power of readily applied measures of protein topology for altering kinetic stability and recommend core engineering as a tractable target for rationally designing kinetic stability that may be widely applicable.

Engineering non-conserved salt bridges in GH11 xylanase from Bacillus pumilus SSP34 for improved thermal stability: an in-silico evaluation

GH11 xylanases are commercially important enzymes for degradation of xylan fibers. We have identified the presence of nine non-conserved and five conserved salt bridges in GH11 xylanase from Bacillus pumilus SSP34. We have designed two sets of mutants viz., (1) substitution mutants in which non-conserved charged amino acid residues have been replaced with appropriate hydrophobic residues based on side chain occupancy and hydrophobicity and (2) deletion mutants where non-conserved charged residues have been deleted. The stability of the mutants has been evaluated in-silico by analyzing the contributions of non-covalent interactions like hydrophobic interaction clusters and salt bridges. The stability of the resultant mutants was evaluated using parameters such as radius of gyration, solvent accessible surface area, root mean square deviation, root mean square fluctuations and protein unfolding measurements using molecular dynamic simulations. The deletion of certain charged residues resulted in mutants having lowered radius of gyration and decreased surface areas. However, RMSD and RMSF measurements indicated lowered stability in comparison to substitution mutants. Of the substitution mutants, the SBM 3 was the most stable mutant as indicated by Rg, SASA, RMSF and simulated protein unfolding measurements. The major contributing factors for improved stability could be strengthening of hydrophobic interactions in the GH11 xylanase from B. pumilus. These in-silico stability measurements of salt bridge mutants may lead to better design of GH11 xylanases for commercial applications.Communicated by Ramaswamy H. Sarma.

Improving the catalytic performance of xylanase from Bacillus circulans through structure-based rational design

Endo-1,4-β-xylanase is one of the most important enzymes employed in biorefineries for obtaining fermentable sugars from hemicellulosic components. Herein, we aimed to improve the catalytic performance of Bacillus circulans xylanase (Bcx) using a structure-guided rational design. A systematic analysis of flexible motions revealed that the R49 component of Bcx (i) constrains the global conformational changes essential for substrate binding and (ii) is involved in modulating flexible motion. Site-saturated mutagenesis of the R49 residue led to the engineering of the active mutants with the trade-off between flexibility and rigidity. The most active mutant R49N improved the catalytic performance, including its catalytic efficiency (7.51-fold), conformational stability (0.7 °C improvement), and production of xylose oligomers (2.18-fold higher xylobiose and 1.72-fold higher xylotriose). The results discussed herein can be applied to enhance the catalytic performance of industrially important enzymes by controlling flexibility.

Aspartic protease-pepstatin A interactions: Structural insights on the thermal inactivation mechanism

Aspartic proteases are the targets for structure-based drug design for their role in physiological processes and pharmaceutical applications. Structural insights into the thermal inactivation mechanism of an aspartic protease in presence and absence of bound pepstatin A have been obtained by kinetics of thermal inactivation, CD, fluorescence spectroscopy and molecular dynamic simulations. The irreversible thermal inactivation of the aspartic protease comprised of loss of tertiary and secondary structures succeeded by the loss of activity, autolysis and aggregation The enthalpy and entropy of thermal inactivation of the enzyme in presence of pepstatin A increased from 81.2 to 148.5 kcal mol^-1, and from 179 to 359 kcal mol^-1 K^-1 respectively. Pepstatin A shifted the mid-point of thermal inactivation of the protease from 58 °C to 77 °C. The association constant (K) for pepstatin A with aspartic protease was 2.5 ± 0.3 × 10 ⁵ M^-1 and ΔG^o value was -8.3 kcal mol^-1. Molecular dynamic simulation studies were able to delineate the role of pepstatin A in stabilizing backbone conformation and side chain interactions. In the Cα-backbone, the short helical segments and the conserved glycines were part of the most unstable segments of the protein. Understanding the mechanism of thermal inactivation has the potential to develop re-engineered thermostable proteases.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Engineering the Cad pathway in Escherichia coli to produce glutarate from L-lysine

For the efficient industrial production of glutarate, an important C5 platform chemical that is widely used in the chemical and pharmaceutical industries, a five-enzyme cascade pathway was designed and reconstructed in vitro to synthesize glutarate from L-lysine. Then, the imbalanced enzyme expression levels of L-lysine decarboxylase from Escherichia coli (EcCA), putrescine aminotransferase (KpcPA) and γ-aminovaleraldehyde dehydrogenase (KpcPD) from Klebsiella pneumoniae, and the poor catalytic efficiency of KpcPA were identified as the rate-limiting bottlenecks. To this end, ribosome binding site regulation was employed to coordinate the enzyme molar ratio of EcCA:KpcPA:KpcPD at approximately 4:8:7 (the optimum ratio obtained in vitro), and volume scanning and hydrophobicity scanning were applied to increase KpcPA activity toward cadaverine from 15.89 ± 0.52 to 75.87 ± 1.51 U·mg^-1. Furthermore, the extracellular accumulation of 5-aminovalerate (5AVA) was considerably reduced by overexpressing gabP encoding the 5AVA importer. Combining these strategies into the engineered strain Glu-02, 77.62 g/L glutarate, the highest titer by E. coli to date, was produced from 100 g/L L-lysine in 42 h, with a yield and productivity of 0.78 g/g L-lysine and 1.85 g/L/h, respectively, at a 5-L scale. The results presented here provide a novel and potential enzymatic process at industrial-scale to produce glutarate from cheaper amino acids. KEY POINTS: • The bioconversion of l-lysine to glutarate using the Cad pathway was first achieved. • Enhancing the conversion efficiency of the Cad route maximizes glutarate in E. coli. • Achieving the highest titer of glutarate by E. coli to date.

Design of mutants of GH11 xylanase from Bacillus pumilus for enhanced stability by amino acid substitutions in the N-terminal region: an in silico analysis

GH11 xylanases are versatile small-molecular-weight single-polypeptide chain monofunctional enzymes. This family of glycoside hydrolases has important applications in food, feed and chemical industries. We designed mutants for improved thermal stability with substitutions in the first six residues of the N-terminal region and evaluated the stability in silico. The first six residues RTITNN of native xylanase have been mutated accordingly to introduce β structure, increase hydrophobic clusters and enhance conformational rigidity in the molecule. To design stable mutants, the approach consisted of constructing root mean square fluctuation (RMSF) plots of both mesophilic and thermophilic xylanases to check the localized backbone displacement maxima, identify the hydrophobic interaction cluster in and around the peaks of interest, construct mutants by substituting appropriate residues based on beta propensity, hydrophobicity, side chain occupancy and conformational rigidity. This resulted in the decreased number of possible substitutions from 19 to 6 residues. Introduction of conformational rigidity by substitution of asparagine residues at 5th and 6th residue position with proline and valine enhanced the stability. Deletion of N-terminal region increased the stability probably by reducing entropic factors. The structure and stability of GH11 xylanase and resultant mutants were analyzed by root mean square deviation, RMSF, radius of gyration and solvent accessible surface area analysis. The stability of the mutants followed the order N-del > Y1P5 >Y1V5 > ATRLM. The contribution of N-terminal end to overall stability of the molecule is significant because of the proximity of the C-terminal end to the N-terminal end which reinforces long-range interactions. Communicated by Ramaswamy H. Sarma.

Tailoring an aldo-keto reductase KmAKR for robust thermostability and catalytic efficiency by stepwise evolution and structure-guided consensus engineering

t-Butyl 6-cyano-(3R,5R)-dihydroxyhexanoate ((3R,5R)-2) is an advanced chiral diol intermediate of the cholesterol-lowering drug atorvastatin. KmAKR_M5 (W297H/Y296W/K29H/Y28A/T63M) constructed in our previous work, displayed good biocatalytic performance on (3R,5R)-2. In the present work, stepwise evolution was applied to further enhance the thermostability and activity of KmAKR_M5. For thermostability enhancement, N109 and S196 located far from the active site were picked out by structure-guided consensus engineering, and mutated by site-directed mutagenesis (SDM). For catalytic efficiency improvement, the residues A30 and T302 adjacent to the substrate-binding pocket were subjected to site-saturation mutagenesis (SSM). As a result, the "best" mutant KmAKR_M9 (W297H/Y296W/K29H/Y28A/T63M/A30P/T302S/N109K/S196C) was developed, of which T₅₀¹⁵ and T_m were 5.0 °C and 8.2 °C higher than those of KmAKR_M5. Moreover, compared to KmAKR_M5, KmAKR_M9 displayed a 1.9-fold (846 vs 2436 min) and 6.7-fold (126 vs 972 min) longer half-lives at 40 and 50 °C, respectively. Structural analysis suggested that beneficial mutations introduced additional hydrophobic interactions and hydrogen bonds, contributing rigidification of the flexible loops and the increase of internal forces, hence increasing the thermostability and activity. 5 g DCW (dry cell weight) L^-1KmAKR_M9 completely reduced 350 g L^-1t-butyl 6-cyano-(5R)-hydroxy-3-oxo-hexanoate ((5R)-1), within 3.7 h at 40 °C, yielding optically pure (3R,5R)-2 (d.e._p > 99.5%) with a space-time yield (STY) of 1.82 kg L^-1 d^-1. Hence, KmAKR_M9 is a robust biocatalyst for the synthesis of (3R,5R)-2.

Enzymatic characterization of a novel thermostable and alkaline tolerant GH10 xylanase and activity improvement by multiple rational mutagenesis strategies

Thermo-alkaline xylanases are widely applied in paper pulping industry. In this study, a novel thermostable and alkaline tolerant GH10 xylanase (Xyn30Y5) gene from alkaliphilic Bacillus sp. 30Y5 was cloned and the surface-layer homology (SLH) domains truncated enzyme (Xyn30Y5-SLH) was expressed in Escherichia coli. The purified Xyn30Y5-SLH was most active at 70 °C and pH 7.0 and showed the highest specific activity of 349.4 U mg^-1. It retained more than 90% activity between pH 6.0 to 9.5 and was stable at pH 6.0-10.0. To improve the activity, 47 mutants were designed based on eight rational strategies and 21 mutants showed higher activity. By combinatorial mutagenesis, the best mutant 3B demonstrated specific activity of 1016.8 U mg^-1 with a doubled catalytic efficiency (k_cat/K_m) and RA60_1/2h value, accompanied by optimal pH shift to 8.0. The molecular dynamics simulation analysis indicated that the increase of flexibility of α5 helix and loop7 located near to the catalytic residues is likely responsible for its activity improvement. And the decrease of flexibility of the most unstable regions is vital for the thermostablity improvement. This work provided not only a novel thermostable and alkaline tolerant xylanase with industrial application potential but also an effective mutagenesis strategy for xylanase activity improvement.

Exploiting the activity-stability trade-off of glucose oxidase from Aspergillus niger using a simple approach to calculate thermostability of mutants

Glucose oxidase (Gox) is a biocatalyst that is widely applied in the food industry, as well as other biotechnological industries. However, the industrial application of Gox is hampered by its low thermostability and activity. Here, we aimed to improve the thermostability of GoxM4 from Aspergillus niger without reducing its activity due to the activity-stability trade-off. A simple and effective approach combining enzyme activity and structure stability was adopted to evaluate the thermostability of GoxM4 and its mutants. After four rounds of computer-aided rational design, the best mutant, GoxM8, was obtained. The melting temperature (T_m) of GoxM8 was increased by 9 °C compared with GoxM4. The catalytic efficiency of GoxM8 was similar to GoxM4, suggesting that the enzyme activity-stability trade-off was counteracted. To explore its mechanism, we performed molecular dynamics simulations of GoxM4 and its mutants. Our findings provided a typical example for researching the enzyme activity-stability trade-off.

Improvement of BsAPA Aspartic Protease Thermostability via Autocatalysis-Resistant Mutation

An aspartic protease gene (Bsapa) was cloned from Bispora sp. MEY-1 and expressed in Pichia pastoris. The recombinant BsAPA showed maximal activity at pH 3.0 and 75 °C and remained stable at 70 °C and below, indicating the thermostable nature of BsAPA. However, heat inactivation still limits the application of BsAPA. To further improve its thermostability, an autocatalysis site (L205-F206) in BsAPA was identified and three mutants (F193W, K204P, and A371V) were generated based on the analysis of the structure neighboring the autocatalysis site. These mutants have improved thermostability, and their half-life at 75 °C increased by 0.5-, 0.2-, and 0.3-fold, respectively. A triple-site mutant (F193W/K204P/A371V) was generated, with 1.5-fold increased half-life at 80 and a 10.7 °C increased T_m, compared with those of the wild-type. These results indicate that autocatalysis of aspartic protease reduces enzyme thermostability. Furthermore, site-directed mutagenesis at regions near the autocatalysis site is an efficient approach to improve aspartic protease thermostability.

Mutagenesis of N-terminal residues confer thermostability on a Penicillium janthinellum MA21601 xylanase

BACKGROUND

A mesophilic xylanase PjxA from Penicillium janthinellum MA21601 has high specific activity under acidic condition and holds great potential for applications in the animal feed industry. To enhance the thermostability of xylanase PjxA, two mutation strategies in the N-terminal region were examined and then integrated into the xylanase to further improvement. The recombinant xylanase PTxA-DB (The meaning of DB is disulfide-bridge.) was constructed by replacement of five residues in the mutated region in TfxA (T10Y, N11H, N12D, Y15F, N30 L), combined with an additional disulfide bridge in the N-terminal region.

RESULTS

The T_m value of mutant PTxA-DB was improved from 21.3 °C to 76.6 °C, and its half-life was found to be 53.6 min at 60 °C, 107-fold higher than the wild type strain. The location of the disulfide bridge (T2C-T29C) was between the irregular loop and the β-strand A2, accounting for most of the improvement in thermostability of PjxA. Further analysis indicated T2C, T29C, N30 L and Y15F lead to increase N-terminal hydrophobicity. Moreover, the specific activity and substrate affinity of PTxA-DB were also enhanced under the acidic pH values.

CONCLUSIONS

These results indicated PTxA-DB could be a prospective additive to industrial animal feeds.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Improvement in the thermostability of a type A feruloyl esterase, AuFaeA, from Aspergillus usamii by iterative saturation mutagenesis

Feruloyl or ferulic acid esterase (Fae, EC 3.1.1.73) catalyzes the hydrolysis of ester bonds between polysaccharides and phenolic acid compounds in xylan side chain. In this study, the thermostability of a type A feruloyl esterase (AuFaeA) from Aspergillus usamii was increased by iterative saturation mutagenesis (ISM). Two amino acids, Ser33 and Asn92, were selected for saturation mutagenesis according to the B-factors analyzed by B-FITTER software and ΔΔG values predicted by PoPMuSiC algorithm. After screening the saturation mutagenesis libraries constructed in Pichia pastoris, 15 promising variants were obtained. The best variant S33E/N92-4 (S33E/N92R) produced a T m value of 44.5 °C, the half-lives (t1/2) of 35 and 198 min at 55 and 50 °C, respectively, corresponding to a 4.7 °C, 2.33- and 3.96-fold improvement compared to the wild type. Additionally, the best S33 variant S33-6 (S33E) was thermostable at 50 °C with a t1/2 of 82 min, which was 32 min longer than that of the wild type. All the screened S33E/N92 variants were more thermostable than the best S33 variant S33-6 (S33E). This work would contribute to the further studies on higher thermostability modification of type A feruloyl esterases, especially those from fungi. The thermostable feruloyl esterase variants were expected to be potential candidates for industrial application in prompting the enzymic degradation of plant biomass materials at elevated temperatures.

Efficient CO2-reducing activity of NAD-dependent formate dehydrogenase from Thiobacillus sp. KNK65MA for formate production from CO2 gas

NAD-dependent formate dehydrogenase (FDH) from Candida boidinii (CbFDH) has been widely used in various CO₂-reduction systems but its practical applications are often impeded due to low CO₂-reducing activity. In this study, we demonstrated superior CO₂-reducing properties of FDH from Thiobacillus sp. KNK65MA (TsFDH) for production of formate from CO₂ gas. To discover more efficient CO₂-reducing FDHs than a reference enzyme, i.e. CbFDH, five FDHs were selected with biochemical properties and then, their CO₂-reducing activities were evaluated. All FDHs including CbFDH showed better CO₂-reducing activities at acidic pHs than at neutral pHs and four FDHs were more active than CbFDH in the CO₂ reduction reaction. In particular, the FDH from Thiobacillus sp. KNK65MA (TsFDH) exhibited the highest CO₂-reducing activity and had a dramatic preference for the reduction reaction, i.e., a 84.2-fold higher ratio of CO₂ reduction to formate oxidation in catalytic efficiency (k_cat/K_B) compared to CbFDH. Formate was produced from CO₂ gas using TsFDH and CbFDH, and TsFDH showed a 5.8-fold higher formate production rate than CbFDH. A sequence and structural comparison showed that FDHs with relatively high CO₂-reducing activities had elongated N- and C-terminal loops. The experimental results demonstrate that TsFDH can be an alternative to CbFDH as a biocatalyst in CO₂ reduction systems.

Structure of a His170Tyr mutant of thermostable pNPPase from Geobacillus stearothermophilus

Using directed evolution based on random mutagenesis and heat-treated selection, a thermostable His170Tyr mutant of Geobacillus stearothermophilus thermostable p-nitrophenylphosphatase (TpNPPase) was obtained. The temperature at which the His170Tyr mutant lost 50% of its activity (T1/2) was found to be 4.40 K higher than that of wild-type TpNPPase, and the melting temperature of the His170Tyr mutant increased by 2.39 K. The crystal structure of the His170Tyr mutant was then determined at 2.0 Å resolution in the presence of a sodium ion and a sulfate ion in the active site. The cap domain of chain B shows a half-closed conformation. The hydrophobic side chain of the mutated residue, the hydroxyphenyl group, forms a hydrophobic contact with the methyl group of Ala166. This hydrophobic interaction was found using the Protein Interactions Calculator (PIC) web server with an interaction distance of 4.6 Å, and might be a key factor in the thermostabilization of the His170Tyr mutant. This study potentially offers a molecular basis for both investigation of the catalytic mechanism and thermostable protein engineering.

Enhancing the activity of Bacillus circulans xylanase by modulating the flexibility of the hinge region

Enzymes undergo multiple conformational changes in solution, and these dynamics are considered to play a critical role in enzyme activity. Hinge-bending motions, resulting from reciprocal movements of dynamical quasi-rigid bodies, are thought to be related to turnover rate and are affected by the physical properties of the hinge regions. In this study, hinge identification and flexibility modification of the regions by mutagenesis were conducted to explore the relationship between hinge flexibility and catalytic activity. Bacillus circulans xylanase was selected for the identification and mutation of the hinge regions. As a result, turnover rate (V(max)) was improved approximately twofold in mutants that have more rigid hinge structure, despite the decrease in K(m) and V(max)/K(m). This result indicates that the rigidly mutated hinge has positive effects on B. circulans xylanase activity.

Quality matters: extension of clusters of residues with good hydrophobic contacts stabilize (hyper)thermophilic proteins

Identifying determinant(s) of protein thermostability is key for rational and data-driven protein engineering. By analyzing more than 130 pairs of mesophilic/(hyper)thermophilic proteins, we identified the quality (residue-wise energy) of hydrophobic interactions as a key factor for protein thermostability. This distinguishes our study from previous ones that investigated predominantly structural determinants. Considering this key factor, we successfully discriminated between pairs of mesophilic/(hyper)thermophilic proteins (discrimination accuracy: ∼80%) and searched for structural weak spots in E. coli dihydrofolate reductase (classification accuracy: 70%).

Computational tools for designing and engineering enzymes

Protein engineering strategies aimed at constructing enzymes with novel or improved activities, specificities, and stabilities greatly benefit from in silico methods. Computational methods can be principally grouped into three main categories: bioinformatics; molecular modelling; and de novo design. Particularly de novo protein design is experiencing rapid development, resulting in more robust and reliable predictions. A recent trend in the field is to combine several computational approaches in an interactive manner and to complement them with structural analysis and directed evolution. A detailed investigation of designed catalysts provides valuable information on the structural basis of molecular recognition, biochemical catalysis, and natural protein evolution.

Improvement in thermostability of metagenomic GH11 endoxylanase (Mxyl) by site-directed mutagenesis and its applicability in paper pulp bleaching process

An attempt has been made for enhancing the thermostability of xylanase (Mxyl) retrieved from a compost-soil-based metagenomic library. The analysis of the structure of xylanase by molecular dynamics simulation revealed more structural fluctuations in β-sheets. When the surface of β-sheets was enriched with arginine residues by substituting serine/threonine by site-directed mutagenesis, the enzyme with four arginine substitutions (MxylM4) exhibited enhanced thermostability at 80 °C. The T 1/2 of MxylM4 at 80 °C, in the presence of birchwood xylan, increased from 130 to 150 min at 80 °C without any alteration in optimum pH and temperature and molecular mass. Improvement in thermostability of MxylM4 was corroborated by increase in T m by 6 °C over that of Mxyl. The K m of MxylM4, however, increased from 8.01 ± 0.56 of Mxyl to 12.5 ± 0.32 mg ml(-1), suggesting a decrease in the affinity as well as specific enzyme activity. The Mxyl as well as MxylM4 liberated chromophores and lignin-derived compounds from kraft pulp, indicating their applicability in pulp bleaching.

Structure of the Aeropyrum pernix L7Ae multifunctional protein and insight into its extreme thermostability

Archaeal ribosomal protein L7Ae is a multifunctional RNA-binding protein that directs post-transcriptional modification of archaeal RNAs. The L7Ae protein from Aeropyrum pernix (Ap L7Ae), a member of the Crenarchaea, was found to have an extremely high melting temperature (>383 K). The crystal structure of Ap L7Ae has been determined to a resolution of 1.56 Å. The structure of Ap L7Ae was compared with the structures of two homologs: hyperthermophilic Methanocaldococcus jannaschii L7Ae and the mesophilic counterpart mammalian 15.5 kD protein. The primary stabilizing feature in the Ap L7Ae protein appears to be the large number of ion pairs and extensive ion-pair network that connects secondary-structural elements. To our knowledge, Ap L7Ae is among the most thermostable single-domain monomeric proteins presently observed.

Structure- and sequence-analysis inspired engineering of proteins for enhanced thermostability

Protein engineering strategies for increasing stability can be improved by replacing random mutagenesis and high-throughput screening by approaches that include bioinformatics and computational design. Mutations can be focused on regions in the structure that are most flexible and involved in the early steps of thermal unfolding. Sequence analysis can often predict the position and nature of stabilizing mutations, and may allow the reconstruction of thermostable ancestral sequences. Various computational tools make it possible to design stabilizing features, such as hydrophobic clusters and surface charges. Different methods for designing chimeric enzymes can also support the engineering of more stable proteins without the need of high-throughput screening.