Do all backbone polar groups in proteins form hydrogen bonds?

Evidence from proteins and peptides supports the conclusion that intrapeptide hydrogen bonds stabilize the folded form of proteins. Paradoxically, evidence from small molecules supports the opposite conclusion, that intrapeptide hydrogen bonds are less favorable than peptide-water hydrogen bonds. A related issue-often lost in this debate about comparing peptide-peptide to peptide- water hydrogen bonds-involves the energetic cost of an unsatisfied hydrogen bond. Here, experiment and theory agree that breaking a hydrogen bond costs between 5 and 6 kcal/mol. Accordingly, the likelihood of finding an unsatisfied hydrogen bond in a protein is insignificant. This realization establishes a powerful rule for evaluating protein conformations.     

Learning about protein hydrogen bonding by minimizing contrastive divergence

Defining the strength and geometry of hydrogen bonds in protein structures has been a challenging task since early days of structural biology. In this article, we apply a novel statistical machine learning technique, known as contrastive divergence, to efficiently estimate both the hydrogen bond strength and the geometric characteristics of strong interpeptide backbone hydrogen bonds, from a dataset of structures representing a variety of different protein folds. Despite the simplifying assumptions of the interatomic energy terms used, we determine the strength of these hydrogen bonds to be between 1.1 and 1.5 kcal/mol, in good agreement with earlier experimental estimates. The geometry of these strong backbone hydrogen bonds features an almost linear arrangement of all four atoms involved in hydrogen bond formation. We estimate that about a quarter of all hydrogen bond donors and acceptors participate in these strong interpeptide hydrogen bonds.     

Contribution of hydrogen bonding to the conformational stability of ribonuclease T1.

For 30 years, the prevailing view has been that the hydrophobic effect contributes considerably more than hydrogen bonding to the conformational stability of globular proteins. The results and reasoning presented here suggest that hydrogen bonding and the hydrophobic effect make comparable contributions to the conformational stability of ribonuclease T1 (RNase T1). When RNase T1 folds, 86 intramolecular hydrogen bonds with an average length of 2.95 A are formed. Twelve mutants of RNase T1 [Tyr----Phe (5), Ser----Ala (3), and Asn----Ala (4)] have been prepared that remove 17 of the hydrogen bonds with an average length of 2.93 A. On the basis of urea and thermal unfolding studies of these mutants, the average decrease in conformational stability due to hydrogen bonding is 1.3 kcal/mol per hydrogen bond. This estimate is in good agreement with results from several related systems. Thus, we estimate that hydrogen bonding contributes about 110 kcal/mol to the conformational stability of RNase T1 and that this is comparable to the contribution of the hydrophobic effect. Accepting the idea that intramolecular hydrogen bonds contribute 1.3 +/- 0.6 kcal/mol to the stability of systems in an aqueous environment makes it easier to understand the stability of the "molten globule" states of proteins, and the alpha-helical conformations of small peptides.     

Prudent modeling of core polar residues in computational protein design

Hydrogen bond interactions were surveyed in a set of protein structures. Compared to surface positions, polar side-chains at core positions form a greater number of intra-molecular hydrogen bonds. Furthermore, the majority of polar side-chains at core positions form at least one hydrogen bond to main-chain atoms that are not involved in hydrogen bonds to other main-chain atoms. Based on this structural survey, hydrogen bond rules were generated for each polar amino acid for use in protein core design. In the context of protein core design, these prudent polar rules were used to eliminate from consideration polar amino acid rotamers that do not form a minimum number of hydrogen bonds. As an initial test, the core of Escherichia coli thioredoxin was selected as a design target. For this target, the prudent polar strategy resulted in a minor increase in computational complexity compared to a strategy that did not allow polar residues. Dead-end elimination was used to identify global minimum energy conformations for the prudent polar and no polar strategies. The prudent polar strategy identified a protein sequence that was thermodynamically stabilized by 2.5 kcal/mol relative to wild-type thioredoxin and 2.2 kcal/mol relative to a thioredoxin variant whose core was designed without polar residues.     

Structural bases of lectin-carbohydrate affinities: comparison with protein-folding energetics.

We have made a comparative structure based analysis of the thermodynamics of lectin-carbohydrate (L-C) binding and protein folding. Examination of the total change in accessible surface area in those processes revealed a much larger decrease in free energy per unit of area buried in the case of L-C associations. According to our analysis, this larger stabilization of L-C interactions arises from a more favorable enthalpy of burying a unit of polar surface area, and from higher proportions of polar areas. Hydrogen bonds present at 14 L-C interfaces were identified, and their overall characteristics were compared to those reported before for hydrogen bonds in protein structures. Three major factors might explain why polar-polar interactions are stronger in L-C binding than in protein folding: (1) higher surface density of hydrogen bonds; (2) better hydrogen-bonding geometry; (3) larger proportion of hydrogen bonds involving charged groups. Theoretically, the binding entropy can be partitioned into three main contributions: entropy changes due to surface desolvation, entropy losses arising from freezing rotatable bonds, and entropic effects that result from restricting translation and overall rotation motions. These contributions were estimated from structural information and added up to give calculated binding entropies. Good correlation between experimental and calculated values was observed when solvation effects were treated according to a parametrization developed by other authors from protein folding studies. Finally, our structural parametrization gave calculated free energies that deviate from experimental values by 1.1 kcal/mol on the average; this amounts to an uncertainty of one order of magnitude in the binding constant.     

Ser45 plays an important role in managing both the equilibrium and transition state energetics of the streptavidin-biotin system

The contribution of the Ser45 hydrogen bond to biotin binding activation and equilibrium thermodynamics was investigated by biophysical and X-ray crystallographic studies. The S45A mutant exhibits a 1,700-fold greater dissociation rate and 907-fold lower equilibrium affinity for biotin relative to wild-type streptavidin at 37 degrees C, indicating a crucial role in binding energetics. The crystal structure of the biotin-bound mutant reveals only small changes from the wild-type bound structure, and the remaining hydrogen bonds to biotin retain approximately the same lengths. No additional water molecules are observed to replace the missing hydroxyl, in contrast to the previously studied D128A mutant. The equilibrium deltaG degrees, deltaH degrees, deltaS degrees, deltaC degrees(p), and activation deltaG++ of S45A at 37 degrees C are 13.7+/-0.1 kcal/mol, -21.1+/-0.5 kcal/mol, -23.7+/-1.8 cal/mol K, -223+/-12 cal/mol K, and 20.0+/-2.5 kcal/mol, respectively. Eyring analysis of the large temperature dependence of the S45A off-rate resolves the deltaH++ and deltaS++ of dissociation, 25.8+/-1.2 kcal/mol and 18.7+/-4.3 cal/mol K. The large increases of deltaH++ and deltaS++ in the mutant, relative to wild-type, indicate that Ser45 could form a hydrogen bond with biotin in the wild-type dissociation transition state, enthalpically stabilizing it, and constraining the transition state entropically. The postulated existence of a Ser45-mediated hydrogen bond in the wild-type streptavidin transition state is consistent with potential of mean force simulations of the dissociation pathway and with molecular dynamics simulations of biotin pullout, where Ser45 is seen to form a hydrogen bond with the ureido oxygen as biotin slips past this residue after breaking the native hydrogen bonds.     

Cooperative hydrogen bond interactions in the streptavidin-biotin system

The thermodynamic and structural cooperativity between the Ser45- and D128-biotin hydrogen bonds was measured by calorimetric and X-ray crystallographic studies of the S45A/D128A double mutant of streptavidin. The double mutant exhibits a binding affinity approximately 2x10(7) times lower than that of wild-type streptavidin at 25 degrees C. The corresponding reduction in binding free energy (DeltaDeltaG) of 10.1 kcal/mol was nearly completely due to binding enthalpy losses at this temperature. The loss of binding affinity is 11-fold greater than that predicted by a linear combination of the single-mutant energetic perturbations (8.7 kcal/mol), indicating that these two mutations interact cooperatively. Crystallographic characterization of the double mutant and comparison with the two single mutant structures suggest that structural rearrangements at the S45 position, when the D128 carboxylate is removed, mask the true energetic contribution of the D128-biotin interaction. Taken together, the thermodynamic and structural analyses support the conclusion that the wild-type hydrogen bond between D128-OD and biotin-N2 is thermodynamically stronger than that between S45-OG and biotin-N1.     

Thermodynamics of amide hydrogen bond formation in polar and apolar solvents

We present the initial findings of a theoretical study of hydrogen bond formation between two formamide molecules in water and in carbon tetrachloride. These systems were chosen as the simplest models for secondary structure formation in the polar environment near the protein surface and the apolar environment of the protein interior. We have employed thermodynamic simulation methods to obtain absolute binding free energies and free energy profiles for the formation of peptide hydrogen bonds in the two solvents. We find that the amide hydrogen bond is stable by 8.4 kcal/mol in CCl4, and by 0.3 kcal/mol in water. Our results indicate also that the hydrogen-bonded dimer is 2.2 kcal/mol more stable in water than it is in CCl4. We compare our results with those from experiment, and discuss their use in interpreting mechanisms of protein folding.     

Probing the role of hydrogen bonds in the stability of base pairs in double-helical DNA

Aromatic stacking and hydrogen bonding between nucleobases are two of the key interactions responsible for stabilization of DNA double-helical structures. The present work aims at defining the specific contributions of these interactions to the stability of individual base pairs in DNA. The two DNA double helices investigated are formed, respectively, by the palindromic base sequences 5'-dCCAACGTTGG-3' and 5'-dCGCAGATCTGCG-3'. The strength of the N==H...N inter-base hydrogen bond in each base pair is characterized from the measurement of the protium-deuterium fractionation factor of the corresponding imino proton using NMR spectroscopy. The structural stability of each base pair is evaluated from the exchange rate of the imino proton, measured by NMR. The results reveal that the fractionation factors of the imino protons in the two DNA double helices investigated fall within a narrow range of values, between 0.92 and 1.0. In contrast, the free energies of structural stabilization for individual base pairs span 3.5 kcal/mol, from 5.2 to 8.7 kcal/mol (at 15 degrees C). These findings indicate that, in the two DNA double helices investigated, the strength of N==H...N inter-base hydrogen bonds does not change significantly depending on the nature or the sequence context of the base pair. Hence, the variations in structural stability detected by proton exchange do not involve changes in the strength of inter-base hydrogen bonds. Instead, the results suggest that the energetic identity of a base pair is determined by the number of inter-base hydrogen bonds, and by the stacking interactions with neighboring base pairs.     

Metal-activated histidine carbon donor hydrogen bonds contribute to metalloprotein folding and function

Carbon donor hydrogen bonds are typically weak interactions that contribute less than 2 kcal/mol, and provide only modest stabilization in proteins. One exception is the class of hydrogen bonds donated by heterocyclic side chain carbons. Histidine is capable of particularly strong interactions through the Cε¹ and Cδ² carbons when the imidazole is protonated or bound to metal. Given the frequent occurrence of metal-bound histidines in metalloproteins, we characterized the energies of these interactions through DFT calculations on model compounds. Imidazole-water hydrogen bonding could vary from −11.0 to −17.0 kcal/mol, depending on the metal identity and oxidation state. A geometric search of metalloprotein structures in the PDB identified a number of candidate His C-H···O hydrogen bonds which may be important for folding or function. DFT calculations on model complexes of superoxide reductase show a carbon donor hydrogen bond positioning a water molecule above the active site.     

The contribution of polar group burial to protein stability is strongly context-dependent

We previously suggested that proteins gain more stability from the burial and hydrogen bonding of polar groups than from the burial of nonpolar groups (Pace, C. N. (2001) Biochemistry 40, 310-313). To study this further, we prepared eight Thr-to-Val mutants of RNase Sa, four in which the Thr side chain is hydrogen-bonded and four in which it is not. We measured the stability of these mutants by analyzing their thermal denaturation curves. The four hydrogen-bonded Thr side chains contribute 1.3 +/- 0.9 kcal/mol to the stability; those that are not still contribute 0.4 +/- 0.9 kcal/mol to the stability. For 40 Thr-to-Val mutants of 11 proteins, the average decrease in stability is 1.0 +/- 1.0 kcal/mol when the Thr side chain is hydrogen-bonded and 0.0 +/- 0.5 kcal/mol when it is not. This is clear evidence that hydrogen bonds contribute favorably to protein stability. In addition, we prepared four Val-to-Thr mutants of RNase Sa, measured their stability, and determined their crystal structures. In all cases, the mutants are less stable than the wild-type protein, with the decreases in stability ranging from 0.5 to 4.4 kcal/mol. For 41 Val-to-Thr mutants of 11 proteins, the average decrease in stability is 1.8 +/- 1.3 kcal/mol and is unfavorable for 40 of 41 mutants. This shows that placing an [bond]OH group at a site designed for a [bond]CH3 group is very unfavorable. So, [bond]OH groups can contribute favorably to protein stability, even if they are not hydrogen-bonded, if the site was selected for an [bond]OH group, but they will make an unfavorable contribution to stability, even if they are hydrogen-bonded, when they are placed at a site selected for a [bond]CH3 group. The contribution that polar groups make to protein stability depends strongly on their environment.     

Binding energetics of phosphorus-containing inhibitors of thermolysin

The importance of a specific hydrogen bond between thermolysin and a phosphonamidate inhibitor, Z-NHCH2-PO(O-)-Leu-Leu (1) [Bartlett, P. A., & Marlowe, C. K. (1987) Science (Washington D.C.) 235, 569-571], has been reevaluated. We have determined the inhibition constants (binding free energies) for thermolysin of phosphonamidate n-hexyl-P(O)(O-)-Leu-Trp-NHMe (4), phosphonate n-hexyl-P-(O)(O-)OCH(iBu)CO-Trp-NHMe (5), and phosphinates n-hexyl-P(O)(O-)CH2CH(iBu)CO-Trp-NHMe (6) and Z-NHCH2PO(O-)CH2CH(iBu)CO-Leu (3). Replacement of the P-NH group by P-CH2 (1----3 and 4----6) weakens the overall binding free energy by about 1.5 kcal/mol. A negligible difference in solvation energy has been measured for these pairs, and the basicity of the P-O- ligand for zinc in each pair remains nearly unchanged as determined by pH titration of their 31P NMR resonances. Therefore, this value of 1.5 kcal/mol can be assigned to the specific hydrogen bond known to exist between the P-NH of 1 and thermolysin [Tronrud, D. E., Holden, H. M., & Matthews, B. W. (1987) Science (Washington, D.C.) 235, 871-574] and inferred to exist between 4 and the enzyme. Substitution of P-O for P-NH (1----2 [Bartlett, P. A., & Marlowe, C. K. (1987) Science (Washington, D.C.) 235, 569-571] and 4----5) weakens the overall binding free energy by 4.1 kcal/mol for each pair as the basicity of the P-O- ligand decreases by about 1.6 pH units. The measured solvation energy difference between 4 and 5 (and by inference between 1 and 2) is negligible.(ABSTRACT TRUNCATED AT 250 WORDS)     

A low-barrier hydrogen bond between histidine of secreted phospholipase A2 and a transition state analog inhibitor

This work describes in-depth NMR characterization of a unique low-barrier hydrogen bond (LBHB) between an active site residue from the enzyme and a bound inhibitor: the complex between secreted phospholipase A(2) (sPLA(2), from bee venom and bovine pancreas) and a transition-state analog inhibitor HK32. A downfield proton NMR resonance, at 17-18 ppm, was observed in the complex but not in the free enzyme. On the basis of site-specific mutagenesis and specific 15N-decoupling, this downfield resonance was assigned to the active site H48, which is part of the catalytic dyad D99-H48. These results led to a hypothesis that the downfield resonance represents the proton (H(epsilon 2) of H48) involved in the H-bonding between D99 and H48, in analogy with serine proteases. However, this was shown not to be the case by use of the bovine enzyme labeled with specific [15N(epsilon 2)]His. Instead, the downfield resonance arises from H(delta1) of H48, which forms a hydrogen bond with a non-bridging phosphonate oxygen of the inhibitor. Further studies showed that this proton displays a fractionation factor of 0.62(+/-0.06), and an exchange rate protection factor of >100 at 285 K and >40 at 298 K, which are characteristic of a LBHB. The pK(a) of the imidazole ring of H48 was shown to be shifted from 5.7 for the free enzyme to an apparent value of 9.0 in the presence of the inhibitor. These properties are very similar to those of the Asp em leader His LBHBs in serine proteases. Possible structural bases and functional consequences for the different locations of the LBHB between these two types of enzymes are discussed. The results also underscore the importance of using specific isotope labeling, rather than extrapolation of NMR results from other enzyme systems, to assign the downfield proton resonance to a specific hydrogen bond. Although our studies did not permit the strength of the LBHB to be accurately measured, the data do not provide support for an unusually strong hydrogen bond strength (i.e. >10 kcal/mol).     

Hydrogen bonding stabilizes globular proteins.

It is clear that intramolecular hydrogen bonds are essential to the structure and stability of globular proteins. It is not clear, however, whether they make a net favorable contribution to this stability. Experimental and theoretical studies are at odds over this important question. Measurements of the change in conformational stability, delta (delta G), for the mutation of a hydrogen bonded residue to one incapable of hydrogen bonding suggest a stabilization of 1.0 kcal/mol per hydrogen bond. If the delta (delta G) values are corrected for differences in side-chain hydrophobicity and conformational entropy, then the estimated stabilization becomes 2.2 kcal/mol per hydrogen bond. These and other experimental studies discussed here are consistent and compelling: hydrogen bonding stabilizes globular proteins.     

Free energy simulations and MM-PBSA analyses on the affinity and specificity of steroid binding to antiestradiol antibody

Antiestradiol antibody 57-2 binds 17beta-estradiol (E2) with moderately high affinity (K(a) = 5 x 10(8) M(-1)). The structurally related natural estrogens estrone and estriol as well synthetic 17-deoxy-estradiol and 17alpha-estradiol are bound to the antibody with 3.7-4.9 kcal mol(-1) lower binding free energies than E2. Free energy perturbation (FEP) simulations and the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were applied to investigate the factors responsible for the relatively low cross-reactivity of the antibody with these four steroids, differing from E2 by the substituents of the steroid D-ring. In addition, computational alanine scanning of the binding site residues was carried out with the MM-PBSA method. Both the FEP and MM-PBSA methods reproduced the experimental relative affinities of the five steroids in good agreement with experiment. On the basis of FEP simulations, the number of hydrogen bonds formed between the antibody and steroids, which varied from 0 to 3 in the steroids studied, determined directly the magnitude of the steroid-antibody interaction free energies. One hydrogen bond was calculated to contribute about 3 kcal mol(-1) to the interaction energy. Because the relative binding free energies of estrone (two antibody-steroid hydrogen bonds), estriol (three hydrogen bonds), 17-deoxy-estradiol (no hydrogen bonds), and 17alpha-estradiol (two hydrogen bonds) are close to each other and clearly lower than that of E2 (three hydrogen bonds), the water-steroid interactions lost upon binding to the antibody make an important contribution to the binding free energies. The MM-PBSA calculations showed that the binding of steroids to the antiestradiol antibody is driven by van der Waals interactions, whereas specificity is solely due to electrostatic interactions. In addition, binding of steroids to the antiestradiol antibody 57-2 was compared to the binding to the antiprogesterone antibody DB3 and antitestosterone antibody 3-C4F5, studied earlier with the MM-PBSA method.     

Chrysophanol anthrone: a theoretical study on the potential energy surface

The potential energy surface (PES) of chrysophanol anthrone, the active component of Goa Powder, was systematically explored and thoroughly scrutinised via density functional theory, in order to gain an understanding of its physicochemical properties. In particular, we focused on the rotations of the two hydroxyl-phenyl dihedral angles. A picture with a four stable rotamers emerged where only the most stable conformer has a planar structure and the less stable conformer has the maximum deviation from planarity. The computed PES shows that the energy barriers for the conformer interconversion are less than 15 kcal/mol. From the analysis of the calculated intramolecular hydrogen bond enthalpy, we conclude that the number of the intramolecular hydrogen bonds governs the conformer stability. Additionally, the conformational equilibrium was pursued by means of an analysis of the energy of OH internal rotation barriers. The total energy changes were decomposed in an electrostatic decomposition scheme in order to gain an insight into the effects governing the torsional barrier and preferred conformations. This analysis shows that the interplay between the repulsive and attractive potentials causes the conformer stability, where the attractive term dominates the conformer stabilisation.     

Stability of a model beta-sheet in water.

We have used molecular dynamics simulations to determine the stability in water of a model beta-sheet formed by two alanine dipeptide molecules with two intermolecular hydrogen bonds in the closely spaced antiparallel arrangement. In this paper we describe our computations of the binding free energy of the model sheet and a portion of the free energy surface as a function of a reaction co-ordinate for sheet formation. We used the free energy surface to identify stable conformations along the reaction co-ordinate. To determine whether or not the model sheet with two hydrogen bonds is more stable than a single amide hydrogen bond in water, we compared the results of the present calculations to results from our earlier study of linear hydrogen bond formation between two formamide molecules (the formamide "dimer"). The free energy surfaces for the sheet and formamide dimer each have two minima corresponding to locally stable hydrogen-bonded and solvent-separated configurations. The binding free energies of the model sheet and the formamide dimer are -5.5 and -0.34 kcal/mol, respectively. Thus, the model sheet with two hydrogen bonds is quite stable while the simple amide hydrogen bond is only marginally stable. To understand the relative stabilities of the model sheet and formamide dimer in terms of solute-solute and solute-water interactions, we decomposed the free energy differences between hydrogen-bonded and solvent-separated conformations into energetic and entropic contributions. The changes in the peptide-peptide energy and the entropy are roughly twice as large for the sheet as they are for the formamide dimer. The magnitude of the peptide-water energy difference for the sheet is less than twice (by about 3.5 kcal/mol) that for the formamide dimer, and this accounts for the stability of the sheet. The presence of the side-chains and/or blocking groups apparently prevents the amide groups in the sheet from being solvated as favorably in the separated arrangement as in the formamide dimer, where the amide groups are completely exposed to the solvent.     

Semiquantitative calculations of catalytic free energies in genetically modified enzymes

The catalytic free energy and binding free energies of the native and the Asn-155----Thr, Asn-155----Leu, and Asn-155----Ala mutants of subtilisin are calculated by the empirical valence bond method and a free energy perturbation method. Two simple procedures are used; one "mutates" the substrate, and the other "mutates" the enzyme. The calculated changes in free energies (delta delta G not equal to cat and delta delta Gbind) between the mutant and native enzymes are within 1 kcal/mol of the corresponding observed values. This indicates that we are approaching a quantitative structure-function correlation. The calculated changes in catalytic free energies are almost entirely due to the electrostatic interaction between the enzyme-water system and the charges of the reacting system. This supports the idea that the electrostatic free energy associated with the changes of charges of the reacting system is the key factor in enzyme catalysis.     

Contribution of a tyrosine side chain to ribonuclease A catalysis and stability.

An intricate architecture of covalent bonds and noncovalent interactions appear to position the side chain of Lys 41 properly within the active site of bovine pancreatic ribonuclease A (RNase A). One of these interactions arises from Tyr 97, which is conserved in all 41 RNase A homologues of known sequence. Tyr 97 has a solvent-inaccessible side chain that donates a hydrogen bond to the main-chain oxygen of Lys 41. Here, the role of Tyr 97 was examined by replacing Tyr 97 with a phenylalanine, alanine, or glycine residue. All three mutant proteins have diminished catalytic activity, with the value of Kcat being perturbed more significantly than that of Km. The free energies with which Y97F, Y97A, and Y97G RNase A bind to the rate-limiting transition state during the cleavage of poly(cytidylic acid) are diminished by 0.74, 3.3, and 3.8 kcal/mol, respectively. These results show that even though Tyr 97 is remote from the active site, its side chain contributes to catalysis. The role of Tyr 97 in the thermal stability of RNase A is large. The conformational free energies of native Y97F, Y97A, and Y97G RNase A are decreased by 3.54, 12.0, and 11.7 kcal/mol, respectively. The unusually large decrease in stability caused by the Tyr-->Phe mutation could result from a decrease in the barrier to isomerization of the Lys 41-Pro 42 peptide bond.     

Genetic selection reveals the role of a buried, conserved polar residue

The burial of nonpolar surface area is known to enhance markedly the conformational stability of proteins. The contribution from the burial of polar surface area is less clear. Here, we report on the tolerance to substitution of Ser75 of bovine pancreatic ribonuclease (RNase A), a residue that has the unusual attributes of being buried, conserved, and polar. To identify variants that retain biological function, we used a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity. Cell growth at 30 degrees C, 37 degrees C, and 44 degrees C correlated with residue size, indicating that the primary attribute of Ser75 is its small size. The side-chain hydroxyl group of Ser75 forms a hydrogen bond with a main-chain nitrogen. The conformational stability of the S75A variant, which lacks this hydrogen bond, was diminished by DeltaDeltaG = 2.5 kcal/mol. Threonine, which can reinstate this hydrogen bond, provided a catalytically active RNase A variant at higher temperatures than did some smaller residues (including aspartate), indicating that a secondary attribute of Ser75 is the ability of its uncharged side chain to accept a hydrogen bond. These results provide insight on the imperatives for the conservation of a buried polar residue.