An online resource for GPCR structure determination and analysis.

G-protein-coupled receptors (GPCRs) transduce physiological and sensory stimuli into appropriate cellular responses and mediate the actions of one-third of drugs. GPCR structural studies have revealed the general bases of receptor activation, signaling, drug action and allosteric modulation, but so far cover only 13% of nonolfactory receptors. We broadly surveyed the receptor modifications/engineering and methods used to produce all available GPCR crystal and cryo-electron microscopy (cryo-EM) structures, and present an interactive resource integrated in GPCRdb ( http://www.gpcrdb.org ) to assist users in designing constructs and browsing appropriate experimental conditions for structure studies.

Identification of Histamine H3 Receptor Ligands Using a New Crystal Structure Fragment-based Method.

Virtual screening offers an efficient alternative to high-throughput screening in the identification of pharmacological tools and lead compounds. Virtual screening is typically based on the matching of target structures or ligand pharmacophores to commercial or in-house compound catalogues. This study provides the first proof-of-concept for our recently reported method where pharmacophores are instead constructed based on the inference of residue-ligand fragments from crystal structures. We demonstrate its unique utility for G protein-coupled receptors, which represent the largest families of human membrane proteins and drug targets. We identified five neutral antagonists and one inverse agonist for the histamine H3 receptor with potencies of 0.7-8.5 µM in a recombinant receptor cell-based inositol phosphate accumulation assay and validated their activity using a radioligand competition binding assay. H3 receptor antagonism is of large therapeutic value and our ligands could serve as starting points for further lead optimisation. The six ligands exhibit four chemical scaffolds, whereof three have high novelty in comparison to the known H3 receptor ligands in the ChEMBL database. The complete pharmacophore fragment library is freely available through the GPCR database, GPCRdb, allowing the successful application herein to be repeated for most of the 285 class A GPCR targets. The method could also easily be adapted to other protein families.

The GPR139 reference agonists 1a and 7c, and tryptophan and phenylalanine share a common binding site.

GPR139 is an orphan G protein-coupled receptor expressed in the brain, in particular in the habenula, hypothalamus and striatum. It has therefore been suggested that GPR139 is a possible target for metabolic disorders and Parkinson's disease. Several surrogate agonist series have been published for GPR139. Two series published by Shi et al. and Dvorak et al. included agonists 1a and 7c respectively, with potencies in the ten-nanomolar range. Furthermore, Isberg et al. and Liu et al. have previously shown that tryptophan (Trp) and phenylalanine (Phe) can activate GPR139 in the hundred-micromolar range. In this study, we produced a mutagenesis-guided model of the GPR139 binding site to form a foundation for future structure-based ligand optimization. Receptor mutants studied in a Ca2+ assay demonstrated that residues F1093×33, H1875×43, W2416×48 and N2717×38, but not E1083×32, are highly important for the activation of GPR139 as predicted by the receptor model. The initial ligand-receptor complex was optimized through free energy perturbation simulations, generating a refined GPR139 model in agreement with experimental data. In summary, the GPR139 reference surrogate agonists 1a and 7c, and the endogenous amino acids L-Trp and L-Phe share a common binding site, as demonstrated by mutagenesis, ligand docking and free energy calculations.

The orphan G protein-coupled receptor GPR139 is activated by the peptides: Adrenocorticotropic hormone (ACTH), α-, and ß-melanocyte stimulating hormone (α-MSH, and ß-MSH), and the conserved core motif HFRW.

GPR139 is an orphan G protein-coupled receptor that is expressed primarily in the brain. Not much is known regarding the function of GPR139. Recently we have shown that GPR139 is activated by the amino acids l-tryptophan and l-phenylalanine (EC50 values of 220 µM and 320 µM, respectively), as well as di-peptides comprised of aromatic amino acids. This led us to hypothesize that GPR139 may be activated by peptides. Sequence alignment of the binding cavities of all class A GPCRs, revealed that the binding pocket of the melanocortin 4 receptor is similar to that of GPR139. Based on the chemogenomics principle "similar targets bind similar ligands", we tested three known endogenous melanocortin 4 receptor agonists; adrenocorticotropic hormone (ACTH) and α- and ß-melanocyte stimulating hormone (α-MSH and ß-MSH) on CHO-k1 cells stably expressing the human GPR139 in a Fluo-4 Ca2+-assay. All three peptides, as well as their conserved core motif HFRW, were found to activate GPR139 in the low micromolar range. Moreover, we found that peptides consisting of nine or ten N-terminal residues of α-MSH activate GPR139 in the submicromolar range. α-MSH1-9 was found to correspond to the product of a predicted cleavage site in the pre-pro-protein pro-opiomelanocortin (POMC). Our results demonstrate that GPR139 is a peptide receptor, activated by ACTH, α-MSH, ß-MSH, the conserved core motif HFRW as well as a potential endogenous peptide α-MSH1-9. Further studies are needed to determine the functional relevance of GPR139 mediated signaling by these peptides.

Novel Agonist Bioisosteres and Common Structure-Activity Relationships for The Orphan G Protein-Coupled Receptor GPR139.

GPR139 is an orphan class A G protein-coupled receptor found mainly in the central nervous system. It has its highest expression levels in the hypothalamus and striatum, regions regulating metabolism and locomotion, respectively, and has therefore been suggested as a potential target for obesity and Parkinson's disease. The two aromatic amino acids L-Trp and L-Phe have been proposed as putative endogenous agonists, and three structurally related benzohydrazide, glycine benzamide, and benzotriazine surrogate agonist series have been published. Herein, we assayed 158 new analogues selected from a pharmacophore model, and identified 12 new GPR139 agonists, containing previously untested bioisosteres. Furthermore, we present the first combined structure-activity relationships, and a refined pharmacophore model to serve as a rationale for future ligand identification and optimization.

Integrating structural and mutagenesis data to elucidate GPCR ligand binding.

G protein-coupled receptors (GPCRs) represent the largest family of human membrane proteins, as well as drug targets. A recent boom in GPCR structural biology has provided detailed images of receptor ligand binding sites and interactions on the molecular level. An ever-increasing number of ligands is reported that exhibit activity through multiple receptors, binding in allosteric sites, and bias towards different intracellular signalling pathways. Furthermore, a wealth of single point mutants has accumulated in literature and public databases. Integrating these structural and mutagenesis data will help elucidate new GPCR ligand binding sites, and ultimately design drugs with tailored pharmacological activity.

GPCRdb: an information system for G protein-coupled receptors.

Recent developments in G protein-coupled receptor (GPCR) structural biology and pharmacology have greatly enhanced our knowledge of receptor structure-function relations, and have helped improve the scientific foundation for drug design studies. The GPCR database, GPCRdb, serves a dual role in disseminating and enabling new scientific developments by providing reference data, analysis tools and interactive diagrams. This paper highlights new features in the fifth major GPCRdb release: (i) GPCR crystal structure browsing, superposition and display of ligand interactions; (ii) direct deposition by users of point mutations and their effects on ligand binding; (iii) refined snake and helix box residue diagram looks; and (iii) phylogenetic trees with receptor classification colour schemes. Under the hood, the entire GPCRdb front- and back-ends have been re-coded within one infrastructure, ensuring a smooth browsing experience and development. GPCRdb is available at http://www.gpcrdb.org/ and it's open source code at https://bitbucket.org/gpcr/protwis.

Selective Negative Allosteric Modulation Of Metabotropic Glutamate Receptors ­ A Structural Perspective of Ligands and Mutants.

The metabotropic glutamate receptors have a wide range of modulatory functions in the central nervous system. They are among the most highly pursued drug targets, with relevance for several neurological diseases, and a number of allosteric modulators have entered clinical trials. However, so far this has not led to a marketed drug, largely because of the difficulties in achieving subtype-selective compounds with desired properties. Very recently the first crystal structures were published for the transmembrane domain of two metabotropic glutamate receptors in complex with negative allosteric modulators. In this analysis, we make the first comprehensive structural comparison of all metabotropic glutamate receptors, placing selective negative allosteric modulators and critical mutants into the detailed context of the receptor binding sites. A better understanding of how the different mGlu allosteric modulator binding modes relates to selective pharmacological actions will be very valuable for rational design of safer drugs.

Creating a specialist protein resource network: a meeting report for the protein bioinformatics and community resources retreat.

During 11-12 August 2014, a Protein Bioinformatics and Community Resources Retreat was held at the Wellcome Trust Genome Campus in Hinxton, UK. This meeting brought together the principal investigators of several specialized protein resources (such as CAZy, TCDB and MEROPS) as well as those from protein databases from the large Bioinformatics centres (including UniProt and RefSeq). The retreat was divided into five sessions: (1) key challenges, (2) the databases represented, (3) best practices for maintenance and curation, (4) information flow to and from large data centers and (5) communication and funding. An important outcome of this meeting was the creation of a Specialist Protein Resource Network that we believe will improve coordination of the activities of its member resources. We invite further protein database resources to join the network and continue the dialogue.

A new crystal structure fragment-based pharmacophore method for G protein-coupled receptors.

We have developed a new method for the building of pharmacophores for G protein-coupled receptors, a major drug target family. The method is a combination of the ligand- and target-based pharmacophore methods and founded on the extraction of structural fragments, interacting ligand moiety and receptor residue pairs, from crystal structure complexes. We describe the procedure to collect a library with more than 250 fragments covering 29 residue positions within the generic transmembrane binding pocket. We describe how the library fragments are recombined and inferred to build pharmacophores for new targets. A validating retrospective virtual screening of histamine H1 and H3 receptor pharmacophores yielded area-under-the-curves of 0.88 and 0.82, respectively. The fragment-based method has the unique advantage that it can be applied to targets for which no (homologous) crystal structures or ligands are known. 47% of the class A G protein-coupled receptors can be targeted with at least four-element pharmacophores. The fragment libraries can also be used to grow known ligands or for rotamer refinement of homology models. Researchers can download the complete fragment library or a subset matching their receptor of interest using our new tool in GPCRDB.

Generic GPCR residue numbers - aligning topology maps while minding the gaps.

Generic residue numbers facilitate comparisons of, for example, mutational effects, ligand interactions, and structural motifs. The numbering scheme by Ballesteros and Weinstein for residues within the class A GPCRs (G protein-coupled receptors) has more than 1100 citations, and the recent crystal structures for classes B, C, and F now call for a community consensus in residue numbering within and across these classes. Furthermore, the structural era has uncovered helix bulges and constrictions that offset the generic residue numbers. The use of generic residue numbers depends on convenient access by pharmacologists, chemists, and structural biologists. We review the generic residue numbering schemes for each GPCR class, as well as a complementary structure-based scheme, and provide illustrative examples and GPCR database (GPCRDB) web tools to number any receptor sequence or structure.

Computer-aided discovery of aromatic l-α-amino acids as agonists of the orphan G protein-coupled receptor GPR139.

GPR139 is an orphan G protein-coupled receptor expressed mainly in the central nervous system. We developed a pharmacophore model based on known GPR139 surrogate agonists which led us to propose aromatic-containing dipeptides as potential ligands. Upon testing, the dipeptides demonstrated agonism in the Gq pathway. Next, in testing all 20 proteinogenic l-α-amino acids, L-tryptophan and l-phenylalanine were found to have EC50 values of 220 and 320 µM, respectively, making them the first putative endogenous agonists for GPR139.

GPCRDB: an information system for G protein-coupled receptors.

For the past 20 years, the GPCRDB (G protein-coupled receptors database; http://www.gpcr.org/7tm/) has been a 'one-stop shop' for G protein-coupled receptor (GPCR)-related data. The GPCRDB contains experimental data on sequences, ligand-binding constants, mutations and oligomers, as well as many different types of computationally derived data, such as multiple sequence alignments and homology models. The GPCRDB also provides visualization and analysis tools, plus a number of query systems. In the latest GPCRDB release, all multiple sequence alignments, and >65,000 homology models, have been significantly improved, thanks to a recent flurry of GPCR X-ray structure data. Tools were introduced to browse X-ray structures, compare binding sites, profile similar receptors and generate amino acid conservation statistics. Snake plots and helix box diagrams can now be custom coloured (e.g. by chemical properties or mutation data) and saved as figures. A series of sequence alignment visualization tools has been added, and sequence alignments can now be created for subsets of sequences and sequence positions, and alignment statistics can be produced for any of these subsets.

Structure-activity relationships of constrained phenylethylamine ligands for the serotonin 5-HT2 receptors.

Serotonergic ligands have proven effective drugs in the treatment of migraine, pain, obesity, and a wide range of psychiatric and neurological disorders. There is a clinical need for more highly 5-HT2 receptor subtype-selective ligands and the most attention has been given to the phenethylamine class. Conformationally constrained phenethylamine analogs have demonstrated that for optimal activity the free lone pair electrons of the 2-oxygen must be oriented syn and the 5-oxygen lone pairs anti relative to the ethylamine moiety. Also the ethyl linker has been constrained providing information about the bioactive conformation of the amine functionality. However, combined 1,2-constriction by cyclization has only been tested with one compound. Here, we present three new 1,2-cyclized phenylethylamines, 9-11, and describe their synthetic routes. Ligand docking in the 5-HT2B crystal structure showed that the 1,2-heterocyclized compounds can be accommodated in the binding site. Conformational analysis showed that 11 can only bind in a higher-energy conformation, which would explain its absent or low affinity. The amine and 2-oxygen interactions with D3.32 and S3.36, respectively, can form but shift the placement of the core scaffold. The constraints in 9-11 resulted in docking poses with the 4-bromine in closer vicinity to 5.46, which is polar only in the human 5-HT2A subtype, for which 9-11 have the lowest affinity. The new ligands, conformational analysis and docking expand the structure-activity relationships of constrained phenethylamines and contributes towards the development of 5-HT2 receptor subtype-selective ligands.

Structure-activity relationships and identification of optmized CC-chemokine receptor CCR1, 5, and 8 metal-ion chelators.

Chemokine receptors are involved in trafficking of leukocytes and represent targets for autoimmune conditions, inflammatory diseases, viral infections, and cancer. We recently published CCR1, CCR8, and CCR5 agonists and positive modulators based on a three metal-ion chelator series: 2,2'-bipyridine, 1,10-phenanthroline, and 2,2';6',2″-terpyridine. Here, we have performed an in-depth structure-activity relationship study and tested eight new optimized analogs. Using density functional theory calculations we demonstrate that the chelator zinc affinities depend on how electron-donating and -withdrawing substituents modulate the partial charges of chelating nitrogens. The zinc affinity was found to constitute the major factor for receptor potency, although the activity of some chelators deviate suggesting favorable or unfavorable interactions. Hydrophobic and halogen substituents are generally better accommodated in the receptors than polar groups. The new analog brominated terpyridine (29) resulted in the highest chelator potencies observed so far CCR1 (EC50: 0.49 µM) and CCR8 (EC50: 0.28 µM). Furthermore, we identified the first selective CCR5 agonist chelator, meta dithiomethylated bipyridine (23). The structure-activity relationships contribute to small-molecule drug development, and the novel chelators constitute valuable tools for studies of structural mechanisms for chemokine receptor activation.

Design, synthesis, and pharmacological characterization of N- and O-substituted 5,6,7,8-tetrahydro-4H-isoxazolo[4,5-d]azepin-3-ol analogues: novel 5-HT(2A)/5-HT(2C) receptor agonists with pro-cognitive properties.

The isoxazol-3-one tautomer of the bicyclic isoxazole, 5,6,7,8-tetrahydro-4H-isoxazolo[4,5-d]azepin-3-ol (THAZ), has previously been shown to be a weak GABA(A) and glycine receptor antagonist. In the present study, the potential in this scaffold has been explored through the synthesis and pharmacological characterization of a series of N- and O-substituted THAZ analogues. The analogues N-Bn-THAZ (3d) and O-Bn-THAZ (4d) were found to be potent agonists of the human 5-HT(2A) and 5-HT(2C) receptors. Judging from an elaborate pharmacological profiling at numerous other CNS targets, the 3d analogue appears to be selective for the two receptors. Administration of 3d substantially improved the cognitive performance of mice in a place recognition Y-maze model, an effect fully reversible by coadministration of the selective 5-HT(2C) antagonist SB242084. In conclusion, as novel bioavailable cognitive enhancers that most likely mediate their effects through 5-HT(2A) and/or 5-HT(2C) receptors, the isoxazoles 3d and 4d constitute interesting leads for further medicinal chemistry development.

G protein- and agonist-bound serotonin 5-HT2A receptor model activated by steered molecular dynamics simulations.

A 5-HT(2A) receptor model was constructed by homology modeling based on the ß(2)-adrenergic receptor and the G protein-bound opsin crystal structures. The 5-HT(2A) receptor model was transferred into an active conformation by an agonist ligand and a G(αq) peptide in four subsequent steered molecular dynamics (MD) simulations. The driving force for the transformation was the addition of several known intermolecular and receptor interhelical hydrogen bonds enforcing the necessary helical and rotameric movements. Subsquent MD simulations without constraints confirmed the stability of the activated receptor model as well as revealed new information about stabilizing residues and bonds. The active 5-HT(2A) receptor model was further validated by retrospective ligand screening of more than 9400 compounds, whereof 182 were known ligands. The results show that the model can be used in drug discovery for virtual screening and structure-based ligand design as well as in GPCR activation studies.