Free energy calculations of the functional selectivity of 5-HT2B G protein-coupled receptor.

G Protein-Coupled Receptors (GPCRs) mediate intracellular signaling in response to extracellular ligand binding and are the target of one-third of approved drugs. Ligand binding modulates the GPCR molecular free energy landscape by preferentially stabilizing active or inactive conformations that dictate intracellular protein recruitment and downstream signaling. We perform enhanced sampling molecular dynamics simulations to recover the free energy surfaces of a thermostable mutant of the GPCR serotonin receptor 5-HT2B in the unliganded form and bound to a lysergic acid diethylamide (LSD) agonist and lisuride antagonist. LSD binding imparts a ∼110 kJ/mol driving force for conformational rearrangement into an active state. The lisuride-bound form is structurally similar to the apo form and only ∼24 kJ/mol more stable. This work quantifies ligand-induced conformational specificity and functional selectivity of 5-HT2B and presents a platform for high-throughput virtual screening of ligands and rational engineering of the ligand-bound molecular free energy landscape.

Biased Ligands Differentially Shape the Conformation of the Extracellular Loop Region in 5-HT2B Receptors.

G protein-coupled receptors are linked to various intracellular transducers, each pathway associated with different physiological effects. Biased ligands, capable of activating one pathway over another, are gaining attention for their therapeutic potential, as they could selectively activate beneficial pathways whilst avoiding those responsible for adverse effects. We performed molecular dynamics simulations with known ß-arrestin-biased ligands like lysergic acid diethylamide and ergotamine in complex with the 5-HT2B receptor and discovered that the extent of ligand bias is directly connected with the degree of closure of the extracellular loop region. Given a loose allosteric coupling of extracellular and intracellular receptor regions, we delineate a concept for biased signaling at serotonin receptors, by which conformational interference with binding pocket closure restricts the signaling repertoire of the receptor. Molecular docking studies of biased ligands gathered from the BiasDB demonstrate that larger ligands only show plausible docking poses in the ergotamine-bound structure, highlighting the conformational constraints associated with bias. This emphasizes the importance of selecting the appropriate receptor conformation on which to base virtual screening workflows in structure-based drug design of biased ligands. As this mechanism of ligand bias has also been observed for muscarinic receptors, our studies provide a general mechanism of signaling bias transferable between aminergic receptors.

Influence of Cholesterol and Its Stereoisomers on Members of the Serotonin Receptor Family.

Despite the ubiquity of cholesterol within the cell membrane, the mechanism by which it influences embedded proteins remains elusive. Numerous G-protein coupled receptors exhibit dramatic responses to membrane cholesterol with regard to the ligand-binding affinity and functional properties, including the 5-HT receptor family. Here, we use over 25 µs of unbiased atomistic molecular dynamics simulations to identify cholesterol interaction sites in the 5-HT1B and 5-HT2B receptors and evaluate their impact on receptor structure. Susceptibility to membrane cholesterol is shown to be subtype dependent and determined by the quality of interactions between the extracellular loops. Charged residues are essential for maintaining the arrangement of the extracellular surface in 5-HT2B; in the absence of such interactions, the extracellular surface of the 5-HT1B is malleable, populating a number of distinct conformations. Elevated cholesterol density near transmembrane helix 4 is considered to be conducive to the conformation of extracellular loop 2. Occupation of this site is also shown to be stereospecific, illustrated by differential behavior of nat-cholesterol isomers, ent- and epi-cholesterol. In simulations containing the endogenous agonist, serotonin, cholesterol binding at transmembrane helix 4 biases bound serotonin molecules toward an unexpected binding mode in the extended binding pocket. The results highlight the capability of membrane cholesterol to influence the mobility of the extracellular surface in the 5-HT1 receptor family and manipulate the architecture of the extracellular ligand-binding pocket.

Phospholipid Membrane Fluidity Alters Ligand Binding Activity of a G Protein-Coupled Receptor by Shifting the Conformational Equilibrium.

The affinity of a ligand for a receptor on the cell surface will be influenced by the membrane composition. Herein, we evaluated the effects of differences in membrane fluidity, controlled by phospholipid composition, on the ligand binding activity of the G protein-coupled receptor human serotonin 2B. Using Nanodisc technology to control membrane properties, we performed biophysical analysis and employed molecular dynamics simulations to demonstrate that increased membrane fluidity shifted the equilibrium toward an active form of the receptor. Our quantitative study will enable development of more realistic in vitro drug discovery assays involving membrane-bound proteins such as G protein-coupled receptors.

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

Molecular Signatures of Cholesterol Interaction with Serotonin Receptors.

The interaction of G protein-coupled receptors (GPCRs) with cholesterol is a hallmark of their function, organization, and structural dynamics. Several cholesterol interaction sites, such as the cholesterol recognition amino acid consensus (CRAC) and cholesterol consensus motif (CCM), have been mapped from crystallography, bioinformatics, and simulation studies. In this article, we characterize common descriptors for cholesterol interaction sites in the serotonin1A receptor from a series of coarse-grain simulations. We have identified a novel interaction mode for cholesterol in which the cholesterol polar headgroup interacts with aromatic amino acid residues, such as tryptophan and tyrosine. The cholesterol rings interact with both aromatic residues and nonpolar residues, thereby constituting a signature aromatic interaction site. In addition, we report a similar binding mode in the crystal structures of the serotonin2B receptor, suggesting that this binding mode could be a general feature of the serotonin receptor family. Interestingly, this signature aromatic interaction site is present along with one of the CRAC motifs in the serotonin1A receptor. Our results represent an important step toward mapping out the diversity of cholesterol-GPCR interaction sites.

Assessment of the transmembrane domain structures in GPCR Dock 2013 models.

The community-wide blind prediction of G-protein coupled receptor (GPCR) structures and ligand docking has been conducted three times and the quality of the models was primarily assessed by the accuracy of ligand binding modes. The seven transmembrane (TM) helices of the receptors were taken as a whole; thus the model quality within the 7TM domains has not been evaluated. Here we evaluate the 7TM domain structures in the models submitted for the last round of prediction - GPCR Dock 2013. Applying the 7 × 7 RMSD matrix analysis described in our prior work, we show that the models vary widely in prediction accuracy of the 7TM structures, exhibiting diverse structural differences from the targets. For the prediction of the 5-hydroxytryptamine receptors, the top 7TM models are rather close to the targets, which however are not ranked top by ligand-docking. On the other hand, notable deviations of the TMs are found in in the previously identified top docking models that closely resemble other receptors. We further reveal reasons of success and failure in ligand docking for the models. This current assessment not only complements the previous assessment, but also provides important insights into the current status of GPCR modeling and ligand docking.

Development of Fluorescent Probes that Target Serotonin 5-HT2B Receptors.

Some 5-HT2B fluorescent probes were obtained by tagging 1-(2,5-dimethoxy-4-iodophenyl)-propan-2-amine (DOI) with a subset of fluorescent amines. Some of the resulting fluorescent ligands showed excellent affinity and selectivity profiles at the 5-HT2B receptors (e.g. 12b), while retain the agonistic functional behaviour of the model ligand (DOI). The study highlighted the most salient features of the structure-activity relationship in this series and these were substantiated by a molecular modelling study based on a receptor-driven docking model constructed on the basis of the crystal structure of the human 5-HT2B receptor. One of the fluorescent ligands developed in this work, compound 12i, specifically labelled CHO-K1 cells expressing 5-HT2B receptors and not parental CHO-K1 cells in a concentration-dependent manner. 12i enables imaging and quantification of specific 5-HT2B receptor labelling in live cells by automated fluorescence microscopy as well as quantification by measurements of fluorescence intensity using a fluorescence plate reader.

Structural insights into the extracellular recognition of the human serotonin 2B receptor by an antibody.

Monoclonal antibodies provide an attractive alternative to small-molecule therapies for a wide range of diseases. Given the importance of G protein-coupled receptors (GPCRs) as pharmaceutical targets, there has been an immense interest in developing therapeutic monoclonal antibodies that act on GPCRs. Here we present the 3.0-Å resolution structure of a complex between the human 5-hydroxytryptamine 2B (5-HT2B) receptor and an antibody Fab fragment bound to the extracellular side of the receptor, determined by serial femtosecond crystallography with an X-ray free-electron laser. The antibody binds to a 3D epitope of the receptor that includes all three extracellular loops. The 5-HT2B receptor is captured in a well-defined active-like state, most likely stabilized by the crystal lattice. The structure of the complex sheds light on the mechanism of selectivity in extracellular recognition of GPCRs by monoclonal antibodies.

Mechanism Exploration of Arylpiperazine Derivatives Targeting the 5-HT2A Receptor by In Silico Methods.

As a G-protein coupled receptor, the 5-hydroxytryptamine 2A (5-HT2A) receptor is known for its critical role in the cognitive, behavioural and physiological functions, and thus is a primary molecular target to treat psychiatric diseases, including especially depression. With purpose to explore the structural traits affecting the inhibitory activity, currently a dataset of 109 arylpiperazine derivatives as promising 5-HT2A antagonists was built, based on which the ligand-based three-dimensional quantitative structure-activity relationship (3D-QSAR) study by using both comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) approaches was carried out. The resultant optimal CoMSIA model displays proper validity and predictability with cross-validated correlation coefficient Q² = 0.587, non-cross-validated correlation coefficient R²ncv = 0.900 and predicted correlation coefficient for the test set of compounds R²pre = 0.897, respectively. Besides, molecular docking was also conducted to investigate the binding mode between these ligands and the active site of the 5-HT2A receptor. Meanwhile, as a docking supplementary tool to study the antagonists' conformation in the binding cavity, molecular dynamics (MD) simulation was also performed, providing further elucidation about the changes in the ligand-receptor complex. Lastly, some new molecules were also newly-designed based on the above results that are potential arylpiperazine antagonists of 5-HT2A receptor. We hope that the present models and derived information may be of help for facilitating the optimization and design of novel potent antagonists as antidepressant drugs as well as exploring the interaction mechanism of 5-HT2A antagonists.

Crystal Structure of an LSD-Bound Human Serotonin Receptor.

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

Selective 5-HT2C receptor agonists: Design and synthesis of pyridazine-fused azepines.

Heterocycle-fused azepines are discussed as potent 5-HT2C receptor agonists with excellent selectivity over 5-HT2B agonism. Synthesis and structure activity relationships are outlined for a series of bicyclic pyridazino[3,4-d]azepines. By comparison with earlier published work, in vitro assays predict a high probability for achieving CNS penetration for a potent and selective compound 15a, a pre-requisite to achieve in vivo efficacy.

Receptor interaction profiles of novel psychoactive tryptamines compared with classic hallucinogens.

The present study investigated interactions between the novel psychoactive tryptamines DiPT, 4-OH-DiPT, 4-OH-MET, 5-MeO-AMT, and 5-MeO-MiPT at monoamine receptors and transporters compared with the classic hallucinogens lysergic acid diethylamide (LSD), psilocin, N,N-dimethyltryptamine (DMT), and mescaline. We investigated binding affinities at human monoamine receptors and determined functional serotonin (5-hydroxytryptamine [5-HT]) 5-HT2A and 5-HT2B receptor activation. Binding at and the inhibition of human monoamine uptake transporters and transporter-mediated monoamine release were also determined. All of the novel tryptamines interacted with 5-HT2A receptors and were partial or full 5-HT2A agonists. Binding affinity to the 5-HT2A receptor was lower for all of the tryptamines, including psilocin and DMT, compared with LSD and correlated with the reported psychoactive doses in humans. Several tryptamines, including psilocin, DMT, DiPT, 4-OH-DiPT, and 4-OH-MET, interacted with the serotonin transporter and partially the norepinephrine transporter, similar to 3,4-methylenedioxymethamphetamine but in contrast to LSD and mescaline. LSD but not the tryptamines interacted with adrenergic and dopaminergic receptors. In conclusion, the receptor interaction profiles of the tryptamines predict hallucinogenic effects that are similar to classic serotonergic hallucinogens but also MDMA-like psychoactive properties.

Application of an Integrated GPCR SAR-Modeling Platform To Explain the Activation Selectivity of Human 5-HT2C over 5-HT2B.

Agonism of the 5-HT2C serotonin receptor has been associated with the treatment of a number of diseases including obesity, psychiatric disorders, sexual health, and urology. However, the development of effective 5-HT2C agonists has been hampered by the difficulty in obtaining selectivity over the closely related 5-HT2B receptor, agonism of which is associated with irreversible cardiac valvulopathy. Understanding how to design selective agonists requires exploration of the structural features governing the functional uniqueness of the target receptor relative to related off targets. X-ray crystallography, the major experimental source of structural information, is a slow and challenging process for integral membrane proteins, and so is currently not feasible for every GPCR or GPCR-ligand complex. Therefore, the integration of existing ligand SAR data with GPCR modeling can be a practical alternative to provide this essential structural insight. To demonstrate this, we integrated SAR data from 39 azepine series 5-HT2C agonists, comprising both selective and unselective examples, with our hierarchical GPCR modeling protocol (HGMP). Through this work we have been able to demonstrate how relatively small differences in the amino acid sequences of GPCRs can lead to significant differences in secondary structure and function, as supported by experimental data. In particular, this study suggests that conformational differences in the tilt of TM7 between 5-HT2B and 5-HT2C, which result from differences in interhelical interactions, may be the major source of selectivity in G-protein activation between these two receptors. Our approach also demonstrates how the use of GPCR models in conjunction with SAR data can be used to explain activity cliffs.

Structure and function of serotonin G protein-coupled receptors.

Serotonin receptors are prevalent throughout the nervous system and the periphery, and remain one of the most lucrative and promising drug discovery targets for disorders ranging from migraine headaches to neuropsychiatric disorders such as schizophrenia and depression. There are 14 distinct serotonin receptors, of which 13 are G protein-coupled receptors (GPCRs), which are targets for approximately 40% of the approved medicines. Recent crystallographic and biochemical evidence has provided a converging understanding of the basic structure and functional mechanics of GPCR activation. Currently, two GPCR crystal structures exist for the serotonin family, the 5-HT1B and 5-HT2B receptor, with the antimigraine and valvulopathic drug ergotamine bound. The first serotonin crystal structures not only provide the first evidence of serotonin receptor topography but also provide mechanistic explanations into functional selectivity or biased agonism. This review will detail the findings of these crystal structures from a molecular and mutagenesis perspective for driving rational drug design for novel therapeutics incorporating biased signaling.

Molecular interactions between general anesthetics and the 5HT2B receptor.

BACKGROUND: Serotonin modulates many processes through a family of seven serotonin receptors. However, no studies have screened for interactions between general anesthetics currently in clinical use and serotonergic G-protein-coupled receptors (GPCRs). Given that both intravenous and inhalational anesthetics have been shown to target other classes of GPCRs, we hypothesized that general anesthetics might interact directly with some serotonin receptors and thus modify their function. METHODS: Radioligand binding assays were performed to screen serotonin receptors for interactions with propofol and isoflurane as well as for affinity determinations. Docking calculations using the crystal structure of 5-HT2B were performed to computationally confirm the binding assay results and locate anesthetic binding sites. RESULTS: The 5-HT2B class of receptors interacted significantly with both propofol and isoflurane in the primary screen. The affinities for isoflurane and propofol were determined to be 7.78 and .95 µM, respectively, which were at or below the clinical concentrations for both anesthetics. The estimated free energy derived from docking calculations for propofol (-6.70 kcal/mol) and isoflurane (-5.10 kcal/mol) correlated with affinities from the binding assay. The anesthetics were predicted to dock at a pharmacologically relevant binding site of 5HT2B. CONCLUSIONS: The molecular interactions between propofol and isoflurane with the 5-HT2B class of receptors were discovered and characterized. This finding implicates the serotonergic GPCRs as potential anesthetic targets.

Multidimensional de†novo design reveals 5-HT2B receptor-selective ligands.

We report a multi-objective de novo design study driven by synthetic tractability and aimed at the prioritization of computer-generated 5-HT2B receptor ligands with accurately predicted target-binding affinities. Relying on quantitative bioactivity models we designed and synthesized structurally novel, selective, nanomolar, and ligand-efficient 5-HT2B modulators with sustained cell-based effects. Our results suggest that seamless amalgamation of computational activity prediction and molecular design with microfluidics-assisted synthesis enables the swift generation of small molecules with the desired polypharmacology.

A dynamic view of molecular switch behavior at serotonin receptors: implications for functional selectivity.

Functional selectivity is a property of G protein-coupled receptors that allows them to preferentially couple to particular signaling partners upon binding of biased agonists. Publication of the X-ray crystal structure of serotonergic 5-HT1B and 5-HT2B receptors in complex with ergotamine, a drug capable of activating G protein coupling and ß-arrestin signaling at the 5-HT1B receptor but clearly favoring ß-arrestin over G protein coupling at the 5-HT2B subtype, has recently provided structural insight into this phenomenon. In particular, these structures highlight the importance of specific residues, also called micro-switches, for differential receptor activation. In our work, we apply classical molecular dynamics simulations and enhanced sampling approaches to analyze the behavior of these micro-switches and their impact on the stabilization of particular receptor conformational states. Our analysis shows that differences in the conformational freedom of helix 6 between both receptors could explain their different G protein-coupling capacity. In particular, as compared to the 5-HT1B receptor, helix 6 movement in the 5-HT2B receptor can be constrained by two different mechanisms. On the one hand, an anchoring effect of ergotamine, which shows an increased capacity to interact with the extracellular part of helices 5 and 6 and stabilize them, hinders activation of a hydrophobic connector region at the center of the receptor. On the other hand, this connector region in an inactive conformation is further stabilized by unconserved contacts extending to the intracellular part of the 5-HT2B receptor, which hamper opening of the G protein binding site. This work highlights the importance of considering receptor capacity to adopt different conformational states from a dynamic perspective in order to underpin the structural basis of functional selectivity.

Synthesis and structure-activity relationships of new carbonyl guanidine derivatives as novel dual 5-HT2B and 5-HT7 receptor antagonists. Part 2.

We previously reported that the novel dual 5-HT2B and 5-HT7 receptor antagonist N-(9-hydroxy-9H-fluorene-2-carbonyl)guanidine (4) exerted a suppressing effect on 5-HT-induced dural protein extravasation in guinea pigs. To develop a synthetic strategy, we performed docking studies of lead compound 4 bound to 5-HT2B and 5-HT7 receptors, and observed that the carbonyl guanidine group forms a tight interaction network with an active center Asp (D135:5-HT2B, D162:5-HT7), Tyr (Y370:5-HT2B, Y374:5-HT7) and aromatic residue (W131:5-HT2B, F158:5-HT7). Based on molecular modeling results, we optimized the substituents at the 5- to 8-position and 9-position of the fluorene ring and identified N-(diaminomethylene)-9-hydroxy-9-methyl-9H-fluorene-2-carboxamide (24a) exhibits potent affinity for 5-HT2B (Ki=4.3 nM) and 5-HT7 receptor (Ki=4.3 nM) with high selectivity over 5-HT2A, 5-HT2C, α1, D2 and M1 receptors. Compound 24a reversed the hypothermic effect of 5-carboxamidotryptamine (5-CT) in mice and also showed a suppressing effect on 5-HT-induced dural protein extravasation in guinea pigs when orally administered at 30 mg/kg. Compound 24a is therefore a promising candidate for a novel class of anti-migraine agent without any adverse effects.

5-Hydroxy-2-(2-phenylethyl)chromone (5-HPEC): a novel non-nitrogenous ligand for 5-HT2B receptor.

Chromones are a class of natural products found in almost every known terrestrial plant with over 4000 naturally occurring derivatives having been isolated and structurally elucidated. Recently, 5-hydroxy-2-(2-phenylethyl)chromone (5-HPEC), isolated from Imperata cylindrical, showed neuroprotective activity against glutamate induced excitotoxicity in primary cultures of rat cortical cells. In comparison to other naturally occurring neuroprotective chromones, 5-HPEC contains fewer hydroxyl groups. Here we report our most recent characterization on this interesting natural product against a number of CNS receptors for the purpose to identify the potential molecular targets that may be related to its biological activity. Based on our studies, including radiobinding assays, calcium flux functional assays and molecular modeling studies, 5-HPEC may represent a type of novel nonnitrogenous ligands to the 5-HT2B receptor.