Comparison of helix interactions in membrane and soluble alpha-bundle proteins

Helix-helix interactions are important for the folding, stability, and function of membrane proteins. Here, two independent and complementary methods are used to investigate the nature and distribution of amino acids that mediate helix-helix interactions in membrane and soluble alpha-bundle proteins. The first method characterizes the packing density of individual amino acids in helical proteins based on the van der Waals surface area occluded by surrounding atoms. We have recently used this method to show that transmembrane helices pack more tightly, on average, than helices in soluble proteins. These studies are extended here to characterize the packing of interfacial and noninterfacial amino acids and the packing of amino acids in the interfaces of helices that have either right- or left-handed crossing angles, and either parallel or antiparallel orientations. We show that the most abundant tightly packed interfacial residues in membrane proteins are Gly, Ala, and Ser, and that helices with left-handed crossing angles are more tightly packed on average than helices with right-handed crossing angles. The second method used to characterize helix-helix interactions involves the use of helix contact plots. We find that helices in membrane proteins exhibit a broader distribution of interhelical contacts than helices in soluble proteins. Both helical membrane and soluble proteins make use of a general motif for helix interactions that relies mainly on four residues (Leu, Ala, Ile, Val) to mediate helix interactions in a fashion characteristic of left-handed helical coiled coils. However, a second motif for mediating helix interactions is revealed by the high occurrence and high average packing values of small and polar residues (Ala, Gly, Ser, Thr) in the helix interfaces of membrane proteins. Finally, we show that there is a strong linear correlation between the occurrence of residues in helix-helix interfaces and their packing values, and discuss these results with respect to membrane protein structure prediction and membrane protein stability.     

Interactions between two beta-sheets. Energetics of beta/beta packing in proteins

The analysis of the interactions between regularly folded segments of the polypeptide chain contributes to an understanding of the energetics of protein folding. Conformational energy-minimization calculations have been carried out to determine the favorable ways of packing two right-twisted beta-sheets. The packing of two five-stranded beta-sheets was investigated, with the strands having the composition CH3CO-(L-Ile)6-NHCH3 in one beta-sheet and CH3CO-(L-Val)6-NHCH3 in the other. Two distinct classes of low-energy packing arrangements were found. In the class with lowest energies, the strands of the two beta-sheets are aligned nearly parallel (or antiparallel) with each other, with a preference for a negative orientation angle, because this arrangement corresponds to the best complementary packing of the two twisted saddle-shaped beta-sheets. In the second class, with higher interaction energies, the strands of the two beta-sheets are oriented nearly perpendicular to each other. While the surfaces of the two beta-sheets are not complementary in this arrangement, there is good packing between the corner of one beta-sheet and the interior part of the surface of the other, resulting in a favorable energy of packing. Both classes correspond to frequently observed orientations of beta-sheets in proteins. In proteins, the second class of packing is usually observed when the two beta-sheets are covalently linked, i.e. when a polypeptide strand passes from one beta-sheet to the other, but we have shown here that a large contribution to the stabilization of this packing arrangement arises from noncovalent interactions.     

Shorter side chains optimize helix-helix packing

A systematic study of helix-helix packing in a comprehensive database of protein structures revealed that the side chains inside helix-helix interfaces on average are shorter than those in the noninterface parts of the helices. The study follows our earlier study of this effect in transmembrane helices. The results obtained on the entire database of protein structures are consistent with those obtained on the transmembrane helices. The difference in the length of interface and noninterface side chains is small but statistically significant. It indicates that helices, if viewed along their main axis, statistically are not circular, but have a flattened interface. This effect brings the helices closer to each other and creates a tighter structural packing. The results provide an interesting insight into the aspects of protein structure and folding.     

Rate of beta-structure formation in polypeptides

An explanation is suggested for why a marginally stable beta-structure folds extremely slowly; it is predicted that even a small increase in stability drastically accelerates beta-folding. According to the theory, this folding is a first-order phase transition, and the rate-limiting step is nucleation. The rate-determining "nucleus" (transition state) is the smallest beta-sheet that is sufficiently large to provide an overall free energy reduction during subsequent folding. If the stability of the beta-structure is low, the nucleus is large and possesses a high free energy due to having a large perimeter. When the net stability of the final beta-structure increases (due to either an increase of the beta-sheet stability or a decrease in stability of the competing structures, e.g., alpha-helices), the size and energy of a nucleus decrease and the rate of folding increases exponentially. This must result in a fast folding of polypeptides enriched by beta-forming residues (e.g., protein chains). The theory is developed for intramolecular beta-structure, but it can also explain the overall features of intermolecular beta-folding; it is applicable both to antiparallel and parallel beta-sheets. The difference in folding of beta-sheets, alpha-helices, and proteins is discussed.     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

Sequence specificity in the dimerization of transmembrane alpha-helices.

While several reports have suggested a role for helix-helix interactions in membrane protein oligomerization, there are few direct biochemical data bearing on this subject. Here, using mutational analysis, we show that dimerization of the transmembrane alpha-helix of glycophorin A in a detergent environment is spontaneous and highly specific. Very subtle changes in the side-chain structure at certain sensitive positions disrupt the helix-helix association. These sensitive positions occur at approximately every 3.9 residues along the helix, consistent with their comprising the interface of a closely fit transmembranous supercoil of alpha-helices. By contrast with other reported cases of interactions between transmembrane helices, the set of interfacial residues in this case contains no highly polar groups. Amino acids with aliphatic side chains define much of the interface, indicating that precise packing interactions between the helices may provide much of the energy for association. These data highlight the potential general importance of specific interactions between the hydrophobic anchors of integral membrane proteins.     

Self-association of transmembrane alpha-helices in model membranes: importance of helix orientation and role of hydrophobic mismatch

Interactions between transmembrane helices play a key role in almost all cellular processes involving membrane proteins. We have investigated helix-helix interactions in lipid bilayers with synthetic tryptophan-flanked peptides that mimic the membrane spanning parts of membrane proteins. The peptides were functionalized with pyrene to allow the self-association of the helices to be monitored by pyrene fluorescence and Trp-pyrene fluorescence resonance energy transfer (FRET). Specific labeling of peptides at either their N or C terminus has shown that helix-helix association occurs almost exclusively between antiparallel helices. Furthermore, computer modeling suggested that antiparallel association arises primarily from the electrostatic interactions between alpha-helix backbone atoms. We propose that such interactions may provide a force for the preferentially antiparallel association of helices in polytopic membrane proteins. Helix-helix association was also found to depend on the lipid environment. In bilayers of dioleoylphosphatidylcholine, in which the hydrophobic length of the peptides approximately matched the bilayer thickness, association between the helices was found to require peptide/lipid ratios exceeding 1/25. Self-association of the helices was promoted by either increasing or decreasing the bilayer thickness, and by adding cholesterol. These results indicate that helix-helix association in membrane proteins can be promoted by unfavorable protein-lipid interactions.     

Modeling microdomains: the surface area of globin helices.

The accessible surface areas of 53 high-resolution globin helices are correlated with molecular weight. The linear fit is assessed for statistical accuracy using a boot-strap analysis, and by comparison to the areas of 13 ideal polyalanine alpha-helices. The accessible area of the unfolded helices is compared with the folded values before helix-helix packing. An analytical physical model is presented to explain the correlation, and to provide an analytical value for the surface area parameter in the diffusion-collision model of protein folding.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

How do helix–helix interactions help determine the folds of membrane proteins? Perspectives from the study of homo‐oligomeric helical bundles

The final, structure-determining step in the folding of membrane proteins involves the coalescence of preformed transmembrane helices to form the native tertiary structure. Here, we review recent studies on small peptide and protein systems that are providing quantitative data on the interactions that drive this process. Gel electrophoresis, analytical ultracentrifugation, and fluorescence resonance energy transfer (FRET) are useful methods for examining the assembly of homo-oligomeric transmembrane helical proteins. These methods have been used to study the assembly of the M2 proton channel from influenza A virus, glycophorin, phospholamban, and several designed membrane proteins-all of which have a single transmembrane helix that is sufficient for association into a transmembrane helical bundle. These systems are being studied to determine the relative thermodynamic contributions of van der Waals interactions, conformational entropy, and polar interactions in the stabilization of membrane proteins. Although the database of thermodynamic information is not yet large, a few generalities are beginning to emerge concerning the energetic differences between membrane and water-soluble proteins: the packing of apolar side chains in the interior of helical membrane proteins plays a smaller, but nevertheless significant, role in stabilizing their structure. Polar, hydrogen-bonded interactions occur less frequently, but, nevertheless, they often provide a strong driving force for folding helix-helix pairs in membrane proteins. These studies are laying the groundwork for the design of sequence motifs that dictate the association of membrane helices.     

Specific sequence combinations at parallel and antiparallel helix-helix interfaces

Orientation of helices at parallel and antiparallel helix-helix interfaces in proteins depends on interacting amino acids from both helices. Particularly important are amino acids at positions analogous to a and d in GCN4 leucine zipper nomenclature, which form hydrophobic core. In this work repeating sequence combinations at a and d positions characteristic for both parallel and antiparallel packing are shown. Layer packing of hydrophobic groups is compared for possible combinations of aliphatic amino acids at all four positions. Correlation between specific position of methyl groups and interhelical angle is found for parallel and antiparallel types of packing.     

GXXXG and AXXXA: common alpha-helical interaction motifs in proteins, particularly in extremophiles

The GXXXG motif is a frequently occurring sequence of residues that is known to favor helix-helix interactions in membrane proteins. Here we show that the GXXXG motif is also prevalent in soluble proteins whose structures have been determined. Some 152 proteins from a non-redundant PDB set contain at least one alpha-helix with the GXXXG motif, 41 +/- 9% more than expected if glycine residues were uniformly distributed in those alpha-helices. More than 50% of the GXXXG-containing alpha-helices participate in helix-helix interactions. In fact, 26 of those helix-helix interactions are structurally similar to the helix-helix interaction of the glycophorin A dimer, where two transmembrane helices associate to form a dimer stabilized by the GXXXG motif. As for the glycophorin A structure, we find backbone-to-backbone atomic contacts of the C alpha-H...O type in each of these 26 helix-helix interactions that display the stereochemical hallmarks of hydrogen bond formation. These glycophorin A-like helix-helix interactions are enriched in the general set of helix-helix interactions containing the GXXXG motif, suggesting that the inferred C alpha-H...O hydrogen bonds stabilize the helix-helix interactions. In addition to the GXXXG motif, some 808 proteins from the non-redundant PDB set contain at least one alpha-helix with the AXXXA motif (30 +/- 3% greater than expected). Both the GXXXG and AXXXA motifs occur frequently in predicted alpha-helices from 24 fully sequenced genomes. Occurrence of the AXXXA motif is enhanced to a greater extent in thermophiles than in mesophiles, suggesting that helical interaction based on the AXXXA motif may be a common mechanism of thermostability in protein structures. We conclude that the GXXXG sequence motif stabilizes helix-helix interactions in proteins, and that the AXXXA sequence motif also stabilizes the folded state of proteins.     

Helix Interaction Tool (HIT): a web-based tool for analysis of helix-helix interactions in proteins

MOTIVATION: In many proteins, helix-helix interactions can be critical to establishing protein conformation (folding) and dynamics, as well as determining associations between protein units. However, the determination of a set of rules that guide helix-helix interaction has been elusive. In order to gain further insight into the helix-helix interface, we have developed a comprehensive package of tools for analyzing helix-helix packing in proteins. These tools are available at http://helix.gersteinlab.org. They include quantitative measures of the helix interaction surface area and helix crossing angle, as well as several methods for visualizing the helical interaction. These methods can be used for analysis of individual protein conformations or to gain insight into dynamic changes in helix interactions. For the latter purpose, a direct interface from entries in the Molecular Motions Database to the HIT site has been provided.     

Secondary structure length as a determinant of folding rate of proteins with two- and three-state kinetics

We present a simple method for determining the folding rates of two- and three-state proteins from the number of residues in their secondary structures (secondary structure length). The method is based on the hypothesis that two- and three-state foldings share a common pattern. Three-state proteins first condense into metastable intermediates, subsequent forming of alpha-helices, turns, and beta-sheets at slow rate-limiting step. The folding rate of such proteins anticorrelate with the length of these beta-secondary structures. It is also assumed that in two-state folding, rapidly folded alpha-helices and turns may facilitate formation of fleeting unobservable intermediates and thus show two-state behavior. There is an inverse relationship between the folding rate and the length of beta-sheets and loops. Our study achieves 94.0 and 88.1% correlations with folding rates determined experimentally for 21 three- and 38 two-state proteins, respectively, suggesting that protein-folding rates are determined by the secondary structure length. The kinetic kinds are selected on the basis of a competitive formation of hydrophobic collapse and alpha-structure in early intermediates.     

Amino acid composition of parallel helix-helix interfaces

Amino acids at helix-helix parallel interfaces influence arrangement of helices and interhelical angles. Parallel interfaces in 79 proteins were considered. Location of amino acids at the positions analogous to a and d in GCN4 leucine zipper nomenclature shows that certain combinations of amino acids characteristic for parallel packing occur more often than could be expected by chance. Repeating sequence combinations occur at a and d positions of parallel helix-helix interfaces with similar values of interhelical angles not only in homologous proteins but also within the same protein and in nonhomologous proteins. Within each group of observed combinations correlation exists between the size of amino acid and magnitude of the interhelical angle.     

Role of stabilization centers in 4 helix bundle proteins

Stabilization centers (SCs) were shown to play an important role in preventing decay of three-dimensional protein structures. These residue clusters, stabilized by cooperative long-range interactions, were proposed to serve as anchoring points for arranging secondary structure elements. In all-alpha proteins, SC elements appear less frequently than in all-beta, alpha/beta, and alpha+beta proteins suggesting that tertiary structure formation of all-alpha proteins is governed by different principles than in other protein classes. Here we analyzed the relation between the formation of stabilization centers and the inter-axial angles (Omega) of alpha-helices in 4 helix bundle proteins. In the distance range, where dipoles have dominant effect on the helix pair arrangement, those helix pairs, where residues from both helices participate in SC elements, appear as parallel more frequently than those helices where no SC elements are present. For SC containing helix pairs, the energetic difference between the parallel and anti-parallel states decreases considerably from 1.1 kcal/mol to 0.4 kcal/mol. Although the observed effect is weak for more distant helices, a competition between the SC element formation and the optimal dipole-dipole interaction of alpha-helices is proposed as a mechanism for tertiary structure formation in 4 helix bundle proteins. The SC-forming potential of different arrangements as well as the pitfalls of the SC definition are also discussed.     

Size and folding in globular proteins

We have modeled protein folding by packing a unified length of regular structural elements (alpha-helices and beta-sheets) into a 'cube'. In a globular protein with m alpha-helices and n beta-strands, this unified length is expressed in units of heptapeptides in alpha-helices, and in units of tripeptides in beta-strands. Calculations using published data show that a 4-helix bundle (m = 4, n = 0) has at least 2 x 2 x 2 helical heptapeptides; the 16-strand beta-barrel of porin (m = 0, n = 16) is at most 4 x 4 x 4 tripeptides in beta-strands. Compact, recurring protein modules with mixed helices and beta-strands are the ones that actually acquire a geometrically quasi-spherical, or cubic, shape.     

Helix packing in polytopic membrane proteins: role of glycine in transmembrane helix association.

The nature and distribution of amino acids in the helix interfaces of four polytopic membrane proteins (cytochrome c oxidase, bacteriorhodopsin, the photosynthetic reaction center of Rhodobacter sphaeroides, and the potassium channel of Streptomyces lividans) are studied to address the role of glycine in transmembrane helix packing. In contrast to soluble proteins where glycine is a noted helix breaker, the backbone dihedral angles of glycine in transmembrane helices largely fall in the standard alpha-helical region of a Ramachandran plot. An analysis of helix packing reveals that glycine residues in the transmembrane region of these proteins are predominantly oriented toward helix-helix interfaces and have a high occurrence at helix crossing points. Moreover, packing voids are generally not formed at the position of glycine in folded protein structures. This suggests that transmembrane glycine residues mediate helix-helix interactions in polytopic membrane proteins in a fashion similar to that seen in oligomers of membrane proteins with single membrane-spanning helices. The picture that emerges is one where glycine residues serve as molecular notches for orienting multiple helices in a folded protein complex.     

Distinct protein interfaces in transmembrane domains suggest an in vivo folding model

Given the known high-resolution structures of alpha-helical transmembrane domains, we show that there are statistically distinct classes of transmembrane interfaces which relate to the folding and oligomerization of transmembrane domains. Distinct types of interfaces have been categorized and refer to those between: the same polypeptide chain, different polypeptide chains, helices that are sequential neighbors, and those that are nonsequential. These different interfaces may reflect different phases in the mechanism of transmembrane domain folding and are consistent with the current experimental evidence pertaining to the folding and oligomerization of transmembrane domains. The classes of helix-helix interfaces have been identified in terms of the numbers and different types of pairwise amino acid interactions. The specific measures used are interaction entropy, the information content of interacting partners compared to a random set of contacts, the amino acid composition of the classes and the abundances of specific amino acid pairs in close contact. Knowledge of the clear differences in the types of helix-helix contacts helps with the derivation of knowledge-based constraints which until now have focused on only the interiors of transmembrane domains as compared to the exterior. Taken together, an in vivo model for membrane protein folding is presented, which is distinct from the familiar two-stage model. The model takes into account the different interfaces of membrane helices defined herein, and the available data regarding folding in the translocation channel.     

Flexibility of beta-sheets: principal component analysis of database protein structures

Protein folds are built primarily from the packing together of two types of structures: alpha-helices and beta-sheets. Neither structure is rigid, and the flexibility of helices and sheets is often important in determining the final fold (e.g., coiled coils and beta-barrels). Recent work has quantified the flexibility of alpha-helices using a principal component analysis (PCA) of database helical structures (J. Mol. Bio. 2003, 327, pp. 229-237). Here, we extend the analysis to beta-sheet flexibility using PCA on a database of beta-sheet structures. For sheets of varying dimension and geometry, we find two dominant modes of flexibility: twist and bend. The distributions of amplitudes for these modes are found to be Gaussian and independent, suggesting that the PCA twist and bend modes can be identified as the soft elastic normal modes of sheets. We consider the scaling of mode eigenvalues with sheet size and find that parallel beta-sheets are more rigid than antiparallel sheets over the entire range studied. Finally, we discuss the application of our PCA results to modeling and design of beta-sheet proteins.