Thermostable Mutant Research Articles

Computational design of thermostable mutants for cephalosporin C acylase from Pseudomonas strain SE83

Computational protein design strategies can be used to increase enzyme stability without the need for high-throughput screening. In this report, computational methods were used to redesign cephalosporin C acylase from Pseudomonas strain SE83 to enhance its stability by repacking the hydrophobic core regions and reconstructing the protein-protein interactions in the segment interface regions. A nine-fold mutant with enhanced catalytic activity in the hydrolysis of cephalosporin C to 7-aminocephalosporanic acid, but with low stability, was used as a starting point. A computational enzyme design strategy was used to identify target regions to increase the protein melting temperature (Tm). Single point mutations Asn2βThr, Asn2βVal, Cys470βSer, Leu154βPhe, and Leu180βPhe in hydrophobic core regions, and Ala100αSer and Ala37βSer in segment-segment interface regions, increased the Tm by 4.7–19.7° C, while combining these confirmed single mutations increased the Tm by up to 20.5° C.

Enhancing the thermostability of fumarase C from Corynebacterium glutamicum via molecular modification

Fumarases have been successfully applied in industry for the production of l-malate. However, the industrialization of fumarases is limited by their low thermostability. In this study, the thermostability of fumarase C from Corynebacterium glutamicum was enhanced through directed evolution, simulated mutagenesis, site-directed mutagenesis and saturated mutagenesis. Mutant 2G (A411V) was initially constructed through directed evolution. Its half-life at 50 °C (t1/2, 50°C) increased from 1 min to 2.2 min, and the T5015 (temperature at which the activity of enzyme decreased by 50% in 15 min) increased from 44.8 °C to 47.2 °C. Besides, several different mutants were obtained by site-directed mutation. Among them, mutant 3G (A227V) showed significant improvement in thermostability with a 3.3-fold improvement of t1/2, 50°C and a 3.6 °C increase in T5015 compared to the wild-type enzyme. Then, 2/3G (A227V, A411V) was obtained by combining the mutant 2G with the mutant 3G, for which the t1/2, 50°C and T5015 increased to more than 768 min and 52.4 °C, respectively. Finally, site-saturated mutagenesis was employed on amino acid residues 175-Glu, 228-Gly, 297-Gly, 320-Lys and 464-Glu to maximize the thermostability of mutant 2/3G. The most thermostable mutant 175G with amino acid substitutions (A227V, A411V, E175K) was isolated. Its t1/2,50°C increased to more than 2700 min while that of wild-type enzyme was only 1 min and T5015 was 9.8 °C higher than the wild-type enzyme. The thermostable mutated enzymes generated without affecting the activity in this study would be an attractive candidate for industrial applications.

Directed -in vitro- evolution of Precambrian and extant Rubiscos

Rubisco is an ancient, catalytically conserved yet slow enzyme, which plays a central role in the biosphere’s carbon cycle. The design of Rubiscos to increase agricultural productivity has hitherto relied on the use of in vivo selection systems, precluding the exploration of biochemical traits that are not wired to cell survival. We present a directed -in vitro- evolution platform that extracts the enzyme from its biological context to provide a new avenue for Rubisco engineering. Precambrian and extant form II Rubiscos were subjected to an ensemble of directed evolution strategies aimed at improving thermostability. The most recent ancestor of proteobacteria -dating back 2.4 billion years- was uniquely tolerant to mutagenic loading. Adaptive evolution, focused evolution and genetic drift revealed a panel of thermostable mutants, some deviating from the characteristic trade-offs in CO2-fixing speed and specificity. Our findings provide a novel approach for identifying Rubisco variants with improved catalytic evolution potential.

Machine Learning Methods for Rapid Prediction of Thermostabilizing Mutants of G‐protein‐coupled receptors in Detergents

G protein‐coupled receptors (GPCRs) are highly dynamic and often denature when extracted in detergents. Deriving thermostable mutants has been a successful strategy to stabilize GPCRs in detergents, but this process is experimentally tedious. We have developed a computational method that scans the entire receptor for alanine or leucine mutations and calculates the stability score for mutation of each position in the GPCR structural model. The method involves generating homology models of the receptors of varying accuracies and an ensemble of conformations by sampling the rigid body degrees of freedom of transmembrane helices. Then, an all‐atom force field function is used to calculate the enthalpy gain, known as the “stability score” upon mutation of every residue, in these receptor structures, to alanine. Inclusion of conformational sampling to account for structural perturbations improves the thermostability predictions compared to using a single structure. We have validated the method against experimentally measured thermostability data for single mutants of the β1‐adrenergic receptor (β1AR), adenosine A2A receptor (A2AR) and neurotensin receptor 1 (NTSR1). We will present some of our recent results on blind predictions of thermostabilizing mutations on NTSR1.Using the computed energy components in our dataset, we have recently tested several machine learning algorithms for rapid and reliable prediction of thermostable and non‐thermostable mutants. We used Matlab to build a binary classification model (thermostable vs. non‐thermostable) trained on data from an ensemble of three proteins (25% of dataset used for holdout cross‐validation), and tested against a fourth protein omitted from the training set. Initial dataset was heavily imbalanced to the “non‐thermostable” class. This was balanced by random replication of the undersampled “thermostable” class. Models trained by decision tree algorithms performed best in our study compared to k‐nearest neighbor (k‐NN) algorithms. Our precision and recall statistics for the best performing ensemble dataset were as follows: Boosted Trees – 62% precision, 34% recall, Simple Trees – 30% precision, 90% recall, Complex Trees – 36% precision, 73% recall, Medium k‐NN – 28% precision, 34% recall.These results are promising for predicting true thermostable mutants, and we are working to improve our methods for reducing the number of false positive predictions. This requires testing different training algorithms such as Support Vector Machines and other k‐NN and Decision Tree models, and potentially ensemble training models combining multiples of these algorithms. Our method is the first step toward a computational method for rapid prediction of thermostable mutants of GPCRs. These tools are critical for reliable and quick prediction of thermostabilizing mutations of GPCRs for drug discovery.Support or Funding InformationFunding for this work was from NIH‐RO1GM097261 to N.V.This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

Engineering a thermostable highly active glucose 6-phosphate dehydrogenase and its application to hydrogen production in vitro.

Glucose 6-phosphate dehydrogenase (G6PDH) is one of the most important dehydrogenases responsible for generating reduced NADPH for anabolism and is also the rate-limiting enzyme in the Entner-Doudoroff pathway. For in vitro biocatalysis, G6PDH must possess both high activity and good thermostability due to requirements of efficient use and low expense of biocatalyst. Here, we used directed evolution to improve thermostability of the highly active G6PDH from Zymomonas mobilis. Four generations of random mutagenesis and Petri-dish-based double-layer screening evolved the thermolabile wild-type enzyme to the thermostable mutant Mut 4-1, which showed a more than 124-fold increase in half-life time (t1/2) at 60°C, a 3.4°C increase in melting temperature (T m ), and a 5°C increase in optimal temperature (Topt), without compromising the specific activity. In addition, the thermostable mutant was conducted to generate hydrogen from maltodextrin via in vitro synthetic biosystems (ivSB), gaining a more than 8-fold improvement of productivity rate with 76% of theoretical yield at 60°C. Thus, the engineered G6PDH has been shown to effectively regenerate NADPH at high temperatures and will be applicable for NAD(P)H regeneration in numerous in vitro biocatalysis applications.

Enhancing the Thermostability of Rhizomucor miehei Lipase with a Limited Screening Library by Rational-Design Point Mutations and Disulfide Bonds.

Rhizomucor miehei lipase (RML), as a kind of eukaryotic protein catalyst, plays an important role in the food, organic chemical, and biofuel industries. However, RML retains its catalytic activity below 50°C, which limits its industrial applications at higher temperatures. Soluble expression of this eukaryotic protein in Escherichia coli not only helps to screen for thermostable mutants quickly but also provides the opportunity to develop rapid and effective ways to enhance the thermal stability of eukaryotic proteins. Therefore, in this study, RML was engineered using multiple computational design methods, followed by filtration via conservation analysis and functional region assessment. We successfully obtained a limited screening library (only 36 candidates) to validate thermostable single point mutants, among which 24 of the candidates showed higher thermostability and 13 point mutations resulted in an apparent melting temperature ([Formula: see text]) of at least 1°C higher. Furthermore, both of the two disulfide bonds predicted from four rational-design algorithms were further introduced and found to stabilize RML. The most stable mutant, with T18K/T22I/E230I/S56C-N63C/V189C-D238C mutations, exhibited a 14.3°C-higher [Formula: see text] and a 12.5-fold increase in half-life at 70°C. The catalytic efficiency of the engineered lipase was 39% higher than that of the wild type. The results demonstrate that rationally designed point mutations and disulfide bonds can effectively reduce the number of screened clones to enhance the thermostability of RML.IMPORTANCER. miehei lipase, whose structure is well established, can be widely applied in diverse chemical processes. Soluble expression of R. miehei lipase in E. coli provides an opportunity to explore efficient methods for enhancing eukaryotic protein thermostability. This study highlights a strategy that combines computational algorithms to predict single point mutations and disulfide bonds in RML without losing catalytic activity. Through this strategy, an RML variant with greatly enhanced thermostability was obtained. This study provides a competitive alternative for wild-type RML in practical applications and further a rapid and effective strategy for thermostability engineering.

Параллельные вычисления в разработке термостабильных мутантов липаз

Параллельные вычисления в разработке термостабильных мутантов липаз

Erratum to “Expression, Purification and Crystallization of Thermostable Mutant of Cutinase Est1 from <i>Thermobifida alba</i>.” [Advances in Bioscience and Biotechnology 9 (2018), 215-223

The original online version of this article (Kitadokoro, K., Matsui, S., Osokoshi, R., Nakata, K. and Kamitani, S. (2018) Expression, Purification and Crystallization of Thermostable Mutant of Cutinase Est1 from Thermobifida alba. Advances in Bioscience and Biotechnology, 9, 215-223. https://doi.org/10.4236/abb.2018.95015) unfortunately contains some mistakes. The author wishes to correct the errors.

Engineering stability in NADPH oxidases: A common strategy for enzyme production.

NADPH oxidases (NOXs) are membrane enzymes whose sole function is the generation of reactive oxygen species. Humans have seven NOX isoenzymes that feature distinct functions in immune response and cell signaling but share the same catalytic core comprising a FAD-binding dehydrogenase domain and a heme-binding transmembrane domain. We previously described a mutation that stabilizes the dehydrogenase domain of a prokaryotic homolog of human NOX5. The thermostable mutant exhibited a large 19 °C increase in the apparent melting temperature (app Tm) and a much tighter binding of the FAD cofactor, which allowed the crystallization and structure determination of the domain holo-form. Here, we analyze the transferability of this mutation onto prokaryotic and eukaryotic full-length NOX enzymes. We found that the mutation exerts a significative stabilizing effect on the full-length NOX5 from both Cylindrospermum stagnale (app Tm increase of 8 °C) and Homo sapiens (app ΔTm of 2 °C). Enhanced thermal stability resulted in more homogeneous preparations of the bacterial NOX5 with less aggregation problems. Moreover, we also found that the mutation increases the overall expression of recombinant human NOX4 and NOX5 in mammalian cells. Such a 2-5-fold increase is mainly due to the lowered cell toxicity, which leads to higher biomasses. Because of the high sequence identity of the catalytic core within this family of enzymes, this strategy can be a general tool to boost the production of all NOXs.

Improved thermostability of a metagenomic glucose-tolerant β-glycosidase based on its X-ray crystal structure.

MeBglD2, a metagenomic β-glycosidase, is stimulated by various saccharides, including D-glucose, D-xylose, and maltose, and it promotes the enzymatic saccharification of plant biomass. To improve the thermostability of MeBglD2, its X-ray crystal structure was analyzed, and the amino acid residues responsible for its thermostability were identified using the structural information. Mutations in His8, Asn59, and Gly295 improved the thermostability of MeBglD2, and the combination of these mutations resulted in the highest thermostability. Compared with wild-type MeBglD2, thermostable MeBglD2 mutants promoted plant biomass saccharification using Trichoderma reesei cellulase. In addition to thermostability, the thermostable mutants exhibited higher tolerance to ethanol, dimethyl sulfoxide, and copper ions, indicating that the MeBglD2 mutants generated in this study were improved in their tolerance to not only high temperature but also to organic solvents and metal ions.

Characterization of a thermostable mutant of Agaricus brasiliensis laccase created by phylogeny-based design

Laccases are enzymes that oxidize various aromatic compounds, and therefore they have attracted much attention from the standpoints of medical and industrial applications. We previously isolated the cDNA that codes for a laccase isozyme (Lac2a) from the medicinal mushroom Agaricus brasiliensis (Matsumoto-Akanuma etal., Int. J. Med. Mushrooms, 16, 375-393, 2014). In this study, we first attempted heterologous expression of the wild-type laccase using a Pichia pastoris secretory expression system. However, the trial was unsuccessful most likely because the enzyme was too unstable and degraded immediately after production. Therefore, we improved the stability of the laccase by using a phylogeny-based design method. We created a mutant laccase in which sixteen original residues were replaced with those found in the phylogenetically inferred ancestral sequence. The resulting mutant protein was successfully produced using the P.pastoris secretory expression system and then purified. The designed laccase showed catalytic properties similar to those of other fungal laccases. Moreover, the laccase is highly thermally stable at acidic and neutral pH and is also stable at alkaline pH at moderate temperatures. We expect that the laccase will serve as a useful tool for enzymatic polymerization of di-phenolic compounds.

Improved genome recovery and integrated cell-size analyses of individual uncultured microbial cells and viral particles

Microbial single-cell genomics can be used to provide insights into the metabolic potential, interactions, and evolution of uncultured microorganisms. Here we present WGA-X, a method based on multiple displacement amplification of DNA that utilizes a thermostable mutant of the phi29 polymerase. WGA-X enhances genome recovery from individual microbial cells and viral particles while maintaining ease of use and scalability. The greatest improvements are observed when amplifying high G+C content templates, such as those belonging to the predominant bacteria in agricultural soils. By integrating WGA-X with calibrated index-cell sorting and high-throughput genomic sequencing, we are able to analyze genomic sequences and cell sizes of hundreds of individual, uncultured bacteria, archaea, protists, and viral particles, obtained directly from marine and soil samples, in a single experiment. This approach may find diverse applications in microbiology and in biomedical and forensic studies of humans and other multicellular organisms.

Thermophilic Enzyme or Mesophilic Enzyme with Enhanced Thermostability: Can We Draw a Line?

Aminoglycoside nucleotidyltransferase 4' (ANT) is a homodimeric enzyme that modifies the C4'-OH site of aminoglycoside antibiotics by nucleotidylation. A few single- and double-residue mutants of this enzyme (T130K, D80Y, and D80Y/T130K) from Bacillus stearothermophilus show increased thermostability. This article investigates how such residue replacements, which are distant from the active site and monomer-monomer interface, result in various changes of the thermostability of the enzyme. In this work, we show that the thermodynamic properties of enzyme-ligand complexes and protein dynamics may be indicators of a thermophilic behavior. Our data suggests that one of the single-site mutants of ANT, D80Y, may be a thermophilic protein and the other thermostable mutant, T130K, is actually a more heat-stable variant of the mesophilic wild type (WT) with a higher Tm. Our data also suggest that T130K and D80Y adopt different global dynamics strategies to achieve different levels of thermostability enhancement and that the differences between the properties of the species can be described in terms of global dynamics rather than in terms of specific structural features. Thermophilicity of the D80Y comes at the cost of less favorable thermodynamic parameters for ligand binding relative to WT. On the other hand, the T130K species exhibits the same affinity to ligands and the same thermodynamic parameters of complex formation as the WT enzyme. These observations suggest that a quantitative characterization of ligand binding and protein dynamics can be used to differentiate thermophilic proteins from their simply more heat-stable mesophilic counterparts.

FireProt: web server for automated design of thermostable proteins.

There is a continuous interest in increasing proteins stability to enhance their usability in numerous biomedical and biotechnological applications. A number of in silico tools for the prediction of the effect of mutations on protein stability have been developed recently. However, only single-point mutations with a small effect on protein stability are typically predicted with the existing tools and have to be followed by laborious protein expression, purification, and characterization. Here, we present FireProt, a web server for the automated design of multiple-point thermostable mutant proteins that combines structural and evolutionary information in its calculation core. FireProt utilizes sixteen tools and three protein engineering strategies for making reliable protein designs. The server is complemented with interactive, easy-to-use interface that allows users to directly analyze and optionally modify designed thermostable mutants. FireProt is freely available at http://loschmidt.chemi.muni.cz/fireprot.

Improving the reversibility of thermal denaturation and catalytic efficiency of Bacillus licheniformis α-amylase through stabilizing a long loop in domain B.

The reversibility of thermal denaturation and catalytic efficiency of Bacillus licheniformis α-amylase were improved through site-directed mutagenesis. By using multiple sequence alignment and PoPMuSiC algorithm, Ser187 and Asn188, which located within a long loop in Domain B of Bacillus licheniformis α-amylase, were selected for mutation. In addition, Ala269, which is adjacent to Ser187 and Asn188, was also investigated. Seven mutants carrying the mutations S187D, N188T, N188S, A269K, A269K/S187D, S187D/N188T, and A269K/S187D/N188T were generated and characterized. The most thermostable mutant, A269K/S187D/N188T, exhibited a 9-fold improvement in half-life at 95°C and pH 5.5, compared with that of the wild-type enzyme. Mutant A269K/S187D/N188T also exhibited improved catalytic efficiency. The catalytic efficiency of mutant A269K/S187D/N188T reached 5.87×103±0.17 g·L-1·s-1 at pH 5.5, which is 1.84-fold larger than the corresponding value determined for the wild-type enzyme. Furthermore, the structure analysis showed that immobilization of the loop containing Ser187 and Asn188 plays a significant role in developing the properties of Bacillus licheniformis α-amylase.

The role of local residue environmental changes in thermostable mutants of the GH11 xylanase from Bacillus subtilis

A thermostable variant of the mesophilic xylanase A from Bacillus subtilis (BsXynA-G3_4x) contains the four mutations Gln7His, Gly13Arg, Ser22Pro, and Ser179Cys. The crystal structure of the BsXynA-G3_4x has been solved, and the local environments around each of these positions investigated by molecular dynamics (MD) simulations at 328K and 348K. The structural and MD simulation results were correlated with thermodynamic data of the wild-type enzyme, the 4 single mutants and the BsXynA-G3_4x. This analysis suggests that the overall stabilizing effect is entropic, and is consistent with solvation of charged residues and reduction of main-chain flexibility. Furthermore, increased protein-protein hydrogen bonding and hydrophobic interactions also contribute to stabilize the BsXynA-G3_4x. The study revealed that a combination of several factors is responsible for increased thermostability of the BsXynA-G3_4x; (i) introduction of backbone rigidity in regions of high flexibility, (ii) solvation effects and (iii) hydrophobic contacts.

Modification of the peroxygenative:peroxidative activity ratio in the unspecific peroxygenase from Agrocybe aegerita by structure-guided evolution.

Unspecific peroxygenase (UPO) is a heme-thiolate peroxidase capable of performing with high-selectivity C-H oxyfunctionalizations of great interest in organic synthesis through its peroxygenative activity. However, the convergence of such activity with an unwanted peroxidative activity encumbers practical applications. In this study, we have modified the peroxygenative:peroxidative activity ratio (P:p ratio) of UPO from Agrocybe aegerita by structure-guided evolution. Several flexible loops (Glu1-Pro35, Gly103-Asp131, Ser226-Gly243, Gln254-Thr276 and Ty293-Arg327) were selected on the basis on their B-factors and ΔΔG values. The full ensemble of segments (43% of UPO sequence) was subjected to focused evolution by the Mutagenic Organized Recombination Process by Homologous IN vivo Grouping (MORPHING) method in Saccharomyces cerevisiae. Five independent mutant libraries were screened in terms of P:p ratio and thermostability. We identified several variants that harbored substitutions at positions 120 and 320 with a strong enhancement in the P:p ratio albeit at the cost of stability. The most thermostable mutant of this process (S226G with an increased T50 of 2°C) was subjected to further combinatorial saturation mutagenesis on Thr120 and Thr320 yielding a collection of variants with modified P:p ratio and recovered stability. Our results seem to indicate the coexistence of several oxidation sites for peroxidative and peroxygenative activities in UPO.

Site‐directed mutagenesis studies of hydrophobic residues in the lid region of T1 lipase

T1 lipase is a potential biocatalyst for industrial application since it is highly thermostable and displays optimal activity at 60–75°C. Structural analysis of T1 lipase shows that its lid region undergoes a spatial displacement along with a distinct secondary structure reorganization upon activation. To study structure/function of this atypical lid, we performed site‐directed mutagenesis on the hydrophobic residues in the lid region. These residues were mutated to hydrophilic ones and biochemical properties of mutants were investigated. Results showed that F181 might be an important residue for enzyme‐substrate binding. Mutants A186S and A190S had 35–50% increase in catalytic efficiencies compared to wild‐type T1, without compromising their functions at high temperatures. In general, mutagenesis did not cause large changes to chain‐length preference in T1 lipase. Mutants A186S and V187N were inactive towards long‐chain pNP esters (p‐Nitrophenyl stearate) and V187N showed lower activity towards long‐chain triacylglycerols than wt T1, which makes it a potential catalyst in dairy industry. Thermostability of mutants were affected at different extent due to the influence of hydrophobic contact between the lid and the protein core. These findings not only shed light on the lipase structure/function relationship but also lay the framework for further engineering to gain more potent, stable, and selective lipases.Practical applications: A thermostable T1 lipase owns an atypical lid and hydrophobic residues in the lid influence its catalytic properties. Here, we screened mutants A186S and A190S with 35–50% increase in catalytic efficiency, which have industrial application in the modification of fats and oils. Moreover, mutant V187N showed lower activity towards long‐chain TAGs than wt T1. This mutant can be a potential biocatalyst in dairy industry which promote generation of short‐chain fatty acids from milk fats to increase flavors in dairy products.The lid region of thermostable T1 lipase undergoes a spatial displacement along with a distinct secondary structure reorganization upon activation (PDB ID: 2DSN, 2W22 [16,17]). Mutagenesis on hydrophobic residues in the lid affect catalytic properties of T1 lipase.

Effects of point mutations on the thermostability of B. subtilis lipase: investigating nonadditivity.

Molecular level understanding of mutational effects on stability and activity of enzymes is challenging particularly when several point mutations are incorporated during the directed evolution experiments. In our earlier study, we have suggested the lack of consistency in the effect of point mutations incorporated during the initial generations of directed evolution experiments, towards conformational stabilization of B. subtilis lipase mutants of later generations. Here, we report that the cumulative point mutations incorporated in mutants 2M (with two point mutations) to 6M (with six point mutations) possibly do not retain their original stabilizing nature in the most thermostable 12M mutant (with 12 point mutations). We have carried out MD simulations using structures incorporating reversal of different sets of point mutations to assess their effect on the conformational stability and activity of 12M. Our analysis has revealed that reversal of certain point mutations in 12M had little effect on its conformational stability, suggesting that these mutations were probably inconsequential towards the thermostability of the 12M mutant. Interestingly these mutations involved evolutionarily conserved residues. On the other hand, some of the other point mutations incorporated in nonconserved regions, appeared to contribute significantly towards the conformational stability and/or activity of 12M. Based on the analysis of dynamics of in silico mutants generated using the consensus sequence, we identified experimentally verifiable residue positions to further increase the conformational stability and activity of the 12M mutant.

Engineering increased thermostability in the GH-10 endo-1,4-β-xylanase from Thermoascus aurantiacus CBMAI 756

The GH10 endo-xylanase from Thermoascus aurantiacus CBMAI 756 (XynA) is industrially attractive due to its considerable thermostability and high specific activity. Considering the possibility of a further improvement in thermostability, eleven mutants were created in the present study via site-directed mutagenesis using XynA as a template. XynA and its mutants were successfully overexpressed in Escherichia coli Rosetta-gami DE3 and purified, exhibiting maximum xylanolytic activity at pH 5 and 65°C. Three of the eleven mutants, Q158R, H209N, and N257D, demonstrated increased thermostability relative to the wild type at 70°C and 75°C.Q158R and N257D were stable in the pH range 5.0–10.0, while WT and H209N were stable from pH 8–10. CD analysis demonstrated that the WT and the three mutant enzymes were expressed in a folded form. H209N was the most thermostable mutant, showing a Tm of 71.3°C. Molecular dynamics modeling analyses suggest that the increase in H209N thermostability may beattributed to a higher number of short helices and salt bridges, which displayed a positive charge in the catalytic core, stabilizing its tertiary structure.