Structure-based discovery of conformationally selective inhibitors of the serotonin transporter.

The serotonin transporter (SERT) removes synaptic serotonin and is the target of anti-depressant drugs. SERT adopts three conformations: outward-open, occluded, and inward-open. All known inhibitors target the outward-open state except ibogaine, which has unusual anti-depressant and substance-withdrawal effects, and stabilizes the inward-open conformation. Unfortunately, ibogaine's promiscuity and cardiotoxicity limit the understanding of inward-open state ligands. We docked over 200 million small molecules against the inward-open state of the SERT. Thirty-six top-ranking compounds were synthesized, and thirteen inhibited; further structure-based optimization led to the selection of two potent (low nanomolar) inhibitors. These stabilized an outward-closed state of the SERT with little activity against common off-targets. A cryo-EM structure of one of these bound to the SERT confirmed the predicted geometry. In mouse behavioral assays, both compounds had anxiolytic- and anti-depressant-like activity, with potencies up to 200-fold better than fluoxetine (Prozac), and one substantially reversed morphine withdrawal effects.

Structural insights into the human D1 and D2 dopamine receptor signaling complexes.

The D1- and D2-dopamine receptors (D1R and D2R), which signal through Gs and Gi, respectively, represent the principal stimulatory and inhibitory dopamine receptors in the central nervous system. D1R and D2R also represent the main therapeutic targets for Parkinson's disease, schizophrenia, and many other neuropsychiatric disorders, and insight into their signaling is essential for understanding both therapeutic and side effects of dopaminergic drugs. Here, we report four cryoelectron microscopy (cryo-EM) structures of D1R-Gs and D2R-Gi signaling complexes with selective and non-selective dopamine agonists, including two currently used anti-Parkinson's disease drugs, apomorphine and bromocriptine. These structures, together with mutagenesis studies, reveal the conserved binding mode of dopamine agonists, the unique pocket topology underlying ligand selectivity, the conformational changes in receptor activation, and potential structural determinants for G protein-coupling selectivity. These results provide both a molecular understanding of dopamine signaling and multiple structural templates for drug design targeting the dopaminergic system.

Discovery of Human Signaling Systems: Pairing Peptides to G Protein-Coupled Receptors.

The peptidergic system is the most abundant network of ligand-receptor-mediated signaling in humans. However, the physiological roles remain elusive for numerous peptides and more than 100 G protein-coupled receptors (GPCRs). Here we report the pairing of cognate peptides and receptors. Integrating comparative genomics across 313 species and bioinformatics on all protein sequences and structures of human class A GPCRs, we identify universal characteristics that uncover additional potential peptidergic signaling systems. Using three orthogonal biochemical assays, we pair 17 proposed endogenous ligands with five orphan GPCRs that are associated with diseases, including genetic, neoplastic, nervous and reproductive system disorders. We also identify additional peptides for nine receptors with recognized ligands and pathophysiological roles. This integrated computational and multifaceted experimental approach expands the peptide-GPCR network and opens the way for studies to elucidate the roles of these signaling systems in human physiology and disease. VIDEO ABSTRACT.

5-HT2C Receptor Structures Reveal the Structural Basis of GPCR Polypharmacology.

Drugs frequently require interactions with multiple targets-via a process known as polypharmacology-to achieve their therapeutic actions. Currently, drugs targeting several serotonin receptors, including the 5-HT2C receptor, are useful for treating obesity, drug abuse, and schizophrenia. The competing challenges of developing selective 5-HT2C receptor ligands or creating drugs with a defined polypharmacological profile, especially aimed at G protein-coupled receptors (GPCRs), remain extremely difficult. Here, we solved two structures of the 5-HT2C receptor in complex with the highly promiscuous agonist ergotamine and the 5-HT2A-C receptor-selective inverse agonist ritanserin at resolutions of 3.0 Å and 2.7 Å, respectively. We analyzed their respective binding poses to provide mechanistic insights into their receptor recognition and opposing pharmacological actions. This study investigates the structural basis of polypharmacology at canonical GPCRs and illustrates how understanding characteristic patterns of ligand-receptor interaction and activation may ultimately facilitate drug design at multiple GPCRs.

Structure of the Nanobody-Stabilized Active State of the Kappa Opioid Receptor.

The κ-opioid receptor (KOP) mediates the actions of opioids with hallucinogenic, dysphoric, and analgesic activities. The design of KOP analgesics devoid of hallucinatory and dysphoric effects has been hindered by an incomplete structural and mechanistic understanding of KOP agonist actions. Here, we provide a crystal structure of human KOP in complex with the potent epoxymorphinan opioid agonist MP1104 and an active-state-stabilizing nanobody. Comparisons between inactive- and active-state opioid receptor structures reveal substantial conformational changes in the binding pocket and intracellular and extracellular regions. Extensive structural analysis and experimental validation illuminate key residues that propagate larger-scale structural rearrangements and transducer binding that, collectively, elucidate the structural determinants of KOP pharmacology, function, and biased signaling. These molecular insights promise to accelerate the structure-guided design of safer and more effective κ-opioid receptor therapeutics.

Structures of the human dopamine D3 receptor-Gi complexes.

The dopamine system, including five dopamine receptors (D1R-D5R), plays essential roles in the central nervous system (CNS), and ligands that activate dopamine receptors have been used to treat many neuropsychiatric disorders. Here, we report two cryo-EM structures of human D3R in complex with an inhibitory G protein and bound to the D3R-selective agonists PD128907 and pramipexole, the latter of which is used to treat patients with Parkinson's disease. The structures reveal agonist binding modes distinct from the antagonist-bound D3R structure and conformational signatures for ligand-induced receptor activation. Mutagenesis and homology modeling illuminate determinants of ligand specificity across dopamine receptors and the mechanisms for Gi protein coupling. Collectively our work reveals the basis of agonist binding and ligand-induced receptor activation and provides structural templates for designing specific ligands to treat CNS diseases targeting the dopaminergic system.

Molecular basis for selective activation of DREADD-based chemogenetics.

Designer receptors exclusively activated by designer drugs (DREADDs) represent a powerful chemogenetic technology for the remote control of neuronal activity and cellular signalling1-4. The muscarinic receptor-based DREADDs are the most widely used chemogenetic tools in neuroscience research. The Gq-coupled DREADD (hM3Dq) is used to enhance neuronal activity, whereas the Gi/o-coupled DREADD (hM4Di) is utilized to inhibit neuronal activity5. Here we report four DREADD-related cryogenic electron microscopy high-resolution structures: a hM3Dq-miniGq complex and a hM4Di-miniGo complex bound to deschloroclozapine; a hM3Dq-miniGq complex bound to clozapine-N-oxide; and a hM3R-miniGq complex bound to iperoxo. Complemented with mutagenesis, functional and computational simulation data, our structures reveal key details of the recognition of DREADD chemogenetic actuators and the molecular basis for activation. These findings should accelerate the structure-guided discovery of next-generation chemogenetic tools.

Synthon-based ligand discovery in virtual libraries of over 11 billion compounds.

Structure-based virtual ligand screening is emerging as a key paradigm for early drug discovery owing to the availability of high-resolution target structures1-4 and ultra-large libraries of virtual compounds5,6. However, to keep pace with the rapid growth of virtual libraries, such as readily available for synthesis (REAL) combinatorial libraries7, new approaches to compound screening are needed8,9. Here we introduce a modular synthon-based approach-V-SYNTHES-to perform hierarchical structure-based screening of a REAL Space library of more than 11 billion compounds. V-SYNTHES first identifies the best scaffold-synthon combinations as seeds suitable for further growth, and then iteratively elaborates these seeds to select complete molecules with the best docking scores. This hierarchical combinatorial approach enables the rapid detection of the best-scoring compounds in the gigascale chemical space while performing docking of only a small fraction (<0.1%) of the library compounds. Chemical synthesis and experimental testing of novel cannabinoid antagonists predicted by V-SYNTHES demonstrated a 33% hit rate, including 14 submicromolar ligands, substantially improving over a standard virtual screening of the Enamine REAL diversity subset, which required approximately 100 times more computational resources. Synthesis of selected analogues of the best hits further improved potencies and affinities (best inhibitory constant (Ki) = 0.9 nM) and CB2/CB1 selectivity (50-200-fold). V-SYNTHES was also tested on a kinase target, ROCK1, further supporting its use for lead discovery. The approach is easily scalable for the rapid growth of combinatorial libraries and potentially adaptable to any docking algorithm.

Bespoke library docking for 5-HT2A receptor agonists with antidepressant activity.

There is considerable interest in screening ultralarge chemical libraries for ligand discovery, both empirically and computationally1-4. Efforts have focused on readily synthesizable molecules, inevitably leaving many chemotypes unexplored. Here we investigate structure-based docking of a bespoke virtual library of tetrahydropyridines-a scaffold that is poorly sampled by a general billion-molecule virtual library but is well suited to many aminergic G-protein-coupled receptors. Using three inputs, each with diverse available derivatives, a one pot C-H alkenylation, electrocyclization and reduction provides the tetrahydropyridine core with up to six sites of derivatization5-7. Docking a virtual library of 75 million tetrahydropyridines against a model of the serotonin 5-HT2A receptor (5-HT2AR) led to the synthesis and testing of 17 initial molecules. Four of these molecules had low-micromolar activities against either the 5-HT2A or the 5-HT2B receptors. Structure-based optimization led to the 5-HT2AR agonists (R)-69 and (R)-70, with half-maximal effective concentration values of 41 nM and 110 nM, respectively, and unusual signalling kinetics that differ from psychedelic 5-HT2AR agonists. Cryo-electron microscopy structural analysis confirmed the predicted binding mode to 5-HT2AR. The favourable physical properties of these new agonists conferred high brain permeability, enabling mouse behavioural assays. Notably, neither had psychedelic activity, in contrast to classic 5-HT2AR agonists, whereas both had potent antidepressant activity in mouse models and had the same efficacy as antidepressants such as fluoxetine at as low as 1/40th of the dose. Prospects for using bespoke virtual libraries to sample pharmacologically relevant chemical space will be considered.

Structural insights into angiotensin receptor signaling modulation by balanced and biased agonists.

The peptide hormone angiotensin II regulates blood pressure mainly through the type 1 angiotensin II receptor AT1 R and its downstream signaling proteins Gq and ß-arrestin. AT1 R blockers, clinically used as antihypertensive drugs, inhibit both signaling pathways, whereas AT1 R ß-arrestin-biased agonists have shown great potential for the treatment of acute heart failure. Here, we present a cryo-electron microscopy (cryo-EM) structure of the human AT1 R in complex with a balanced agonist, Sar1 -AngII, and Gq protein at 2.9 Å resolution. This structure, together with extensive functional assays and computational modeling, reveals the molecular mechanisms for AT1 R signaling modulation and suggests that a major hydrogen bond network (MHN) inside the receptor serves as a key regulator of AT1 R signal transduction from the ligand-binding pocket to both Gq and ß-arrestin pathways. Specifically, we found that the MHN mutations N1113.35 A and N2947.45 A induce biased signaling to Gq and ß-arrestin, respectively. These insights should facilitate AT1 R structure-based drug discovery for the treatment of cardiovascular diseases.

Structures of the σ2 receptor enable docking for bioactive ligand discovery.

The σ2 receptor has attracted intense interest in cancer imaging1, psychiatric disease2, neuropathic pain3-5 and other areas of biology6,7. Here we determined the crystal structure of this receptor in complex with the clinical candidate roluperidone2 and the tool compound PB288. These structures templated a large-scale docking screen of 490 million virtual molecules, of which 484 compounds were synthesized and tested. We identified 127 new chemotypes with affinities superior to 1 µM, 31 of which had affinities superior to 50 nM. The hit rate fell smoothly and monotonically with docking score. We optimized three hits for potency and selectivity, and achieved affinities that ranged from 3 to 48 nM, with up to 250-fold selectivity versus the σ1 receptor. Crystal structures of two ligands bound to the σ2 receptor confirmed the docked poses. To investigate the contribution of the σ2 receptor in pain, two potent σ2-selective ligands and one potent σ1/σ2 non-selective ligand were tested for efficacy in a mouse model of neuropathic pain. All three ligands showed time-dependent decreases in mechanical hypersensitivity in the spared nerve injury model9, suggesting that the σ2 receptor has a role in nociception. This study illustrates the opportunities for rapid discovery of in vivo probes through structure-based screens of ultra large libraries, enabling study of underexplored areas of biology.

Structure, function and pharmacology of human itch GPCRs.

The MRGPRX family of receptors (MRGPRX1-4) is a family of mas-related G-protein-coupled receptors that have evolved relatively recently1. Of these, MRGPRX2 and MRGPRX4 are key physiological and pathological mediators of itch and related mast cell-mediated hypersensitivity reactions2-5. MRGPRX2 couples to both Gi and Gq in mast cells6. Here we describe agonist-stabilized structures of MRGPRX2 coupled to Gi1 and Gq in ternary complexes with the endogenous peptide cortistatin-14 and with a synthetic agonist probe, respectively, and the development of potent antagonist probes for MRGPRX2. We also describe a specific MRGPRX4 agonist and the structure of this agonist in a complex with MRGPRX4 and Gq. Together, these findings should accelerate the structure-guided discovery of therapeutic agents for pain, itch and mast cell-mediated hypersensitivity.

Virtual discovery of melatonin receptor ligands to modulate circadian rhythms.

The neuromodulator melatonin synchronizes circadian rhythms and related physiological functions through the actions of two G-protein-coupled receptors: MT1 and MT2. Circadian release of melatonin at night from the pineal gland activates melatonin receptors in the suprachiasmatic nucleus of the hypothalamus, synchronizing the physiology and behaviour of animals to the light-dark cycle1-4. The two receptors are established drug targets for aligning circadian phase to this cycle in disorders of sleep5,6 and depression1-4,7-9. Despite their importance, few in vivo active MT1-selective ligands have been reported2,8,10-12, hampering both the understanding of circadian biology and the development of targeted therapeutics. Here we docked more than 150 million virtual molecules to an MT1 crystal structure, prioritizing structural fit and chemical novelty. Of these compounds, 38 high-ranking molecules were synthesized and tested, revealing ligands with potencies ranging from 470 picomolar to 6 micromolar. Structure-based optimization led to two selective MT1 inverse agonists-which were topologically unrelated to previously explored chemotypes-that acted as inverse agonists in a mouse model of circadian re-entrainment. Notably, we found that these MT1-selective inverse agonists advanced the phase of the mouse circadian clock by 1.3-1.5 h when given at subjective dusk, an agonist-like effect that was eliminated in MT1- but not in MT2-knockout mice. This study illustrates the opportunities for modulating melatonin receptor biology through MT1-selective ligands and for the discovery of previously undescribed, in vivo active chemotypes from structure-based screens of diverse, ultralarge libraries.

Ligand recognition and allosteric modulation of the human MRGPRX1 receptor.

The human MAS-related G protein-coupled receptor X1 (MRGPRX1) is preferentially expressed in the small-diameter primary sensory neurons and involved in the mediation of nociception and pruritus. Central activation of MRGPRX1 by the endogenous opioid peptide fragment BAM8-22 and its positive allosteric modulator ML382 has been shown to effectively inhibit persistent pain, making MRGPRX1 a promising target for non-opioid pain treatment. However, the activation mechanism of MRGPRX1 is still largely unknown. Here we report three high-resolution cryogenic electron microscopy structures of MRGPRX1-Gαq in complex with BAM8-22 alone, with BAM8-22 and ML382 simultaneously as well as with a synthetic agonist compound-16. These structures reveal the agonist binding mode for MRGPRX1 and illuminate the structural requirements for positive allosteric modulation. Collectively, our findings provide a molecular understanding of the activation and allosteric modulation of the MRGPRX1 receptor, which could facilitate the structure-based design of non-opioid pain-relieving drugs.

Publisher Correction: Structural basis of ligand recognition at the human MT1 melatonin receptor.

Change history: In this Letter, the rotation signs around 90°, 135° and 15° were missing and in the HTML, Extended Data Tables 2 and 3 were the wrong tables; these errors have been corrected online.

XFEL structures of the human MT2 melatonin receptor reveal the basis of subtype selectivity.

The human MT1 and MT2 melatonin receptors1,2 are G-protein-coupled receptors (GPCRs) that help to regulate circadian rhythm and sleep patterns3. Drug development efforts have targeted both receptors for the treatment of insomnia, circadian rhythm and mood disorders, and cancer3, and MT2 has also been implicated in type 2 diabetes4,5. Here we report X-ray free electron laser (XFEL) structures of the human MT2 receptor in complex with the agonists 2-phenylmelatonin (2-PMT) and ramelteon6 at resolutions of 2.8 Å and 3.3 Å, respectively, along with two structures of function-related mutants: H2085.46A (superscripts represent the Ballesteros-Weinstein residue numbering nomenclature7) and N862.50D, obtained in complex with 2-PMT. Comparison of the structures of MT2 with a published structure8 of MT1 reveals that, despite conservation of the orthosteric ligand-binding site residues, there are notable conformational variations as well as differences in [3H]melatonin dissociation kinetics that provide insights into the selectivity between melatonin receptor subtypes. A membrane-buried lateral ligand entry channel is observed in both MT1 and MT2, but in addition the MT2 structures reveal a narrow opening towards the solvent in the extracellular part of the receptor. We provide functional and kinetic data that support a prominent role for intramembrane ligand entry in both receptors, and suggest that there might also be an extracellular entry path in MT2. Our findings contribute to a molecular understanding of melatonin receptor subtype selectivity and ligand access modes, which are essential for the design of highly selective melatonin tool compounds and therapeutic agents.

Structural basis of ligand recognition at the human MT1 melatonin receptor.

Melatonin (N-acetyl-5-methoxytryptamine) is a neurohormone that maintains circadian rhythms1 by synchronization to environmental cues and is involved in diverse physiological processes2 such as the regulation of blood pressure and core body temperature, oncogenesis, and immune function3. Melatonin is formed in the pineal gland in a light-regulated manner4 by enzymatic conversion from 5-hydroxytryptamine (5-HT or serotonin), and modulates sleep and wakefulness5 by activating two high-affinity G-protein-coupled receptors, type 1A (MT1) and type 1B (MT2)3,6. Shift work, travel, and ubiquitous artificial lighting can disrupt natural circadian rhythms; as a result, sleep disorders affect a substantial population in modern society and pose a considerable economic burden7. Over-the-counter melatonin is widely used to alleviate jet lag and as a safer alternative to benzodiazepines and other sleeping aids8,9, and is one of the most popular supplements in the United States10. Here, we present high-resolution room-temperature X-ray free electron laser (XFEL) structures of MT1 in complex with four agonists: the insomnia drug ramelteon11, two melatonin analogues, and the mixed melatonin-serotonin antidepressant agomelatine12,13. The structure of MT2 is described in an accompanying paper14. Although the MT1 and 5-HT receptors have similar endogenous ligands, and agomelatine acts on both receptors, the receptors differ markedly in the structure and composition of their ligand pockets; in MT1, access to the ligand pocket is tightly sealed from solvent by extracellular loop 2, leaving only a narrow channel between transmembrane helices IV and V that connects it to the lipid bilayer. The binding site is extremely compact, and ligands interact with MT1 mainly by strong aromatic stacking with Phe179 and auxiliary hydrogen bonds with Asn162 and Gln181. Our structures provide an unexpected example of atypical ligand entry for a non-lipid receptor, lay the molecular foundation of ligand recognition by melatonin receptors, and will facilitate the design of future tool compounds and therapeutic agents, while their comparison to 5-HT receptors yields insights into the evolution and polypharmacology of G-protein-coupled receptors.

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.

De Novo Design of κ-Opioid Receptor Antagonists Using a Generative Deep-Learning Framework.

Likely effective pharmacological interventions for the treatment of opioid addiction include attempts to attenuate brain reward deficits during periods of abstinence. Pharmacological blockade of the κ-opioid receptor (KOR) has been shown to abolish brain reward deficits in rodents during withdrawal, as well as to reduce the escalation of opioid use in rats with extended access to opioids. Although KOR antagonists represent promising candidates for the treatment of opioid addiction, very few potent selective KOR antagonists are known to date and most of them exhibit significant safety concerns. Here, we used a generative deep-learning framework for the de novo design of chemotypes with putative KOR antagonistic activity. Molecules generated by models trained with this framework were prioritized for chemical synthesis based on their predicted optimal interactions with the receptor. Our models and proposed training protocol were experimentally validated by binding and functional assays.

Structural optimizations and bioevaluation of N-H aporphine analogues as Gq-biased and selective serotonin 5-HT2C receptor agonists.

The concept of subtype selectivity and functional bias has recently reshaped current GPCR drug discovery for G protein-coupled receptors. A series of new N-H aporphines with A-ring modifications have been synthesized, and their efficacy on 5-HT2 receptor activation was evaluated. SAR analysis led to the discovery of several more potent and selective 5-HT2C receptor agonists (e.g., 11b and 11f) with low nanomolar activity. Molecular docking analysis of this series of aporphines was in accordance with our SAR results. The functional selectivity of specific compounds was tested via both Gq-mediated calcium flux and ß-arrestin recruitment assays, which revealed that these compounds exhibited no ß-arrestin recruitment activity. Further ADMET study combined with behavioral assessment using a methamphetamine-induced hyperactivity model identified compound 11b and 11f possessing promising drug-like profiles and having antipsychotic potential. These agonists with an exclusive bias toward Gq signaling may serve as valuable pharmacological probes to facilitate the elucidation of therapeutically relevant 5-HT2C signaling pathways and the development of alternative antipsychotic medications.