Enhanced Thermostability and Catalytic Activity of Streptomyces mobaraenesis Transglutaminase by Rationally Engineering Its Flexible Regions

Navigating the landscape of enzyme design: from molecular simulations to machine learning

Global environmental issues and sustainable development call for new technologies for fine chemical synthesis and waste valorization. Biocatalysis has attracted great attention as the alternative to the traditional organic synthesis. However, it is challenging to navigate the vast sequence space to identify those proteins with admirable biocatalytic functions. The recent development of deep-learning based structure prediction methods such as AlphaFold2 reinforced by different computational simulations or multiscale calculations has largely expanded the 3D structure databases and enabled structure-based design. While structure-based approaches shed light on site-specific enzyme engineering, they are not suitable for large-scale screening of potential biocatalysts. Effective utilization of big data using machine learning techniques opens up a new era for accelerated predictions. Here, we review the approaches and applications of structure-based and machine-learning guided enzyme design. We also provide our view on the challenges and perspectives on effectively employing enzyme design approaches integrating traditional molecular simulations and machine learning, and the importance of database construction and algorithm development in attaining predictive ML models to explore the sequence fitness landscape for the design of admirable biocatalysts.

The Frequently Used Industrial Food Process Additive, Microbial Transglutaminase: Boon or Bane

Microbial transglutaminase (mTG) is a frequently consumed processed food additive, and use of its cross-linked complexes is expanding rapidly. It was designated as a processing aid and was granted the generally recognized as safe (GRAS) classification decades ago, thus avoiding thorough assessment according to current criteria of toxicity and public health safety. In contrast to the manufacturer's declarations and claims, mTG and/or its transamidated complexes are proinflammatory, immunogenic, allergenic, pathogenic, and potentially toxic, hence raising concerns for public health. Being a member of the transglutaminase family and functionally imitating the tissue transglutaminase, mTG was recently identified as a potential inducer of celiac disease. Microbial transglutaminase and its docked complexes have numerous detrimental effects. Those harmful aspects are denied by the manufacturers, who claim the enzyme is deactivated when heated or by gastric acidity, and that its covalently linked isopeptide bonds are safe. The present narrative review describes the potential side effects of mTG, highlighting its thermostability and activity over a broad pH range, thus, challenging the manufacturers' and distributers' safety claims. The national food regulatory authorities and the scientific community are urged to reevaluate mTG's GRAS status, prioritizing public health protection against the possible risks associated with this enzyme and its health-damaging consequences.

Rapid aggregation of amyloid-like protein enhanced by mTGase to prepare functional wool fabrics for efficient and sustainable remove heavy metals from wastewater

To counteract the increasing severity of water pollution and purify water sources, wastewater treatment materials are essential. In particular, it is necessary to improve the bonding strength between the adsorption material and the substrate in a long-term humid environment, and resist the invasion of microorganisms to prolong the service life. In this study, an amyloid-like aggregation method of lysozyme catalyzed by microbial transglutaminase (mTGase). Lysozyme self-assembles into an amyloid-like phase-transited lysozyme (PTL) in the presence of a reducing agent. Simultaneously, mTGase catalyzes acyl transfer reactions within lysozyme molecules or between lysozyme and keratin molecules, and driving PTL assembly on the wool fiber (TG-PTL@wool). This process enhances the grafting amount and fastness of PTL on the wool. Moreover, the tensile strength of wool fabric increased to 523 N. TG-PTL@wool achieves a 97.32 % removal rate of heavy metals, maintaining a removal rate of over 95 % after 5 cycles. TG-PTL@wool has excellent antibacterial property (99 %), and it remains above 90 % after 50 times of circulating washing. This study proved that mTGase can enhance the amyloid aggregation of lysozyme and enhance the bonding strength between PTL coating and substrate. Moreover, TG-PTL@wool provides a sustainable, efficient and cleaner solution for removing heavy metals from water.

The two-step strategy for enhancing the specific activity and thermostability of alginate lyase AlyG2 with mechanism for improved thermostability

To overcome the trade-off challenge encountered in the engineering of alginate lyase AlyG2 from Seonamhaeicola algicola Gy8T and to expand its potential industrial applications, we devised a two-step strategy encompassing activity enhancement followed by thermal stability engineering. To enhance the specific activity of efficient AlyG2, we strategically substituted residues with bulky steric hindrance proximal to the active pocket with glycine or alanine. This led to the generation of three promising positive mutants, with particular emphasis on the T91S mutant, exhibiting a 1.91-fold specific activity compared to the wild type. To mitigate the poor thermal stability of T91S, mutants with negative ΔΔG values in the thermal flexibility region were screened out. Notably, the S72Ya mutant not only displayed 17.96 % further increase in specific activity but also exhibited improved stability compared to T91S, manifesting as a remarkable 30.97 % increase in relative activity following a 1-hour incubation at 42 °C. Furthermore, enhanced kinetic stability was observed. To gain deeper insights into the mechanism underlying the enhanced thermostability of the S72Ya mutant, we conducted molecular dynamics simulations, principal component analysis (PCA), dynamic cross-correlation map (DCCM), and free energy landscape (FEL) analysis. The results unveiled a reduction in the flexibility of the surface loop, a stronger correlation dynamic and a narrower motion subspace in S72Ya system, along with the formation of more stable hydrogen bonds. Collectively, our findings suggest amino acids substitutions resulting in smaller side chains proximate to the active site can positively impact enzyme activity, while reducing the flexibility of surface loops emerges as a pivotal factor in conferring thermal stability. These insights offer valuable guidance and a framework for the engineering of other enzyme types.

Functional Characterization and Catalytic Activity Improvement of Borneol Acetyltransferase from Wurfbainia longiligularis

In plant secondary metabolite biosynthesis, acylation is a diverse physiological process, with BAHD acyltransferases playing an essential role. Borneol acetyltransferase (BAT) is an alcohol acetyltransferase, which catalyzes borneol and acetyl-CoA to synthesize bornyl acetate (BA). However, the enzymes involved in the biosynthesis of BA have so far only been characterized in Wurfbainia villosa, the studies on the WvBATs have only been conducted in vitro, and the catalytic activity was relatively low. In this research, three genes (WlBAT1, WlBAT2, and WlBAT3) have been identified to encode BATs that are capable of acetylating borneol to synthesize BA in vitro. We also determined that WlBAT1 has the highest catalytic efficiency for borneol-type substrates, including (+)-borneol, (-)-borneol, and isoborneol. Furthermore, we found that BATs could catalyze a wide range of substrate types in vitro, but in vivo, they exclusively catalyzed borneol-type substrates. Through molecular simulations and site-directed mutagenesis, it was revealed that residues D32, N36, H168, N297, N355, and H384 are crucial for the catalytic activity of WlBAT1, while the R382I-D385R double mutant of WlBAT1 exhibited an increasing acylation efficiency for borneol-type substrates in vitro and in vivo. These findings offer key genetic elements for the metabolic engineering of plants and synthetic biology to produce BA.

Combining Flexible Region Design and Automatic Design to Enhance the Thermal Stability and Catalytic Efficiency of Leucine Dehydrogenase

Leucine dehydrogenase (LeuDH, EC 1.4.1.9) can reversibly catalyze the oxidative deamination of l-leucine and some other specific α-amino acids to form the corresponding α-ketoacids. This reaction has great significance in the field of food additives and the pharmaceutical industry. The LeuDH from Exiguobacterium sibiricum (EsLeuDH) has high catalytic efficiency but limited thermal stability, hindering its widespread industrial application. In this study, a mutant N5F/I12L/A352Y of EsLeuDH (referred to as M2) was developed with enhanced thermal stability and catalytic activity through rational modification. The M2 mutant exhibits a half-life at 60 °C (t1/2(60 °C)) of 975.7 min and a specific activity of 69.6 U mg-1, which is 5.4 and 2.1 times higher than those of EsLeuDH, respectively. This research may facilitate the utilization of EsLeuDH at elevated temperatures, enhancing its potential for industrial applications. The findings offer a practical and efficient approach for optimizing LeuDH and other industrial enzymes.

Reshaping Phosphatase Substrate Preference for Controlled Biosynthesis Using a "Design-Build-Test-Learn" Framework

Biosynthesis is the application of enzymes in microbial cell factories and has emerged as a promising alternative to chemical synthesis. However, natural enzymes with limited catalytic performance often need to be engineered to meet specific needs through a time-consuming trial-and-error process. This study presents a quantum mechanics (QM)-incorporated design-build-test-learn (DBTL) framework to rationally design phosphatase BT4131, an enzyme with an ambiguous substrate spectrum involved in N-acetylglucosamine (GlcNAc) biosynthesis. First, mutant M1 (L129Q) is designed using force field-based methods, resulting in a 1.4-fold increase in substrate preference (kcat/Km) toward GlcNAc-6-phosphate (GlcNAc6P). QM calculations indicate that the shift in substrate preference is caused by a 13.59 kcal mol-1 reduction in activation energy. Furthermore, an iterative computer-aided design is conducted to stabilize the transition state. As a result, mutant M4 (I49Q/L129Q/G172L) with a 9.5-fold increase in kcat-GlcNAc6P/Km-GlcNAc6P and a 59% decrease in kcat-Glc6P/Km-Glc6P is highly desirable compared to the wild type in the GlcNAc-producing chassis. The GlcNAc titer increases to 217.3 g L-1 with a yield of 0.597 g (g glucose)-1 in a 50-L bioreactor, representing the highest reported level. Collectively, this DBTL framework provides an easy yet fascinating approach to the rational design of enzymes for industrially viable biocatalysts.

Enhancing Acid Resistance of Aspergillus niger Pectin Lyase through Surface Charge Design for Improved Application in Juice Clarification

Pectin lyases (PNLs) can enhance juice clarity and flavor by degrading pectin in highly esterified fruits, but their inadequate acid resistance leads to rapid activity loss in juice. This study aimed to improve the acid resistance of Aspergillus niger PNL pelA through surface charge design. A modification platform was established by fusing pelA with a protein tag and expressing the fusion enzyme in Escherichia coli. Four single-point mutants were identified to increase the surface charge using computational tools. Moreover, the combined mutant M6 (S514D/S538E) exhibited 99.8% residual activity at pH 3.0. The M6 gene was then integrated into the A. niger genome using a multigene integration system to obtain the recombinant PNL AM6. Notably, AM6 improved the light transmittance of orange juice to 45.3%, which was 8.39 times higher than that of pelA. In conclusion, AM6 demonstrated the best-reported acid resistance, making it a promising candidate for industrial juice clarification.

Computer-Aided Semi-Rational Design to Enhance the Activity of l-Sorbosone Dehydrogenase from Gluconobacter oxidans WSH-004

The titer of the microbial fermentation products can be increased by enzyme engineering. l-Sorbosone dehydrogenase (SNDH) is a key enzyme in the production of 2-keto-l-gulonic acid (2-KLG), which is the precursor of vitamin C. Enhancing the activity of SNDH may have a positive impact on 2-KLG production. In this study, a computer-aided semirational design of SNDH was conducted. Based on the analysis of SNDH's substrate pocket and multiple sequence alignment, three modification strategies were established: (1) expanding the entrance of SNDH's substrate pocket, (2) engineering the residues within the substrate pocket, and (3) enhancing the electron transfer of SNDH. Finally, mutants S453A, L460V, and E471D were obtained, whose specific activity was increased by 20, 100, and 10%, respectively. In addition, the ability of Gluconobacter oxidans WSH-004 to synthesize 2-KLG was improved by eliminating H2O2. This study provides mutant enzymes and metabolic engineering strategies for the microbial-fermentation-based production of 2-KLG.

Efficient Production of a Thermostable Mutant of Transglutaminase by Streptomyces mobaraensis

The transglutaminase (TGase) from Streptomyces mobaraensis is widely used to improve the texture of protein-based foods. However, wild-type TGase is not heat-resistant, which is unfavorable for its application. In this study, we successfully constructed a S. mobaraensis strain that can efficiently produce TGm2, a thermostable mutant of S. mobaraensis TGase. First, S. mobaraensis DSM40587 was subjected to atmospheric room temperature plasma mutagenesis, generating mutant smY2022 with a 12.2-fold increase in TGase activity. Then, based on the double-crossover recombination, we replaced the coding sequence of the TGase with that of TGm2 in smY2022, obtaining the strain smY2022-TGm2. The extracellular TGase activity of smY2022-TGm2 reached 61.7 U/mL, 147% higher than that of smY2022. Finally, the catalytic properties of TGm2 were characterized. The half-life time at 60 °C and specific activity of TGm2 reached 64 min and 71.15 U/mg, 35.6- and 2.9-fold higher than those of the wild-type TGase, respectively. As indicated by SDS-PAGE analysis, TGm2 exhibited demonstrably better protein cross-linking ability than the wild-type TGase at 70 °C, although both enzymes shared a similar ability at 40 °C. With improved enzyme production and thermal stability, smY2022-TGm2 could be a competitive strain for the industrial production of transglutaminase.

In silico design of multipoint mutants for enhanced performance of Thermomyces lanuginosus lipase for efficient biodiesel production

BACKGROUND

Biodiesel, an emerging sustainable and renewable clean energy, has garnered considerable attention as an alternative to fossil fuels. Although lipases are promising catalysts for biodiesel production, their efficiency in industrial-scale application still requires improvement.

RESULTS

In this study, a novel strategy for multi-site mutagenesis in the binding pocket was developed via FuncLib (for mutant enzyme design) and Rosetta Cartesian_ddg (for free energy calculation) to improve the reaction rate and yield of lipase-catalyzed biodiesel production. Thermomyces lanuginosus lipase (TLL) with high activity and thermostability was obtained using the Pichia pastoris expression system. The specific activities of the mutants M11 and M21 (each with 5 and 4 mutations) were 1.50- and 3.10-fold higher, respectively, than those of the wild-type (wt-TLL). Their corresponding melting temperature profiles increased by 10.53 and 6.01 °C, [Formula: see text] (the temperature at which the activity is reduced to 50% after 15 min incubation) increased from 60.88 to 68.46 °C and 66.30 °C, and the optimum temperatures shifted from 45 to 50 °C. After incubation in 60% methanol for 1 h, the mutants M11 and M21 retained more than 60% activity, and 45% higher activity than that of wt-TLL. Molecular dynamics simulations indicated that the increase in thermostability could be explained by reduced atomic fluctuation, and the improved catalytic properties were attributed to a reduced binding free energy and newly formed hydrophobic interaction. Yields of biodiesel production catalyzed by mutants M11 and M21 for 48 h at an elevated temperature (50 °C) were 94.03% and 98.56%, respectively, markedly higher than that of the wt-TLL (88.56%) at its optimal temperature (45 °C) by transesterification of soybean oil.

CONCLUSIONS

An integrating strategy was first adopted to realize the co-evolution of catalytic efficiency and thermostability of lipase. Two promising mutants M11 and M21 with excellent properties exhibited great potential for practical applications for in biodiesel production.

Improving the thermal stability and catalytic activity of ulvan lyase by the combination of FoldX and KnowVolution campaign

Thermal stability is one of the most important properties of ulvan lyases for their application in algae biomass degradation. The Knowledge gaining directed eVolution (KnowVolution) protein engineering strategy could be employed to improve thermostability of ulvan lyase with less screening effort. Herein, the unfolding free energies (ΔΔG) of the loop region were calculated using FoldX and four sites (D103, G104, T113, Q229) were selected for saturation mutagenesis, resulting in the identification of a favorable single-site mutant Q229M. Subsequently, iteration mutation was carried out with the mutant N57P (previously obtained by our group) to further enhance the performance of ulvan lyase. The results showed that the most beneficial variant N57P/Q229M exhibited a 1.67-fold and 2-fold increase in residual activity compared to the wild type after incubation at 40 °C and 50 °C for 1 h, respectively. In addition, the variant produced 1.06 mg/mL of reducing sugar in 2 h, which was almost four times as much as the wild type. Molecular dynamics simulations revealed that N57P/Q229M mutant enhanced the structural rigidity by augmenting intramolecular hydrogen bonds. Meanwhile, the shorter proton transmission distance between the general base of the enzyme and the substrate contributed to the glycosidic bond breakage. Our research showed that in silico saturation mutagenesis using position scan module in FoldX allowed for faster screening of mutants with improved thermal stability, and combining it with KnowVolution enabled a balanced effect of thermal stability and enzyme activity in protein engineering.

Rational Redesign of Chitosanase to Enhance Thermostability and Catalytic Activity to Produce Chitooligosaccharides with a Relatively High Degree of Polymerization

Chitooligosaccharides (hdpCOS) with a high degree of polymerization (hdp, DP 4-10) generally have greater biological activities than those of low-DP (ldp, DP 2-3) COS. Chitosanase from Bacillus amyloliquefaciens KCP2 (Csn46) can degrade chitosan to more hdpCOS at high temperature (70 °C), but low thermal stability at this temperature makes it unsuitable for industrial application; the wild-type enzyme can only produce COS (DP 2-4) at lower temperatures. Several thermostable mutants were obtained by modifying chitosanase using a comprehensive strategy based on a computer-aided mutant design. A combination of four beneficial single-point mutations (A129L/T175 V/K70T/D34G) to Csn46 was selected to obtain a markedly improved mutant, Mut4, with a half-life at 60 °C extended from 34.31 to 690.80 min, and the specific activity increased from 1671.73 to 3528.77 U/mg. Mut4 produced COS with DPs of 2-4 and 2-7 at 60 and 70 °C, respectively. Therefore, Mut4 has the potential to be applied to the industrial-scale preparation of hdpCOS with high biological activity.

Dynamics of loops surrounding the active site architecture in GH5_2 subfamily TfCel5A for cellulose degradation

BACKGROUND

Lignocellulose is the most abundant natural biomass resource for the production of biofuels and other chemicals. The efficient degradation of cellulose by cellulases is a critical step for the lignocellulose bioconversion. Understanding the structure-catalysis relationship is vital for rational design of more stable and highly active enzymes. Glycoside hydrolase (GH) family 5 is the largest and most functionally diverse group of cellulases, with a conserved TIM barrel structure. The important roles of the various loop regions of GH5 enzymes in catalysis, however, remain poorly understood.

RESULTS

In the present study, we investigated the relationship between the loops surrounding active site architecture and its catalytic efficiency, taking TfCel5A, an enzyme from GH5_2 subfamily of Thermobifida fusca, as an example. Large-scale computational simulations and site-directed mutagenesis experiments revealed that three loops (loop 8, 3, and 7) around active cleft played diverse roles in substrate binding, intermediate formation, and product release, respectively. The highly flexible and charged residue triad of loop 8 was responsible for capturing the ligand into the active cleft. Severe fluctuation of loop 3 led to the distortion of sugar conformation at the - 1 subsite. The wobble of loop 7 might facilitate product release, and the enzyme activity of the mutant Y361W in loop 7 was increased by approximately 40%.

CONCLUSION

This study unraveled the vital roles of loops in active site architecture and provided new insights into the catalytic mechanism of the GH5_2 cellulases.

Combined Computer-Aided Predictors to Improve the Thermostability of Nattokinase

Food-derived nattokinase has strong thrombolytic activity and few side effects. In the field of medicine, nattokinase has been developed as an adjuvant drug for the treatment of thrombosis, and nattokinase-rich beverages and health foods have also shown great potential in the field of food development. At present, the poor thermostability of nattokinase limits its industrial production and application. In this study, we used several thermostability-prediction algorithms to predict nattokinase from Bacillus mojavensis LY-06 (AprY), and screened two variants S33T and T174V with increased thermostability and fibrinolytic activity. The t1/2 of S33T and T174V were 8.87-fold and 2.51-fold those of the wild type AprY, respectively, and their enzyme activities were also increased (1.17-fold and 1.28-fold, respectively). Although the thermostability of N218L was increased by 2.7 times, the fibrinolytic activity of N218L was only 73.3% of that of wild type AprY. The multiple-point mutation results showed that S33T-N218L and S33T-T174V-N218L variants lost their activity, and the T174V-N218L variant did not show any significant change in catalytic performance, while S33T-T174V increased its thermostability and activity by 21.3% and 24.8%, respectively. Although the S33T-T174V variant did not show the additive effect of thermostability, it combined the excellent transient thermostability of S33T with the better thrombolytic activity of T174V. Bioinformatics analysis showed that the overall structure of S33T and T174V variants tended to be stable, while the structure of S33T-T174V variant was more flexible. Local structure analysis showed that the increased rigidity of the active center region (positions 64-75) and the key loop region (positions 129-130, 155-163, 187-192, 237-241, and 268-270) determined the increased thermostability of all variants. In addition, the enhanced flexibility of S33T-T174V variant in the Ca1 binding region (positions 1-4, 75-82) and the peripheral region of the catalytic pocket (positions 210-216) may account for the inability to superpose its thermostability. We explored the effective strategy to enhance the thermostability of nattokinase, and the resulting variants have potential industrial production and application.