Protein Stability Engineering Research Articles

Correlation Analysis of Key Residue Sites between Computational-Aided Design Thermostability d-Amino Acid Oxidase and Ancestral Enzymes.

The d-amino acid oxidase (DAAO) from Rhodotorula taiwanensis has proven to have great potential for applications due to its excellent catalytic kinetic parameters. However, its poor thermal stability has limited its performance in biocatalysis. Herein, starting from the variant SHVG of RtwDAAO, this study employed a comprehensive computational design approach for protein stability engineering, resulting in positive substitutions at specific sites (A43S, T45M, C234L, E195Y). The generated variant combination, SHVG/SMLY, exhibited a significant synergistic effect, leading to an extension of the half-life and Tmapp. The ancestral sequence reconstruction revealed the conservation of the variant sites. The association of the variant sites with the highly stable ancestral enzyme was further explored. After determining the contribution of the variant sites to thermal stability, it was applied to other homologous sequences and validated. Molecular dynamics simulations indicated that the increased hydrophobicity of the variant SHVG/SMLY was a key factor for the increased stability, with strengthened intersubunit interactions playing an important role. In addition, the physical properties of the amino acids themselves were identified as crucial factors for thermal stability generality in homologous enzymes, which is important for the rapid acquisition of a series of stable enzymes.

Enhanced Thermostability of Candida Ketoreductase by Computation-Based Cross-Regional Combinatorial Mutagenesis

The flexible regions represented by loops have been the focus of protein stability engineering. However, residues located in specific rigid regions are also crucial for protein stability. The role of dynamic correlations between mutations in these two regions on protein stabilization remains inadequately understood. In this work, a ketoreductase CpKR originated from Candida parapsilosis was rationally engineered to improve its stability by comprehensively redesigning the flexible and rigid regions. Hot spots were successively identified and validated in the rigid regions with low B-factor or root-mean-square fluctuation values, as well as in relatively flexible areas of the protein. Based on the non-additive and cooperative mutational effects, a combinatorial variant (M2) with beneficial mutations in these two distinct regions showed a 4630-fold increased half-life at 40 °C and a 12 °C higher apparent melting temperature as compared with those of the wild type. The variant M2 also showed improved pH compatibility, better organic substance tolerance, and enhanced performance in asymmetric reduction of methyl 2-chloro-3-(4-methoxyphenyl)-3-oxo-propanoate (1a), the precursor of cardiovascular drug Diltiazem. The crystal structures and molecular dynamics simulations demonstrated that additional interaction networks in both the rigid and flexible regions modulated protein stability, and the dynamic correlation between these two regions mediated the observed synergistic combination. These results showed that both the flexible and rigid regions, especially the cross-correlated areas, are crucial for enhancing protein stability, and comprehensively redesigning these two regions is a powerful and attractive approach for acquiring stable enzymes.

General method to stabilize mesophilic proteins in hyperthermal water

SummaryThe stability of protein structures and biological functions at normal temperature is closely linked with the universal aqueous environment of organisms. Preserving bioactivities of proteins in hyperthermia water would expand their functional capabilities beyond those in native environments. However, only a limited number of proteins derived from hyperthermophiles are thermostable at elevated temperatures. Triggered by this, here we describe a general method to stabilize mesophilic proteins in hyperthermia water. The mesophilic proteins, protected by amphiphilic polymers with multiple binding sites, maintain their secondary and tertiary structures after incubation even in boiling water. This approach, outside the conventional environment for bioactivities of mesophilic proteins, provides a general strategy to dramatically increase the Tm (melting temperature) of mesophilic proteins without any changes to amino sequences of the native proteins. Current work offers a new insight with protein stability engineering for potential application, including vaccine storage and enzyme engineering.

Protein stability engineering insights revealed by domain-wide comprehensive mutagenesis

The accurate prediction of protein stability upon sequence mutation is an important but unsolved challenge in protein engineering. Large mutational datasets are required to train computational predictors, but traditional methods for collecting stability data are either low-throughput or measure protein stability indirectly. Here, we develop an automated method to generate thermodynamic stability data for nearly every single mutant in a small 56-residue protein. Analysis reveals that most single mutants have a neutral effect on stability, mutational sensitivity is largely governed by residue burial, and unexpectedly, hydrophobics are the best tolerated amino acid type. Correlating the output of various stability-prediction algorithms against our data shows that nearly all perform better on boundary and surface positions than for those in the core and are better at predicting large-to-small mutations than small-to-large ones. We show that the most stable variants in the single-mutant landscape are better identified using combinations of 2 prediction algorithms and including more algorithms can provide diminishing returns. In most cases, poor in silico predictions were tied to compositional differences between the data being analyzed and the datasets used to train the algorithm. Finally, we find that strategies to extract stabilities from high-throughput fitness data such as deep mutational scanning are promising and that data produced by these methods may be applicable toward training future stability-prediction tools.

DeepDDG: Predicting the Stability Change of Protein Point Mutations Using Neural Networks.

Accurately predicting changes in protein stability due to mutations is important for protein engineering and for understanding the functional consequences of missense mutations in proteins. We have developed DeepDDG, a neural network-based method, for use in the prediction of changes in the stability of proteins due to point mutations. The neural network was trained on more than 5700 manually curated experimental data points and was able to obtain a Pearson correlation coefficient of 0.48-0.56 for three independent test sets, which outperformed 11 other methods. Detailed analysis of the input features shows that the solvent accessible surface area of the mutated residue is the most important feature, which suggests that the buried hydrophobic area is the major determinant of protein stability. We expect this method to be useful for large-scale design and engineering of protein stability. The neural network is freely available to academic users at http://protein.org.cn/ddg.html .

RiSLnet: Rapid identification of smart mutant libraries using protein structure network. Application to thermal stability enhancement

A key point of protein stability engineering is to identify specific target residues whose mutations can stabilize the protein structure without negatively affecting the function or activity of the protein. Here, we propose a method called RiSLnet (Rapid identification of Smart mutant Library using residue network) to identify such residues by combining network analysis for protein residue interactions, identification of conserved residues, and evaluation of relative solvent accessibility. To validate its performance, the method was applied to four proteins, that is, T4 lysozyme, ribonuclease H, barnase, and cold shock protein B. Our method predicted beneficial mutations in thermal stability with ~62% average accuracy when the thermal stability of the mutants was compared with the ones in the Protherm database. It was further applied to lysine decarboxylase (CadA) to experimentally confirm its accuracy and effectiveness. RiSLnet identified mutations increasing the thermal stability of CadA with the accuracy of ~60% and significantly reduced the number of candidate residues (~99%) for mutation. Finally, combinatorial mutations designed by RiSLnet and in silico saturation mutagenesis yielded a thermally stable triple mutant with the half-life (T 1/2 ) of 114.9 min at 58°C, which is approximately twofold higher than that of the wild-type.

Deciphering the factors responsible for the stability of a GFP variant resistant to alkaline pH using molecular dynamics simulations

Charged amino acids having ionizable side chains play crucial roles in maintaining the solubility and stability of a protein. These charged amino acids are mostly exposed on protein surface and participate in electrostatic interactions with neighboring charged amino acids as well as with solvent. Therefore, the change in the solvent pH affects the protein stability in most cases. Previously, we reported a GFP variant, GFP14R having 14 surface lysines replaced with arginines, that showed enhanced stability under alkaline pH. Here, we analyzed the factors that contribute to the stability of the GFP14R under alkaline pH quantitatively using molecular dynamics simulations. Protonation state of the charged amino acids of GFP14R and control GFP under neutral pH and alkaline pH were modeled, and molecular dynamics simulations were performed. This comparative analysis revealed that the GFP14R with more arginine frequency on the surface maintained the stability under both pH conditions without much change in their salt-bridge interactions as well as the hydrogen bond interactions with solvent. On the other hand, these interactions were significantly reduced for the control GFP under alkaline pH due to the deprotonated lysine side chains. These results suggest that the advantageous property of arginine over lysine can be considered one of the parameter for the protein stability engineering under alkaline pH conditions.

Variations in the stability of NCR ene reductase by rational enzyme loop modulation

The engineering of protein stability is of major importance for the application of enzymes in a wide range of industrial applications. Here we study the determinants of the thermo- and solvent stability of the Zymomonas mobilis ene reductase NCR using a rational protein engineering approach based on analyses of structural and sequence data. We designed and created two loop mutants with the aim to increase their overall stability. They all retained catalytic activity but exhibited altered thermostability relative to the wild-type enzyme. The modulation of one specific loop segment near the active site of NCR showed an increased tolerance to organic solvents along with an enhanced thermostability.

Consistent mutational paths predict eukaryotic thermostability

BackgroundProteomes of thermophilic prokaryotes have been instrumental in structural biology and successfully exploited in biotechnology, however many proteins required for eukaryotic cell function are absent from bacteria or archaea. With Chaetomium thermophilum, Thielavia terrestris and Thielavia heterothallica three genome sequences of thermophilic eukaryotes have been published.ResultsStudying the genomes and proteomes of these thermophilic fungi, we found common strategies of thermal adaptation across the different kingdoms of Life, including amino acid biases and a reduced genome size. A phylogenetics-guided comparison of thermophilic proteomes with those of other, mesophilic Sordariomycetes revealed consistent amino acid substitutions associated to thermophily that were also present in an independent lineage of thermophilic fungi. The most consistent pattern is the substitution of lysine by arginine, which we could find in almost all lineages but has not been extensively used in protein stability engineering. By exploiting mutational paths towards the thermophiles, we could predict particular amino acid residues in individual proteins that contribute to thermostability and validated some of them experimentally. By determining the three-dimensional structure of an exemplar protein from C. thermophilum (Arx1), we could also characterise the molecular consequences of some of these mutations.ConclusionsThe comparative analysis of these three genomes not only enhances our understanding of the evolution of thermophily, but also provides new ways to engineer protein stability.

A Study on the Effect of Surface Lysine to Arginine Mutagenesis on Protein Stability and Structure Using Green Fluorescent Protein

Two positively charged basic amino acids, arginine and lysine, are mostly exposed to protein surface, and play important roles in protein stability by forming electrostatic interactions. In particular, the guanidinium group of arginine allows interactions in three possible directions, which enables arginine to form a larger number of electrostatic interactions compared to lysine. The higher pKa of the basic residue in arginine may also generate more stable ionic interactions than lysine. This paper reports an investigation whether the advantageous properties of arginine over lysine can be utilized to enhance protein stability. A variant of green fluorescent protein (GFP) was created by mutating the maximum possible number of lysine residues on the surface to arginines while retaining the activity. When the stability of the variant was examined under a range of denaturing conditions, the variant was relatively more stable compared to control GFP in the presence of chemical denaturants such as urea, alkaline pH and ionic detergents, but the thermal stability of the protein was not changed. The modeled structure of the variant indicated putative new salt bridges and hydrogen bond interactions that help improve the rigidity of the protein against different chemical denaturants. Structural analyses of the electrostatic interactions also confirmed that the geometric properties of the guanidinium group in arginine had such effects. On the other hand, the altered electrostatic interactions induced by the mutagenesis of surface lysines to arginines adversely affected protein folding, which decreased the productivity of the functional form of the variant. These results suggest that the surface lysine mutagenesis to arginines can be considered one of the parameters in protein stability engineering.

3DM: Systematic analysis of heterogeneous superfamily data to discover protein functionalities

Ten years of experience with molecular class-specific information systems (MCSIS) such as with the hand-curated G protein-coupled receptor database (GPCRDB) or the semiautomatically generated nuclear receptor database has made clear that a wide variety of questions can be answered when protein-related data from many different origins can be flexibly combined. MCSISes revolve around a multiple sequence alignment (MSA) that includes "all" available sequences from the entire superfamily, and it has been shown at many occasions that the quality of these alignments is the most crucial aspect of the MCSIS approach. We describe here a system called 3DM that can automatically build an entire MCSIS. 3DM bases the MSA on a multiple structure alignment, which implies that the availability of a large number of superfamily members with a known three-dimensional structure is a requirement for 3DM to succeed well. Thirteen MCSISes were constructed and placed on the Internet for examination. These systems have been instrumental in a large series of research projects related to enzyme activity or the understanding and engineering of specificity, protein stability engineering, DNA-diagnostics, drug design, and so forth.

Consensus-Based Engineering of Protein Stability: From Intrabodies to Thermostable Enzymes

Publisher Summary This chapter elaborates the consensus-based engineering of protein stability. Protein stability is an experimental metric of the information encoded in a protein sequence. This approach has allowed to engineer hyperstable immunoglobulin domains that can be expressed as intrabodies. Its application by a number of groups to a diverse set of proteins has demonstrated this procedure to be the most reliable strategy for the rational generation of stabilizing mutations. The ability to generate synthetic enzymatically active proteins from consensus sequences emphasizes the statistical nature of natural protein sequences as products of neutral drift and selective pressure. Stabilization is not the only objective for which the approach is useful, and others could include the selection of residues for nondisruptive destabilizing mutations to reduce the lifetime of proteins, or to generate proteins that can be used to screen for second-site revertants in evolutionary engineering strategies. It is found that the final steps of mutagenesis, characterization, and combination of successful residues are specific for each experimental system at hand and should be straightforward.

Stabilization of GroEL minichaperones by core and surface mutations

We report the crystal structures of two hexa-substituted mutants of a GroEL minichaperone that are more stable than wild-type by 7.0 and 6.1 kcal mol−1. Their structures imply that the increased stability results from multiple factors including improved hydrophobic packing, optimised hydrogen bonding and favourable structural rearrangements. It is commonly believed that protein core residues are immutable and generally optimized for energy, while on the contrary, surface residues are variable and hence unimportant for stability. But, it is now becoming clear that mutations of both core and surface residues can increase protein stability, and that protein cores are more flexible and thus more tolerant to mutation than expected. Sequence comparison of homologous proteins has provided a way to pinpoint the residues that contribute constructively to stability and to guide the engineering of protein stability. Stabilizing mutations identified by this approach are most frequently located at protein surfaces but with a few found in protein cores. In the latter case, local flexibility in the hydrophobic core is the key factor that allows the energetically favourable burial of larger hydrophobic side-chains without undue energetic penalties from steric clashes.