Satisfaction of hydrogen-bonding potential influences the conservation of polar sidechains

Although polar amino acids tend to be found on the surface of proteins due to their hydrophilic nature, their important roles within the core of proteins are now becoming better recognized. It has long been understood that a significant number of mainchain functions will not achieve hydrogen bond satisfaction through the formation of secondary structures; in these circumstances, it is generally buried polar residues that provide hydrogen bond satisfaction. Here, we describe an analysis of the hydrogen-bonding of polar amino acids in a set of structurally aligned protein families. This allows us not only to calculate the conservation of each polar residue but also to assess whether conservation is correlated with the hydrogen-bonding potential of polar sidechains. We show that those polar sidechains whose hydrogen-bonding potential is satisfied tend to be more conserved than their unsatisfied or nonhydrogen-bonded counterparts, particularly when buried. Interestingly, these buried and satisfied polar residues are significantly more conserved than buried hydrophobic residues. Forming hydrogen bonds to mainchain amide atoms also influences conservation, with those satisfied buried polar residues that form two hydrogen bonds to mainchain amides being significantly more conserved than those that form only one or none. These results indicate that buried polar residues whose hydrogen-bonding potential is satisfied are likely to have important roles in maintaining protein structure.     

Hydrogen bonding increases packing density in the protein interior

The contribution of hydrogen bonds and the burial of polar groups to protein stability is a controversial subject. Theoretical studies suggest that burying polar groups in the protein interior makes an unfavorable contribution to the stability, but experimental studies show that burying polar groups, especially those that are hydrogen bonded, contributes favorably to protein stability. Understanding the factors that are not properly accounted for by the theoretical models would improve the models so that they more accurately describe experimental results. It has been suggested that hydrogen bonds may contribute to protein stability, in part, by increasing packing density in the protein interior, and thereby increasing the contribution of van der Waals interactions to protein stability. To investigate the influence of hydrogen bonds on packing density, we analyzed 687 crystal structures and determined the volume of buried polar groups as a function of their extent of hydrogen bonding. Our findings show that peptide groups and polar side chains that form hydrogen bonds occupy a smaller volume than the same groups when they do not form hydrogen bonds. For example, peptide groups in which both polar groups are hydrogen bonded occupy a volume, on average, 5.2 A3 less than a peptide group that is not hydrogen bonded.     

Analysing the ability to retain sidechain hydrogen-bonds in mutant proteins

MOTIVATION: Hydrogen bonds are one of the most important inter-atomic interactions in biology. Previous experimental, theoretical and bioinformatics analyses have shown that the hydrogen bonding potential of amino acids is generally satisfied and that buried unsatisfied hydrogen-bond-capable residues are destabilizing. When studying mutant proteins, or introducing mutations to residues involved in hydrogen bonding, one needs to know whether a hydrogen bond can be maintained. Our aim, therefore, was to develop a rapid method to evaluate whether a sidechain can form a hydrogen-bond. RESULTS: A novel knowledge-based approach was developed in which the conformations accessible to the residues involved are taken into account. Residues involved in hydrogen bonds in a set of high resolution crystal structures were analyzed and this analysis is then applied to a given protein. The program was applied to assess mutations in the tumour-suppressor protein, p53. This raised the number of distinct mutations identified as disrupting sidechain-sidechain hydrogen bonding from 181 in our previous analysis to 202 in this analysis.     

Statistical analysis of protein structures suggests that buried ionizable residues in proteins are hydrogen bonded or form salt bridges

It is well known that nonpolar residues are largely buried in the interior of proteins, whereas polar and ionizable residues tend to be more localized on the protein surface where they are solvent exposed. Such a distribution of residues between surface and interior is well understood from a thermodynamic point: nonpolar side chains are excluded from the contact with the solvent water, whereas polar and ionizable groups have favorable interactions with the water and thus are preferred at the protein surface. However, there is an increasing amount of information suggesting that polar and ionizable residues do occur in the protein core, including at positions that have no known functional importance. This is inconsistent with the observations that dehydration of polar and in particular ionizable groups is very energetically unfavorable. To resolve this, we performed a detailed analysis of the distribution of fractional burial of polar and ionizable residues using a large set of ?2600 nonhomologous protein structures. We show that when ionizable residues are fully buried, the vast majority of them form hydrogen bonds and/or salt bridges with other polar/ionizable groups. This observation resolves an apparent contradiction: the energetic penalty of dehydration of polar/ionizable groups is paid off by favorable energy of hydrogen bonding and/or salt bridge formation in the protein interior. Our conclusion agrees well with the previous findings based on the continuum models for electrostatic interactions in proteins. Proteins 2011; © 2011 Wiley‐Liss, Inc.     

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

Conformational transitions in eosinophil cationic protein: a molecular dynamics study in aqueous environment

Extensive molecular dynamics simulations have been performed on eosinophil cationic protein (ECP). The two structures found in the crystallographic dimer (ECPA and ECPB) have been independently simulated. A small difference in the pattern of the sidechain hydrogen bonds in the starting structure has resulted in interesting differences in the conformations accessed during the simulations. In one simulation (ECPB), a stable equilibrium conformation was obtained and in the other (ECPA), conformational transitions at the level of sidechain interactions were observed. The conformational transitions exhibit the involvement of the solvent (water) molecules with a pore-like construct in the equilibrium conformation and an opening for a large number of water molecules during the transition phase. The details of these transitions are examined in terms of intra-protein hydrogen bonds, protein-water networks and the residence times of water molecules on the polar atoms of the protein. These properties show some significant differences in the region between the N-terminal helix and the loop before the C-terminal strand as a function of different conformations accessed during the simulations. However, the stable hydrogen bonds, the protein-water networks, and the hydration patterns in most part of the protein including the active site are very much similar in both the simulations, indicating the fact that these are intrinsic properties of proteins.     

Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.     

Contribution of polar groups in the interior of a protein to the conformational stability

It has been generally believed that polar residues are usually located on the surface of protein structures. However, there are many polar groups in the interior of the structures in reality. To evaluate the contribution of such buried polar groups to the conformational stability of a protein, nonpolar to polar mutations (L8T, A9S, A32S, I56T, I59T, I59S, A92S, V93T, A96S, V99T, and V100T) in the interior of a human lysozyme were examined. The thermodynamic parameters for denaturation were determined using a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. If a polar group had a heavy energy cost to be buried, a mutant protein would be remarkably destabilized. However, the stability (Delta G) of the Ala to Ser and Val to Thr mutant human lysozymes was comparable to that of the wild-type protein, suggesting a low-energy penalty of buried polar groups. The structural analysis showed that all polar side chains introduced in the mutant proteins were able to find their hydrogen bond partners, which are ubiquitous in protein structures. The empirical structure-based calculation of stability change (Delta Delta G) [Takano et al. (1999) Biochemistry 38, 12698--12708] revealed that the mutant proteins decreased the hydrophobic effect contributing to the stability (Delta G(HP)), but this destabilization was recovered by the hydrogen bonds newly introduced. The present study shows the favorable contribution of polar groups with hydrogen bonds in the interior of protein molecules to the conformational stability.     

Are aromatic carbon donor hydrogen bonds linear in proteins?

Proteins fold and maintain structure through the collective contributions of a large number of weak, noncovalent interactions. The hydrogen bond is one important category of forces that acts on very short distances. As our knowledge of protein structure continues to expand, we are beginning to appreciate the role that weak carbon-donor hydrogen bonds play in structure and function. One property that differentiates hydrogen bonds from other packing forces is propensity for forming a linear donor-hydrogen-acceptor orientation. To ascertain if carbon-donor hydrogen bonds are able to direct acceptor linearity, we surveyed the geometry of interactions specifically involving aromatic sidechain ring carbons in a data set of high resolution protein structures. We found that while donor-acceptor distances for most carbon donor hydrogen bonds were tighter than expected for van der Waals packing, only the carbons of histidine showed a significant bias for linear geometry. By categorizing histidines in the data set into charged and neutral sidechains, we found only the charged subset of histidines participated in linear interactions. B3LYP/6-31G**++ level optimizations of imidazole and indole-water interactions at various fixed angles demonstrates a clear orientation dependence of hydrogen bonding capacity for both charged and neutral sidechains. We suggest that while all aromatic carbons can participate in hydrogen bonding, only charged histidines are able to overcome protein packing forces and enforce linear interactions. The implications for protein modeling and design are discussed.     

Thermodynamic consequences of burial of polar and non-polar amino acid residues in the protein interior

Effects of amino acid substitutions at four fully buried sites of the ubiquitin molecule on the thermodynamic parameters (enthalpy, Gibbs energy) of unfolding were evaluated experimentally using differential scanning calorimetry. The same set of substitutions has been incorporated at each of four sites. These substitutions have been designed to perturb packing (van der Waals) interactions, hydration, and/or hydrogen bonding. From the analysis of the thermodynamic parameters for these ubiquitin variants we conclude that: (i) packing of non-polar groups in the protein interior is favorable and is largely defined by a favorable enthalpy of van der Waals interactions. The removal of one methylene group from the protein interior will destabilize a protein by approximately 5 kJ/mol, and will decrease the enthalpy of a protein by 12 kJ/mol. (ii) Burial of polar groups in the non-polar interior of a protein is highly destabilizing, and the degree of destabilization depends on the relative polarity of this group. For example, burial of Thr side-chain in the non-polar interior will be less destabilizing than burial of Asn side-chain. This decrease in stability is defined by a large enthalpy of dehydration of polar groups upon burial. (iii) The destabilizing effect of dehydration of polar groups upon burial can be compensated if these buried polar groups form hydrogen bonding. The enthalpy of this hydrogen bonding will compensate for the unfavorable dehydration energy and as a result the effect will be energetically neutral or even slightly stabilizing.     

The response of T4 lysozyme to large-to-small substitutions within the core and its relation to the hydrophobic effect.

To further examine the structural and thermodynamic basis of hydrophobic stabilization in proteins, all of the bulky non-polar residues that are buried or largely buried within the core of T4 lysozyme were substituted with alanine. In 25 cases, including eight reported previously, it was possible to determine the crystal structures of the variants. The structures of four variants with double substitutions were also determined. In the majority of cases the "large-to-small" substitutions lead to internal cavities. In other cases declivities or channels open to the surface were formed. In some cases the structural changes were minimal (mainchain shifts < or = 0.3 A); in other cases mainchain atoms moved up to 2 A. In the case of Ile 29 --> Ala the structure collapsed to such a degree that the volume of the putative cavity was zero. Crystallographic analysis suggests that the occupancy of the engineered cavities by solvent is usually low. The mutants Val 149 --> Ala (V149A) and Met 6 --> Ala (M6A), however, are exceptions and have, respectively, one and two well-ordered water molecules within the cavity. The Val 149 --> Ala substitution allows the solvent molecule to hydrogen bond to polar atoms that are occluded in the wild-type molecule. Similarly, the replacement of Met 6 with alanine allows the two solvent molecules to hydrogen bond to each other and to polar atoms on the protein. Except for Val 149 --> Ala the loss of stability of all the cavity mutants can be rationalized as a combination of two terms. The first is a constant for a given class of substitution (e.g., -2.1 kcal/mol for all Leu --> Ala substitutions) and can be considered as the difference between the free energy of transfer of leucine and alanine from solvent to the core of the protein. The second term can be considered as the energy cost of forming the cavity and is consistent with a numerical value of 22 cal mol(-1) A(-3). Physically, this term is due to the loss of van der Waal''s interactions between the bulky sidechain that is removed and the atoms that form the wall of the cavity. The overall results are consistent with the prior rationalization of Leu --> Ala mutants in T4 lysozyme by Eriksson et al. (Eriksson et al., 1992, Science 255:178-183).     

Tolerance to the substitution of buried apolar residues by charged residues in the homologous protein structures

Occurrence and accommodation of charged amino acid residues in proteins that are structurally equivalent to buried non-polar residues in homologues have been investigated. Using a dataset of 1,852 homologous pairs of crystal structures of proteins available at 2A or better resolution, 14,024 examples of apolar residues in the structurally conserved regions replaced by charged residues in homologues have been identified. Out of 2,530 cases of buried apolar residues, 1,677 of the equivalent charged residues in homologues are exposed and the rest of the charged residues are buried. These drastic substitutions are most often observed in homologous protein pairs with low sequence identity (<30%) and in large protein domains (>300 residues). Such buried charged residues in the large proteins are often located in the interface of sub-domains or in the interface of structural repeats, Beyond 7A of residue depth of buried apolar residues, or less than 4% of solvent accessibility, almost all the substituting charged residues are buried. It is also observed that acidic sidechains have higher preference to get buried than the positively charged residues. There is a preference for buried charged residues to get accommodated in the interior by forming hydrogen bonds with another sidechain than the main chain. The sidechains interacting with a buried charged residue are most often located in the structurally conserved regions of the alignment. About 50% of the observations involving hydrogen bond between buried charged sidechain and another sidechain correspond to salt bridges. Among the buried charged residues interacting with the main chain, positively charged sidechains form hydrogen bonds commonly with main chain carbonyls while the negatively charged residues are accommodated by hydrogen bonding with the main chain amides. These carbonyls and amides are usually located in the loops that are structurally variable among homologous proteins.     

Aggregation patterns of bile salts: crystal structure of calcium cholate chloride heptahydrate

Crystals of calcium cholate chloride heptahydrate, CaC24H39O7Cl . 7H2O, are monoclinic, space group P2(1), with a = 11.918(2), b = 8.636(1), c = 15.302(3) A, beta = 97.93(3) degrees, V = 1559.9(8) A3, and Z = 2. A trial structure was obtained by Patterson and Fourier techniques and was refined by full-matrix least-squares calculations using absorption corrected CuK-alpha diffractometer data. The final R index is 0.047. The crystal structure contains bilayer-type arrangements, with hydrophobic portions of cholate rings sandwiched between layers of polar groups that are interacting with calcium ions and water molecules. The calcium ion is coordinated to five water molecules and to the two carboxylate oxygen atoms of the cholate residue. Two additional water molecules are involved only in crystal packing through the formation of hydrogen bonds. Cholate-cholate hydrophobic interactions involve contacts between the hydrocarbon portions of the carboxylate sidechains and the A and B rings. This results in a staggered packing pattern that is nearly identical to that found in crystals of sodium cholate and rubidium deoxycholate. Similar bilayer aggregation patterns may also be involved in the formation of bile salt micelles in aqueous media. The characteristic bilayer packing arrangement can accommodate a variety of cation-binding patterns, as evidenced by the finding that calcium, sodium, and rubidium ions interact with the polar faces of the bilayers in different ways. The carboxylate sidechain displays two different conformations in the crystal structure of calcium cholate chloride heptahydrate. Variation in sidechain conformation may be of importance in the adjustment required to accommodate different cation coordination schemes.     

Structural and thermodynamic analysis of the binding of solvent at internal sites in T4 lysozyme

To investigate the structural and thermodynamic basis of the binding of solvent at internal sites within proteins a number of mutations were constructed in T4 lysozyme. Some of these were designed to introduce new solvent-binding sites. Others were intended to displace solvent from preexisting sites. In one case Val-149 was replaced with alanine, serine, cysteine, threonine, isoleucine, and glycine. Crystallographic analysis shows that, with the exception of isoleucine, each of these substitutions results in the binding of solvent at a polar site that is sterically blocked in the wild-type enzyme. Mutations designed to perturb or displace a solvent molecule present in the native enzyme included the replacement of Thr-152 with alanine, serine, cysteine, valine, and isoleucine. Although the solvent molecule was moved in some cases by up to 1.7 A, in no case was it completely removed from the folded protein. The results suggest that hydrogen bonds from the protein to bound solvent are energy neutral. The binding of solvent to internal sites within proteins also appears to be energy neutral except insofar as the bound solvent may prevent a loss of energy due to potential hydrogen bonding groups that would otherwise be unsatisfied. The introduction of a solvent-binding site appears to require not only a cavity to accommodate the water molecule but also the presence of polar groups to help satisfy its hydrogen-bonding potential. It may be easier to design a site to accommodate two or more water molecules rather than one as the solvent molecules can then hydrogen-bond to each other. For similar reasons it is often difficult to design a point mutation that will displace a single solvent molecule from the core of a protein.     

Replica exchange simulation of reversible folding/unfolding of the Trp-cage miniprotein in explicit solvent: on the structure and possible role of internal water

We simulate the folding/unfolding equilibrium of the 20-residue miniprotein Trp-cage. We use replica exchange molecular dynamics simulations of the AMBER94 atomic detail model of the protein explicitly solvated by water, starting from a completely unfolded configuration. We employ a total of 40 replicas, covering the temperature range between 280 and 538 K. Individual simulation lengths of 100 ns sum up to a total simulation time of about 4 micros. Without any bias, we observe the folding of the protein into the native state with an unfolding-transition temperature of about 440 K. The native state is characterized by a distribution of root mean square distances (RMSD) from the NMR data that peaks at 1.8A, and is as low as 0.4A. We show that equilibration times of about 40 ns are required to yield convergence. A folded configuration in the entire extended ensemble is found to have a lifetime of about 31 ns. In a clamp-like motion, the Trp-cage opens up during thermal denaturation. In line with fluorescence quenching experiments, the Trp-residue sidechain gets hydrated when the protein opens up, roughly doubling the number of water molecules in the first solvation shell. We find the helical propensity of the helical domain of Trp-cage rather well preserved even at very high temperatures. In the folded state, we can identify states with one and two buried internal water molecules interconnecting parts of the Trp-cage molecule by hydrogen bonds. The loss of hydrogen bonds of these buried water molecules in the folded state with increasing temperature is likely to destabilize the folded state at elevated temperatures.     

Free‐energy function based on an all‐atom model for proteins

We have developed a free‐energy function based on an all‐atom model for proteins. It comprises two components, the hydration entropy (HE) and the total dehydration penalty (TDP). Upon a transition to a more compact structure, the number of accessible configurations arising from the translational displacement of water molecules in the system increases, leading to a water‐entropy gain. To fully account for this effect, the HE is calculated using a statistical‐mechanical theory applied to a molecular model for water. The TDP corresponds to the sum of the hydration energy and the protein intramolecular energy when a fully extended structure, which possesses the maximum number of hydrogen bonds with water molecules and no intramolecular hydrogen bonds, is chosen as the standard one. When a donor and an acceptor (e.g., N and O, respectively) are buried in the interior after the break of hydrogen bonds with water molecules, if they form an intramolecular hydrogen bond, no penalty is imposed. When a donor or an acceptor is buried with no intramolecular hydrogen bond formed, an energetic penalty is imposed. We examine all the donors and acceptors for backbone‐backbone, backbone‐side chain, and side chain‐side chain intramolecular hydrogen bonds and calculate the TDP. Our free‐energy function has been tested for three different decoy sets. It is better than any other physics‐based or knowledge‐based potential function in terms of the accuracy in discriminating the native fold from misfolded decoys and the achievement of high Z‐scores. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Significant role of electrostatic interactions for stabilization of protein assemblies

Contribution of electrostatic interactions to stability of BPTI orthorhombic, pig-insulin cubic crystals, and horse L ferritin crystals was evaluated with numerical calculation of Poisson-Boltzmann equation based on a dielectric model. The stability of a ferritin molecule (24-mer) composed of 24 subunits was also evaluated. It was found that the surface charge-charge interactions at separation distances (< 5 Å) were insensitive to variations in the ionic strength, and thus stabilized assembled states of the proteins (i.e., crystalline state and oligomeric state). It was also revealed that the charge density and the packing of the protein crystals were largely responsible for the ionic strength dependence of the crystal stability. The stability of the 5PTI crystalline state with a high charge density drastically increased as the concentration of the solvent ions increased. In contrast, that of the insulin crystal with a low charge density and large solvent region was insensitive to changes in the ionic concentration. The electrostatic interaction between ferritin 24-mers was attributed to two salt bridges mediated by Cd ion. For the stability of the ferritin 24-mer, which is evolutionally designed, the electrostatic stabilization between the subunits was attributed to polar bonds such as buried salt bridges or hydrogen bonds, which occasionally yielded more than 5 kcal/mol and were numerous and very strong compared with the bonds between molecules in the 5PTI and 9INS crystals.By analyzing the atomic charge-charge interactions in detail, it was found that charge pairs separated by less than 3 Å, such as hydrogen bonds, dominantly stabilize the assembled states, and that pairs 3 to 5 Å apart were also important. The stability of the assembled states evaluated by the total EET was determined by the fine balance between the two competing contributions arising from the stabilizing atoms and the destabilizing atoms.Changes of the ASA and hydration free energy were also evaluated in accordance with the process of the subunit assembly. The change of hydration free energy, which was very large (i.e., ~+ 100 kcal/mol/subunit) and unfavorable for the assembly, was proportional to the electrostatic hydration energy (i.e., Born energy change in the hydration process). Hydrophobic groups were likely to appear more frequently than hydrophilic groups at the interfaces.This study offers a method which can improve the stability of protein crystals by introducing polar or charged residues that are properly designed to form specific hydrogen bonds or salt bridges between neighboring protein molecules. This method is also applicable to crystallography, because it improves refinement of protein structures in crystals by taking the inter-protein interactions into account.     

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.     

Why are polar residues within the membrane core evolutionary conserved?

Here, we present a study of polar residues within the membrane core of alpha-helical membrane proteins. As expected, polar residues are less frequent in the membrane than expected. Further, most of these residues are buried within the interior of the protein and are only rarely exposed to lipids. However, the polar groups often border internal water filled cavities, even if the rest of the sidechain is buried. A survey of their functional roles in known structures showed that the polar residues are often directly involved in binding of small compounds, especially in channels and transporters, but other functions including proton transfer, catalysis, and selectivity have also been attributed to these proteins. Among the polar residues histidines often interact with prosthetic groups in photosynthetic- and oxidoreductase-related proteins, whereas prolines often are required for conformational changes of the proteins. Indeed, the polar residues in the membrane core are more conserved than other residues in the core, as well as more conserved than polar residues outside the membrane. The reason is twofold; they are often (i) buried in the interior of the protein and (ii) directly involved in the function of the proteins. Finally, a method to identify which polar residues are present within the membrane core directly from protein sequences was developed. Applying the method to the set of all human membrane proteins the prediction indicates that polar residues were most frequent among active transporter proteins and GPCRs, whereas infrequent in families with few transmembrane regions, such as non-GPCR receptors.     

Helix-capping interaction in λ cro protein: A free energy simulation analysis

The stability mutant Tyr-26 → Asp was studied in the Cro protein from bacteriophage λ using free energy molecular dynamics simulations. The mutant was calculated to be more stable than the wild type by 3.0 ± 1.7 kcal/mol/monomer, in reasonable agreement with experiment (1.4 kcal/mol/monomer). Moreover, the aspartic acid in the mutant was found to form a capping interation with the amino terminus of the third α-helix of Cro. The simulations were analyzed to understand better the source of the stability of this helix-capping interaction and to examine the results in light of previous explanations of stabilizing helix caps-namely, a model of local unsatisfied hydrogen bonds at the helix termini and the helix macro dipole model. Analysis of the simulations shows that the stabilizing effect of this charged helical cap is due both to favorable hydrogen bonds with backbone NH groups at the helix terminus and to favorable electrostatic interactions (but not hydrogen bonds) with their carbonyls (effectively the next row of local dipoles in the helix). However, electrostatic interactions are weak or negligible with backbone dipolar groups in the helix further away from the terminus. Moreover, the importance of other local electrostatic interactions with polar side chains near the helix terminus, which are neglected in most treatments of this effect, are shown to be important. Thus, the results support a model that is intermediate between the two previous explanations: both unsatisfied hydrogen bonds at the helix terminus and other, local preoriented dipolar groups stabilize the helix cap. These findings suggest that similar interactions with preoriented dipolar groups may be important for cooperativity in other charge–dipole interactions and may be employed to advantage for molecular design. © 1994 Wiley-Liss, Inc.