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.

Pan-cancer ion transport signature reveals functional regulators of glioblastoma aggression.

Ion channels, transporters, and other ion-flux controlling proteins, collectively comprising the "ion permeome", are common drug targets, however, their roles in cancer remain understudied. Our integrative pan-cancer transcriptome analysis shows that genes encoding the ion permeome are significantly more often highly expressed in specific subsets of cancer samples, compared to pan-transcriptome expectations. To enable target selection, we identified 410 survival-associated IP genes in 33 cancer types using a machine-learning approach. Notably, GJB2 and SCN9A show prominent expression in neoplastic cells and are associated with poor prognosis in glioblastoma, the most common and aggressive brain cancer. GJB2 or SCN9A knockdown in patient-derived glioblastoma cells induces transcriptome-wide changes involving neuron projection and proliferation pathways, impairs cell viability and tumor sphere formation in vitro, perturbs tunneling nanotube dynamics, and extends the survival of glioblastoma-bearing mice. Thus, aberrant activation of genes encoding ion transport proteins appears as a pan-cancer feature defining tumor heterogeneity, which can be exploited for mechanistic insights and therapy development.

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.

Low-dose metformin targets the lysosomal AMPK pathway through PEN2.

Metformin, the most prescribed antidiabetic medicine, has shown other benefits such as anti-ageing and anticancer effects1-4. For clinical doses of metformin, AMP-activated protein kinase (AMPK) has a major role in its mechanism of action4,5; however, the direct molecular target of metformin remains unknown. Here we show that clinically relevant concentrations of metformin inhibit the lysosomal proton pump v-ATPase, which is a central node for AMPK activation following glucose starvation6. We synthesize a photoactive metformin probe and identify PEN2, a subunit of γ-secretase7, as a binding partner of metformin with a dissociation constant at micromolar levels. Metformin-bound PEN2 forms a complex with ATP6AP1, a subunit of the v-ATPase8, which leads to the inhibition of v-ATPase and the activation of AMPK without effects on cellular AMP levels. Knockout of PEN2 or re-introduction of a PEN2 mutant that does not bind ATP6AP1 blunts AMPK activation. In vivo, liver-specific knockout of Pen2 abolishes metformin-mediated reduction of hepatic fat content, whereas intestine-specific knockout of Pen2 impairs its glucose-lowering effects. Furthermore, knockdown of pen-2 in Caenorhabditis elegans abrogates metformin-induced extension of lifespan. Together, these findings reveal that metformin binds PEN2 and initiates a signalling route that intersects, through ATP6AP1, the lysosomal glucose-sensing pathway for AMPK activation. This ensures that metformin exerts its therapeutic benefits in patients without substantial adverse effects.

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.

Acetylation coordinates the crosstalk between carbon metabolism and ammonium assimilation in Salmonella enterica.

Enteric bacteria use up to 15% of their cellular energy for ammonium assimilation via glutamine synthetase (GS)/glutamate synthase (GOGAT) and glutamate dehydrogenase (GDH) in response to varying ammonium availability. However, the sensory mechanisms for effective and appropriate coordination between carbon metabolism and ammonium assimilation have not been fully elucidated. Here, we report that in Salmonella enterica, carbon metabolism coordinates the activities of GS/GDH via functionally reversible protein lysine acetylation. Glucose promotes Pat acetyltransferase-mediated acetylation and activation of adenylylated GS. Simultaneously, glucose induces GDH acetylation to inactivate the enzyme by impeding its catalytic centre, which is reversed upon GDH deacetylation by deacetylase CobB. Molecular dynamics (MD) simulations indicate that adenylylation is required for acetylation-dependent activation of GS. We show that acetylation and deacetylation occur within minutes of "glucose shock" to promptly adapt to ammonium/carbon variation and finely balance glutamine/glutamate synthesis. Finally, in a mouse infection model, reduced S. enterica growth caused by the expression of adenylylation-mimetic GS is rescued by acetylation-mimicking mutations. Thus, glucose-driven acetylation integrates signals from ammonium assimilation and carbon metabolism to fine-tune bacterial growth control.

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.

Paralytic, the Drosophila voltage-gated sodium channel, regulates proliferation of neural progenitors.

Proliferating cells, typically considered "nonexcitable," nevertheless, exhibit regulation by bioelectric signals. Notably, voltage-gated sodium channels (VGSC) that are crucial for neuronal excitability are also found in progenitors and up-regulated in cancer. Here, we identify a role for VGSC in proliferation of Drosophila neuroblast (NB) lineages within the central nervous system. Loss of paralytic (para), the sole gene that encodes Drosophila VGSC, reduces neuroblast progeny cell number. The type II neuroblast lineages, featuring a population of transit-amplifying intermediate neural progenitors (INP) similar to that found in the developing human cortex, are particularly sensitive to para manipulation. Following a series of asymmetric divisions, INPs normally exit the cell cycle through a final symmetric division. Our data suggests that loss of Para induces apoptosis in this population, whereas overexpression leads to an increase in INPs and overall neuroblast progeny cell numbers. These effects are cell autonomous and depend on Para channel activity. Reduction of Para expression not only affects normal NB development, but also strongly suppresses brain tumor mass, implicating a role for Para in cancer progression. To our knowledge, our studies are the first to identify a role for VGSC in neural progenitor proliferation. Elucidating the contribution of VGSC in proliferation will advance our understanding of bioelectric signaling within development and disease states.

Directed differentiation of human hindbrain neuroepithelial stem cells recapitulates cerebellar granule neurogenesis.

Cerebellar granule neurons (CGNs) are the most abundant neurons in the human brain. Dysregulation of their development underlies movement disorders and medulloblastomas. It is suspected that these disorders arise in progenitor states of the CGN lineage, for which human models are lacking. Here, we have differentiated human hindbrain neuroepithelial stem (hbNES) cells to CGNs in vitro using soluble growth factors, recapitulating key progenitor states in the lineage. We show that hbNES cells are not lineage committed and retain rhombomere 1 regional identity. Upon differentiation, hbNES cells transit through a rhombic lip (RL) progenitor state at day 7, demonstrating human specific sub-ventricular cell identities. This RL state is followed by an ATOH1+ CGN progenitor state at day 14. By the end of a 56-day differentiation procedure, we obtain functional neurons expressing CGN markers GABAARα6 and vGLUT2. We show that sonic hedgehog promotes GABAergic lineage specification and CGN progenitor proliferation. Our work presents a new model with which to study development and diseases of the CGN lineage in a human context.

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.

Indentation induces instantaneous nuclear stiffening and unfolding of nuclear envelope wrinkles.

The nuclear envelope (NE) separates genomic DNA from the cytoplasm and regulates transport between the cytosol and the nucleus in eukaryotes. Nuclear stiffening enables the cell nucleus to protect itself from extensive deformation, loss of NE integrity, and genome instability. It is known that the reorganization of actin, lamin, and chromatin can contribute to nuclear stiffening. In this work, we show that structural alteration of NE also contributes to instantaneous nuclear stiffening under indentation. In situ mechanical characterization of cell nuclei in intact cells shows that nuclear stiffening and unfolding of NE wrinkles occur simultaneously at the indentation site. A positive correlation between the initial state of NE wrinkles, the unfolding of NE wrinkles, and the stiffening ratio (stiffness fold-change) is found. Additionally, NE wrinkles unfold throughout the nucleus outside the indentation site. Finite element simulation, which involves the purely passive process of structural unfolding, shows that unfolding of NE wrinkles alone can lead to an increase in nuclear stiffness and a reduction in stress and strain levels. Together, these results provide a perspective on how cell nucleus adapts to mechanical stimuli through structural alteration of the NE.

The photomorphogenic repressors BBX28 and BBX29 integrate light and brassinosteroid signaling to inhibit seedling development in Arabidopsis.

B-box containing proteins (BBXs) integrate light and various hormonal signals to regulate plant growth and development. Here, we demonstrate that the photomorphogenic repressors BBX28 and BBX29 positively regulate brassinosteroid (BR) signaling in Arabidopsis thaliana seedlings. Treatment with the BR brassinolide stabilized BBX28 and BBX29, which partially depended on BR INSENSITIVE1 (BRI1) and BIN2. bbx28 bbx29 seedlings exhibited larger cotyledon aperture than the wild-type when treated with brassinazole in the dark, which partially suppressed the closed cotyledons of brassinazole resistant 1-1D (bzr1-1D). Consistently, overexpressing BBX28 and BBX29 partially rescued the short hypocotyls of bri1-5 and bin2-1 in both the dark and light, while the loss-of-function of BBX28 and BBX29 partially suppressed the long hypocotyls of bzr1-1D in the light. BBX28 and BBX29 physically interacted with BR-ENHANCED EXPRESSION1 (BEE1), BEE2, and BEE3 and enhanced their binding to and activation of their target genes. Moreover, BBX28 and BBX29 as well as BEE1, BEE2, and BEE3 increased BZR1 accumulation to promote the BR signaling pathway. Therefore, both BBX28 and BBX29 interact with BEE1, BEE2, and BEE3 to orchestrate light and BR signaling by facilitating the transcriptional activity of BEE target genes. Our study provides insights into the pivotal roles of BBX28 and BBX29 as signal integrators in ensuring normal seedling development.