A Systematic Analysis of Buried Polar Side Chains and Their Interactions in α-Helical Proteins

Examination of 80 -helical proteins and domains demonstrates that they contain from 1 to more than 20 completely buried (water-inaccessible) polar side chains. As a rule the latter have partners for H-bonding but the resulting H-bond system is often not saturating. Basing on statistical analysis, we determined the optimal number of H-bonds for every type of polar side chain, and discuss the structural role of vacant donors and acceptors. About half of the H-bonds formed by buried side chains pertain to interhelix contacts of the (side chain)–(side chain) and (side chain)–(main chain) types. Such interactions appear to be a most important factor determining the mutual arrangement of -helices in proteins. Analysis of the frequency of occurrence of various interacting pairs reveals that these interactions are selective.     

Analysis of interactions of buried polar side chains in beta-proteins

It was shown in qualitative and quantitative analyses of polar side chains inaccessible for water molecules as well as their interactions in 100 globular beta-structural proteins that completely buried polar side chains are widespread in beta-proteins, their vast majority being involved in "side chain-side chain" or "side chain-main chain" interactions. The analysis of the occurrence of different "side chain-partner" pairs permitted us to demonstrate that such interactions are selective. The results were compared with similar data obtained previously for alpha-helical proteins.     

Analysis of Interactions of Buried Polar Side Chains in β-Proteins

Qualitative and quantitative analysis of polar side chains inaccessible to water molecules, as well as their interactions in 100 globular β-sheet proteins, was performed. It was shown that completely buried polar side chains are widespread in β-proteins, with their vast majority being involved in side chain-side chain or side chain-main chain interactions. An analysis of frequency of occurrence of different side chain-partner pairs demonstrated that these interactions are selective. The results were compared with similar data obtained earlier for α-helical proteins.     

Polar residues in membrane domains of proteins: molecular basis for helix-helix association in a mutant CFTR transmembrane segment

Polar side chains constitute over 20% of residues in the transmembrane (TM) helices of membrane proteins, where they may serve as hydrogen bond interaction sites for phenotypic polar mutations that arise in membrane protein-related diseases. To systematically explore the structural consequences of H-bonds between TM helices, we focused on TM4 of the cystic fibrosis conductance regulator (CFTR) and its cystic fibrosis- (CF-) phenotypic mutation, V232D, as a model system. Synthetic peptides corresponding to wild-type (TM4-wt) (residues 219-242: LQASAFCGLGFLIVLALFQAGLGR) and mutant (TM4-V232D) sequences both adopt helical structures in SDS micelles and display dimer bands on SDS-PAGE arising from disulfide bond formation via wild-type residue Cys-225. However, the TM4-V232D peptide additionally forms a ladder of noncovalent oligomers, including tetramers, hexamers, and octamers, mediated by a hydrogen bond network involving Asp-Gln side chain-side chain interactions. Ala-scanning mutagenesis of the TM4 sequence indicated that ladder formation minimally required the simultaneous presence of the Cys-225, Asp-232, and Gln-237 residues. As random hydrophobic sequences containing these three residues at TM4 equivalent positions did not oligomerize, specific van der Waals packing interactions between helix side chains were also shown to play a crucial role. Overall, the results suggest that polar mutations in membrane domains, in conjunction with critically positioned polar partner residues, potentially constitute a source of aberrant helix interactions that could contribute to loss of function when they arise in protein transmembrane domains.     

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.     

Position-dependence of stabilizing polar interactions of asparagine in transmembrane helical bundles

Recent studies with model peptides and statistical analyses of the crystal structures of membrane proteins have shown that buried polar interactions contribute significantly to the stabilization of the three-dimensional structures of membrane proteins. Here, we probe how the location of these polar groups along the transmembrane helices affect their free energies of interaction. Asn residues were placed singly and in pairs at three positions within a model transmembrane helix, which had previously been shown to support the formation of trimers in micelles. The model helix was designed to form a transmembrane coiled coil, with Val side chains at the "a" positions of the heptad repeat. Variants of this peptide were prepared in which an Asn residue was introduced at one or more of the "a" positions, and their free energies of association were determined by analytical ultracentrifugation. When placed near the middle of the transmembrane helix, the formation of trimers was stabilized by at least -2.0 kcal/mol per Asn side chain. When the Asn was placed at the interface between the hydrophobic and polar regions of the peptide, the substitution was neither stabilizing nor destabilizing (0.0 +/- 0.5 kcal/mol of monomer). Finally, it has previously been shown that a Val-for-Asn mutation in a water-soluble coiled coil destabilizes the structure by approximately 1.5 kcal/mol of monomer [Acharya, A., et al. (2002) Biochemistry 41, 14122-14131]. Thus, the headgroup region of a micelle appears to have a conformational impact intermediate between that of bulk water and the apolar region of micelle. A similarly large dependence on the location of the polar residues was found in a statistical survey of helical transmembrane proteins. The tendency of different types of residues to be buried in the interiors versus being exposed to lipids was analyzed. Asn and Gln show a very strong tendency to be buried when they are located near the middle of a transmembrane helix. However, when placed near the ends of transmembrane helices, they show little preference for the surface versus the interior of the protein. These data show that Asn side chains within the apolar region of the transmembrane helix provide a significantly larger driving force for association than Asn residues near the apolar/polar interface. Thus, although polar interactions are able to strongly stabilize the folding of membrane proteins, the energetics of association depend on their location within the hydrophobic region of a transmembrane helix.     

Statistical analysis of atomic contacts at RNA-protein interfaces

Forty-five crystals of complexes between proteins and RNA molecules from the Protein Data Bank have been statistically surveyed for the number of contacts between RNA components (phosphate, ribose and the four bases) and amino acid side chains. Three groups of complexes were defined: the tRNA synthetases; the ribosomal complexes; and a third group containing a variety of complexes. The types of atomic contacts were a priori classified into ionic, neutral H-bond, C-H...O H-bond, or van der Waals interaction. All the contacts were organized into a relational database which allows for statistical analysis. The main conclusions are the following: (i) in all three groups of complexes, the most preferred amino acids (Arg, Asn, Ser, Lys) and the less preferred ones (Ala, Ile, Leu, Val) are the same; Trp and Cys are rarely observed (respectively 15 and 5 amino acids in the ensemble of interfaces); (ii) of the total number of amino acids located at the interfaces 22% are hydrophobic, 40% charged (positive 32%, negative 8%), 30% polar and 8% are Gly; (iii) in ribosomal complexes, phosphate is preferred over ribose, which is preferred over the bases, but there is no significant preference in the other two groups; (iv) there is no significant prevalence of a base type at protein-RNA interfaces, but specifically Arg and Lys display a preference for phosphate over ribose and bases; Pro and Asn prefer bases over ribose and phosphate; Met, Phe and Tyr prefer ribose over phosphate and bases. Further, Ile, Pro, Ser prefer A over the others; Leu prefers C; Asp and Gly prefer G; and Asn prefers U. Considering the contact types, the following conclusions could be drawn: (i) 23% of the contacts are via potential H-bonds (including CH...O H-bonds and ionic interactions), 72% belong to van der Waals interactions and 5% are considered as short contacts; (ii) of all potential H-bonds, 54% are standard, 33% are of the C-H...O type and 13% are ionic; (iii) the Watson-Crick sites of G, O6(G) and principally N2(G) and the hydroxyl group O2' is more often involved in H-bonds than expected; the protein main chain is involved in 32% and the side chains in 68% of the H-bonds; considering the neutral and ionic H-bonds, the following couples are more frequent than expected-base A-Ser, base G-Asp/Glu, base U-Asn. The RNA CH groups interact preferentially with oxygen atoms (62% on the main chain and 19% on the side chains); (iv) the bases are involved in 38% of all H-bonds and more than 26% of the H-bonds have the H donor group on the RNA; (v) the atom O2' is involved in 21% of all H-bonds, a number greater than expected; (vi) amino acids less frequently in direct contact with RNA components interact frequently via their main chain atoms through water molecules with RNA atoms; in contrast, those frequently observed in direct contact, except Ser, use instead their side chain atoms for water bridging interactions.     

Use of a potential of mean force to analyze free energy contributions in protein folding.

A method for calculation of the free energy of residues as a function of residue burial is proposed. The method is based on the potential of mean force, with a reaction coordinate expressed by residue burial. Residue burials are calculated from high-resolution protein structures. The largest individual contributions to the free energy of a residue are found to be due to the hydrophobic interactions of the nonpolar atoms, interactions of the main chain polar atoms, and interactions of the charged groups of residues Arg and Lys. The contribution to the free energy of folding due to the uncharged side chain polar atoms is small. The contribution to the free energy of folding due to the main chain polar atoms is favorable for partially buried residues and less favorable or unfavorable for fully buried residues. Comparison of the accessible surface areas of proteins and model spheres shows that proteins deviate considerably from a spherical shape and that the deviations increase with the size of a protein. The implications of these results for protein folding are also discussed.     

Alpha helix capping in synthetic model peptides by reciprocal side chain–main chain interactions: Evidence for an N terminal “capping box”

A significant fraction of the amino acids in proteins are alpha helical in conformation. Alpha helices in globular proteins are short, with an average length of about twelve residues, so that residues at the ends of helices make up an important fraction of all helical residues. In the middle of a helix, H-bonds connect the NH and CO groups of each residue to partners four residues along the chain. At the ends of a helix, the H-bond potential of the main chain remains unfulfilled, and helix capping interactions involving bonds from polar side chains to the NH or CO of the backbone have been proposed and detected. In a study of synthetic helical peptides, we have found that the sequence Ser-Glu-Asp-Glu stabilizes the alpha helix in a series of helical peptides with consensus sequences. Following the report by Harper and Rose, which identifies SerXaaXaaGlu as a member of a class of common motifs at the N termini of alpha helices in proteins that they refer to as “capping boxes,” we have reexamined the side chain–main chain interactions in a varient sequence using ¹H NMR, and find that the postulated reciprocal side chain-backbone bonding between the first Ser and last Glu side chains and their peptide NH partners can be resolved: Deletion of two residues N terminal to the Ser-Glu-Asp-Glu sequence in these peptides has no effect on the initiation of helical structure, as defined by two-dimensional (2D) NMR experiments on this variant. Thus the capping box sequence Ser-Glu-Asp-Glu inhibits N terminal fraying of the N terminus of alpha helix in these peptides, and shows the side chain–main chain interactions proposed by Harper and Rose. It thus acts as a helix initiating signal. Since normal a helix cannot propagate beyond the N terminus of this structure, the box acts as a termination signal in this direction as well. © 1994 John Wiley & Sons, Inc.     

Side chain dynamics and alternative hydrogen bonding in the mechanism of protein thermostabilization

To elucidate the mechanism of protein thermostabilization, the thermodynamic properties of small monomeric proteins from mesophilic and thermophilic organisms have been analyzed. Molecular dynamics simulations were employed in the study of dynamic features of charged and polar side chains of amino acid residues. The basic conclusion has been made: surface charged and polar side chains with high conformational mobility can form alternative hydrogen bonded (H-bonded) donor-acceptor pairs. The correlation between the quantitative content of alternative H-bonds per residue and the temperature of maximal thermostability of proteins has been found. The proposed mechanism of protein thermostabilization suggests continuous disruption of the primary H-bonds and formation of alternative ones, which maintain constant the enthalpy value in the native state and prevent a rapid increase of the conformational entropy with the rising temperature. The analysis of the results show that the more residues located in the N- and C-terminal regions and in the extended loops that are capable of forming alternative longer-range H-bonded pairs, the higher the protein thermostability.     

Asparagine-mediated self-association of a model transmembrane helix

In membrane proteins, the extent to which polarity, hydrogen bonding, and van der Waals packing interactions of the buried, internal residues direct protein folding and association of transmembrane segments is poorly understood. The energetics associated with these various interactions should differ substantially between membrane versus water-soluble proteins. To help evaluate these energetics, we have altered a water-soluble, two-stranded coiled-coil peptide to render its sequence soluble in membranes. The membrane-soluble peptide associates in a monomer-dimer-trimer equilibrium, in which the trimer predominates at the highest peptide/detergent ratios. The oligomers are stabilized by a buried Asn side chain. Mutation of this Asn to Val essentially eliminates oligomerization of the membrane-soluble peptide. Thus, within a membrane-like environment, interactions involving a polar Asn side chain provide a strong thermodynamic driving force for membrane helix association.     

Stabilization of alpha-helix structure by polar side-chain interactions: complex salt bridges,cation-pi interactions,and C-H em leader O H-bonds

It is generally understood that helical proteins are stabilized by a combination of hydrophobic and packing interactions, together with H-bonds and electrostatic interactions. Here we show that polar side-chain interactions on the surface can play an important role in helix formation and stability. We review studies on model helical peptides that reveal the effect of weak interactions between side chains on helix stability, focusing on some nonclassical side-chain-side-chain interactions: complex salt bridges, cation-pi, and C-H em leader O H-bonding interactions. Each of these can be shown to contribute to helix stability, and thus must be included in a comprehensive catalogue of helix stabilizing effects. The issue of the structure of the unfolded states of helical peptides is also discussed, in the light of recent experiments showing that these contain substantial amounts of polyproline II conformation.     

Relationship between intramolecular hydrogen bonding and solvent accessibility of side-chain donors and acceptors in proteins

This study shows that intramolecular hydrogen bonding in proteins depends on the accessibility of donors and acceptors to water molecules. The frequency of occurrence of H-bonded side chains in proteins is inversely proportional to the solvent accessibility of their donors and acceptors. Estimates of the notional free energy of hydrogen bonding suggest that intramolecular hydrogen-bonding interactions of buried and half-buried donors and acceptors can contribute favorably to the stability of a protein, whereas those of solvent-exposed polar atoms become less favorable or unfavorable.     

Analysis of tryptophan surface accessibility in proteins by MALDI-TOF mass spectrometry

Surface accessible amino acids can play an important role in proteins. They can participate in enzyme's active center structure or in specific intermolecular interactions. Thus, the information about selected amino acids' surface accessibility can contribute to the understanding of protein structure and function. In this paper, we present a simple method for surface accessibility mapping of tryptophan side chains by their chemical modification and identification by MALDI-TOF mass spectrometry. The reaction with 2-hydroxy-5-nitrobenzyl bromide, a common and highly specific covalent modification of tryptophan, seems to be very useful for this purpose. The method was tested on four model proteins with known spatial structure. In the native proteins (1) only surface accessible tryptophan side chains were found to react with the modification agent and (2) no buried one was found to react at lower reagent concentrations. These results indicate that the described method can be a potent tool for identification of surface-located tryptophan side chain in a protein of unknown conformation.     

Non-native interhelical hydrogen bonds in the cystic fibrosis transmembrane conductance regulator domain modulated by polar mutations

Polar residues comprise about 15% of the transmembrane (TM) domains of proteins, where they can stabilize structure via native side chain-side chain interhelical hydrogen bonds between TM helices. However, non-native H-bonds may be implicated in disease states, through limiting protein dynamics during transport and/or misfolding the protein by inducing non-native rotational positions about TM helical axes. Here we have undertaken an investigation of the presence and strength of H-bond interactions within a series of helix-loop-helix ("hairpin") constructs derived from TM helices 3 and 4 (italic) of the cystic fibrosis transmembrane conductance regulator (CFTR) (prototypic sequence G(194)LALAHFVWIAPLQ(207)VALLMGLIWELLQASAFAGLGFLIV(232)LALFQ(237)AGLG(241)) in which wild-type Q207 in TM3 forms an interhelical H-bond with CF-phenotypic mutant V232D in TM4 [Therien, A. G., Grant, F. E., and Deber, C. M. (2001) Nat. Struct. Biol 8, 597-601]. In the present work, a library of 21 TM3/4 constructs was prepared, where Asp residues were placed individually at TM4 positions 221-241. Using gel shift assays-in which H-bond-linked hairpins (closed conformation) migrate faster than the elongated forms (open conformation)-we found that Q207 in TM3 is able to "capture" all 21 TM4 D mutations into measurable populations of interhelical H-bonds. A similar library of TM4 D mutants-but also containing Q207L-reverted to wild-type migration rates, confirming Q207 as the polar partner for TM4 D residues. In view of the broad capture range of Q207, these results emphasize the potential consequences to folding and dynamics of introducing polar mutations into the TM domains of membrane proteins in the vicinity of a native polar TM residue.     

Structural and conformational changes concomitant with the E1-E2 transition in H(+)K(+)-ATPase: a comparative protein modeling study

Comparative modeling studies on conserved regions of the gastric H(+)K(+)-ATPase reveal that the E1-E2 conformational transition induces significant tertiary structural changes while conserving the secondary structure. The residues 516-530 of the cytoplasmic domain and TM10 within the transmembrane (TM) regions undergo maximum tertiary structural changes. The luminal regions exhibit comparatively lesser tertiary structural deviations. Residues 249-304 show maximum secondary structural deviation in the conformational transition. The Cys-815 and Cys-323 residues involved in inhibitor binding are found to have smaller buried side chain areas in the E1 conformation compared to E2. Retention of activity correlates well with the buried side chain area when selected amino acid residues in TM6 are mutated using modeling techniques with bulkier amino acid residues. Conformational specificity for ion binding is corroborated with the fraction of side chains exposed to polar atoms of the residues E345, D826, V340, A341, V343, and E822.     

Lipid-protein interactions of integral membrane proteins: a comparative simulation study

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.     

Helix propensities of the amino acids measured in alanine-based peptides without helix-stabilizing side-chain interactions.

Helix propensities of the amino acids have been measured in alanine-based peptides in the absence of helix-stabilizing side-chain interactions. Fifty-eight peptides have been studied. A modified form of the Lifson-Roig theory for the helix-coil transition, which includes helix capping (Doig AJ, Chakrabartty A, Klingler TM, Baldwin RL, 1994, Biochemistry 33:3396-3403), was used to analyze the results. Substitutions were made at various positions of homologous helical peptides. Helix-capping interactions were found to contribute to helix stability, even when the substitution site was not at the end of the peptide. Analysis of our data with the original Lifson-Roig theory, which neglects capping effects, does not produce as good a fit to the experimental data as does analysis with the modified Lifson-Roig theory. At 0 degrees C, Ala is a strong helix former, Leu and Arg are helix-indifferent, and all other amino acids are helix breakers of varying severity. Because Ala has a small side chain that cannot interact significantly with other side chains, helix formation by Ala is stabilized predominantly by the backbone ("peptide H-bonds"). The implication for protein folding is that formation of peptide H-bonds can largely offset the unfavorable entropy change caused by fixing the peptide backbone. The helix propensities of most amino acids oppose folding; consequently, the majority of isolated helices derived from proteins are unstable, unless specific side-chain interactions stabilize them.     

Extension of the fragment method to calculate amino acid zwitterion and side chain partition coefficients

The fragment method of calculating partition coefficients (P) has been extended to include the common amino acids (AAs). The results indicate that polar and charged side chains influence the hydrophobicity of atoms in the side chain in a predictable manner. Field effects, as evidenced through polar proximity factors and bond factors, need to be considered for accurate estimation of transfer phenomena. The calculated log P and delta G degree ' values of the 20 AAs agree well with the observed values. Pro calculates to be more hydrophilic than the observed log P. Hydrophobicity scales for peptide side chain residues are compared and evaluated in terms of suitability. Calculated pi values for nonpolar side chain residues agree well with the observed values; calculated values for uncharged polar side chain residues deviate by about 0.6 log units except for Gln and Cys; and polar side chain residues with charged side chains calculate as too hydrophilic. Reasons for the differences are explored. We also suggest that tightly bound water to polar moieties in amino acids and peptides may be transferred into the octanol phase during partitioning experiments. A quantitative methodology is presented which characterizes the thermodynamic partitioning of groups and individual atoms in amino acids and proteins.     

One- and two-body decomposable Poisson-Boltzmann methods for protein design calculations

Successfully modeling electrostatic interactions is one of the key factors required for the computational design of proteins with desired physical, chemical, and biological properties. In this paper, we present formulations of the finite difference Poisson-Boltzmann (FDPB) model that are pairwise decomposable by side chain. These methods use reduced representations of the protein structure based on the backbone and one or two side chains in order to approximate the dielectric environment in and around the protein. For the desolvation of polar side chains, the two-body model has a 0.64 kcal/mol RMSD compared to FDPB calculations performed using the full representation of the protein structure. Screened Coulombic interaction energies between side chains are approximated with an RMSD of 0.13 kcal/mol. The methods presented here are compatible with the computational demands of protein design calculations and produce energies that are very similar to the results of traditional FDPB calculations.