Crystal structure of the human 5-HT1B serotonin receptor bound to an inverse agonist

Structural insights into the unexpected agonism of tetracyclic antidepressants at serotonin receptors 5-HT1eR and 5-HT1FR

Serotonin [5-hydroxytryptamine (5-HT)] acts via 13 different receptors in humans. Of these receptor subtypes, all but 5-HT1eR have confirmed roles in native tissue and are validated drug targets. Despite 5-HT1eR's therapeutic potential and plausible druggability, the mechanisms of its activation remain elusive. To illuminate 5-HT1eR's pharmacology in relation to the highly homologous 5-HT1FR, we screened a library of aminergic receptor ligands at both receptors and observe 5-HT1eR/5-HT1FR agonism by multicyclic drugs described as pan-antagonists at 5-HT receptors. Potent agonism by tetracyclic antidepressants mianserin, setiptiline, and mirtazapine suggests a mechanism for their clinically observed antimigraine properties. Using cryo-EM and mutagenesis studies, we uncover and characterize unique agonist-like binding poses of mianserin and setiptiline at 5-HT1eR distinct from similar drug scaffolds in inactive-state 5-HTR structures. Together with computational studies, our data suggest that these binding poses alongside receptor-specific allosteric coupling in 5-HT1eR and 5-HT1FR contribute to the agonist activity of these antidepressants.

Structural Insights into the Unexpected Agonism of Tetracyclic Antidepressants at Serotonin Receptors 5-HT1eR and 5-HT1FR

Serotonin (5-hydroxytryptamine, 5-HT) acts via 13 different receptors in humans. Of these receptor subtypes, all but 5-HT1eR have confirmed roles in native tissue and are validated drug targets. Despite 5-HT1eR's therapeutic potential and plausible druggability, the mechanisms of its activation remain elusive. To illuminate 5-HT1eR's pharmacology in relation to the highly homologous 5-HT1FR, we screened a library of aminergic receptor ligands at both receptors and observe 5-HT1e/1FR agonism by multicyclic drugs described as pan-antagonists at 5-HT receptors. Potent agonism by tetracyclic antidepressants mianserin, setiptiline, and mirtazapine suggests a mechanism for their clinically observed anti-migraine properties. Using cryoEM and mutagenesis studies, we uncover and characterize unique agonist-like binding poses of mianserin and setiptiline at 5-HT1eR distinct from similar drug scaffolds in inactive-state 5-HTR structures. Together with computational studies, our data suggest that these binding poses alongside receptor-specific allosteric coupling in 5-HT1eR and 5-HT1FR contribute to the agonist activity of these antidepressants.

Structural studies of serotonin receptor family

Serotonin receptors, also known as 5-HT receptors, belong to the G protein-coupled receptors (GPCRs) superfamily. They mediate the effects of serotonin, a neurotransmitter that plays a key role in a wide range of functions including mood regulation, cognition and appetite. The functions of serotonin are mediated by a family of 5-HT receptors including 12 GPCRs belonging to six major families: 5-HT1, 5-HT2, 5-HT4, 5-HT5, 5-HT6 and 5-HT7. Despite their distinct characteristics and functions, these receptors' subtypes share common structural features and signaling mechanisms. Understanding the structure, functions and pharmacology of the serotonin receptor family is essential for unraveling the complexities of serotonin signaling and developing targeted therapeutics for neuropsychiatric disorders. However, developing drugs that selectively target specific receptor subtypes is challenging due to the structural similarities in their orthosteric binding sites. This review focuses on the recent advancements in the structural studies of 5-HT receptors, highlighting the key structural features of each subtype and shedding light on their potential as targets for mental health and neurological disorders (such as depression, anxiety, schizophrenia, and migraine) drugs. [BMB Reports 2023; 56(10): 527-536].

Identification, characterization, and transcription of serotonin receptors in rainbow trout (Oncorhynchus mykiss) in response to bacterial infection and salinity changes

Serotonergic system is involved in the regulation of physiological functions and behavioral traits including cognition, memory, aggression, stress coping, appetite and immunomodulation. Serotonin exerts its functions via binding distinct serotonin receptors which are classified into 7 groups. Salmonid exhibits expanded functional gene copies due to salmonid-specific whole genome duplication. However, serotonin receptor (htr) repertoire is not fully identified in rainbow trout (Oncorhynchus mykiss). In this study, we identified 39 htr genes, including 14 htr1, 4 htr2, 4 htr2 like, 3 htr3, 4 htr4, 2 htr5, 2 htr6, and 6 htr7 subtypes. We investigated physiological functions of serotonin receptors in response to bacterial pathogens exposure and salinity changes. We showed htr1, htr2, htr4 and htr7 subtypes were associated with immunomodulation in response to Vibrio anguillarum or Aeromonas salmonicida infection. Saltwater (salinity of 15) transfer significantly altered htr1, htr2, htr4, and htr7 subtypes, suggesting trout Htr was associated with osmoregulation. We further showed residues interacted with inverse agonist (methiothepin) and serotonin analogue (5-Carboxamidotryptamine) were conserved between trout and human, suggesting exogenous ligands targeting human HTRs might have a role in aquaculture. This study showed duplicated trout Htrs might be physiologically neofunctionalized and potentially exhibit pleiotropic effects in regulating immunomodulation and osmoregulation.

Synergistic Interaction of 5-HT1B and 5-HT2B Receptors in Cytoplasmic Ca2+ Regulation in Human Umbilical Vein Endothelial Cells: Possible Involvement in Pathologies

The aim of this work was to explore the involvement of 5-HT1B and 5-HT2B receptors (5-HT1BR and 5-HT2BR) in the regulation of free cytoplasmic calcium concentration ([Ca2+]i) in human umbilical vein endothelial cells (HUVEC). We have shown by quantitative PCR analysis, that 5-HT1BR and 5-HT2BR mRNAs levels are almost equal in HUVEC. Immunofluorescent staining demonstrated, that 5-HT1BR and 5-HT2BR are expressed both in plasma membrane and inside the cells. Intracellular 5-HT1BR are localized mainly in the nuclear region, whereas 5-HT2BR receptors are almost evenly distributed in HUVEC. 5-HT, 5-HT1BR agonist CGS12066B, or 5-HT2BR agonist BW723C86 added to HUVEC caused a slight increase in [Ca2+]i, which was much lower than that of histamine, ATP, or SFLLRN, an agonist of protease-activated receptors (PAR1). However, activation of 5-HT1BR with CGS12066B followed by activation of 5-HT2BR with BW723C86 manifested a synergism of response, since several-fold higher rise in [Ca2+]i occurred. CGS12066B caused more than a 5-fold increase in [Ca2+]i rise in HUVEC in response to 5-HT. This 5-HT induced [Ca2+]i rise was abolished by 5-HT2BR antagonist RS127445, indicating that extracellular 5-HT acts through 5-HT2BR. Synergistic [Ca2+]i rise in response to activation of 5-HT1BR and 5-HT2BR persisted in a calcium-free medium. It was suppressed by the phospholipase C inhibitor U73122 and was not inhibited by the ryanodine and NAADP receptors antagonists dantrolene and NED-19. [Ca2+]i measurements in single cells demonstrated that activation of 5-HT2BR alone by BW723C86 caused single asynchronous [Ca2+]i oscillations in 19.8 ± 4.2% (n = 3) of HUVEC that occur with a long delay (66.1 ± 4.3 s, n = 71). On the contrary, histamine causes a simultaneous and almost immediate increase in [Ca2+]i in all the cells. Pre-activation of 5-HT1BR by CGS12066B led to a 3-4 fold increase in the number of HUVEC responding to BW723C86, to synchronization of their responses with a delay shortening, and to the bursts of [Ca2+]i oscillations in addition to single oscillations. In conclusion, to get a full rise of [Ca2+]i in HUVEC in response to 5-HT, simultaneous activation of 5-HT1BR and 5-HT2BR is required. 5-HT causes an increase in [Ca2+]i via 5-HT2BR while 5-HT1BR could be activated by the membrane-permeable agonist CGS12066B. We hypothesized that CGS12066B acts via intracellular 5-HT1BR inaccessible to extracellular 5-HT. Intracellular 5-HT1BR might be activated by 5-HT which could be accumulated in EC under certain pathological conditions.

Structure-based pharmacophore modeling 2. Developing a novel framework for structure-based pharmacophore model generation and selection

Pharmacophore models are three-dimensional arrangements of molecular features required for biological activity that are used in ligand identification efforts for many biological targets, including G protein-coupled receptors (GPCR). Though GPCR are integral membrane proteins of considerable interest as targets for drug development, many of these receptors lack known ligands or experimentally determined structures necessary for ligand- or structure-based pharmacophore model generation, respectively. Thus, we here present a structure-based pharmacophore modeling approach that uses fragments placed with Multiple Copy Simultaneous Search (MCSS) to generate high-performing pharmacophore models in the context of experimentally determined, as well as modeled GPCR structures. Moreover, we have addressed the oft-neglected topic of pharmacophore model selection via development of a cluster-then-predict machine learning workflow. Herein score-based pharmacophore models were generated in experimentally determined and modeled structures of 13 class A GPCR and resulted in pharmacophore models exhibiting high enrichment factors when used to search a database containing 569 class A GPCR ligands. In addition, classification of pharmacophore models with the best performing cluster-then-predict logistic regression classifier resulted in positive predictive values (PPV) of 0.88 and 0.76 for selecting high enrichment pharmacophore models from among those generated in experimentally determined and modeled structures, respectively.

Migraine drugs

According to recent studies, migraine affects more than 1 billion people worldwide, making it one of the world’s most prevalent diseases. Although this highly debilitating illness has been known since ancient times, the first therapeutic drugs to treat migraine, ergotamine (Gynergen) and dihydroergotamine (Dihydergot), did not appear on the market until 1921 and 1946, respectively. Both drugs originated from Sandoz, the world’s leading pharmaceutical company in ergot alkaloid research at the time. Historically, ergot alkaloids had been primarily used in obstetrics, but with methysergide (1-methyl-lysergic acid 1′-hydroxy-butyl-(2S)-amide), it became apparent that they also held some potential in migraine treatment. Methysergide was the first effective prophylactic drug developed specifically to prevent migraine attacks in 1959. On the basis of significantly improved knowledge of migraine pathophysiology and the discovery of serotonin and its receptors, Glaxo was able to launch sumatriptan in 1992. It was the first member from the class of triptans, which are selective 5-HT1B/1D receptor agonists. Recent innovations in acute and preventive migraine therapy include lasmiditan, a selective 5-HT1F receptor agonist from Eli Lilly, the gepants, which are calcitonin gene-related peptide (CGRP) receptor antagonists discovered at Merck & Co and BMS, and anti-CGRP/receptor monoclonal antibodies from Amgen, Pfizer, Eli Lilly, and others.

Graphical abstract

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.

Molecular Dynamics Refinement of Open State Serotonin 5-HT3A Receptor Structures

Pentameric ligand-gated ion channels play an important role in mediating fast neurotransmissions. As a member of this receptor family, cation-selective 5-HT₃ receptors are a clinical target for treating nausea and vomiting associated with chemotherapy and radiation therapy (Thompson and Lummis, 2006). Multiple cryo-electron microscopy (cryo-EM) structures of 5-HT₃ receptors have been determined in distinct functional states (e.g., open, closed, etc.) (Basak et al., 2018; Basak et al., 2018; Polovinkin et al., 2018; Zhang et al., 2015). However, recent work has shown that the transmembrane pores of the open 5-HT₃ receptor structures rapidly collapse and become artificially asymmetric in molecular dynamics (MD) simulations. To avoid this hydrophobic collapse, Dämgen and Biggin developed an equilibration protocol that led to a stable open state structure of the glycine receptor in MD simulations (Dämgen and Biggin, 2020). However, the protocol failed to yield open-like structures of the 5-HT₃ receptor in our simulations. Here, we present a refined equilibration protocol that involves the rearrangement of the transmembrane helices to achieve stable open state structures of the 5-HT₃ receptor that allow both water and ion permeation through the channel. Notably, channel gating is mediated through collective movement of the transmembrane helices, involving not only pore lining M2 helices but also their cross-talk with the adjacent M1 and M3 helices. Thus, the successful application of our refined equilibration protocol underscores the importance of the conformational coupling between the transmembrane helices in stabilizing open-like structures of the 5-HT₃ receptor.

The Bright Side of Psychedelics: Latest Advances and Challenges in Neuropharmacology

The need to identify effective therapies for the treatment of psychiatric disorders is a particularly important issue in modern societies. In addition, difficulties in finding new drugs have led pharmacologists to review and re-evaluate some past molecules, including psychedelics. For several years there has been growing interest among psychotherapists in psilocybin or lysergic acid diethylamide for the treatment of obsessive-compulsive disorder, of depression, or of post-traumatic stress disorder, although results are not always clear and definitive. In fact, the mechanisms of action of psychedelics are not yet fully understood and some molecular aspects have yet to be well defined. Thus, this review aims to summarize the ethnobotanical uses of the best-known psychedelic plants and the pharmacological mechanisms of the main active ingredients they contain. Furthermore, an up-to-date overview of structural and computational studies performed to evaluate the affinity and binding modes to biologically relevant receptors of ibogaine, mescaline, N,N-dimethyltryptamine, psilocin, and lysergic acid diethylamide is presented. Finally, the most recent clinical studies evaluating the efficacy of psychedelic molecules in some psychiatric disorders are discussed and compared with drugs already used in therapy.

Structural identification of lysophosphatidylcholines as activating ligands for orphan receptor GPR119

Lysophosphatidylcholine (LPC) is an essential mediator in human lipid metabolism and is associated with a variety of diseases, but the exact identity of LPC receptors remains controversial. Through extensive biochemical and structural analyses, we have identified the orphan receptor GPR119 as the receptor for LPC. The structure of the GPR119-G-protein complex without any added ligands reveals a density map that fits well with LPC, which is further confirmed by mass spectrometry and functional studies. As LPCs are abundant on the cell membrane, their preoccupancy in the receptor may lead to 'constitutive activity' of GPR119. The structure of GPR119 bound to APD668, a clinical drug candidate for type 2 diabetes, reveals an exceedingly similar binding mode to LPC. Together, these data highlight structural evidence for LPC function in regulating glucose-dependent insulin secretion through direct binding and activation of GPR119, and provide structural templates for drug design targeting GPR119.

Evolutionary association of receptor-wide amino acids with G protein-coupling selectivity in aminergic GPCRs

Evolutionary analyses of aminergic G protein–coupled receptors reveal receptor-wide potential determinants of G protein–coupling selectivity.

G protein-coupled receptors (GPCRs) induce signal transduction pathways through coupling to four main subtypes of G proteins (G_s, G_i, G_q, and G_12/13), selectively. However, G protein selective activation mechanisms and residual determinants in GPCRs have remained obscure. Herein, we performed extensive phylogenetic analysis and identified specifically conserved residues for the aminergic receptors having similar coupling profiles. By integrating our methodology of differential evolutionary conservation of G protein–specific amino acids with structural analyses, we identified specific activation networks for G_s, G_i1, G_o, and G_q. To validate that these networks could determine coupling selectivity we further analyzed G_s-specific activation network and its association with G_s selectivity. Through molecular dynamics simulations, we showed that previously uncharacterized Glycine at position 7x41 plays an important role in receptor activation and it may determine G_s coupling selectivity by facilitating a larger TM6 movement. Finally, we gathered our results into a comprehensive model of G protein selectivity called “sequential switches of activation” describing three main molecular switches controlling GPCR activation: ligand binding, G protein selective activation mechanisms, and G protein contact.

Structural insights into the ligand binding and Gi coupling of serotonin receptor 5-HT5A

5-hydroxytryptamine receptor 5A (5-HT_5A) belongs to the 5-HT receptor family and signals through the G_i/o protein. It is involved in nervous system regulation and an attractive target for the treatment of psychosis, depression, schizophrenia, and neuropathic pain. 5-HT_5A is the only G_i/o-coupled 5-HT receptor subtype lacking a high-resolution structure, which hampers the mechanistic understanding of ligand binding and G_i/o coupling for 5-HT_5A. Here we report a cryo-electron microscopy structure of the 5-HT_5A–G_i complex bound to 5-Carboxamidotryptamine (5-CT). Combined with functional analysis, this structure reveals the 5-CT recognition mechanism and identifies the receptor residue at 6.55 as a determinant of the 5-CT selectivity for G_i/o-coupled 5-HT receptors. In addition, 5-HT_5A shows an overall conserved G_i protein coupling mode compared with other G_i/o-coupled 5-HT receptors. These findings provide comprehensive insights into the ligand binding and G protein coupling of G_i/o-coupled 5-HT receptors and offer a template for the design of 5-HT_5A-selective drugs.

Spectroscopic Evaluation of the Potential Neurotoxic Effects of a New Candidate for Anti-Seizure Medication-TP-315 during Chronic Administration (In Vivo)

The aim of this study was to investigate the potential neurotoxic effect of the new anti-seizure medication candidate—5-(3-chlorophenyl)-4-hexyl-2,4-dihydro-3H-1,2,4-triazole-3-thione (TP-315), after chronic administration to mice. TP-315 was administered to mice intraperitoneally for 14 days. At 24 h post the last injection, animals were decapitated, their brains were acquired, flash-frozen in liquid nitrogen and cut into 10 µm slices. The FT-IR chemical imaging technique was used for the investigation of the potential neurotoxic effect in the cerebral cortex and hippocampus. The effect on the lipidomic and proteomic profile and on oxidative stress was investigated. The results showed no statistically significant changes in the above-mentioned parameters. TP-315 seems to pose no neurotoxic effect on the mouse brain after chronic use, therefore, its use should be safe.

Binding and Activation of Serotonergic G-Protein Coupled Receptors by the Multimodal Antidepressant Vortioxetine

G-protein coupled receptors (GPCRs) are important pharmacological targets. Despite substantial progress, important questions still remain concerning the details of activation: how can a ligand act as an agonist in one receptor but as an antagonist in a homologous receptor, and how can agonists activate a receptor despite lacking polar functional groups able to interact with helix 5 as is the case for the related adrenergic receptors? Studying vortioxetine (VXT), an important multimodal antidepressant drug, may elucidate both questions. Herein, we present a thorough in silico analysis of VXT binding to 5-HT_1A, 5-HT_1B, and 5-HT₇ receptors and compare it with available experimental data. We are able to rationalize the differential mode of action of VXT at different receptors, but also, in the case of the 5-HT_1A receptor, we observe the initial steps of activation that inform about an activation mechanism that does not involve polar interaction with helix 5. The results extend our current understanding of agonist and antagonist action at aminergic GPCRs.

In silico study of dimethyltryptamine analogues against 5-HT1B receptor: Molecular docking, dynamic simulations and ADMET prediction

Introduction: The 5-HT1B receptor has a potential role in various psychiatric disorders such as depression, anxiety, and post-traumatic stress disorder. The objective of this study was to perform docking and molecular dynamics simulation to evaluate at atomic level the behavior of N,N-dimethyltryptamine (DMT) on 5-HT1B receptor.Methods: In this study, initially, a search for DMT was performed using the PubChem database. Subsequently, molecular docking was executed using AutoDock Vina based in PyRx 0.8 with a 95% analogy. Additionally, ergotamine (ERG) and serotonin were used as control. Then, it ran a total of 100 ns molecular dynamics simulations on 5-HT1B bound with DMT, serotonin, 112814775, and ERG. Finally, pharmacokinetic prediction and IV acute toxicity for analogues and DMT were performed.Results: It was possible to show that 112814775 had the lowest binding energy with the receptor. In addition, 112814775 presented great conformational stability, low mobility, and stiffness compared to the control ligands: ERG, serotonin, and DMT subsequent dynamic analysis. With respect to the free energy calculation, contributions such as Van der Waals, electrostatics, and nonpolar interactions for all systems, were highlighted.Conclusion: 112814775 showed affinities with 5-HT1B receptor and evidenced notable behavior by molecular dynamic simulation according to root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), solvent-accessible surface area (SASA), the radius of gyration, number of hydrogen bond, and free energy calculated. These results established the possible relevance of in-silico studies in search of DMT analogues against the 5-HT1B receptor, which may be associated with alterations such as depression and anxiety, and may become future study molecules for the treatment of this type of disorder.

Ligand-Dependent Modulation of the Dynamics of Intracellular Loops Dictates Functional Selectivity of 5-HT2AR

The serotonin 2A receptor (5-HT_2AR) subtype of the G protein-coupled receptor (GPCR) family is involved in a plethora of neuromodulatory functions (e.g., neurogenesis, sleep, and cognitive processes). 5-HT_2AR is the target of pharmacologically distinct classes of ligands, binding of which either activate or inactivate the receptor. Although high-resolution structures of 5-HT_2AR as well as several other 5-HT GPCRs provided snapshots of both active and inactive conformational states, these structures, representing a truncated form of the receptor, cannot fully explain the mechanism of conformational transitions during their function. Importantly, biochemical studies have suggested the importance of intracellular loops in receptor functions. In our previous study, a model of the ligand-free form of 5-HT_2AR with the third intracellular loop (ICL3) has been meticulously built. Here, we have investigated the functional regulation of 5-HT_2AR with intact intracellular loops in ligand-free and five distinct ligand-bound configurations using unbiased atomistic molecular dynamics (MD) simulations. The selected ligands belong to either of the full, partial, or inverse agonist classes, which exert distinct pharmacological responses. We have observed significant structural, dynamic, and thermodynamic differences within ligand-bound complexes. Our results revealed, for the first time, that either activation or inactivation of the receptor upon specific ligand binding is primarily achieved through conformational transitions of its second and third intracellular loops (ICL2 and ICL3). A remarkable allosteric cross-talk between the ligand-binding site and the distal intracellular parts of the receptor, where binding of a specific ligand thermodynamically controls (either stabilizes or destabilizes) the intracellular region, consisting of crucial dynamic elements ICL2 and ICL3, and differential conformational transitions of these loops determine ligand-dependent functional selectivity.

Binding of SEP-363856 within TAAR1 and the 5HT1A receptor: implications for the design of novel antipsychotic drugs

Current medications for schizophrenia typically modulate dopaminergic neurotransmission. While affecting positive symptoms, antipsychotic drugs have little clinical effect on negative symptoms and cognitive impairment. Moreover, newer 'atypical' antipsychotic drugs also have significant metabolic adverse-effects. The recent positive clinical trial of the novel drug candidate SEP-363856, which targets non-dopamine receptors (trace amine-associated receptor and the 5HT_1A receptor), is a potentially promising development for the management of schizophrenia. In this perspective, we briefly overview the role of TAAR1 and the 5HT_1A receptor in schizophrenia and explore the specific binding characteristics of SEP-363856 at these receptors. Molecular dynamics simulations (MDS) indicate that SEP-363856 interacts with a small, common set of conserved residues within the TAAR1 and 5HT_1A ligand-binding domain. The primary interaction of SEP-363856 involves binding to the negatively charged aspartate residue (Asp103^3.32, TAAR1; Asp116^3.32, 5HT_1A). In general, the binding of SEP-363856 within TAAR1 involves a greater number of aromatic contacts compared to 5HT_1A. MDS provides important insights into the molecular basis of binding site interactions of SEP-363856 with TAAR1 and the 5HT_1A receptor, which will be beneficial for understanding the pharmacological uniqueness of SEP-363856 and for the design of novel drug candidates for these newly targeted receptors in the treatment of schizophrenia and related disorders.

Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

G protein-coupled receptors (GPCRs) constitute the largest family of membrane proteins in the human genome and are important therapeutic targets. During the last decade, the number of atomic-resolution structures of GPCRs has increased rapidly, providing insights into drug binding at the molecular level. These breakthroughs have created excitement regarding the potential of using structural information in ligand design and initiated a new era of rational drug discovery for GPCRs. The molecular docking method is now widely applied to model the three-dimensional structures of GPCR-ligand complexes and screen for chemical probes in large compound libraries. In this review article, we first summarize the current structural coverage of the GPCR superfamily and the understanding of receptor-ligand interactions at atomic resolution. We then present the general workflow of structure-based virtual screening and strategies to discover GPCR ligands in chemical libraries. We assess the state of the art of this research field by summarizing prospective applications of virtual screening based on experimental structures. Strategies to identify compounds with specific efficacy and selectivity profiles are discussed, illustrating the opportunities and limitations of the molecular docking method. Our overview shows that structure-based virtual screening can discover novel leads and will be essential in pursuing the next generation of GPCR drugs. SIGNIFICANCE STATEMENT: Extraordinary advances in the structural biology of G protein-coupled receptors have revealed the molecular details of ligand recognition by this large family of therapeutic targets, providing novel avenues for rational drug design. Structure-based docking is an efficient computational approach to identify novel chemical probes from large compound libraries, which has the potential to accelerate the development of drug candidates.

Vibrational spectroscopy analysis of ligand efficacy in human M2 muscarinic acetylcholine receptor (M2R)

The intrinsic efficacy of ligand binding to G protein-coupled receptors (GPCRs) reflects the ability of the ligand to differentially activate its receptor to cause a physiological effect. Here we use attenuated total reflection-Fourier transform infrared (ATR-FTIR) spectroscopy to examine the ligand-dependent conformational changes in the human M2 muscarinic acetylcholine receptor (M2R). We show that different ligands affect conformational alteration appearing at the C=O stretch of amide-I band in M2R. Notably, ATR-FTIR signals strongly correlated with G-protein activation levels in cells. Together, we propose that amide-I band serves as an infrared probe to distinguish the ligand efficacy in M2R and paves the path to rationally design ligands with varied efficacy towards the target GPCR.

Cholesterol footprint in high-resolution structures of serotonin receptors: Where are we now and what does it mean?

An emerging feature of several high-resolution GPCR structures is the presence of closely bound cholesterol molecules. In this Perspective, we share the excitement of the recent advancements in GPCR structural biology. We further highlight our laboratory's journey in comprehensively elucidating functional sensitivity of GPCRs (using the serotonin_1A receptor as a representative neurotransmitter GPCR) to membrane cholesterol and validation using a variety of assays and molecular dynamics simulations. Although high-resolution structures of many GPCRs have been reported in the last few years, the structure of the serotoin_1A receptor proved to be elusive for a long time. Very recently the cryo-EM structure of the serotoin_1A receptor displaying 10 bound cholesterol molecules has been reported. We conclude by providing a critical analysis of caveats involved in GPCR structure determination.

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-HT_2A 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.

Molecular Dynamics Simulations for the Determination of the Characteristic Structural Differences between Inactive and Active States of Wild Type and Mutants of the Orexin2 Receptor

The orexin2 receptor (OX2R), which is classified as a class A G protein-coupled receptor (GPCR), is the target of our study. We performed over 20 several-microsecond-scale molecular dynamics simulations of the wild type and mutants of OX2R to extract the characteristics of the structural changes taking place in the active state. We introduced mutations that exhibited the stable inactive state and the constitutively active state in class A GPCRs. In these simulations, significant characteristic structural changes were observed in the V309^6.40Y mutant, which corresponded to a constitutively active mutant. These conformational changes include the outward movement of the transmembrane helix 6 (TM6) and the inward movement of TM7, which are common structural changes in the activation of GPCRs. In addition, we extracted a suitable index for the quantitative evaluation of the active and inactive states of GPCRs, namely, the inter-atomic distance of Cα atoms between x(3.46) and Y(7.53). The structures of the inactive and active states solved by X-ray crystallography and cryo-electron microscopy can be classified using the inter-atomic distance. Furthermore, we clarified that the inward movement of TM7 requires the swapping of M305^6.36 on TM6 and L367^7.56 on TM7. Finally, we discussed the structural advantages of TM7 inward movement for GPCR activation.

Structural insights into the lipid and ligand regulation of serotonin receptors

Serotonin, or 5-hydroxytryptamine (5-HT), is an important neurotransmitter^1,2 that activates the largest subtype family of G-protein-coupled receptors³. Drugs that target 5-HT_1A, 5-HT_1D, 5-HT_1E and other 5-HT receptors are used to treat numerous disorders⁴. 5-HT receptors have high levels of basal activity and are subject to regulation by lipids, but the structural basis for the lipid regulation and basal activation of these receptors and the pan-agonism of 5-HT remains unclear. Here we report five structures of 5-HT receptor-G-protein complexes: 5-HT_1A in the apo state, bound to 5-HT or bound to the antipsychotic drug aripiprazole; 5-HT_1D bound to 5-HT; and 5-HT_1E in complex with a 5-HT_1E- and 5-HT_1F-selective agonist, BRL-54443. Notably, the phospholipid phosphatidylinositol 4-phosphate is present at the G-protein-5-HT_1A interface, and is able to increase 5-HT_1A-mediated G-protein activity. The receptor transmembrane domain is surrounded by cholesterol molecules-particularly in the case of 5-HT_1A, in which cholesterol molecules are directly involved in shaping the ligand-binding pocket that determines the specificity for aripiprazol. Within the ligand-binding pocket of apo-5-HT_1A are structured water molecules that mimic 5-HT to activate the receptor. Together, our results address a long-standing question of how lipids and water molecules regulate G-protein-coupled receptors, reveal how 5-HT acts as a pan-agonist, and identify the determinants of drug recognition in 5-HT receptors.

Identification of a Potent and Selective 5-HT1A Receptor Agonist with In Vitro and In Vivo Antinociceptive Activity

Opioids are the gold standard drugs for the treatment of acute and chronic severe pain, although their serious side effects constitute a big limitation. In the search for new and safer drugs, 5-HT_1AR agonists are emerging as potential candidates in pain relief therapy. In this work, we evaluated the affinity and activity of enantiomers of the two newly synthesized, potent 5-HT_1AR agonists N-[(2,2-diphenyl-1,3-dioxolan-4-yl)methyl]-2-[2-(pyridin-4-yl)phenoxy]ethan-1-ammonium hydrogenoxalate (rac-1) and N-((2,2-diphenyl-1,3-dioxolan-4-yl)methyl)-2-(2-(1-methyl-1H-imidazol-5-yl)phenoxy)ethan-1-ammonium hydrogenoxalate (rac-2) in vitro and in vivo. The role of chirality in the interaction with 5-HT_1AR was evaluated by molecular docking. The activity of the rac-1 was tested in mouse models of acute pain (hot plate) and severe tonic nociceptive stimulation (intraplantar formalin test). Rac-1 was active in the formalin test with a reduction in paw licking in both phases at 10 mg/kg, and its effect was abolished by the selective 5-HT_1AR antagonist, WAY-100635. The eutomer (S)-1, but not the racemate, was active during the hot plate test at 10 and 20 mg/kg, and this effect was abolished by 30 min treatment with WAY-100635 at 30 min. Similarly to 8-OH-DPAT, (S)-1 evoked a slow outward current and depressed spontaneous glutamatergic transmission in superficial dorsal horn neurons, more effectively than rac-1. The eutomer (S)-1 showed promising developability properties, such as high selectivity over 5-HT subtypes, no interaction with the μ receptors, and low hepato- and cardiotoxicity. Therefore, (S)-1 may represent a potential candidate for the treatment of acute and chronic pain without having the adverse effects that are commonly associated with the classic opioid drugs.

Structure, dynamics and lipid interactions of serotonin receptors: excitements and challenges

Serotonin (5-hydroxytryptamine, 5-HT) is an intrinsically fluorescent neurotransmitter found in organisms spanning a wide evolutionary range. Serotonin exerts its diverse actions by binding to distinct cell membrane receptors which are classified into many groups. Serotonin receptors are involved in regulating a diverse array of physiological signaling pathways and belong to the family of either G protein-coupled receptors (GPCRs) or ligand-gated ion channels. Serotonergic signaling appears to play a key role in the generation and modulation of various cognitive and behavioral functions such as sleep, mood, pain, anxiety, depression, aggression, and learning. Serotonin receptors act as drug targets for a number of diseases, particularly neuropsychiatric disorders. The signaling mechanism and efficiency of serotonin receptors depend on their amazing ability to rapidly access multiple conformational states. This conformational plasticity, necessary for the wide variety of functions displayed by serotonin receptors, is regulated by binding to various ligands. In this review, we provide a succinct overview of recent developments in generating and analyzing high-resolution structures of serotonin receptors obtained using crystallography and cryo-electron microscopy. Capturing structures of distinct conformational states is crucial for understanding the mechanism of action of these receptors, which could provide important insight for rational drug design targeting serotonin receptors. We further provide emerging information and insight from studies on interactions of membrane lipids (such as cholesterol) with serotonin receptors. We envision that a judicious combination of analysis of high-resolution structures and receptor-lipid interaction would allow a comprehensive understanding of GPCR structure, function and dynamics, thereby leading to efficient drug discovery.

In silico decryption of serotonin-receptor binding: local non-covalent interactions and long-range conformational changes

Serotonin-receptor binding is the key step in the process behind serotonin functionality, including several psychological and physiological behaviours. This study is focused on identifying the main non-covalent interactions controlling the stability of serotonin-receptor complexes as well as the main conformational changes in the receptor due to serotonin-receptor binding using classical molecular dynamics simulations and quantum chemical calculations. A qualitative analysis based on two order parameters ((i) the centre of mass distance and (ii) the angle between the surface normals of each aromatic residue and serotonin in the binding site) on the serotonin-receptor complex trajectory suggests that the T-type stacking interaction is predominant in the binding site. Quantum chemical calculations of the stacking interaction energy provide the quantitative contributions of important aromatic residues to the stabilization of the complex. Furthermore, a three body stacking interaction (named 'L'-type) was observed and likely contributes to the stability of the complex. Direct and water-mediated hydrogen bonding between the residues in the binding site and serotonin contributes to the complex stability. Principal component analysis of the molecular dynamics simulation trajectory of the serotonin-receptor complex and the apo-receptor in water indicates that the whole receptor is significantly stabilized due to serotonin binding. An analysis based on the dynamic cross correlation function reflects the strong correlation between trans-membrane (TM)3, TM5, TM6 (containing residues responsible for the stacking interaction and hydrogen bonding) and mini-G0 which may participate in signal transduction leading to the functionality of serotonin.

Anxiolytic activity of Psidium guajava in mice subjected to chronic restraint stress and effect on neurotransmitters in brain

The anxiolytic activity of Psidium guajava L. leaf ethanolic extract (PLE) and its effect on neurotransmitter systems was investigated. PLE, extracted using Soxhlet apparatus, was subjected to preliminary qualitative and quantitative (flavonoids and phenols) analyses. The anxiolytic activity at 100, 200, and 400 mg/Kg doses were assessed in mice using elevated plus maze (EPM) and light/dark transition (LDT) test models on days 1 and 16. Neurotransmitters such as monoamines (serotonin, norepinephrine, and dopamine), γ-aminobutyric acid (GABA), and glutamate were estimated in different regions of the brain (cortex, hippocampus, and cerebellum and brain stem). Phytoconstituents identified using gas chromatography-mass spectrometry and liquid chromatography-mass spectrometry were analyzed in silico to evaluate their potential binding mode(s) to GABA_A and 5-HT_1A receptors. Phytochemical studies showed the presence of alkaloids, tannins, flavonoids, saponins, anthraquinone glycosides, carbohydrates, and proteins, whereas total flavonoid and phenol contents were estimated to be 64.96 ± 0.95 and 206.58 ± 1.60 mg/g of dried extract, respectively. PLE treatment significantly enhanced exploratory activity of mice in EPM and LDT models with significant effects on monoamines, GABA and glutamate levels in the brain. The in silico studies suggested the interaction(s) of PLE component(s) with GABA_A /5-HT_1A receptors as a potential mechanism of its anxiolytic activity.

Molecular Mechanism for Ligand Recognition and Subtype Selectivity of α2C Adrenergic Receptor

Adrenergic G-protein-coupled receptors (GPCRs) mediate different cellular signaling pathways in the presence of endogenous catecholamines and play important roles in both physiological and pathological conditions. Extensive studies have been carried out to investigate the structure and function of β adrenergic receptors (βARs). However, the structure of α adrenergic receptors (αARs) remains to be determined. Here, we report the structure of the human α_2C adrenergic receptor (α_2CAR) with the non-selective antagonist, RS79948, at 2.8 Å. Our structure, mutations, modeling, and functional experiments indicate that a α_2CAR-specific D206^ECL2-R409^ECL3-Y405^6.58 network plays a role in determining α₂ adrenergic subtype selectivity. Furthermore, our results show that a specific loosened helix at the top of TM4 in α_2CAR is involved in receptor activation. Together, our structure of human α_2CAR-RS79948 provides key insight into the mechanism underlying the α₂ adrenergic receptor activation and subtype selectivity.

Neuropharmacological assessment in mice and molecular docking of piperazine derivative LQFM212

Piperazine derivatives are an attractive class of chemical compounds for the treatment of various mental illness. Herein, we demonstrated the synthesis of LQFM212, a piperazine derivative, behavioral evaluation in mice and computational studies. In neuropharmacological assessment, LQFM212 treatment at doses of 18, 54 or 162 μmol/kg increased the sleep duration in sodium pentobarbital-induced sleep test. LQFM212 at dose of 162 μmol/kg increased climbing time in the chimney test and decreased the number of squares crossed in the open field test, suggesting that LQFM212 in high doses reduces spontaneous movement. However, LQFM212 treatment at the doses of 18 or 54 μmol/kg increased the preference for the center of field which could be indicative of anxiolytic-like effects. In elevated plus maze and light-dark box tests, LQFM212 treatment altered all parameters observed that demonstrate anxiolytic-like activity. These effects were reversed by flumazenil, mecamylamine, WAY-100635 and PCPA, but not with ketanserin, showing that anxiolytic-like activity involve benzodiazepine site of GABA_A receptor, nicotinic and serotonergic pathways. Molecular docking of LQFM212 showed that the ligand has more interactions with GABA_A receptor than with 5-HT_1A receptor. Despite the involvement of benzodiazepine site on anxiolytic-like effect of LQFM212, treatment with this compound did not alter cognitive function in the step-down avoidance test. In this sense, this piperazine derivative is a good prototype for treating anxiety disorders with putative mechanism of action.

Anxiolytic-like effect of chalcone N-{4'[(2E)-3-(3-nitrophenyl)-1-(phenyl)prop-2-en-1-one]} acetamide on adult zebrafish (Danio rerio): Involvement of the 5-HT system

The action of anxiolytic compounds that act on selective serotonin receptors (SSRIs) have been scarcely evaluated. Serotonergic drugs have been shown to be effective in treating anxiety without presenting adverse effects as benzodiazepines. However, the anxiolytic effects take days to occur. This study aimed to evaluate the anxiolytic effect of the synthetic chalcone, 4'-[(2E) -3- (3-nitrophenyl) -1- (phenyl) prop-2-en-1-one] acetamide (PAAMNBA), and its possible mechanism of action in adult zebrafish (Danio rerio). PAAMNBA was synthesized with a yield of 51.3% and its chemical structure was determined by ¹H and ¹³C NMR. Initially, PAAPMNBA was intraperitoneally administered to zebrafish (n = 6/group) at doses of 4, 12, or 40 mg/kg, and the animals were subsequently subjected to acute and open field toxicity tests. PAAMNBA was administered to the other groups (n = 6/group) for analyzing its effect in the light and dark test. The involvement of the serotonergic (5HT) system was also evaluated using 5-HTR 1, 5-HTR 2A/2C, and 5-HTR 3A/3B receptor antagonists, namely, pizotifeo, granizetron, and ciproeptadina, respectively. Molecular coupling was performed using the 5-HT1 receptor. PAAMNBA was found to be non-toxic, reduced the locomotor activity, and had an anxiolytic effect in adult zebrafish. The effect was reduced by pretreatment with pizotifene and was not reversed by treatment with granizetron and cyproeptadine. A previous in vivo molecular coupling study indicated that chalcones interact with the 5-HT1 receptor. The results suggested that the chalcone, PAAPMNBA, has anxiolytic activity, that is mediated by the serotonergic system via the 5-HT1 receptor. The interaction of PAAPMNBA with the 5-HT1 receptor was confirmed by molecular docking studies.

Molecular mechanisms of fentanyl mediated β-arrestin biased signaling

The development of novel analgesics with improved safety profiles to combat the opioid epidemic represents a central question to G protein coupled receptor structural biology and pharmacology: What chemical features dictate G protein or β-arrestin signaling? Here we use adaptively biased molecular dynamics simulations to determine how fentanyl, a potent β-arrestin biased agonist, binds the μ-opioid receptor (μOR). The resulting fentanyl-bound pose provides rational insight into a wealth of historical structure-activity-relationship on its chemical scaffold. Following an in-silico derived hypothesis we found that fentanyl and the synthetic opioid peptide DAMGO require M153 to induce β-arrestin coupling, while M153 was dispensable for G protein coupling. We propose and validate an activation mechanism where the n-aniline ring of fentanyl mediates μOR β-arrestin through a novel M153 “microswitch” by synthesizing fentanyl-based derivatives that exhibit complete, clinically desirable, G protein biased coupling. Together, these results provide molecular insight into fentanyl mediated β-arrestin biased signaling and a rational framework for further optimization of fentanyl-based analgesics with improved safety profiles.

The global opioid crisis has drawn significant attention to the risks associated with over-use of synthetic opioids. Despite the public attention, and perhaps in-line with the profit-based incentives of the pharmaceutical industry, there is no public structure of mu-opioid receptor bound to fentanyl or fentanyl derivatives. A publicly available structure of the complex would allow open-source development of safer painkillers and synthetic antagonists. Current overdose antidotes, antagonists, require natural products in their synthesis which persists a sizable barrier to market and develop better antidotes. In this work we use advance molecular dynamics techniques to obtain the bound geometry of mu-opioid receptor with fentanyl (and derivatives) and uncovered a novel molecular switch involved in receptor activation. Based on our in-silico structure, we synthesized and tested novel compounds to validate our predicted structure. Herein we report the bound state of several dangerous fentanyl derivatives and introduce new derivatives with signaling profiles that may lead to lower risk of respiratory depression.

Binding Interactions of Ergotamine and Dihydroergotamine to 5-Hydroxytryptamine Receptor 1B (5-HT1b) Using Molecular Dynamics Simulations and Dynamic Network Analysis

Ergotamine (ERG) and dihydroergotamine (DHE), common migraine drugs, have small structural differences but lead to clinically important distinctions in their pharmacological profiles. For example, DHE is less potent than ERG by about 10-fold at the 5-hydroxytrptamine receptor 1B (5-HT_1B). Although the high-resolution crystal structures of the 5-HT_1B receptor with both ligands have been solved, the high similarity between these two complex structures does not sufficiently explain their activity differences and the activation mechanism of the receptor. Hence, an examination of the dynamic motion of both drugs with the receptor is required. In this study, we ran a total of 6.0 μs molecular dynamics simulations on each system. Our simulation data show the subtle variations between the two systems in terms of the ligand-receptor interactions and receptor secondary structures. More importantly, the ligand and protein root-mean-square fluctuations (RMSFs) for the two systems were distinct, with ERG having a trend of lower RMSF values, indicating it to be bound tighter to 5-HT_1B with less fluctuations. The molecular mechanism-general born surface area (MM-GBSA) binding energies illustrate this further, proving ERG has an overall stronger MM-GBSA binding energy. Analysis of several different microswitches has shown that the 5-HT_1B-ERG complex is in a more active conformation state than 5-HT_1B-DHE, which is further supported by the dynamic network model, with reference to mutagenesis data with the critical nodes and the first three low-energy modes from the normal mode analysis. We also identify Trp327^6.48 and Phe331^6.52 as key residues involved in the active state 5-HT_1B for both ligands. Using the detailed dynamic information from our analysis, we made predictions for possible modifications to DHE and ERG that yielded five derivatives that might have more favorable binding energies and reduced structural fluctuations.

Potential for imaging the high-affinity state of the 5-HT1B receptor: a comparison of three PET radioligands with differing intrinsic activity

Background

Over the last decade, a few radioligands have been developed for PET imaging of brain 5-HT_1B receptors. The 5-HT_1B receptor is a G-protein-coupled receptor (GPCR) that exists in two different agonist affinity states. An agonist ligand is expected to be more sensitive towards competition from another agonist, such as endogenous 5-HT, than an antagonist ligand. It is of interest to know whether the intrinsic activity of a PET radioligand for the 5-HT_1B receptor impacts on its ability to detect changes in endogenous synaptic 5-HT density. Three high-affinity ¹¹C-labeled 5-HT_1B PET radioligands with differing intrinsic activity were applied to PET measurements in cynomolgus monkey to evaluate their sensitivity to be displaced within the brain by endogenous 5-HT. For these experiments, fenfluramine was pre-administered at two different doses (1.0 and 5.0 mg/kg, i.v.) to induce synaptic 5-HT release.

Results

A dose-dependent response to fenfluramine was detected for all three radioligands. At the highest dose of fenfluramine (5.0 mg/kg, i.v.), reductions in specific binding in the occipital cortex increased with radioligand agonist efficacy, reaching 61% for [¹¹C]3. The most antagonistic radioligand showed the lowest reduction in specific binding.

Conclusions

Three 5-HT_1B PET radioligands were identified with differing intrinsic activity that could be used in imaging high- and low-affinity states of 5-HT_1B receptors using PET. From this limited study, radioligand sensitivity to endogenous 5-HT appears to depend on agonist efficacy. More extensive studies are required to substantiate this suggestion.

Down-regulation effects of IFN-α on p11, 5-htr1b and 5-HTR4 protein levels were affected by NH4CL or MG132 treatment in SH-sy5y cells

In previous studies, we found interferon-α (IFN-α) could reduce protein levels of p11, 5-hydroxytryptamine receptor 1b (5-HT1b) and 5-hydroxytryptamine receptor 4 (5-HT4), but does not influence their messenger RNA levels in SH-sy5y cells. Thus, we investigated the post-transcriptional modulation of these molecules by IFN-α. SH-sy5y cells were treated with IFN-α, NH₄Cl or MG132 alone or in combination, and then the protein levels of p11, 5-HT1b and 5-HT4 were analyzed by western blots. The regulatory effects of p11 on 5-HT1b and 5-HT4 were also determined in p11 knock-down cells. NH4Cl but not MG132 could reverse the protein level of p11 in IFN-α-treated SH-sy5y cells. MG132 could recover the protein levels of 5-HT1b and 5-HT4 in p11 knock-down cells. The down-regulation effects of IFN-α on p11, 5-HT1b and 5-HT4 were associated with the lysosome and ubiquitin-proteasome-mediated pathways. p11 was identified as a potent regulator to modulate the ubiquitination of 5-HT1b and 5-HT4. Therefore, it could be potential target therapies in IFN-ainduced depression.

Functional Characterization of the Prejunctional Receptors Mediating the Inhibition by Ergotamine of the Rat Perivascular Sensory Peptidergic Drive

Calcitonin gene-related peptide (α-CGRP) released from perivascular sensory nerves induces decreases in diastolic blood pressure (DBP). Experimentally, this can be shown by spinal thoracic (T₉-T₁₂) electrical stimulation of these afferent fibers. Because ergotamine inhibits these neurogenic vascular responses and displays affinity for monoaminergic receptors that inhibit neurotransmitter release, we investigated whether this ergotamine-induced inhibition results from activation of serotonin 5-HT_1B/1D, dopamine D₂-like, and α₂-adrenergic receptors. Wistar rats were pithed and, under autonomic ganglion blockade, received intravenous infusions of methoxamine followed by ergotamine (0.1-3.1 μg kg^-1 min^-1). Thoracic T₉-T₁₂ electrical stimulation or an intravenous bolus of α-CGRP resulted in decreases in DBP. Ergotamine inhibited the electrically induced, but not α-CGRP-induced, responses. The vasodilator sensory inhibition by 3.1 μg of ergotamine kg^-1 min^-1 was resistant to simultaneous blockade of 5-HT_1B/1D, D₂-like, and α₂-adrenergic receptors upon addition of antagonists GR127935, haloperidol, and rauwolscine. Moreover, the inhibition by 0.31 μg of ergotamine kg^-1 min^-1 was unaltered by GR127935 and haloperidol, partly blocked by GR127935 and rauwolscine or rauwolscine and haloperidol, and abolished by GR127935, haloperidol, and rauwolscine. These findings imply that prejunctional 5-HT_1B/1D, D₂-like, and α₂-adrenergic receptors mediate the sensory inhibition induced by 0.31 μg of ergotamine kg^-1 min^-1, whereas larger doses may involve other receptors. Thus, ergotamine's ability to inhibit the perivascular sensory peptidergic drive may result in facilitation of its systemic vasoconstrictor properties.

Structural basis of ligand recognition at the human MT1 melatonin receptor

Melatonin (N-acetyl-5-methoxytryptamine) is a neurohormone that maintains circadian rhythms¹ by synchronization to environmental cues and is involved in diverse physiological processes² such as the regulation of blood pressure and core body temperature, oncogenesis, and immune function³. Melatonin is formed in the pineal gland in a light-regulated manner⁴ by enzymatic conversion from 5-hydroxytryptamine (5-HT or serotonin), and modulates sleep and wakefulness⁵ by activating two high-affinity G-protein-coupled receptors, type 1A (MT₁) and type 1B (MT₂)^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 burden⁷. Over-the-counter melatonin is widely used to alleviate jet lag and as a safer alternative to benzodiazepines and other sleeping aids^8,9, and is one of the most popular supplements in the United States¹⁰. Here, we present high-resolution room-temperature X-ray free electron laser (XFEL) structures of MT₁ in complex with four agonists: the insomnia drug ramelteon¹¹, two melatonin analogues, and the mixed melatonin-serotonin antidepressant agomelatine^12,13. The structure of MT₂ is described in an accompanying paper¹⁴. Although the MT₁ 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 MT₁, 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 MT₁ 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.

Structures of the 5-HT2A receptor in complex with the antipsychotics risperidone and zotepine

Many drugs target the serotonin 2A receptor (5-HT_2AR), including second-generation antipsychotics that also target the dopamine D₂ receptor (D₂R). These drugs often produce severe side effects due to non-selective binding to other aminergic receptors. Here, we report the structures of human 5-HT_2AR in complex with the second-generation antipsychotics risperidone and zotepine. These antipsychotics effectively stabilize the inactive conformation by forming direct contacts with the residues at the bottom of the ligand-binding pocket, the movements of which are important for receptor activation. 5-HT_2AR is structurally similar to 5-HT_2CR but possesses a unique side-extended cavity near the orthosteric binding site. A docking study and mutagenic studies suggest that a highly 5-HT_2AR-selective antagonist binds the side-extended cavity. The conformation of the ligand-binding pocket in 5-HT_2AR significantly differs around extracellular loops 1 and 2 from that in D₂R. These findings are beneficial for the rational design of safer antipsychotics and 5-HT_2AR-selective drugs.

Structure-activity relationship of serotonin derived tocopherol lipids

Tocopherol-based lipids are widely used for nucleic acid delivery. Using tocopherol molecules, we designed and synthesized 5-HT functionalized lipids by tethering 5-hydroxytryptamine (5-HT), a small molecule ligand as the head group to a natural amphiphilic molecule namely α-tocopherol (Vitamin E). This is with the aim of delivering nucleic acids specifically into cells expressing the serotonin receptors (5-hydroxytryptamine[5-HT]) which are abundant in the central nervous system. In order to achieve target recognition, we adopted an approach wherein two structurally different lipid molecules having serotonin as the head group was conjugated to tocopherol via different linkers thus generating lipids with either free -NH₂ or -OH moiety. The corresponding lipids designated as Lipid A (Tocopheryl carbonate serotonin-NH₂) and Lipid B (Tocopheryl 2-hydroxy propyl ammonium serotonin-OH), were formulated with co-lipids 1,2-dioleoyl-sn-glycero-3-phosphatidyl-ethanolamine (DOPE) and 1,2-dioleoyl-sn-glycero-sn-3-phosphatidylcholine (DOPC) and evaluated for their ability to deliver plasmid DNA through reporter gene expression assays in vitro. Furthermore, the physicochemical characteristics and cellular interactions of the formulations were examined using serotonin-receptor enriched cells in order to distinguish the structural and functional attributes of both lipids. Cell-based gene expression studies reveal that in comparison to Lipid A, a formulation of Lipid B prepared with DOPE as the co-lipid, contributes to efficient uptake leading to significant enhancement in transfection. Specific interactions explored by molecular docking studies suggests the role of the hydroxyl moiety and the enantiospecific significance of serotonin- conjugated tocopherol lipids in recognizing these receptors thus signifying a promising lipid-based approach to target the serotonin receptors in the central nervous system.

Three-dimensional descriptors for aminergic GPCRs: dependence on docking conformation and crystal structure

Three-dimensional descriptors are often used to search for new biologically active compounds, in both ligand- and structure-based approaches, capturing the spatial orientation of molecules. They frequently constitute an input for machine learning-based predictions of compound activity or quantitative structure-activity relationship modeling; however, the distribution of their values and the accuracy of depicting compound orientations might have an impact on the power of the obtained predictive models. In this study, we analyzed the distribution of three-dimensional descriptors calculated for docking poses of active and inactive compounds for all aminergic G protein-coupled receptors with available crystal structures, focusing on the variation in conformations for different receptors and crystals. We demonstrated that the consistency in compound orientation in the binding site is rather not correlated with the affinity itself, but is more influenced by other factors, such as the number of rotatable bonds and crystal structure used for docking studies. The visualizations of the descriptors distributions were prepared and made available online at http://chem.gmum.net/vischem_stability , which enables the investigation of chemical structures referring to particular data points depicted in the figures. Moreover, the performed analysis can assist in choosing crystal structure for docking studies, helping in selection of conditions providing the best discrimination between active and inactive compounds in machine learning-based experiments.

Aminergic GPCR-Ligand Interactions: A Chemical and Structural Map of Receptor Mutation Data

The aminergic family of G protein-coupled receptors (GPCRs) plays an important role in various diseases and represents a major drug discovery target class. Structure determination of all major aminergic subfamilies has enabled structure-based ligand design for these receptors. Site-directed mutagenesis data provides an invaluable complementary source of information for elucidating the structural determinants of binding of different ligand chemotypes. The current study provides a comparative analysis of 6692 mutation data points on 34 aminergic GPCR subtypes, covering the chemical space of 540 unique ligands from mutagenesis experiments and information from experimentally determined structures of 52 distinct aminergic receptor-ligand complexes. The integrated analysis enables detailed investigation of structural receptor-ligand interactions and assessment of the transferability of combined binding mode and mutation data across ligand chemotypes and receptor subtypes. An overview is provided of the possibilities and limitations of using mutation data to guide the design of novel aminergic receptor ligands.

Oncotoxic Properties of Serotonin Transporter Inhibitors and 5-HT1A Receptor Ligands

The cytotoxic activity of several serotonin transporter (SERT) inhibitors and subtype of serotonin receptor 1A (5-HT_1A receptor) ligands have been examined in androgen-insensitive human PC-3 prostate and neuroblastoma SH-SY5Y cancer cells. Almost all of the studied compounds (except 5-HT_1A receptor agonist (2R)-(+)-8-Hydroxy-2-(di-n-propylamino)tetralin hydrobromide (8-OH-DPAT)) exhibited absolute cytotoxic activity against the examined cancer cells. The compound 4-Fluoro-N-[2-[4-(7-methoxy-1-naphthalenyl)-1-piperazinyl]ethyl]benzamide hydrochloride (S14506) that showed highest activity against neuroblastoma tumors was the 5-HT_1A receptor agonist (although not alike other 5-HT_1A receptor agonists). On the other hand, the compound 6-nitro-2-(4-undecylpiperazin-1-yl)quinoline hydrochloride (AZ07) that had the highest activity against PC-3 prostate cancer cells was a compound exhibiting antagonistic activity against the 5-HT_1A receptor. Thus, compounds of oncotoxic properties S14506 and AZ07 should be evaluated further for their potential use in the prevention and treatment of cancer. Most of the 15 compounds tested exhibited either agonistic or antagonistic activity for both the cyclic adenosine monophosphate (cAMP) and extracellular signal-regulated kinase 1 and 2 (ERK1/2) pathways in human embryonic kidney 293 (HEK293) cells that overexpress the 5HT1AR gene. However, compounds paroxetine, N-Ac-paroxetine and 2-[4-(cyclobutylmethyl)piperazin-1-yl]-6-nitroquinoline hydrochloride (AB22) simultaneously exhibited antagonistic activity on the cAMP pathway and agonistic activity on the ERK1/2 pathway. Fluoxetine relative to compound AZ07 had almost three times lower cytotoxic activity against PC-3 prostate cancer cells. However, the proapoptotic activity of fluoxetine compared to compound AZ07 is almost two times higher which would suggest that the cytotoxic activity of both compounds may be dependent on different cell death mechanisms. Compound S14506 was found to be an antagonist of the serine-threonine protein kinase B (Akt) pathway. Prosurvival Akt activity may be reversed by Akt antagonists. Therefore, the antagonistic activity of S14506 on the Akt pathway may evoke caspase-3 expression and cytotoxicity. It appears that one should not expect a straightforward relationship between the activation of particular serotonergic pathways by selective serotonin reuptake inhibitors (SSRIs) and 5-HT_1A receptor ligands and their cytotoxic or cytoprotective activity. Additionally, nuclear transcription factor κB (NF-κB), which may be involved in 5-HT-dependent biochemical pathways by coordinating different subunits in the formation of a dimer, may regulate the transcription of different transduction pathways. Therefore, it can be suggested that the mechanism of the cytotoxic activity of certain compounds (serotonergic against nonserotonergic) may depend on the compound and cancer type being examined. Docking studies showed that S14506, buspirone and spiperone bind in similar ways in the 5-HT_1A receptor model and interacted with similar 5-HT_1A receptor residues. S14506 and spiperone were found to be located closer to both phenylalanines in TM6 than buspirone, thus exhibiting more antagonist binding modes.

Structural determinants of 5-HT2B receptor activation and biased agonism

Serotonin (5-hydroxytryptamine; 5-HT) receptors modulate a variety of physiological processes ranging from perception, cognition and emotion to vascular and smooth muscle contraction, platelet aggregation, gastrointestinal function and reproduction. Drugs that interact with 5-HT receptors effectively treat diseases as diverse as migraine headaches, depression and obesity. Here we present four structures of a prototypical serotonin receptor-the human 5-HT_2B receptor-in complex with chemically and pharmacologically diverse drugs, including methysergide, methylergonovine, lisuride and LY266097. A detailed analysis of these structures complemented by comprehensive interrogation of signaling illuminated key structural determinants essential for activation. Additional structure-guided mutagenesis experiments revealed binding pocket residues that were essential for agonist-mediated biased signaling and β-arrestin2 translocation. Given the importance of 5-HT receptors for a large number of therapeutic indications, insights derived from these studies should accelerate the design of safer and more effective medications.