Crosstalk in G protein-coupled receptors: changes at the transmembrane homodimer interface determine activation

Multiple Ecological Axes Drive Molecular Evolution of Cone Opsins in Beloniform Fishes

Ecological and evolutionary transitions offer an excellent opportunity to examine the molecular basis of adaptation. Fishes of the order Beloniformes include needlefishes, flyingfishes, halfbeaks, and allies, and comprise over 200 species occupying a wide array of habitats-from the marine epipelagic zone to tropical rainforest rivers. These fishes also exhibit a diversity of diets, including piscivory, herbivory, and zooplanktivory. We investigated how diet and habitat affected the molecular evolution of cone opsins, which play a key role in bright light and colour vision and are tightly linked to ecology and life history. We analyzed a targeted-capture dataset to reconstruct the evolutionary history of beloniforms and assemble cone opsin sequences. We implemented codon-based clade models of evolution to examine how molecular evolution was affected by habitat and diet. We found high levels of positive selection in medium- and long-wavelength beloniform opsins, with piscivores showing increased positive selection in medium-wavelength opsins and zooplanktivores showing increased positive selection in long-wavelength opsins. In contrast, short-wavelength opsins showed purifying selection. While marine/freshwater habitat transitions have an effect on opsin molecular evolution, we found that diet plays a more important role. Our study suggests that evolutionary transitions along ecological axes produce complex adaptive interactions that affect patterns of selection on genes that underlie vision.

The World of GPCR dimers - Mapping dopamine receptor D2 homodimers in different activation states and configuration arrangements

G protein-coupled receptors (GPCRs) are known to dimerize, but the molecular and structural basis of GPCR dimers is not well understood. In this study, we developed a computational framework to generate models of symmetric and asymmetric GPCR dimers using different monomer activation states and identified their most likely interfaces with molecular details. We chose the dopamine receptor D2 (D2R) homodimer as a case study because of its biological relevance and the availability of structural information. Our results showed that transmembrane domains 4 and 5 (TM4 and TM5) are mostly found at the dimer interface of the D2R dimer and that these interfaces have a subset of key residues that are mostly nonpolar from TM4 and TM5, which was in line with experimental studies. In addition, TM2 and TM3 appear to be relevant for D2R dimers. In some cases, the inactive configuration is unaffected by the partnered protomer, whereas in others, the active protomer adopts the properties of an inactive receptor. Additionally, the β-arrestin configuration displayed the properties of an active receptor in the absence of an agonist, suggesting that a switch to another meta-state during dimerization occurred. Our findings are consistent with the experimental data, and this method can be adapted to study heterodimers and potentially extended to include additional proteins such as G proteins or β-arrestins. In summary, this approach provides insight into the impact of the conformational status of partnered protomers on the overall quaternary GPCR macromolecular structure and dynamics.

The transmembrane domains of GPCR dimers as targets for drug development

G-protein-coupled receptors (GPCRs) can form homodimers or heterodimers that modulate specific signal transduction pathways to regulate a wide range of physiological and pathological functions. As such, GPCR dimers are novel drug targets for disorders including depression, hypertension, diabetes, and vascular dementia. The interaction between two receptors in a GPCR dimer involves a conformational change in the transmembrane domain (TMD). It has been demonstrated that the TMD has an important role in GPCR dimer formation and stability in vitro and in vivo. Moreover, increasing evidence shows that the TMD of GPCRs affects the function of dimers. Therefore, the TMD of GPCRs is an emerging target for the development of drugs to treat diseases that involve GPCR dimerization.

Recent advances in dopamine D2 receptor ligands in the treatment of neuropsychiatric disorders

Dopamine is a biologically active amine synthesized in the central and peripheral nervous system. This biogenic monoamine acts by activating five types of dopamine receptors (D_1-5 Rs), which belong to the G protein-coupled receptor family. Antagonists and partial agonists of D₂ Rs are used to treat schizophrenia, Parkinson's disease, depression, and anxiety. The typical pharmacophore with high D₂ R affinity comprises four main areas, namely aromatic moiety, cyclic amine, central linker and aromatic/heteroaromatic lipophilic fragment. From the literature reviewed herein, we can conclude that 4-(2,3-dichlorophenyl), 4-(2-methoxyphenyl)-, 4-(benzo[b]thiophen-4-yl)-1-substituted piperazine, and 4-(6-fluorobenzo[d]isoxazol-3-yl)piperidine moieties are critical for high D₂ R affinity. Four to six atoms chains are optimal for D₂ R affinity with 4-butoxyl as the most pronounced one. The bicyclic aromatic/heteroaromatic systems are most frequently occurring as lipophilic appendages to retain high D₂ R affinity. In this review, we provide a thorough overview of the therapeutic potential of D₂ R modulators in the treatment of the aforementioned disorders. In addition, this review summarizes current knowledge about these diseases, with a focus on the dopaminergic pathway underlying these pathologies. Major attention is paid to the structure, function, and pharmacology of novel D₂ R ligands, which have been developed in the last decade (2010-2021), and belong to the 1,4-disubstituted aromatic cyclic amine group. Due to the abundance of data, allosteric D₂ R ligands and D₂ R modulators from patents are not discussed in this review.

Biased Activation Mechanism Induced by GPCR Heterodimerization: Observations from μOR/δOR Dimers

GPCRs regulate multiple intracellular signaling cascades. Biasedly activating one signaling pathway over the others provides additional clinical utility to optimize GPCR-based therapies. GPCR heterodimers possess different functions from their monomeric states, including their selectivity to different transducers. However, the biased signaling mechanism induced by the heterodimerization remains unclear. Motivated by the issue, we select an important GPCR heterodimer (μOR/δOR heterodimer) as a case and use microsecond Gaussian accelerated molecular dynamics simulation coupled with potential of mean force and protein structure network (PSN) to probe mechanisms regarding the heterodimerization-induced constitutive β-arrestin activity and efficacy change of the agonist DAMGO. The results show that only the lowest energy state of the μOR/δOR heterodimer, which adopts a slightly outward shift of TM6 and an ICL2 conformation close to the receptor core, can selectively accommodate β-arrestins. PSN further reveals important roles of H8, ICL1, and ICL2 in regulating the constitutive β-arrestin-biased activity for the apo μOR/δOR heterodimer. In addition, the heterodimerization can allosterically alter the binding mode of DAMGO mainly by means of W7.35. Consequently, DAMGO transmits the structural signal mainly through TM6 and TM7 in the dimer, rather than TM3 similar to the μOR monomer, thus changing the efficacy of DAMGO from a balanced agonist to the β-arrestin-biased one. On the other side, the binding of DAMGO to the heterodimer can stabilize μOR/δOR heterodimers through a stronger interaction of TM1/TM1 and H8/H8, accordingly enhancing the interaction of μOR with δOR and the binding affinity of the dimer to the β-arrestin. The agonist DAMGO does not change main compositions of the regulation network from the dimer interface to the transducer binding pocket of the μOR protomer, but induces an increase in the structural communication of the network, which should contribute to the enhanced β-arrestin coupling. Our observations, for the first time, reveal the molecular mechanism of the biased signaling induced by the heterodimerization for GPCRs, which should be beneficial to more comprehensively understand the GPCR bias signaling.

Computationally designed GPCR quaternary structures bias signaling pathway activation

Communication across membranes controls critical cellular processes and is achieved by receptors translating extracellular signals into selective cytoplasmic responses. While receptor tertiary structures can be readily characterized, receptor associations into quaternary structures are challenging to study and their implications in signal transduction remain poorly understood. Here, we report a computational approach for predicting receptor self-associations, and designing receptor oligomers with various quaternary structures and signaling properties. Using this approach, we designed chemokine receptor CXCR4 dimers with reprogrammed binding interactions, conformations, and abilities to activate distinct intracellular signaling proteins. In agreement with our predictions, the designed CXCR4s dimerized through distinct conformations and displayed different quaternary structural changes upon activation. Consistent with the active state models, all engineered CXCR4 oligomers activated the G protein Gi, but only specific dimer structures also recruited β-arrestins. Overall, we demonstrate that quaternary structures represent an important unforeseen mechanism of receptor biased signaling and reveal the existence of a bias switch at the dimer interface of several G protein-coupled receptors including CXCR4, mu-Opioid and type-2 Vasopressin receptors that selectively control the activation of G proteins vs β-arrestin-mediated pathways. The approach should prove useful for predicting and designing receptor associations to uncover and reprogram selective cellular signaling functions.

Class A and C GPCR Dimers in Neurodegenerative Diseases

Neurodegenerative diseases affect over 30 million people worldwide with an ascending trend. Most individuals suffering from these irreversible brain damages belong to the elderly population, with onset between 50 and 60 years. Although the pathophysiology of such diseases is partially known, it remains unclear upon which point a disease turns degenerative. Moreover, current therapeutics can treat some of the symptoms but often have severe side effects and become less effective in long-term treatment. For many neurodegenerative diseases, the involvement of G proteincoupled receptors (GPCRs), which are key players of neuronal transmission and plasticity, has become clearer and holds great promise in elucidating their biological mechanism. With this review, we introduce and summarize class A and class C GPCRs, known to form heterodimers or oligomers to increase their signalling repertoire. Additionally, the examples discussed here were shown to display relevant alterations in brain signalling and had already been associated with the pathophysiology of certain neurodegenerative diseases. Lastly, we classified the heterodimers into two categories of crosstalk, positive or negative, for which there is known evidence.

Molecular evolution and depth-related adaptations of rhodopsin in the adaptive radiation of cichlid fishes in Lake Tanganyika

The visual sensory system is essential for animals to perceive their environment and is thus under strong selection. In aquatic environments, light intensity and spectrum differ primarily along a depth gradient. Rhodopsin (RH1) is the only opsin responsible for dim‐light vision in vertebrates and has been shown to evolve in response to the respective light conditions, including along a water depth gradient in fishes. In this study, we examined the diversity and sequence evolution of RH1 in virtually the entire adaptive radiation of cichlid fishes in Lake Tanganyika, focusing on adaptations to the environmental light with respect to depth. We show that Tanganyikan cichlid genomes contain a single copy of RH1. The 76 variable amino acid sites detected in RH1 across the radiation were not uniformly distributed along the protein sequence, and 31 of these variable sites show signals of positive selection. Moreover, the amino acid substitutions at 15 positively selected sites appeared to be depth‐related, including three key tuning sites that directly mediate shifts in the peak spectral sensitivity, one site involved in protein stability and 11 sites that may be functionally important on the basis of their physicochemical properties. Among the strongest candidate sites for deep‐water adaptations are two known key tuning sites (positions 292 and 299) and three newly identified variable sites (37, 104 and 290). Our study, which is the first comprehensive analysis of RH1 evolution in a massive adaptive radiation of cichlid fishes, provides novel insights into the evolution of RH1 in a freshwater environment.

Molecular evolution and convergence of the rhodopsin gene in Gymnogobius, a goby group having diverged into coastal to freshwater habitats

Adaptive evolution of vision-related genes has been frequently observed in the process of invasion of new environments in a wide range of animal taxa. The typical example is that of the molecular evolution of rhodopsin associated with habitat changes in aquatic animals. However, few studies have investigated rhodopsin evolution during adaptive radiation across various habitats. In the present study, we examined the link between molecular evolutionary patterns in the rhodopsin gene and macroscopic habitat changes in Gymnogobius species (Gobiidae), which have adaptively radiated to diverse aquatic habitats including the sea, brackish waters, rivers and lakes. Analysis of amino acid substitutions in rhodopsin in the phylogenetic framework revealed convergent substitutions in 4-5 amino acids in three groups (four species), including two spectral tuning amino acid sites known to change rhodopsin's absorption wavelength. Positive selection was detected in the basal branches of each of these three groups, suggesting adaptive molecular convergence of rhodopsin. However, no significant correlation was observed between amino acid substitutions and the species' habitat changes, suggesting molecular adaptation to some unidentified micro-ecological environments. Taken together, these results emphasize the importance of considering not only macroscopic habitats but also micro-ecological environments when elucidating the driving forces of adaptive evolution of the visual system.

BRET Analysis of GPCR Dimers in Neurons and Non-Neuronal Cells: Evidence for Inactive, Agonist, and Constitutive Conformations

G-protein-coupled receptors (GPCRs) are dimeric proteins, but the functional consequences of the process are still debated. Active GPCR conformations are promoted either by agonists or constitutive activity. Inverse agonists decrease constitutive activity by promoting inactive conformations. The histamine H₃ receptor (H₃R) is the target of choice for the study of GPCRs because it displays high constitutive activity. Here, we study the dimerization of recombinant and brain H₃R and explore the effects of H₃R ligands of different intrinsic efficacy on dimerization. Co-immunoprecipitations and Western blots showed that H₃R dimers co-exist with monomers in transfected HEK 293 cells and in rodent brains. Bioluminescence energy transfer (BRET) analysis confirmed the existence of spontaneous H₃R dimers, not only in living HEK 293 cells but also in transfected cortical neurons. In both cells, agonists and constitutive activity of the H₃R decreased BRET signals, whereas inverse agonists and GTPγS, which promote inactive conformations, increased BRET signals. These findings show the existence of spontaneous H₃R dimers not only in heterologous systems but also in native tissues, which are able to adopt a number of allosteric conformations, from more inactive to more active states.

Ligand-Dependent Conformational Transitions in Molecular Dynamics Trajectories of GPCRs Revealed by a New Machine Learning Rare Event Detection Protocol

Central among the tools and approaches used for ligand discovery and design are Molecular Dynamics (MD) simulations, which follow the dynamic changes in molecular structure in response to the environmental condition, interactions with other proteins, and the effects of ligand binding. The need for, and successes of, MD simulations in providing this type of essential information are well documented, but so are the challenges presented by the size of the resulting datasets encoding the desired information. The difficulty of extracting information on mechanistically important state-to-state transitions in response to ligand binding and other interactions is compounded by these being rare events in the MD trajectories of complex molecular machines, such as G-protein-coupled receptors (GPCRs). To address this problem, we have developed a protocol for the efficient detection of such events. We show that the novel Rare Event Detection (RED) protocol reveals functionally relevant and pharmacologically discriminating responses to the binding of different ligands to the 5-HT_2AR orthosteric site in terms of clearly defined, structurally coherent, and temporally ordered conformational transitions. This information from the RED protocol offers new insights into specific ligand-determined functional mechanisms encoded in the MD trajectories, which opens a new and rigorously reproducible path to understanding drug activity with application in drug discovery.

Discovery of a true bivalent dopamine D2 receptor agonist

Employing two different alkyne-modified dopamine agonists to construct bivalent compounds via click chemistry resulted in the identification of a bivalent ligand (11c) for dopamine D₂ receptor homodimer, which, compared to its parent monomeric alkyne, showed a 16-fold higher binding affinity for the dopamine D₂ receptor and a 5-fold higher potency in a cAMP assay in HEK 293T cells stably expressing D₂R. Molecular modeling revealed that 11c can indeed bridge the orthosteric binding sites of a D₂R homodimer in a relaxed conformation via the TM5-TM6 interface and allows to largely rationalize the results of the receptor assays.

PAR4 activation involves extracellular loop 3 and transmembrane residue Thr153

Protease-activated receptor 4 (PAR4) mediates sustained thrombin signaling in platelets and is required for a stable thrombus. PAR4 is activated by proteolysis of the N terminus to expose a tethered ligand. The structural basis for PAR4 activation and the location of its ligand binding site (LBS) are unknown. Using hydrogen/deuterium exchange (H/D exchange), computational modeling, and signaling studies, we determined the molecular mechanism for tethered ligand-mediated PAR4 activation. H/D exchange identified that the LBS is composed of transmembrane 3 (TM3) domain and TM7. Unbiased computational modeling further predicted an interaction between Gly48 from the tethered ligand and Thr153 from the LBS. Mutating Thr153 significantly decreased PAR4 signaling. H/D exchange and modeling also showed that extracellular loop 3 (ECL3) serves as a gatekeeper for the interaction between the tethered ligand and LBS. A naturally occurring sequence variant (P310L, rs2227376) and 2 experimental mutations (S311A and P312L) determined that the rigidity conferred by prolines in ECL3 are essential for PAR4 activation. Finally, we examined the role of the polymorphism at position 310 in venous thromboembolism (VTE) using the International Network Against Venous Thrombosis (INVENT) consortium multi-ancestry genome-wide association study (GWAS) meta-analysis. Individuals with the PAR4 Leu310 allele had a 15% reduction in relative risk for VTE (odds ratio, 0.85; 95% confidence interval, 0.77-0.94) compared with the Pro310 allele. These data are consistent with our H/D exchange, molecular modeling, and signaling studies. In conclusion, we have uncovered the structural basis for PAR4 activation and identified a previously unrecognized role for PAR4 in VTE.

Structural Complexity and Plasticity of Signaling Regulation at the Melanocortin-4 Receptor

The melanocortin-4 receptor (MC4R) is a class A G protein-coupled receptor (GPCR), essential for regulation of appetite and metabolism. Pathogenic inactivating MC4R mutations are the most frequent cause of monogenic obesity, a growing medical and socioeconomic problem worldwide. The MC4R mediates either ligand-independent or ligand-dependent signaling. Agonists such as α-melanocyte-stimulating hormone (α-MSH) induce anorexigenic effects, in contrast to the endogenous inverse agonist agouti-related peptide (AgRP), which causes orexigenic effects by suppressing high basal signaling activity. Agonist action triggers the binding of different subtypes of G proteins and arrestins, leading to concomitant induction of diverse intracellular signaling cascades. An increasing number of experimental studies have unraveled molecular properties and mechanisms of MC4R signal transduction related to physiological and pathophysiological aspects. In addition, the MC4R crystal structure was recently determined at 2.75 Å resolution in an inactive state bound with a peptide antagonist. Underpinned by structural homology models of MC4R complexes simulating a presumably active-state conformation compared to the structure of the inactive state, we here briefly summarize the current understanding and key players involved in the MC4R switching process between different activity states. Finally, these perspectives highlight the complexity and plasticity in MC4R signaling regulation and identify gaps in our current knowledge.

Update of the GRIP web service

G protein-coupled receptors (GPCRs) can form homodimers, heterodimers, or higher-order molecular complexes (oligomers). The reports on the change of functions through the oligomerization have been accumulated. Inhibition of GPCR oligomerization without affecting the protomer's overall structure would clarify the oligomer-specific functions although inhibition experiments are costly and require accurate information about the interface location. Unfortunately, the number of experimentally determined interfaces is limited. The precise prediction of the oligomerization interfaces is, therefore, useful for inhibition experiments to examine the oligomer-specific functions, which would accelerate investigations of the GPCR signaling. However, interface prediction for GPCR oligomerization is difficult because different GPCR subtypes belonging to the same subfamily often use different structural regions as their interfaces. We previously developed a high-performance method to predict the interfaces for GPCR oligomerization, by identifying the conserved surfaces with the sequence and structure information. Then, the structural characteristic of a GPCR structure is regarded to be a thick-tube like conformation that is approximately perpendicular to the membrane plane. Our method had successfully predicted all of the interfaces available on that day. We had launched a web server for our interface prediction of GPCRs (GRIP). We have improved the previous version of GRIP server and enhanced its usability. First, we discarded the approximation of the GPCR structure as the thick-tube-like conformation. This improvement increased the number of structures for the prediction. Second, the FUGUE-based template recommendation service was introduced to facilitate the choice of an appropriate structure for the prediction. The new prediction server is available at http://grip.b.dendai.ac.jp/∼grip/.

Orexin/Hypocretin Type 2 Receptor (HCRTR2) Gene as a Candidate Gene in Sertraline-Associated Insomnia in Depressed Patients

BACKGROUND

Selective serotonin reuptake inhibitors (SSRIs) are considered as first-line drugs for treating depressive disorders. Among the adverse effects reported with sertraline is sleep disturbances; however, the etiology lying beneath is obscure. Orexin, the most recently discovered neurotransmitter, is involved in the sleep cycle. It exerts its physiological actions through orexin or hypocretin type 1 and 2 receptors (HCRTR1 and HCRTR2). Dysfunction of the orexin system contributes to various psychiatric, neurologic and neuropsychiatric disorders. Thus, our study aimed to assess the possible association of genetic variation of HCRTR2 G1246A with hypersomnia reported with sertraline in a group of major depressive disorder (MDD) patients.

PATIENTS AND METHODS

Ninety-six newly diagnosed MDD patients were enrolled in our cohort study. MDD was assessed using DSM-V criteria. Insomnia Severity Index (ISI) was used to assess insomnia at baseline (week 0) and week 4. Blood samples were collected for further genotyping of HCRTR2 G1246A (rs2653349) using polymerase chain reaction-restriction fragment length polymorphism.

RESULTS

A significant association between G1264A polymorphism of HCRTR2 and insomnia was observed. Insomnia with sertraline happens by 2.5-fold (P=0.022; odds ratio (OR)=2.5; 95% confidence interval (CI): 1.1-5.7) in patients having GG genotype. Patients with G allele experience insomnia by 2.1-fold more than A allele carriers (P=0.022; OR=2.1; 95% CI= 1.1-4.0). Subgroup analysis showed a significant association between GG genotype as well as the G allele and insomnia only in female MDD patients (P=0.011; OR=4.0; 95% CI=1.3-12.0 and P=0.033; OR=2.4; 95% CI=1.02-5.7, respectively).

CONCLUSION

In conclusion, the G1246A variant might be a predictor for insomnia in MDD patients treated with sertraline. Our findings support the idea that some variants of the HCRTR might contribute to inter-individual variability in the sleep pattern of patients receiving antidepressants.

The effect of dimerization on the activation and conformational dynamics of adenosine A1 receptor

The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.

Cryo-EM structure of the native rhodopsin dimer in nanodiscs

Imaging of rod photoreceptor outer-segment disc membranes by atomic force microscopy and cryo-electron tomography has revealed that the visual pigment rhodopsin, a prototypical class A G protein-coupled receptor (GPCR), can organize as rows of dimers. GPCR dimerization and oligomerization offer possibilities for allosteric regulation of GPCR activity, but the detailed structures and mechanism remain elusive. In this investigation, we made use of the high rhodopsin density in the native disc membranes and of a bifunctional cross-linker that preserves the native rhodopsin arrangement by covalently tethering rhodopsins via Lys residue side chains. We purified cross-linked rhodopsin dimers and reconstituted them into nanodiscs for cryo-EM analysis. We present cryo-EM structures of the cross-linked rhodopsin dimer as well as a rhodopsin dimer reconstituted into nanodiscs from purified monomers. We demonstrate the presence of a preferential 2-fold symmetrical dimerization interface mediated by transmembrane helix 1 and the cytoplasmic helix 8 of rhodopsin. We confirmed this dimer interface by double electron-electron resonance measurements of spin-labeled rhodopsin. We propose that this interface and the arrangement of two protomers is a prerequisite for the formation of the observed rows of dimers. We anticipate that the approach outlined here could be extended to other GPCRs or membrane receptors to better understand specific receptor dimerization mechanisms.

Coevolution of Residues Provides Evidence of a Functional Heterodimer of 5-HT2AR and 5-HT2CR Involving Both Intracellular and Extracellular Domains

Serotonin is a neurotransmitter that plays a role in regulating activities such as sleep, appetite, mood and substance abuse disorders; serotonin receptors 5-HT2AR and 5-HT2CR are active within pathways associated with substance abuse. It has been suggested that 5-HT2AR and 5-HT2CR may form a dimer that affects behavioral processes. Here we study the coevolution of residues in 5-HT2AR and 5-HT2CR to identify potential interactions between residues in both proteins. Coevolution studies can detect protein interactions, and since the thus uncovered interactions are subject to evolutionary pressure, they are likely functional. We assessed the significance of the 5-HT2AR/5-HT2CR interactions using randomized phylogenetic trees and found the coevolution significant (p-value = 0.01). We also discuss how co-expression of the receptors suggests the predicted interaction is functional. Finally, we analyze how several single nucleotide polymorphisms for the 5-HT2AR and 5-HT2CR genes affect their interaction. Our findings are the first to characterize the binding interface of 5-HT2AR/5-HT2CR and indicate a correlation between this interface and location of SNPs in both proteins.

Computational prediction of GPCR oligomerization

There has been a recent and prolific expansion in the number of GPCR crystal structures being solved: in both active and inactive forms and in complex with ligand, with G protein and with each other. Despite this, there is relatively little experimental information about the precise configuration of GPCR oligomers during these different biologically relevant states. While it may be possible to identify the experimental conditions necessary to crystallize a GPCR preferentially in a specific structural conformation, computational approaches afford a potentially more tractable means of describing the probability of formation of receptor dimers and higher order oligomers. Ensemble-based computational methods based on structurally determined dimers, coupled with a computational workflow that uses quantum mechanical methods to analyze the chemical nature of the molecular interactions at a GPCR dimer interface, will generate the reproducible and accurate predictions needed to predict previously unidentified GPCR dimers and to inform future advances in our ability to understand and begin to precisely manipulate GPCR oligomers in biological systems. It may also provide information needed to achieve an increase in the number of experimentally determined oligomeric GPCR structures.

GPCR interaction as a possible way for allosteric control between receptors

For more than twenty years now, GPCR dimers and larger oligomers have been the subject of intense debates. Evidence for a role of such complexes in receptor trafficking to and from the plasma membrane have been provided. However, one main issue is of course to determine whether or not such a phenomenon can be responsible for an allosteric and reciprocal control (allosteric control) of the subunits. Such a possibility would indeed add to the possible ways a cell integrates various signals targeting GPCRs. Among the large GPCR family, the class C receptors that include mGlu and GABA_B receptors, represent excellent models to examine such a possibility as they are mandatory dimers. In the present review, we will report on the observed allosteric interaction between the subunits of class C GPCRs, both mGluRs and GABA_BRs, and on the structural bases of these interactions. We will then discuss these findings for other GPCR types such as the rhodopsin-like class A receptors. We will show that many of the observations made with class C receptors have also been reported with class A receptors, suggesting that the mechanisms involved in the allosteric control between subunits in GPCR dimers may not be unique to class C GPCRs.

Signal Transduction and Pathogenic Modifications at the Melanocortin-4 Receptor: A Structural Perspective

The melanocortin-4 receptor (MC4R) can be endogenously activated by binding of melanocyte-stimulating hormones (MSH), which mediates anorexigenic effects. In contrast, the agouti-related peptide (AgRP) acts as an endogenous inverse agonist and suppresses ligand-independent basal signaling activity (orexigenic effects). Binding of ligands to MC4R leads to the activation of different G-protein subtypes or arrestin and concomitant signaling pathways. This receptor is a key protein in the hypothalamic regulation of food intake and energy expenditure and naturally-occurring inactivating MC4R variants are the most frequent cause of monogenic obesity. In general, obesity is a growing problem on a global scale and is of social, medical, and economic relevance. A significant goal is to develop optimized pharmacological tools targeting MC4R without adverse effects. To date, this has not been achieved because of inter alia non-selective ligands across the five functionally different MCR subtypes (MC1-5R). This motivates further investigation of (i) the three-dimensional MC4R structure, (ii) binding mechanisms of various ligands, and (iii) the molecular transfer process of signal transduction, with the aim of understanding how structural features are linked with functional-physiological aspects. Unfortunately, experimentally elucidated structural information is not yet available for the MC receptors, a group of class A G-protein coupled receptors (GPCRs). We, therefore, generated MC4R homology models and complexes with interacting partners to describe approximate structural properties associated with signaling mechanisms. In addition, molecular insights from pathogenic mutations were incorporated to discriminate more precisely their individual malfunction of the signal transfer mechanism.

Synthesis toward Bivalent Ligands for the Dopamine D2 and Metabotropic Glutamate 5 Receptors

In this study, we designed and synthesized heterobivalent ligands targeting heteromers consisting of the metabotropic glutamate 5 receptor (mGluR5) and the dopamine D₂ receptor (D₂R). Bivalent ligand 22a with a linker consisting of 20 atoms showed 4-fold increase in affinity for cells coexpressing D₂R and mGluR5 compared to cells solely expressing D₂R. Likewise, the affinity of 22a for mGluR5 increased 2-fold in the coexpressing cells. Additionally, 22a exhibited a 5-fold higher mGluR5 affinity than its monovalent precursor 21a in cells coexpressing D₂R and mGluR5. These results indicate that 22a is able to bridge binding sites on both receptors constituting the heterodimer. Likewise, cAMP assays revealed that 22a had a 4-fold higher potency in stable D₂R and mGluR5 coexpressing cell lines than 1. Furthermore, molecular modeling reveals that 22a is able to simultaneously bind both receptors by passing between the TM5-TM6 interface and establishing six protein-ligand H-bonds.

Mapping the Interface of a GPCR Dimer: A Structural Model of the A2A Adenosine and D2 Dopamine Receptor Heteromer

The A_2A adenosine (A_2AR) and D₂ dopamine (D₂R) receptors form oligomers in the cell membrane and allosteric interactions across the A_2AR–D₂R heteromer represent a target for development of drugs against central nervous system disorders. However, understanding of the molecular determinants of A_2AR–D₂R heteromerization and the allosteric antagonistic interactions between the receptor protomers is still limited. In this work, a structural model of the A_2AR–D₂R heterodimer was generated using a combined experimental and computational approach. Regions involved in the heteromer interface were modeled based on the effects of peptides derived from the transmembrane (TM) helices on A_2AR–D₂R receptor–receptor interactions in bioluminescence resonance energy transfer (BRET) and proximity ligation assays. Peptides corresponding to TM-IV and TM-V of the A_2AR blocked heterodimer interactions and disrupted the allosteric effect of A_2AR activation on D₂R agonist binding. Protein–protein docking was used to construct a model of the A_2AR–D₂R heterodimer with a TM-IV/V interface, which was refined using molecular dynamics simulations. Mutations in the predicted interface reduced A_2AR–D₂R interactions in BRET experiments and altered the allosteric modulation. The heterodimer model provided insights into the structural basis of allosteric modulation and the technique developed to characterize the A_2AR–D₂R interface can be extended to study the many other G protein-coupled receptors that engage in heteroreceptor complexes.

The Heptahelical Domain of the Sweet Taste Receptor T1R2 Is a New Allosteric Binding Site for the Sweet Taste Modulator Amiloride That Modulates Sweet Taste in a Species-Dependent Manner

The activity of sweet taste receptor (heterodimeric T1R2 and T1R3) can be modulated by sweet regulators. The compound amiloride can inhibit the sweet sensitivity of the human sweet taste receptor. This study describes the species-dependent regulation of the response of sweet taste receptors by this sweet inhibitor. Amiloride inhibited the sweet taste response of humans and mice but not that of squirrel monkeys. Using human/squirrel monkey/mouse chimeric T1R2 and T1R3 receptors as well as the agonist perillartine (which can activate the single heptahelical domain of T1R2), we found that the heptahelical domain of T1R2 is the molecular determinant that mediates the species-dependent sensitivity to this sweet regulator. Compared to the sweet inhibitor lactisole (which acts on T1R3), amiloride has a different allosteric binding site on the sweet receptor, which is important new information for the design of novel sweet taste modulators that act on T1R2.

Dynamic tuneable G protein-coupled receptor monomer-dimer populations

G protein-coupled receptors (GPCRs) are the largest class of membrane receptors, playing a key role in the regulation of processes as varied as neurotransmission and immune response. Evidence for GPCR oligomerisation has been accumulating that challenges the idea that GPCRs function solely as monomeric receptors; however, GPCR oligomerisation remains controversial primarily due to the difficulties in comparing evidence from very different types of structural and dynamic data. Using a combination of single-molecule and ensemble FRET, double electron–electron resonance spectroscopy, and simulations, we show that dimerisation of the GPCR neurotensin receptor 1 is regulated by receptor density and is dynamically tuneable over the physiological range. We propose a “rolling dimer” interface model in which multiple dimer conformations co-exist and interconvert. These findings unite previous seemingly conflicting observations, provide a compelling mechanism for regulating receptor signalling, and act as a guide for future physiological studies.

Evidence suggests oligomerisation of G protein-coupled receptors in membranes, but this is controversial. Here, authors use single-molecule and ensemble FRET, and spectroscopy to show that the neurotensin receptor 1 forms multiple dimer conformations that interconvert - “rolling” interfaces.

Structural Investigation of the Dopamine-2 Receptor Agonist Bromocriptine Binding to Dimeric D2HighR and D2LowR States

The active (D2^HighR) and inactive (D2^LowR) states of dimeric dopamine D2 receptor (D2R) models were investigated to clarify the binding mechanisms of the dopamine agonist bromocriptine, using Molecular Dynamics (MD) simulation. The aim of this comprehensive study was to investigate the critical effects of bromocriptine binding on each distinct receptor conformation. The different binding modes of the bromocriptine ligand in the active and inactive states have a significant effect on the conformational changes of the receptor. Based on the MM/GBSA approach, the calculated binding enthalpies of bromocriptine demonstrated selectivity toward the D2^HighR active state. There is good agreement between the calculated and experimentally measured D2^HighR selectivity. In the ligand-binding site, the key amino acids identified for D2^HighR were Asp114(3.32) and Glu95(2.65), and for D2^LowR, it was Ser193(5.42). Moreover, analysis of replicate MD trajectories demonstrated that the bromocriptine structure was more rigid at the D2^HighR state and more flexible at the D2^LowR state. However, the side chains of the ligand-receptor complex of D2^HighR showed larger variations relative to the corresponding regions of D2^LowR. The present study is part of an ongoing research program to study D2R conformational changes during ligand activation and to evaluate the conformational state selectivity for ligand binding.

Glutamate heteroreceptor complexes in the brain

The existence of mGluR, NMDAR, AMPAR and putative KAR heteroreceptor complexes in synaptic and extrasynaptic regions of brain glutamate synapses represents a major integrative mechanism. Our aim in the current article is to analyze if the formation of the different types glutamate hetereceptor complexes involves the contribution of triplet amino acid homologies (protriplets) in a postulated receptor interface based on the triplet puzzle theory. Seven main sets (lists) of receptor pairs in databases were used containing various sets (lists) of human receptor heteromers and nonheteromers obtained from the available scientific publications including the publically available GPCR-hetnet database. Brain mGluR1-mGluR5 and mGluR2-mGluR4 isoreceptor complexes were demonstrated with a predominant extrasynaptic localization at a post- and prejunctional localization. The existence of putative mGluR4-mGluR7 heteroreceptor complexes in the basal ganglia is proposed. Metabotropic glutamate receptor subtypes also participated in the formation of a large number of heteroreceptor complexes like mGluR1-A1R, mGluR5-A2AR, mGluR5-D2R and D2R-A2AR-mGluR5, located in relation to glutamate synapses, especially in the basal ganglia. A putative mGluR1-GABAB1/2 heterocomplex may also exist. NMDAR heteroreceptor complexes were also demonstrated as a fundamental integrative mechanism in the glutamate synapse and its extrasynaptic membranes. It represented fundamental work on inter alia NMDAR-mGluR5, NMDAR-D1R and NMDAR-D2R heteroreceptor complexes involving both antagonistic and facilitatory allosteric receptor-receptor interactions. As to AMPA receptors, a heterocomplex was found for the interaction between IFNgR1 and the AMPAR mediated via the subunit GluA1 which may be of relevance for neuroinflammation. AMPAR-D2R heteroreceptor complexes were also demonstrated. Besides glutamate heteroreceptor complexes and their allosteric receptor-receptor interactions, a significant mechanism for the functional crosstalk can also be phosphorylation and/or reorganization of adapter proteins with dynamic binding to the two receptors modulating the allosteric receptor mechanism.

Evidence for functional pre-coupled complexes of receptor heteromers and adenylyl cyclase

G protein-coupled receptors (GPCRs), G proteins and adenylyl cyclase (AC) comprise one of the most studied transmembrane cell signaling pathways. However, it is unknown whether the ligand-dependent interactions between these signaling molecules are based on random collisions or the rearrangement of pre-coupled elements in a macromolecular complex. Furthermore, it remains controversial whether a GPCR homodimer coupled to a single heterotrimeric G protein constitutes a common functional unit. Using a peptide-based approach, we here report evidence for the existence of functional pre-coupled complexes of heteromers of adenosine A_2A receptor and dopamine D₂ receptor homodimers coupled to their cognate Gs and Gi proteins and to subtype 5 AC. We also demonstrate that this macromolecular complex provides the necessary frame for the canonical Gs-Gi interactions at the AC level, sustaining the ability of a Gi-coupled GPCR to counteract AC activation mediated by a Gs-coupled GPCR.

It is unclear whether GPCRs, G proteins and adenylyl cyclase (AC) associate through random collisions or defined pre-coupling mechanisms. Using a peptide-based approach, the authors show that heteromers of adenosine A_2A and dopamine D₂ receptors form pre-coupled complexes with their cognate G proteins and AC5.

Capturing the multiscale dynamics of membrane protein complexes with all-atom, mixed-resolution, and coarse-grained models

The structures and dynamics of protein complexes are often challenging to model in heterogeneous environments such as biological membranes. Herein, we meet this fundamental challenge at attainable cost with all-atom, mixed-resolution, and coarse-grained models of vital membrane proteins. We systematically simulated five complex models formed by two distinct G protein-coupled receptors (GPCRs) in the lipid-bilayer membrane on the ns-to-μs timescales. These models, which suggest the swinging motion of an intracellular loop, for the first time, provide the molecular details for the regulatory role of such a loop. For the models at different resolutions, we observed consistent structural stability but various levels of speed-ups in protein dynamics. The mixed-resolution and coarse-grained models show two and four times faster protein diffusion than the all-atom models, in addition to a 4- and 400-fold speed-up in the simulation performance. Furthermore, by elucidating the strengths and challenges of combining all-atom models with reduced resolution models, this study can serve as a guide to simulating other complex systems in heterogeneous environments efficiently.

Probing the cooperative mechanism of the μ–δ opioid receptor heterodimer by multiscale simulation

The activation-cooperativity of the μ–δ opioid receptor heterodimer was probed by multiscale simulation coupled with a protein structure network.

An engineered opsin monomer scrambles phospholipids

The G protein-coupled receptor opsin is a phospholipid scramblase that facilitates rapid transbilayer phospholipid exchange in liposomes. The mechanism by which opsin scrambles lipids is unknown. It has been proposed that lipid translocation may occur at protein-protein interfaces of opsin dimers. To test this possibility, we rationally engineered QUAD opsin by tryptophan substitution of four lipid-facing residues in transmembrane helix 4 (TM4) that is known to be important for dimerization. Atomistic molecular dynamics simulations of wild type and QUAD opsins combined with continuum modeling revealed that the tryptophan substitutions lower the energetically unfavorable residual hydrophobic mismatch between TM4 and the membrane, reducing the drive of QUAD opsin to dimerize. We purified thermostable wild type and QUAD opsins, with or without a SNAP tag for fluorescence labeling. Single molecule fluorescence measurements of purified SNAP-tagged constructs revealed that both proteins are monomers. Fluorescence-based activity assays indicated that QUAD opsin is a fully functional scramblase. However, unlike wild type opsin which dimerizes en route to insertion into phospholipid vesicles, QUAD opsin reconstitutes as a monomer. We conclude that an engineered opsin monomer can scramble phospholipids, and that the lipid-exposed face of TM4 is unlikely to contribute to transbilayer phospholipid exchange.

The Class-A GPCR Dopamine D2 Receptor Forms Transient Dimers Stabilized by Agonists: Detection by Single-Molecule Tracking

Whether class-A G-protein coupled receptors (GPCRs) exist and work as monomers or dimers has drawn extensive attention. A class-A GPCR dopamine D2 receptor (D2R) is involved in many physiological and pathological processes and diseases, indicating its critical role in proper functioning of neuronal circuits. In particular, D2R homodimers might play key roles in schizophrenia development and amphetamine-induced psychosis. Here, using single-molecule imaging, we directly tracked single D2R molecules in the plasma membrane at a physiological temperature of 37 °C, and unequivocally determined that D2R forms transient dimers with a lifetime of 68 ms in its resting state. Agonist addition prolonged the dimer lifetime by a factor of ~1.5, suggesting the possibility that transient dimers might be involved in signaling.

Spatial Intensity Distribution Analysis: Studies of G Protein-Coupled Receptor Oligomerisation

Spatial intensity distribution analysis (SpIDA) is a recently developed approach for determining quaternary structure information on fluorophore-labelled proteins of interest in situ. It can be applied to live or fixed cells and native tissue. Using confocal images, SpIDA generates fluorescence intensity histograms that are analysed by super-Poissonian distribution functions to obtain density and quantal brightness values of the fluorophore-labelled protein of interest. This allows both expression level and oligomerisation state of the protein to be determined. We describe the application of SpIDA to investigate the oligomeric state of G protein-coupled receptors (GPCRs) at steady state and following cellular challenge, and consider how SpIDA may be used to explore GPCR quaternary organisation in pathophysiology and to stratify medicines.

GPCRs may exist and function as monomers: however, abundant evidence suggests they can form dimers/oligomers.

This concept has implications for drug discovery as it may offer opportunities to modulate the effects of known pharmaceuticals or identify new drug therapies.

A variety of approaches have been applied to this issue from traditional biochemical techniques, via resonance energy transfer approaches to recently developed image analysis-based techniques such as SpIDA. This uses mathematical analysis of confocal microscopy images to generate quantal brightness and density information for a fluorophore-tagged receptor.

SpIDA can be applied to live or fixed cells and native tissue.

SpIDA has been applied to GPCRs from each of the major subfamilies to explore their oligomerisation status at steady state and their regulation by receptor density and ligand binding.

Allosteric control of an asymmetric transduction in a G protein-coupled receptor heterodimer

GPCRs play critical roles in cell communication. Although GPCRs can form heteromers, their role in signaling remains elusive. Here we used rat metabotropic glutamate (mGlu) receptors as prototypical dimers to study the functional interaction between each subunit. mGluRs can form both constitutive homo- and heterodimers. Whereas both mGlu2 and mGlu4 couple to G proteins, G protein activation is mediated by mGlu4 heptahelical domain (HD) exclusively in mGlu2-4 heterodimers. Such asymmetric transduction results from the action of both the dimeric extracellular domain, and an allosteric activation by the partially-activated non-functional mGlu2 HD. G proteins activation by mGlu2 HD occurs if either the mGlu2 HD is occupied by a positive allosteric modulator or if mGlu4 HD is inhibited by a negative modulator. These data revealed an oriented asymmetry in mGlu heterodimers that can be controlled with allosteric modulators. They provide new insight on the allosteric interaction between subunits in a GPCR dimer.

Analysis of the Glutamate Agonist LY404,039 Binding to Nonstatic Dopamine Receptor D2 Dimer Structures and Consensus Docking

Dopamine receptor D2 (D2R) plays an important role in the human central nervous system and is a focal target of antipsychotic agents. The D2^HighR and D2^LowR dimeric models previously developed by our group are used to investigate the prediction of binding affinity of the LY404,039 ligand and its binding mechanism within the catalytic domain. The computational data obtained using molecular dynamics simulations fit well with the experimental results. The calculated binding affinities of LY404,039 using MM/PBSA for the D2^HighR and D2^LowR targets were -12.04 and -9.11 kcal/mol, respectively. The experimental results suggest that LY404,039 binds to D2^HighR and D2^LowR with binding affinities (K_i) of 8.2 and 1640 nM, respectively. The high binding affinity of LY404,039 in terms of binding to [³H]domperidone was inhibited by the presence of a guanine nucleotide, indicating an agonist action of the drug at D2^HighR. The interaction analysis demonstrated that while Asp114 was among the most critical amino acids for D2^HighR binding, residues Ser193 and Ser197 were significantly more important within the binding cavity of D2^LowR. Molecular modeling analyses are extended to ensemble docking as well as structure-based pharmacophore model (E-pharmacophore) development using the bioactive conformation of LY404,039 at the binding pocket as a template and screening of small-molecule databases with derived pharmacophore models.

A Molecular Basis for Selective Antagonist Destabilization of Dopamine D3 Receptor Quaternary Organization

The dopamine D3 receptor (D3R) is a molecular target for both first-generation and several recently-developed antipsychotic agents. Following stable expression of this mEGFP-tagged receptor, Spatial Intensity Distribution Analysis indicated that a substantial proportion of the receptor was present within dimeric/oligomeric complexes and that increased expression levels of the receptor favored a greater dimer to monomer ratio. Addition of the antipsychotics, spiperone or haloperidol, resulted in re-organization of D3R quaternary structure to promote monomerization. This action was dependent on ligand concentration and reversed upon drug washout. By contrast, a number of other antagonists with high affinity at the D3R, did not alter the dimer/monomer ratio. Molecular dynamics simulations following docking of each of the ligands into a model of the D3R derived from the available atomic level structure, and comparisons to the receptor in the absence of ligand, were undertaken. They showed that, in contrast to the other antagonists, spiperone and haloperidol respectively increased the atomic distance between reference α carbon atoms of transmembrane domains IV and V and I and II, both of which provide key interfaces for D3R dimerization. These results offer a molecular explanation for the distinctive ability of spiperone and haloperidol to disrupt D3R dimerization.

Biological Insights of the Dopaminergic Stabilizer ACR16 at the Binding Pocket of Dopamine D2 Receptor

The dopamine D2 receptor (D2R) plays an important part in the human central nervous system and it is considered to be a focal target of antipsychotic agents. It is structurally modeled in active and inactive states, in which homodimerization reaction of the D2R monomers is also applied. The ASP2314 (also known as ACR16) ligand, a D2R stabilizer, is used in tests to evaluate how dimerization and conformational changes may alter the ligand binding space and to provide information on alterations in inhibitory mechanisms upon activation. The administration of the D2R agonist ligand ACR16 [³H](+)-4-propyl-3,4,4a,5,6,10b-hexahydro-2H-naphtho[1,2-b][1,4]oxazin-9-ol ((+)PHNO) revealed K_i values of 32 nM for the D2^highR and 52 μM for the D2^lowR. The calculated binding affinities of ACR16 with post processing molecular dynamics (MD) simulations analyses using MM/PBSA for the monomeric and homodimeric forms of the D2^highR were -9.46 and -8.39 kcal/mol, respectively. The data suggests that the dimerization of the D2R leads negative cooperativity for ACR16 binding. The dimerization reaction of the D2^highR is energetically favorable by -22.95 kcal/mol. The dimerization reaction structurally and thermodynamically stabilizes the D2^highR conformation, which may be due to the intermolecular forces formed between the TM4 of each monomer, and the result strongly demonstrates dimerization essential for activation of the D2R.

Structural-Functional Features of the Thyrotropin Receptor: A Class A G-Protein-Coupled Receptor at Work

The thyroid-stimulating hormone receptor (TSHR) is a member of the glycoprotein hormone receptors, a sub-group of class A G-protein-coupled receptors (GPCRs). TSHR and its endogenous ligand thyrotropin (TSH) are of essential importance for growth and function of the thyroid gland and proper function of the TSH/TSHR system is pivotal for production and release of thyroid hormones. This receptor is also important with respect to pathophysiology, such as autoimmune (including ophthalmopathy) or non-autoimmune thyroid dysfunctions and cancer development. Pharmacological interventions directly targeting the TSHR should provide benefits to disease treatment compared to currently available therapies of dysfunctions associated with the TSHR or the thyroid gland. Upon TSHR activation, the molecular events conveying conformational changes from the extra- to the intracellular side of the cell across the membrane comprise reception, conversion, and amplification of the signal. These steps are highly dependent on structural features of this receptor and its intermolecular interaction partners, e.g., TSH, antibodies, small molecules, G-proteins, or arrestin. For better understanding of signal transduction, pathogenic mechanisms such as autoantibody action and mutational modifications or for developing new pharmacological strategies, it is essential to combine available structural data with functional information to generate homology models of the entire receptor. Although so far these insights are fragmental, in the past few decades essential contributions have been made to investigate in-depth the involved determinants, such as by structure determination via X-ray crystallography. This review summarizes available knowledge (as of December 2016) concerning the TSHR protein structure, associated functional aspects, and based on these insights we suggest several receptor complex models. Moreover, distinct TSHR properties will be highlighted in comparison to other class A GPCRs to understand the molecular activation mechanisms of this receptor comprehensively. Finally, limitations of current knowledge and lack of information are discussed highlighting the need for intensified efforts toward TSHR structure elucidation.

Beta2-adrenergic receptor homodimers: Role of transmembrane domain 1 and helix 8 in dimerization and cell surface expression

Even though there are hundreds of reports in the published literature supporting the hypothesis that G protein-coupled receptors (GPCR) form and function as dimers this remains a highly controversial area of research and mechanisms governing homodimer formation are poorly understood. Crystal structures revealing homodimers have been reported for many different GPCR. For adrenergic receptors, a potential dimer interface involving transmembrane domain 1 (TMD1) and helix 8 (H8) was identified in crystal structures of the beta₁-adrenergic (β₁-AR) and β₂-AR. The purpose of this study was to investigate a potential role for TMD1 and H8 in dimerization and plasma membrane expression of functional β₂-AR. Charged residues at the base of TMD1 and in the distal portion of H8 were replaced, singly and in combination, with non-polar residues or residues of opposite charge. Wild type and mutant β₂-AR, tagged with YFP and expressed in HEK293 cells, were evaluated for plasma membrane expression and function. Homodimer formation was evaluated using bioluminescence resonance energy transfer, bimolecular fluorescence complementation, and fluorescence correlation spectroscopy. Amino acid substitutions at the base of TMD1 and in the distal portion of H8 disrupted homodimer formation and caused receptors to be retained in the endoplasmic reticulum. Mutations in the proximal region of H8 did not disrupt dimerization but did interfere with plasma membrane expression. This study provides biophysical evidence linking a potential TMD1/H8 interface with ER export and the expression of functional β₂-AR on the plasma membrane. This article is part of a Special Issue entitled: Interactions between membrane receptors in cellular membranes edited by Kalina Hristova.

Impact of Lipid Composition and Receptor Conformation on the Spatio-temporal Organization of μ-Opioid Receptors in a Multi-component Plasma Membrane Model

The lipid composition of cell membranes has increasingly been recognized as playing an important role in the function of various membrane proteins, including G Protein-Coupled Receptors (GPCRs). For instance, experimental and computational evidence has pointed to lipids influencing receptor oligomerization directly, by physically interacting with the receptor, and/or indirectly, by altering the bulk properties of the membrane. While the exact role of oligomerization in the function of class A GPCRs such as the μ-opioid receptor (MOR) is still unclear, insight as to how these receptors oligomerize and the relevance of the lipid environment to this phenomenon is crucial to our understanding of receptor function. To examine the effect of lipids and different MOR conformations on receptor oligomerization we carried out extensive coarse-grained molecular dynamics simulations of crystal structures of inactive and/or activated MOR embedded in an idealized mammalian plasma membrane composed of 63 lipid types asymmetrically distributed across the two leaflets. The results of these simulations point, for the first time, to specific direct and indirect effects of the lipids, as well as the receptor conformation, on the spatio-temporal organization of MOR in the plasma membrane. While sphingomyelin-rich, high-order lipid regions near certain transmembrane (TM) helices of MOR induce an effective long-range attractive force on individual protomers, both long-range lipid order and interface formation are found to be conformation dependent, with a larger number of different interfaces formed by inactive MOR compared to active MOR.

The μ-opioid receptor (MOR) is an important pharmaceutical target in the treatment of pain. In order to develop novel pain therapies, devoid of the serious side-effects of present opioid analgesics, we need to understand the fundamentals of how MOR works on the molecular level. While some studies suggest that oligomers of MOR could play a role in signaling, how MOR forms dimers, which interfaces form, and the exact role of oligomers in MOR function remain unclear. While research has shown that the membrane environment can affect membrane protein function, most previous computational work to study oligomerization has been performed in a very simple membrane. Here, we use molecular dynamics simulations of MOR in a heterogeneous plasma membrane model (comprising 63 lipid types) to investigate how the presence of the protein modulates its lipid environment, affecting species distribution and sculpting characteristic order and thickness profiles around the receptors. Such modulations, in turn, induce long-range interactions between the proteins and favor the formation of specific dimeric conformations.

Insights into Basal Signaling Regulation, Oligomerization, and Structural Organization of the Human G-Protein Coupled Receptor 83

The murine G-protein coupled receptor 83 (mGPR83) is expressed in the hypothalamus and was previously suggested to be involved in the regulation of metabolism. The neuropeptide PEN has been recently identified as a potent GPR83 ligand. Moreover, GPR83 constitutes functionally relevant hetero-oligomers with other G-protein coupled receptors (GPCR) such as the ghrelin receptor (GHSR) or GPR171. Previous deletion studies also revealed that the long N-terminal extracellular receptor domain (eNDo) of mGPR83 may act as an intra-molecular ligand, which participates in the regulation of basal signaling activity, which is a key feature of GPCR function. Here, we investigated particular amino acids at the eNDo of human GPR83 (hGPR83) by side-directed mutagenesis to identify determinants of the internal ligand. These studies were accompanied by structure homology modeling to combine functional insights with structural information. The capacity for hetero-oligomer formation of hGPR83 with diverse family A GPCRs such as the melanocortin-4 receptor (MC4R) was also investigated, with a specific emphasis on the impact of the eNDo on oligomerization and basal signaling properties. Finally, we demonstrate that hGPR83 exhibits an unusual basal signaling for different effectors, which also supports signaling promiscuity. hGPR83 interacts with a variety of hypothalamic GPCRs such as the MC4R or GHSR. These interactions are not dependent on the ectodomain and most likely occur at interfaces constituted in the transmembrane regions. Moreover, several amino acids at the transition between the eNDo and transmembrane helix 1 were identified, where mutations lead also to biased basal signaling modulation.

Dynamic Cholesterol-Conditioned Dimerization of the G Protein Coupled Chemokine Receptor Type 4

G protein coupled receptors (GPCRs) allow for the transmission of signals across biological membranes. For a number of GPCRs, this signaling was shown to be coupled to prior dimerization of the receptor. The chemokine receptor type 4 (CXCR4) was reported before to form dimers and their functionality was shown to depend on membrane cholesterol. Here, we address the dimerization pattern of CXCR4 in pure phospholipid bilayers and in cholesterol-rich membranes. Using ensembles of molecular dynamics simulations, we show that CXCR4 dimerizes promiscuously in phospholipid membranes. Addition of cholesterol dramatically affects the dimerization pattern: cholesterol binding largely abolishes the preferred dimer motif observed for pure phospholipid bilayers formed mainly by transmembrane helices 1 and 7 (TM1/TM5-7) at the dimer interface. In turn, the symmetric TM3,4/TM3,4 interface is enabled first by intercalating cholesterol molecules. These data provide a molecular basis for the modulation of GPCR activity by its lipid environment.

Membrane-Mediated Oligomerization of G Protein Coupled Receptors and Its Implications for GPCR Function

The dimerization or even oligomerization of G protein coupled receptors (GPCRs) causes ongoing, controversial debates about its functional role and the coupled biophysical, biochemical or biomedical implications. A continously growing number of studies hints to a relation between oligomerization and function of GPCRs and strengthens the assumption that receptor assembly plays a key role in the regulation of protein function. Additionally, progress in the structural analysis of GPCR-G protein and GPCR-ligand interactions allows to distinguish between actively functional and non-signaling complexes. Recent findings further suggest that the surrounding membrane, i.e., its lipid composition may modulate the preferred dimerization interface and as a result the abundance of distinct dimeric conformations. In this review, the association of GPCRs and the role of the membrane in oligomerization will be discussed. An overview of the different reported oligomeric interfaces is provided and their capability for signaling discussed. The currently available data is summarized with regard to the formation of GPCR oligomers, their structures and dependency on the membrane microenvironment as well as the coupling of oligomerization to receptor function.