Identification of Two Serine Residues Involved in Agonist Activation of the β-Adrenergic Receptor

Cathodoluminescent and Characteristic X-Ray-Emissive Rare-Earth-Doped Core/Shell Protein Labels for Spectromicroscopic Analysis of Cell Surface Receptors

Understanding the localization and the interactions of biomolecules at the nanoscale and in the cellular context remains challenging. Electron microscopy (EM), unlike light-based microscopy, gives access to the cellular ultrastructure yet results in grey-scale images and averts unambiguous (co-)localization of biomolecules. Multimodal nanoparticle-based protein labels for correlative cathodoluminescence electron microscopy (CCLEM) and energy-dispersive X-ray spectromicroscopy (EDX-SM) are presented. The single-particle STEM-cathodoluminescence (CL) and characteristic X-ray emissivity of sub-20 nm lanthanide-doped nanoparticles are exploited as unique spectral fingerprints for precise label localization and identification. To maximize the nanoparticle brightness, lanthanides are incorporated in a low-phonon host lattice and separated from the environment using a passivating shell. The core/shell nanoparticles are then functionalized with either folic (terbium-doped) or caffeic acid (europium-doped). Their potential for (protein-)labeling is successfully demonstrated using HeLa cells expressing different surface receptors that bind to folic or caffeic acid, respectively. Both particle populations show single-particle CL emission along with a distinctive energy-dispersive X-ray signal, with the latter enabling color-based localization of receptors within swift imaging times well below 2 min perμ m $\umu\text{m}$ 2 while offering high resolution with a pixel size of 2.78 nm. Taken together, these results open a route to multi-color labeling based on electron spectromicroscopy.

Unraveling Activation-Related Rearrangements and Intrinsic Divergence from Ligand-Specific Conformational Changes of the Dopamine D3 and D2 Receptors

Effective rational drug discovery hinges on understanding the functional states of the target protein and distinguishing it from homologues. However, for the G protein coupled receptors, both activation-related conformational changes (ACCs) and intrinsic divergence among receptors can be misled or obscured by ligand-specific conformational changes (LCCs). Here, we unraveled ACCs and intrinsic divergence from LCCs of the dopamine D3 and D2 receptors (D3R and D2R), by analyzing their experimentally determined structures and the molecular dynamics (MD) simulation results of the receptors bound with various ligands. In addition to the ACCs common to other aminergic receptors, we revealed unique ACCs for these two receptors, including the extracellular portion of TM5 (TM5e) and TM6e shifting away from TM2e and TM3e, with a subtle rotation of TM5e. In identifying intrinsic divergence, we found more outward tilting of TM6e in the D2R compared to the D3R in both the experimental structures and simulations bound with ligands in different scaffolds. However, this difference was drastically reduced in the simulations bound with nonselective agonist quinpirole, suggesting a misleading effect of LCCs. Further, in the quinpirole-bound simulations, TM1 showed a greater disparity between these receptors, indicating that LCCs may also obscure intrinsic divergence. Importantly, our MD simulations revealed divergence in the dynamics of these receptors. Specifically, the D2R exhibited heightened flexibility compared to the D3R in the extracellular loops and TMs 5e, 6e, and 7e, associated with its greater ligand binding site plasticity. Our results lay the groundwork for crafting ligands specifically targeting the D2R and D3R with more precise pharmacological profiles.

A review on the effect of COX-2-mediated mechanisms on development and progression of gastric cancer induced by nicotine

Smoking is a documented risk factor for cancer, e.g., gastric cancer. Nicotine, the principal tobacco alkaloid, would exert its role of contribution to gastric cancer development and progression through nicotinic acetylcholine receptors (nAChRs) and β-adrenergic receptors (β-ARs), which then promote cancer cell proliferation, migration and invasion. As a key isoenzyme in conversion of arachidonic acid to prostaglandins, cyclooxygenase-2 (COX-2) has been demonstrated to have a wide range of effects in carcinogenesis and tumor development. At present, many studies have reported the effect of nicotine on gastric cancer by binding to nAChR, as well as indirectly stimulating β-AR to mediate COX-2-related pathways. This review summarizes these studies, and also proposes more potential COX-2-mediated mechanisms. These events might contribute to the growth and progression of gastric cancer exposed to nicotine through tobacco smoke or cigarette substitutes. Also, this review article has therefore the potential not only to make a significant contribution to the treatment and prognosis of gastric cancer for smokers but also to the clinical application of COX-2 antagonists. In addition, this work also discusses the considerable challenges of this field with special reference to the future perspective of COX-2-mediated mechanisms in development and progression of gastric cancer induced by nicotine.

Inside the Biology of the β3-Adrenoceptor

Since the first discovery in 1989, the β3-adrenoceptor (β3-AR) has gained great attention because it showed the ability to regulate many physiologic and metabolic activities, such as thermogenesis and lipolysis in brown and white adipose tissue, respectively (BAT, WAT), negative inotropic effects in cardiomyocytes, and relaxation of the blood vessels and the urinary bladder. The β3-AR has been suggested as a potential target for cancer treatment, both in adult and pediatric tumors, since under hypoxia its upregulation in the tumor microenvironment (TME) regulates stromal cell differentiation, tumor growth and metastases, signifying that its agonism/antagonism could be useful for clinical benefits. Promising results in cancer research have proposed the β3-AR being targeted for the treatment of many conditions, with some drugs, at present, undergoing phase II and III clinical trials. In this review, we report the scientific journey followed by the research from the β3-Ars' discovery, with focus on the β3-Ars' role in cancer initiation and progression that elects it an intriguing target for novel antineoplastic approaches. The overview highlights the great potential of the β3-AR, both in physiologic and pathologic conditions, with the intention to display the possible benefits of β3-AR modulation in cancer reality.

Predicting the duration of action of β2-adrenergic receptor agonists: Ligand and structure-based approaches

Agonists of the β2 adrenergic receptor (ADRB2) are an important class of medications used for the treatment of respiratory diseases. They can be classified as short acting (SABA) or long acting (LABA), with each class playing a different role in patient management. In this work we explored both ligand-based and structure-based high-throughput approaches to classify β2-agonists based on their duration of action. A completely in-silico prediction pipeline using an AlphaFold generated structure was used for structure-based modelling. Our analysis identified the ligands' 3D structure and lipophilicity as the most relevant features for the prediction of the duration of action. Interaction-based methods were also able to select ligands with the desired duration of action, incorporating the bias directly in the structure-based drug discovery pipeline without the need for further processing.

Unraveling activation-related rearrangements and intrinsic divergence from ligand-induced conformational changes of the dopamine D3 and D2 receptors

Effective rational drug discovery targeting a specific protein hinges on understanding their functional states and distinguishing it from homologs. However, for the G protein coupled receptors, both the activation-related conformational changes (ACCs) and the intrinsic divergence among receptors can be misled or obscured by ligand-induced conformational changes (LCCs). Here, we unraveled ACCs and intrinsic divergence from LCCs of the dopamine D3 and D2 receptors (D3R and D2R), by analyzing their experimentally determined structures and the molecular dynamics simulation results of the receptors bound with different ligands. In addition to the ACCs common to other aminergic receptors, we revealed unique ACCs for these two receptors including TM5e and TM6e shifting away from TM2e and TM3e, with a subtle rotation of TM5e. In identifying intrinsic divergence, we found pronounced outward tilting of TM6e in the D2R compared to the D3R in both experimental structures and simulations with ligands in different scaffolds. This tilting was drastically reduced in the simulations of the receptors bound with nonselective full agonist quinpirole, suggesting a misleading impact of LCCs. Further, in the quinpirole-bound simulations, TM1 showed a greater disparity between these receptors, indicating that LCCs may obscure intrinsic divergence. In addition, our analysis showed that the impact of the nonconserved TM1 propagated to conserved Trp7.40 and Glu2.65, both are ligand binding residues. We also found that the D2R exhibited heightened flexibility compared to the D3R in the extracellular portions of TMs 5, 6, and 7, potentially associated with its greater ligand binding site plasticity. Our results lay the groundwork for crafting ligands specifically targeting D2R or D3R with more precise pharmacological profiles.

Activation/Deactivation Free-Energy Profiles for the β2-Adrenergic Receptor: Ligand Modes of Action

We use enhanced-sampling simulations with an effective collective variable to study the activation of the β2-adrenergic receptor in the presence of ligands with different efficacy. The free-energy profiles are computed for the ligand-free (apo) receptor and binary (apo-receptor + G-protein α-subunit and receptor + ligand) and ternary complexes. The results are not only compatible with available experiments but also allow unprecedented structural insight into the nature of GPCR conformations along the activation pathway and their role in the activation mechanism. In particular, the simulations reveal an unexpected mode of action of partial agonists such as salmeterol and salbutamol that arises already in the binary complex without the G-protein. Specific differences in the polar interactions with residues in TM5, which are required to stabilize an optimal TM6 conformation that facilitates G-protein binding and receptor activation, play a major role in differentiating them from full agonists.

A Novel "Activation Switch" Motif Common to All Aminergic Receptors

Aminergic receptors are G protein-coupled receptors (GPCRs) that transduce signals from small endogenous biogenic amines to regulate intracellular signaling pathways. Agonist binding in the ligand binding pocket on the extracellular side opens and prepares a cavity on the intracellular face of the receptors to interact with and activate G proteins and β-arrestins. Here, by reviewing and analyzing all available aminergic receptor structures, we seek to identify activation-related conformational changes that are independent of the specific scaffold of the bound agonist, which we define as "activation conformational changes" (ACCs). While some common intracellular ACCs have been well-documented, identifying common extracellular ACCs, including those in the ligand binding pocket, is complicated by local adjustments to different ligand scaffolds. Our analysis shows no common ACCs at the extracellular ends of the transmembrane helices. Furthermore, the restricted access to the ligand binding pocket identified previously in some receptors is not universal. Notably, the Trp6.48 toggle switch and the Pro5.50-Ile3.40-Phe6.44 (PIF) motif at the bottom of the ligand binding pocket have previously been proposed to mediate the conformational consequences of ligand binding to the intracellular side of the receptors. Our analysis shows that common ACCs in the ligand binding pocket are associated with the PIF motif and nearby residues, including Trp6.48, but fails to support a shared rotamer toggle associated with activation. However, we identify two common rearrangements between the extracellular and middle subsegments, and propose a novel "activation switch" motif common to all aminergic receptors. This motif includes the middle subsegments of transmembrane helices 3, 5, and 6 and integrates both the PIF motif and Trp6.48.

General Metadynamics Protocol To Simulate Activation/Deactivation of Class A GPCRs: Proof of Principle for the Serotonin Receptor

We present a generally applicable metadynamics protocol for characterizing the activation free-energy profiles of class A G-protein coupled receptors and a proof-of-principle study for the 5HT_1A-receptor. The almost universal A¹⁰⁰ activation index, which depends on five inter-helix distances, is used as the single collective variable in well-tempered multiple-walker metadynamics simulations. Here, we show free-energy profiles for the serotonin receptor as binary (apo-receptor + G-protein-α-subunit and receptor + ligand) and ternary complexes with two prototypical orthosteric ligands: the full agonist serotonin and the partial agonist aripiprazole. Our results are not only compatible with previously reported experimental and computational data, but they also allow differences between active and inactive conformations to be determined in unprecedented atomic detail, and with respect to the so-called microswitches that have been suggested as determinants of activation, giving insight into their role in the activation mechanism.

Unconventional interactions of the TRPV4 ion channel with beta-adrenergic receptor ligands

The transient receptor potential vanilloid 4 (TRPV4) ion channel is present in different tissues including those of the airways. This channel is activated in response to stimuli such as changes in temperature, hypoosmotic conditions, mechanical stress, and chemicals from plants, lipids, and others. TRPV4's overactivity and/or dysfunction has been associated with several diseases, such as skeletal dysplasias, neuromuscular disorders, and lung pathologies such as asthma and cardiogenic lung edema and COVID-19-related respiratory malfunction. TRPV4 antagonists and blockers have been described; nonetheless, the mechanisms involved in achieving inhibition of the channel remain scarce, and the search for safe use of these molecules in humans continues. Here, we show that the widely used bronchodilator salbutamol and other ligands of β-adrenergic receptors inhibit TRPV4's activation. We also demonstrate that inhibition of TRPV4 by salbutamol is achieved through interaction with two residues located in the outer region of the pore and that salbutamol leads to channel closing, consistent with an allosteric mechanism. Our study provides molecular insights into the mechanisms that regulate the activity of this physiopathologically important ion channel.

L-DOPA and Droxidopa: From Force Field Development to Molecular Docking into Human β2-Adrenergic Receptor

The increasing interest in the molecular mechanism of the binding of different agonists and antagonists to β₂-adrenergic receptor (β₂AR) inactive and active states has led us to investigate protein–ligand interactions using molecular docking calculations. To perform this study, the 3.2 Å X-ray crystal structure of the active conformation of human β₂AR in the complex with the endogenous agonist adrenaline has been used as a template for investigating the binding of two exogenous catecholamines to this adrenergic receptor. Here, we show the derivation of L-DOPA and Droxidopa OPLS all atom (AA) force field (FF) parameters via quantum mechanical (QM) calculations, molecular dynamics (MD) simulations in aqueous solutions of the two catecholamines and the molecular docking of both ligands into rigid and flexible β₂AR models. We observe that both ligands share with adrenaline similar experimentally observed binding anchor sites, which are constituted by Asp113/Asn312 and Ser203/Ser204/Ser207 side chains. Moreover, both L-DOPA and Droxidopa molecules exhibit binding affinities comparable to that predicted for adrenaline, which is in good agreement with previous experimental and computational results. L-DOPA and Droxidopa OPLS AA FFs have also been tested by performing MD simulations of these ligands docked into β₂AR proteins embedded in lipid membranes. Both hydrogen bonds and hydrophobic interaction networks observed over the 1 μs MD simulation are comparable with those derived from molecular docking calculations and MD simulations performed with the CHARMM FF.

Evaluating the conformational space of the active site of D2 dopamine receptor. Scope and limitations of the standard docking methods

We report here for the first time the potential energy surfaces (PES) of phenyletilamine (PEA) and meta-tyramine (m-OH-PEA) at the D₂ dopamine receptor (D2DR) binding site. PESs not only allow us to observe all the critical points of the surface (minimums, maximums, and transition states), but also to note the ease or difficulty that each local minima have for their conformational inter-conversions and therefore know the conformational flexibility that these ligands have in their active sites. Taking advantage of possessing this valuable information, we analyze how accurate a standard docking study is in these cases. Our results indicate that although we have to be careful in how to carry out this type of study and to consider performing some extra-simulations, docking calculations can be satisfactory. In order to analyze in detail the different molecular interactions that are stabilizing the different ligand-receptor (L-R) complexes, we carried out quantum theory of atoms in molecules (QTAIM) computations and NMR shielding calculations. Although some of these techniques are a bit tedious and require more computational time, our results demonstrate the importance of performing computational simulations using different types of combined techniques (docking/MD/hybrid QM-MM/QTAIM and NMR shielding calculations) in order to obtain more accurate results. Our results allow us to understand in details the molecular interactions stabilizing and destabilizing the different L-R complexes reported here. Thus, the different activities observed for dopamine (DA), m-OH-PEA, and PEA can be clearly explained at molecular level.

Analysis of L-DOPA and droxidopa binding to human β2-adrenergic receptor

Over the last two decades, an increasing number of studies has been devoted to a deeper understanding of the molecular process involved in the binding of various agonists and antagonists to active and inactive conformations of β2-adrenergic receptor (β2AR). The 3.2 Å x-ray crystal structure of human β2AR active state in combination with the endogenous low affinity agonist adrenaline offers an ideal starting structure for studying the binding of various catecholamines to adrenergic receptors. We show that molecular docking of levodopa (L-DOPA) and droxidopa into rigid and flexible β2AR models leads for both ligands to binding anchor sites comparable to those experimentally reported for adrenaline, namely D113/N312 and S203/S204/S207 side chains. Both ligands have a hydrogen bond network that is extremely similar to those of noradrenaline and dopamine. Interestingly, redocking neutral and protonated versions of adrenaline to rigid and flexible β2AR models results in binding poses that are more energetically stable and distinct from the x-ray crystal structure. Similarly, lowest energy conformations of noradrenaline and dopamine generated by docking into flexible β2AR models had binding free energies lower than those of best poses in rigid receptor models. Furthermore, our findings show that L-DOPA and droxidopa molecules have binding affinities comparable to those predicted for adrenaline, noradrenaline, and dopamine, which are consistent with previous experimental and computational findings and supported by the molecular dynamics simulations of β2AR-ligand complexes performed here.

Molecular and pharmacological characterization of biogenic amine receptors from the diamondback moth, Plutella xylostella

BACKGROUD

Insect biogenic amines play important roles in mediating behavioral and physiological processes. They exert their effects by binding to biogenic amine receptors (BARs), which are specific receptor proteins in the G-protein-coupled receptor superfamily. BAR genes have been cloned and characterized from multiple model insects, including Drosophila melanogaster, Anopheles gambiae, Bombyx mori, Apis mellifera and Tribolium castaneum. However, relatively little work has addressed the molecular properties, expression profiles, and pharmacological characterization of BARs from other insects, including important pests.

RESULTS

In this study, we cloned 17 genes encoding putative biogenic amine receptor proteins from Plutella xylostella, a global pest of Brassica crops. These PxBAR genes were five octopamine receptors (PxOA1, PxOA2B1, PxOA2B2, PxOA2B3, and PxOA3), three tyramine receptors (PxTAR1A, PxTAR1B, and PxTAR2), four dopamine receptors (PxDOP1, PxDOP2, PxDOP3, and PxDopEcR), and five serotonin receptors (Px5-HT_1A , Px5-HT_1B , Px5-HT_2A , Px5-HT_2B , and Px5-HT₇ ). All PxBARs showed considerable sequence identity with orthologous BARs, and phylogenetic analysis clustered the receptors within their respective groups while preserving organismal evolutionary relationships. We investigated their molecular properties and expression profiles, and pharmacologically characterized the dopamine receptor, PxDOP2.

CONCLUSIONS

Our study provides important information and resources on biogenic amine receptors from P. xylostella, which suggests potential target sites for controlling this pest species. © 2021 Society of Chemical Industry.

Dopamine receptors in emesis: Molecular mechanisms and potential therapeutic function

Dopamine is a member of the catecholamine family and is associated with multiple physiological functions. Together with its five receptor subtypes, dopamine is closely linked to neurological disorders such as schizophrenia, Parkinson's disease, depression, attention deficit-hyperactivity, and restless leg syndrome. Unfortunately, several dopamine receptor-based agonists used to treat some of these diseases cause nausea and vomiting as impending side-effects. The high degree of cross interactions of dopamine receptor ligands with many other targets including G-protein coupled receptors, transporters, enzymes, and ion-channels, add to the complexity of discovering new targets for the treatment of nausea and vomiting. Using activation status of signaling cascades as mechanism-based biomarkers to foresee drug sensitivity combined with the development of dopamine receptor-based biased agonists may hold great promise and seems as the next step in drug development for the treatment of such multifactorial diseases. In this review, we update the present knowledge on dopamine and dopamine receptors and their potential roles in nausea and vomiting. The pre- and clinical evidence provided in this review supports the implication of both dopamine and dopamine receptor agonists in the incidence of emesis. Besides the conventional dopaminergic antiemetic drugs, potential novel antiemetic targeting emetic protein signaling cascades may offer superior selectivity profile and potency.

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

Activation of the α2B adrenoceptor by the sedative sympatholytic dexmedetomidine

The α₂ adrenergic receptors (α₂ARs) are G protein-coupled receptors (GPCRs) that respond to adrenaline and noradrenaline and couple to the Gi/o family of G proteins. α₂ARs play important roles in regulating the sympathetic nervous system. Dexmedetomidine is a highly selective α₂AR agonist used in post-operative patients as an anxiety-reducing, sedative medicine that decreases the requirement for opioids. As is typical for selective αAR agonists, dexmedetomidine consists of an imidazole ring and a substituted benzene moiety lacking polar groups, which is in contrast to βAR-selective agonists, which share an ethanolamine group and an aromatic system with polar, hydrogen-bonding substituents. To better understand the structural basis for the selectivity and efficacy of adrenergic agonists, we determined the structure of the α_2BAR in complex with dexmedetomidine and Go at a resolution of 2.9 Å by single-particle cryo-EM. The structure reveals the mechanism of α₂AR-selective activation and provides insights into Gi/o coupling specificity.

Crystal structure of dopamine receptor D4 bound to the subtype selective ligand, L745870

Multiple subtypes of dopamine receptors within the GPCR superfamily regulate neurological processes through various downstream signaling pathways. A crucial question about the dopamine receptor family is what structural features determine the subtype-selectivity of potential drugs. Here, we report the 3.5-angstrom crystal structure of mouse dopamine receptor D4 (DRD4) complexed with a subtype-selective antagonist, L745870. Our structure reveals a secondary binding pocket extended from the orthosteric ligand-binding pocket to a DRD4-specific crevice located between transmembrane helices 2 and 3. Additional mutagenesis studies suggest that the antagonist L745870 prevents DRD4 activation by blocking the relative movement between transmembrane helices 2 and 3. These results expand our knowledge of the molecular basis for the physiological functions of DRD4 and assist new drug design.

β-blockers augment L-type Ca2+ channel activity by targeting spatially restricted β2AR signaling in neurons

G protein-coupled receptors (GPCRs) transduce pleiotropic intracellular signals in mammalian cells. Here, we report neuronal excitability of β-blockers carvedilol and alprenolol at clinically relevant nanomolar concentrations. Carvedilol and alprenolol activate β₂AR, which promote G protein signaling and cAMP/PKA activities without action of G protein receptor kinases (GRKs). The cAMP/PKA activities are restricted within the immediate vicinity of activated β₂AR, leading to selectively enhance PKA-dependent phosphorylation and stimulation of endogenous L-type calcium channel (LTCC) but not AMPA receptor in rat hippocampal neurons. Moreover, we have engineered a mutant β₂AR that lacks the catecholamine binding pocket. This mutant is preferentially activated by carvedilol but not the orthosteric agonist isoproterenol. Carvedilol activates the mutant β₂AR in mouse hippocampal neurons augmenting LTCC activity through cAMP/PKA signaling. Together, our study identifies a mechanism by which β-blocker-dependent activation of GPCRs promotes spatially restricted cAMP/PKA signaling to selectively target membrane downstream effectors such as LTCC in neurons.

Ligand-binding affinity of alternative conformers of human β2 -adrenergic receptor in the presence of intracellular loop 3 (ICL3) and their potential use in virtual screening studies

This study investigates the structural distinctiveness of orthosteric ligand-binding sites of several human β₂ adrenergic receptor (β₂ -AR) conformations that have been obtained from a set of independent molecular dynamics (MD) simulations in the presence of intracellular loop 3 (ICL3). A docking protocol was established in order to classify each receptor conformation via its binding affinity to selected ligands with known efficacy. This work's main goal was to reveal many subtle features of the ligand-binding site, presenting alternative conformations, which might be considered as either active- or inactive-like but mostly specific for that ligand. Agonists, inverse agonists, and antagonists were docked to each MD conformer with distinct binding pockets, using different docking tools and scoring functions. Mostly favored receptor conformation persistently observed in all docking/scoring evaluations was classified as active or inactive based on the type of ligand's biological effect. Classified MD conformers were further tested for their ability to discriminate agonists from inverse agonists/antagonists, and several conformers were proposed as important targets to be used in virtual screening experiments that were often limited to a single X-ray structure.

A critical role for Dop1-mediated dopaminergic signaling in the plasticity of behavioral and neuronal responses to sex pheromone in a moth

Most animal species, including insects, are able to modulate their responses to sexual chemosignals and this flexibility originates from the remodeling of olfactory areas under the influence of dopaminergic system. In the moth Agrotis ipsilon, the behavioral response of males to the female-emitted sex pheromone increases throughout adult life and after a prior exposure to pheromone signal and this change is accompanied by an increase in neuronal sensitivity within the primary olfactory centers, the antennal lobes (ALs). To identify the underlying neuromodulatory mechanisms, we examined whether this age- and experience-dependent olfactory plasticity is mediated by dopamine (DA) through the Dop1 receptor, an ortholog of the vertebrate D1-type dopamine receptors, which is positively coupled to adenylyl cyclase. We cloned A. ipsilon Dop1 (AiDop1) which is expressed predominantly in brain and especially in ALs and its knockdown induced decreased AL cAMP amounts and altered sex pheromone-orientated flight. The levels of DA, AiDop1 expression and cAMP in ALs increased from the third day of adult life and at 24h and 48h following pre-exposure to sex pheromone and the dynamic of these changes correlated with the increased responsiveness to sex pheromone. These results demonstrate that Dop1 is required for the display of male sexual behavior and that age- and experience-related neuronal and behavioral changes are sustained by DA-Dop1 signaling that operates within ALs probably through cAMP-dependent mechanisms in A. ipsilon. Thus, this study expands our understanding of the neuromodulatory mechanisms underlying olfactory plasticity, mechanisms that appear to be highly conserved between insects and mammals.

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.

Pharmacological characterization of dopamine receptors in the rice striped stem borer, Chilo suppressalis

Dopamine is an important neurotransmitter and neuromodulator in both vertebrates and invertebrates and is the most abundant monoamine present in the central nervous system of insects. A complement of functionally distinct dopamine receptors mediate the signal transduction of dopamine by modifying intracellular Ca²⁺ and cAMP levels. In the present study, we pharmacologically characterized three types of dopamine receptors, CsDOP1, CsDOP2 and CsDOP3, from the rice striped stem borer, Chilo suppressalis. All three receptors show considerable sequence identity with orthologous dopamine receptors. The phylogenetic analysis also clusters the receptors within their respective groups. Transcript levels of CsDOP1, CsDOP2 and CsDOP3 were all expressed at high levels in the central nervous system, indicating their important roles in neural processes. After heterologous expression in HEK 293 cells, CsDOP1, CsDOP2 and CsDOP3 were dose-dependently activated by dopamine and synthetic dopamine receptor agonists. They can also be blocked by different series of antagonists. This study offers important information on three dopamine receptors from C. suppressalis that will provide the basis for forthcoming studies investigating their roles in behaviors and physiology, and facilitate the development of new insecticides for pest control.

Functional roles of T3.37 and S5.46 in the activation mechanism of the dopamine D1 receptor

The activation mechanism of dopamine receptors is unknown. The amino acids S5.42, S5.43, and S5.46 located in helix 5 appear to be crucial, but their specific roles in receptor activation have not been studied. We modeled the D1 dopamine receptor using the crystal structures of the D3 dopamine and β2 adrenergic receptors. Molecular dynamics simulations show that the interaction of dopamine with the D1 receptor leads to the formation of a hydrogen-bond network with its catechol group and helices 3, 5, and 6, including water molecules. The para hydroxyl group of dopamine binds directly to S5.42 and N6.55, the latter also interacting with S5.43. Unexpectedly, S5.46 does not interact directly with the catechol; instead, it interacts through a water molecule with S5.42 and directly with T3.37. The formation of this hydrogen-bond network, part of which was previously observed in docking studies with dopamine agonists, triggers the opening of the E6.30-R3.60 ionic lock associated with the activation of GPCRs. These changes do not occur in the unbonded (apo) receptor or when it is in a complex with the antagonist 3-methoxy-5,6,7,8,9,14-hexahydrodibenz[d,g]azecine. Our results provide valuable insight into the T3.37-S5.46-water-S5.43-ligand interaction, which may be crucial to the activation of the D1 dopamine receptor and should be considered during the design of novel agonists. Graphical Abstract General representation of the relationship between the formation of the HBN and the opening of the R3.50-E6.30 ionic lock.

Confirming therapeutic target of protopine using immobilized β2 -adrenoceptor coupled with site-directed molecular docking and the target-drug interaction by frontal analysis and injection amount-dependent method

Drug-protein interaction analysis is pregnant in designing new leads during drug discovery. We prepared the stationary phase containing immobilized β₂ -adrenoceptor (β₂ -AR) by linkage of the receptor on macroporous silica gel surface through N,N'-carbonyldiimidazole method. The stationary phase was applied in identifying antiasthmatic target of protopine guided by the prediction of site-directed molecular docking. Subsequent application of immobilized β₂ -AR in exploring the binding of protopine to the receptor was realized by frontal analysis and injection amount-dependent method. The association constants of protopine to β₂ -AR by the 2 methods were (1.00 ± 0.06) × 10⁵ M^-1 and (1.52 ± 0.14) × 10⁴ M^-1 . The numbers of binding sites were (1.23 ± 0.07) × 10^-7 M and (9.09 ± 0.06) × 10^-7 M, respectively. These results indicated that β₂ -AR is the specific target for therapeutic action of protopine in vivo. The target-drug binding occurred on Ser¹⁶⁹ in crystal structure of the receptor. Compared with frontal analysis, injection amount-dependent method is advantageous to drug saving, improvement of sampling efficiency, and performing speed. It has grave potential in high-throughput drug-receptor interaction analysis.

Conserved Mechanism of Conformational Stability and Dynamics in G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are transmembrane receptors involved in diverse biological functions. Despite the diversity in their amino acid sequences, class A GPCRs exhibit a conserved structural topology and possibly a common mechanism of receptor activation. To understand how this high sequence diversity translates to a conserved functional mechanism, we have compared the dynamic behavior of eight class A GPCRs comprised of six biogenic amine receptors, adenosine A_2A, and the peptide receptor protease-activated receptor 1. Starting from the crystal structures of the inactive state of these receptors bound to inverse agonists or antagonists, we have performed multiple all-atom MD simulations adding up to several microseconds of simulation. We elucidated the similarities and differences in the dynamic behavior and the conformational ensembles sampled by these eight class A GPCRs. Among the six biogenic amine receptors studied here, β₂-adrenergic receptor shows the highest level of fluctuation in the sixth and seventh transmembrane helices, possibly explaining its high basal activity. In contrast, the muscarinic acetylcholine receptors show the lowest fluctuations as well as tight packing and low hydration of the transmembrane domain. All eight GPCRs show several conserved allosteric communication pipelines from the residues in the agonist binding site with the G-protein interface. Positions of the residues along these pipelines that serve as major hubs of allosteric communication are conserved in their respective structures. These findings have important implications in understanding the dynamics and allosteric mechanism of communication in class A GPCRs and hence are useful for designing conformation-specific drugs.

How To Minimize Artifacts in Atomistic Simulations of Membrane Proteins, Whose Crystal Structure Is Heavily Engineered: β₂-Adrenergic Receptor in the Spotlight

Atomistic molecular dynamics (MD) simulations are used extensively to elucidate membrane protein properties. These simulations are based on three-dimensional protein structures that in turn are often based on crystallography. The protein structures resolved in crystallographic studies typically do not correspond to pristine proteins, however. Instead the crystallized proteins are commonly engineered, including structural modifications (mutations, replacement of protein sequences by antibodies, bound ligands, etc.) whose impact on protein structure and dynamics is largely unknown. Here we explore this issue through atomistic MD simulations (∼5 μs in total), focusing on the β2-adrenergic receptor (β2AR) that is one of the most studied members of the G-protein coupled receptor superfamily. Starting from an inactive-state crystal structure of β2AR, we remove the many modifications in β2AR systematically one at a time, in six consecutive steps. After each step, we equilibrate the system and simulate it quite extensively. The results of this step-by-step approach highlight that the structural modifications used in crystallization can affect ligand and G-protein binding sites, packing at the transmembrane-helix interface region, and the dynamics of connecting loops in β2AR. When the results of the systematic step-by-step approach are compared to an all-at-once technique where all modifications done on β2AR are removed instantaneously at the same time, it turns out that the step-by-step method provides results that are superior in terms of maintaining protein structural stability. The results provide compelling evidence that for membrane proteins whose 3D structure is based on structural engineering, the preparation of protein structure for atomistic MD simulations is a delicate and sensitive process. The results show that most valid results are found when the structural modifications are reverted slowly, one at a time.

Molecular Evolution and Functional Divergence of Trace Amine-Associated Receptors

Trace amine-associated receptors (TAARs) are a member of the G-protein-coupled receptor superfamily and are known to be expressed in olfactory sensory neurons. A limited number of molecular evolutionary studies have been done for TAARs so far. To elucidate how lineage-specific evolution contributed to their functional divergence, we examined 30 metazoan genomes. In total, 493 TAAR gene candidates (including 84 pseudogenes) were identified from 26 vertebrate genomes. TAARs were not identified from non-vertebrate genomes. An ancestral-type TAAR-like gene appeared to have emerged in lamprey. We found four therian-specific TAAR subfamilies (one eutherian-specific and three metatherian-specific) in addition to previously known nine subfamilies. Many species-specific TAAR gene duplications and losses contributed to a large variation of TAAR gene numbers among mammals, ranging from 0 in dolphin to 26 in flying fox. TAARs are classified into two groups based on binding preferences for primary or tertiary amines as well as their sequence similarities. Primary amine-detecting TAARs (TAAR1-4) have emerged earlier, generally have single-copy orthologs (very few duplication or loss), and have evolved under strong functional constraints. In contrast, tertiary amine-detecting TAARs (TAAR5-9) have emerged more recently and the majority of them experienced higher rates of gene duplications. Protein members that belong to the tertiary amine-detecting TAAR group also showed the patterns of positive selection especially in the area surrounding the ligand-binding pocket, which could have affected ligand-binding activities and specificities. Expansions of the tertiary amine-detecting TAAR gene family may have played important roles in terrestrial adaptations of therian mammals. Molecular evolution of the TAAR gene family appears to be governed by a complex, species-specific, interplay between environmental and evolutionary factors.

Pharmacological Analysis and Structure Determination of 7-Methylcyanopindolol-Bound β1-Adrenergic Receptor

Comparisons between structures of the β1-adrenergic receptor (AR) bound to either agonists, partial agonists, or weak partial agonists led to the proposal that rotamer changes of Ser(5.46), coupled to a contraction of the binding pocket, are sufficient to increase the probability of receptor activation. (RS)-4-[3-(tert-butylamino)-2-hydroxypropoxy]-1H-indole-2-carbonitrile (cyanopindolol) is a weak partial agonist of β1AR and, based on the hypothesis above, we predicted that the addition of a methyl group to form 4-[(2S)-3-(tert-butylamino)-2-hydroxypropoxy]-7-methyl-1H-indole-2-carbonitrile (7-methylcyanopindolol) would dramatically reduce its efficacy. An eight-step synthesis of 7-methylcyanopindolol was developed and its pharmacology was analyzed. 7-Methylcyanopindolol bound with similar affinity to cyanopindolol to both β1AR and β2AR. As predicted, the efficacy of 7-methylcyanopindolol was reduced significantly compared with cyanopindolol, acting as a very weak partial agonist of turkey β1AR and an inverse agonist of human β2AR. The structure of 7-methylcyanopindolol-bound β1AR was determined to 2.4-Å resolution and found to be virtually identical to the structure of cyanopindolol-bound β1AR. The major differences in the orthosteric binding pocket are that it has expanded by 0.3 Å in 7-methylcyanopindolol-bound β1AR and the hydroxyl group of Ser(5.46) is positioned 0.8 Å further from the ligand, with respect to the position of the Ser(5.46) side chain in cyanopindolol-bound β1AR. Thus, the molecular basis for the reduction in efficacy of 7-methylcyanopindolol compared with cyanopindolol may be regarded as the opposite of the mechanism proposed for the increase in efficacy of agonists compared with antagonists.

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.

Pharmacological and signalling properties of a D2-like dopamine receptor (Dop3) in Tribolium castaneum

Dopamine is an important neurotransmitter in the central nervous system of vertebrates and invertebrates. Despite their evolutionary distance, striking parallels exist between deuterostomian and protostomian dopaminergic systems. In both, signalling is achieved via a complement of functionally distinct dopamine receptors. In this study, we investigated the sequence, pharmacology and tissue distribution of a D2-like dopamine receptor from the red flour beetle Tribolium castaneum (TricaDop3) and compared it with related G protein-coupled receptors in other invertebrate species. The TricaDop3 receptor-encoding cDNA shows considerable sequence similarity with members of the Dop3 receptor class. Real time qRT-PCR showed high expression in both the central brain and the optic lobes, consistent with the role of dopamine as neurotransmitter. Activation of TricaDop3 expressed in mammalian cells increased intracellular Ca(2+) signalling and decreased NKH-477 (a forskolin analogue)-stimulated cyclic AMP levels in a dose-dependent manner. We studied the pharmacological profile of the TricaDop3 receptor and demonstrated that the synthetic vertebrate dopamine receptor agonists, 2 - amino- 6,7 - dihydroxy - 1,2,3,4 - tetrahydronaphthalene hydrobromide (6,7-ADTN) and bromocriptine acted as agonists. Methysergide was the most potent of the antagonists tested and showed competitive inhibition in the presence of dopamine. This study offers important information on the Dop3 receptor from Tribolium castaneum that will facilitate functional analyses of dopamine receptors in insects and other invertebrates.

Characterization of the Anopheles gambiae octopamine receptor and discovery of potential agonists and antagonists using a combined computational-experimental approach

BACKGROUND

Octopamine receptors (OARs) perform key functions in the biological pathways of primarily invertebrates, making this class of G-protein coupled receptors (GPCRs) a potentially good target for insecticides. However, the lack of structural and experimental data for this insect-essential GPCR family has promoted the development of homology models that are good representations of their biological equivalents for in silico screening of small molecules.

METHODS

Two Anopheles gambiae OARs were cloned, analysed and functionally characterized using a heterologous cell reporter system. Four antagonist- and four agonist-binding homology models were generated and virtually screened by docking against compounds obtained from the ZINC database. Resulting compounds from the virtual screen were tested experimentally using an in vitro reporter assay and in a mosquito larvicide bioassay.

RESULTS

Six An. gambiae OAR/tyramine receptor genes were identified. Phylogenetic analysis revealed that the OAR (AGAP000045) that encodes two open reading frames is an α-adrenergic-like receptor. Both splice variants signal through cAMP and calcium. Mutagenesis analysis revealed that D100 in the TM3 region and S206 and S210 in the TM5 region are important to the activation of the GPCR. Some 2,150 compounds from the virtual screen were structurally analysed and 70 compounds were experimentally tested against AgOAR45B expressed in the GloResponse™CRE-luc2P HEK293 reporter cell line, revealing 21 antagonists, 17 weak antagonists, 2 agonists, and 5 weak agonists.

CONCLUSION

Reported here is the functional characterization of two An. gambiae OARs and the discovery of new OAR agonists and antagonists based on virtual screening and molecular dynamics simulations. Four compounds were identified that had activity in a mosquito larva bioassay, three of which are imidazole derivatives. This combined computational and experimental approach is appropriate for the discovery of new and effective insecticides.

Strategies for improved modeling of GPCR-drug complexes: blind predictions of serotonin receptors bound to ergotamine

The recent increase in the number of atomic-resolution structures of G protein-coupled receptors (GPCRs) has contributed to a deeper understanding of ligand binding to several important drug targets. However, reliable modeling of GPCR-ligand complexes for the vast majority of receptors with unknown structure remains to be one of the most challenging goals for computer-aided drug design. The GPCR Dock 2013 assessment, in which researchers were challenged to predict the crystallographic structures of serotonin 5-HT(1B) and 5-HT(2B) receptors bound to ergotamine, provided an excellent opportunity to benchmark the current state of this field. Our contributions to GPCR Dock 2013 accurately predicted the binding mode of ergotamine with RMSDs below 1.8 Å for both receptors, which included the best submissions for the 5-HT(1B) complex. Our models also had the most accurate description of the binding sites and receptor-ligand contacts. These results were obtained using a ligand-guided homology modeling approach, which combines extensive molecular docking screening with incorporation of information from multiple crystal structures and experimentally derived restraints. In this work, we retrospectively analyzed thousands of structures that were generated during the assessment to evaluate our modeling strategies. Major contributors to accuracy were found to be improved modeling of extracellular loop two in combination with the use of molecular docking to optimize the binding site for ligand recognition. Our results suggest that modeling of GPCR-drug complexes has reached a level of accuracy at which structure-based drug design could be applied to a large number of pharmaceutically relevant targets.

Immobilised histidine tagged β2-adrenoceptor oriented by a diazonium salt reaction and its application in exploring drug-protein interaction using ephedrine and pseudoephedrine as probes

A new oriented method using a diazonium salt reaction was developed for linking β₂-adrenoceptor (β₂-AR) on the surface of macroporous silica gel. Stationary phase containing the immobilised receptor was used to investigate the interaction between β₂-AR and ephedrine plus pseudoephedrine by zonal elution. The isotherms of the two drugs best fit the Langmuir model. Only one type of binding site was found for ephedrine and pseudoephedrine targeting β₂-AR. At 37 °C, the association constants during the binding were (5.94±0.05)×10³/M for ephedrine and (3.80±0.02) ×10³/M for pseudoephedrine, with the binding sites of (8.92±0.06) ×10⁻⁴ M. Thermodynamic studies showed that the binding of the two compounds to β₂-AR was a spontaneous reaction with exothermal processes. The ΔG^θ, ΔH^θ and ΔS^θ for the interaction between ephedrine and β₂-AR were −(22.33±0.04) kJ/mol, −(6.51±0.69) kJ/mol and 50.94±0.31 J/mol·K, respectively. For the binding of pseudoephedrine to the receptor, these values were −(21.17±0.02) kJ/mol, −(7.48±0.56) kJ/mol and 44.13±0.01 J/mol·K. Electrostatic interaction proved to be the driving force during the binding of the two drugs to β₂-AR. The proposed immobilised method will have great potential for attaching protein to solid substrates and realizing the interactions between proteins and drugs.

Beyond standard molecular dynamics: investigating the molecular mechanisms of G protein-coupled receptors with enhanced molecular dynamics methods

The majority of biological processes mediated by G Protein-Coupled Receptors (GPCRs) take place on timescales that are not conveniently accessible to standard molecular dynamics (MD) approaches, notwithstanding the current availability of specialized parallel computer architectures, and efficient simulation algorithms. Enhanced MD-based methods have started to assume an important role in the study of the rugged energy landscape of GPCRs by providing mechanistic details of complex receptor processes such as ligand recognition, activation, and oligomerization. We provide here an overview of these methods in their most recent application to the field.

Molecular pharmacology and ligand docking studies reveal a single amino acid difference between mouse and human serotonin 5-HT2A receptors that impacts behavioral translation of novel 4-phenyl-2-dimethylaminotetralin ligands

During translational studies to develop 4-phenyl-2-dimethylaminotetralin (PAT) compounds for neuropsychiatric disorders, the (2R,4S)-trans-(+)- and (2S,4R)-trans-(-)-enantiomers of the analog 6-hydroxy-7-chloro-PAT (6-OH-7-Cl-PAT) demonstrated unusual pharmacology at serotonin (5-HT) 5-HT2 G protein-coupled receptors (GPCRs). The enantiomers had similar affinities (Ki) at human (h) 5-HT2A receptors (≈ 70 nM). In an in vivo mouse model of 5-HT2A receptor activation [(±)-(2,5)-dimethoxy-4-iodoamphetamine (DOI)-elicited head twitch], however, (-)-6-OH-7-Cl-PAT was about 5-fold more potent than the (+)-enantiomer at attenuating the DOI-elicited response. It was discovered that (+)-6-OH-7-Cl-PAT (only) had ≈ 40-fold-lower affinity at mouse (m) compared with h5-HT2A receptors. Molecular modeling and computational ligand docking studies indicated that the 6-OH moiety of (+)- but not (-)-6-OH-7-Cl-PAT could form a hydrogen bond with serine residue 5.46 of the h5-HT2A receptor. The m5-HT2A as well as m5-HT2B, h5-HT2B, m5-HT2C, and h5-HT2C receptors have alanine at position 5.46, obviating this interaction; (+)-6-OH-7-Cl-PAT also showed ≈ 50-fold lower affinity than (-)-6-OH-7-Cl-PAT at m5-HT2C and h5-HT2C receptors. Mutagenesis studies confirmed that 5-HT2A S5.46 is critical for (+)- but not (-)-6-OH-7-Cl-PAT binding, as well as function. The (+)-6-OH-7-Cl-PAT enantiomer showed partial agonist effects at h5-HT2A wild-type (WT) and m5-HT2A A5.46S point-mutated receptors but did not activate m5-HT2A WT and h5-HT2A S5.46A point-mutated receptors, or h5-HT2B, h5-HT2C, and m5-HT2C receptors; (-)-6-OH-7-Cl-PAT did not activate any of the 5-HT2 receptors. Experiments also included the (2R,4S)-trans-(+)- and (2S,4R)-trans-(-)-enantiomers of 6-methoxy-7-chloro-PAT to validate hydrogen bonding interactions proposed for the corresponding 6-OH analogs. Results indicate that PAT ligand three-dimensional structure impacts target receptor binding and translational outcomes, supporting the hypothesis that GPCR ligand structure governs orthosteric binding pocket molecular determinants and resulting pharmacology.

Adrenaline-activated structure of β2-adrenoceptor stabilized by an engineered nanobody

G protein-coupled receptors (GPCRs) are integral membrane proteins that play an essential role in human physiology, yet the molecular processes through which they bind to their endogenous agonists and activate effector proteins remain poorly understood. Thus far, it has not been possible to capture an active-state GPCR bound to its native neurotransmitter. Crystal structures of agonist-bound GPCRs have relied on the use of either exceptionally high-affinity agonists^1,2 or receptor stabilization by mutagenesis^3-5. Many natural agonists such as adrenaline, which activates the β₂-adrenoceptor (β₂AR), bind with relatively low affinity, and they are often chemically unstable. Using directed evolution, we engineered a high-affinity camelid antibody fragment that stabilizes the active state of the β₂AR, and used this to obtain crystal structures of the activated receptor bound to multiple ligands. Here, we present structures of active-state β₂AR bound to three chemically distinct agonists: the ultra high-affinity agonist BI167107, the high-affinity catecholamine agonist hydroxybenzyl isoproterenol, and the low-affinity endogenous agonist adrenaline. The crystal structures reveal a highly conserved overall ligand recognition and activation mode despite diverse ligand chemical structures and affinities that range from 100 nM to ~80 pM. The adrenaline-bound receptor structure is overall similar to the others, but shows substantial rearrangements in extracellular loop three and the extracellular tip of transmembrane helix 6. These structures also reveal a water-mediated hydrogen bond between two conserved tyrosines, which appears to stabilize the active state of the β₂AR and related GPCRs.

Molecular interactions between fenoterol stereoisomers and derivatives and the β₂-adrenergic receptor binding site studied by docking and molecular dynamics simulations

The β2 adrenergic receptor (β2-AR) has become a model system for studying the ligand recognition process and mechanism of the G protein coupled receptors activation. In the present study stereoisomers of fenoterol and some of its derivatives (N = 94 molecules) were used as molecular probes to identify differences in stereo-recognition interactions between β2-AR and structurally similar agonists. The present study aimed at determining the 3D molecular models of the fenoterol derivative-β2-AR complexes. Molecular models of β2-AR have been developed by using the crystal structure of the human β2-AR T4 lysozyme fusion protein with bound (S)-carazolol (PDB ID: 2RH1) and more recently reported structure of a nanobody-stabilized active state of the β2-AR with the bound full agonist BI-167107 (PDB ID: 3P0G). The docking procedure allowed us to study the similarities and differences in the recognition binding site(s) for tested ligands. The agonist molecules occupied the same binding region, between TM III, TM V, TM VI and TM VII. The residues identified by us during docking procedure (Ser203, Ser207, Asp113, Lys305, Asn312, Tyr308, Asp192) were experimentally indicated in functional and biophysical studies as being very important for the agonist-receptor interactions. Moreover, the additional space, an extension of the orthosteric pocket, was identified and described. Furthermore, the molecular dynamics simulations were used to study the molecular mechanism of interaction between ligands ((R,R')- and (S,S')-fenoterol) and β2-AR. Our research offers new insights into the ligand stereoselective interaction with one of the most important GPCR member. This study may also facilitate the design of improved selective medications, which can be used to treat, prevent and control heart failure symptoms.

Determinants involved in subtype-specific functions of rat trace amine-associated receptors 1 and 4

AIMS

The trace amine-associated receptor (Taar) family displays high species- and subtype-specific pharmacology. Several trace amines such as β-phenylethylamine (β-PEA), p-tyramine and tryptamine are agonists at TA(1) but poorly activate rat and mouse Taar4.

PRINCIPAL RESULTS

Using rat TA(1) and Taar4 chimera, we identified determinants in transmembrane helices 3 and 6, which, when replaced by the corresponding portion of rat TA(1) , can rescue cell surface expression of rat Taar4. When expressed at the cell surface, rat Taar4 pharmacology was very similar to that of TA(1) and coupled to the Gα(s) -protein/AC pathway. Our data suggest that binding pockets of Taar for surrogate agonists overlap between paralogs.

CONCLUSIONS

This implicates that the repertoire of Taar ensures functional redundancy, tissue- and cell-specific expression and/or different downstream signalling rather than different agonist specificity.

The importance of interactions with helix 5 in determining the efficacy of β-adrenoceptor ligands

Structures of the inactive state of the thermostabilized β1-adrenoceptor have been determined bound to eight different ligands, including full agonists, partial agonists, inverse agonists and biased agonists. Comparison of the structures shows distinct differences within the binding pocket that correlate with the pharmacological properties of the ligands. These data suggest that full agonists stabilize a structure with a contracted binding pocket and a rotamer change of serine (5.46) compared with when antagonists are bound. Inverse agonists may prevent both of these occurrences, whereas partial agonists stabilize a contraction of the binding pocket but not the rotamer change of serine (5.46). It is likely that subtle changes in the interactions between transmembrane helix 5 (H5) and H3/H4 on agonist binding promote the formation of the activated state.

Conformation guides molecular efficacy in docking screens of activated β-2 adrenergic G protein coupled receptor

A prospective, large library virtual screen against an activated β2-adrenergic receptor (β2AR) structure returned potent agonists to the exclusion of inverse-agonists, providing the first complement to the previous virtual screening campaigns against inverse-agonist-bound G protein coupled receptor (GPCR) structures, which predicted only inverse-agonists. In addition, two hits recapitulated the signaling profile of the co-crystal ligand with respect to the G protein and arrestin mediated signaling. This functional fidelity has important implications in drug design, as the ability to predict ligands with predefined signaling properties is highly desirable. However, the agonist-bound state provides an uncertain template for modeling the activated conformation of other GPCRs, as a dopamine D2 receptor (DRD2) activated model templated on the activated β2AR structure returned few hits of only marginal potency.

Structure of β-adrenergic receptors

β-Adrenergic receptors (βARs) control key physiological functions by transducing signals encoded in catecholamine hormones and neurotransmitters to activate intracellular signaling pathways. As members of the large family of G protein-coupled receptors (GPCRs), βARs have a seven-transmembrane helix topology and signal via G protein- and arrestin-dependent pathways. Until 2007, three-dimensional structural information of GPCRs activated by diffusible ligands, including βARs, was limited to homology models that used the related photoreceptor rhodopsin as a template. Over many years, several labs have developed strategies that have finally allowed the structures of the turkey β(1)AR and the human β(2)AR to be determined experimentally. The challenges to overcome included heterologous receptor overexpression, design of stabilized and crystallizable modified receptor constructs, ligand-affinity purification of active receptor and the development of novel techniques in crystallization and microcrystallography. The structures of βARs in complex with inverse agonists, antagonists, and agonists have revealed the binding mode of ligands with different efficacies, have allowed to obtain insights into ligand selectivity, and have provided better templates for drug design. Also, the structures of β(2)AR in complex with a G protein and a G protein-mimicking nanobody have provided important insights into the mechanism of receptor activation and G protein coupling. This chapter summarizes the strategies and methods that have been successfully applied to the structural studies of βARs. These are exemplified with detailed protocols toward the structure determination of stabilized turkey β(1)AR-ligand complexes. We also discuss the spectacular insights into adrenergic receptor function that were obtained from the structures.

Identifying ligand binding conformations of the β2-adrenergic receptor by using its agonists as computational probes

Recently available G-protein coupled receptor (GPCR) structures and biophysical studies suggest that the difference between the effects of various agonists and antagonists cannot be explained by single structures alone, but rather that the conformational ensembles of the proteins need to be considered. Here we use an elastic network model-guided molecular dynamics simulation protocol to generate an ensemble of conformers of a prototypical GPCR, β(2)-adrenergic receptor (β(2)AR). The resulting conformers are clustered into groups based on the conformations of the ligand binding site, and distinct conformers from each group are assessed for their binding to known agonists of β(2)AR. We show that the select ligands bind preferentially to different predicted conformers of β(2)AR, and identify a role of β(2)AR extracellular region as an allosteric binding site for larger drugs such as salmeterol. Thus, drugs and ligands can be used as "computational probes" to systematically identify protein conformers with likely biological significance.

The essential role for aromatic cluster in the β3 adrenergic receptor

AIM

To explore the function of the conserved aromatic cluster F213(5.47), F308(6.51), and F309(6.52) in human β3 adrenergic receptor (hβ3AR).

METHODS

Point mutation technology was used to produce plasmid mutations of hβ3AR. HEK-293 cells were transiently co-transfected with the hβ3AR (wild-type or mutant) plasmids and luciferase reporter vector pCRE-luc. The expression levels of hβ3AR in the cells were determined by Western blot analysis. The constitutive signalling and the signalling induced by the β3AR selective agonist, BRL (BRL37344), were then evaluated. To further explore the interaction mechanism between BRL and β3AR, a three-dimensional complex model of β3AR and BRL was constructed by homology modelling and molecular docking.

RESULTS

For F308(6.51), Ala and Leu substitution significantly decreased the constitutive activities of β3AR to approximately 10% of that for the wild-type receptor. However, both the potency and maximal efficacy were unchanged by Ala substitution. In the F308(6.51)L construct, the EC(50) value manifested as a "right shift" of approximately two orders of magnitude with an increased E(max). Impressively, the molecular pharmacological phenotype was similar to the wild-type receptor for the introduction of Tyr at position 308(6.51), though the EC(50) value increased by approximately five-fold for the mutant. For F309(6.52), the constitutive signalling for both F309(6.52)A and F309(6.52)L constructs were strongly impaired. In the F309(6.52)A construct, BRL-stimulated signalling showed a normal E(max) but reduced potency. Leu substitution of F309(6.52) reduced both the E(max) and potency. When F309(6.52) was mutated to Tyr, the constitutive activity was decreased approximately three-fold, and BRL-stimulated signalling was significantly impaired. Furthermore, the double mutant (F308(6.51)A_F309(6.52)A) caused the total loss of β3AR function. The predicted binding mode between β3AR and BRL revealed that both F308(6.51) and F309(6.52) were in the BRL binding pocket of β3AR, while F213(5.47) and W305(6.48) were distant from the binding site.

CONCLUSION

These results revealed that aromatic residues, especially F308(6.51) and F309(6.52), play essential roles in the function of β3AR. Aromatic residues maintained the receptor in a partially activated state and significantly contributed to ligand binding. The results supported the common hypothesis that the aromatic cluster F[Y]5.47/F[Y]6.52/F[Y]6.51 conserved in class A G protein-coupled receptor (GPCR) plays an important role in the structural stability and activation of GPCRs.