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.     

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.     

TM-MOTIF: an alignment viewer to annotate predicted transmembrane helices and conserved motifs in aligned set of sequences

Multiple sequence alignments become biologically meaningful only if conserved and functionally important residues and secondary structural elements preserved can be identified at equivalent positions. This is particularly important for transmembrane proteins like G-protein coupled receptors (GPCRs) with seven transmembrane helices. TM-MOTIF is a software package and an effective alignment viewer to identify and display conserved motifs and amino acid substitutions (AAS) at each position of the aligned set of homologous sequences of GPCRs. The key feature of the package is to display the predicted membrane topology for seven transmembrane helices in seven colours (VIBGYOR colouring scheme) and to map the identified motifs on its respective helices /loop regions. It is an interactive package which provides options to the user to submit query or pre-aligned set of GPCR sequences to align with a reference sequence, like rhodopsin, whose structure has been solved experimentally. It also provides the possibility to identify the nearest homologue from the available inbuilt GPCR or Olfactory Receptor cluster dataset whose association is already known for its receptor type. AVAILABILITY: The database is available for free at mini@ncbs.res.in.     

Simulation studies on bacteriorhodopsin bundle of transmembrane α segments

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Seven transmembrane BR helical segments are subjected to simulation studies in order to investigate the packing process of transmembrane helices. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Helix packing is optimized according to a semi-empirical potential mainly composed of six components: a bilayer potential, a crossing angle potential, a helix dipole potential, a helix-helix distance potential, a helix orientation potential and a helix-helix distance restraint potential (a loop potential). Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. The structures from the simulations are compared with the experimentally determined structures in terms of geometry. The structures generated show similar shapes to the experimentally suggested structure even without the helix-helix distance restraint potential. However, the relative locations of individual helices were reproduced only when the helix-helix distance restraint potential was used with restraint conditions. Our results suggest that transmembrane helix bundles resembling those observed experimentally may be generated by simulations using simple potentials. Received: 19 April 1999 / Revised version: 6 September 1999 / Accepted: 17 September 1999     

Simulation studies on bacteriorhodopsin α-helices

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Transmembrane BR helical segments are subjected to simulation studies in order to investigate the effect of bilayer environment in various simulation conditions. A bilayer potential is introduced to the system to mimic the lipid membrane. The structures from the simulations are compared with the experimentally determined structures in terms of geometrical properties. Electrostatic contribution to the helix packing is also investigated. The simulation results show that the packing geometry of the transmembrane helices is highly affected by the bilayer potential. The results obtained from the simulations may be used for further simulation studies and analysis in investigating transmembrane helix packing. Received: 5 July 1999 / Revised version: 5 October 1999 / Accepted: 15 October 1999     

Statistical sequence analyses of G-protein-coupled receptors: structural and functional characteristics viewed with periodicities of entropy, hydrophobicity, and volume

G-protein-coupled receptors (GPCRs) play a crucial role in signal transduction and receive a wide variety of ligands. GPCRs are a major target in drug design, as nearly 50% of all contemporary medicines act on GPCRs. GPCRs are membrane proteins possessing a common structural feature, seven transmembrane helices. In order to design an effective drug to act on a GPCR, knowledge of the three-dimensional (3D) structure of the target GPCR is indispensable. However, as GPCRs are membrane bound, their 3D structures are difficult to obtain. Thus we conducted statistical sequence analyses to find information about 3D structure and ligand binding using the receptors' primary sequences. We present statistical sequence analyses of 270 human GPCRs with regard to entropy (Shannon entropy in sequence alignment), hydrophobicity and volume, which are associated with the alpha-helical periodicity of the accessibility to the surrounding lipid. We found periodicity such that the phase changes once in the middle of each transmembrane region, both in the entropy plot and in the hydrophobicity plot. The phase shift in the entropy plot reflects the variety of ligands and the generality of the mechanism of signal transduction. The two periodic regions in the hydrophobicity plot indicate the regions facing the hydrophobic lipid chain and the polar phospholipid headgroup. We also found a simple periodicity in the plot of volume deviation, which suggests conservation of the stable structural packing among the transmembrane helices.     

Role of conserved prolines in the structure and function of the Na+/dicarboxylate cotransporter 1, NaDC1

The Na+/dicarboxylate cotransporter 1 (NaDC1) is a low-affinity transporter for citric acid cycle intermediates such as succinate and citrate. The sequence of NaDC1 contains a number of conserved proline residues in predicted transmembrane helices (TMs) 7 and 10. These transmembrane domains are of particular importance because they may be involved in determining the substrate or cation-binding affinity in NaDC1. Four conserved proline residues in TMs 7 and 10 of rabbit NaDC1 were replaced with alanine to promote ideal alpha helix or glycine to promote free conformation, and the mutant transporters were expressed in the HRPE cell line. Mutations of prolines in TM 10 produced decreased protein expression and activity, whereas mutations of prolines in TM 7 completely abolished protein expression and activity. The chemical chaperone glycerol was found to improve the expression of the Pro-351 mutants in TM 7, suggesting that these mutants had defects in trafficking. The inactive mutant transporters at position 351 could also be rescued by the addition of a proline at a second site. For example, the P351A-F347P mutant had restored activity, although its substrate specificity was altered. We conclude that, in TM 7, Pro-327 may be of particular importance in the function of the transporter, whereas Pro-351 may affect protein targeting. The prolines in TM 10, at positions 523 and 524, may not be directly involved in the transporter function but may be necessary for maintaining structure.     

Proline residues in transmembrane helices of channel and transport proteins: a molecular modelling study.

Proline residues are commonly found in putative transbilayer helices of many integral membrane proteins which act as transporters, channels and receptors. Intramembranous prolines are often conserved between homologous proteins. It has been suggested that such intrahelical prolines provide liganding sites for cations via exposure of the backbone carbonyl oxygen atoms of residues i-3 and i-4 (relative to the proline). Molecular modelling studies have been carried out to evaluate this proposal. Bundles of parallel proline-kinked helices are considered as simplified models of ion channels. The energetics of K+ ion-helix bundle interactions are explored. It is shown that carbonyl oxygens exposed by the proline-induced kink and at the C-terminus of the helices may provide cation-liganding sites. 'Hybrid' bundles of antiparallel helices, only some of which contain proline residues, are considered as models of transport proteins. Again, proline-exposed carbonyl oxygens are shown to be capable of liganding cations. The roles of alpha-helix dipoles and of the geometry of helix packing are considered in relation to cation-bundle interactions. Implications with respect to modelling of ion channel and transport proteins are discussed.     

Comparison of class A and D G protein-coupled receptors: common features in structure and activation

All G protein-coupled receptors (GPCRs) share a common seven TM helix architecture and the ability to activate heterotrimeric G proteins. Nevertheless, these receptors have widely divergent sequences with no significant homology. We present a detailed structure-function comparison of the very divergent Class A and D receptors to address whether there is a common activation mechanism across the GPCR superfamily. The Class A and D receptors are represented by the vertebrate visual pigment rhodopsin and the yeast alpha-factor pheromone receptor Ste2, respectively. Conserved amino acids within each specific receptor class and amino acids where mutation alters receptor function were located in the structures of rhodopsin and Ste2 to assess whether there are functionally equivalent positions or regions within these receptors. We find several general similarities that are quite striking. First, strongly polar amino acids mediate helix interactions. Their mutation generally leads to loss of function or constitutive activity. Second, small and weakly polar amino acids facilitate tight helix packing. Third, proline is essential at similar positions in transmembrane helices 6 and 7 of both receptors. Mapping the specific location of the conserved amino acids and sites of constitutively active mutations identified conserved microdomains on transmembrane helices H3, H6, and H7, suggesting that there are underlying similarities in the mechanism of the widely divergent Class A and Class D receptors.     

Why are polar residues within the membrane core evolutionary conserved?

Here, we present a study of polar residues within the membrane core of alpha-helical membrane proteins. As expected, polar residues are less frequent in the membrane than expected. Further, most of these residues are buried within the interior of the protein and are only rarely exposed to lipids. However, the polar groups often border internal water filled cavities, even if the rest of the sidechain is buried. A survey of their functional roles in known structures showed that the polar residues are often directly involved in binding of small compounds, especially in channels and transporters, but other functions including proton transfer, catalysis, and selectivity have also been attributed to these proteins. Among the polar residues histidines often interact with prosthetic groups in photosynthetic- and oxidoreductase-related proteins, whereas prolines often are required for conformational changes of the proteins. Indeed, the polar residues in the membrane core are more conserved than other residues in the core, as well as more conserved than polar residues outside the membrane. The reason is twofold; they are often (i) buried in the interior of the protein and (ii) directly involved in the function of the proteins. Finally, a method to identify which polar residues are present within the membrane core directly from protein sequences was developed. Applying the method to the set of all human membrane proteins the prediction indicates that polar residues were most frequent among active transporter proteins and GPCRs, whereas infrequent in families with few transmembrane regions, such as non-GPCR receptors.     

Proline-induced distortions of transmembrane helices

Proline residues in the transmembrane (TM) alpha-helices of integral membrane proteins have long been suspected to play a key role for helix packing and signal transduction by inducing regions of helix distortion and/or dynamic flexibility (hinges). In this study we try to characterise the effect of proline on the geometric properties of TM alpha-helices. We have examined 199 transmembrane alpha-helices from polytopic membrane proteins of known structure. After examining the location of proline residues within the amino acid sequences of TM helices, we estimated the helix axes either side of a hinge and hence identified a hinge residue. This enabled us to calculate helix kink and swivel angles. The results of this analysis show that proline residues occur with a significant concentration in the centre of sequences of TM alpha-helices. In this location, they may induce formation of molecular hinges, located on average about four residues N-terminal to the proline residue. A superposition of proline-containing TM helices structures shows that the distortion induced is anisotropic and favours certain relative orientations (defined by helix kink and swivel angles) of the two helix segments.     

Molecular dynamics simulations reveal insights into key structural elements of adenosine receptors

The crystal structure of the human A(2A) adenosine receptor, a member of the G protein-coupled receptor (GPCR) family, is used as a starting point for the structural characterization of the conformational equilibrium around the inactive conformation of the human A(2) (A(2A) and A(2B)) adenosine receptors (ARs). A homology model of the closely related A(2B)AR is reported, and the two receptors were simulated in their apo form through all-atom molecular dynamics (MD) simulations. Different conditions were additionally explored in the A(2A)AR, including the protonation state of crucial histidines or the presence of the cocrystallized ligand. Our simulations reveal the role of several conserved residues in the ARs in the conformational equilibrium of the receptors. The "ionic lock" absent in the crystal structure of the inactive A(2A)AR is rapidly formed in the two simulated receptors, and a complex network of interacting residues is presented that further stabilizes this structural element. Notably, the observed rotameric transition of Trp6.48 ("toggle switch"), which is thought to initiate the activation process in GPCRs, is accompanied by a concerted rotation of the conserved residue of the A(2)ARs, His6.52. This new conformation is further stabilized in the two receptors under study by a novel interaction network involving residues in transmembrane (TM) helices TM5 (Asn5.42) and TM3 (Gln3.37), which resemble the conformational changes recently observed in the agonist-bound structure of β-adrenoreceptors. Finally, the interaction between Glu1.39 and His7.43, a pair of conserved residues in the family of ARs, is found to be weaker than previously thought, and the role of this interaction in the structure and dynamics of the receptor is thoroughly examined. All these findings suggest that, despite the commonalities with other GPCRs, the conformational equilibrium of ARs is also modulated by specific residues of the family.     

Dimerization in aminergic G-protein-coupled receptors: application of a hidden-site class model of evolution

G-Protein-coupled receptors (GPCRs) are an important superfamily of transmembrane proteins involved in cellular communication. Recently, it has been shown that dimerization is a widely occurring phenomenon in the GPCR superfamily, with likely important physiological roles. Here we use a novel hidden-site class model of evolution as a sequence analysis tool to predict possible dimerization interfaces in GPCRs. This model aims to simulate the evolution of proteins at the amino acid level, allowing the analysis of their sequences in an explicitly evolutionary context. Applying this model to aminergic GPCR sequences, we first validate the general reasoning behind the model. We then use the model to perform a family specific analysis of GPCRs. Accounting for the family structure of these proteins, this approach detects different evolutionarily conserved and accessible patches on transmembrane (TM) helices 4-6 in different families. On the basis of these findings, we propose an experimentally testable dimerization mechanism, involving interactions among different combinations of these helices in different families of aminergic GPCRs.     

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.     

NMR structure of the second intracellular loop of the alpha 2A adrenergic receptor: evidence for a novel cytoplasmic helix

A major, unresolved question in signal transduction by G protein coupled receptors (GPCRs) is to understand how, at atomic resolution, a GPCR activates a G protein. A step toward answering this question was made with the determination of the high-resolution structure of rhodopsin; we now know the intramolecular interactions that characterize the resting conformation of a GPCR. To what degree does this structure represent a structural paradigm for other GPCRs, especially at the cytoplasmic surface where GPCR-G protein interaction occurs and where the sequence homology is low among GPCRs? To address this question, we performed NMR studies on approximately 35-residue-long peptides including the critical second intracellular loop (i2) of the alpha 2A adrenergic receptor (AR) and of rhodopsin. To stabilize the secondary structure of the peptide termini, 4-12 residues from the adjacent transmembrane helices were included and structures determined in dodecylphosphocholine micelles. We also characterized the effects on an alpha 2A AR peptide of a D130I mutation in the conserved DRY motif. Our results show that in contrast to the L-shaped loop in the i2 of rhodopsin, the i2 of the alpha 2A AR is predominantly helical, supporting the hypothesis that there is structural diversity within GPCR intracellular loops. The D130I mutation subtly modulates the helical structure. The spacing of nonpolar residues in i2 with helical periodicity is a predictor of helical versus loop structure. These data should lead to more accurate models of the intracellular surface of GPCRs and of receptor-mediated G protein activation.     

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.     

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.     

Prediction of the odorant binding site of olfactory receptor proteins by human-mouse comparisons

Olfactory receptors (ORs) are a large family of proteins involved in the recognition and discrimination of numerous odorants. These receptors belong to the G-protein coupled receptor (GPCR) hyperfamily, for which little structural data are available. In this study we predict the binding site residues of OR proteins by analyzing a set of 1441 OR protein sequences from mouse and human. The central insight utilized is that functional contact residues would be conserved among pairs of orthologous receptors, but considerably less conserved among paralogous pairs. Using judiciously selected subsets of 218 ortholog pairs and 518 paralog pairs, we have identified 22 sequence positions that are both highly conserved among the putative orthologs and variable among paralogs. These residues are disposed on transmembrane helices 2 to 7, and on the second extracellular loop of the receptor. Strikingly, although the prediction makes no assumption about the location of the binding site, these amino acid positions are clustered around a pocket in a structural homology model of ORs, mostly facing the inner lumen. We propose that the identified positions constitute the odorant binding site. This conclusion is supported by the observation that all but one of the predicted binding site residues correspond to ligand-contact positions in other rhodopsin-like GPCRs.     

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.     

Residues from transmembrane helices 3 and 5 participate in leukotriene B4 binding to BLT1

Leukotrienes are inflammatory mediators that bind to seven transmembrane, G-protein-coupled receptors (GPCRs). Here we examine residues from transmembrane helices 3 and 5 of the leukotriene B4 (LTB4) receptor BLT1 to elucidate how these residues are involved in ligand binding. We have selected these residues on the basis of (1) amino acid sequence analysis, (2) receptor binding and activation studies with a variety of leukotriene-like ligands and recombinant BLT1 receptors, (3) previously published recombinant BLT1 mutants, and (4) a computed model of the active structure of the BLT1 receptor. We propose that LTB4 binds with the polar carboxylate group of LTB4 near the extracellular surface of BLT1 and with the hydrophobic LTB4 tail pointing into the transmembrane regions of the receptor protein. The carboxylate group and the two hydroxyls of LTB4 interact with Arg178 and Glu185 in transmembrane helix 5. Residues from transmembrane helix 3, Val105 and Ile108, also line the pocket deeper inside the receptor. LTB4 is becoming increasingly important as an immunomodulator during a number of pathologies, including atherosclerosis. Detailed information about the LTB4 binding mechanism, and the receptor residues involved, will hopefully aid in the design of new immunomodulatory drugs.