Knowledge-guided laboratory evolution of protein thermolability

A general temperature-guided language model to design proteins of enhanced stability and activity

Designing protein mutants with both high stability and activity is a critical yet challenging task in protein engineering. Here, we introduce PRIME, a deep learning model, which can suggest protein mutants with improved stability and activity without any prior experimental mutagenesis data for the specified protein. Leveraging temperature-aware language modeling, PRIME demonstrated superior predictive ability compared to current state-of-the-art models on the public mutagenesis dataset across 283 protein assays. Furthermore, we validated PRIME's predictions on five proteins, examining the impact of the top 30 to 45 single-site mutations on various protein properties, including thermal stability, antigen-antibody binding affinity, and the ability to polymerize nonnatural nucleic acid or resilience to extreme alkaline conditions. More than 30% of PRIME-recommended mutants exhibited superior performance compared to their premutation counterparts across all proteins and desired properties. We developed an efficient and effective method based on PRIME to rapidly obtain multisite mutants with enhanced activity and stability. Hence, PRIME demonstrates broad applicability in protein engineering.

Impact of single-residue mutations on protein thermal stability: The case of threonine 83 of BC2L-CN lectin

The thermal stability of trimeric lectin BC2L-CN was investigated and found to be considerably altered when mutating residue 83, originally a threonine, located at the fucose-binding loop. Mutants were analyzed using differential scanning calorimetry and isothermal microcalorimetry. Although most mutations decreased the affinity of the protein for oligosaccharide H type 1, six mutations increased the melting temperature (Tm) by >5 °C; one mutation, T83P, increased the Tm value by 18.2 °C(T83P, Tm = 96.3 °C). In molecular dynamic simulations, the investigated thermostable mutants, T83P, T83A, and T83S, had decreased fluctuations in the loop containing residue 83. In the T83S mutation, the side-chain hydroxyl group of serine formed a hydrogen bond with a nearby residue, suggesting that the restricted movement of the side-chain resulted in fewer fluctuations and enhanced thermal stability. Residue 83 is located at the interface and near the upstream end of the equivalent loop in a different protomer; therefore, fluctuations by this residue likely propagate throughout the loop. Our study of the dramatic change in thermal stability by a single amino acid mutation provides useful insights into the rational design of protein structures, especially the structures of oligomeric proteins.

Structure of an Alkaline Pectate Lyase and Rational Engineering with Improved Thermo-Alkaline Stability for Efficient Ramie Degumming

Alkaline pectate lyases have biotechnological applications in plant fiber processing, such as ramie degumming. Previously, we characterized an alkaline pectate lyase from Bacillus clausii S10, named BacPelA, which showed potential for enzymatic ramie degumming because of its high cleavage activity toward methylated pectins in alkaline conditions. However, BacPelA displayed poor thermo-alkaline stability. Here, we report the 1.78 Å resolution crystal structure of BacPelA in apo form. The enzyme has the characteristic right-handed β-helix fold of members of the polysaccharide lyase 1 family and shows overall structural similarity to them, but it displays some differences in the details of the secondary structure and Ca²⁺-binding site. On the basis of the structure, 10 sites located in flexible regions and showing high B-factor and positive ΔT_m values were selected for mutation, aiming to improve the thermo-alkaline stability of the enzyme. Following site-directed saturation mutagenesis and screening, mutants A238C, R150G, and R216H showed an increase in the T₅₀¹⁵ value at pH 10.0 of 3.0 °C, 6.5 °C, and 7.0 °C, respectively, compared with the wild-type enzyme, interestingly accompanied by a 24.5%, 46.6%, and 61.9% increase in activity. The combined mutant R150G/R216H/A238C showed an 8.5 °C increase in the T₅₀¹⁵ value at pH 10.0, and an 86.1% increase in the specific activity at 60 °C, with approximately doubled catalytic efficiency, compared with the wild-type enzyme. Moreover, this mutant retained 86.2% activity after incubation in ramie degumming conditions (4 h, 60 °C, pH 10.0), compared with only 3.4% for wild-type BacPelA. The combined mutant increased the weight loss of ramie fibers in degumming by 30.2% compared with wild-type BacPelA. This work provides a thermo-alkaline stable, highly active pectate lyase with great potential for application in the textile industry, and also illustrates an effective strategy for rational design and improvement of pectate lyases.

Recombination of Compatible Substitutions by 2GenReP and InSiReP

The CompassR rule enables to identify the beneficial substitutions, which can be recombined in directed evolution with gradually improving the enzymatic properties. However, the question of how to efficiently explore the protein sequence space when ten or more beneficial substitutions are identified has not yet been addressed. Two recombination strategies 2GenReP and InSiReP employing CompassR are systematically investigated to minimize experimental efforts and maximize possible improvements. Here we describe the details of the 2GenReP and InSiReP procedure with an example of recombining 15 substitutions and discuss some important practical issues that should be considered for the application of 2GenReP and InSiReP, such as placing the substitutions into subsets. The core part of the protocol (Step1 to Step5) is transferable to other enzymes and any recombination of potential substitutions.

Research Advances of d-allulose: An Overview of Physiological Functions, Enzymatic Biotransformation Technologies, and Production Processes

d-allulose has a significant application value as a sugar substitute, not only as a food ingredient and dietary supplement, but also with various physiological functions, such as improving insulin resistance, anti-obesity, and regulating glucolipid metabolism. Over the decades, the physiological functions of d-allulose and the corresponding mechanisms have been studied deeply, and this product has been applied to various foods to enhance food quality and prolong shelf life. In recent years, biotransformation technologies for the production of d-allulose using enzymatic approaches have gained more attention. However, there are few comprehensive reviews on this topic. This review focuses on the recent research advances of d-allulose, including (1) the physiological functions of d-allulose; (2) the major enzyme families used for the biotransformation of d-allulose and their microbial origins; (3) phylogenetic and structural characterization of d-allulose 3-epimerases, and the directed evolution methods for the enzymes; (4) heterologous expression of d-allulose ketose 3-epimerases and biotransformation techniques for d-allulose; and (5) production processes for biotransformation of d-allulose based on the characterized enzymes. Furthermore, the future trends on biosynthesis and applications of d-allulose in food and health industries are discussed and evaluated in this review.

Significantly improving the thermostability of a hyperthermophilic GH10 family xylanase XynAF1 by semi-rational design

Xylanases have a broad range of applications in industrial biotechnologies, which require the enzymes to resist the high-temperature environments. The majority of xylanases have maximum activity at moderate temperatures, which limited their potential applications in industries. In this study, a thermophilic GH10 family xylanase XynAF1 from the high-temperature composting strain Aspergillus fumigatus Z5 was characterized and engineered to further improve its thermostability. XynAF1 has the optimal reaction temperature of 90 °C. The crystal structure of XynAF1 was obtained by X-ray diffraction after heterologous expression, purification, and crystallization. The high-resolution X-ray crystallographic structure of the protein-product complex was obtained by soaking the apo-state crystal with xylotetraose. Structure analysis indicated that XynAF1 has a rigid skeleton, which helps to maintain the hyperthermophilic characteristic. The homologous structure analysis and the catalytic center mutant construction of XynAF1 indicated the conserved catalytic center contributed to the high optimum catalytic temperature. The amino acids in the surface of xylanase XynAF1 which might influence the enzyme thermostability were identified by the structure analysis. Combining the rational design with the saturation mutation at the high B-value regions, the integrative mutant XynAF1-AC with a 6-fold increase of thermostability was finally obtained. This study efficiently improved the thermostability of a GH10 family xylanase by semi-rational design, which provided a new biocatalyst for high-temperature biotechnological applications. KEY POINTS: • Obtained the crystal structure of GH10 family hyperthermophilic xylanase XynAF1. • Shed light on the understanding of the GH10 family xylanase thermophilic mechanism. • Constructed a 6-fold increased thermostability recombinant xylanase.

CompassR Yields Highly Organic-Solvent-Tolerant Enzymes through Recombination of Compatible Substitutions

The CompassR (computer-assisted recombination) rule enables, among beneficial substitutions, the identification of those that can be recombined in directed evolution. Herein, a recombination strategy is systematically investigated to minimize experimental efforts and maximize possible improvements. In total, 15 beneficial substitutions from Bacillus subtilis lipase A (BSLA), which improves resistance to the organic cosolvent 1,4-dioxane (DOX), were studied to compare two recombination strategies, the two-gene recombination process (2GenReP) and the in silico guided recombination process (InSiReP), employing CompassR. Remarkably, both strategies yielded a highly DOX-resistant variant, M4 (I12R/Y49R/E65H/N98R/K122E/L124K), with up to 14.6-fold improvement after screening of about 270 clones. M4 has a remarkably enhanced resistance in 60 % (v/v) acetone (6.0-fold), 30 % (v/v) ethanol (2.1-fold), and 60 % (v/v) methanol (2.4-fold) compared with wild-type BSLA. Molecular dynamics simulations revealed that attracting water molecules by charged surface substitutions is the main driver for increasing the DOX resistance of BSLA M4. Both strategies and obtained molecular knowledge can likely be used to improve the properties of other enzymes with a similar α/β-hydrolase fold.

Computer-Assisted Recombination (CompassR) Teaches us How to Recombine Beneficial Substitutions from Directed Evolution Campaigns

A main remaining challenge in protein engineering is how to recombine beneficial substitutions. Systematic recombination studies show that poorly performing variants are usually obtained after recombination of 3 to 4 beneficial substitutions. This limits researchers in exploiting nature's potential in generating better enzymes. The Computer-assisted Recombination (CompassR) strategy provides a selection guide for beneficial substitutions that can be recombined to gradually improve enzyme performance by analysis of the relative free energy of folding (ΔΔGfold ). The performance of CompassR was evaluated by analysis of 84 recombinants located on 13 positions of Bacillus subtilis lipase A. The finally obtained variant F17S/V54K/D64N/D91E had a 2.7-fold improved specific activity in 18.3 % (v/v) 1-butyl-3-methylimidazolium chloride ([BMIM][Cl]). In essence, the deducted CompassR rule allows recombination of beneficial substitutions in an iterative manner and empowers researchers to generate better enzymes in a time-efficient manner.

RankProt: A multi criteria-ranking platform to attain protein thermostabilizing mutations and its in vitro applications - Attribute based prediction method on the principles of Analytical Hierarchical Process

Attaining recombinant thermostable proteins is still a challenge for protein engineering. The complexity is the length of time and enormous efforts required to achieve the desired results. Present work proposes a novel and economic strategy of attaining protein thermostability by predicting site-specific mutations at the shortest possible time. The success of the approach can be attributed to Analytical Hierarchical Process and the outcome was a rationalized thermostable mutation(s) prediction tool- RankProt. Briefly the method involved ranking of 17 biophysical protein features as class predictors, derived from 127 pairs of thermostable and mesostable proteins. Among the 17 predictors, ionic interactions and main-chain to main-chain hydrogen bonds were the highest ranked features with eigen value of 0.091. The success of the tool was judged by multi-fold in silico validation tests and it achieved the prediction accuracy of 91% with AUC 0.927. Further, in vitro validation was carried out by predicting thermostabilizing mutations for mesostable Bacillus subtilis lipase and performing the predicted mutations by multi-site directed mutagenesis. The rationalized method was successful to render the lipase thermostable with optimum temperature stability and T_m increase by 20°C and 7°C respectively. Conclusively it can be said that it was the minimum number of mutations in comparison to the number of mutations incorporated to render Bacillus subtilis lipase thermostable, by directed evolution techniques. The present work shows that protein stabilizing mutations can be rationally designed by balancing the biophysical pleiotropy of proteins, in accordance to the selection pressure.

Rational and Semirational Protein Design

This mini review gives an overview over different design approaches and methodologies applied in rational and semirational enzyme engineering. The underlying principles for engineering novel activities, enantioselectivity, substrate specificity, stability, and pH optimum are summarized.

Catalytic stimulation by restrained active-site floppiness--the case of high density lipoprotein-bound serum paraoxonase-1

Despite the abundance of membrane-associated enzymes, the mechanism by which membrane binding stabilizes these enzymes and stimulates their catalysis remains largely unknown. Serum paraoxonase-1 (PON1) is a lipophilic lactonase whose stability and enzymatic activity are dramatically stimulated when associated with high-density lipoprotein (HDL) particles. Our mutational and structural analyses, combined with empirical valence bond simulations, reveal a network of hydrogen bonds that connect HDL binding residues with Asn168--a key catalytic residue residing >15Å from the HDL contacting interface. This network ensures precise alignment of N168, which, in turn, ligates PON1's catalytic calcium and aligns the lactone substrate for catalysis. HDL binding restrains the overall motion of the active site and particularly of N168, thus reducing the catalytic activation energy barrier. We demonstrate herein that disturbance of this network, even at its most far-reaching periphery, undermines PON1's activity. Membrane binding thus immobilizes long-range interactions via second- and third-shell residues that reduce the active site's floppiness and pre-organize the catalytic residues. Although this network is critical for efficient catalysis, as demonstrated here, unraveling these long-rage interaction networks is challenging, let alone their implementation in artificial enzyme design.

Enhancing thermostability and the structural characterization of Microbacterium saccharophilum K-1 β-fructofuranosidase

A β-fructofuranosidase from Microbacterium saccharophilum K-1 (formerly known as Arthrobacter sp. K-1) is useful for producing the sweetener lactosucrose (4(G)-β-D-galactosylsucrose). Thermostability of the β-fructofuranosidase was enhanced by random mutagenesis and saturation mutagenesis. Clones with enhanced thermostability included mutations at residues Thr47, Ser200, Phe447, Phe470, and Pro500. In the highest stability mutant, T47S/S200T/F447P/F470Y/P500S, the half-life at 60 °C was 182 min, 16.5-fold longer than the wild-type enzyme. A comparison of the crystal structures of the full-length wild-type enzyme and three mutants showed that various mechanisms appear to be involved in thermostability enhancement. In particular, the replacement of Phe447 with Val or Pro induced a conformational change in an adjacent residue His477, which results in the formation of a new hydrogen bond in the enzyme. Although the thermostabilization mechanisms of the five residue mutations were explicable on the basis of the crystal structures, it appears to be difficult to predict which amino acid residues should be modified to obtain thermostabilized enzymes.

Engineering lipase A from mesophilic Bacillus subtilis for activity at low temperatures

Loops or unordered regions of a protein are structurally dynamic and are strongly implicated in activity, stability and proteolytic susceptibility of proteins. Diminished activity of proteins at lower temperatures is considered to be due to compromised dynamics of the protein at lower temperatures. To evolve an active mesophilic lipase (Bacillus subtilis) at low temperatures, we subjected all the loop residues (n = 88) to site saturation mutagenesis (SSM). Based on a three-level screening protocol, we identified 14 substitutions, among 16,000 mutant population, which contributed to a substantial increase in activity at 5 °C. Based on the preliminary activity of recombinants at several temperatures, 5 substitutions among the 14 were found to be beneficial. A recombinant of these five mutations, named as 5CR, exhibited 7-fold higher catalytic efficiency than wild-type (WT) lipase at 10 °C. All the mutants, individually and in a recombinant (5CR), were characterized by substrate-binding parameters, melting temperatures and secondary structure. 5CR was similar to WT in substrate preferences and showed a significant improvement in activity at both lower and higher temperatures compared with the WT. To establish the contribution of mutations on the dynamics of the protein, we performed 100-ns molecular dynamics (MD) simulations on the WT and mutant lipase at 10 and 37 °C. The root mean square fluctuations (RMSFs) indeed showed that the mutations enhance the protein dynamics locally in the loop region having a catalytic residue, which may help in improved activities at lower temperatures.

Structural analysis on mutation residues and interfacial water molecules for human TIM disease understanding

Background

Human triosephosphate isomerase (HsTIM) deficiency is a genetic disease caused often by the pathogenic mutation E104D. This mutation, located at the side of an abnormally large cluster of water in the inter-subunit interface, reduces the thermostability of the enzyme. Why and how these water molecules are directly related to the excessive thermolability of the mutant have not been investigated in structural biology.

Results

This work compares the structure of the E104D mutant with its wild type counterparts. It is found that the water topology in the dimer interface of HsTIM is atypical, having a "wet-core-dry-rim" distribution with 16 water molecules tightly packed in a small deep region surrounded by 22 residues including GLU104. These water molecules are co-conserved with their surrounding residues in non-archaeal TIMs (dimers) but not conserved across archaeal TIMs (tetramers), indicating their importance in preserving the overall quaternary structure. As the structural permutation induced by the mutation is not significant, we hypothesize that the excessive thermolability of the E104D mutant is attributed to the easy propagation of atoms' flexibility from the surface into the core via the large cluster of water. It is indeed found that the B factor increment in the wet region is higher than other regions, and, more importantly, the B factor increment in the wet region is maintained in the deeply buried core. Molecular dynamics simulations revealed that for the mutant structure at normal temperature, a clear increase of the root-mean-square deviation is observed for the wet region contacting with the large cluster of interfacial water. Such increase is not observed for other interfacial regions or the whole protein. This clearly suggests that, in the E104D mutant, the large water cluster is responsible for the subunit interface flexibility and overall thermolability, and it ultimately leads to the deficiency of this enzyme.

Conclusions

Our study reveals that a large cluster of water buried in protein interfaces is fragile and high-maintenance, closely related to the structure, function and evolution of the whole protein.

Improvement of biocatalysts for industrial and environmental purposes by saturation mutagenesis

Laboratory evolution techniques are becoming increasingly widespread among protein engineers for the development of novel and designed biocatalysts. The palette of different approaches ranges from complete randomized strategies to rational and structure-guided mutagenesis, with a wide variety of costs, impacts, drawbacks and relevance to biotechnology. A technique that convincingly compromises the extremes of fully randomized vs. rational mutagenesis, with a high benefit/cost ratio, is saturation mutagenesis. Here we will present and discuss this approach in its many facets, also tackling the issue of randomization, statistical evaluation of library completeness and throughput efficiency of screening methods. Successful recent applications covering different classes of enzymes will be presented referring to the literature and to research lines pursued in our group. The focus is put on saturation mutagenesis as a tool for designing novel biocatalysts specifically relevant to production of fine chemicals for improving bulk enzymes for industry and engineering technical enzymes involved in treatment of waste, detoxification and production of clean energy from renewable sources.

A multi-factors rational design strategy for enhancing the thermostability of Escherichia coli AppA phytase

Despite recent advances in our understanding of the importance of protein surface properties for protein thermostability,there are seldom studies on multi-factors rational design strategy, so a more scientific, simple and effective rational strategy is urgent for protein engineering. Here, we first attempted to use a three-factors rational design strategy combining three common structural features, protein flexibility, protein surface, and salt bridges. Escherichia coli AppA phytase was used as a model enzyme to improve its thermostability. Moreover, the structure and enzyme features of the thermostable mutants designed by our strategy were analyzed roundly. For the single mutants, two (Q206E and Y311K), in five exhibited thermostable property with a higher success rate of prediction (40 %). For the multiple mutants, the themostable sites were combined with another site, I427L, we obtained by directed evolution, Q206E/I427L, Y311K/I427L, and Q206E/Y311K/I427L, all exhibited thermostable property. The Y311K/I427L doubled thermostability (61.7 %, and was compared to 30.97 % after being heated at 80 °C for 10 min) and catalytic efficiency (4.46 was compared to 2.37) improved more than the wild-type AppA phytase almost without hampering catalytic activity. These multi-factors of rational design strategy can be applied practically as a thermostabilization strategy instead of the conventional single-factor approach.

Probing protein stability and proteolytic resistance by loop scanning: a comprehensive mutational analysis

Improvement in protein thermostability was often found to be associated with increase in its proteolytic resistance as revealed by comparative studies of homologous proteins from extremophiles or mutational studies. Structural elements of protein responsible for this association are not firmly established although loops are implicated indirectly due to their structural role in protein stability. To get a better insight, a detailed study of protein wide mutants and their influence on stability and proteolytic resistance would be helpful. To generate such a data set, a model protein, Bacillus subtilis lipase was subjected to loop scanning site-saturation mutagenesis on 86 positions spanning all loops including termini. Upon screening of ~16,000 clones, 17 single mutants with improved thermostability were identified with increment in apparent melting temperature (Tm(app) ) by 1-6°C resulting in an increase in free energy of unfolding (ΔG(unf) ) by 0.04-1.16 kcal/mol. Proteolytic resistance of all single mutants upon incubation with nonspecific protease, Subtilisin A, was determined. Upon comparison, post-proteolysis residual activities as well as kinetics of proteolysis of mutants showed excellent correlation with ΔG(unf) , (r > 0.9), suggesting that proteolysis was strongly correlated with the global stability of this protein. This significant correlation in this set, with least possible sequence changes (single aa substitution), while covering >60% of protein surface strongly argues for the covariance of these two variables. Compared to studies from extremophiles, with large sequence heterogeneity, the observed correlation in such a narrow sequence space (ΔΔG(unf) = 1.57 kcal⁻¹) justifies the robustness of this relation.

Determination of ensemble-average pairwise root mean-square deviation from experimental B-factors

Root mean-square deviation (RMSD) after roto-translational least-squares fitting is a measure of global structural similarity of macromolecules used commonly. On the other hand, experimental x-ray B-factors are used frequently to study local structural heterogeneity and dynamics in macromolecules by providing direct information about root mean-square fluctuations (RMSF) that can also be calculated from molecular dynamics simulations. We provide a mathematical derivation showing that, given a set of conservative assumptions, a root mean-square ensemble-average of an all-against-all distribution of pairwise RMSD for a single molecular species, <RMSD(2)>(1/2), is directly related to average B-factors (<B>) and <RMSF(2)>(1/2). We show this relationship and explore its limits of validity on a heterogeneous ensemble of structures taken from molecular dynamics simulations of villin headpiece generated using distributed-computing techniques and the Folding@Home cluster. Our results provide a basis for quantifying global structural diversity of macromolecules in crystals directly from x-ray experiments, and we show this on a large set of structures taken from the Protein Data Bank. In particular, we show that the ensemble-average pairwise backbone RMSD for a microscopic ensemble underlying a typical protein x-ray structure is approximately 1.1 A, under the assumption that the principal contribution to experimental B-factors is conformational variability.