Increasing protein stability: Importance of ΔCp and the denatured state

Increasing the conformational stability of proteins is an important goal for both basic research and industrial applications. In vitro selection has been used successfully to increase protein stability, but more often site‐directed mutagenesis is used to optimize the various forces that contribute to protein stability. In previous studies, we showed that improving electrostatic interactions on the protein surface and improving the β‐turn sequences were good general strategies for increasing protein stability, and used them to increase the stability of RNase Sa. By incorporating seven of these mutations in RNase Sa, we increased the stability by 5.3 kcal/mol. Adding one more mutation, D79F, gave a total increase in stability of 7.7 kcal/mol, and a melting temperature 28°C higher than the wild‐type enzyme. Surprisingly, the D79F mutation lowers the change in heat capacity for folding, ΔC_p, by 0.6 kcal/mol/K. This suggests that this mutation stabilizes structure in the denatured state ensemble. We made other mutants that give some insight into the structure present in the denatured state. Finally, the thermodynamics of folding of these stabilized variants of RNase Sa are compared with those observed for proteins from thermophiles.     

Biophysical characterization of mutants of Bacillus subtilis lipase evolved for thermostability: factors contributing to increased activity retention

Previously, Lipase A from Bacillus subtilis was subjected to in vitro directed evolution using iterative saturation mutagenesis, with randomization sites chosen on the basis of the highest B-factors available from the crystal structure of the wild-type (WT) enzyme. This provided mutants that, unlike WT enzyme, retained a large part of their activity after heating above 65 °C and cooling down. Here, we subjected the three best mutants along with the WT enzyme to biophysical and biochemical characterization. Combining thermal inactivation profiles, circular dichroism, X-ray structure analyses and NMR experiments revealed that mutations of surface amino acid residues counteract the tendency of Lipase A to undergo precipitation under thermal stress. Reduced precipitation of the unfolding intermediates rather than increased conformational stability of the evolved mutants seems to be responsible for the activity retention.     

Engineering deamidation-susceptible asparagines leads to improved stability to thermal cycling in a lipase

At high temperatures, protein stability is influenced by chemical alterations; most important among them is deamidation of asparagines. Deamidation kinetics of asparagines depends on the local sequence, solvent, pH, temperature, and the tertiary structure. Suitable replacement of deamidated asparagines could be a viable strategy to improve deamidation-mediated loss in protein properties, specifically protein thermostability. In this study, we have used nano RP-HPLC coupled ESI MS/MS approach to identify residues susceptible to deamidation in a lipase (6B) on heat treatment. Out of 15 asparagines and six glutamines in 6B, only five asparagines were susceptible to deamidation at temperatures higher than 75°C. These five positions were subjected to site saturation mutagenesis followed by activity screen to identify the most suitable substitutions. Only three of the five asparagines were found to be tolerant to substitutions. Best substitutions at these positions were combined into a mutant. The resultant lipase (mutC) has near identical secondary structure and improved thermal tolerance as compared to its parent. The triple mutant has shown almost two-fold higher residual activity compared to 6B after four cycles at 90°C. MutC has retained more than 50% activity even after incubation at 100°C. Engineering asparagines susceptible to deamidation would be a potential strategy to improve proteins to withstand very high temperatures.     

Origins of protein denatured state compactness and hydrophobic clustering in aqueous urea: inferences from nonpolar potentials of mean force

Free energies of pairwise hydrophobic association are simulated in aqueous solutions of urea at concentrations ranging from 0-8 M. Consistent with the expectation that hydrophobic interactions are weakened by urea, the association of relatively large nonpolar solutes is destabilized by urea. However, the association of two small methane-sized nonpolar solutes in water has the opposite tendency of being slightly strengthened by the addition of urea. Such size effects and the dependence of urea-induced stability changes on the configuration of nonpolar solutes are not predicted by solvent accessible surface area approaches based on energetic parameters derived from bulk-phase solubilities of model compounds. Thus, to understand hydrophobic interactions in proteins, it is not sufficient to rely solely on transfer experiment data that effectively characterize a single nonpolar solute in an aqueous environment but not the solvent-mediated interactions among two or more nonpolar solutes. We find that the m-values for the rate of change of two-methane association free energy with respect to urea concentration is a dramatically nonmonotonic function of the spatial separation between the two methanes, with a distance-dependent profile similar to the corresponding two-methane heat capacity of association in pure water. Our results rationalize the persistence of residual hydrophobic contacts in some proteins at high urea concentrations and explain why the heat capacity signature (DeltaC(P)) of a compact denatured state can be similar to DeltaC(P) values calculated by assuming an open random-coil-like unfolded state.     

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.     

Mutational analysis demonstrates that specific electrostatic interactions can play a key role in the denatured state ensemble of proteins

The nature of the denatured state ensemble has been controversial for decades owing, in large part, to the difficulty in characterizing the structure and energetics of denatured state interactions. There is increasing evidence for relatively non-specific hydrophobic clustering in the denatured states of some proteins but other types of interactions are much less well characterized. Here, we report the characterization of highly specific electrostatic interactions in the denatured state of a small alpha-beta protein, the N-terminal domain of the ribosomal protein L9 (NTL9). Mutation of Lys12 to Met has been shown to increase the stability of NTL9 significantly through the disruption of denatured state interactions. Here, we describe the analysis of the pH-dependent stability of 13 mutants designed to probe the nature of the Lys12 denatured state interaction. Lys12 is located in a lysine-rich region of the protein but analysis of a set of Lys to Met mutants shows that it plays a unique role in the denatured state. Analysis of mutants of all of the acidic residues in NTL9 shows that Lys12 forms a specific non-native electrostatic interaction with Asp8 in the denatured state ensemble. Thus the distribution of charge-charge interactions in the denatured state ensemble of NTL9 appears to be biased by few key interactions and is very different from that expected in a random coil. We propose that these interactions are not encoded by local sequence effects but rather reflect interactions among residues more distant in sequence. These results demonstrate that electrostatic as well as hydrophobic interactions can play an important role in the denatured state ensemble.     

Mutants provide evidence of the importance of glycosydic chains in the activation of lipase 1 from Candida rugosa

Sequence analysis of Candida rugosa lipase 1 (LIP1) predicts the presence of three N-linked glycosylation sites at asparagine 291, 314, 351. To investigate the relevance of sugar chains in the activation and stabilization of LIP1, we directed site mutagenesis to replace the above mentioned asparagine with glutamine residues. Comparison of the activity of mutants with that of the wild-type (wt) lipase indicates that both 314 and 351 Asn to Gln substitutions influence, although at a different extent, the enzyme activity both in hydrolysis and esterification reactions, but they do not alter the enzyme water activity profiles in organic solvents or temperature stability. Introduction of Gln to replace Asn351 is likely to disrupt a stabilizing interaction between the sugar chain and residues of the inner side of the lid in the enzyme active conformation. The effect of deglycosylation at position 314 is more difficult to explain and might suggest a more general role of the sugar moiety for the structural stability of lipase 1. Conversely, Asn291Gln substitution does not affect the lipolytic or the esterase activity of the mutant that behaves essentially as the wt enzyme. This observation supports the hypothesis that changes in activity of Asn314Gln and Asn351Gln mutants are specifically due to deglycosylation.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Reversal of coenzyme specificity of 2,3‐butanediol dehydrogenase from Saccharomyces cerevisae and in vivo functional analysis

Saccharomyces cerevisiae NAD(H)‐dependent 2,3‐butanediol dehydrogenase (Bdh1), a medium chain dehydrogenase/reductase is the main enzyme catalyzing the reduction of acetoin to 2,3‐butanediol. In this work we focused on altering the coenzyme specificity of Bdh1 from NAD(H) to NADP(H). Based on homology studies and the crystal structure of the NADP(H)‐dependent yeast alcohol dehydrogenase Adh6, three adjacent residues (Glu²²¹, Ile²²², and Ala²²³) were predicted to be involved in the coenzyme specificity of Bdh1 and were altered by site‐directed mutagenesis. Coenzyme reversal of Bdh1 was obtained with double Glu221Ser/Ile222Arg and triple Glu221Ser/Ile222Arg/Ala223Ser mutants. The performance of the triple mutant for NADPH was close to that of native Bdh1 for NADH. The three engineered mutants were able to restore the growth of a phosphoglucose isomerase deficient strain (pgi), which cannot grow on glucose unless an alternative NADPH oxidizing system is provided, thus demonstrating their in vivo functionality. These mutants are interesting tools to reduce the excess of acetoin produced by engineered brewing or wine yeasts overproducing glycerol. In addition, they represent promising tools for the manipulation of the NADP(H) metabolism and for the development of a powerful catalyst in biotransformations requiring NADPH regeneration. Biotechnol. Bioeng. 2009; 104: 381–389 © 2009 Wiley Periodicals, Inc.     

Three steps in the thermal unfolding of F-actin: an experimental evidence

The development of differential scanning calorimetry has resulted in an increased interest in studies of the unfolding process in proteins with the aim of identifying domains and interactions with ligands or other proteins. Several of these studies were done with actin and showed that the thermal unfolding of F-actin occurs in at least three steps; this was interpreted as the denaturation of independent domains. In the present work, we have followed the thermal unfolding of F-actin using differential scanning calorimetry (DSC), CD spectroscopy, and probe fluorescence. We found that the three steps revealed through DSC are not the denaturation of independent domains. These three steps are a change in the environment of cys 374 at 49.5 degrees C; a modification at the nucleotide-binding site at 55 degrees C; and the unfolding of the peptide chain at 64 degrees C. Previous interpretations of the thermograms of F-actin were thus erroneous. Since DSC is now widely used to study proteins, our experimental approach and conclusions may also be relevant in denaturation studies of proteins in general.     

An analysis approach to identify specific functional sites in orthologous proteins using sequence and structural information: Application to neuroserpin reveals regions that differentially regulate inhibitory activity

The analysis of sequence conservation is commonly used to predict functionally important sites in proteins. We have developed an approach that first identifies highly conserved sites in a set of orthologous sequences using a weighted substitution‐matrix‐based conservation score and then filters these conserved sites based on the pattern of conservation present in a wider alignment of sequences from the same family and structural information to identify surface‐exposed sites. This allows us to detect specific functional sites in the target protein and exclude regions that are likely to be generally important for the structure or function of the wider protein family. We applied our method to two members of the serpin family of serine protease inhibitors. We first confirmed that our method successfully detected the known heparin binding site in antithrombin while excluding residues known to be generally important in the serpin family. We next applied our sequence analysis approach to neuroserpin and used our results to guide site‐directed polyalanine mutagenesis experiments. The majority of the mutant neuroserpin proteins were found to fold correctly and could still form inhibitory complexes with tissue plasminogen activator (tPA). Kinetic analysis of tPA inhibition, however, revealed altered inhibitory kinetics in several of the mutant proteins, with some mutants showing decreased association with tPA and others showing more rapid dissociation of the covalent complex. Altogether, these results confirm that our sequence analysis approach is a useful tool that can be used to guide mutagenesis experiments for the detection of specific functional sites in proteins. Proteins 2015; 83:135–152. © 2014 Wiley Periodicals, Inc.