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.     

Helix-helix interfaces and ligand binding

Helix-helix parallel interfaces can be characterized by certain combinations of amino acids, which repeatedly occur at core positions a and d (leucine zipper nomenclature) in homologous and nonhomologous proteins and influence interhelical angles. Applied for the prediction of interhelical angles in glutathione S-transferase, intracellular chloride channel and annexin molecules from various sources, correct results were achieved in 58 out of 62 proteins. Interhelical angles are found to correlate with the conformation of the glutathione S-transferase ligands glutathione, s-hexylglutathione, glutathione sulfonic acid, and glutathione-s-dinitrobenzene.     

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.     

A parallel coiled-coil tetramer with offset helices

Specific helix-helix interactions are fundamental in assembling the native state of proteins and in protein-protein interfaces. Coiled coils afford a unique model system for elucidating principles of molecular recognition between alpha helices. The coiled-coil fold is specified by a characteristic seven amino acid repeat containing hydrophobic residues at the first (a) and fourth (d) positions. Nonpolar side chains spaced three and four residues apart are referred to as the 3-4 hydrophobic repeat. The presence of apolar amino acids at the e or g positions (corresponding to a 3-3-1 hydrophobic repeat) can provide new possibilities for close-packing of alpha-helices that includes examples such as the lac repressor tetramerization domain. Here we demonstrate that an unprecedented coiled-coil interface results from replacement of three charged residues at the e positions in the dimeric GCN4 leucine zipper by nonpolar valine side chains. Equilibrium circular dichroism and analytical ultracentrifugation studies indicate that the valine-containing mutant forms a discrete alpha-helical tetramer with a significantly higher stability than the parent leucine-zipper molecule. The 1.35 A resolution crystal structure of the tetramer reveals a parallel four-stranded coiled coil with a three-residue interhelical offset. The local packing geometry of the three hydrophobic positions in the tetramer conformation is completely different from that seen in classical tetrameric structures yet bears resemblance to that in three-stranded coiled coils. These studies demonstrate that distinct van der Waals interactions beyond the a and d side chains can generate a diverse set of helix-helix interfaces and three-dimensional supercoil structures.     

Ionic Interactions Promote Transmembrane Helix-Helix Association Depending on Sequence Context

Folding and oligomerization of integral membrane proteins frequently depend on specific interactions of transmembrane helices. Interacting amino acids of helix-helix interfaces may form complex motifs and exert different types of molecular forces. Here, a set of strongly self-interacting transmembrane domains (TMDs), as isolated from a combinatorial library, was found to contain basic and acidic residues, in combination with polar nonionizable amino acids and C-terminal GxxxG motifs. Mutational analyses of selected sequences and reconstruction of high-affinity interfaces confirmed the cooperation of these residues in homotypic interactions. Probing heterotypic interaction indicated the presence of interhelical charge-charge interactions. Furthermore, simple motifs of an ionizable residue and GxxxG are significantly overrepresented in natural TMDs, and a specific combination of these motifs exhibits high-affinity heterotypic interaction. We conclude that intramembrane charge-charge interactions depend on sequence context. Moreover, they appear important for homotypic and heterotypic interactions of numerous natural TMDs.     

Driving Forces for Transmembrane α-Helix Oligomerization

We present what we believe to be a novel statistical contact potential based on solved structures of transmembrane (TM) α-helical bundles, and we use this contact potential to investigate the amino acid likelihood of stabilizing helix-helix interfaces. To increase statistical significance, we have reduced the full contact energy matrix to a four-flavor alphabet of amino acids, automatically determined by our methodology, in which we find that polarity is a more dominant factor of group identity than is size, with charged or polar groups most often occupying the same face, whereas polar/apolar residue pairs tend to occupy opposite faces. We found that the most polar residues strongly influence interhelical contact formation, although they occur rarely in TM helical bundles. Two-body contact energies in the reduced letter code are capable of determining native structure from a large decoy set for a majority of test TM proteins, at the same time illustrating that certain higher-order sequence correlations are necessary for more accurate structure predictions.     

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.     

Analysis and prediction of helix-helix interactions in membrane channels and transporters

Membrane proteins span a large variety of different functions such as cell-surface receptors, redox proteins, ion channels, and transporters. Proteins with functional pores show different characteristics of helix-helix packing as other helical membrane proteins. We found that the helix-helix contacts of 13 nonhomologous high-resolution structures of membrane channels and transporters are mainly accomplished by weakly polar amino acids (G > S > T > F) that preferably create contacts every fourth residue, typical for right-handed helix crossings. There is a strong correlation between the now available biological hydrophobicity scale and the propensities of the weakly polar and hydrophobic residues to be buried at helix-helix interfaces or to be exposed to the lipids in membrane channels and transporters. The polar residues, however, make no major contribution towards the packing of their transmembrane helices, and are therefore subsumed to be primarily exposed to the polar milieu during the folding process. The contact formation of membrane channels and transporters is therefore ruled by the solubility of the residues, which we suppose to be the driving force for the assembly of their transmembrane helices. By contrast, in 14 nonhomologous high-resolution structures of other membrane protein coils, also large and polar amino acids (D > S > M > Q) create characteristic contacts every 3.5th residues, which is a signature for left-handed helix crossings. Accordingly, it seems that dependent on the function, different concepts of folding and stabilization are realized for helical membrane proteins. Using a sequence-based matrix prediction method these differences are exploited to improve the prediction of buried and exposed residues of transmembrane helices significantly. When the sequence motifs typical for membrane channels and transporters were applied for the prediction of helix-helix contacts the quality of prediction rises by 16% to an average value of 76%, compared to the same approach when only single amino acid positions are taken into account.     

Complex Patterns of Histidine, Hydroxylated Amino Acids and the GxxxG Motif Mediate High-affinity Transmembrane Domain Interactions

Specific interactions of transmembrane helices play a pivotal role in the folding and oligomerization of integral membrane proteins. The helix-helix interfaces frequently depend on specific amino acid patterns. In this study, a heptad repeat pattern was randomized with all naturally occurring amino acids to uncover novel sequence motifs promoting transmembrane domain interactions. Self-interacting transmembrane domains were selected from the resulting combinatorial library by means of the ToxR/POSSYCCAT system. A comparison of the amino acid composition of high-and low-affinity sequences revealed that high-affinity transmembrane domains exhibit position-specific enrichment of histidine. Further, sequences containing His preferentially display Gly, Ser, and/or Thr residues at flanking positions and frequently contain a C-terminal GxxxG motif. Mutational analysis of selected sequences confirmed the importance of these residues in homotypic interaction. Probing heterotypic interaction indicated that His interacts in trans with hydroxylated residues. Reconstruction of minimal interaction motifs within the context of an oligo-Leu sequence confirmed that His is part of a hydrogen bonded cluster that is brought into register by the GxxxG motif. Notably, a similar motif contributes to self-interaction of the BNIP3 transmembrane domain.     

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.     

Higher-order interhelical spatial interactions in membrane proteins

Higher-order interactions are important for protein folding and assembly. We introduce the concept of interhelical three-body interactions as derived from Delaunay triangulation and alpha shapes of protein structures. In addition to glycophorin A, where triplets are strongly correlated with protein stability, we found that tight interhelical triplet interactions exist extensively in other membrane proteins, where many types of triplets occur far more frequently than in soluble proteins. We developed a probabilistic model for estimating the value of membrane helical interaction triplet (MHIT) propensity. Because the number of known structures of membrane proteins is limited, we developed a bootstrap method for determining the 95% confidence intervals of estimated MHIT values. We identified triplets that have high propensity for interhelical interactions and are unique to membrane proteins, e.g. AGF, AGG, GLL, GFF and others. A significant fraction (32%) of triplet types contains triplets that may be involved in interhelical hydrogen bond interactions, suggesting the prevalent and important roles of H-bond in the assembly of TM helices. There are several well-defined spatial conformations for triplet interactions on helices with similar parallel or antiparallel orientations and with similar right-handed or left-handed crossing angles. Often, they contain small residues and correspond to the regions of the closest contact between helices. Sequence motifs such as GG4 and AG4 can be part of the three-body interactions that have similar conformations, which in turn can be part of a higher-order cooperative four residue spatial motif observed in helical pairs from different proteins. In many cases, spatial motifs such as serine zipper and polar clamp are part of triplet interactions. On the basis of the analysis of the archaeal rhodopsin family of proteins, tightly packed triplet interactions can be achieved with several different choices of amino acid residues.     

Importance of potential interhelical salt-bridges involving interior residues for coiled-coil stability and quaternary structure

Coiled coils are formed by two or more alpha-helices that align in a parallel or an antiparallel relative orientation. Polar interactions involving residues at the interior a and d positions are important for determining the quaternary structure of coiled coils. In the model heterodimeric coiled-coil Acid-a1-Base-a1, a buried a-d' Asn-Asn interaction is sufficient to specify both a dimeric structure and an antiparallel relative helix orientation. Although the equivalent a-a' interaction is found in parallel coiled coils, there is no example of an a-d' Asn-Asn interaction in structurally characterized, naturally occurring antiparallel coiled coils. Instead, interior charged residues form interhelical salt-bridges with residues at the adjacent e or g positions. Using a model coiled-coil heterodimer, we have explored the role of a potential interhelical interaction between an Arg at an interior d position and a Glu at the adjacent g' position. Our results demonstrate that this potentially attractive interhelical Coulombic interaction has little or no influence on helix orientation. Instead, we show that burying a single Arg residue at an interior position is sufficient to specify a dimeric state at a significantly lower thermodynamic cost than burial of two interacting Asn residues.     

The structure of interfaces between subunits of dimeric and tetrameric proteins

The structures of the interfaces of nine dimeric and nine tetrameric proteins have been analyzed and have been seen to follow general principles. These interfaces are combinations of four structural motifs, which resemble features of monomeric proteins. These are: (i) extended beta sheet; (ii) helix-helix packing; (iii) sheet-sheet packing; and (iv) loop interactions. Other common structural features in the interfaces studied are two-fold symmetry, charged hydrogen bonds and channel formation (found only in tetramers). Monomer-monomer interfaces are intermediate in hydrophobicity and charge between the interfaces between secondary structures of monomeric proteins and the exteriors of monomeric proteins. A typical interface has one of the first three of the structural motifs at its centre and loop interactions around the outside, where most of the charge resides.     

A thermodynamic scale for leucine zipper stability and dimerization specificity: e and g interhelical interactions.

The leucine zipper is a dimeric coiled-coil protein structure composed of two amphipathic alpha-helices with the hydrophobic surfaces interacting to create the dimer interface. This structure has been found to mediate the dimerization of two abundant classes of DNA binding proteins: the bZIP and bHLH-Zip proteins. Several workers have reported that amino acids in the e and g positions of the coiled coil can modulate dimerization stability and specificity. Using the bZIP protein VBP as a host molecule, we report a thermodynamic scale (delta delta G) for 27 interhelical interactions in 35 proteins between amino acids in the g and the following e positions (g<==>e') of a leucine zipper coiled coil. We have examined the four commonly occurring amino acids in the e and g positions of bZIP proteins, lysine (K), arginine (R), glutamine (Q), glutamic acid (E), as well as the only other remaining charged amino acid aspartic acid (D), and finally alanine (A) as a reference amino acid. These results indicate that E<==>R is the most stable interhelical pair, being 0.35 kcal/mol more stable than E<==>K. A thermodynamic cycle analysis shows that the E<==>R pair is 1.33 kcal/mol more stable than A<==>A with -1.14 kcal/mol of coupling energy (delta delta Gint) coming from the interaction of E with R. The E<==>K coupling energy is only -0.14 kcal/mol. E interacts with more specificity than Q. The R<==>R pair is less stable than the K<==>K by 0.24 kcal/mol. R interacts with more specificity than K. Q forms more stable pairs with the basic amino acids K and R rather than with E. Changing amino acids in the e position to A creates bZIP proteins that form tetramers.     

Mutagenic mapping of helical structures in the transmembrane segments of the yeast alpha-factor receptor

The alpha-mating pheromone receptor encoded by the yeast STE2 gene is a G protein coupled receptor that initiates signaling via a MAP kinase pathway that prepares haploid cells for mating. To establish the range of allowed amino acid substitutions within transmembrane segments of this receptor, we conducted extensive random mutagenesis of receptors followed by screening for receptor function. A total of 157 amino acid positions in seven different mutagenic libraries corresponding to the seven predicted transmembrane segments were analyzed, yielding 390 alleles that retain at least 60 % of normal signaling function. These alleles contained a total of 576 unique amino acid substitutions, including 61 % of all the possible amino acid changes that can arise from single base substitutions. The receptor exhibits a surprising tolerance for amino acid substitutions. Every amino acid in the mutagenized regions of the transmembrane regions could be substituted by at least one other residue. Polar amino acids were tolerated in functional receptors at 115 different positions (73 % of the total). Hydrophobic amino acids were tolerated in functional receptors at all mutagenized positions. Substitutions introducing proline residues were recovered at 53 % of all positions where they could be brought about by single base changes. Residues with charged side-chains could also be tolerated at 53 % of all positions where they were accessible through single base changes. The spectrum of allowed amino acid substitutions was characterized in terms of the hydrophobicity, radius of gyration, and charge of the allowed substitutions and mapped onto alpha-helical structures. By comparing the patterns of allowed substitutions with the recently determined structure of rhodopsin, structural features indicative of helix-helix interactions can be discerned in spite of the extreme sequence divergence between these two proteins.     

Interhelical angle and distance preferences in globular proteins

Orientational preferences between interacting helices within globular proteins have been studied extensively over the years. A number of classical structural models such as "knobs into holes" and "ridges into grooves" were developed decades ago to explain perceived preferences in interhelical angle distributions. In contrast, relatively recent works have examined statistical biases in angular distributions which result from spherical geometric effects. Those works have concluded that the predictions of classical models are due in large part to these biases. In this article we perform an analysis on the largest set of helix-helix interactions within high-resolution structures of nonhomologous proteins studied to date. We examine the interhelical angle distribution as a function of spatial distance between helix pairs. We show that previous efforts to normalize angle distribution data did not include two important effects: 1), helices can interact with each other in three distinct ways which we refer to as "line-on-line," "endpoint-to-line," and "endpoint-to-endpoint," and each of these interactions has its own geometric effects which must be included in the proper normalization of data; and 2), all normalizations that depend on geometric parameters such as interhelical angle must occur before the data is binned to avoid artifacts of bin size from biasing the conclusions. Taking these two points into account, we find that there are very pronounced preferences for helices to interact at angles of approximately +/-160 and +/-20 degrees in the line-on-line case. This pattern persists when the closest alpha-carbons in the helices vary from 4 to 12 A. The endpoint-to-line and endpoint-to-endpoint cases also exhibit distinct preferences when the data is normalized properly. Analysis of the local structural interactions which give rise to these preferences has not been studied here and is left for future work.     

Substitution rates in α-helical transmembrane proteins

It has been shown previously that some membrane proteins have a conserved core of amino acid residues. This idea not only serves to orient helices during model building exercises but may also provide insight into the structural role of residues mediating helix-helix interactions. Using experimentally determined high-resolution structures of alpha-helical transmembrane proteins we show that, of the residues within the hydrophobic transmembrane spans, the residues at lipid and subunit interfaces are more evolutionarily variable than those within the lipid-inaccessible core of a polypeptide's transmembrane domain. This supports the idea that helix-helix interactions within the same polypeptide chain and those at the interface between different polypeptide chains may arise in distinct ways. To show this, we use a new method to estimate the substitution rate of an amino acid residue given an alignment and phylogenetic tree of closely related proteins. This method gives better sensitivity in the otherwise-conserved transmembrane domains than a conventional similarity analysis and is relatively insensitive to the sequences used.     

Side-chain conformations in 4-alpha-helical bundles.

The distribution of the chi(1), chi(2) dihedral angles in a dataset consisting of 12 unrelated 4-alpha-helical bundle proteins was determined and qualitatively compared with that observed in globular proteins. The analysis suggests that the 4-alpha-helical bundle motif could occasionally impose steric constraints on side chains: (i) the side-chain conformations are limited to only a subset of the conformations observed in globular proteins and for some amino acids they are sterically more constrained than those in helical regions of globular proteins; (ii) aspartic acid and asparagine occasionally adopt rotamers that have not been previously reported for globular or helical proteins; (iii) some rotamers of tyrosine and isoleucine are predominantly or exclusively associated with hydrophobic core positions (a, d); (iv) mutations in the hydrophobic core occur preferentially between residue types which among other physicochemical properties also share a predominant rotamer.     

Molecular dynamics simulations of the E1/E2 transmembrane domain of the Semliki Forest virus

Transmembrane (TM) helix-helix interactions are important for virus budding and fusion. We have developed a simulation strategy that reveals the main features of the helical packing between the TM domains of the two glycoproteins E1 and E2 of the alpha-virus Semliki Forest virus and that can be extrapolated to sketch TM helical packing in other alpha-viruses. Molecular dynamics simulations were performed in wild-type and mutant peptides, both isolated and forming E1/E2 complexes. The simulations revealed that the isolated wild-type E1 peptide formed a more flexible helix than the rest of peptides and that the wild-type E1/E2 complex consists of two helices that intimately pack their N-terminals. The residues located at the interhelical interface displayed the typical motif of the left-handed coiled-coils. These were small and medium residues as Gly, Ala, Ser, and Leu, which also had the possibility to form interhelical Calpha-H...O hydrogen bonds. Results from the mutant complexes suggested that correct packing is a compromise between these residues at both E1 and E2 interhelical interfaces. This compromise allowed prediction of E1-E2 contact residues in the TM spanning domain of other alphaviruses even though the sequence identity of E2 peptides is low in this domain.