Average assignment method for predicting the stability of protein mutants

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and it will be helpful for designing stable mutants. In this work, we have analyzed the stability of protein mutants using three different data sets of 1791, 1396, and 2204 mutants, respectively, for thermal stability (DeltaTm), free energy change due to thermal (DeltaDeltaG), and denaturant denaturations (DeltaDeltaGH2O), obtained from the ProTherm database. We have classified the mutants into 380 possible substitutions and assigned the stability of each mutant using the information obtained with similar type of mutations. We observed that this assignment could distinguish the stabilizing and destabilizing mutants to an accuracy of 70-80% at different measures of stability. Further, we have classified the mutants based on secondary structure and solvent accessibility (ASA) and observed that the classification significantly improved the accuracy of prediction. The classification of mutants based on helix, strand, and coil distinguished the stabilizing/destabilizing mutants at an average accuracy of 82% and the correlation is 0.56; information about the location of residues at the interior, partially buried, and surface regions of a protein correctly identified the stabilizing/destabilizing residues at an average accuracy of 81% and the correlation is 0.59. The nine subclassifications based on three secondary structures and solvent accessibilities improved the accuracy of assigning stabilizing/destabilizing mutants to an accuracy of 84-89% for the three data sets. Further, the present method is able to predict the free energy change (DeltaDeltaG) upon mutations within a deviation of 0.64 kcal/mol. We suggest that this method could be used for predicting the stability of protein mutants.     

Importance of mutant position in Ramachandran plot for predicting protein stability of surface mutations

Understanding the mechanisms by which mutations affect protein stability is one of the most important problems in molecular biology. In this work, we analyzed the relationship between changes in protein stability caused by surface mutations and changes in 49 physicochemical, energetic, and conformational properties of amino acid residues. We found that the hydration entropy was the major contributor to the stability of surface mutations in helical segments; other properties responsible for size and volume of molecule also correlated significantly with stability. Classification of coil mutations based on their locations in the (phi-psi) map improved the correlation significantly, demonstrating the existence of a relationship between stability and strain energy, which indicates that the role of strain energy is very important for the stability of surface mutations. We observed that the inclusion of sequence and structural information raised the correlation, indicating the influence of surrounding residues on the stability of surface mutations. Further, we examined the previously reported "inverse relationship" between stability and hydrophobicity, and observed that the inverse hydrophobic effect was generally applicable only to coil mutations. The present study leads to a simple method for predicting protein stability changes caused by amino acid substitutions, which will be useful for protein engineering in designing novel proteins with increased stability and altered function.     

Revisiting “reverse hydrophobic effect”: Applicable only to coil mutations at the surface

In a seminal paper, Pakula and Sauer (Nature, 1990, 344, 363–364) demonstrated that the increase in side‐chain hydrophobicity has a reverse relationship with protein stability. We have addressed this problem with several examples of mutants that span at different locations in protein structure based on secondary structure and solvent accessibility. We confirmed that the stability change upon single coil mutation at exposed region is reversely correlated with hydrophobicity with a single exception. In addition, we found the existence of such relationship in partially buried coil mutants. The stability of exposed helical mutants is governed by conformational properties. In buried and partially buried helical and strand mutants properties reflecting hydrophobicity have direct relationship with stability, whereas an opposite relationship was obtained with entropy and flexibility. The structural analysis of partially buried/exposed mutants showed that the surrounding residues are important for the stability change upon mutation. These results provide insights to understand the general behavior for the stability of proteins upon amino acid substitutions. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 591–599, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Role of structural and sequence information in the prediction of protein stability changes: comparison between buried and partially buried mutations.

Predicting mutation-induced changes in protein stability is one of the greatest challenges in molecular biology. In this work, we analyzed the correlation between stability changes caused by buried and partially buried mutations and changes in 48 physicochemical, energetic and conformational properties. We found that properties reflecting hydrophobicity strongly correlated with stability of buried mutations, and there was a direct relation between the property values and the number of carbon atoms. Classification of mutations based on their location within helix, strand, turn or coil segments improved the correlation of mutations with stability. Buried mutations within beta-strand segments correlated better than did those in alpha-helical segments, suggesting stronger hydrophobicity of the beta-strands. The stability changes caused by partially buried mutations in ordered structures (helix, strand and turn) correlated most strongly and were mainly governed by hydrophobicity. Due to the disordered nature of coils, the mechanism underlying their stability differed from that of the other secondary structures: the stability changes due to mutations within the coil were mainly influenced by the effects of entropy. Further classification of mutations within coils, based on their hydrogen-bond forming capability, led to much stronger correlations. Hydrophobicity was the major factor in determining the stability of buried mutations, whereas hydrogen bonds, other polar interactions and hydrophobic interactions were all important determinants of the stability of partially buried mutations. Information about local sequence and structural effects were more important for the prediction of stability changes caused by partially buried mutations than for buried mutations; they strengthened correlations by an average of 27% among all data sets.     

Prediction of protein mutant stability using classification and regression tool

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and the solving of which would help for designing stable mutants. In this work, we have analyzed the stability of protein mutants using two different datasets of 1396 and 2204 mutants obtained from ProTherm database, respectively for free energy change due to thermal (DeltaDeltaG) and denaturant denaturations (DeltaDeltaG(H(2)O)). We have used a set of 48 physical, chemical energetic and conformational properties of amino acid residues and computed the difference of amino acid properties for each mutant in both sets of data. These differences in amino acid properties have been related to protein stability (DeltaDeltaG and DeltaDeltaG(H(2)O)) and are used to train with classification and regression tool for predicting the stability of protein mutants. Further, we have tested the method with 4 fold, 5 fold and 10 fold cross validation procedures. We found that the physical properties, shape and flexibility are important determinants of protein stability. The classification of mutants based on secondary structure (helix, strand, turn and coil) and solvent accessibility (buried, partially buried, partially exposed and exposed) distinguished the stabilizing/destabilizing mutants at an average accuracy of 81% and 80%, respectively for DeltaDeltaG and DeltaDeltaG(H(2)O). The correlation between the experimental and predicted stability change is 0.61 for DeltaDeltaG and 0.44 for DeltaDeltaG(H(2)O). Further, the free energy change due to the replacement of amino acid residue has been predicted within an average error of 1.08 kcal/mol and 1.37 kcal/mol for thermal and chemical denaturation, respectively. The relative importance of secondary structure and solvent accessibility, and the influence of the dataset on prediction of protein mutant stability have been discussed.     

Strain in protein structures as viewed through nonrotameric side chains: I. their position and interaction.

We studied the relative spatial positioning of nonrotameric side chains with atypical and strained dihedral angles in well-refined protein tertiary structures. The analysis was confined to buried protein cores, which are less error prone to side-chain positioning. More than half of the proteins with two or more nonrotameric residues displayed clusters of two or more (and up to five) nonrotameric residues. The clusters exhibited lower average crystallographic temperature factors compared with isolated nonrotameric residues. Nonrotameric clusters showed significantly tighter packing than corresponding rotameric clusters and had distinct residue compositions that did not correlate with amino acid characteristics such as size, hydrophobicity, turn preference, and the like. Such nonrotameric residue biases would suggest that spatially concentrated strain in protein folds would be minimized by lowered vibrational energy. Furthermore, nonrotameric residues avoided helices and strands and mostly preferred coil regions. If they were in the helical conformation, then they preferred to be within N-terminal segments. Proteins 1999;37:30-43.     

Role of medium- and long-range interactions to the stability of the mutants of T4 lysozyme

Inter-residue interactions play an important role to the folding and stability of protein molecules. In this work, we analyze the role of medium- and long-range interactions to the stability of T4 lysozyme mutants. We found that, in buried mutations, the increase in long-range contacts upon mutations destabilizes the protein, whereas, in surface mutations, the increase in long-range contacts increases the stability, indicating the importance of surrounding polar residues to the stability of surface mutations. Further, the increase in medium-range contacts decreases the stability of buried and surface mutations and a direct relationship is observed between the increase of medium-range contacts and increase in stability for partially buried/exposed mutations. Moreover, the relationship between amino acid properties and stability of T4 lysozyme mutants at positions Ile3, Phe53, and Leu99 showed that the effect of medium- and long-range contacts is less for buried mutations and the inter-residue contacts have significant correlation with the stability of partially buried mutations.     

Turn scanning by site-directed mutagenesis: application to the protein folding problem using the intestinal fatty acid binding protein

We have systematically mutated residues located in turns between beta-strands of the intestinal fatty acid binding protein (IFABP), and a glycine in a half turn, to valine and have examined the stability, refolding rate constants and ligand dissociation constants for each mutant protein. IFABP is an almost all beta-sheet protein exhibiting a topology comprised of two five-stranded sheets surrounding a large cavity into which the fatty acid ligand binds. A glycine residue is located in seven of the eight turns between the antiparallel beta-strands and another in a half turn of a strand connecting the front and back sheets. Mutations in any of the three turns connecting the last four C-terminal strands slow the folding and decrease stability with the mutation between the last two strands slowing folding dramatically. These data suggest that interactions between the last four C-terminal strands are highly cooperative, perhaps triggered by an initial hydrophobic collapse. We suggest that this trigger is collapse of the highly hydrophobic cluster of amino acids in the D and E strands, a region previously shown to also affect the last stage of the folding process (Kim et al., 1997). Changing the glycine in the strand between the front and back sheets also results in a unstable, slow folding protein perhaps disrupting the D-E strand interactions. For most of the other turn mutations there was no apparent correlation between stability and refolding rate constants. In some turns, the interaction between strands, rather than the turn type, appears to be critical for folding while in others, turn formation itself appears to be a rate limiting step. Although there is no simple correlation between turn formation and folding kinetics, we propose that turn scanning by mutagenesis will be a useful tool for issues related to protein folding.     

A study of the preferred environment of amino acid residues in globular proteins

In the native folded conformation of a globular protein, amino acid residues distant along the polypeptide chain come together to form the compact structure. This spatial structure is such that most of the polar residues are on the surface and have contact with the solvent medium and the nonpolar residues buried in the interior which have contact with similar nonpolar side chains. This cooperativity and mutual interaction among the randomly aligned amino acid residues suggest that each type of residue may prefer to have a specific environment. To gain more insight into this aspect of residue-residue cooperativity, a detailed analysis of the preferred environment associated with each of the 20 different amino acid residues in a number of protein crystals has been carried out. The variation of nonpolar nature computed for different sizes of spheres shows that the spatial region between radii of 6 and 8 Å is more favored for hydrophobic interactions and indicates that the influence of each residue over the surrounding medium extends predominantly up to a distance of 8 Å. The analysis of the surrounding amino acid residues associated with each type of residue shows that there is a definite tendency for each type of residue to have association with specific residues. The variation in environment is found even within the polar group as well as in the nonpolar group of residues. The surrounding residues associated with isoleucine, leucine, and valine are purely nonpolar. Proline, a nonpolar residue, is often surrounded by polar residues. The surrounding nonpolar nature of the tryptophan and tyrosine residues implies that even a single polar atom in a nonpolar side chain is sufficient to reduce their hydrophobic environment. There exists a high degree of mutual residue-residue cooperativity between the pairs glutamic acid-lysine, methionine-arginine, asparagine-tryptophan, and glutamine-proline, and the mutual residue-residue noncooperativity is high for the pairs methionine-aspartic acid, cysteine-glutamic acid, histidine-glutamine, and leucine-asparagine. The formation of secondary and tertiary structures is discussed in terms of the preferred environment and mutual cooperativity among various types of amino acid residues.     

Context-dependence of amino acid residue pairing in antiparallel beta-sheets.

In an effort to understand the driving forces behind antiparallel beta-sheet assembly, we have investigated the mutational tolerance of four pairs of residues in CspA, the major cold shock protein of E. coli. Two buried pairs and two exposed pairs of neighboring amino acids were separately randomized and the corresponding effects on protein stability were assessed using a protein expression screen. The thermal denaturation of a subset of the recovered proteins was measured by circular dichroism spectroscopy in order to determine the range of stabilities sampled by the expressed mutants. As anticipated, buried sites are substantially less tolerant of substitutions than exposed sites with more than half of the exposed residue combinations giving rise to stably folded proteins. The two exposed residue pairs, however, display different degrees of tolerance to substitution and accept different residue pair combinations. Except for the prohibition of proline from interior strand positions, no obvious correlations of mutant stability with any single parameter such as beta-sheet propensity or hydrophobicity can be detected. Mutant combinations recovered in both orientations (e.g. XY and YX) at a given exposed pair site often show markedly different stabilities, indicating that the local environment plays a substantial role in modulating the pairing preferences of residues in beta-sheets.     

Suggestions for "safe" residue substitutions in site-directed mutagenesis

The conserved topological structure observed in various molecular families such as globins or cytochromes c allows structural equivalencing of residues in every homologous structure and defines in a coherent way a global alignment in each sequence family. A search was performed for equivalent residue pairs in various topological families that were buried in protein cores or exposed at the protein surface and that had mutated but maintained similar unmutated environments. Amino acid residues with atoms in contact with the mutated residue pairs defined the environment. Matrices of preferred amino acid exchanges were then constructed and preferred or avoided amino acid substitutions deduced. Given the conserved atomic neighborhoods, such natural in vivo substitutions are subject to similar constrains as point mutations performed in site-directed mutagenesis experiments. The exchange matrices should provide guidelines for "safe" amino acid substitutions least likely to disturb the protein structure, either locally or in its overall folding pathway, and most likely to allow probing the structural and functional significance of the substituted site.     

Thermal stabilities of mutant Escherichia coli tryptophan synthase alpha subunits.

Random chemical mutagenesis, in vitro, of the 5' portion of the Escherichia coli trpA gene has yielded 66 mutant alpha subunits containing single amino acid substitutions at 49 different residue sites within the first 121 residues of the protein; this portion of the alpha subunit contains four of the eight alpha helices and three of the eight beta strands in the protein. Sixty-two of the subunits were examined for their heat stabilities by sensitivity to enzymatic inactivation (52 degrees C for 20 min) in crude extracts and by differential scanning calorimetry (DSC) with 29 purified proteins. The enzymatic activities of mutant alpha subunits that contained amino acid substitutions within the alpha and beta secondary structures were more heat labile than the wild-type alpha subunit. Alterations only in three regions, at or immediately C-terminal to the first three beta strands, were stability neutral or stability enhancing with respect to enzymatic inactivation. Enzymatic thermal inactivation appears to be correlated with the relative accessibility of the substituted residues; stability-neutral mutations are found at accessible residual sites, stability-enhancing mutations at buried sites. DSC analyses showed a similar pattern of stabilization/destabilization as indicated by inactivation studies. Tm differences from the wild-type alpha subunit varied +/- 7.6 degrees C. Eighteen mutant proteins containing alterations in helical and sheet structures had Tm's significantly lower (-1.6 to -7.5 degrees C) than the wild-type Tm (59.5 degrees C). In contrast, 6 mutant alpha subunits with alterations in the regions following beta strands 1 and 3 had increased Tm's (+1.4 to +7.6 degrees C). Because of incomplete thermal reversibilities for many of the mutant alpha subunits, most likely due to identifiable aggregated forms in the unfolded state, reliable differences in thermodynamic stability parameters are not possible. The availability of this group of mutant alpha subunits which clearly contain structural alterations should prove useful in defining the roles of certain residues or sequences in the unfolding/folding pathway for this protein when examined by urea/guaninidine denaturation kinetic analysis.     

Computation of non-covalent interactions in TNF proteins and interleukins

The roles played by the non-covalent interactions have been investigated for a set of six TNF proteins and nine Interleukins. The stabilizing residues have been identified by a consensus approach using the concepts of available surface area, medium and long-range interactions and conservation of amino acid residues. The cation-pi interactions have been computed based on a geometric approach such as distance and energy criteria. We identified an average of 1 energetically significant cation-pi interactions in every 94 residues in TNF proteins and 1 in every 62 residues in Interleukins. In TNF proteins, the cationic groups Lys preferred to be in helix while Arg preferred to be in strand regions while in Interleukins the Arg residues preferred to be in helix and Lys preferred to be in strand regions. From the available surface area calculations, we found that, almost all the cation and pi residues in TNF proteins and Interleukins were either in buried or partially buried regions and none of them in the exposed regions. Medium and long-range interactions were predominant in both TNF proteins and Interleukins. It was observed that the percentage of stabilizing centers were more in TNF proteins as compared to the Interleukins, while the percentage of conserved residues were more in Interleukins than in TNF proteins. In the stabilizing residues Lys was observed to be a stabilizing residue in both TNF proteins and Interleukins. Among the aromatic group, Phe was seen to be a stabilizing residue in both TNF and Interleukins. We suggest that this study on the computation of cation-pi interactions in TNF proteins and Interleukins would be very helpful in further understanding the structure, stability and functional similarity of these proteins.     

Possible role of two phenylalanine residues in the active site of human cytidine deaminase

Site-directed mutagenesis on human cytidine deaminase (CDA) was employed to mutate specifically two highly conserved phenylalanine residues, F36 and F137, to tryptophan; at the same time, the unique tryptophan residue present in the sequence at position 113 was mutated to phenylalanine. These double mutations were performed in order to have for each protein a single tryptophan signal for fluorescence studies relative to position 36 or 137. The mutant enzymes thus obtained, W113F, F36W/W113F and F137W/W113F, showed by circular dicroism and thermal stability an overall structure not greatly affected by the mutations. The titration of Trp residues by N-bromosuccinimide (NBS) suggested that residue W113 of the wild-type CDA and W36 of mutant F36W/W113F are buried in the tertiary structure of the enzyme, whereas the residue W137 of mutant F137W/W113F is located near the surface of the molecule. Kinetic experiments and equilibrium experiments with FZEB showed that the residue W113 seems not to be part of the active site of the enzyme whereas the Phe/Trp substitution in F36W/W113F and F137W/W113F mutant enzymes had a negative effect on substrate binding and catalysis, suggesting that F137 and F36 of the wild-type CDA are involved in a stabilizing interaction between ligand and enzyme.     

Dramatic stabilization of an SH3 domain by a single substitution: roles of the folded and unfolded states

The N-terminal SH3 domain of the Drosophila drk protein (drkN SH3) exists in equilibrium between folded and unfolded states under non-denaturing buffer conditions. In order to examine the origins of this instability, we have made mutations in the domain and characterized the thermodynamics and kinetics of folding. Results of substitutions of negatively charged residues to neutral amino acid residues suggest that the large electrostatic potential of the domain does not play a dominant role in the instability of the domain. Sequence alignment of a large number of SH3 domains reveals that the drkN SH3 domain has a threonine (T22) at a position corresponding to an otherwise highly conserved glycine residue in the diverging beta-turn connecting the beta3 and beta4 strands. Mutation of T22 to glycine results in significant stabilization of the drkN SH3 domain by 2.5 kcal/mole. To further characterize the basis for the stabilization of the T22 mutant relative to wild-type, we made additional mutant proteins with substitutions of residue T22. A strong correlation is seen between protein stability or folding rate and propensity for native beta-turn structure at this position. Correlation of folding rates with AGADIR predictions of non-native helical structure in the diverging turn region, along with our previous NMR evidence for non-native structure in this region of the unfolded state of the drkN SH3 domain, suggests that the free energy of the unfolded state also plays a role in stability. This result highlights the importance of both folded and unfolded states for understanding protein stability.     

Analyses of the Sequence and Structural Properties Corresponding to Pentapeptide and Large Palindromes in Proteins

The analyses of 3967 representative proteins selected from the Protein Data Bank revealed the presence of 2803 pentapeptide and large palindrome sequences with known secondary structure conformation. These represent 2014 unique palindrome sequences. 60% palindromes are not associated with any regular secondary structure and 28% are in helix conformation, 11% in strand conformation and 1% in the coil conformation. The average solvent accessibility values are in the range between 0–155.28 Ų suggesting that the palindromes in proteins can be either buried, exposed to the solvent or share an intermittent property. The number of residue neighborhood contacts defined by interactions ≤ 3.2 Ǻ is in the range between 0–29 residues. Palindromes of the same length in helix, strand and coil conformation are associated with different amino acid residue preferences at the individual positions. Nearly, 20% palindromes interact with catalytic/active site residues, ligand or metal ions in proteins and may therefore be important for function in the corresponding protein. The average hydrophobicity values for the pentapeptide and large palindromes range between -4.3 to +4.32 and the number of palindromes is almost equally distributed between the negative and positive hydrophobicity values. The palindromes represent 107 different protein families and the hydrolases, transferases, oxidoreductases and lyases contain relatively large number of palindromes.     

Molecular dissection,tissue localization and Ca2+ binding of the ryanodine receptor of Caenorhabditis elegans

We investigated the possible role of residues at the Ccap position in an alpha-helix on protein stability. A set of 431 protein alpha-helices containing a C'-Gly from the Protein Data Bank (PDB) was analyzed, and the normalized frequencies for finding particular residues at the Ccap position, the average fraction of buried surface area, and the hydrogen bonding patterns of the Ccap residue side-chain were calculated. We found that on average the Ccap position is 70% buried and noted a significant correlation (R=0.8) between the relative burial of this residue and its hydrophobicity as defined by the Gibbs energy of transfer from octanol or cyclohexane to water. Ccap residues with polar side-chains are commonly involved in hydrogen bonding. The hydrogen bonding pattern is such that, the longer side-chains of Glu, Gln, Arg, Lys, His form hydrogen bonds with residues distal (>+/-4) in sequence, while the shorter side-chains of Asp, Asn, Ser, Thr exhibit hydrogen bonds with residues close in sequence (<+/-4), mainly involving backbone atoms. Experimentally we determined the thermodynamic propensities of residues at the Ccap position using the protein ubiquitin as a model system. We observed a large variation in the stability of the ubiquitin variants depending on the nature of the Ccap residue. Furthermore, the measured changes in stability of the ubiquitin variants correlate with the hydrophobicity of the Ccap residue. The experimental results, together with the statistical analysis of protein structures from the PDB, indicate that the key hydrophobic capping interactions between a helical residue (C3 or C4) and a residue outside the helix (C", C3' or C4') are frequently enhanced by the hydrophobic interactions with Ccap residues.     

Single tryptophanyl substitutions affect the structure of apomyoglobin

Mammalian myoglobins contain two tryptophanyl residues at the invariant positions A-5 (W7) and A-12 (W14) in the N-terminal region (A helix) of the protein molecule. To determine the contribution of each tryptophanyl residue to the structure and stability of myoglobin, recombinant proteins with single indole residue, i.e., W7 or W14, were obtained by site-directed mutagenesis. The mutant proteins, expressed in Escherichia coli, were found correctly folded, the far ultraviolet circular dichroism of both mutants as well as the Soret absorption being superimposed to that of wild type protein. The removal of the prosthetic group from mutant proteins determined a loss of helical content much larger than that observed in the case of wild type myoglobin. These results suggest that tryptophanyl residues can play a crucial role on globin folding and structure.     

The role of strand 1 of the C β-sheet in the structure and function of α1-antitrypsin

Serpins inhibit cognate serine proteases involved in a number of important processes including blood coagulation and inflammation. Consequently, loss of serpin function or stability results in a number of disease states. Many of the naturally occurring mutations leading to disease are located within strand 1 of the C beta-sheet of the serpin. To ascertain the structural and functional importance of each residue in this strand, which constitutes the so-called distal hinge of the reactive center loop of the serpin, an alanine scanning study was carried out on recombinant alpha(1)-antitrypsin Pittsburgh mutant (P1 = Arg). Mutation of the P10' position had no effect on its inhibitory properties towards thrombin. Mutations to residues P7' and P9' caused these serpins to have an increased tendency to act as substrates rather than inhibitors, while mutations at P6' and P8' positions caused the serpin to behave almost entirely as a substrate. Mutations at the P6' and P8' residues of the C beta-sheet, which are buried in the hydrophobic core in the native structure, caused the serpin to become highly unstable and polymerize much more readily. Thus, P6' and P8' mutants of alpha(1)-antitrypsin had melting temperatures 14 degrees lower than wild-type alpha(1)-antitrypsin. These results indicate the importance of maintaining the anchoring of the distal hinge to both the inhibitory mechanism and stability of serpins, the inhibitory mechanism being particularly sensitive to any perturbations in this region. The results of this study allow more informed analysis of the effects of mutations found at these positions in disease-associated serpin variants.     

Hydrogen-exchange stabilities of RNase T1 and variants with buried and solvent-exposed Ala --> Gly mutations in the helix

Hydrogen-exchange rates were measured for RNase T1 and three variants with Ala --> Gly substitutions at a solvent-exposed (residue 21) and a buried (residue 23) position in the helix: A21G, G23A, and A21G + G23A. These results were used to measure the stabilities of the proteins. The hydrogen-exchange stabilities (DeltaG(HX)) for the most stable residues in each variant agree with the equilibrium conformational stability measured by urea denaturation (DeltaG(U)), if the effects of D(2)O and proline isomerization are included [Huyghues-Despointes, B. M. P., Scholtz, J. M., and Pace, C. N. (1999) Nat. Struct. Biol. 6, 210-212]. These residues also show similar changes in DeltaG(HX) upon Ala --> Gly mutations (DeltaDeltaG(HX)) as compared to equilibrium measurements (DeltaDeltaG(U)), indicating that the most stable residues are exchanging from the globally unfolded ensemble. Alanine is stabilizing compared to glycine by 1 kcal/mol at a solvent-exposed site 21 as seen by other methods for the RNase T1 protein and peptide helix [Myers, J. K., Pace, C. N., and Scholtz, J. M. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 3833-2837], while it is destabilizing at the buried site 23 by the same amount. For the A21G variant, only local NMR chemical shift perturbations are observed compared to RNase T1. For the G23A variant, large chemical shift changes are seen throughout the sequence, although X-ray crystal structures of the variant and RNase T1 are nearly superimposable. Ala --> Gly mutations in the helix of RNase T1 at both helical positions alter the native-state hydrogen-exchange stabilities of residues throughout the sequence.