Identifying a Folding Nucleus for the Lysozyme/α-Lactalbumin Family from Sequence Conservation Clusters

Four subfamilies of c-type lysozyme and one subfamily of α-lactalbumin are defined from 78 sequences, and their folding nucleus is identified with a method based on conserved residues and native structural contacts between pairs of conserved residues. One large cluster of 19 conserved residues is found which is mostly nonpolar, buried, and nonfunctional. It can be subdivided into three subclusters: (1) conserved residues in four helices; (2) conserved residues that stabilize the connector between the α and the β domains; and (3) a β-turn, sitting in the middle of a bowl of α-helix residues. It is proposed that this folding nucleus initiates four helices, A, B, C, and D, three β sheets, and the connector, which corresponds closely to the nucleation of the so-called fast folding track pathway. As the secondary structures propagate, nonconserved residues and functionally conserved residues would form additional contacts. The conserved residues are selected with a phylogenetic scheme in which single members of subfamilies are selected. Subfamilies are then equally weighted to obtain the consensus conservation. Received: 11 June 2001 / Accepted: 28 August 2001     

In silico activation of KcsA K+ channel by lateral forces applied to the C-termini of inner helices

Crystallographic studies of K(+) channels in the closed (KcsA) and open (MthK) states suggest that Gly(99) (KcsA numbering) in the inner helices serves as a gating hinge during channel activation. However, some P-loop channels have larger residues in the corresponding position. The comparison of x-ray structures of KcsA and MthK shows that channel activation alters backbone torsions and helical H-bonds in residues 95-105. Importantly, the changes in Gly(99) are not the largest ones. This raises questions about the mechanism of conformational changes upon channel gating. In this work, we have built a model of the open KcsA using MthK as a template and simulated opening and closing of KcsA by constraining C-ends of the inner helices at a gradually changing distance from the pore axis without restraining mobility of the helices along the axis. At each imposed distance, the energy was Monte Carlo-minimized. The channel-opening and channel-closing trajectories arrived to the structures in which the backbone geometry was close to that seen in MthK and KcsA, respectively. In the channel-opening trajectory, the constraints-induced lateral forces caused kinks at midpoints of the inner helices between Val(97) and Gly(104) but did not destroy interdomain contacts, the pore helices, and the selectivity filter. The simulated activation of the Gly(99)Ala mutant yielded essentially similar results. Analysis of interresidue energies shows that the N-terminal parts of the inner helices form strong attractive contacts with the pore helices and the outer helices. The lateral forces induce kinks at the position where the helix-breaking torque is maximal and the intersegment contacts vanish. This mechanism may be conserved in different P-loop channels.     

Solid-phase synthesis of short α-helices stabilized by the hydrogen bond surrogate approach

Stabilized α-helices and nonpeptidic helix mimetics have emerged as powerful molecular scaffolds for the discovery of protein-protein interaction inhibitors. Protein-protein interactions often involve large contact areas, which are often difficult for small molecules to target with high specificity. The hypothesis behind the design of stabilized helices and helix mimetics is that these medium-sized molecules may pursue their targets with higher specificity because of a larger number of contacts. This protocol describes an optimized synthetic strategy for the preparation of stabilized α-helices that feature a carbon-carbon linkage in place of the characteristic N-terminal main-chain hydrogen bond of canonical helices. Formation of the carbon-carbon bond is enabled by a microwave-assisted ring-closing metathesis reaction between two terminal olefins on the peptide chain. The outlined strategy allows the synthesis and purification of a hydrogen bond surrogate (HBS) α-helix in ～ 1 week.     

Non-equivalent role of TM2 gating hinges in heteromeric Kir4.1/Kir5.1 potassium channels

Comparison of the crystal structures of the KcsA and MthK potassium channels suggests that the process of opening a K⁺ channel involves pivoted bending of the inner pore-lining helices at a highly conserved glycine residue. This bending motion is proposed to splay the transmembrane domains outwards to widen the gate at the “helix-bundle crossing”. However, in the inwardly rectifying (Kir) potassium channel family, the role of this “hinge” residue in the second transmembrane domain (TM2) and that of another putative glycine gating hinge at the base of TM2 remain controversial. We investigated the role of these two positions in heteromeric Kir4.1/Kir5.1 channels, which are unique amongst Kir channels in that both subunits lack a conserved glycine at the upper hinge position. Contrary to the effect seen in other channels, increasing the potential flexibility of TM2 by glycine substitutions at the upper hinge position decreases channel opening. Furthermore, the contribution of the Kir4.1 subunit to this process is dominant compared to Kir5.1, demonstrating a non-equivalent contribution of these two subunits to the gating process. A homology model of heteromeric Kir4.1/Kir5.1 shows that these upper “hinge” residues are in close contact with the base of the pore α-helix that supports the selectivity filter. Our results also indicate that the highly conserved glycine at the “lower” gating hinge position is required for tight packing of the TM2 helices at the helix-bundle crossing, rather than acting as a hinge residue.     

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.     

Intramolecular steric effects and hydrogen bonding in regular conformations of polyamino acids

A survey has been made, by using computer methods, of the types of helices which polypeptide chains can form, taking into account steric requirements and intramolecular hydrogen-bonding interactions. The influence on these two requirements, of small variations in the bond angles of the peptide residues, or of small changes in the overall dimensions of the helix (pitch and residues per turn), have been assessed for the special case of the α-helix. Criteria for the formation of acceptable hydrogen bonds have also been applied to helices of other types, viz., the 3, γ^?, ω^?, and π-helices. It was shown that the N? H … O and H … O? C angles in hydrogen bonds are sensitive to changes in either the NC^αC′ bond angle or in the rotational angles about the N? C^α and C^α? C′ bonds. However, the variants of the α-helix observed experimentally in myoglobin can all be constructed without distortion of the hydrogen bonds. For α-helices, the steric and hydrogen bonding requirements are more easily fulfilled with an NC^αC′ bond angle of 111°, rather than 109.5°. The decreased stability observed for the left-handed α-helix relative to the right-handed one for L -amino acids is due essentially only to interactions of the C^β atom of the side chains with atoms in adjacent peptide units in the backbone, and interactions with atoms in adjacent turns of the helical backbone are not significantly different in the two helices. Restrictions in the freedom of rotation of bulky side chains may have significant kinetic effects during the formation of the α-helix from the “random coil” state.     

Dynamics and hydration of the alpha-helices of apolipophorin III

Apolipophorin III (apoLp-III) is an exchangeable apolipoprotein whose structure is represented as a bundle of five amphipathic alpha-helices. In order to study the properties of the helical domains of apolipophorin III, we designed and obtained five single-tryptophan mutants of Locusta migratoria apoLp-III. The proteins were studied by UV absorption spectroscopy, time-resolved and steady-state fluorescence spectroscopy, and circular dichroism. Fluorescence anisotropy, near-UV CD and solute fluorescence quenching studies indicate that the Trp residues in helices 1 (N-terminal) and 5 (C-terminal) have the highest conformational flexibility. These two residues also showed the highest degree of hydration. Trp residues in helices 3 and 4 display the lowest mobility, as assessed by fluorescence anisotropy and near UV CD. The Trp residue in helix 2 is protected from the solvent but shows high mobility. As inferred from the properties of the Trp residues, helices 1 and 5 appear to have the highest conformational flexibility. Helix 2 has an intermediate mobility, whereas helices 3 and 4 appear to constitute a highly ordered domain. From the configuration of the helices in the tertiary structure of the protein, we estimated the relative strength of the five interhelical interactions of apoLp-III. These interactions can be ordered according to their apparent stabilizing strengths as: helix 3-helix 4 > helix 2-helix 3 > helix 4-helix 1 approximately helix 2-helix 5 > helix 1-helix 5. A new model for the conformational change that is expected to occur upon binding of the apolipoprotein to lipid is proposed. This model is significantly different from the currently accepted model (Breiter, D. R., Kanost, M. R., Benning, M. M., Wesemberg, G., Law, J. H., Wells, M. A., Rayment, I., and Holden, M. (1991) Biochemistry 30, 603-608). The model presented here predicts that the relaxation of the tertiary structure and the concomitant exposure of the hydrophobic core take place through the disruption of the weak interhelical contacts between helices 1 and 5. To some extent, the weakness of the helix 1-helix 5 interaction would be due to the parallel arrangement of these helices.     

Solution structure of P01, a natural scorpion peptide structurally analogous to scorpion toxins specific for apamin-sensitive potassium channel

The venom of the North African scorpion Androctonus mauretanicus mauretanicus possesses numerous highly active neurotoxins that specifically bind to various ion channels. One of these, P05, has been found to bind specifically to calcium-activated potassium channels and also to compete with apamin, a toxin extracted from bee venom. Besides the highly potent ones, several of these peptides (including that of P01) have been purified and been found to possess only a very weak, although significant, activity in competition with apamin. The amino acid sequence of P01 shows that it is shorter than P05 by two residues. This deletion occurs within an α-helix stretch (residues 5–12). This α-helix has been shown to be involved in the interaction of P05 with its receptor via two arginine residues. These two arginines are absent in the P01 sequence. Furthermore, a proline residue in position 7 of the P01 sequence may act as an α-helix breaker. We have determined the solution structure of P01 by conventional two-dimensional ¹H nuclear magnetic resonance and show that 1) the proline residue does not disturb the α-helix running from residues 5 to 12; 2) the two arginines are topologically replaced by two acidic residues, which explains the drop in activity; 3) the residual binding activity may be due to the histidine residue in position 9; and 4) the overall secondary structure is conserved, i.e., an α-helix running from residues 5 to 12, two antiparallel stretches of β-sheet (residues 15–20 and 23–27) connected by a type I′ β-turn, and three disulfide bridges connecting the α-helix to the β-sheet.     

A conserved π-helix plays a key role in thermoadaptation of catalysis in the glycoside hydrolase family 4

Here, we characterize the role of a π-helix in the molecular mechanisms underlying thermoadaptation in the glycoside hydrolase family 4 (GH4). The interspersed π-helix present in a subgroup is evolutionarily related to a conserved α-helix in other orthologs by a single residue insertion/deletion event. The insertional residue, Phe407, in a hyperthermophilic α-glucuronidase, makes specific interactions across the inter-subunit interface. In order to establish the sequence-structure-stability implications of the π-helix, the wild-type and the deletion variant (Δ407) were characterized. The variant showed a significant lowering of melting temperature and optimum temperature for the highest activity. Crystal structures of the proteins show a transformation of the π-helix to a continuous α-helix in the variant, identical to that in orthologs lacking this insertion. Thermodynamic parameters were determined from stability curves representing the temperature dependence of unfolding free energy. Though the proteins display maximum stabilities at similar temperatures, a higher melting temperature in the wild-type is achieved by a combination of higher enthalpy and lower heat capacity of unfolding. Comparisons of the structural changes, and the activity and thermodynamic profiles allow us to infer that specific non-covalent interactions, and the existence of residual structure in the unfolded state, are crucial determinants of its thermostability. These features permit the enzyme to balance the preservation of structure at a higher temperature with the thermodynamic stability required for optimum catalysis.     

Dynamic Motions of the HIV-1 Frameshift Site RNA

The HIV-1 frameshift site (FS) plays a critical role in viral replication. During translation, the HIV-1 FS transitions from a 3-helix to a 2-helix junction RNA secondary structure. The 2-helix junction structure contains a GGA bulge, and purine-rich bulges are common motifs in RNA secondary structure. Here, we investigate the dynamics of the HIV-1 FS 2-helix junction RNA. Interhelical motions were studied under different ionic conditions using NMR order tensor analysis of residual dipolar couplings. In 150 mM potassium, the RNA adopts a 43°(±4°) interhelical bend angle (β) and displays large amplitude, anisotropic interhelical motions characterized by a 0.52(±0.04) internal generalized degree of order (GDO_int) and distinct order tensor asymmetries for its two helices (η = 0.26(±0.04) and 0.5(±0.1)). These motions are effectively quenched by addition of 2 mM magnesium (GDO_int = 0.87(±0.06)), which promotes a near-coaxial conformation (β = 15°(±6°)) of the two helices. Base stacking in the bulge was investigated using the fluorescent purine analog 2-aminopurine. These results indicate that magnesium stabilizes extrahelical conformations of the bulge nucleotides, thereby promoting coaxial stacking of helices. These results are highly similar to previous studies of the HIV transactivation response RNA, despite a complete lack of sequence similarity between the two RNAs. Thus, the conformational space of these RNAs is largely determined by the topology of their interhelical junctions.     

High resolution structure of BipD: an invasion protein associated with the type III secretion system of Burkholderia pseudomallei

Burkoldheria pseudomallei is a Gram-negative bacterium that possesses a protein secretion system similar to those found in Salmonella and Shigella. Recent work has indicated that the protein encoded by the BipD gene of B. pseudomallei is an important secreted virulence factor. BipD is similar in sequence to IpaD from Shigella and SipD from Salmonella and is therefore likely to be a translocator protein in the type-III secretion system of B. pseudomallei. The crystal structure of BipD has been solved at a resolution of 2.1 A revealing the detailed tertiary fold of the molecule. The overall structure is appreciably extended and consists of a bundle of antiparallel alpha-helical segments with two small beta-sheet regions. The longest helices of the molecule form a four-helix bundle and most of the remaining secondary structure elements (three helices and two three-stranded beta-sheets) are formed by the region linking the last two helices of the four-helix bundle. The structure suggests that the biologically active form of the molecule may be a dimer formed by contacts involving the C-terminal alpha-helix, which is the most strongly conserved part of the protein. Comparison of the structure of BipD with immunological and other data for IpaD indicates that the C-terminal alpha-helix is also involved in contacts with other proteins that form the translocon.     

Effects of single-base bulges on intercalator binding to small RNA and DNA hairpins and a ribosomal RNA fragment

The way in which a single-base bulge might affect the structure of an RNA helix has been examined by preparing a series of six RNA hairpins, all with seven base pairs and a four-nucleotide loop. Five of the hairpins have single-base bulges at different positions. The intercalating cleavage reagent (methidiumpropyl)-EDTA-Fe(II) [MPE-Fe(II)] binds preferentially at a CpG sequence in the helix lacking a bulge and in four of the five hairpins with bulges. Hairpins with a bulge one or two bases to the 3' side of the CpG sequence bind ethidium 4-5-fold more strongly than the others. V1 RNase, which is sensitive to RNA backbone conformation in helices, detects a conformational change in all of the helices when ethidium binds; the most dramatic changes, involving the entire hairpin stem, are in one of the two hairpins with enhanced ethidium affinity. Only a slight conformational change is detected in the hairpin lacking a bulge. A bulge adjacent to a CpG sequence in a 100-nucleotide ribosomal RNA fragment enhances MPE-Fe(II) binding by an order of magnitude. These results extend our previous observations of bulges at a single position in an RNA hairpin [White, S. A., & Draper, D.E. (1987) Nucleic Acids Res. 15, 4049] and show that (1) a structural change in an RNA helix may be propagated for several base pairs, (2) bulges tend to increase the number of conformations available to a helix, and (3) the effects observed in small RNA hairpins are relevant to larger RNAs with more extensive structure. A bulge in a DNA hairpin identical in sequence with the RNA hairpins does not enhance MPE-Fe(II) binding affinity, relative to a control DNA hairpin. The effects of bulges on ethidium intercalation are evidently modulated by helix structure.     

Evolutionary origin of a secondary structure: π-helices as cryptic but widespread insertional variations of α-helices that enhance protein functionality

Formally annotated π-helices are rare in protein structures but have been correlated with functional sites. Here, we analyze protein structures to show that π-helices are the same as structures known as α-bulges, α-aneurisms, π-bulges, and looping outs, and are evolutionarily derived by the insertion of a single residue into an α-helix. This newly discovered evolutionary origin explains both why π-helices are cryptic, being rarely annotated despite occurring in 15% of known proteins, and why they tend to be associated with function. An analysis of π-helices in the diverse ferritin-like superfamily illustrates their tendency to be conserved in protein families and identifies a putative π-helix-containing primordial precursor, a “missing link” intermediary form of the ribonucleotide reductase family, vestigial π-helices, and a novel function for π-helices that we term a “peristaltic-like shift.” This new understanding of π-helices paves the way for this generally overlooked motif to become a noteworthy feature that will aid in tracing the evolution of many protein families, guide investigations of protein and π-helix functionality, and contribute additional tools to the protein engineering toolkit.     

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.     

Helix-helix interactions and their impact on protein motifs and assemblies

Protein secondary structure elements are arranged in distinct structural motifs such as four-α-helix bundle, 8α/8β TIM-barrel, Rossmann dinucleotide binding fold, assembly of a helical rod. Each structural motif is characterized by a particular type of helix-helix interactions. A unique pattern of contacts is formed by interacting helices of the structural motif. In each type of fold, edges of the helix surface, which participate in the formation of helix-helix contacts with preceding and following helices, differ. This work shows that circular arrangements of the four, eight, and sixteen α-helices, which are found in the four-α-helical motif, TIM-barrel 8α/8β fold, and helical rod of 16.3¯ helices per turn correspondingly, can be associated with the mutual positioning of the edges of the helix surfaces. Edges (i, i+1)−(i+1, i+2) of the helix surface are central for the interhelical contacts in a four-α-helix bundle. Edges (i, i+1)−(i+2, i+3) are involved in the assembly of four-α-helix subunits into helical rod of a tobacco mosaic virus and a three-helix fragment of a Rossmann fold. In 8α/8β TIM-barrel fold, edges (i, i+1)−(i+5, i+6) are involved in the octagon arrangement. Approximation of a cross section of each motif with a polygon (n-gon, n=4, 8, 16) shows that a good correlation exists between polygon interior angles and angles formed by the edges of helix surfaces.     

Possible roles of exceptionally conserved residues around the selectivity filters of sodium and calcium channels

In the absence of x-ray structures of sodium and calcium channels their homology models are used to rationalize experimental data and design new experiments. A challenge is to model the outer-pore region that folds differently from potassium channels. Here we report a new model of the outer-pore region of the NaV1.4 channel, which suggests roles of highly conserved residues around the selectivity filter. The model takes from our previous study (Tikhonov, D. B., and Zhorov, B. S. (2005) Biophys. J. 88, 184-197) the general disposition of the P-helices, selectivity filter residues, and the outer carboxylates, but proposes new intra- and inter-domain contacts that support structural stability of the outer pore. Glycine residues downstream from the selectivity filter are proposed to participate in knob-into-hole contacts with the P-helices and S6s. These contacts explain the adapted tetrodotoxin resistance of snakes that feed on toxic prey through valine substitution of isoleucine in the P-helix of repeat IV. Polar residues five positions upstream from the selectivity filter residues form H-bonds with the ascending-limb backbones. Exceptionally conserved tryptophans are engaged in inter-repeat H-bonds to form a ring whose π-electrons would facilitate passage of ions from the outer carboxylates to the selectivity filter. The outer-pore model of CaV1.2 derived from the NaV1.4 model is also stabilized by the ring of exceptionally conservative tryptophans and H-bonds between the P-helices and ascending limbs. In this model, the exceptionally conserved aspartate downstream from the selectivity-filter glutamate in repeat II facilitates passage of calcium ions to the selectivity-filter ring through the tryptophan ring. Available experimental data are discussed in view of the models.     

Modeling the ion channel structure of cecropin.

Atomic-scale computer models were developed for how cecropin peptides may assemble in membranes to form two types of ion channels. The models are based on experimental data and physiochemical principles. Initially, cecropin peptides, in a helix-bend-helix motif, were arranged as antiparallel dimers to position conserved residues of adjacent monomers in contact. The dimers were postulated to bind to the membrane with the NH2-terminal helices sunken into the head-group layer and the COOH-terminal helices spanning the hydrophobic core. This causes a thinning of the top lipid layer of the membrane. A collection of the membrane bound dimers were then used to form the type I channel structure, with the pore formed by the transmembrane COOH-terminal helices. Type I channels were then assembled into a hexagonal lattice to explain the large number of peptides that bind to the bacterium. A concerted conformational change of a type I channel leads to the larger type II channel, in which the pore is formed by the NH2-terminal helices. By having the dimers move together, the NH2-terminal helices are inserted into the hydrophobic core without having to desolvate the charged residues. It is also shown how this could bring lipid head-groups into the pore lining.     

Structural understanding of the transmembrane domains of inositol triphosphate receptors and ryanodine receptors towards calcium channeling.

Inositol 1,4,5-triphosphate receptors (Insp(3)Rs) and ryanodine receptors (ryRs) act as cationic channels transporting calcium ions from the endoplasmic reticulum to cytosol by forming tetramers and are proteins localized to the endoplasmic reticulum (ER). Despite the absence of classical calcium-binding motifs, calcium channeling occurs at the transmembrane domain. We have investigated putative calcium binding motifs in these sequences. Prediction methods indicate the presence of six transmembrane helices in the C-terminal domain, one of the three domains conserved between Insp(3)R and ryR receptors. The recently identified crystal structure of the K(+) channel, which also forms tetramers, revealed that two transmembrane helices, an additional pore helix and a selectivity filter are responsible for selective K(+) ion channeling. The last three TM helices of Insp(3)R and ryR are particularly well conserved and we found analogous pore helix and selectivity filter motif in these sequences. We obtained a three-dimensional structural model for the transmembrane tetramer by extrapolating the distant structural similarity to the K(+) channels.     

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.