Common activation mechanism of class A GPCRs

Distinct activation mechanisms regulate subtype selectivity of Cannabinoid receptors

Design of cannabinergic subtype selective ligands is challenging because of high sequence and structural similarities of cannabinoid receptors (CB₁ and CB₂). We hypothesize that the subtype selectivity of designed selective ligands can be explained by the ligand binding to the conformationally distinct states between cannabinoid receptors. Analysis of ~ 700 μs of unbiased simulations using Markov state models and VAMPnets identifies the similarities and distinctions between the activation mechanism of both receptors. Structural and dynamic comparisons of metastable intermediate states allow us to observe the distinction in the binding pocket volume change during CB₁ and CB₂ activation. Docking analysis reveals that only a few of the intermediate metastable states of CB₁ show high affinity towards CB₂ selective agonists. In contrast, all the CB₂ metastable states show a similar affinity for these agonists. These results mechanistically explain the subtype selectivity of these agonists by deciphering the activation mechanism of cannabinoid receptors.

G protein-coupled receptors in neurodegenerative diseases and psychiatric disorders

Neuropsychiatric disorders are multifactorial disorders with diverse aetiological factors. Identifying treatment targets is challenging because the diseases are resulting from heterogeneous biological, genetic, and environmental factors. Nevertheless, the increasing understanding of G protein-coupled receptor (GPCR) opens a new possibility in drug discovery. Harnessing our knowledge of molecular mechanisms and structural information of GPCRs will be advantageous for developing effective drugs. This review provides an overview of the role of GPCRs in various neurodegenerative and psychiatric diseases. Besides, we highlight the emerging opportunities of novel GPCR targets and address recent progress in GPCR drug development.

Computational insights into ligand-induced G protein and β-arrestin signaling of the dopamine D1 receptor

The dopamine D1 receptor (D1R), is a class A G protein coupled-receptor (GPCR) which has been a promising drug target for psychiatric and neurological disorders such as Parkinson's disease (PD). Previous studies have suggested that therapeutic effects can be realized by targeting the β-arrestin signaling pathway of dopamine receptors, while overactivation of the G protein-dependent pathways leads to side effects, such as dyskinesias. Therefore, it is highly desirable to develop a D1R ligand that selectively regulates the β-arrestin pathway. Currently, most D1R agonists are signaling-balanced and stimulate both G protein and β-arrestin pathways, with a few reports of G protein biased ligands. However, identification and characterization of β-arrestin biased D1R agonists has been a challenge thus far. In this study, we implemented Gaussian accelerated molecular dynamics (GaMD) simulations to provide valuable computational insights into the possible underlying molecular mechanism of the different signaling properties of two catechol and two non-catechol D1R agonists that are either G protein biased or signaling-balanced. Dynamic network analysis further identified critical residues in the allosteric signaling network of D1R for each ligand at different conformational or binding states. Some of these residues are crucial for G protein or arrestin signals of GPCRs based on previous studies. Finally, we provided a molecular design strategy which can be utilized by medicinal chemists to develop potential β-arrestin biased D1R ligands. The proposed hypotheses are experimentally testable and can guide the development of safer and more effective medications for a variety of CNS disorders.

Structural insights into constitutive activity of 5-HT6 receptor

While most therapeutic research on G-protein-coupled receptors (GPCRs) focuses on receptor activation by (endogenous) agonists, significant therapeutic potential exists through agonist-independent intrinsic constitutive activity that can occur in various physiological and pathophysiological settings. For example, inhibiting the constitutive activity of 5-HT6R-a receptor that is found almost exclusively in the brain and mediates excitatory neurotransmission-has demonstrated a therapeutic effect on cognitive/memory impairment associated with several neuropsychiatric disorders. However, the structural basis of such constitutive activity remains unclear. Here, we present a cryo-EM structure of serotonin-bound human 5-HT6R-Gs heterotrimer at 3.0-Å resolution. Detailed analyses of the structure complemented by comprehensive interrogation of signaling illuminate key structural determinants essential for constitutive 5-HT6R activity. Additional structure-guided mutagenesis leads to a nanobody mimic Gαs for 5-HT6R that can reduce its constitutive activity. Given the importance of 5-HT6R for a large number of neuropsychiatric disorders, insights derived from these studies will accelerate the design of more effective medications, and shed light on the molecular basis of constitutive activity.

Activation mechanism of the human Smoothened receptor

Smoothened (SMO) is a membrane protein of the class F subfamily of G protein-coupled receptors (GPCRs) and maintains homeostasis of cellular differentiation. SMO undergoes conformational change during activation, transmitting the signal across the membrane, making it amenable to bind to its intracellular signaling partner. Receptor activation has been studied at length for class A receptors, but the mechanism of class F receptor activation remains unknown. Agonists and antagonists bound to SMO at sites in the transmembrane domain (TMD) and the cysteine-rich domain have been characterized, giving a static view of the various conformations SMO adopts. Although the structures of the inactive and active SMO outline the residue-level transitions, a kinetic view of the overall activation process remains unexplored for class F receptors. We describe SMO's activation process in atomistic detail by performing 300 μs of molecular dynamics simulations and combining it with Markov state model theory. A molecular switch, conserved across class F and analogous to the activation-mediating D-R-Y motif in class A receptors, is observed to break during activation. We also show that this transition occurs in a stage-wise movement of the transmembrane helices: TM6 first, followed by TM5. To see how modulators affect SMO activity, we simulated agonist and antagonist-bound SMO. We observed that agonist-bound SMO has an expanded hydrophobic tunnel in SMO's core TMD, whereas antagonist-bound SMO shrinks this tunnel, further supporting the hypothesis that cholesterol travels through a tunnel inside Smoothened to activate it. In summary, this study elucidates the distinct activation mechanism of class F GPCRs and shows that SMO's activation process rearranges the core TMD to open a hydrophobic conduit for cholesterol transport.

Neurotensin Receptor Allosterism Revealed in Complex with a Biased Allosteric Modulator

The NTSR1 neurotensin receptor (NTSR1) is a G protein-coupled receptor (GPCR) found in the brain and peripheral tissues with neurotensin (NTS) being its endogenous peptide ligand. In the brain, NTS modulates dopamine neuronal activity, induces opioid-independent analgesia, and regulates food intake. Recent studies indicate that biasing NTSR1 toward β-arrestin signaling can attenuate the actions of psychostimulants and other drugs of abuse. Here, we provide the cryoEM structures of NTSR1 ternary complexes with heterotrimeric Gq and GoA with and without the brain-penetrant small-molecule SBI-553. In functional studies, we discovered that SBI-553 displays complex allosteric actions exemplified by negative allosteric modulation for G proteins that are Gα subunit selective and positive allosteric modulation and agonism for β-arrestin translocation at NTSR1. Detailed structural analysis of the allosteric binding site illuminated the structural determinants for biased allosteric modulation of SBI-553 on NTSR1.

Prediction of dynamic allostery for the transmembrane domain of the sweet taste receptor subunit, TAS1R3

The sweet taste receptor plays an essential role as an energy sensor by detecting carbohydrates. However, the dynamic mechanisms of receptor activation remain unclear. Here, we describe the interactions between the transmembrane domain of the G protein-coupled sweet receptor subunit, TAS1R3, and allosteric modulators. Molecular dynamics simulations reproduced species-specific sensitivity to ligands. We found that a human-specific sweetener, cyclamate, interacted with the mouse receptor as a negative allosteric modulator. Agonist-induced allostery during receptor activation was found to destabilize the intracellular part of the receptor, which potentially interfaces with the Gα subunit, through ionic lock opening. A common human variant (R757C) of the TAS1R3 exhibited a reduced response to sweet taste, in support of our predictions. Furthermore, histidine residues in the binding site acted as pH-sensitive microswitches to modulate the sensitivity to saccharin. This study provides important insights that may facilitate the prediction of dynamic activation mechanisms for other G protein-coupled receptors.

Structural insights into the human niacin receptor HCA2-Gi signalling complex

The hydroxycarboxylic acid receptor 2 (HCA2) agonist niacin has been used as treatment for dyslipidemia for several decades albeit with skin flushing as a common side-effect in treated individuals. Extensive efforts have been made to identify HCA2 targeting lipid lowering agents with fewer adverse effects, despite little being known about the molecular basis of HCA2 mediated signalling. Here, we report the cryo-electron microscopy structure of the HCA2-G_i signalling complex with the potent agonist MK-6892, along with crystal structures of HCA2 in inactive state. These structures, together with comprehensive pharmacological analysis, reveal the ligand binding mode and activation and signalling mechanisms of HCA2. This study elucidates the structural determinants essential for HCA2 mediated signalling and provides insights into ligand discovery for HCA2 and related receptors.

Molecular simulations of SSTR2 dynamics and interaction with ligands

The cyclic peptide hormone somatostatin regulates physiological processes involved in growth and metabolism, through its binding to G-protein coupled somatostatin receptors. The isoform 2 (SSTR2) is of particular relevance for the therapy of neuroendocrine tumours for which different analogues to somatostatin are currently in clinical use. We present an extensive and systematic computational study on the dynamics of SSTR2 in three different states: active agonist-bound, inactive antagonist-bound and apo inactive. We exploited the recent burst of SSTR2 experimental structures to perform μs-long multi-copy molecular dynamics simulations to sample conformational changes of the receptor and rationalize its binding to different ligands (the agonists somatostatin and octreotide, and the antagonist CYN154806). Our findings suggest that the apo form is more flexible compared to the holo ones, and confirm that the extracellular loop 2 closes upon the agonist octreotide but not upon the antagonist CYN154806. Based on interaction fingerprint analyses and free energy calculations, we found that all peptides similarly interact with residues buried into the binding pocket. Conversely, specific patterns of interactions are found with residues located in the external portion of the pocket, at the basis of the extracellular loops, particularly distinguishing the agonists from the antagonist. This study will help in the design of new somatostatin-based compounds for theranostics of neuroendocrine tumours.

Cryo-EM structure of the human chemerin receptor 1-Gi protein complex bound to the C-terminal nonapeptide of chemerin

Chemerin is a processed protein that acts on G protein-coupled receptors (GPCRs) for its chemotactic and adipokine activities. The biologically active chemerin (chemerin 21-157) results from proteolytic cleavage of prochemerin and uses its C-terminal peptide containing the sequence YFPGQFAFS for receptor activation. Here we report a high-resolution cryo-electron microscopy (cryo-EM) structure of human chemerin receptor 1 (CMKLR1) bound to the C-terminal nonapeptide of chemokine (C9) in complex with Gi proteins. C9 inserts its C terminus into the binding pocket and is stabilized through hydrophobic interactions involving its Y1, F2, F6, and F8, as well as polar interactions between G4, S9, and several amino acids lining the binding pocket of CMKLR1. Microsecond scale molecular dynamics simulations support a balanced force distribution across the whole ligand-receptor interface that enhances thermodynamic stability of the captured binding pose of C9. The C9 interaction with CMKLR1 is drastically different from chemokine recognition by chemokine receptors, which follow a two-site two-step model. In contrast, C9 takes an "S"-shaped pose in the binding pocket of CMKLR1 much like angiotensin II in the AT1 receptor. Our mutagenesis and functional analyses confirmed the cryo-EM structure and key residues in the binding pocket for these interactions. Our findings provide a structural basis for chemerin recognition by CMKLR1 for the established chemotactic and adipokine activities.

Photodynamic Activation of Cholecystokinin 1 Receptor Is Conserved in Mammalian and Avian Pancreatic Acini

Cholecystokinin 1 receptor (CCK1R) is the only G protein coupled receptor that is activated in type II photodynamic action, but whether this is a property common to both mammalian and avian species is not known. In this work, pancreatic acini were isolated from the rat, mouse, and Peking duck, and photodynamic CCK1R activation was examined. Isolated pancreatic acini were exposed to photosensitizer sulphonated aluminum phthalocyanine (SALPC) and photodynamic action elicited by a brief light-emitting diode (LED 675 nm) pulse (1.5 min); photodynamic CCK1R activation was assessed by Fura-2 fluorescent calcium imaging. Photodynamic action was found to induce persistent calcium oscillations in rat, mouse, and Peking duck pancreatic acini, with the sensitivity order of mouse > rat > Peking duck. Photodynamically-activated CCK1R could be inhibited reversibly by CCK1R antagonist devazepide (1 μM); photodynamic CCK1R activation was blocked by pre-incubation with 1O2 quencher Trolox C (300 µM). The sensitivity of photodynamic CCK1R activation was correlated with the increasing size of the disordered region in intracellular loop 3. These data suggest that photodynamic CCK1R activation is conserved in both mammalian and avian species, as evidenced by the presence of the photodynamic activation motif “YFM” in transmembrane domain 3.

Photodynamic Activation of Cholecystokinin 1 Receptor Is Conserved in Mammalian and Avian Pancreatic Acini

Cholecystokinin 1 receptor (CCK1R) is the only G protein coupled receptor that is activated in type II photodynamic action, but whether this is a property common to both mammalian and avian species is not known. In this work, pancreatic acini were isolated from the rat, mouse, and Peking duck, and photodynamic CCK1R activation was examined. Isolated pancreatic acini were exposed to photosensitizer sulphonated aluminum phthalocyanine (SALPC) and photodynamic action elicited by a brief light-emitting diode (LED 675 nm) pulse (1.5 min); photodynamic CCK1R activation was assessed by Fura-2 fluorescent calcium imaging. Photodynamic action was found to induce persistent calcium oscillations in rat, mouse, and Peking duck pancreatic acini, with the sensitivity order of mouse > rat > Peking duck. Photodynamically-activated CCK1R could be inhibited reversibly by CCK1R antagonist devazepide (1 μM); photodynamic CCK1R activation was blocked by pre-incubation with 1O2 quencher Trolox C (300 µM). The sensitivity of photodynamic CCK1R activation was correlated with the increasing size of the disordered region in intracellular loop 3. These data suggest that photodynamic CCK1R activation is conserved in both mammalian and avian species, as evidenced by the presence of the photodynamic activation motif “YFM” in transmembrane domain 3.

Structural basis of peptide recognition and activation of endothelin receptors

Endothelin system comprises three endogenous 21-amino-acid peptide ligands endothelin-1, -2, and -3 (ET-1/2/3), and two G protein-coupled receptor (GPCR) subtypes-endothelin receptor A (ET_AR) and B (ET_BR). Since ET-1, the first endothelin, was identified in 1988 as one of the most potent endothelial cell-derived vasoconstrictor peptides with long-lasting actions, the endothelin system has attracted extensive attention due to its critical role in vasoregulation and close relevance in cardiovascular-related diseases. Here we present three cryo-electron microscopy structures of ET_AR and ET_BR bound to ET-1 and ET_BR bound to the selective peptide IRL1620. These structures reveal a highly conserved recognition mode of ET-1 and characterize the ligand selectivity by ETRs. They also present several conformation features of the active ETRs, thus revealing a specific activation mechanism. Together, these findings deepen our understanding of endothelin system regulation and offer an opportunity to design selective drugs targeting specific ETR subtypes.

Allosteric modulation of a human odorant receptor

Odor perception is first determined by how the myriad of environmental volatiles are detected at the periphery of the olfactory system. The combinatorial activation of dedicated odorant receptors generates enough encoding power for the discrimination of tens of thousands of odorants. Recent studies have revealed that odorant receptors undergo widespread inhibitory modulation of their activity when presented with mixtures of odorants, a property likely required to maintain discrimination and ensure sparsity of the code for complex mixtures. Here, we establish the role of human OR5AN1 in the detection of musks and identify distinct odorants capable of enhancing its activity in binary mixtures. Chemical and pharmacological characterization indicate that specific α-β unsaturated aliphatic aldehydes act as positive allosteric modulators. Sensory experiments show decreased odor detection threshold in humans, suggesting that allosteric modulation of odorant receptors is perceptually relevant and likely adds another layer of complexity to how odors are encoded in the peripheral olfactory system.

Alleviation of Cocaine Withdrawal and Pertinent Interactions between Salvinorin-Based Antagonists and Kappa Opioid Receptor

The kappa opioid receptor (KOR) is involved in the regulation of both the reward and mood processes. Recent reports find that the use of drugs of abuse increases the production of dynorphin and the overall activation of KOR. Long-acting KOR antagonists, such as norbinaltorphimine (nor-BNI), JDTic, and 5'-guanidinonaltrindole (GNTI), have been shown to stop depressive and anxiety-related disorders, which are the common side effects of withdrawal that can lead to a relapse in drug use. Unfortunately, these prototypical KOR antagonists are known to induce selective KOR antagonism that is delayed by hours and extremely prolonged, and their use in humans comes with serious safety concerns because they possess a large window for potential drug-drug interactions. Furthermore, their persistent pharmacodynamic activities can hinder the ability to reverse unanticipated side effects immediately. Herein, we report our studies of the lead selective, salvinorin-based KOR antagonist (1) as well as nor-BNI on C57BL/6N male mice for spontaneous cocaine withdrawal. Assessment of pharmacokinetics showed that 1 is a short-acting compound with an average half-life of 3.75 h across different compartments (brain, spinal cord, liver, and plasma). Both 1 (5 mg/kg) and nor-BNI (5 mg/kg) were shown to reduce spontaneous withdrawal behavior in mice, with 1 producing additional anti-anxiety-like behavior in a light-dark transition test (however, no mood-related effects of 1 or nor-BNI were observed at the current dosing in an elevated plus maze or a tail suspension test). Our results support the study of selective, short-acting KOR antagonists for the treatment of psychostimulant withdrawal and the associated negative mood states that contribute to relapse. Furthermore, we identified pertinent interactions between 1 and KOR via computational studies, including induced-fit docking, mutagenesis, and molecular dynamics simulations, to gain insight into the design of future selective, potent, and short-acting salvinorin-based KOR antagonists.

FSHR activation through small molecule modulators: Mechanistic insights from MD simulations

Follicle-stimulating hormone receptor (FSHR) is a glycoprotein hormone receptor that plays a vital role in reproduction, cancer progression and osteoporosis. Owing to its therapeutic importance, several small molecule modulators have been identified by researchers through high throughput studies that usually include virtual screening of chemical libraries followed by in vitro validation through radio-ligand binding assays, cAMP accumulation and luciferase-based luminescence assays. The binding site of these modulators and structural changes that accompany modulator binding remains elusive. Here, we address these aspects through molecular docking and MD simulations on well-studied FSHR modulators and comparing the domain motions between agonist/FSH bound and antagonist bound FSHR structures. It was observed that agonist and antagonist modulators bind to the same site, but interact with distinct residues in transmembrane domain(TMD). FSHR(TMD) residues Ile522, Ala595, Ile602 and Val604 were found to interact only with agonist. Notably, these residues are conserved in the close homolog luteinizing hormone/choriogonadotropin receptor (LHCGR) and participate in interaction with its agonist Org43553. We observed distinctly prominent domain motions and conformational changes in TM helices 3, 4 and 6 for agonist bound FSHR structure. These structural changes have also been reported for LHCGR, and few GPCR members suggesting an important and well conserved mechanism of GPHR activation that could be exploited for design of novel modulators.

Structural basis of odorant recognition by a human odorant receptor

Our sense of smell enables us to navigate a vast space of chemically diverse odour molecules. This task is accomplished by the combinatorial activation of approximately 400 odorant G protein-coupled receptors encoded in the human genome1-3. How odorants are recognized by odorant receptors remains unclear. Here we provide mechanistic insight into how an odorant binds to a human odorant receptor. Using cryo-electron microscopy, we determined the structure of the active human odorant receptor OR51E2 bound to the fatty acid propionate. Propionate is bound within an occluded pocket in OR51E2 and makes specific contacts critical to receptor activation. Mutation of the odorant-binding pocket in OR51E2 alters the recognition spectrum for fatty acids of varying chain length, suggesting that odorant selectivity is controlled by tight packing interactions between an odorant and an odorant receptor. Molecular dynamics simulations demonstrate that propionate-induced conformational changes in extracellular loop 3 activate OR51E2. Together, our studies provide a high-resolution view of chemical recognition of an odorant by a vertebrate odorant receptor, providing insight into how this large family of G protein-coupled receptors enables our olfactory sense.

Structural Basis of the Binding Mode of the Antineoplastic Compound Motixafortide (BL-8040) in the CXCR4 Chemokine Receptor

Modulation of the CXCL12-CXCR4 signaling axis is of the utmost importance due to its central involvement in several pathological disorders, including inflammatory diseases and cancer. Among the different currently available drugs that inhibit CXCR4 activation, motixafortide-a best-in-class antagonist of this GPCR receptor-has exhibited promising results in preclinical studies of pancreatic, breast, and lung cancers. However, detailed information on the interaction mechanism of motixafortide is still lacking. Here, we characterize the motixafortide/CXCR4 and CXCL12/CXCR4 protein complexes by using computational techniques including unbiased all-atom molecular dynamics simulations. Our microsecond-long simulations of the protein systems indicate that the agonist triggers changes associated with active-like GPCR conformations, while the antagonist favors inactive conformations of CXCR4. Detailed ligand-protein analysis indicates the importance of motixafortide's six cationic residues, all of which established charge-charge interactions with acidic CXCR4 residues. Furthermore, two synthetic bulky chemical moieties of motixafortide work in tandem to restrict the conformations of important residues associated with CXCR4 activation. Our results not only elucidate the molecular mechanism by which motixafortide interacts with the CXCR4 receptor and stabilizes its inactive states, but also provide essential information to rationally design CXCR4 inhibitors that preserve the outstanding pharmacological features of motixafortide.

Application of computational methods for class A GPCR Ligand discovery

G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for drug development due to their role in transmitting cellular signals in a multitude of biological processes. Of the six classes categorizing GPCR (A, B, C, D, E, and F), class A contains the largest number of therapeutically relevant GPCR. Despite their importance as drug targets, many challenges exist for the discovery of novel class A GPCR ligands serving as drug precursors. Though knowledge of the structural and functional characteristics of GPCR has grown significantly over the past 20 years, a large portion of GPCR lack reported, experimentally determined structures. Furthermore, many GPCR have no known endogenous and/or synthetic ligands, limiting further exploration of their biochemical, cellular, and physiological roles. While many successes in GPCR ligand discovery have resulted from experimental high-throughput screening, computational methods have played an increasingly important role in GPCR ligand identification in the past decade. Here we discuss computational techniques applied to GPCR ligand discovery. This review summarizes class A GPCR structure/function and provides an overview of many obstacles currently faced in GPCR ligand discovery. Furthermore, we discuss applications and recent successes of computational techniques used to predict GPCR structure as well as present a summary of ligand- and structure-based methods used to identify potential GPCR ligands. Finally, we discuss computational hit list generation and refinement and provide comprehensive workflows for GPCR ligand identification.

All-Atom Molecular Dynamics Simulations Indicated the Involvement of a Conserved Polar Signaling Channel in the Activation Mechanism of the Type I Cannabinoid Receptor

The type I cannabinoid G protein-coupled receptor (CB1, GPCR) is an intensely investigated pharmacological target, owing to its involvement in numerous physiological functions as well as pathological processes such as cancers, neurodegenerative diseases, metabolic disorders and neuropathic pain. In order to develop modern medications that exert their effects through binding to the CB1 receptor, it is essential to understand the structural mechanism of activation of this protein. The pool of atomic resolution experimental structures of GPCRs has been expanding rapidly in the past decade, providing invaluable information about the function of these receptors. According to the current state of the art, the activity of GPCRs involves structurally distinct, dynamically interconverting functional states and the activation is controlled by a cascade of interconnecting conformational switches in the transmembrane domain. A current challenge is to uncover how different functional states are activated and what specific ligand properties are responsible for the selectivity towards those different functional states. Our recent studies of the μ-opioid and β₂-adrenergic receptors (MOP and β₂AR, respectively) revealed that the orthosteric binding pockets and the intracellular surfaces of these receptors are connected through a channel of highly conserved polar amino acids whose dynamic motions are in high correlation in the agonist- and G protein-bound active states. This and independent literature data led us to hypothesize that, in addition to consecutive conformational transitions, a shift of macroscopic polarization takes place in the transmembrane domain, which is furnished by the rearrangement of polar species through their concerted movements. Here, we examined the CB1 receptor signaling complexes utilizing microsecond scale, all-atom molecular dynamics (MD) simulations in order to see if our previous assumptions could be applied to the CB1 receptor too. Apart from the identification of the previously proposed general features of the activation mechanism, several specific properties of the CB1 have been indicated that could possibly be associated with the signaling profile of this receptor.

Mechanism of activation and biased signaling in complement receptor C5aR1

The complement system plays an important role in the innate immune response to invading pathogens. The complement fragment C5a is one of its important effector components and exerts diverse physiological functions through activation of the C5a receptor 1 (C5aR1) and associated downstream G protein and β-arrestin signaling pathways. Dysfunction of the C5a-C5aR1 axis is linked to numerous inflammatory and immune-mediated diseases, but the structural basis for activation and biased signaling of C5aR1 remains elusive. Here, we present cryo-electron microscopy structures of the activated wild-type C5aR1-G_i protein complex bound to each of the following: C5a, the hexapeptidic agonist C5a^pep, and the G protein-biased agonist BM213. The structures reveal the landscape of the C5a-C5aR1 interaction as well as a common motif for the recognition of diverse orthosteric ligands. Moreover, combined with mutagenesis studies and cell-based pharmacological assays, we deciphered a framework for biased signaling using different peptide analogs and provided insight into the activation mechanism of C5aR1 by solving the structure of C5aR1^I116A mutant-G_i signaling activation complex induced by C089, which exerts antagonism on wild-type C5aR1. In addition, unusual conformational changes in the intracellular end of transmembrane domain 7 and helix 8 upon agonist binding suggest a differential signal transduction process. Collectively, our study provides mechanistic understanding into the ligand recognition, biased signaling modulation, activation, and G_i protein coupling of C5aR1, which may facilitate the future design of therapeutic agents.

Allosteric modulation of G protein-coupled receptor signaling

G protein-coupled receptors (GPCRs), the largest family of transmembrane proteins, regulate a wide array of physiological processes in response to extracellular signals. Although these receptors have proven to be the most successful class of drug targets, their complicated signal transduction pathways (including different effector G proteins and β-arrestins) and mediation by orthosteric ligands often cause difficulties for drug development, such as on- or off-target effects. Interestingly, identification of ligands that engage allosteric binding sites, which are different from classic orthosteric sites, can promote pathway-specific effects in cooperation with orthosteric ligands. Such pharmacological properties of allosteric modulators offer new strategies to design safer GPCR-targeted therapeutics for various diseases. Here, we explore recent structural studies of GPCRs bound to allosteric modulators. Our inspection of all GPCR families reveals recognition mechanisms of allosteric regulation. More importantly, this review highlights the diversity of allosteric sites and presents how allosteric modulators control specific GPCR pathways to provide opportunities for the development of new valuable agents.

Structural Insights into ATP-Sensitive Potassium Channel Mechanics: A Role of Intrinsically Disordered Regions

Commonly used techniques, such as CryoEM or X-ray, are not able to capture the structural reorganizations of disordered regions of proteins (IDR); therefore, it is difficult to assess their functions in proteins based exclusively on experiments. To fill this gap, we used computational molecular dynamics (MD) simulation methods to capture IDR dynamics and trace biological function-related interactions in the Kir6.2/SUR1 potassium channel. This ATP-sensitive octameric complex, one of the critical elements in the insulin secretion process in human pancreatic β-cells, has four to five large, disordered fragments. Using unique MD simulations of the full Kir6.2/SUR1 channel complex, we present an in-depth analysis of the dynamics of the disordered regions and discuss the possible functions they could have in this system. Our MD results confirmed the crucial role of the N-terminus of the Kir6.2 fragment and the L0-loop of the SUR1 protein in the transfer of mechanical signals between domains that trigger insulin release. Moreover, we show that the presence of IDRs affects natural ligand binding. Our research takes us one step further toward understanding the action of this vital complex.

Structural Basis for Agonistic Activity and Selectivity toward Melatonin Receptors hMT1 and hMT2

Glaucoma, a major ocular neuropathy originating from a progressive degeneration of retinal ganglion cells, is often associated with increased intraocular pressure (IOP). Daily IOP fluctuations are physiologically influenced by the antioxidant and signaling activities of melatonin. This endogenous modulator has limited employment in treating altered IOP disorders due to its low stability and bioavailability. The search for low-toxic compounds as potential melatonin agonists with higher stability and bioavailability than melatonin itself could start only from knowing the molecular basis of melatonergic activity. Thus, using a computational approach, we studied the melatonin binding toward its natural macromolecular targets, namely melatonin receptors 1 (MT1) and 2 (MT2), both involved in IOP signaling regulation. Besides, agomelatine, a melatonin-derivative agonist and, at the same time, an atypical antidepressant, was also included in the study due to its powerful IOP-lowering effects. For both ligands, we evaluated both stability and ligand positioning inside the orthosteric site of MTs, mapping the main molecular interactions responsible for receptor activation. Affinity values in terms of free binding energy (ΔG_bind) were calculated for the selected poses of the chosen compounds after stabilization through a dynamic molecular docking protocol. The results were compared with experimental in vivo effects, showing a higher potency and more durable effect for agomelatine with respect to melatonin, which could be ascribed both to its higher affinity for hMT2 and to its additional activity as an antagonist for the serotonin receptor 5-HT2c, in agreement with the in silico results.

Correlated evolution of social organization and lifespan in mammals

Discerning the relationship between sociality and longevity would permit a deeper understanding of how animal life history evolved. Here, we perform a phylogenetic comparative analysis of ~1000 mammalian species on three states of social organization (solitary, pair-living, and group-living) and longevity. We show that group-living species generally live longer than solitary species, and that the transition rate from a short-lived state to a long-lived state is higher in group-living than non-group-living species, altogether supporting the correlated evolution of social organization and longevity. The comparative brain transcriptomes of 94 mammalian species identify 31 genes, hormones and immunity-related pathways broadly involved in the association between social organization and longevity. Further selection features reveal twenty overlapping pathways under selection for both social organization and longevity. These results underscore a molecular basis for the influence of the social organization on longevity.

Mechanism of hormone and allosteric agonist mediated activation of follicle stimulating hormone receptor

Follicle stimulating hormone (FSH) is an essential glycoprotein hormone for human reproduction, which functions are mediated by a G protein-coupled receptor, FSHR. Aberrant FSH-FSHR signaling causes infertility and ovarian hyperstimulation syndrome. Here we report cryo-EM structures of FSHR in both inactive and active states, with the active structure bound to FSH and an allosteric agonist compound 21 f. The structures of FSHR are similar to other glycoprotein hormone receptors, highlighting a conserved activation mechanism of hormone-induced receptor activation. Compound 21 f formed extensive interactions with the TMD to directly activate FSHR. Importantly, the unique residue H615^7.42 in FSHR plays an essential role in determining FSHR selectivity for various allosteric agonists. Together, our structures provide a molecular basis of FSH and small allosteric agonist-mediated FSHR activation, which could inspire the design of FSHR-targeted drugs for the treatment of infertility and controlled ovarian stimulation for in vitro fertilization.

Ligand recognition mechanism of the human relaxin family peptide receptor 4 (RXFP4)

Members of the insulin superfamily regulate pleiotropic biological processes through two types of target-specific but structurally conserved peptides, insulin/insulin-like growth factors and relaxin/insulin-like peptides. The latter bind to the human relaxin family peptide receptors (RXFPs). Here, we report three cryo-electron microscopy structures of RXFP4-G_i protein complexes in the presence of the endogenous ligand insulin-like peptide 5 (INSL5) or one of the two small molecule agonists, compound 4 and DC591053. The B chain of INSL5 adopts a single α-helix that penetrates into the orthosteric pocket, while the A chain sits above the orthosteric pocket, revealing a peptide-binding mode previously unknown. Together with mutagenesis and functional analyses, the key determinants responsible for the peptidomimetic agonism and subtype selectivity were identified. Our findings not only provide insights into ligand recognition and subtype selectivity among class A G protein-coupled receptors, but also expand the knowledge of signaling mechanisms in the insulin superfamily.

Structural Understanding of Peptide-Bound G Protein-Coupled Receptors: Peptide-Target Interactions

The activation of G protein-coupled receptors (GPCRs) is triggered by ligand binding to their orthosteric sites, which induces ligand-specific conformational changes. Agonists and antagonists bound to GPCR orthosteric sites provide detailed information on ligand-binding modes. Among these, peptide ligands play an instrumental role in GPCR pharmacology and have attracted increased attention as therapeutic drugs. The recent breakthrough in GPCR structural biology has resulted in the remarkable availability of peptide-bound GPCR complexes. Despite the several structural similarities shared by these receptors, they exhibit distinct features in terms of peptide recognition and receptor activation. From this perspective, we have summarized the current status of peptide-bound GPCR structural complexes, largely focusing on the interactions between the receptor and its peptide ligand at the orthosteric site. In-depth structural investigations have yielded valuable insights into the molecular mechanisms underlying peptide recognition. This study would contribute to the discovery of GPCR peptide drugs with improved therapeutic effects.

Cryo-EM structure of orphan G protein-coupled receptor GPR21

GPR21 belongs to class A orphan G protein-coupled receptor (GPCR). The endogenous ligands for human GPR21 remain unidentified. GPR21 expression is associated with developing type 2 diabetes (T2DM), a multifactorial metabolic disease caused by pancreatic β-cell dysfunction, decreasing insulin production, insulin resistance, and obesity. Animal studies suggested that GPR21 is a potential therapeutic target for T2DM treatment. The underlying mechanisms leading to GPR21 self-activation remain unknown. In our co-expression analysis, we noted that GPR21 could also form a stable complex with an unreported Gα protein subtype, Gαs. To gain further insights into the structural mechanisms of GPR21 activation, we employed cryo-electron microscopy (cryo-EM) and single-particle analysis to resolve the high-resolution structure of GPR21-Gαs complexes. The clear electron density map of the GPR21-Gαs provided direct evidence that GPR21 could couple to Gαs protein at physiological conditions. Thus, GPR21 might mediate previously unexplored pathways in normal or pathological conditions, which warrants further investigation. Structure-guided mutagenesis and biochemical analysis revealed that extracellular loop 2 (ECL2) of GPR21 is essential for the receptor transducing intracellular signal via cAMP. Together, the new structure data reveal a novel signaling cascade of human GPR21 mediated by ECL2 and provide fundamental information for future structure-based drug development.

Structural elucidation of ETHR-A and ETHR-B from Plutella xylostella and insight into non-conservative mutations in transmembrane helix-6

The development of Diamondback moth (DBM) depends on the ecdysis triggering hormone receptor (ETHR); a neuronal membrane G-protein coupled receptor (GPCR) connected to the metamorphosis cascade. Lepidopteran insect DBM is an infamous pest of cruciferous plants. This study examined the full-length coding sequences (CDS) of PxETHR-A and PxETHR-B from the DBM genome. The three-dimensional (3 D) models of both receptors were generated in an inactive state. The behaviour and stability of receptors were examined using molecular dynamics simulations in a lipid membrane system for 300 ns and established a GPCR family-based view. Secondary interactions within receptors were studied to know more about factors contributing to their stability. Multiple sequence alignment revealed conserved features of insect ETHRs those compared to the GPCR family proteins. These features were helpful during the evaluation of the molecular models of both receptors. Side-chain orientation of conserved residues, non-conserved and conserved hydrogen-bond networks (HBN) and hydrophobic clusters were examined in the structures of both receptors. The non-conserved residues L6.35, T6.39, C/S6.43, and L6.48, are present in a conserved position on the transmembrane helix-6 (TM6) of ETHRs. In TM6, PxETHR-A and PxETHR-B differ at positions C/S6.43 and Y/F6.51, both being part of the HBN.Communicated by Ramaswamy H. Sarma.

Cryo-EM structures of orphan GPR21 signaling complexes

GPR21 is a class-A orphan G protein-coupled receptor (GPCR) and a potential therapeutic target for type 2 diabetes and other metabolic disorders. This receptor shows high basal activity in coupling to multiple G proteins in the absence of any known endogenous agonist or synthetic ligand. Here, we present the structures of ligand-free human GPR21 bound to heterotrimeric miniGs and miniG15 proteins, respectively. We identified an agonist-like motif in extracellular loop 2 (ECL2) that occupies the orthosteric pocket and promotes receptor activation. A side pocket that may be employed as a new ligand binding site was also uncovered. Remarkably, G protein binding is accommodated by a flexible cytoplasmic portion of transmembrane helix 6 (TM6) which adopts little or undetectable outward movement. These findings will enable the design of modulators for GPR21 for understanding its signal transduction and exploring opportunity for deorphanization.

Thermodynamic architecture and conformational plasticity of GPCRs

G-protein-coupled receptors (GPCRs) are ubiquitous integral membrane proteins involved in diverse cellular signaling processes. Here, we carry out a large-scale ensemble thermodynamic study of 45 ligand-free GPCRs employing a structure-based statistical mechanical framework. We find that multiple partially structured states co-exist in the GPCR native ensemble, with the TM helices 1, 6 and 7 displaying varied folding status, and shaping the conformational landscape. Strongly coupled residues are anisotropically distributed, accounting for only 13% of the residues, illustrating that a large number of residues are inherently dynamic. Active-state GPCRs are characterized by reduced conformational heterogeneity with altered coupling-patterns distributed throughout the structural scaffold. In silico alanine-scanning mutagenesis reveals that extra- and intra-cellular faces of GPCRs are coupled thermodynamically, highlighting an exquisite structural specialization and the fluid nature of the intramolecular interaction network. The ensemble-based perturbation methodology presented here lays the foundation for understanding allosteric mechanisms and the effects of disease-causing mutations in GCPRs.

AlphaFold2 versus experimental structures: evaluation on G protein-coupled receptors

As important drug targets, G protein-coupled receptors (GPCRs) play pivotal roles in a wide range of physiological processes. Extensive efforts of structural biology have been made on the study of GPCRs. However, a large portion of GPCR structures remain unsolved due to structural instability. Recently, AlphaFold2 has been developed to predict structure models of many functionally important proteins including all members of the GPCR family. Herein we evaluated the accuracy of GPCR structure models predicted by AlphaFold2. We revealed that AlphaFold2 could capture the overall backbone features of the receptors. However, the predicted models and experimental structures were different in many aspects including the assembly of the extracellular and transmembrane domains, the shape of the ligand-binding pockets, and the conformation of the transducer-binding interfaces. These differences impeded the use of predicted structure models in the functional study and structure-based drug design of GPCRs, which required reliable high-resolution structural information.

Cryo-EM structure of the human adenosine A2B receptor-Gs signaling complex

The human adenosine A_2B receptor (A_2BR) is a class A G protein-coupled receptor that is involved in several major physiological and pathological processes throughout the body. A_2BR recognizes its ligands adenosine and NECA with relatively low affinity, but the detailed mechanism for its ligand recognition and signaling is still elusive. Here, we present two structures determined by cryo-electron microscopy of A_2BR bound to its agonists NECA and BAY60-6583, each coupled to an engineered G_s protein. The structures reveal conserved orthosteric binding pockets with subtle differences, whereas the selectivity or specificity can mainly be attributed to regions extended from the orthosteric pocket. We also found that BAY60-6583 occupies a secondary pocket, where residues V250^6.51 and N273^7.36 were two key determinants for its selectivity against A_2BR. This study offers a better understanding of ligand selectivity for the adenosine receptor family and provides a structural template for further development of A_2BR ligands for related diseases.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Insights into divalent cation regulation and G13-coupling of orphan receptor GPR35

Endogenous ions play important roles in the function and pharmacology of G protein-coupled receptors (GPCRs) with limited atomic evidence. In addition, compared with G protein subtypes G_s, G_i/o, and G_q/11, insufficient structural evidence is accessible to understand the coupling mechanism of G_12/13 protein by GPCRs. Orphan receptor GPR35, which is predominantly expressed in the gastrointestinal tract and is closely related to inflammatory bowel diseases (IBDs), stands out as a prototypical receptor for investigating ionic modulation and G₁₃ coupling. Here we report a cryo-electron microscopy structure of G₁₃-coupled GPR35 bound to an anti-allergic drug, lodoxamide. This structure reveals a novel divalent cation coordination site and a unique ionic regulatory mode of GPR35 and also presents a highly positively charged binding pocket and the complementary electrostatic ligand recognition mode, which explain the promiscuity of acidic ligand binding by GPR35. Structural comparison of the GPR35-G₁₃ complex with other G protein subtypes-coupled GPCRs reveals a notable movement of the C-terminus of α5 helix of the Gα₁₃ subunit towards the receptor core and the least outward displacement of the cytoplasmic end of GPR35 TM6. A featured 'methionine pocket' contributes to the G₁₃ coupling by GPR35. Together, our findings provide a structural basis for divalent cation modulation, ligand recognition, and subsequent G₁₃ protein coupling of GPR35 and offer a new opportunity for designing GPR35-targeted drugs for the treatment of IBDs.

Pan-cancer functional analysis of somatic mutations in G protein-coupled receptors

G Protein-coupled receptors (GPCRs) are the most frequently exploited drug target family, moreover they are often found mutated in cancer. Here we used a dataset of mutations found in patient samples derived from the Genomic Data Commons and compared it to the natural human variance as exemplified by data from the 1000 genomes project. We explored cancer-related mutation patterns in all GPCR classes combined and individually. While the location of the mutations across the protein domains did not differ significantly in the two datasets, a mutation enrichment in cancer patients was observed among class-specific conserved motifs in GPCRs such as the Class A "DRY" motif. A Two-Entropy Analysis confirmed the correlation between residue conservation and cancer-related mutation frequency. We subsequently created a ranking of high scoring GPCRs, using a multi-objective approach (Pareto Front Ranking). Our approach was confirmed by re-discovery of established cancer targets such as the LPA and mGlu receptor families, but also discovered novel GPCRs which had not been linked to cancer before such as the P2Y Receptor 10 (P2RY10). Overall, this study presents a list of GPCRs that are amenable to experimental follow up to elucidate their role in cancer.

Building the Chordata Olfactory Receptor Database using more than 400,000 receptors annotated by Genome2OR

Olfactory receptors are poorly annotated for most genome-sequenced chordates. To address this deficiency, we developed a nhmmer-based olfactory receptor annotation tool Genome2OR ( https://github.com/ToHanwei/Genome2OR.git ), and used it to process 1,695 sequenced chordate genomes in the NCBI Assembly database as of January, 2021. In total, 765,248 olfactory receptor genes were annotated, with 404,426 functional genes and 360,822 pseudogenes, which represents a four-fold increase in the number of annotated olfactory receptors. Based on the annotation data, we built a database called Chordata Olfactory Receptor Database (CORD, https://cord.ihuman.shanghaitech.edu.cn ) for archiving, analysing and disseminating the data. Beyond the primary data, we offer derivative information, including pictures of species, cross references to public databases, structural models, sequence similarity networks and sequence profiles in the CORD. Furthermore, we did brief analyses on these receptors, including building a huge protein sequence similarity network covering all receptors in the database, and clustering them into 20 communities, classifying the 20 communities into three categories based on their presences/absences in ray-finned fish and/or lobe-finned fish. We infer that olfactory receptors should have unique activation and desensitization mechanisms by analysing their sequences and structural models. We believe the CORD can benefit the researchers and the general public who are interested in olfaction.

Superconserved receptors expressed in the brain: Expression, function, motifs and evolution of an orphan receptor family

GPR27, GPR85 and GPR173 constitute a small family of G protein-coupled receptors (GPCR) that share the distinctive characteristics of being highly conserved throughout vertebrate evolution and predominantly expressed in the brain. Accordingly, they have been coined as "Superconserved Receptors Expressed in the Brain" (SREB), although their expression profile is more complex than what was originally thought. SREBs have no known validated endogenous ligands and are thus labeled as "orphan" receptors. The investigation of this particular category of uncharacterized receptors holds great promise both in terms of physiology and drug development. In the largest GPCR family, the Rhodopsin-like or Class A, around 100 receptors are considered orphans. Because GPCRs are the most successful source of drug targets, the discovery of a novel function or ligand most likely will lead to significant breakthroughs for the discovery of innovative therapies. The high level of conservation is one of the characteristic features of the SREBs. We propose herein a detailed analysis of the putative evolutionary origin of this family. We highlight the properties that distinguish SREBs from other rhodopsin-like GPCRs. We present the current evidence for these receptors downstream signaling pathways and functions. We discuss the pharmacological challenge for the identification of natural or synthetic ligands of orphan receptors like SREBs. The different SREB-related scientific questions are presented with a highlight on what should be addressed in the near future, including the confirmation of published evidence and their validation as drug targets. In particular, we discuss in which pathological conditions these receptors may be of great relevance to solve unmet medical needs.

Mechanism of enhanced sensitivity of mutated β-adrenergic-like octopamine receptor to amitraz in honeybee Apis mellifera: An insight from MD simulations

BACKGROUND

Amitraz is one of the critical acaricides/insecticides for effective control of pest infestation of Varroa destructor mite, a devastating parasite of Apis mellifera, because of its low toxicity to honeybees. Previous assays verified that a typical G protein-coupled receptor, β-adrenergic-like octopamine receptor (Octβ2R), is the unique target of amitraz, but the honeybee Octβ2R resists to amitraz. However, the underlying molecular mechanism of the enhanced sensitivity or toxicity of amitraz to mutated honeybee Octβ2R^{E208V/I335T/I350V} is not fully understood. Here, molecular dynamics simulations are employed to explore the implied mechanism of the enhanced sensitivity to amitraz in mutant honeybee Octβ2R.

RESULTS

We found that amitraz binding stabilized the structure of Octβ2R, particularly the intracellular loop 3 associated with the Octβ2R signaling. Then, it was further demonstrated that both mutations and ligand binding resulted in a more rigid and compact amitraz binding site, as well as the outward movement of the transmembrane helix 6, which was a prerequisite for G protein coupling and activation. Moreover, mutations were found to promote the binding between Octβ2R and amitraz. Finally, community analysis illuminated that mutations and amitraz strengthened the residue-residue communication within the transmembrane domain, which might facilitate the allosteric signal propagation and activation of Octβ2R.

CONCLUSION

Our results unveiled structural determinants of improved sensitivity in the Octβ2R-amitraz complex and may contribute to further structure-based drug design for safer and less toxic selective insecticides. © 2022 Society of Chemical Industry.

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.

A Single Point Mutation Blocks the Entrance of Ligands to the Cannabinoid CB2 Receptor via the Lipid Bilayer

Molecular dynamic (MD) simulations have become a common tool to study the pathway of ligand entry to the orthosteric binding site of G protein-coupled receptors. Here, we have combined MD simulations and site-directed mutagenesis to study the binding process of the potent JWH-133 agonist to the cannabinoid CB₂ receptor (CB₂R). In CB₂R, the N-terminus and extracellular loop 2 fold over the ligand binding pocket, blocking access to the binding cavity from the extracellular environment. We, thus, hypothesized that the binding pathway is a multistage process consisting of the hydrophobic ligand diffusing in the lipid bilayer to contact a lipid-facing vestibule, from which the ligand enters an allosteric site inside the transmembrane bundle through a tunnel formed between TMs 1 and 7 and finally moving from the allosteric to the orthosteric binding cavity. This pathway was experimentally validated by the Ala282^7.36Phe mutation that blocks the entrance of the ligand, as JWH-133 was not able to decrease the forskolin-induced cAMP levels in cells expressing the mutant receptor. This proposed ligand entry pathway defines transient binding sites that are potential cavities for the design of synthetic modulators.

Multi-state modeling of G-protein coupled receptors at experimental accuracy

The family of G-protein coupled receptors (GPCRs) is one of the largest protein families in the human genome. GPCRs transduct chemical signals from extracellular to intracellular regions via a conformational switch between active and inactive states upon ligand binding. While experimental structures of GPCRs remain limited, high-accuracy computational predictions are now possible with AlphaFold2. However, AlphaFold2 only predicts one state and is biased toward either the active or inactive conformation depending on the GPCR class. Here, a multi-state prediction protocol is introduced that extends AlphaFold2 to predict either active or inactive states at very high accuracy using state-annotated templated GPCR databases. The predicted models accurately capture the main structural changes upon activation of the GPCR at the atomic level. For most of the benchmarked GPCRs (10 out of 15), models in the active and inactive states were closer to their corresponding activation state structures. Median RMSDs of the transmembrane regions were 1.12 Å and 1.41 Å for the active and inactive state models, respectively. The models were more suitable for protein-ligand docking than the original AlphaFold2 models and template-based models. Finally, our prediction protocol predicted accurate GPCR structures and GPCR-peptide complex structures in GPCR Dock 2021, a blind GPCR-ligand complex modeling competition. We expect that high accuracy GPCR models in both activation states will promote understanding in GPCR activation mechanisms and drug discovery for GPCRs. At the time, the new protocol paves the way towards capturing the dynamics of proteins at high-accuracy via machine-learning methods.

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.