The second extracellular loop of the dopamine D2 receptor lines the binding-site crevice

1-(Phenylselanyl)-2-(p-tolyl)indolizine: A selenoindolizine with potential antidepressant-like activity in mice mediated by the modulation of dopaminergic and noradrenergic systems

1-(Phenylselanyl)-2-(p-tolyl)indolizine (MeSeI) is a selenoindolizine with an antidepressant-like effect in mice by regulation of the serotonergic system. This study investigated the involvement of dopaminergic and noradrenergic systems in the antidepressant-like action of MeSeI. For this purpose, Swiss male mice were pretreated with different antagonists, after 15 min, the MeSeI was administrated by intragastric (i.g.) via; after 30 min, the mouse behavior was assessed in the forced swimming test (FST). The action of MeSeI on the activity of monoamine oxidase (MAO) was determined. The pretreatment of mice with haloperidol (0.05 mg/kg, intraperitoneally, i.p.; non-selective dopamine receptor antagonist), sulpiride (50 mg/kg, i.p.; D2 receptor antagonist), yohimbine (1 mg/kg, i.p.; α2 receptor antagonist), and propranolol (2 mg/kg, i.p.; non-selective β receptor antagonist), inhibited the anti-immobility action of MeSeI (50 mg/kg, i.g.) in the FST. This blocking effect was not observed when SCH23390 (0.01 mg/kg, i.p.; D1 receptor antagonist), and prazosin (1 mg/kg, i.p.; α1 receptor antagonist) were administered. The coadministration of subeffective doses of bupropion (3 mg/kg. i.g.; dopamine and noradrenaline reuptake inhibitor) and MeSeI (0.5 mg/kg. i.g.) reduced the immobility time in the FST. Furthermore, MeSeI inhibited MAO-A and B activities in vitro and ex vivo tests. These results suggest that MeSeI exerts its antidepressant-like effect via regulation of the D2, α2, and β1 receptors and the inhibition of MAO-A and B activities. Molecular docking investigations corroborated these results. This study provides comprehensive insights into the antidepressant-like mechanism of MeSeI in mice, suggesting its potential as a novel antidepressant candidate.

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.

Dynamics of the Second Extracellular Loop Control Transducer Coupling of Peptide-Activated GPCRs

Many peptide-activated rhodopsin-like GPCRs share a β-hairpin folding motif in the extracellular loop 2 (ECL2), which interacts with the peptide ligand while at the same time being connected to transmembrane helix 3 (TM3) via a highly conserved disulfide bond. Currently, it remains unknown whether the coupling of the specifically shaped ECL2 to TM3 influences the activation of peptide-activated GPCRs. We investigated this possibility in a selection of peptide GPCRs with known structures. Most of the receptors with cysteine to alanine mutations folded like the respective wild-type and resided in the cell membrane, challenging pure folding stabilization by the disulfide bridge. G-protein signaling of the disulfide mutants was retained to a greater extent in secretin-like GPCRs than in rhodopsin-like GPCRs, while recruitment of arrestin was completely abolished in both groups, which may be linked to alterations in ligand residence time. We found a correlation between receptor activity of the neuropeptide Y2 receptor and alterations in ECL2 dynamics using engineered disulfide bridges or site-directed spin labeling and EPR spectroscopy. These data highlight the functional importance of the TM3-ECL2 link for the activation of specific signaling pathways in peptide-activated GPCRs, which might have implications for future drug discovery.

Collagen a1(XI) structure prediction by Alphafold 2

Collagen α1(XI) is a minor fibrillar collagen involved in the critical regulation of collagen fibrils such as nucleation, assembly, and regulation of fibril diameter. The amino propeptide domain of the collagen α1(XI) is retained on the surface of the collagen fibril for an extended period of time and may play a crucial role in the interaction with extracellular matrix glycosaminoglycans and other proteins during the process of fibrillogenesis. Understanding the mechanism of action of this protein will ultimately help us understand the organization and assembly of the extracellular matrix that underlies the structural integrity of connective tissues.

Benchmarking the Accuracy of AlphaFold 2 in Loop Structure Prediction

The inhibition of protein-protein interactions is a growing strategy in drug development. In addition to structured regions, many protein loop regions are involved in protein-protein interactions and thus have been identified as potential drug targets. To effectively target such regions, protein structure is critical. Loop structure prediction is a challenging subgroup in the field of protein structure prediction because of the reduced level of conservation in protein sequences compared to the secondary structure elements. AlphaFold 2 has been suggested to be one of the greatest achievements in the field of protein structure prediction. The AlphaFold 2 predicted protein structures near the X-ray resolution in the Critical Assessment of protein Structure Prediction (CASP 14) competition in 2020. The purpose of this work is to survey the performance of AlphaFold 2 in specifically predicting protein loop regions. We have constructed an independent dataset of 31,650 loop regions from 2613 proteins (deposited after the AlphaFold 2 was trained) with both experimentally determined structures and AlphaFold 2 predicted structures. With extensive evaluation using our dataset, the results indicate that AlphaFold 2 is a good predictor of the structure of loop regions, especially for short loop regions. Loops less than 10 residues in length have an average Root Mean Square Deviation (RMSD) of 0.33 Å and an average the Template Modeling score (TM-score) of 0.82. However, we see that as the number of residues in a given loop increases, the accuracy of AlphaFold 2's prediction decreases. Loops more than 20 residues in length have an average RMSD of 2.04 Å and an average TM-score of 0.55. Such a correlation between accuracy and length of the loop is directly linked to the increase in flexibility. Moreover, AlphaFold 2 does slightly over-predict α-helices and β-strands in proteins.

Structural and Energetic Insights Into the Interaction of Niacin With the GPR109A Receptor

The transmembrane G-protein coupled receptor GPR109A has been previously shown to function as a receptor for niacin in mediating antilipolytic effects. Although administration of high doses of niacin has shown beneficial effects on lipid metabolism, however, it is often accompanied by disturbing side effects such as flushing, liver damage, glucose intolerance, or gastrointestinal problems. Thus, it is important to understand niacin-GPR109A interactions, which can be beneficial for the development of alternate drugs having antilipolytic effects with less or no side effects. To get into the structural insights on niacin binding to GPR109A, we have performed 100 nanoseconds long all-atom MD simulations of five niacin-GPR109A complexes (automatically docked pose 0, and randomly placed niacin in poses 1 to 4 in the receptor crevice) and analyzed using binding free energy calculations and H-bond analysis. Steered MD simulations were used to get an average force for niacin translocation between the bulk and the external crevice of the wild type and mutant (N86Y, W91 S, S178I, and triple mutant of all three residues) GPR109A receptors, as well as GPR109B (as a control that does not bind niacin). The H-bond analysis revealed that TMH3 residue R111 interacts with niacin in a total of 4 (poses 0 to 3) complexes, while residues C177, S178, and S179 contact niacin in complex pose 4, and all these complexes were energetically stable. According to steered MD simulations, all the GPR109A mutants and GPR109B required greater force than that of wild-type GPR109A to translocate in the external crevice, suggesting increased sterical obstacles. Thus, the residues N86 (at the junction of TMH2/ECL2), W91 (ECL2), R111 (TMH3), and ECL3 residues (C177, S178, S179) play an important role for optimal routing of niacin entry and to bind GPR109A.

Neuronal Dopamine D3 Receptors: Translational Implications for Preclinical Research and CNS Disorders

Dopamine (DA), as one of the major neurotransmitters in the central nervous system (CNS) and periphery, exerts its actions through five types of receptors which belong to two major subfamilies such as D1-like (i.e., D1 and D5 receptors) and D2-like (i.e., D2, D3 and D4) receptors. Dopamine D3 receptor (D3R) was cloned 30 years ago, and its distribution in the CNS and in the periphery, molecular structure, cellular signaling mechanisms have been largely explored. Involvement of D3Rs has been recognized in several CNS functions such as movement control, cognition, learning, reward, emotional regulation and social behavior. D3Rs have become a promising target of drug research and great efforts have been made to obtain high affinity ligands (selective agonists, partial agonists and antagonists) in order to elucidate D3R functions. There has been a strong drive behind the efforts to find drug-like compounds with high affinity and selectivity and various functionality for D3Rs in the hope that they would have potential treatment options in CNS diseases such as schizophrenia, drug abuse, Parkinson’s disease, depression, and restless leg syndrome. In this review, we provide an overview and update of the major aspects of research related to D3Rs: distribution in the CNS and periphery, signaling and molecular properties, the status of ligands available for D3R research (agonists, antagonists and partial agonists), behavioral functions of D3Rs, the role in neural networks, and we provide a summary on how the D3R-related drug research has been translated to human therapy.

Structure of the dopamine D2 receptor in complex with the antipsychotic drug spiperone

In addition to the serotonin 5-HT_2A receptor (5-HT_2AR), the dopamine D₂ receptor (D₂R) is a key therapeutic target of antipsychotics for the treatment of schizophrenia. The inactive state structures of D₂R have been described in complex with the inverse agonists risperidone (D₂R_ris) and haloperidol (D₂R_hal). Here we describe the structure of human D₂R in complex with spiperone (D₂R_spi). In D₂R_spi, the conformation of the extracellular loop (ECL) 2, which composes the ligand-binding pocket, was substantially different from those in D₂R_ris and D₂R_hal, demonstrating that ECL2 in D₂R is highly dynamic. Moreover, D₂R_spi exhibited an extended binding pocket to accommodate spiperone’s phenyl ring, which probably contributes to the selectivity of spiperone to D₂R and 5-HT_2AR. Together with D₂R_ris and D₂R_hal, the structural information of D₂R_spi should be of value for designing novel antipsychotics with improved safety and efficacy.

The dopamine D₂ receptor (D₂R) is a GPCR and an important drug target for schizophrenia treatment. Here, the authors present the crystal structure of human D₂R in complex with the antipsychotic drug spiperone, which is of interest for designing antipsychotics with improved receptor selectivity.

Genetic and functional analysis of a Pacific hagfish opioid system

The actions of endogenous opioids and nociceptin/orphanin FQ are mediated by four homologous G protein-coupled receptors that constitute the opioid receptor family. However, little is known about opioid systems in cyclostomes (living jawless fish) and how opioid systems might have evolved from invertebrates. Here, we leveraged de novo transcriptome and low-coverage whole-genome assembly in the Pacific hagfish (Eptatretus stoutii) to identify and characterize the first full-length coding sequence for a functional opioid receptor in a cyclostome. Additionally, we define two novel endogenous opioid precursors in this species that predict several novel opioid peptides. Bioinformatic analysis shows no closely related opioid receptor genes in invertebrates with regard either to the genomic organization or to conserved opioid receptor-specific sequences that are common in all vertebrates. Furthermore, no proteins analogous to vertebrate opioid precursors could be identified by genomic searches despite previous claims of protein or RNA-derived sequences in several invertebrate species. The presence of an expressed orthologous receptor and opioid precursors in the Pacific hagfish confirms that a functional opioid system was likely present in the common ancestor of all extant vertebrates some 550 million years ago, earlier than all previous authenticated accounts. We discuss the premise that the cyclostome and vertebrate opioid systems evolved from invertebrate systems concerned with antimicrobial defense and speculate that the high concentrations of opioid precursors in tissues such as the testes, gut, and activated immune cells are key remnants of this evolutionary role.

Structure of a D2 dopamine receptor-G-protein complex in a lipid membrane

The D₂ dopamine receptor (DRD2) is a therapeutic target for Parkinson’s disease¹ and antipsychotic drugs². DRD2 is activated by the endogenous neurotransmitter dopamine and synthetic agonist drugs such as bromocriptine³, leading to stimulation of G_i and inhibition of adenylyl cyclase. We used cryo-electron microscopy to elucidate the structure of an agonist-bound activated DRD2-G_i complex reconstituted into a phospholipid membrane. The extracellular ligand binding site of DRD2 is remodeled in response to agonist binding, with conformational changes in extracellular loop 2 (ECL2), transmembrane domain 5 (TM5), TM6, and TM7 propagating to opening of the intracellular G_i binding site. The DRD2-G_i structure represents the first experimental model of a GPCR-G protein complex embedded in a phospholipid bilayer, which serves as a benchmark to validate the interactions seen in previous detergent-bound structures. The structure also reveals interactions that are unique to the membrane-embedded complex, including helix 8 burial in the inner leaflet, ordered lysine and arginine sidechains in the membrane interfacial regions, and lipid anchoring of the G protein in the membrane. Our model of the activated DRD2 will help inform the design of subtype-selective DRD2 ligands for multiple human CNS disorders.

DaReUS-Loop: a web server to model multiple loops in homology models

Loop regions in protein structures often have crucial roles, and they are much more variable in sequence and structure than other regions. In homology modeling, this leads to larger deviations from the homologous templates, and loop modeling of homology models remains an open problem. To address this issue, we have previously developed the DaReUS-Loop protocol, leading to significant improvement over existing methods. Here, a DaReUS-Loop web server is presented, providing an automated platform for modeling or remodeling loops in the context of homology models. This is the first web server accepting a protein with up to 20 loop regions, and modeling them all in parallel. It also provides a prediction confidence level that corresponds to the expected accuracy of the loops. DaReUS-Loop facilitates the analysis of the results through its interactive graphical interface and is freely available at http://bioserv.rpbs.univ-paris-diderot.fr/services/DaReUS-Loop/.

Cell adhesion molecule close homolog of L1 binds to the dopamine receptor D2 and inhibits the internalization of its short isoform

Cell adhesion molecule close homolog of L1 (CHL1) and the dopamine receptor D2 (DRD2) are associated with psychiatric and mental disorders. We here show that DRD2 interacts with CHL1 in mouse brain, as evidenced by co-immunostaining, proximity ligation assay, co-immunoprecipitation, and pull-down assay with recombinant extracellular CHL1 domain fused to Fc (CHL1-Fc). Direct binding of CHL1-Fc to the first extracellular loop of DRD2 was shown by ELISA. Using HEK cells transfected to co-express CHL1 and the short (DRD2-S) or long (DRD2-L) DRD2 isoforms, co-localization of CHL1 and both isoforms was observed by immunostaining and proximity ligation assay. Moreover, CHL1 inhibited agonist-triggered internalization of DRD2-S. Proximity ligation assay showed close interaction between CHL1 and DRD2 in neurons expressing dopamine- and cAMP-regulated phosphoprotein of 32 kDa (DARPP32) or tyrosine hydroxylase (TH) in tissue sections of adult mouse striatum. In cultures of striatum or ventral midbrain, CHL1 was also closely associated with DRD2 in DARPP32- or TH-immunopositive cells, respectively. In the dorsal striatum of CHL1-deficient mice, lower levels of DRD2 and phosphorylated TH were observed, when compared to wild-type littermates. In the ventral striatum of CHL1-deficient mice, levels of phosphorylated DARPP32 were reduced. We propose that CHL1 regulates DRD2-dependent presynaptic and postsynaptic functions.

Distinct inactive conformations of the dopamine D2 and D3 receptors correspond to different extents of inverse agonism

By analyzing and simulating inactive conformations of the highly homologous dopamine D2 and D3 receptors (D2R and D3R), we find that eticlopride binds D2R in a pose very similar to that in the D3R/eticlopride structure but incompatible with the D2R/risperidone structure. In addition, risperidone occupies a sub-pocket near the Na+ binding site, whereas eticlopride does not. Based on these findings and our experimental results, we propose that the divergent receptor conformations stabilized by Na+-sensitive eticlopride and Na+-insensitive risperidone correspond to different degrees of inverse agonism. Moreover, our simulations reveal that the extracellular loops are highly dynamic, with spontaneous transitions of extracellular loop 2 from the helical conformation in the D2R/risperidone structure to an extended conformation similar to that in the D3R/eticlopride structure. Our results reveal previously unappreciated diversity and dynamics in the inactive conformations of D2R. These findings are critical for rational drug discovery, as limiting a virtual screen to a single conformation will miss relevant ligands.

Pharmacological characterization, cellular localization and expression profile of NPY receptor subtypes Y2 and Y7 in large yellow croaker, Larimichthys crocea

Neuropeptide Y (NPY) receptors are suggested to mediate the multi-physiological functions of NPY family peptides, such as food intake, in teleost fish. However, the structure and signaling of fish NPY receptors are yet to be fully elucidated. In this study, we report the cloning and characterization of two neuropeptide Y receptor subtypes, Y2 (NPY2R) and Y7 (NPY7R), in yellow croaker Larimichthys crocea (L. crocea) (LcNPY2R, LcNPY7R). The gene structure, pharmacological characterization, cell location, and tissue expression of these two receptors were explored. The phylogenetic results showed that LcNPY2R and LcNPY7R had typical G protein-coupled receptor profiles, associated with the Y2 subfamily, with coding sequences that are highly conserved in vertebrates. The expression of both LcNPY2R and LcNPY7R could be activated by LcNPY in HEK293 cells. However, truncated LcNPY_18-36 was only able to activate LcNPY2R at the same level as full length LcNPY. Expression analysis revealed that LcNPY2R mRNA was predominantly expressed in the intestine and liver, whereas LcNPY7R was expressed in the stomach, which indicated that both receptors were related to the digestive system. Overall, our data establishes a molecular basis to determine the actions of LcNPY2R and LcNPY7R, which could be used to elucidate the conserved roles of these receptor-ligand pairs in vertebrates.

Molecular Determinants of the Intrinsic Efficacy of the Antipsychotic Aripiprazole

Partial agonists of the dopamine D₂ receptor (D₂R) have been developed to treat the symptoms of schizophrenia without causing the side effects elicited by antagonists. The receptor-ligand interactions that determine the intrinsic efficacy of such drugs, however, are poorly understood. Aripiprazole has an extended structure comprising a phenylpiperazine primary pharmacophore and a 1,2,3,4-tetrahydroquinolin-2-one secondary pharmacophore. We combined site-directed mutagenesis, analytical pharmacology, ligand fragments, and molecular dynamics simulations to identify the D₂R-aripiprazole interactions that contribute to affinity and efficacy. We reveal that an interaction between the secondary pharmacophore of aripiprazole and a secondary binding pocket defined by residues at the extracellular portions of transmembrane segments 1, 2, and 7 determines the intrinsic efficacy of aripiprazole. Our findings reveal a hitherto unappreciated mechanism for fine-tuning the intrinsic efficacy of D₂R agonists.

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.

Molecular dynamics simulation of biased agonists at the dopamine D2 receptor suggests the mechanism of receptor functional selectivity

The dopamine D2 receptor (D2R) is the primary target for antipsychotic drugs. Besides schizophrenia, this receptor is linked to dementia, Parkinson's disease, and depression. Recent studies have shown that β-arrestin biased agonists at this receptor treat schizophrenia with less side effects. Although the high resolution structure of this receptor exists, the mechanism of biased agonism at the receptor is unknown. In this study, dopamine, the endogenous unbiased G-protein agonist, MLS1547, a G-protein biased agonist, and UNC9975, a G-protein antagonist and a β-arrestin biased agonist, were docked to a homology model of the whole D2R including all flexible loops, and molecular dynamics simulations were conducted to study the potential mechanisms of biased agonism. Our thorough analysis on the protein-ligand interaction, secondary structure, tertiary structure, structure dynamics, and molecular switches of all three systems indicates that ligand binding to transmembrane 3 might be essential for G-protein recruitment, while ligand binding to transmembrane 6 might be essential for β-arrestin recruitment. Our analysis also suggests changes in both the secondary and the tertiary structures of TM5 and TM7, molecular switches and ICL3 flexibility are important in biased signaling. Communicated by Ramaswamy H. Sarma.

DaReUS-Loop: accurate loop modeling using fragments from remote or unrelated proteins

Despite efforts during the past decades, loop modeling remains a difficult part of protein structure modeling. Several approaches have been developed in the framework of crystal structures. However, for homology models, the modeling of loops is still far from being solved. We propose DaReUS-Loop, a data-based approach that identifies loop candidates mining the complete set of experimental structures available in the Protein Data Bank. Candidate filtering relies on local conformation profile-profile comparison, together with physico-chemical scoring. Applied to three different template-based test sets, DaReUS-Loop shows significant increase in the number of high-accuracy loops, and significant enhancement for modeling long loops. A special advantage is that our method proposes a prediction confidence score that correlates well with the expected accuracy of the loops. Strikingly, over 50% of successful loop models are derived from unrelated proteins, indicating that fragments under similar constraints tend to adopt similar structure, beyond mere homology.

Dynamic States of the Ligand-Free Class A G Protein-Coupled Receptor Extracellular Side

G protein-coupled receptors (GPCRs) make up the largest family of drug targets. The second extracellular loop (ECL2) and extracellular end of the third transmembrane helix (TM3) are basic structural elements of the GPCR ligand binding site. Currently, the disulfide bond between the two conserved cysteines in the ECL2 and TM3 is considered to be a basic GPCR structural feature. This disulfide bond has a significant effect on receptor dynamics and ligand binding. Here, molecular dynamics simulations and experimental results show that the two cysteines are distant from one another in the highest-population conformational state of ligand-free class A GPCRs and do not form a disulfide bond, indicating that the dynamics of the GPCR extracellular side are different from our conventional understanding. These surprising dynamics should have important effects on the drug binding process. On the basis of the two distinct ligand-free states, we suggest two kinetic processes for binding of ligands to GPCRs. These results challenge our commonly held beliefs regarding both GPCR structural features and ligand binding.

Structural Investigation of the Dopamine-2 Receptor Agonist Bromocriptine Binding to Dimeric D2HighR and D2LowR States

The active (D2^HighR) and inactive (D2^LowR) states of dimeric dopamine D2 receptor (D2R) models were investigated to clarify the binding mechanisms of the dopamine agonist bromocriptine, using Molecular Dynamics (MD) simulation. The aim of this comprehensive study was to investigate the critical effects of bromocriptine binding on each distinct receptor conformation. The different binding modes of the bromocriptine ligand in the active and inactive states have a significant effect on the conformational changes of the receptor. Based on the MM/GBSA approach, the calculated binding enthalpies of bromocriptine demonstrated selectivity toward the D2^HighR active state. There is good agreement between the calculated and experimentally measured D2^HighR selectivity. In the ligand-binding site, the key amino acids identified for D2^HighR were Asp114(3.32) and Glu95(2.65), and for D2^LowR, it was Ser193(5.42). Moreover, analysis of replicate MD trajectories demonstrated that the bromocriptine structure was more rigid at the D2^HighR state and more flexible at the D2^LowR state. However, the side chains of the ligand-receptor complex of D2^HighR showed larger variations relative to the corresponding regions of D2^LowR. The present study is part of an ongoing research program to study D2R conformational changes during ligand activation and to evaluate the conformational state selectivity for ligand binding.

The structural determinants of the bitopic binding mode of a negative allosteric modulator of the dopamine D2 receptor

SB269652 is a negative allosteric modulator of the dopamine D₂ receptor (D₂R) yet possesses structural similarity to ligands with a competitive mode of interaction. In this study, we aimed to understand the ligand-receptor interactions that confer its allosteric action. We combined site-directed mutagenesis with molecular dynamics simulations using both SB269652 and derivatives from our previous structure activity studies. We identify residues within the conserved orthosteric binding site (OBS) and a secondary binding pocket (SBP) that determine affinity and cooperativity. Our results indicate that interaction with the SBP is a requirement for allosteric pharmacology, but that both competitive and allosteric derivatives of SB269652 can display sensitivity to the mutation of a glutamate residue (E95^2.65) within the SBP. Our findings provide the molecular basis for the differences in affinity between SB269652 derivatives, and reveal how changes to interactions made by the primary pharmacophore of SB269652 in the orthosteric pocket can confer changes in the interactions made by the secondary pharmacophore in the SBP. Our insights provide a structure-activity framework towards rational optimization of bitopic ligands for D₂R with tailored competitive versus allosteric properties.

Serotonin 2A receptor disulfide bridge integrity is crucial for ligand binding to different signalling states but not for its homodimerization

The serotonin 2_A (5-HT_2A) receptor is a G-protein coupled receptor (GPCR) with a conserved disulfide bridge formed by Cys¹⁴⁸ (transmembrane helix 3, TM3) and Cys²²⁷ (extracellular loop 2, ECL-2). We hypothesized that disulfide bridges may determine serotonin 5-HT_2A receptor functions such as receptor activation, functional selectivity and ligand recognition. We used the reducing agent dithiothreitol (DTT) to determine how the reduction of disulfide bridges affects radioligand binding, second messenger mobilization and receptor dimerization. A DTT-induced decrease in the number of binding sites (1190 ± 63.55 fmol/mg protein for control cells compared with 921.2 ± 60.84 fmol/mg protein for DTT-treated cells) as well as in the efficacy of both signalling pathways characterized was observed, although the affinity and potency were unchanged. Bioluminiscence resonance energy transfer (BRET) assays revealed the DTT treatment did not modify the homodimeric nature of serotonin 5-HT_2A receptors. In molecular dynamic simulations, the ECL-2 of the receptor with a broken cysteine bond adopts a wider variety of conformations, some of which protrude deeper into the receptor orthosteric binding pocket leading to collapse of the pocket. A shrunken binding pocket would be incapable of accommodating lysergic acid diethylamide (LSD). Our findings suggest that the decrease of efficacy may be due to disruption of disulfide bridge between TM3 and ECL-2. This reveals the integrity of the ECL-2 epitope, which should be explored in the development of novel ligands acting as allosteric modulators of serotonin 5-HT_2A receptors.

Ligand-binding characterization of simulated β-adrenergic-like octopamine receptor in Schistocerca gregaria via progressive structure simulation

It is important to design insecticides having both low drug resistance and less undesirable toxicity for desert locust control. Specific GPCRs of Schistocerca gregaria, especially β-adrenergic-like octopamine receptor (SgOctβR), can be considered as its potential effective insecticide targets. However, either the unavailability of SgOctβR's structure or the inadequate capability of its sequence lead the development of insecticide for Schistocerca gregaria meets its plateau. To relax this difficulty, this paper develops a promising progressive structure simulation from SgOctβR's sequence, to its predicted structure of SgOctβR in vacuum, to its conformation as well as its complex with endogenous ligand octopamine in a solvent-membrane system. The combined approach of multiple sequence alignment, static structural characterization, and dynamic process of conformational change during binding octopamine reveal three important aspects. The first one is the characterization of SgOctβR's active pocket, including the attending secondary structure elements, its hydrophobic residues and nonpolar surface. The second one is the interaction with octopamine, especially the involved hydrogen bonds and an aromatic stacking of pi-pi interactions. The third one is the potential binding sites, including six highly conserved residues and one highly variable residue for locust insecticide design. This work is definitely helpful for the further structure-based drug design for efficient and eco-friendly insecticides, as well as site-directed mutagenesis biochemical research of SgOctβR.

Genetic variants of dopamine D2 receptor impact heterodimerization with dopamine D1 receptor

BACKGROUND

The human dopamine D2 receptor gene has three polymorphic variants that alter its amino acid sequence: alanine substitution by valine in position 96 (V96A), proline substitution by serine in position 310 (P310S) and serine substitution by cysteine in position 311 (S311C). Their functional role has never been the object of extensive studies, even though there is some evidence that their occurrence correlates with schizophrenia.

METHODS

The HEK293 cell line was transfected with dopamine D1 and D2 receptors (or genetic variants of the D2 receptor), coupled to fluorescent proteins which allowed us to measure the extent of dimerization of these receptors, using a highly advanced biophysical approach (FLIM-FRET). Additionally, Fluoro-4 AM was used to examine changes in the level of calcium release after ligand stimulation of cells expressing different combinations of dopamine receptors.

RESULTS

Using FLIM-FRET experiments we have shown that in HEK 293 expressing dopamine receptors, polymorphic mutations in the D2 receptor play a role in dimmer formation with the dopamine D1 receptor. The association level of dopamine receptors is affected by ligand administration, with variable effects depending on polymorphic variant of the D2 dopamine receptor. We have found that the level of heteromer formation is reflected by calcium ion release after ligand stimulation and have observed variations of this effect dependent on the polymorphic variant and the ligand.

CONCLUSION

The data presented in this paper support the hypothesis on the role of calcium signaling regulated by the D1-D2 heteromer which may be of relevance for schizophrenia etiology.

Extracellular Loop 2 of the Adenosine A1 Receptor Has a Key Role in Orthosteric Ligand Affinity and Agonist Efficacy

The adenosine A₁ G protein-coupled receptor (A₁AR) is an important therapeutic target implicated in a wide range of cardiovascular and neuronal disorders. Although it is well established that the A₁AR orthosteric site is located within the receptor's transmembrane (TM) bundle, prior studies have implicated extracellular loop 2 (ECL2) as having a significant role in contributing to orthosteric ligand affinity and signaling for various G protein-coupled receptors (GPCRs). We thus performed extensive alanine scanning mutagenesis of A₁AR-ECL2 to explore the role of this domain on A₁AR orthosteric ligand pharmacology. Using quantitative analytical approaches and molecular modeling, we identified ECL2 residues that interact either directly or indirectly with orthosteric agonists and antagonists. Discrete mutations proximal to a conserved ECL2-TM3 disulfide bond selectively affected orthosteric ligand affinity, whereas a cluster of five residues near the TM4-ECL2 juncture influenced orthosteric agonist efficacy. A combination of ligand docking, molecular dynamics simulations, and mutagenesis results suggested that the orthosteric agonist 5'-N-ethylcarboxamidoadenosine binds transiently to an extracellular vestibule formed by ECL2 and the top of TM5 and TM7, prior to entry into the canonical TM bundle orthosteric site. Collectively, this study highlights a key role for ECL2 in A₁AR orthosteric ligand binding and receptor activation.

Structural modeling of G-protein coupled receptors: An overview on automatic web-servers

Despite the significant efforts and discoveries during the last few years in G protein-coupled receptor (GPCR) expression and crystallization, the receptors with known structures to date are limited only to a small fraction of human GPCRs. The lack of experimental three-dimensional structures of the receptors represents a strong limitation that hampers a deep understanding of their function. Computational techniques are thus a valid alternative strategy to model three-dimensional structures. Indeed, recent advances in the field, together with extraordinary developments in crystallography, in particular due to its ability to capture GPCRs in different activation states, have led to encouraging results in the generation of accurate models. This, prompted the community of modelers to render their methods publicly available through dedicated databases and web-servers. Here, we present an extensive overview on these services, focusing on their advantages, drawbacks and their role in successful applications. Future challenges in the field of GPCR modeling, such as the predictions of long loop regions and the modeling of receptor activation states are presented as well.

New paradigms in chemokine receptor signal transduction: Moving beyond the two-site model

Chemokine receptor (CKR) signaling forms the basis of essential immune cellular functions, and dysregulated CKR signaling underpins numerous disease processes of the immune system and beyond. CKRs, which belong to the seven transmembrane domain receptor (7TMR) superfamily, initiate signaling upon binding of endogenous, secreted chemokine ligands. Chemokine-CKR interactions are traditionally described by a two-step/two-site mechanism, in which the CKR N-terminus recognizes the chemokine globular core (i.e. site 1 interaction), followed by activation when the unstructured chemokine N-terminus is inserted into the receptor TM bundle (i.e. site 2 interaction). Several recent studies challenge the structural independence of sites 1 and 2 by demonstrating physical and allosteric links between these supposedly separate sites. Others contest the functional independence of these sites, identifying nuanced roles for site 1 and other interactions in CKR activation. These developments emerge within a rapidly changing landscape in which CKR signaling is influenced by receptor PTMs, chemokine and CKR dimerization, and endogenous non-chemokine ligands. Simultaneous advances in the structural and functional characterization of 7TMR biased signaling have altered how we understand promiscuous chemokine-CKR interactions. In this review, we explore new paradigms in CKR signal transduction by considering studies that depict a more intricate architecture governing the consequences of chemokine-CKR interactions.

Mutation-Guided Unbiased Modeling of the Fat Sensor GPR119 for High-Yield Agonist Screening

Recent benchmark studies have demonstrated the difficulties in obtaining accurate predictions of ligand binding conformations to comparative models of G-protein-coupled receptors. We have developed a data-driven optimization protocol, which integrates mutational data and structural information from multiple X-ray receptor structures in combination with a fully flexible ligand docking protocol to determine the binding conformation of AR231453, a small-molecule agonist, in the GPR119 receptor. Resulting models converge to one conformation that explains the majority of data from mutation studies and is consistent with the structure-activity relationship for a large number of AR231453 analogs. Another key property of the refined models is their success in separating active ligands from decoys in a large-scale virtual screening. These results demonstrate that mutation-guided receptor modeling can provide predictions of practical value for describing receptor-ligand interactions and drug discovery.

Structure-Function Basis of Attenuated Inverse Agonism of Angiotensin II Type 1 Receptor Blockers for Active-State Angiotensin II Type 1 Receptor

Ligand-independent signaling by the angiotensin II type 1 receptor (AT1R) can be activated in clinical settings by mechanical stretch and autoantibodies as well as receptor mutations. Transition of the AT1R to the activated state is known to lower inverse agonistic efficacy of clinically used AT1R blockers (ARBs). The structure-function basis for reduced efficacy of inverse agonists is a fundamental aspect that has been understudied not only in relation to the AT1R but also regarding other homologous receptors. Here, we demonstrate that the active-state transition in the AT1R indeed attenuates an inverse agonistic effect of four biphenyl-tetrazole ARBs through changes in specific ligand-receptor interactions. In the ground state, tight interactions of four ARBs with a set of residues (Ser109(TM3), Phe182(ECL2), Gln257(TM6), Tyr292(TM7), and Asn295(TM7)) results in potent inverse agonism. In the activated state, the ARB-AT1R interactions shift to a different set of residues (Val108(TM3), Ser109(TM3), Ala163(TM4), Phe182(ECL2), Lys199(TM5), Tyr292(TM7), and Asn295(TM7)), resulting in attenuated inverse agonism. Interestingly, V108I, A163T, N295A, and F182A mutations in the activated state of the AT1R shift the functional response to the ARB binding toward agonism, but in the ground state the same mutations cause inverse agonism. Our data show that the second extracellular loop is an important regulator of the functional states of the AT1R. Our findings suggest that the quest for discovering novel ARBs, and improving current ARBs, fundamentally depends on the knowledge of the unique sets of residues that mediate inverse agonistic potency in the two states of the AT1R.

What can crystal structures of aminergic receptors tell us about designing subtype-selective ligands?

G protein-coupled receptors (GPCRs) are integral membrane proteins that represent an important class of drug targets. In particular, aminergic GPCRs interact with a significant portion of drugs currently on the market. However, most drugs that target these receptors are associated with undesirable side effects, which are due in part to promiscuous interactions with close homologs of the intended target receptors. Here, based on a systematic analysis of all 37 of the currently available high-resolution crystal structures of aminergic GPCRs, we review structural elements that contribute to and can be exploited for designing subtype-selective compounds. We describe the roles of secondary binding pockets (SBPs), as well as differences in ligand entry pathways to the orthosteric binding site, in determining selectivity. In addition, using the available crystal structures, we have identified conformational changes in the SBPs that are associated with receptor activation and explore the implications of these changes for the rational development of selective ligands with tailored efficacy.

FSH receptor-specific residues L501 and I505 in extracellular loop 2 are essential for its function

The extracellular loop 2 (EL2) of FSH receptor (FSHR) plays a pivotal role in various events downstream of FSH stimulation. Because swapping the six FSHR-specific residues in EL2 (chimeric EL2M) with those from LH/choriogonadotropin receptor resulted in impaired internalization of FSH-FSHR complex and low FSH-induced cAMP production, six substitution mutants of EL2 were generated to ascertain the contribution of individual amino acids to the effects shown by chimeric EL2M. Results revealed that L(501)F mainly and I(505)V to a lesser extent contribute to the diminished receptor function in chimeric EL2M. HEK293 cells stably expressing WT and chimeric EL2M FSHR were generated to track the fate of the receptors post FSH induction. The chimeric EL2M FSHR stable clone showed weak internalization and cAMP response similar to transiently transfected cells. Furthermore, reduced FSH-induced ERK phosphorylation was also observed. The interaction of activated chimeric EL2M and L(501)F FSHR with β-arrestins was weak compared with WT FSHR, thus explaining the impaired internalization of chimeric EL2M and corroborating the indispensable role of EL2 in receptor function.

The extracellular loop 2 (ECL2) of the human histamine H4 receptor substantially contributes to ligand binding and constitutive activity

In contrast to the corresponding mouse and rat orthologs, the human histamine H4 receptor (hH4R) shows extraordinarily high constitutive activity. In the extracellular loop (ECL), replacement of F169 by V as in the mouse H4R significantly reduced constitutive activity. Stabilization of the inactive state was even more pronounced for a double mutant, in which, in addition to F169V, S179 in the ligand binding site was replaced by M. To study the role of the FF motif in ECL2, we generated the hH4R-F168A mutant. The receptor was co-expressed in Sf9 insect cells with the G-protein subunits Gαi2 and Gβ1γ2, and the membranes were studied in [3H]histamine binding and functional [35S]GTPγS assays. The potency of various ligands at the hH4R-F168A mutant decreased compared to the wild-type hH4R, for example by 30- and more than 100-fold in case of the H4R agonist UR-PI376 and histamine, respectively. The high constitutive activity of the hH4R was completely lost in the hH4R-F168A mutant, as reflected by neutral antagonism of thioperamide, a full inverse agonist at the wild-type hH4R. By analogy, JNJ7777120 was a partial inverse agonist at the hH4R, but a partial agonist at the hH4R-F168A mutant, again demonstrating the decrease in constitutive activity due to F168A mutation. Thus, F168 was proven to play a key role not only in ligand binding and potency, but also in the high constitutive activity of the hH4R.

In silico characterization of binding mode of CCR8 inhibitor: homology modeling, docking and membrane based MD simulation study

Human CC-chemokine receptor 8 (CCR8) is a crucial drug target in asthma that belongs to G-protein-coupled receptor superfamily, which is characterized by seven transmembrane helices. To date, there is no X-ray crystal structure available for CCR8; this hampers active research on the target. Molecular basis of interaction mechanism of antagonist with CCR8 remains unclear. In order to provide binding site information and stable binding mode, we performed modeling, docking and molecular dynamics (MD) simulation of CCR8. Docking study of biaryl-ether-piperidine derivative (13C) was performed inside predefined CCR8 binding site to get the representative conformation of 13C. Further, MD simulations of receptor and complex (13C-CCR8) inside dipalmitoylphosphatidylcholine lipid bilayers were performed to explore the effect of lipids. Results analyses showed that the Gln91, Tyr94, Cys106, Val109, Tyr113, Cys183, Tyr184, Ser185, Lys195, Thr198, Asn199, Met202, Phe254, and Glu286 were conserved in both docking and MD simulations. This indicated possible role of these residues in CCR8 antagonism. However, experimental mutational studies on these identified residues could be effective to confirm their importance in CCR8 antagonism. Furthermore, calculated Coulombic interactions represented the crucial roles of Glu286, Lys195, and Tyr113 in CCR8 antagonism. Important residues identified in this study overlap with the previous non-peptide agonist (LMD-009) binding site. Though, the non-peptide agonist and currently studied inhibitor (13C) share common substructure, but they differ in their effects on CCR8. So, to get more insight into their agonist and antagonist effects, further side-by-side experimental studies on both agonist (LMD-009) and antagonist (13C) are suggested.

A mechanistic role of Helix 8 in GPCRs: Computational modeling of the dopamine D2 receptor interaction with the GIPC1-PDZ-domain

Helix-8 (Hx8) is a structurally conserved amphipathic helical motif in class-A GPCRs, adjacent to the C-terminal sequence that is responsible for PDZ-domain-recognition. The Hx8 segment in the dopamine D2 receptor (D2R) constitutes the C-terminal segment and we investigate its role in the function of D2R by studying the interaction with the PDZ-containing GIPC1 using homology models based on the X-ray structures of very closely related analogs: the D3R for the D2R model, and the PDZ domain of GIPC2 for GIPC1-PDZ. The mechanism of this interaction was investigated with all-atom unbiased molecular dynamics (MD) simulations that reveal the role of the membrane in maintaining the helical fold of Hx8, and with biased MD simulations to elucidate the energy drive for the interaction with the GIPC1-PDZ. We found that it becomes more favorable energetically for Hx8 to adopt the extended conformation observed in all PDZ-ligand complexes when it moves away from the membrane, and that C-terminus palmitoylation of D2R enhanced membrane penetration by the Hx8 backbone. De-palmitoylation enables Hx8 to move out into the aqueous environment for interaction with the PDZ domain. All-atom unbiased MD simulations of the full D2R-GIPC1-PDZ complex in sphingolipid/cholesterol membranes show that the D2R carboxyl C-terminus samples the region of the conserved GFGL motif located on the carboxylate-binding loop of the GIPC1-PDZ, and the entire complex distances itself from the membrane interface. Together, these results outline a likely mechanism of Hx8 involvement in the interaction of the GPCR with PDZ-domains in the course of signaling.

Docking and Virtual Screening Strategies for GPCR Drug Discovery

Progress in structure determination of G protein-coupled receptors (GPCRs) has made it possible to apply structure-based drug design (SBDD) methods to this pharmaceutically important target class. The quality of GPCR structures available for SBDD projects fall on a spectrum ranging from high resolution crystal structures (<2 Å), where all water molecules in the binding pocket are resolved, to lower resolution (>3 Å) where some protein residues are not resolved, and finally to homology models that are built using distantly related templates. Each GPCR project involves a distinct set of opportunities and challenges, and requires different approaches to model the interaction between the receptor and the ligands. In this review we will discuss docking and virtual screening to GPCRs, and highlight several refinement and post-processing steps that can be used to improve the accuracy of these calculations. Several examples are discussed that illustrate specific steps that can be taken to improve upon the docking and virtual screening accuracy. While GPCRs are a unique target class, many of the methods and strategies outlined in this review are general and therefore applicable to other protein families.

Identification, tissue distribution and functional characterization of the ghrelin receptor in West African lungfish, Protopterus annectens

We identified two ghrelin receptor isoforms, the ghrelin receptor type-1a (GHS-R1a) and its alternative splice form (GHS-R1b) for West African lungfish, Protopterus annectens. Lungfish GHS-R1a and 1b comprised 361 and 281 amino acids, respectively. Lungfish GHS-R1a showed the highest identity to coelacanth GHS-R1a (80.4%). The highest expression of GHS-R1a mRNAs was seen in the brain, liver, ovary, heart, intestine, and gills. GHS-R1b mRNAs were also detected in the same tissues with GHS-R1a, but their expression level was 1/20 that of GHS-R1a. In human embryonic kidney 293 cells transiently expressing lungfish GHS-R1a, rat and bullfrog ghrelin, and two GHS-R1a agonists, GHRP-6 and hexarelin, increased intracellular Ca(2+) concentrations. The intensity of the Ca(2+) increases induced by GHS-R1a agonists was twice when compared to that induced by ghrelin, although the median effective doses (ED50) were similar, suggesting a long-lasting effect of GHS-R1a agonists with similar affinity. We also examined changes in the GHS-R gene expression during an eight-week estivation. Body weight was slightly lowered, but plasma sodium and glucose concentrations decreased; plasma urea concentration increased significantly 4weeks after the start of estivation. Overall, expression of GHS-R1a mRNA decreased, but changes in GHS-R1b mRNA expression were inconsistent with those of GHS-R1a during estivation, suggesting an involvement of GHS-R in energy homeostasis, as seen in mammals. Our results suggest that the ghrelin-GHS-R1a system is present in this lungfish although ghrelin has not yet been found. The structure of GHS-R1a is closer to that of tetrapods than Actinopterygian fish, indicating a process of evolution that follows the Crossopterygii such as coelacanth.

Cloning, expression, and ligand-binding characterization of two neuropeptide Y receptor subtypes in orange-spotted grouper, Epinephelus coioides

As one of the most important multifunctional peptides, neuropeptide Y (NPY) performs its physiological functions through different subtype receptors. In this study, full-length cDNAs of two NPY receptors (YRs) in orange-spotted grouper (Epinephelus coioides) were cloned and named npy8br (y8b) and npy2r (y2). Phylogenetic analysis indicated that the Y8b receptor is an ortholog of the teleostean Y8b receptor, which belongs to the Y1 subfamily, and the Y2 receptor is an ortholog of the teleostean Y2 receptor, which belongs to the Y2 subfamily. Both of the YRs have G protein-coupled receptor family profiles. Multiple alignments demonstrate that the extracellular loop regions of YRs have distinctive residues of each species. Expression profile analysis revealed that the grouper Y8b receptor mRNA is primarily expressed in the brain, stomach and intestine, while the grouper Y2 receptor mRNA is primarily expressed in the brain, ovary, liver and heart. Double immunofluorescence analysis determined that the grouper YRs interact with the grouper NPY around the human embryonic kidney 293T cell surface. Furthermore, site-directed mutagenesis in a phage display system revealed that Asp(6.59) might be a common NPY-binding site, while Asp(2.68) of the Y8b receptor and Glu(5.24) of the Y2 receptor could be likely involved in subtype-specific binding. Combining the expression profile and ligand-binding feature, the grouper Y8b receptor could be involved in regulating food intake via the brain-gut axis and the grouper Y2 receptor might play a role in balancing the regulatory activity of the Y8b receptor and participate in metabolism in the liver and ovary.

New structural insights into the apelin receptor: identification of key residues for apelin binding

Apelin is the endogenous ligand of the orphan 7-transmembrane domain GPCR APJ, now named the apelin receptor (ApelinR). Apelin plays a prominent role in body fluid and cardiovascular homeostasis. To better understand the structural organization of the ApelinR, we built 3 homology 3-dimensional (3D) models of the human ApelinR using the validated cholecystokinin receptor-1 3D model or the X-ray structures of the β2-adrenergic and CXCR4 receptors as templates. Docking of the pyroglutamyl form of apelin 13 (pE13F) into these models revealed the conservation at the bottom of the binding site of a hydrophobic cavity in which the C-terminal Phe of pE13F was embedded. In contrast, at the top of the binding site, depending on the model, different interactions were visualized between acidic residues of the ApelinR and the basic residues of pE13F. Using site-directed mutagenesis, we showed that Asp 92, Glu 172, and Asp 282 of rat ApelinR are key residues in apelin binding by interacting with Lys 8, Arg 2, and Arg 4 of pE13F, respectively. These residues are only seen in the CXCR4-based ApelinR 3D model, further validating this model. These findings bring new insights into the structural organization of the ApelinR and the mode of apelin binding.

Structure-based virtual screening of the nociceptin receptor: hybrid docking and shape-based approaches for improved hit identification

The antagonist-bound crystal structure of the nociceptin receptor (NOP), from the opioid receptor family, was recently reported along with those of the other opioid receptors bound to opioid antagonists. We recently reported the first homology model of the ‘active-state’ of the NOP receptor, which when docked with ‘agonist’ ligands showed differences in the TM helices and residues, consistent with GPCR activation after agonist binding. In this study, we explored the use of the active-state NOP homology model for structure-based virtual screening to discover NOP ligands containing new chemical scaffolds. Several NOP agonist and antagonist ligands previously reported are based on a common piperidine scaffold. Given the structure–activity relationships for known NOP ligands, we developed a hybrid method that combines a structure-based and ligand-based approach, utilizing the active-state NOP receptor as well as the pharmacophoric features of known NOP ligands, to identify novel NOP binding scaffolds by virtual screening. Multiple conformations of the NOP active site including the flexible second extracellular loop (EL2) loop were generated by simulated annealing and ranked using enrichment factor (EF) analysis and a ligand–decoy dataset containing known NOP agonist ligands. The enrichment factors were further improved by combining shape-based screening of this ligand–decoy dataset and calculation of consensus scores. This combined structure-based and ligand-based EF analysis yielded higher enrichment factors than the individual methods, suggesting the effectiveness of the hybrid approach. Virtual screening of the CNS Permeable subset of the ZINC database was carried out using the above-mentioned hybrid approach in a tiered fashion utilizing a ligand pharmacophore-based filtering step, followed by structure-based virtual screening using the refined NOP active-state models from the enrichment analysis. Determination of the NOP receptor binding affinity of a selected set of top-scoring hits resulted in identification of several compounds with measurable binding affinity at the NOP receptor, one of which had a new chemotype for NOP receptor binding. The hybrid ligand-based and structure-based methodology demonstrates an effective approach for virtual screening that leverages existing SAR and receptor structure information for identifying novel hits for NOP receptor binding. The refined active-state NOP homology models obtained from the enrichment studies can be further used for structure-based optimization of these new chemotypes to obtain potent and selective NOP receptor ligands for therapeutic development.

A structural chemogenomics analysis of aminergic GPCRs: lessons for histamine receptor ligand design

BACKGROUND AND PURPOSE

Chemogenomics focuses on the discovery of new connections between chemical and biological space leading to the discovery of new protein targets and biologically active molecules. G-protein coupled receptors (GPCRs) are a particularly interesting protein family for chemogenomics studies because there is an overwhelming amount of ligand binding affinity data available. The increasing number of aminergic GPCR crystal structures now for the first time allows the integration of chemogenomics studies with high-resolution structural analyses of GPCR-ligand complexes.

EXPERIMENTAL APPROACH

In this study, we have combined ligand affinity data, receptor mutagenesis studies, and amino acid sequence analyses to high-resolution structural analyses of (hist)aminergic GPCR-ligand interactions. This integrated structural chemogenomics analysis is used to more accurately describe the molecular and structural determinants of ligand affinity and selectivity in different key binding regions of the crystallized aminergic GPCRs, and histamine receptors in particular.

KEY RESULTS

Our investigations highlight interesting correlations and differences between ligand similarity and ligand binding site similarity of different aminergic receptors. Apparent discrepancies can be explained by combining detailed analysis of crystallized or predicted protein-ligand binding modes, receptor mutation studies, and ligand structure-selectivity relationships that identify local differences in essential pharmacophore features in the ligand binding sites of different receptors.

CONCLUSIONS AND IMPLICATIONS

We have performed structural chemogenomics studies that identify links between (hist)aminergic receptor ligands and their binding sites and binding modes. This knowledge can be used to identify structure-selectivity relationships that increase our understanding of ligand binding to (hist)aminergic receptors and hence can be used in future GPCR ligand discovery and design.

Allosteric and biased g protein-coupled receptor signaling regulation: potentials for new therapeutics

G protein-coupled receptors (GPCRs) are seven-transmembrane proteins that participate in many aspects of the endocrine function and are important targets for drug development. They transduce signals mainly, but not exclusively, via hetero-trimeric G proteins, leading to a diversity of intracellular signaling cascades. Ligands binding at the hormone orthosteric sites of receptors have been classified as agonists, antagonists, and/or inverse agonists based on their ability to mainly modulate G protein signaling. Accumulating evidence also indicates that such ligands, alone or in combination with other ones such as those acting outside the orthosteric hormone binding sites (e.g., allosteric modulators), have the ability to selectively engage subsets of signaling responses as compared to the natural endogenous ligands. Such modes of functioning have been variously referred to as "functional selectivity" or "ligand-biased signaling." In this review, we provide an overview of the current knowledge regarding GPCR-biased signaling and their functional regulation with a focus on the evolving concept that receptor domains can also be targeted to allosterically bias signaling, and discuss the usefulness of such modes of regulation for the design of more efficient therapeutics.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Recent structural advances of β1 and β2 adrenoceptors yield keys for ligand recognition and drug design

Because they represent attractive drug targets, adrenoceptors have been widely studied. Recent progress in structural data of β-adrenoceptors allows us to understand and predict key interactions in ligand recognition and receptor activation. Nevertheless, an important aspect of this process has only begun to be explored: the stabilization of a conformational state of these receptors upon contact with a ligand and the capacity of a ligand to influence receptor conformation through allosteric modulation, biased signaling, and selectivity. The aim of the present Perspective is to identify the well-defined orthosteric binding site and possible allosteric sites and to analyze the importance of the ligand-receptor interaction in the stabilization of certain receptor conformations. For this purpose, we have reviewed recent advances made through the use of X-ray data from ligand-β-adrenoceptor (including ADRB1 and ADRB2) crystal structures. Most importantly, implications in the medicinal chemistry field are explored in relation to drug design.

Global fold of human cannabinoid type 2 receptor probed by solid-state 13C-, 15N-MAS NMR and molecular dynamics simulations

The global fold of human cannabinoid type 2 (CB2 ) receptor in the agonist-bound active state in lipid bilayers was investigated by solid-state (13)C- and (15)N magic-angle spinning (MAS) NMR, in combination with chemical-shift prediction from a structural model of the receptor obtained by microsecond-long molecular dynamics (MD) simulations. Uniformly (13)C- and (15)N-labeled CB2 receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. (13)C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of Cα, Cβ, and C=O bands of superimposed resonances with predictions from protein structures generated by MD. The experimental NMR spectra matched the calculated spectra reasonably well indicating agreement of the global fold of the protein between experiment and simulations. In particular, the (13) C chemical shift distribution of Cα resonances was shown to be very sensitive to both the primary amino acid sequence and the secondary structure of CB2. Thus the shape of the Cα band can be used as an indicator of CB2 global fold. The prediction from MD simulations indicated that upon receptor activation a rather limited number of amino acid residues, mainly located in the extracellular Loop 2 and the second half of intracellular Loop 3, change their chemical shifts significantly (≥ 1.5 ppm for carbons and ≥ 5.0 ppm for nitrogens). Simulated two-dimensional (13) Cα(i)-(13)C=O(i) and (13)C=O(i)-(15)NH(i + 1) dipolar-interaction correlation spectra provide guidance for selective amino acid labeling and signal assignment schemes to study the molecular mechanism of activation of CB2 by solid-state MAS NMR.