A chemogenomic analysis of the transmembrane binding cavity of human G-protein-coupled receptors

Sequence, structure and ligand binding evolution of rhodopsin-like G protein-coupled receptors: a crystal structure-based phylogenetic analysis

G protein-coupled receptors (GPCRs) form the largest family of membrane receptors in the human genome. Advances in membrane protein crystallization so far resulted in the determination of 24 receptors available as high-resolution atomic structures. We performed the first phylogenetic analysis of GPCRs based on the available set of GPCR structures. We present a new phylogenetic tree of known human rhodopsin-like GPCR sequences based on this structure set. We can distinguish the three separate classes of small-ligand binding GPCRs, peptide binding GPCRs, and olfactory receptors. Analyzing different structural subdomains, we found that small molecule binding receptors most likely have evolved from peptide receptor precursors, with a rhodopsin/S1PR1 ancestor, most likely an ancestral opsin, forming the link between both classes. A light-activated receptor therefore seems to be the origin of the small molecule hormone receptors of the central nervous system. We find hints for a common evolutionary path of both ligand binding site and central sodium/water binding site. Surprisingly, opioid receptors exhibit both a binding cavity and a central sodium/water binding site similar to the one of biogenic amine receptors instead of peptide receptors, making them seemingly prone to bind small molecule ligands, e.g. opiates. Our results give new insights into the relationship and the pharmacological properties of rhodopsin-like GPCRs.

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.

D1-like receptors distinguishing thieno-azecine regioisomers

Design of novel azecine derivatives with modulated dopaminergic receptor selectivity and affinity profiles.

Polypharmacology modelling using proteochemometrics (PCM): recent methodological developments, applications to target families, and future prospects

Proteochemometric (PCM) modelling is a computational method to model the bioactivity of multiple ligands against multiple related protein targets simultaneously.

Fishing for accessory binding sites at GPCRs with 'loop-hooks' - an approach towards selectivity? Part I

Receptor-subtype selectivity is an important issue in medicinal chemistry and can become very difficult to achieve if the actual binding pockets of the respective receptors are highly conserved. For such cases, known unselective ligands could be equipped with a spacer that sticks outside the actual orthosteric binding pocket towards the extracellular loops. The end of the spacer bears certain functional groups to enable specific or unspecific interactions with the receptor residues outside the binding cavity. Our experiments indicated that it is possible to achieve selectivity within the dopamine D1 family with such 'loop-hooks'.

What can we learn from molecular dynamics simulations for GPCR drug design?

Recent years have seen a tremendous progress in the elucidation of experimental structural information for G-protein coupled receptors (GPCRs). Although for the vast majority of pharmaceutically relevant GPCRs structural information is still accessible only by homology models the steadily increasing amount of structural information fosters the application of structure-based drug design tools for this important class of drug targets. In this article we focus on the application of molecular dynamics (MD) simulations in GPCR drug discovery programs. Typical application scenarios of MD simulations and their scope and limitations will be described on the basis of two selected case studies, namely the binding of small molecule antagonists to the human CC chemokine receptor 3 (CCR3) and a detailed investigation of the interplay between receptor dynamics and solvation for the binding of small molecules to the human muscarinic acetylcholine receptor 3 (hM3R).

Insights into the binding modes of CC chemokine receptor 4 (CCR4) inhibitors: a combined approach involving homology modelling, docking, and molecular dynamics simulation studies

CC chemokine receptor 4 (CCR4), a G protein-coupled receptor (GPCR), plays a vital role in the progression of asthma, T-cell lymphoma, inflammation, and Alzheimer's disease. To date, the structure of CCR4 has not been determined. Therefore, the nature of the interactions between inhibitors and CCR4 is not well known. In this study, we used CCR5 as a template to model the structure of CCR4. Docking studies were performed for four naphthalene-sulphonamide derivatives and crucial ligand-protein interactions were analysed. Molecular dynamics (MD) simulations of these complexes (100 ns each) were carried out to gain insights into the interactions between ligands and CCR4. MD simulations revealed that the residues identified by the docking were displaced and new residues were inserted near the ligands. Results of a principal component analysis (PCA) suggested that CCR4 unfolds at the extracellular site surrounding the ligands. Our simulations identified crucial residues involved in CCR4 antagonism, which were supported by previous mutational studies. Additionally, we identified Ser3.29, Leu3.33, Ser5.39, Phe6.47, Ile7.35, Thr7.38, Thr7.40, and Ala7.42 as residues that play crucial roles in CCR4 antagonism. Mutational studies will help elucidate the significance of these residues in CCR4 antagonism. An understanding of ligand-CCR4 interactions might aid in the design of novel CCR4 inhibitors.

Evolutionary history of the neuropeptide S receptor/neuropeptide S system

The neuropeptide S receptor (NPSR) belongs to the G protein-coupled receptor (GPCR) superfamily and is activated by the neuropeptide S (NPS). Although recently discovered, the vertebrate NPSR-NPS system has been established as an important signaling system in the central nervous system and is involved in physiological processes such as locomotor activity, wakefulness, asthma pathogenesis, anxiety and food intake. The availability of a large number of genome sequences from multiple bilaterian lineages has provided an opportunity to establish the evolutionary history of the system. This review describes the origin and the molecular evolution of the NPSR-NPS system using data derived primarily from comparative genomic analyses. These analyses indicate that the NPSR-NPS system and the vasopressin-like receptor-vasopressin/oxytocin peptide (VPR-VP/OT) system originated from a single system in an ancestral bilaterian. Multiple duplications of this ancestral system gave rise to the bilaterian VPR-VP/OT system and to the protostomian cardioacceleratory peptide receptor-cardioacceleratory peptide (CCAPR-CCAP) system and to the NPSR-NPS system in the deuterostomes. Gene structure features of the receptors were consistent with the orthology annotations derived from phylogenetic analyses. The orthology of the peptide precursors closely paralleled that of the receptors suggesting an ancient coevolution of the receptor-peptide pair. An important challenge for the coevolution hypothesis will be to establish the molecular and structural basis of the divergence between orthologous receptor-ligand pairs in this system.

Computational studies to predict or explain G protein coupled receptor polypharmacology

Since G protein-coupled receptors (GPCRs) belong to a very large superfamily of evolutionarily related receptors (>800 members in humans), and due to the rapid progress on their structural biology, they are ideal candidates for polypharmacology studies. Broad screening and bioinformatics/chemoinformatics have been applied to understanding off-target effects of GPCR ligands. It is now feasible to approach the question of GPCR polypharmacology using molecular modeling and the available X-ray GPCR structures. As an example, large and sterically constrained adenosine derivatives (potent adenosine receptor ligands with low conformational freedom and multiple extended substituents) were screened for binding at diverse receptors. Unanticipated off-target interactions, including at biogenic amine receptors, were then modeled using a structure-based approach to provide a consistent understanding of recognition. A conserved Asp in TM3 changed its role from counterion for biogenic amines to characteristic H-bonding to adenosine. The same systematic approach could potentially be applied to many GPCRs or other receptors using other sets of congeneric ligands.

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.

Virtual screening: an in silico tool for interlacing the chemical universe with the proteome

In silico screening both in the forward (traditional virtual screening) and reverse sense (inverse virtual screening (IVS)) are helpful techniques for interlacing the chemical universe of small molecules with the proteome. The former, which is using a protein structure and a large chemical database, is well-known by the scientific community. We have chosen here to provide an overview on the latter, focusing on validation and target prioritization strategies. By comparing it to complementary or alternative wet-lab approaches, we put IVS in the broader context of chemical genomics, target discovery and drug design. By giving examples from the literature and an own example on how to validate the approach, we provide guidance on the issues related to IVS.

C(X)CR in silico: Computer-aided prediction of chemokine receptor-ligand interactions

This review will focus on the construction, refinement, and validation of chemokine receptor models for the purpose of structure-based virtual screening and ligand design. The review will present a comparative analysis of ligand binding pockets in chemokine receptors, including a review of the recently released CXCR4 X-ray structures, and their implication on chemokine receptor (homology) modeling. The recommended protein-ligand modeling procedure as well as the use of experimental anchors to steer the modeling procedure is discussed and an overview of several successful structure-based ligand discovery and design studies is provided. This review shows that receptor models, despite structural inaccuracies, can be efficiently used to find novel ligands for chemokine receptors.:

Sequence-structure based phylogeny of GPCR Class A Rhodopsin receptors

Current methods of G protein coupled receptors (GPCRs) phylogenetic classification are sequence based and therefore inappropriate for highly divergent sequences, sharing low sequence identity. In this study, sequence structure profile based alignment generated by PROMALS3D was used to understand the GPCR Class A Rhodopsin superfamily evolution using the MEGA 5 software. Phylogenetic analysis included a combination of Neighbor-Joining method and Maximum Likelihood method, with 1000 bootstrap replicates. Our study was able to identify potential ligand association for Class A Orphans and putative/unclassified Class A receptors with no cognate ligand information: GPR21 and GPR52 with fatty acids; GPR75 with Neuropeptide Y; GPR82, GPR18, GPR141 with N-arachidonylglycine; GPR176 with Free fatty acids, GPR10 with Tachykinin & Neuropeptide Y; GPR85 with ATP, ADP & UDP glucose; GPR151 with Galanin; GPR153 and GPR162 with Adrenalin, Noradrenalin; GPR146, GPR139, GPR142 with Neuromedin, Ghrelin, Neuromedin U-25 & Thyrotropin-releasing hormone; GPR171 with ATP, ADP & UDP Glucose; GPR88, GPR135, GPR161, GPR101with 11-cis-retinal; GPR83 with Tackykinin; GPR148 with Prostanoids, GPR109b, GPR81, GPR31with ATP & UTP and GPR150 with GnRH I & GnRHII. Furthermore, we suggest that this study would prove useful in re-classification of receptors, selecting templates for homology modeling and identifying ligands which may show cross reactivity with other GPCRs as signaling via multiple ligands play a significant role in disease modulation.

Probing the Binding Site of Bile Acids in TGR5

TGR5 is a G-protein-coupled receptor (GPCR) mediating cellular responses to bile acids (BAs). Although some efforts have been devoted to generate homology models of TGR5 and draw structure-activity relationships of BAs, none of these studies has hitherto described how BAs bind to TGR5. Here, we present an integrated computational, chemical, and biological approach that has been instrumental to determine the binding mode of BAs to TGR5. As a result, key residues have been identified that are involved in mediating the binding of BAs to the receptor. Collectively, these results provide new hints to design potent and selective TGR5 agonists.

Mutagenesis and computational modeling of human G-protein-coupled receptor Y2 for neuropeptide Y and peptide YY

Neuropeptide Y and peptide YY receptor type 2 (Y2) is involved in appetite regulation and several other physiological processes. We have investigated the structure of the human Y2 receptor. Computational modeling of receptor-agonist interactions was used as a guide to design a series of receptor mutants, followed by binding assays using full-length and truncated peptide agonists and the Y2-specific antagonist BIIE0246. Our model suggested a hydrogen bond network among highly conserved residues Thr2.61, Gln3.32, and His7.39, which could play roles in ligand binding and/or receptor structure. In addition, the C-terminus of the peptide could make contact with residues Tyr5.38 and Leu6.51. Mutagenesis of all these positions, followed by binding assays, provides experimental support for our computational model: most of the mutants for the residues forming the proposed hydrogen bond network displayed reduced peptide agonist affinities as well as reduced hPYY3-36 potency in a functional assay. The Ala and Leu mutants of Gln3.32 and His7.39 disrupted membrane expression of the receptor. Combined with the modeling, the experimental results support roles for these hydrogen bond network residues in peptide binding as well as receptor architecture. The reduced agonist affinity for mutants of Tyr5.38 and Leu6.51 supports their role in a binding pocket surrounding the invariant tyrosine at position 36 of the peptide ligands. The results for antagonist BIIE0246 suggest several differences in interactions compared to those of the peptides. Our results lead to a new structural model for NPY family receptors and peptide binding.

Novel insights on thyroid-stimulating hormone receptor signal transduction

The TSH receptor (TSHR) is a member of the glycoprotein hormone receptors, a subfamily of family A G protein-coupled receptors. The TSHR is of great importance for the growth and function of the thyroid gland. The TSHR and its endogenous ligand TSH are pivotal proteins with respect to a variety of physiological functions and malfunctions. The molecular events of TSHR regulation can be summarized as a process of signal transduction, including signal reception, conversion, and amplification. The steps during signal transduction from the extra- to the intracellular sites of the cell are not yet comprehensively understood. However, essential new insights have been achieved in recent years on the interrelated mechanisms at the extracellular region, the transmembrane domain, and intracellular components. This review contains a critical summary of available knowledge of the molecular mechanisms of signal transduction at the TSHR, for example, the key amino acids involved in hormone binding or in the structural conformational changes that lead to G protein activation or signaling regulation. Aspects of TSHR oligomerization, signaling promiscuity, signaling selectivity, phenotypes of genetic variations, and potential extrathyroidal receptor activity are also considered, because these are relevant to an understanding of the overall function of the TSHR, including physiological, pathophysiological, and pharmacological perspectives. Directions for future research are discussed.

Benchmarking of protein descriptor sets in proteochemometric modeling (part 2): modeling performance of 13 amino acid descriptor sets

Background

While a large body of work exists on comparing and benchmarking descriptors of molecular structures, a similar comparison of protein descriptor sets is lacking. Hence, in the current work a total of 13 amino acid descriptor sets have been benchmarked with respect to their ability of establishing bioactivity models. The descriptor sets included in the study are Z-scales (3 variants), VHSE, T-scales, ST-scales, MS-WHIM, FASGAI, BLOSUM, a novel protein descriptor set (termed ProtFP (4 variants)), and in addition we created and benchmarked three pairs of descriptor combinations. Prediction performance was evaluated in seven structure-activity benchmarks which comprise Angiotensin Converting Enzyme (ACE) dipeptidic inhibitor data, and three proteochemometric data sets, namely (1) GPCR ligands modeled against a GPCR panel, (2) enzyme inhibitors (NNRTIs) with associated bioactivities against a set of HIV enzyme mutants, and (3) enzyme inhibitors (PIs) with associated bioactivities on a large set of HIV enzyme mutants.

Results

The amino acid descriptor sets compared here show similar performance (<0.1 log units RMSE difference and <0.1 difference in MCC), while errors for individual proteins were in some cases found to be larger than those resulting from descriptor set differences ( > 0.3 log units RMSE difference and >0.7 difference in MCC). Combining different descriptor sets generally leads to better modeling performance than utilizing individual sets. The best performers were Z-scales (3) combined with ProtFP (Feature), or Z-Scales (3) combined with an average Z-Scale value for each target, while ProtFP (PCA8), ST-Scales, and ProtFP (Feature) rank last.

Conclusions

While amino acid descriptor sets capture different aspects of amino acids their ability to be used for bioactivity modeling is still – on average – surprisingly similar. Still, combining sets describing complementary information consistently leads to small but consistent improvement in modeling performance (average MCC 0.01 better, average RMSE 0.01 log units lower). Finally, performance differences exist between the targets compared thereby underlining that choosing an appropriate descriptor set is of fundamental for bioactivity modeling, both from the ligand- as well as the protein side.

Design and pharmacological characterization of VUF14480, a covalent partial agonist that interacts with cysteine 98(3.36) of the human histamine H₄ receptor

BACKGROUND AND PURPOSE

The recently proposed binding mode of 2-aminopyrimidines to the human (h) histamine H₄ receptor suggests that the 2-amino group of these ligands interacts with glutamic acid residue E182(5.46) in the transmembrane (TM) helix 5 of this receptor. Interestingly, substituents at the 2-position of this pyrimidine are also in close proximity to the cysteine residue C98(3.36) in TM3. We hypothesized that an ethenyl group at this position will form a covalent bond with C98(3.36) by functioning as a Michael acceptor. A covalent pyrimidine analogue will not only prove this proposed binding mode, but will also provide a valuable tool for H4 receptor research.

EXPERIMENTAL APPROACH

We designed and synthesized VUF14480, and pharmacologically characterized this compound in hH4 receptor radioligand binding, G protein activation and β-arrestin2 recruitment experiments. The ability of VUF14480 to act as a covalent binder was assessed both chemically and pharmacologically.

KEY RESULTS

VUF14480 was shown to be a partial agonist of hH4 receptor-mediated G protein signalling and β-arrestin2 recruitment. VUF14480 bound covalently to the hH₄ receptor with submicromolar affinity. Serine substitution of C98(3.36) prevented this covalent interaction.

CONCLUSION AND IMPLICATIONS

VUF14480 is thought to bind covalently to the hH₄ receptor-C98(3.36) residue and partially induce hH₄ receptor-mediated G protein activation and β-arrestin2 recruitment. Moreover, these observations confirm our previously proposed binding mode of 2-aminopyrimidines. VUF14480 will be a useful tool to stabilize the receptor into an active confirmation and further investigate the structure of the active hH₄ receptor.

Discovery of highly potent and selective D4 ligands by interactive SAR study

A series of thienylmethylphenylpiperazins was synthesized and tested for affinity towards the five subtypes of dopaminergic receptors. Compound 5f showed more than 1000 folds selectivity to D4 receptors; analogue 5e showed the highest affinity to D4 receptors with Ki 3.9 nM. An interactive SAR approach was adopted and lead to compound 14a with Ki (D4) as low as 0.03 nM. Molecular docking studies showed a potential, first to report arene cation interaction between the D4 unique residue Arg-186 and the ligands' arene moiety, explaining the importance of having a strong negative electrostatic potential at this area of the compound structure.

Chemical genomics approach for GPCR-ligand interaction prediction and extraction of ligand binding determinants

Chemical genomics research has revealed that G-protein coupled receptors (GPCRs) interact with a variety of ligands and that a large number of ligands are known to bind GPCRs even with low transmembrane (TM) sequence similarity. It is crucial to extract informative binding region propensities from large quantities of bioactivity data. To address this issue, we propose a machine learning approach that enables identification of both chemical substructures and amino acid properties that are associated with ligand binding, which can be applied to virtual ligand screening on a GPCR-wide scale. We also address the question of how to select plausible negative noninteraction pairs based on a statistical approach in order to develop reliable prediction models for GPCR-ligand interactions. The key interaction sites estimated by our approach can be of great use not only for screening of active compounds but also for modification of active compounds with the aim of improving activity or selectivity.

International Union of Basic and Clinical Pharmacology. [corrected]. LXXXVII. Complement peptide C5a, C4a, and C3a receptors

The activation of the complement cascade, a cornerstone of the innate immune response, produces a number of small (74-77 amino acid) fragments, originally termed anaphylatoxins, that are potent chemoattractants and secretagogues that act on a wide variety of cell types. These fragments, C5a, C4a, and C3a, participate at all levels of the immune response and are also involved in other processes such as neural development and organ regeneration. Their primary function, however, is in inflammation, so they are important targets for the development of antiinflammatory therapies. Only three receptors for complement peptides have been found, but there are no satisfactory antagonists as yet, despite intensive investigation. In humans, there is a single receptor for C3a (C3a receptor), no known receptor for C4a, and two receptors for C5a (C5a₁ receptor and C5a₂ receptor). The most recently characterized receptor, the C5a₂ receptor (previously known as C5L2 or GPR77), has been regarded as a passive binding protein, but signaling activities are now ascribed to it, so we propose that it be formally identified as a receptor and be given a name to reflect this. Here, we describe the complex biology of the complement peptides, introduce a new suggested nomenclature, and review our current knowledge of receptor pharmacology.

Membrane environment imposes unique selection pressures on transmembrane domains of G protein-coupled receptors

We have investigated the influence of the plasma membrane environment on the molecular evolution of G protein-coupled receptors (GPCRs), the largest receptor family in Metazoa. In particular, we have analyzed the site-specific rate variation across the two primary structural partitions, transmembrane (TM) and extramembrane (EM), of these membrane proteins. We find that TM domains evolve more slowly than do EM domains, though TM domains display increased rate heterogeneity relative to their EM counterparts. Although the majority of residues across GPCRs experience strong to weak purifying selection, many GPCRs experience positive selection at both TM and EM residues, albeit with a slight bias towards the EM. Further, a subset of GPCRs, chemosensory receptors (including olfactory and taste receptors), exhibit increased rates of evolution relative to other GPCRs, an effect which is more pronounced in their TM spans. Although it has been previously suggested that the TM's low evolutionary rate is caused by their high percentage of buried residues, we show that their attenuated rate seems to stem from the strong biophysical constraints of the membrane itself, or by functional requirements. In spite of the strong evolutionary constraints acting on the TM spans of GPCRs, positive selection and high levels of evolutionary rate variability are common. Thus, biophysical constraints should not be presumed to preclude a protein's ability to evolve.

Target based virtual screening by docking into automatically generated GPCR models

Target based virtual screening (VS) combined with high-throughput measurements is an extremely useful tool to identify small molecule hits for proteins and in particular for G-protein coupled receptors (GPCRs). However, this is a quite difficult process for GPCRs due to the paucity of 3D structural information on these receptors. Therefore, the only possibility for target based VS is to build a structural model of the GPCR to be used for docking. However, GPCR model building is a very time consuming process, if the model should be able to explain all experimental findings and this investment is not always justified, if the model is only used for VS. Thus, a fully automated workflow is presented here, where a large number of GPCR models is built, and the best model is identified to be used for docking. The workflow leads to moderate enrichments with a very low effort. The inputs required are the sequence of the targeted GPCR, a reference ligand with experimental information and a database of small molecules to be used for docking. Manual intervention is recommended at various points, but it is strictly speaking not necessary.

Identification of novel amino acid derived CCK-2R antagonists as potential antiulcer agent: homology modeling, design, synthesis, and pharmacology

The present study revisited the three-dimensional (3D) homology model of CCK-2R using human A(2a) adenosine receptor and the resolved NMR based structure of the third extracellular loop of the CCK-2R as templates. Further in order to identify novel antiulcer agents, rational designing have been performed utilizing the substructure of a well-known CCK-2R antagonist benzotript as a lead molecule and submitted to the combined docking and simulation studies. This led to the understanding of the essential structure requirement as well as variation of binding mode among conformational isomers of small molecule CCK-2R antagonists. In the next step, preparation of each configurational isomer of these molecules was carried out and submitted for their in vitro activity followed by in vivo screening into antiulcer rat model. The biological screening of these compounds has not only validated the developed homology model of CCK-2R but also led to the identification of highly potent CCK-2R antagonist 6a as an orally active and safe candidate molecule having better antiulcer properties than the well-known drug benzotript.

Mapping histamine H4 receptor–ligand binding modes

Computational prediction of ligand binding modes in G protein-coupled receptors (GPCRs) remains a challenging task. Systematic consideration of different protein modelling templates, ligand binding poses, and ligand protonation states in extensive molecular dynamics (MD) simulation studies enabled the prediction of ligand-specific mutation effects in the histamine H₄ receptor, a key player in inflammation.

Structural evolution of G-protein-coupled receptors: a sequence space approach

Class A G-protein-coupled receptors (GPCRs) provide a fascinating example of evolutionary success. In this review, we discuss how metric multidimensional scaling (MDS), a multivariate analysis method, complements traditional tree-based phylogenetic methods and helps decipher the mechanisms that drove the evolution of class A GPCRs. MDS provides low-dimensional representations of a distance matrix. Applied to a multiple sequence alignment, MDS represents the sequences in a Euclidean space as points whose interdistances are as close as possible to the distances in the alignment (the so-called sequence space). We detail how to perform the MDS analysis of a multiple sequence alignment and how to analyze and interpret the resulting sequence space. We also show that the projection of supplementary data (a property of the MDS method) can be used to straightforwardly monitor the evolutionary drift of specific subfamilies. The sequence space of class A GPCRs reveals the key role of mutations at the level of the TM2 and TM5 proline residues in the evolution of class A GPCRs.

Small and colorful stones make beautiful mosaics: fragment-based chemogenomics

Smaller stones with a wide variety of colors make a higher resolution mosaic. In much the same way, smaller chemical entities that are structurally diverse are better able to interrogate protein binding sites. This feature article describes the construction of a diverse fragment library and an analysis of the screening of six representative protein targets belonging to three diverse target classes (G protein-coupled receptors ADRB2, H1R, H3R, and H4R, the ligand-gated ion channel 5-HT3R, and the kinase PKA) using chemogenomics approaches. The integration of experimentally determined bioaffinity profiles across related and unrelated protein targets and chemogenomics analysis of fragment binding and protein structure allow the identification of: (i) unexpected similarities and differences in ligand binding properties, and (ii) subtle ligand affinity and selectivity cliffs. With a wealth of fragment screening data being generated in industry and academia, such approaches will contribute to a more detailed structural understanding of ligand-protein interactions.

Genomic perspectives in inter-individual adverse responses following nanomedicine administration: The way forward

The underlying mechanism of intravenous infusion-related adverse reactions inherent to regulatory-approved nanomedicines still remains elusive. There are substantial inter-individual differences in observed adverse reactions, which may include cardiovascular, broncho-pulmonary, muco-cutaneous, neuro-psychosomatic and autonomic manifestations. Although nanomedicine-mediated triggering of complement activation has been suggested to be a significant contributing factor to these adverse events, complement activation may still proceed in non-responders. Whether these reactions share similar immunological mechanisms and underpinning genetic factors with drug hypersensitivity syndrome remains to be investigated. Genetic association studies could be a powerful tool to dissect causative factors and reveal the multiple molecular pathways that induce infusion related adverse reactions. It is envisaged that such research may lead to the design of reliable in vitro profiling tests for risk assessment and treatment decisions, thereby revolutionizing the practice of medicine with nanopharmaceuticals. Such procedures may further improve regulatory approval processes for nanomedicines currently in the pipeline and decrease the overall cost of health care. Here we discuss some key innate immunity genes and their polymorphisms in relation to nanomedicine infusion-mediated symptomatic responses.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Evolutionary Analysis of Functional Divergence among Chemokine Receptors, Decoy Receptors, and Viral Receptors

Chemokine receptors (CKRs) function in the inflammatory response and in vertebrate homeostasis. Decoy and viral receptors are two types of CKR homologs with modified functions from those of the typical CKRs. The decoy receptors are able to bind ligands without signaling. On the other hand, the viral receptors show constitutive signaling without ligands. We examined the sites related to the functional difference. At first, the decoy and viral receptors were each classified into five groups, based on the molecular phylogenetic analysis. A multiple amino acid sequence alignment between each group and the CKRs was then constructed. The difference in the amino acid composition between the group and the CKRs was evaluated as the Kullback-Leibler (KL) information value at each alignment site. The KL information value is considered to reflect the difference in the functional constraints at the site. The sites with the top 5% of KL information values were selected and mapped on the structure of a CKR. The comparisons with decoy receptor groups revealed that the detected sites were biased on the intracellular side. In contrast, the sites detected from the comparisons with viral receptor groups were found on both the extracellular and intracellular sides. More sites were found in the ligand binding pocket in the analyses of the viral receptor groups, as compared to the decoy receptor groups. Some of the detected sites were located in the GPCR motifs. For example, the DRY motif of the decoy receptors was often degraded, although the motif of the viral receptors was basically conserved. The observations for the viral receptor groups suggested that the constraints in the pocket region are loose and that the sites on the intracellular side are different from those for the decoy receptors, which may be related to the constitutive signaling activity of the viral receptors.

Pharmacological modulation of chemokine receptor function

G protein-coupled chemokine receptors and their peptidergic ligands are interesting therapeutic targets due to their involvement in various immune-related diseases, including rheumatoid arthritis, multiple sclerosis, inflammatory bowel disease, chronic obstructive pulmonary disease, HIV-1 infection and cancer. To tackle these diseases, a lot of effort has been focused on discovery and development of small-molecule chemokine receptor antagonists. This has been rewarded by the market approval of two novel chemokine receptor inhibitors, AMD3100 (CXCR4) and Maraviroc (CCR5) for stem cell mobilization and treatment of HIV-1 infection respectively. The recent GPCR crystal structures together with mutagenesis and pharmacological studies have aided in understanding how small-molecule ligands interact with chemokine receptors. Many of these ligands display behaviour deviating from simple competition and do not interact with the chemokine binding site, providing evidence for an allosteric mode of action. This review aims to give an overview of the evidence supporting modulation of this intriguing receptor family by a range of ligands, including small molecules, peptides and antibodies. Moreover, the computer-assisted modelling of chemokine receptor-ligand interactions is discussed in view of GPCR crystal structures. Finally, the implications of concepts such as functional selectivity and chemokine receptor dimerization are considered.

Farnesyl pyrophosphate is an endogenous antagonist to ADP-stimulated P2Y₁₂ receptor-mediated platelet aggregation

Farnesyl pyrophosphate (FPP) is an intermediate in cholesterol biosynthesis, and it has also been reported to activate platelet LPA (lysophosphatidic acid) receptors. The aim of this study was to investigate the role of extracellular FPP in platelet aggregation. Human platelets were studied with light transmission aggregometry, flow cytometry and [³⁵S]GTPγS binding assays. As shown previously, FPP could potentiate LPA-stimulated shape change. Surprisingly, FPP also acted as a selective insurmountable antagonist to ADP-induced platelet aggregation. FPP inhibited ADP-induced expression of P-selectin and the activated glycoprotein (Gp)IIb/IIIa receptor. FPP blocked ADP-induced inhibition of cAMP accumulation and [³⁵S]GTPγS binding in platelets. In Chinese hamster ovary cells expressing the P2Y₁₂ receptor, FPP caused a rightward shift of the [³⁵S]GTPγS binding curve. In Sf9 insect cells expressing the human P2Y₁₂ receptor, FPP showed a concentration-dependent, although incomplete inhibition of [³H]PSB-0413 binding. Docking of FPP in a P2Y₁₂ receptor model revealed molecular similarities with ADP and a good fit into the binding pocket for ADP. In conclusion, FPP is an insurmountable antagonist of ADP-induced platelet aggregation mediated by the P2Y₁₂ receptor. It could be an endogenous antithrombotic factor modulating the strong platelet aggregatory effects of ADP in a manner similar to the use of clopidogrel, prasugrel or ticagrelor in the treatment of ischaemic heart disease.

Life beyond kinases: structure-based discovery of sorafenib as nanomolar antagonist of 5-HT receptors

Of great interest in recent years has been computationally predicting the novel polypharmacology of drug molecules. Here, we applied an "induced-fit" protocol to improve the homology models of 5-HT(2A) receptor, and we assessed the quality of these models in retrospective virtual screening. Subsequently, we computationally screened the FDA approved drug molecules against the best induced-fit 5-HT(2A) models and chose six top scoring hits for experimental assays. Surprisingly, one well-known kinase inhibitor, sorafenib, has shown unexpected promiscuous 5-HTRs binding affinities, K(i) = 1959, 56, and 417 nM against 5-HT(2A), 5-HT(2B), and 5-HT(2C), respectively. Our preliminary SAR exploration supports the predicted binding mode and further suggests sorafenib to be a novel lead compound for 5HTR ligand discovery. Although it has been well-known that sorafenib produces anticancer effects through targeting multiple kinases, carefully designed experimental studies are desirable to fully understand whether its "off-target" 5-HTR binding activities contribute to its therapeutic efficacy or otherwise undesirable side effects.

A prospective cross-screening study on G-protein-coupled receptors: lessons learned in virtual compound library design

We present the systematic prospective evaluation of a protein-based and a ligand-based virtual screening platform against a set of three G-protein-coupled receptors (GPCRs): the β-2 adrenoreceptor (ADRB2), the adenosine A(2A) receptor (AA2AR), and the sphingosine 1-phosphate receptor (S1PR1). Novel bioactive compounds were identified using a consensus scoring procedure combining ligand-based (frequent substructure ranking) and structure-based (Snooker) tools, and all 900 selected compounds were screened against all three receptors. A striking number of ligands showed affinity/activity for GPCRs other than the intended target, which could be partly attributed to the fuzziness and overlap of protein-based pharmacophore models. Surprisingly, the phosphodiesterase 5 (PDE5) inhibitor sildenafil was found to possess submicromolar affinity for AA2AR. Overall, this is one of the first published prospective chemogenomics studies that demonstrate the identification of novel cross-pharmacology between unrelated protein targets. The lessons learned from this study can be used to guide future virtual ligand design efforts.

Determination of different putative allosteric binding pockets at the lutropin receptor by using diverse drug-like low molecular weight ligands

The lutropin/choriogonadotrophin receptor (LHCGR) is a family A G protein-coupled receptor (GPCR) which binds the endogenous hormone-ligands at the large extracellular domain. In contrast, several drug-like low-molecular-weight ligands (LMWs) have been reported to interact allosterically within the seven transmembrane domain (7TMD) of the LHCGR. Here, we were interested to study the putative allosteric LHCGR binding region with focus on the determination of two pockets for LMW ligands. A library of compounds was screened for their ability to modify the binding of an allosteric radiolabeled LMW agonist [³H]Org 43553. Further experimental and computational studies revealed that the putative binding pocket for a newly identified allosteric enhancer (LUF5419) and a previously described allosteric inhibitor (LUF5771) are overlapping and that this site is different from the Org 43553 binding site. The present study showed that these compounds are useful tools to reveal details on different allosteric binding sites located within the 7TMD of the LHCGR.

Molecular evolution of the neuropeptide S receptor

The neuropeptide S receptor (NPSR) is a recently deorphanized member of the G protein-coupled receptor (GPCR) superfamily and is activated by the neuropeptide S (NPS). NPSR and NPS are widely expressed in central nervous system and are known to have crucial roles in asthma pathogenesis, locomotor activity, wakefulness, anxiety and food intake. The NPS-NPSR system was previously thought to have first evolved in the tetrapods. Here we examine the origin and the molecular evolution of the NPSR using in-silico comparative analyses and document the molecular basis of divergence of the NPSR from its closest vertebrate paralogs. In this study, NPSR-like sequences have been identified in a hemichordate and a cephalochordate, suggesting an earlier emergence of a NPSR-like sequence in the metazoan lineage. Phylogenetic analyses revealed that the NPSR is most closely related to the invertebrate cardioacceleratory peptide receptor (CCAPR) and the group of vasopressin-like receptors. Gene structure features were congruent with the phylogenetic clustering and supported the orthology of NPSR to the invertebrate NPSR-like and CCAPR. A site-specific analysis between the vertebrate NPSR and the well studied paralogous vasopressin-like receptor subtypes revealed several putative amino acid sites that may account for the observed functional divergence between them. The data can facilitate experimental studies aiming at deciphering the common features as well as those related to ligand binding and signal transduction processes specific to the NPSR.