Making structural sense of dimerization interfaces of delta opioid receptor homodimers

Opioid receptors, like other members of the G protein-coupled receptor (GPCR) family, have been shown to associate to form dimers and/or oligomers at the plasma membrane. Whether this association is stable or transient is not known. Recent compelling evidence suggests that at least some GPCRs rapidly associate and dissociate. We have recently calculated binding affinities from free energy estimates to predict transient association between mouse delta opioid receptor (DOR) protomers at a symmetric interface involving the fourth transmembrane (TM4) helix (herein termed "4" dimer). Here we present disulfide cross-linking experiments with DOR constructs with cysteines substituted at the extracellular ends of TM4 or TM5 that confirm the formation of DOR complexes involving these helices. Our results are consistent with the involvement of TM4 and/or TM5 at the DOR homodimer interface, but possibly with differing association propensities. Coarse-grained (CG) well-tempered metadynamics simulations of two different dimeric arrangements of DOR involving TM4 alone or with TM5 (herein termed "4/5" dimer) in an explicit lipid-water environment confirmed the presence of two structurally and energetically similar configurations of the 4 dimer, as previously assessed by umbrella sampling calculations, and revealed a single energetic minimum of the 4/5 dimer. Additional CG umbrella sampling simulations of the 4/5 dimer indicated that the strength of association between DOR protomers varies depending on the protein region at the interface, with the 4 dimer being more stable than the 4/5 dimer.     

G-protein-coupled receptor dynamics: dimerization and activation models compared with experiment

Our previously derived models of the active state of the β2-adrenergic receptor are compared with recently published X-ray crystallographic structures of activated GPCRs (G-protein-coupled receptors). These molecular dynamics-based models using experimental data derived from biophysical experiments on activation were used to restrain the receptor to an active state that gave high enrichment for agonists in virtual screening. The β2-adrenergic receptor active model and X-ray structures are in good agreement over both the transmembrane region and the orthosteric binding site, although in some regions the active model is more similar to the active rhodopsin X-ray structures. The general features of the microswitches were well reproduced, but with minor differences, partly because of the unexpected X-ray results for the rotamer toggle switch. In addition, most of the interacting residues between the receptor and the G-protein were identified. This analysis of the modelling has also given important additional insight into GPCR dimerization: re-analysis of results on photoaffinity analogues of rhodopsin provided additional evidence that TM4 (transmembrane helix 4) resides at the dimer interface and that ligands such as bivalent ligands may pass between the mobile helices. A comparison, and discussion, is also carried out between the use of implicit and explicit solvent for active-state modelling.     

Molecular Insights into the Transmembrane Domain of the Thyrotropin Receptor

The thyrotropin receptor (TSHR) is a G protein-coupled receptor (GPCR) that is member of the leucine-rich repeat subfamily (LGR). In the absence of crystal structure, the success of rational design of ligands targeting the receptor internal cavity depends on the quality of the TSHR models built. In this subfamily, transmembrane helices (TM) 2 and 5 are characterized by the absence of proline compared to most receptors, raising the question of the structural conformation of these helices. To gain insight into the structural properties of these helices, we carried out bioinformatics and experimental studies. Evolutionary analysis of the LGR family revealed a deletion in TM5 but provided no information on TM2. Wild type residues at positions 2.58, 2.59 or 2.60 in TM2 and/or at position 5.50 in TM5 were substituted to proline. Depending on the position of the proline substitution, different effects were observed on membrane expression, glycosylation, constitutive cAMP activity and responses to thyrotropin. Only proline substitution at position 2.59 maintained complex glycosylation and high membrane expression, supporting occurrence of a bulged TM2. The TSHR transmembrane domain was modeled by homology with the orexin 2 receptor, using a protocol that forced the deletion of one residue in the TM5 bulge of the template. The stability of the model was assessed by molecular dynamics simulations. TM5 straightened during the equilibration phase and was stable for the remainder of the simulations. Our data support a structural model of the TSHR transmembrane domain with a bulged TM2 and a straight TM5 that is specific of glycoprotein hormone receptors.     

A helix heterodimer in a lipid bilayer: prediction of the structure of an integrin transmembrane domain via multiscale simulations

Dimerization of transmembrane (TM) α helices of membrane receptors plays a key role in signaling. We show that molecular dynamics simulations yield models of integrin TM helix heterodimers, which agree well with available NMR structures. We use?a multiscale simulation approach, combining coarse-grained and subsequent atomistic simulation, to model the dimerization of wild-type (WT) and mutated sequences of the αIIb and β3 integrin TM helices. The WT helices formed a stable, right-handed dimer with the same helix-helix interface as in the published NMR structure (PDB: 2K9J). In contrast, the presence of disruptive mutations perturbed the interface between the helices, altering the conformational stability of the dimer. The αIIb/β3 interface was more flexible than that of, e.g., glycophorin A. This is suggestive of a role for alternative packing modes of the TM helices in transbilayer signaling.     

Cholesterol Modulates the Dimer Interface of the β2-Adrenergic Receptor via Cholesterol Occupancy Sites

The β₂-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β₂-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β₂-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.     

Cholesterol Modulates the Dimer Interface of the β2-Adrenergic Receptor via Cholesterol Occupancy Sites

The β₂-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β₂-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β₂-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.     

Stabilization of the human beta2-adrenergic receptor TM4-TM3-TM5 helix interface by mutagenesis of Glu122(3.41), a critical residue in GPCR structure

G protein-coupled receptor (GPCR) instability represents one of the most profound obstacles in the structural study of GPCRs that bind diffusible ligands. The introduction of targeted mutations at nonconserved residues that lie proximal to helix interfaces has the potential to enhance the fold stability of the receptor helix bundle while maintaining wild-type receptor function. To test this hypothesis, we studied the effect of amino acid substitutions at Glu122^3.41 in the well-studied β₂-adrenergic receptor (β₂AR), which was predicted from sequence conservation to lie at a position equivalent to a tryptophan residue in rhodopsin at the 3,4,5 helix interface among transmembrane (TM) domains 3, 4, and 5. Replacement of Glu122^3.41 with bulky hydrophobic residues, such as tryptophan, tyrosine, and phenylalanine, increases the yield of functionally folded β₂AR by as much as 5-fold. Receptor stability in detergent solution was studied by isothermal denaturation, and it was found that the E122W and E122Y mutations enhanced the β₂AR thermal half-life by 9.3- and 6.7-fold, respectively, at 37 °C. The β₁AR was also stabilized by the introduction of tryptophan at Glu147^3.41, and the effect on protein behavior was similar to the rescue of the unstable wild-type receptor by the antagonist propranolol. Molecular modeling of the E122W and E122Y mutants revealed that the tryptophan ring edge and tyrosine hydroxyl are positioned proximal to the helical break in TM5 introduced by the conserved Pro211^5.50 and may stabilize the helix by interacting favorably with the unpaired carbonyl oxygen of Val206^5.45. Conformational flexibility of TM5 is likely to be a general property of class A GPCRs; therefore, engineering of the TM4-TM3-TM5 interface at the 3.41 position may provide a general strategy for the stabilization of other receptors.     

The extreme C-terminal region of Gαs differentially couples to the luteinizing hormone and beta2-adrenergic receptors

The mechanisms of G protein coupling to G protein-coupled receptors (GPCR) share general characteristics but may exhibit specific interactions unique for each GPCR/G protein partnership. The extreme C terminus (CT) of G protein α-subunits has been shown to be important for association with GPCR. Hypothesizing that the extreme CT of Gα(s) is an essential component of the molecular landscape of the GPCR, human LH receptor (LHR), and β(2)-adrenergic receptor (β(2)-AR), a model cell system was created for the expression and manipulation of Gα(s) subunits in LHR(+) s49 ck cells that lack endogenous Gα(s). On the basis of studies involving truncations, mutations, and chain extensions of Gα(s), the CT was found to be necessary for LHR and β(2)-AR signaling. Some general similarities were found for the responses of the two receptors, but significant differences were also noted. Computational modeling was performed with a combination of comparative modeling, molecular dynamics simulations, and rigid body docking. The resulting models, focused on the Gα(s) CT, are supported by the experimental observations and are characterized by the interaction of the four extreme CT amino acid residues of Gα(s) with residues in LHR and β(2)-AR helix 3, (including R of the DRY motif), helix 6, and intracellular loop 2. This portion of Gα(s) recognizes the same regions of the two GPCR, although with differences in the details of selected interactions. The predicted longer cytosolic extensions of helices 5 and 6 of β(2)-AR are expected to contribute significantly to differences in Gα(s) recognition by the two receptors.     

The extracellular N-terminal domain and transmembrane domains 1 and 2 mediate oligomerization of a yeast G protein-coupled receptor

G protein-coupled receptors (GPCRs) can form homodimers/oligomers and/or heterodimers/oligomers. The mechanisms used to form specific GPCR oligomers are poorly understood because the domains that mediate such interactions and the step(s) in the secretory pathway where oligomerization occurs have not been well characterized. Here we have used subcellular fractionation and fluorescence resonance energy transfer (FRET) experiments to show that oligomerization of a GPCR (alpha-factor receptor; STE2 gene product) of the yeast Saccharomyces cerevisiae occurs in the endoplasmic reticulum. To identify domains of this receptor that mediate oligomerization, we used FRET and endocytosis assays of oligomerization in vivo to analyze receptor deletion mutants. A mutant lacking the N-terminal extracellular domain and transmembrane (TM) domain 1 was expressed at the cell surface but did not self-associate. In contrast, a receptor fragment containing only the N-terminal extracellular domain and TM1 could self-associate and heterodimerize with wild type receptors. Analysis of other mutants suggested that oligomerization is facilitated by the N-terminal extracellular domain and TM2. Therefore, the N-terminal extracellular domain, TM1, and TM2 appear to stabilize alpha-factor receptor oligomers. These domains may form an interface in contact or domain-swapped oligomers. Similar domains may mediate dimerization of certain mammalian GPCRs.     

Characterization of the homodimerization interface and functional hotspots of the CXCR4 chemokine receptor

The recent crystallographic structures of the human chemokine CXC Receptor 4 (CXCR4) provide experimental evidence of a human G Protein-Coupled Receptor (GPCR) dimer in atomic detail. The CXCR4 homodimers reveal an unexpected dimerization mode involving transmembrane helices TM5 and TM6, which is examined here using all-atom molecular dynamics (MD) simulations in the physiological environment of a lipid bilayer. The bacteriophage T4 lysozyme (T4L), which was fused to the crystallized protein but absent in our simulations, is found to slightly affect the observed relative position of the protomers in the two dimers studied here, and consequently some rearrangements of the dimerization interface are proposed. In addition, the simulations provide further evidence about the role of the two stabilizing single point mutations introduced to crystallize the receptor. Finally, this work analyzes the structural and dynamic role of key residues involved both in ligand binding and in the infection process of HIV. In particular, the different side chain conformations of His113(3.39) are found to influence the dynamics of the surrounding functional hotspot region being evaluated both in the presence and in the absence of the co-crystallized ligand IT1t. The analysis reported here adds valuable knowledge for future structure-based drug design (SBDD) efforts on this pharmacological target.     

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.     

Differential regulation of two palmitoylation sites in the cytoplasmic tail of the beta1-adrenergic receptor

S-Palmitoylation of G protein-coupled receptors (GPCRs) is a prevalent modification, contributing to the regulation of receptor function. Despite its importance, the palmitoylation status of the β(1)-adrenergic receptor, a GPCR critical for heart function, has never been determined. We report here that the β(1)-adrenergic receptor is palmitoylated on three cysteine residues at two sites in the C-terminal tail. One site (proximal) is adjacent to the seventh transmembrane domain and is a consensus site for GPCRs, and the other (distal) is downstream. These sites are modified in different cellular compartments, and the distal palmitoylation site contributes to efficient internalization of the receptor following agonist stimulation. Using a bioorthogonal palmitate reporter to quantify palmitoylation accurately, we found that the rates of palmitate turnover at each site are dramatically different. Although palmitoylation at the proximal site is remarkably stable, palmitoylation at the distal site is rapidly turned over. This is the first report documenting differential dynamics of palmitoylation sites in a GPCR. Our results have important implications for function and regulation of the clinically important β(1)-adrenergic receptor.     

Disruption of Rhodopsin Dimerization with Synthetic Peptides Targeting an Interaction Interface

Although homo- and heterodimerizations of G protein-coupled receptors (GPCRs) are well documented, GPCR monomers are able to assemble in different ways, thus causing variations in the interactive interface between receptor monomers among different GPCRs. Moreover, the functional consequences of this phenomenon, which remain to be clarified, could be specific for different GPCRs. Synthetic peptides derived from transmembrane (TM) domains can interact with a full-length GPCR, blocking dimer formation and affecting its function. Here we used peptides corresponding to TM helices of bovine rhodopsin (Rho) to investigate the Rho dimer interface and functional consequences of its disruption. Incubation of Rho with TM1, TM2, TM4, and TM5 peptides in rod outer segment (ROS) membranes shifted the resulting detergent-solubilized protein migration through a gel filtration column toward smaller molecular masses with a reduced propensity for dimer formation in a cross-linking reaction. Binding of these TM peptides to Rho was characterized by both mass spectrometry and a label-free assay from which dissociation constants were calculated. A BRET (bioluminescence resonance energy transfer) assay revealed that the physical interaction between Rho molecules expressed in membranes of living cells was blocked by the same four TM peptides identified in our in vitro experiments. Although disruption of the Rho dimer/oligomer had no effect on the rates of G protein activation, binding of G_t to the activated receptor stabilized the dimer. However, TM peptide-induced disruption of dimer/oligomer decreased receptor stability, suggesting that Rho supramolecular organization could be essential for ROS stabilization and receptor trafficking.     

Structural models for dimerization of G-protein coupled receptors: the opioid receptor homodimers

Among the most exciting functional features of G-protein coupled receptors (GPCRs) that are coming into focus lately are those relating to the role and structural characteristics of their oligomerization (mostly homo- and heterodimers). The structural underpinnings of these novel functional insights are still not clear, as current experimental techniques have not yet succeeded in identifying the dimerization interfaces between GPCR monomers. Two computational approaches have recently been designed in our lab to provide reasonable three-dimensional (3D) molecular models of the transmembrane (TM) regions of GPCR dimers based on a combination of the structural information of receptor monomers and analyses of correlated mutations in receptor families. The modeling of GPCR heterodimers has been described recently. We present here a related approach for modeling of GPCR homodimers that identifies the interfaces in the most likely configurations of the complexes. The approach is illustrated for the three cloned opioid receptor subtypes (OPRD, OPRM, and OPRK).     

Modeling GPCR active state conformations: the β(2)-adrenergic receptor

The recent publication of several G protein-coupled receptor (GPCR) structures has increased the information available for homology modeling inactive class A GPCRs. Moreover, the opsin crystal structure shows some active features. We have therefore combined information from these two sources to generate an extensively validated model of the active conformation of the β(2)-adrenergic receptor. Experimental information on fully active GPCRs from zinc binding studies, site-directed spin labeling, and other spectroscopic techniques has been used in molecular dynamics simulations. The observed conformational changes reside mainly in transmembrane helix 6 (TM6), with additional small but significant changes in TM5 and TM7. The active model has been validated by manual docking and is in agreement with a large amount of experimental work, including site-directed mutagenesis information. Virtual screening experiments show that the models are selective for β-adrenergic agonists over other GPCR ligands, for (R)- over (S)-β-hydroxy agonists and for β(2)-selective agonists over β(1)-selective agonists. The virtual screens reproduce interactions similar to those generated by manual docking. The C-terminal peptide from a model of the stimulatory G protein, readily docks into the active model in a similar manner to which the C-terminal peptide from transducin, docks into opsin, as shown in a recent opsin crystal structure. This GPCR-G protein model has been used to explain site-directed mutagenesis data on activation. The agreement with experiment suggests a robust model of an active state of the β(2)-adrenergic receptor has been produced. The methodology used here should be transferable to modeling the active state of other GPCRs.     

Dopamine D2 receptors form higher order oligomers at physiological expression levels

G-protein-coupled receptors are generally thought to be organized as dimers; whether they form higher order oligomers is a topic of much controversy. We combined bioluminescence/fluorescence complementation and energy transfer to demonstrate that at least four dopamine D2 receptors are located in close molecular proximity in living mammalian cells, consistent with their organization as higher order oligomers at the plasma membrane. This implies the existence of multiple receptor interfaces. In addition to the symmetrical interface in the fourth transmembrane segment (TM4) we identified previously by cysteine (Cys) crosslinking, we now show that a patch of residues at the extracellular end of TM1 forms a second symmetrical interface. Crosslinking of D2 receptor with Cys substituted simultaneously into both TM1 and TM4 led to higher order species, consistent with our novel biophysical results. Remarkably, the rate and extent of crosslinking at both interfaces were unaltered over a 100-fold range of receptor expression. Thus, at physiological levels of expression, the receptor is organized in the plasma membrane into a higher order oligomeric structure.     

Evidence that specific interactions play a role in the cholesterol sensitivity of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are known to be modulated by membrane cholesterol levels, but whether or not the effects are caused by specific receptor-cholesterol interactions or cholesterol's general effects on the membrane is not well-understood. We performed coarse-grained molecular dynamics (CGMD) simulations coupled with structural bioinformatics approaches on the β₂-adrenergic receptor (β₂AR) and the cholecystokinin (CCK) receptor subfamily. The β₂AR has been shown to be sensitive to membrane cholesterol and cholesterol molecules have been clearly resolved in numerous β₂AR crystal structures. The two CCK receptors are highly homologous and preserve similar cholesterol recognition motifs but despite their homology, CCK₁R shows functional sensitivity to membrane cholesterol while CCK₂R does not. Our results offer new insights into how cholesterol modulates GPCR function by showing cholesterol interactions with β₂AR that agree with previously published data; additionally, we observe differential and specific cholesterol binding in the CCK receptor subfamily while revealing a previously unreported Cholesterol Recognition Amino-acid Consensus (CRAC) sequence that is also conserved across 38% of class A GPCRs. A thermal denaturation assay (LCP-T_m) shows that mutation of a conserved CRAC sequence on TM7 of the β₂AR affects cholesterol stabilization of the receptor in a lipid bilayer. The results of this study provide a better understanding of receptor-cholesterol interactions that can contribute to novel and improved therapeutics for a variety of diseases.     

Ivermectin Interaction with transmembrane helices reveals widespread rearrangements during opening of P2X receptor channels

P2X receptors are trimeric cation channels that open in response to binding of extracellular ATP. Each subunit contains a large extracellular ligand binding domain and two flanking transmembrane (TM) helices that form the pore, but the extent of gating motions of the TM helices is unclear. We probed these motions using ivermectin (IVM), a macrocyclic lactone that stabilizes the open state of P2X(4) receptor channels. We find that IVM partitions into lipid membranes and that transfer of the TM regions of P2X(4) receptors is sufficient to convey sensitivity to the lactone, suggesting that IVM interacts most favorably with the open conformation of the two TM helices at the protein-lipid interface. Scanning mutagenesis of the two TMs identifies residues that change environment between closed and open states, and substitutions at a subset of these positions weaken IVM binding. The emerging patterns point to widespread rearrangements of the TM helices during opening of P2X receptor channels.     

Identifying interactions between transmembrane helices from the adenosine A2A receptor

The human adenosine A(2A) receptor (A(2A)R) is an integral membrane protein and a member of the G-protein-coupled receptor (GPCR) superfamily, characterized by seven transmembrane (TM) helices. Although helix-helix association in the lipid bilayer is known to be an essential step in the folding of GPCRs, the determinants of their structures, folding, and assembly in the cell membrane are poorly understood. Previous studies in our group showed that while peptides corresponding to all seven TM domains of A(2A)R form stable helical structures in detergent micelles and lipid vesicles, they display significant variability in their helical propensity. This finding suggested to us that some TM domains might need to interact with other domains to properly insert and fold in hydrophobic environments. In this study, we assessed the ability of TM peptides to interact in pairwise combinations. We analyzed peptide interactions in hydrophobic milieus using circular dichroism spectroscopy and F?rster resonance energy transfer. We find that specific interactions between TM helices occur, leading to additional helical content, especially in weakly helical TM domains, suggesting that some TM domains need a partner for proper folding in the membrane. The approach developed in this study will enable complete analysis of the TM domain interactions and the modeling of a folding pathway for A(2A)R.