Tyrosine hydrogen bonds make a large contribution to protein stability

The aim of this study was to gain a better understanding of the contribution of hydrogen bonds by tyrosine -OH groups to protein stability. The amino acid sequences of RNases Sa and Sa3 are 69 % identical and each contains eight Tyr residues with seven at equivalent structural positions. We have measured the stability of the 16 tyrosine to phenylalanine mutants. For two equivalent mutants, the stability increases by 0.3 kcal/mol (RNase Sa Y30F) and 0.5 kcal/mol (RNase Sa3 Y33F) (1 kcal=4.184 kJ). For all of the other mutants, the stability decreases with the greatest decrease being 3.6 kcal/mol for RNase Sa Y52F. Seven of the 16 tyrosine residues form intramolecular hydrogen bonds and the average decrease in stability for these is 2.0(+/-1.0) kcal/mol. For the nine tyrosine residues that do not form intramolecular hydrogen bonds, the average decrease in stability is 0.4(+/-0.6) kcal/mol. Thus, most tyrosine -OH groups contribute favorably to protein stability even if they do not form intramolecular hydrogen bonds. Generally, the stability changes for equivalent positions in the two proteins are remarkably similar. Crystal structures were determined for two of the tyrosine to phenylalanine mutants of RNase Sa: Y80F (1.2 A), and Y86F (1.7 A). The structures are very similar to that of wild-type RNase Sa, and the hydrogen bonding partners of the tyrosine residues always form intermolecular hydrogen bonds to water in the mutants. These results provide further evidence that the hydrogen bonding and van der Waals interactions of polar groups in the tightly packed interior of folded proteins are more favorable than similar interactions with water in the unfolded protein, and that polar group burial makes a substantial contribution to protein stability.     

Buried, charged, non-ion-paired aspartic acid 76 contributes favorably to the conformational stability of ribonuclease T1.

The side-chain carboxyl of Asp 76 in ribonuclease T1 (RNase T1) is buried, charged, non-ion-paired, and forms three good intramolecular hydrogen bonds (2.63, 2.69, and 2.89 A) and a 2.66 A hydrogen bond to a buried, conserved water molecule. When Asp 76 was replaced by Asn, Ser, and Ala, the conformational stability of the protein decreased by 3.1, 3.2, and 3.7 kcal/mol, respectively. The stability was measured as a function of pH for wild-type RNase T1 and the D76N mutant and showed that the pH dependence below pH 3 was almost entirely due to Asp 76. The pK of Asp 76 is 0.5 in the native state and 3.7 in the denatured state. Thus, the hydrogen bonding of the carboxyl group of Asp 76 contributes more than half of the net stability of RNase T1 at pH 7. In addition, the charged carboxyl of Asp 76 stabilizes structure in the denatured states of RNase T1 that is not present in D76N, D76S, and D76A.     

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.     

Contribution of hydrogen bonds to protein stability

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.     

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.     

Increasing protein stability by altering long-range coulombic interactions.

It is difficult to increase protein stability by adding hydrogen bonds or burying nonpolar surface. The results described here show that reversing the charge on a side chain on the surface of a protein is a useful way of increasing stability. Ribonuclease T1 is an acidic protein with a pI approximately 3.5 and a net charge of approximately -6 at pH 7. The side chain of Asp49 is hyperexposed, not hydrogen bonded, and 8 A from the nearest charged group. The stability of Asp49Ala is 0.5 kcal/mol greater than wild-type at pH 7 and 0.4 kcal/mol less at pH 2.5. The stability of Asp49His is 1.1 kcal/mol greater than wild-type at pH 6, where the histidine 49 side chain (pKa = 7.2) is positively charged. Similar results were obtained with ribonuclease Sa where Asp25Lys is 0.9 kcal/mol and Glu74Lys is 1.1 kcal/mol more stable than the wild-type enzyme. These results suggest that protein stability can be increased by improving the coulombic interactions among charged groups on the protein surface. In addition, the stability of RNase T1 decreases as more hydrophobic aromatic residues are substituted for Ala49, indicating a reverse hydrophobic effect.     

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a –CH₂− group on folding contributes 1.1 ± 0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1 ± 0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.     

The severity of osteogenesis imperfecta: a comparison to the relative free energy differences of collagen model peptides

Molecular dynamics simulations were carried out to calculate free energy differences between the folded and unfolded states of wild type and mutant collagen model peptides. The calculated stability of the collagen models was compared with the severity of osteogenesis imperfecta. Free energy differences of Gly → Xaa (Xaa: Ser, Cys, Glu, and Asp) mutations between the wild type and the mutants at position 15 of the model peptide were 3.8, 4.2, 5.6, and 8.8 kcal/mol, respectively. The corresponding free energy differences of a second Gly mutation at the same position in different chains were, on average, 1.3, 1.5, 2.9, and 5.4 kcal/mol, respectively. Free energy simulations were also performed to estimate the relative stability between an oxidized form and a reduced form of the mutants containing two Cys residues, which indicated that the mutant of the collagen-like peptide containing an intramolecular disulfide bond was more stable than the mutant containing one Cys residue but less stable than the wild type. The calculated free energy differences between an oxidized and a reduced form of the mutants containing two Cys residues are 0.8 and 2.6 kcal/mol for the disulfide bonds between Chains A and B and between Chains A and C, respectively.     

Conformational stability and activity of ribonuclease T1 with zero, one, and two intact disulfide bonds

Ribonuclease T1 has two disulfide bonds linking cysteine residues 2-10 and 6-103. We have prepared a derivative of ribonuclease T1 in which one disulfide bond is broken and the cysteine residues carboxymethylated, (2-10)-RCM-T1, and three derivatives in which both disulfides are broken and the cysteine residues reduced, R-T1, carboxamidomethylated, RCAM-T1, or carboxymethylated, RCM-T1. The RNA hydrolyzing activity of these proteins has been measured, and urea and thermal denaturation studies have been used to determine conformational stability. The activity, melting temperature, and conformational stability of the proteins are: ribonuclease T1 (100%, 59.3 degrees C, 10.2 kcal/mol), (2-10)-RCM-T1 (86%, 53.3 degrees C, 6.8 kcal/mol), R-T1 (53%, 27.2 degrees C, 3.0 kcal/mol), RCAM-T1 (43%, 21.2 degrees C, 1.5 kcal/mol), and RCM-T1 (35%, 16.6 degrees C, 0.9 kcal/mol). Thus, the conformational stability is decreased by 3.4 kcal/mol when one disulfide bond is broken and by 7.2-9.3 kcal/mol when both disulfide bonds are broken. It is quite remarkable that RNase T1 can fold and function with both disulfide bonds broken and the cysteine residues carboxymethylated. The large decrease in the stability is due mainly to an increase in the conformational entropy of the unfolded protein which results when the constraints of the disulfide bonds on the flexibility are removed. We propose a new equation for predicting the effect of a cross-link on the conformational entropy of a protein: delta Sconf = -2.1 - (3/2)R 1n n, where n is the number of residues between the side chains which are cross-linked. This equation gives much better agreement with experimental results than other forms of this equation which have been used previously.     

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6 ± 0.3 kcal/mol per -CH₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH₂- group on folding contributes, on average, 1.1 ± 0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG_tr values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60 ± 4% and hydrogen bonds contribute 40 ± 4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.     

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.     

An energetic scale for equilibrium H/D fractionation factors illuminates hydrogen bond free energies in proteins

Equilibrium H/D fractionation factors have been extensively employed to qualitatively assess hydrogen bond strengths in protein structure, enzyme active sites, and DNA. It remains unclear how fractionation factors correlate with hydrogen bond free energies, however. Here we develop an empirical relationship between fractionation factors and free energy, allowing for the simple and quantitative measurement of hydrogen bond free energies. Applying our empirical relationship to prior fractionation factor studies in proteins, we find: [1] Within the folded state, backbone hydrogen bonds are only marginally stronger on average in α‐helices compared to β‐sheets by ～0.2 kcal/mol. [2] Charge‐stabilized hydrogen bonds are stronger than neutral hydrogen bonds by ～2 kcal/mol on average, and can be as strong as –7 kcal/mol. [3] Changes in a few hydrogen bonds during an enzyme catalytic cycle can stabilize an intermediate state by –4.2 kcal/mol. [4] Backbone hydrogen bonds can make a large overall contribution to the energetics of conformational changes, possibly playing an important role in directing conformational changes. [5] Backbone hydrogen bonding becomes more uniform overall upon ligand binding, which may facilitate participation of the entire protein structure in events at the active site. Our energetic scale provides a simple method for further exploration of hydrogen bond free energies.     

Contribution of the hydrophobic effect to globular protein stability.

The decrease in conformational stability, delta(delta G), has been measured for 72 aliphatic side-chain mutants from four proteins in which a larger side-chain is replaced by a smaller side-chain so that steric effects are minimal. When these delta(delta G) values are corrected to the same accessibility, namely 100% buried, then the following -delta(delta G) values per -CH2- group (in kcal/mol) are obtained: Ile----Val (1.26), Ala (1.26), Gly (1.26); Leu----Ala (1.16), Gly (1.21); Val----Ala (1.23), Gly (1.53). The average of these values is 1.27(+/- 0.07) kcal/mol. The 72 individual values range from 0 to 2.4 kcal/mol with an average value of 1.27(+/- 0.51) (standard deviation) kcal/mol. When the delta Gtr values from n-octanol to water are corrected for the difference in volume between the solutes and the solvents, the average value for the same substitutions is 1.25(+/- 0.05) kcal/mol. This suggests that proteins gain 1.3(+/- 0.5) kcal/mol in stability for each -CH2- group buried in folding, and, furthermore, that the volume corrected delta Gtr values for n-octanol for the amino acid side-chains provide good estimates of the contribution of the hydrophobic effect to globular protein stability.     

Contribution of histidine residues to the conformational stability of ribonuclease T1 and mutant Glu-58----Ala

The pK values of the histidine residues in ribonuclease T1 (RNase T1) are unusually high: 7.8 (His-92), 7.9 (His-40), and 7.3 (His-27) [Inagaki et al. (1981) J. Biochem. 89, 1185-1195]. In the RNase T1 mutant Glu-58----Ala, the first two pK values are reduced to 7.4 (His-92) and 7.1 (His-40). These lower pKs were expected since His-92 (5.5 A) and His-40 (3.7 A) are in close proximity to Glu-58 at the active site. The conformational stability of RNase T1 increases by over 4 kcal/mol between pH 9 and 5, and this can be entirely accounted for by the greater affinity for protons by the His residues in the folded protein (average pK = 7.6) than in the unfolded protein (pk approximately 6.6). Thus, almost half of the net conformational stability of RNase T1 results from a difference between the pK values of the histidine residues in the folded and unfolded conformations. In the Glu-58----Ala mutant, the increase in stability between pH 9 and 5 is halved (approximately 2 kcal/mol), as expected on the basis of the lower pK values for the His residues in the folded protein (average pK = 7.1). As a consequence, RNase T1 is more stable than the mutant below pH 7.5, and less stable above pH 7.5. These results emphasize the importance of measuring the conformational stability as a function of pH when comparing proteins differing in structure.     

Completely buried, non-ion-paired glutamic acid contributes favorably to the conformational stability of pyrrolidone carboxyl peptidases from hyperthermophiles

Pyrrolidone carboxyl peptidases (PCPs) from hyperthermophiles have a structurally conserved and completely buried Glu192 in the hydrophobic core; in contrast, the corresponding residue in the mesophile protein is a hydrophobic residue, Ile. Does the buried ionizable residue contribute to stabilization or destabilization of hyperthermophile PCPs? To elucidate the role of the buried glutamic acid in stabilizing PCP from hyperthermophiles, we constructed five Glu192 mutants of PCP-0SH (C142S/C188S, Cys-free double mutant of PCP) from Pyrococcus furiosus and examined their thermal and pH-induced unfolding and crystal structures and compared them with those of PCP-0SH. The stabilities of apolar (E192A/I/V) and polar (E192D/Q) mutants were less than PCP-0SH at acidic pH values. In the alkaline region, the mutant proteins, except for E192D, were more stable than PCP-0SH. The thermal stability data and theoretical calculations indicated an apparent pKa value > or = 7.3 for Glu192. Present results confirmed that the protonated Glu192 in PCP-0SH forms strong hydrogen bonds with the carbonyl oxygen and peptide nitrogen of Pro168. New intermolecular hydrogen bonds in the E --> A/D mutants were formed by a water molecule introduced into the cavity created around position 192, whereas the hydrogen bonds disappeared in the E --> I/V mutants. Structure-based empirical stability of mutant proteins was in good agreement with the experimental results. The results indicated that (1) completely buried Glu192 contributes to the stabilization of PCP-0SH because of the formation of strong intramolecular hydrogen bonds and (2) the hydrogen bonds by the nonionized and buried Glu can contribute more than the burial of hydrophobic groups to the conformational stability of proteins.     

Destabilizing effect of proline substitutions in two helical regions of T4 lysozyme: leucine 66 to proline and leucine 91 to proline.

A class of temperature-sensitive (ts) mutants of T4 lysozyme with reduced activity at 30 degrees C and no activity at 43 degrees C has been selected. These mutants, designated "tight" ts mutants, differ from most other T4 lysozyme mutants that are active at 43 degrees C, but only manifest their ts lesion by a reduced halo size around phage plaques after exposure of the growth plates to chloroform vapors. For example, in the series of T4 lysozyme mutants at position 157, the original randomly selected mutant, T1571, is the least stable of the series, yet, apart from the halo assay and subsequent in vitro protein stability measurements, this mutant is indistinguishable from wild type (WT) even at 43 degrees C. Two mutants were identified: L91P and L66P. Both insert proline residues into alpha-helical regions of the WT protein structure. The stabilities (delta delta G) as determined by urea denaturation are 8.2 kcal/mol for L91P and 7.1 kcal/mol for L66P. CD spectra indicate that no major conformational changes have occurred in the mutant structures. The structures of the mutants were modeled with a 40-ps molecular dynamics simulation using explicit solvent. For L91P, the reduction of stability appears to be due to an unsatisfied hydrogen bond in the alpha-helix and to a new buried cavity. For L66P, the reduction of stability appears to be due to a disruption of the interdomain alpha-helix, at least two unsatisfied hydrogen bonds, and a newly formed solvent-filled pocket that protrudes into the hydrophobic core, possibly reducing the stabilizing contribution of a partially buried intrachain salt bridge.     

Contribution of hydrogen bonds to the conformational stability of human lysozyme: calorimetry and X-ray analysis of six Ser --> Ala mutants.

To further examine the contribution of hydrogen bonds to the conformational stability of the human lysozyme, six Ser to Ala mutants were constructed. The thermodynamic parameters for denaturation of these six Ser mutant proteins were investigated by differential scanning calorimetry (DSC), and the crystal structures were determined by X-ray analysis. The denaturation Gibbs energy (DeltaG) of the Ser mutant proteins was changed from 2.0 to -5.7 kJ/mol, compared to that of the wild-type protein. With an analysis in which some factors that affected the stability due to mutation were considered, the contribution of hydrogen bonds to the stability (Delta DeltaGHB) was extracted on the basis of the structures of the mutant proteins. The results showed that hydrogen bonds between protein atoms and between a protein atom and a water bound with the protein molecule favorably contribute to the protein stability. The net contribution of one intramolecular hydrogen bond to protein stability (DeltaGHB) was 8.9 +/- 2.6 kJ/mol on average. However, the contribution to the protein stability of hydrogen bonds between a protein atom and a bound water molecule was smaller than that for a bond between protein atoms.     

Hydrogen bonding markedly reduces the pK of buried carboxyl groups in proteins

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.     

Energy of hydrogen bonds probed by the adhesion of functionalized lipid layers

It is now well admitted that hydrophobic interactions and hydrogen bonds are the main forces driving protein folding and stability. However, because of the complex structure of a protein, it is still difficult to separate the different energetic contributions and have a reliable estimate of the hydrogen bond part. This energy can be quantified on simpler systems such as surfaces bearing hydrogen-bonding groups. Using the surface force apparatus, we have directly measured the interaction energy between monolayers of lipids whose headgroups can establish hydrogen bonds in water: nitrilotriacetate, adenosine, thymidine, and methylated thymidine lipids. From the adhesion energy between the surfaces, we have deduced the energy of a single hydrogen bond in water. We found in each case an energy of 0.5 kcal/mol. This result is in good agreement with recent experimental and theoretical studies made on protein systems showing that intramolecular hydrogen bonds make a positive contribution to protein stabilization.     

Stabilization of apoflavodoxin by replacing hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or Gln

Knowledge of protein stability principles provides a means to increase protein stability in a rational way. Here we explore the feasibility of stabilizing proteins by replacing solvent-exposed hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or GLN: The rationale behind this is a previous observation that, in some cases, neutral hydrogen bonds may be more stable that charged ones. We identified, in the apoflavodoxin from Anabaena PCC 7119, three surface-exposed aspartate or glutamate residues involved in hydrogen bonding with a single partner and we mutated them to asparagine or glutamine, respectively. The effect of the mutations on apoflavodoxin stability was measured by both urea and temperature denaturation. We observed that the three mutant proteins are more stable than wild-type (on average 0.43 kcal/mol from urea denaturation and 2.8 degrees C from a two-state analysis of fluorescence thermal unfolding data). At high ionic strength, where potential electrostatic repulsions in the acidic apoflavodoxin should be masked, the three mutants are similarly more stable (on average 0.46 kcal/mol). To rule out further that the stabilization observed is due to removal of electrostatic repulsions in apoflavodoxin upon mutation, we analysed three control mutants and showed that, when the charged residue mutated to a neutral one is not hydrogen bonded, there is no general stabilizing effect. Replacing hydrogen-bonded charged Asp or Glu residues by Asn or Gln, respectively, could be a straightforward strategy to increase protein stability.