Structure-based discovery of CFTR potentiators and inhibitors.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a crucial ion channel whose loss of function leads to cystic fibrosis, whereas its hyperactivation leads to secretory diarrhea. Small molecules that improve CFTR folding (correctors) or function (potentiators) are clinically available. However, the only potentiator, ivacaftor, has suboptimal pharmacokinetics and inhibitors have yet to be clinically developed. Here, we combine molecular docking, electrophysiology, cryo-EM, and medicinal chemistry to identify CFTR modulators. We docked ∼155 million molecules into the potentiator site on CFTR, synthesized 53 test ligands, and used structure-based optimization to identify candidate modulators. This approach uncovered mid-nanomolar potentiators, as well as inhibitors, that bind to the same allosteric site. These molecules represent potential leads for the development of more effective drugs for cystic fibrosis and secretory diarrhea, demonstrating the feasibility of large-scale docking for ion channel drug discovery.

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.

Structure of a Hallucinogen-Activated Gq-Coupled 5-HT2A Serotonin Receptor.

Hallucinogens like lysergic acid diethylamide (LSD), psilocybin, and substituted N-benzyl phenylalkylamines are widely used recreationally with psilocybin being considered as a therapeutic for many neuropsychiatric disorders including depression, anxiety, and substance abuse. How psychedelics mediate their actions-both therapeutic and hallucinogenic-are not understood, although activation of the 5-HT2A serotonin receptor (HTR2A) is key. To gain molecular insights into psychedelic actions, we determined the active-state structure of HTR2A bound to 25-CN-NBOH-a prototypical hallucinogen-in complex with an engineered Gαq heterotrimer by cryoelectron microscopy (cryo-EM). We also obtained the X-ray crystal structures of HTR2A complexed with the arrestin-biased ligand LSD or the inverse agonist methiothepin. Comparisons of these structures reveal determinants responsible for HTR2A-Gαq protein interactions as well as the conformational rearrangements involved in active-state transitions. Given the potential therapeutic actions of hallucinogens, these findings could accelerate the discovery of more selective drugs for the treatment of a variety of neuropsychiatric disorders.

The Psychiatric Cell Map Initiative: A Convergent Systems Biological Approach to Illuminating Key Molecular Pathways in Neuropsychiatric Disorders.

Although gene discovery in neuropsychiatric disorders, including autism spectrum disorder, intellectual disability, epilepsy, schizophrenia, and Tourette disorder, has accelerated, resulting in a large number of molecular clues, it has proven difficult to generate specific hypotheses without the corresponding datasets at the protein complex and functional pathway level. Here, we describe one path forward-an initiative aimed at mapping the physical and genetic interaction networks of these conditions and then using these maps to connect the genomic data to neurobiology and, ultimately, the clinic. These efforts will include a team of geneticists, structural biologists, neurobiologists, systems biologists, and clinicians, leveraging a wide array of experimental approaches and creating a collaborative infrastructure necessary for long-term investigation. This initiative will ultimately intersect with parallel studies that focus on other diseases, as there is a significant overlap with genes implicated in cancer, infectious disease, and congenital heart defects.

Crystal Structure of an LSD-Bound Human Serotonin Receptor.

The prototypical hallucinogen LSD acts via serotonin receptors, and here we describe the crystal structure of LSD in complex with the human serotonin receptor 5-HT2B. The complex reveals conformational rearrangements to accommodate LSD, providing a structural explanation for the conformational selectivity of LSD's key diethylamide moiety. LSD dissociates exceptionally slow from both 5-HT2BR and 5-HT2AR-a major target for its psychoactivity. Molecular dynamics (MD) simulations suggest that LSD's slow binding kinetics may be due to a "lid" formed by extracellular loop 2 (EL2) at the entrance to the binding pocket. A mutation predicted to increase the mobility of this lid greatly accelerates LSD's binding kinetics and selectively dampens LSD-mediated ß-arrestin2 recruitment. This study thus reveals an unexpected binding mode of LSD; illuminates key features of its kinetics, stereochemistry, and signaling; and provides a molecular explanation for LSD's actions at human serotonin receptors. PAPERCLIP.

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.

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.

Profiling the proximal proteome of the activated µ-opioid receptor.

The µ-opioid receptor (µOR) represents an important target of therapeutic and abused drugs. So far, most understanding of µOR activity has focused on a subset of known signal transducers and regulatory molecules. Yet µOR signaling is coordinated by additional proteins in the interaction network of the activated receptor, which have largely remained invisible given the lack of technologies to interrogate these networks systematically. Here we describe a proteomics and computational approach to map the proximal proteome of the activated µOR and to extract subcellular location, trafficking and functional partners of G-protein-coupled receptor (GPCR) activity. We demonstrate that distinct opioid agonists exert differences in the µOR proximal proteome mediated by endocytosis and endosomal sorting. Moreover, we identify two new µOR network components, EYA4 and KCTD12, which are recruited on the basis of receptor-triggered G-protein activation and might form a previously unrecognized buffering system for G-protein activity broadly modulating cellular GPCR signaling.

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.