Interaction of transmembrane helices by a knobs-into-holes packing characteristic of soluble coiled coils

Membrane-embedded protein domains frequently exist as α-helical bundles, as exemplified by photosynthetic reaction centers, bacteriorhodopsin, and cytochrome C oxidase. The sidechain packing between their transmembrane helices was investigated by a nearest-neighbor analysis which identified sets of interfacial residues for each analyzed helix–helix interface. For the left-handed helix–helix pairs, the interfacial residues almost exclusively occupy positions a, d, e, or g within a heptad motif (abcdefg) which is repeated two to three times for each interacting helical surface. The connectivity between the interfacial residues of adjacent helices conforms to the knobs-into-holes type of sidechain packing known from soluble coiled coils. These results demonstrate on a quantitative basis that the geometry of sidechain packing is similar for left-handed helix–helix pairs embedded in membranes and coiled coils of soluble proteins. The transmembrane helix–helix interfaces studied are somewhat less compact and regular as compared to soluble coiled coils and tolerate all hydrophobic amino acid types to similar degrees. The results are discussed with respect to previous experimental findings which demonstrate that specific interactions between transmembrane helices are important for membrane protein folding and/or oligomerization. Proteins 31:150–159, 1998. © 1998 Wiley-Liss, Inc.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

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.     

Robust Driving Forces for Transmembrane Helix Packing

The packing structures of transmembrane helices are traditionally attributed to patterns in residues along the contact surface. In this view, besides keeping the helices confined in the membrane, the bilayer has only a minor effect on the helices structure. Here, we use two different approaches to show that the lipid environment has a crucial effect in determining the cross-angle distribution of packed helices. We analyzed structural data of a membrane proteins database. We show that the distribution of cross angles of helix pairs in this database is statistically indistinguishable from the cross-angle distribution of two noninteracting helices imbedded in the membrane. These results suggest that the cross angle is, to a large extent, determined by the tilt angle of the individual helices. We test this hypothesis using molecular simulations of a coarse-grained model that contains no specific residue interactions. These simulations reproduce the same cross-angle distribution as found in the database. As the tilt angle of a helix is dominated by hydrophobic mismatch between the protein and surrounding lipids, our results indicate that hydrophobic mismatch is the dominant factor guiding the transmembrane helix packing. Other short-range forces might then fine-tune the structure to its final configuration.     

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.     

Current status of membrane protein structure classification

For over 2 decades, continuous efforts to organize the jungle of available protein structures have been underway. Although a number of discrepancies between different classification approaches for soluble proteins have been reported, the classification of membrane proteins has so far not been comparatively studied because of the limited amount of available structural data. Here, we present an analysis of α‐helical membrane protein classification in the SCOP and CATH databases. In the current set of 63 α‐helical membrane protein chains having between 1 and 13 transmembrane helices, we observed a number of differently classified proteins both regarding their domain and fold assignment. The majority of all discrepancies affect single transmembrane helix, two helix hairpin, and four helix bundle domains, while domains with more than five helices are mostly classified consistently between SCOP and CATH. It thus appears that the structural constraints imposed by the lipid bilayer complicate the classification of membrane proteins with only few membrane‐spanning regions. This problem seems to be specific for membrane proteins as soluble four helix bundles, not restrained by the membrane, are more consistently classified by SCOP and CATH. Our findings indicate that the structural space of small membrane helix bundles is highly continuous such that even minor differences in individual classification procedures may lead to a significantly different classification. Membrane proteins with few helices and limited structural diversity only seem to be reasonably classifiable if the definition of a fold is adapted to include more fine‐grained structural features such as helix–helix interactions and reentrant regions. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Hydrogen-bonding and packing features of membrane proteins: functional implications

The recent structural elucidation of about one dozen channels (in which we include transporters) has provided further evidence that these membrane proteins typically undergo large movements during their function. However, it is still not well understood how these proteins achieve the necessary trade-off between stability and mobility. To identify specific structural properties of channels, we compared the helix-packing and hydrogen-bonding patterns of channels with those of membrane coils; the latter is a class of membrane proteins whose structures are expected to be more rigid. We describe in detail how in channels, helix pairs are usually arranged in packing motifs with large crossing angles (|τ| ≈ 40°), where the (small) side chains point away from the packing core and the backbones of the two helices are in close contact. We found that this contributes to a significant enrichment of Cα-H…O bonds and to a packing geometry where right-handed parallel (τ = −40° ± 10°) and antiparallel (τ = +140° ± 25°) arrangements are equally preferred. By sharp contrast, the interdigitation and hydrogen bonding of side chains in helix pairs of membrane coils results in narrowly distributed left-handed antiparallel arrangements with crossing angles τ = −160° ± 10° (|τ| ≈ 20°). In addition, we show that these different helix-packing modes of the two types of membrane proteins correspond to specific hydrogen-bonding patterns. In particular, in channels, three times as many of the hydrogen-bonded helix pairs are found in parallel right-handed motifs than are non-hydrogen-bonded helix pairs. Finally, we discuss how the presence of weak hydrogen bonds, water-containing cavities, and right-handed crossing angles may facilitate the required conformational flexibility between helix pairs of channels while maintaining sufficient structural stability.     

Interhelical hydrogen bonding drives strong interactions in membrane proteins

Polar residues in transmembrane alpha-helices may strongly influence the folding or association of integral membrane proteins. To test whether a motif that promotes helix association in a soluble protein could do the same within a membrane, we designed a model transmembrane helix based on the GCN4 leucine zipper. We found in both detergent micelles and biological membranes that helix association is driven strongly by asparagine, independent of the rest of the hydrophobic leucine and/or valine sequence. Hydrogen bonding between membrane helices gives stronger associations than the packing of surfaces in glycophorin A helices, creating an opportunity to stabilize structures, but also implying a danger that non-specific interactions might occur. Thus, membrane proteins may fold to avoid exposure of strongly hydrogen bonding groups at their lipid exposed surfaces.     

Structural role of glycine in amyloid fibrils formed from transmembrane alpha-helices

Amyloid fibrils associated with diseases such as Alzheimer's are often derived from the transmembrane helices of membrane proteins. It is known that the fibrils have a cross-beta-sheet structure where main chain hydrogen bonding occurs between beta-strands in the direction of the fibril axis. However, the structural basis for how the membrane-spanning helix is converted into a beta-sheet or how protofibrils associate into fibrils is not known. Here, we use a model peptide corresponding to a portion of the single transmembrane helix of glycophorin A to investigate the structural role of glycine in amyloid-like fibrils formed from transmembrane helices. Glycophorin A contains a GxxxG motif that is found in many transmembrane sequences including that of the amyloid precursor protein and prion protein. We propose that glycine, which mediates helix interactions in membrane proteins, also provides key packing motifs when it occurs in beta-sheets. We show that glycines in the glycophorin A transmembrane helix promote extended beta-strand formation when the helix partitions into aqueous environments and stabilize the packing of beta-sheets in the formation of amyloid-like fibrils. We demonstrate that fibrillization can be disrupted with a new class of inhibitors that target the molecular grooves created by glycine.     

Fluorescence resonance energy transfer shows a close helix-helix distance in the transmembrane M13 procoat protein

During the membrane insertion process the major coat protein of bacteriophage M13 assumes a conformation in which two transmembrane helices corresponding to the leader sequence and the anchor region in the mature part of the protein coming into close contact with each other. Previous studies on the molecular mechanism of membrane insertion of M13 procoat protein have shown that this interaction between the two helices might drive the actual translocation process. We investigated the intramolecular distance between the two helices of the transmembrane procoat protein by measuring fluorescence resonance energy transfer (FRET) between the donor (Tyr) placed in one helix and the acceptor (Trp) placed in the other helix. Various mutant procoat proteins with differently positioned donor-acceptor pairs were generated, purified, and reconstituted into artificial lipid bilayers. The results obtained from the FRET measurements, combined with molecular modeling, show that the transmembrane helices are in close contact on the order of 1-1.5 nm. The present approach might be of general interest for determining the topology and the folding of membrane proteins.     

Close packing of helices 3 and 12 of 16 S rRNA is required for the normal ribosome function

The along-groove packing motif is a quasi-reciprocal arrangement of two RNA double helices in which a backbone of each helix is closely packed within the minor groove of the other helix. At the center of the inter-helix contact, a GU base pair in one helix packs against a Watson-Crick base pair in the other helix. Here, based on in vivo selection from a combinatorial gene library of 16 S rRNA and on functional characterization of the selected clones, we demonstrate that the normal ribosome performance requires that helices 3 and 12 be closely packed. In some clones the Watson-Crick and GU base pairs exchange in their positions between the two helices, which affects neither the quality of the helix packing, nor the ribosome function. On the other hand, perturbations in the close packing usually lead to a substantial drop in the ribosome activity. The functionality of the clones containing such perturbations may depend on the presence of particular elements in the vicinity of the area of contact between helices 3 and 12. Such cases do not exist in natural 16 S rRNA, and their selection enriches our knowledge of the constraints imposed on the structure of ribosomal RNA in functional ribosomes.     

High cationic charge and bilayer interface-binding helices in a regulatory lipid glycosyltransferase

In the prokaryote Acholeplasma laidlawii, membrane bilayer properties are sensed and regulated by two interface glycosyltransferases (GTs), synthesizing major nonbilayer- (alMGS GT) and bilayer-prone glucolipids. These enzymes are of similar structure, as many soluble GTs, but are sensitive to lipid charge and curvature stress properties. Multivariate and bioinformatic sequence analyses show that such interface enzymes, in relation to soluble ones of similar fold, are characterized by high cationic charge, certain distances between small and cationic amino acids, and by amphipathic helices. Varying surface contents of Lys/Arg pairs and Trp indicate different membrane-binding subclasses. A predicted potential (cationic) binding helix from alMGS was structurally verified by solution NMR and CD. The helix conformation was induced by a zwitterionic as well as anionic lipid environment, and the peptide was confined to the bilayer interface. Bilayer affinity of the peptide, analyzed by surface plasmon resonance, was higher than that for soluble membrane-seeking proteins/peptides and rose with anionic lipid content. Interface intercalation was supported by phase equilibria in membrane lipid mixtures, analyzed by 31P NMR and DSC. An analogous, potentially binding helix has a similar location in the structurally determined Escherichia coli cell wall precursor GT MurG. These two helices have little sequence conservation in alMGS and MurG homologues but maintain their amphipathic character. The evolutionary modification of the alMGS binding helix and its location close to the acceptor substrate site imply a functional importance in enzyme catalysis, potentially providing a mechanism by which glycolipid synthesis will be sensitive to membrane surface charge and intrinsic curvature strain.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Structural comparison and classification of alpha‐helical transmembrane domains based on helix interaction patterns

Structural classification of membrane proteins is still in its infancy due to the relative paucity of available three‐dimensional structures compared with soluble proteins. However, recent technological advances in protein structure determination have led to a significant increase in experimentally known membrane protein folds, warranting exploration of the structural universe of membrane proteins. Here, a new and completely membrane protein specific structural classification system is introduced that classifies α‐helical membrane proteins according to common helix architectures. Each membrane protein is represented by a helix interaction graph depicting transmembrane helices with their pairwise interactions resulting from individual residue contacts. Subsequently, proteins are clustered according to similarities among these helix interaction graphs using a newly developed structural similarity score called HISS. As HISS scores explicitly disregard structural properties of loop regions, they are more suitable to capture conserved transmembrane helix bundle architectures than other structural similarity scores. Importantly, we are able to show that a classification approach based on helix interaction similarity closely resembles conventional structural classification databases such as SCOP and CATH implying that helix interactions are one of the major determinants of α‐helical membrane protein folds. Furthermore, the classification of all currently available membrane protein structures into 20 recurrent helix architectures and 15 singleton proteins demonstrates not only an impressive variability of membrane helix bundles but also the conservation of common helix interaction patterns among proteins with distinctly different sequences. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Helix kinks are equally prevalent in soluble and membrane proteins

Helix kinks are a common feature of α‐helical membrane proteins, but are thought to be rare in soluble proteins. In this study we find that kinks are a feature of long α‐helices in both soluble and membrane proteins, rather than just transmembrane α‐helices. The apparent rarity of kinks in soluble proteins is due to the relative infrequency of long helices (≥20 residues) in these proteins. We compare length‐matched sets of soluble and membrane helices, and find that the frequency of kinks, the role of Proline, the patterns of other amino acid around kinks (allowing for the expected differences in amino acid distributions between the two types of protein), and the effects of hydrogen bonds are the same for the two types of helices. In both types of protein, helices that contain Proline in the second and subsequent turns are very frequently kinked. However, there are a sizeable proportion of kinked helices that do not contain a Proline in either their sequence or sequence homolog. Moreover, we observe that in soluble proteins, kinked helices have a structural preference in that they typically point into the solvent. Proteins 2014; 82:1960–1970. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Proximity relationships between helices I and XI or XII in the lactose permease of Escherichia coli determined by site-directed thiol cross-linking.

The lactose permease of Escherichia coli was expressed in two fragments (split permease), each with a Cys residue, and cross-linking was studied. Split permease with a discontinuity in either loop II/III (N2C10permease) or loop VI/VII (N6C6permease) was used. Proximity of multiple pairs of Cys residues in helices I and XI or XII was examined by using three homobifunctional thiol-specific cross-linking reagents of different lengths and flexibilities (6 A, rigid; 10 A, rigid; 16 A, flexible) or iodine. Cys residues in the periplasmic half of helix I cross-link to Cys residues in the periplasmic half of helix XI. In contrast, no cross-linking is evident with paired Cys residues near the cytoplasmic ends of helices I and XI. Therefore, the periplasmic halves of helices I and XI are in close proximity, and the helices tilt away from each other towards the cytoplasmic face of the membrane. Cross-linking is also found with paired Cys residues near the middle of helices I and XII, but not with paired Cys residues near either end of the helices. Thus, helices I and XII are in close proximity only in the approximate middle of the membrane. Based on the findings, a modified helix packing model is proposed.     

General architecture of the alpha-helical globule

A model is presented for the arrangement of alpha-helices in globular proteins. In the model, helices are placed on certain ribs of "quasi-spherical" polyhedra. The polyhedra are chosen so as to allow the close packing of helices around a hydrophobic core and to stress the collective interactions of the individual helices. The model predicts a small set of stable architectures for alpha-helices in globular proteins and describes the geometries of the helix packings. Some of the predicted helix arrangements have already been observed in known protein structures; others are new. An analysis of the three-dimensional structures of all proteins for which co-ordinates are available shows that the model closely approximates the arrangements and packing of helices actually observed. The average deviations of the real helix axes from those in the model polyhedra is +/- 20 degrees in orientation and +/- 2 A in position (1 A = 0.1 nm). We also show that for proteins that are not homologous, but whose helix arrangements are described by the same polyhedron, the root-mean-square difference in the position of the C alpha atoms in the helices is 1.6 to 3.0 A.     

Global secondary structure packing angle bias in proteins

While studies of secondary structure interactions have focused on local interacting features, there is a need for a more global characterization of packing-induced aligned packing of secondary structures. This study presents an analysis of the distribution of globally sampled secondary structures within selected subunits of a selected set of multimeric proteins. Comparisons are made between the distribution of the cosines of angles between triplets of linear segments associated to secondary structures and a theoretically obtained distribution for triplets of random uniformly distributed unit vectors. We show that, among all pairs of helix or strand segments, planar configurations appear more frequently than expected for uniformly distributed vectors, and alignment is strongly preferred compared to that expected for uniformly distributed vector triplets. Among all secondary structure triplets, pairs of angle cosines between helix strand segments deviate from uniformity corresponding to alignment and anti-alignment. Furthermore, among all helix or strand segments, including non-interacting secondary structures, the distribution of a single angle cosine indicates a strong preference for alignment and anti-alignment. Selection for interactive triplets shows results consistent with prior studies. Lastly, angle pairs are not statistically independent, indicating that alignment between two helix or strand segments is more likely if another helix or strand is aligned with either of the first two helices or strands. Selection for interactive segment triplets shows results consistent with prior studies.