Structure and ligand-binding site characteristics of the human P2Y11 nucleotide receptor deduced from computational modelling and mutational analysis

Primary cilia, A-kinase anchoring proteins and constitutive activity at the orphan G protein-coupled receptor GPR161: A tale about a tail

Primary cilia are non-motile antennae-like structures responsible for sensing environmental changes in most mammalian cells. Ciliary signalling is largely mediated by the Sonic Hedgehog (Shh) pathway, which acts as a master regulator of ciliary protein transit and is essential for normal embryonic development. One particularly important player in primary cilia is the orphan G protein-coupled receptor, GPR161. In this review, we introduce GPR161 in the context of Shh signalling and describe the unique features on its C-terminus such as PKA phosphorylation sites and an A-kinase anchoring protein motif, which may influence the function of the receptor, cAMP compartmentalisation and/or trafficking within primary cilia. We discuss the recent putative pairing of GPR161 and spexin-1, highlighting the additional steps needed before GPR161 could be considered 'deorphanised'. Finally, we speculate that the marked constitutive activity and unconventional regulation of GPR161 may indicate that the receptor may not require an endogenous ligand.

Molecular pharmacology of P2Y receptor subtypes

Professor Geoffrey Burnstock proposed the concept of purinergic signaling via P1 and P2 receptors. P2Y receptors are G-protein-coupled receptors (GPCRs) for extracellular adenine and uracil nucleotides. Eight mammalian P2Y receptor subtypes have been identified. They are divided into two subgroups (P2Y₁, P2Y₂, P2Y₄, P2Y₆, and P2Y₁₁) and (P2Y₁₂, P2Y₁₃, and P2Y₁₄). P2Y receptors are found in almost all cells and mediate responses in physiology and pathophysiology including pain and inflammation. The antagonism of platelet P2Y₁₂ receptors by cangrelor, ticagrelor or active metabolites of the thienopyridine compounds ticlopidine, clopidogrel and prasugrel reduces the ADP-induced platelet aggregation in patients with thrombotic complications of vascular diseases. The nucleotide agonist diquafosol acting at P2Y₂ receptors is used for the treatment of the dry eye syndrome. Structural information obtained by crystallography of the human P2Y₁ and P2Y₁₂ receptor proteins, site-directed mutagenesis and molecular modeling will facilitate the rational design of novel selective drugs.

Molecular Modeling Applied to the Discovery of New Lead Compounds for P2 Receptors Based on Natural Sources

P2 receptors are a family of transmembrane receptors activated by nucleotides and nucleosides. Two classes have been described in mammals, P2X and P2Y, which are implicated in various diseases. Currently, only P2Y12 has medicines approved for clinical use as antiplatelet agents and natural products have emerged as a source of new drugs with action on P2 receptors due to the diversity of chemical structures. In drug discovery, in silico virtual screening (VS) techniques have become popular because they have numerous advantages, which include the evaluation of thousands of molecules against a target, usually proteins, faster and cheaper than classical high throughput screening (HTS). The number of studies using VS techniques has been growing in recent years and has led to the discovery of new molecules of natural origin with action on different P2X and P2Y receptors. Using different algorithms it is possible to obtain information on absorption, distribution, metabolism, toxicity, as well as predictions on biological activity and the lead-likeness of the selected hits. Selected biomolecules may then be tested by molecular dynamics and, if necessary, rationally designed or modified to improve their interaction for the target. The algorithms of these in silico tools are being improved to permit the precision development of new drugs and, in the future, this process will take the front of drug development against some central nervous system (CNS) disorders. Therefore, this review discusses the methodologies of in silico tools concerning P2 receptors, as well as future perspectives and discoveries, such as the employment of artificial intelligence in drug discovery.

P2Y Purinergic Receptors, Endothelial Dysfunction, and Cardiovascular Diseases

Purinergic G-protein-coupled receptors are ancient and the most abundant group of G-protein-coupled receptors (GPCRs). The wide distribution of purinergic receptors in the cardiovascular system, together with the expression of multiple receptor subtypes in endothelial cells (ECs) and other vascular cells demonstrates the physiological importance of the purinergic signaling system in the regulation of the cardiovascular system. This review discusses the contribution of purinergic P2Y receptors to endothelial dysfunction (ED) in numerous cardiovascular diseases (CVDs). Endothelial dysfunction can be defined as a shift from a “calm” or non-activated state, characterized by low permeability, anti-thrombotic, and anti-inflammatory properties, to a “activated” state, characterized by vasoconstriction and increased permeability, pro-thrombotic, and pro-inflammatory properties. This state of ED is observed in many diseases, including atherosclerosis, diabetes, hypertension, metabolic syndrome, sepsis, and pulmonary hypertension. Herein, we review the recent advances in P2Y receptor physiology and emphasize some of their unique signaling features in pulmonary endothelial cells.

P2Y11 Receptors: Properties, Distribution and Functions

The P2Y₁₁ receptor is a G protein-coupled receptor that is stimulated by endogenous purine nucleotides, particularly ATP. Amongst P2Y receptors it has several unique properties; (1) it is the only human P2Y receptor gene that contains an intron in the coding sequence; (2) the gene does not appear to be present in the rodent genome; (3) it couples to stimulation of both phospholipase C and adenylyl cyclase. Its absence in mice and rats, along with a limited range of selective pharmacological tools, has hampered the development of our knowledge and understanding of its properties and functions. Nonetheless, through a combination of careful use of the available tools, suppression of receptor expression using siRNA and genetic screening for SNPs, possible functions of native P2Y₁₁ receptors have been identified in a variety of human cells and tissues. Many are in blood cells involved in inflammatory responses, consistent with extracellular ATP being a damage-associated signalling molecule in the immune system. Thus proposed potential therapeutic applications relate, in the main, to modulation of acute and chronic inflammatory responses.

Human P2Y11 Expression Level Affects Human P2X7 Receptor-Mediated Cell Death

Adenosine triphosphate (ATP) is known to induce cell death in T lymphocytes at high extracellular concentrations. CD4⁺ and CD8⁺ T lymphocytes have a differential response to ATP, which in mice is due to differences in the P2X7 receptor expression levels. By contrast, we observed that the difference in human CD4⁺ and CD8⁺ T lymphocyte response toward the synthetic ATP-analog BzATP is not explained by a difference in human P2X7 receptor expression. Rather, the BzATP-induced human P2X7 receptor response in naïve and immune-activated lymphocyte subtypes correlated with the expression of another ATP-binding receptor: the human P2Y₁₁ receptor. In a recombinant expression system, the coexpression of the human P2Y₁₁ receptor counteracted BzATP-induced human P2X7 receptor-driven lactate dehydrogenase release (a marker of cell death) and pore formation independent of calcium signaling. A mutated non-signaling human P2Y₁₁ receptor had a similar human P2X7 receptor-inhibitory effect on pore formation, thus demonstrating that the human P2X7 receptor interference was not caused by human P2Y₁₁ receptor signaling. In conclusion, we demonstrate an important species difference in the ATP-mediated cell death between mice and human cells and show that in human T lymphocytes, the expression of the human P2Y₁₁ receptor correlates with human P2X7 receptor-driven cell death following BzATP stimulation.

Identification of Highly Promising Antioxidants/Neuroprotectants Based on Nucleoside 5'-Phosphorothioate Scaffold. Synthesis, Activity, and Mechanisms of Action

With a view to identify novel and biocompatible neuroprotectants, we designed nucleoside 5'-thiophosphate analogues, 6-11. We identified 2-SMe-ADP(α-S), 7A, as a most promising neuroprotectant. 7A reduced ROS production in PC12 cells under oxidizing conditions, IC50 of 0.08 vs 21 μM for ADP. Furthermore, 7A rescued primary neurons subjected to oxidation, EC50 of 0.04 vs 19 μM for ADP. 7A is a most potent P2Y1-R agonist, EC50 of 0.0026 μM. Activity of 7A in cells involved P2Y1/12-R as indicated by blocking P2Y12-R or P2Y1-R. Compound 7A inhibited Fenton reaction better than EDTA, IC50 of 37 vs 54 μM, due to radical scavenging, IC50 of 12.5 vs 30 μM for ADP, and Fe(II)-chelation, IC50 of 80 vs >200 μM for ADP (ferrozine assay). In addition, 7A was stable in human blood serum, t1/2 of 15 vs 1.5 h for ADP, and resisted hydrolysis by NPP1/3, 2-fold vs ADP. Hence, we propose 7A as a highly promising neuroprotectant.

Mutational mapping of the transmembrane binding site of the G-protein coupled receptor TGR5 and binding mode prediction of TGR5 agonists

TGR5 (Gpbar-1, M-Bar) is a class A G-protein coupled bile acid-sensing receptor predominately expressed in brain, liver and gastrointestinal tract, and a promising drug target for the treatment of metabolic disorders. Due to the lack of a crystal structure of TGR5, the development of TGR5 agonists has been guided by ligand-based approaches so far. Three binding mode models of bile acid derivatives have been presented recently. However, they differ from one another in terms of overall orientation or with respect to the location and interactions of the cholane scaffold, or cannot explain all results from mutagenesis experiments. Here, we present an extended binding mode model based on an iterative and integrated computational and biological approach. An alignment of 68 TGR5 agonists based on this binding mode leads to a significant and good structure-based 3D QSAR model, which constitutes the most comprehensive structure-based 3D-QSAR study of TGR5 agonists undertaken so far and suggests that the binding mode model is a close representation of the "true" binding mode. The binding mode model is further substantiated in that effects predicted for eight mutations in the binding site agree with experimental analyses on the impact of these TGR5 variants on receptor activity. In the binding mode, the hydrophobic cholane scaffold of taurolithocholate orients towards the interior of the orthosteric binding site such that rings A and B are in contact with TM5 and TM6, the taurine side chain orients towards the extracellular opening of the binding site and forms a salt bridge with R79(EL1), and the 3-hydroxyl group forms hydrogen bonds with E169(5.44) and Y240(6.51). The binding mode thus differs in important aspects from the ones recently presented. These results are highly relevant for the development of novel, more potent agonists of TGR5 and should be a valuable starting point for the development of TGR5 antagonists, which could show antiproliferative effects in tumor cells.

Nucleotides Acting at P2Y Receptors: Connecting Structure and Function

Eight G protein-coupled P2Y receptor (P2YR) subtypes are important physiologic mediators. The human P2YRs are fully activated by ATP (P2Y2 and P2Y11), ADP (P2Y1, P2Y12, and P2Y13), UTP (P2Y2 and P2Y4), UDP (P2Y6 and P2Y14), and UDP glucose (P2Y14). Their structural elucidation is progressing rapidly. The X-ray structures of three ligand complexes of the Gi-coupled P2Y12R and two of the Gq-coupled P2Y1Rs were recently determined and will be especially useful in structure-based ligand design at two P2YR subfamilies. These high-resolution structures, which display unusual binding site features, complement mutagenesis studies for probing ligand recognition and activation. The structural requirements for nucleotide agonist recognition at P2YRs are relatively permissive with respect to the length of the phosphate moiety, but less so with respect to base recognition. Nucleotide-like antagonists and partial agonists are also known for P2Y1, P2Y2, P2Y4, and P2Y12Rs. Each P2YR subtype has the ability to be activated by structurally bifunctional agonists, such as dinucleotides, typically, dinucleoside triphosphates or tetraphosphates, and nucleoside polyphosphate sugars (e.g., UDP glucose) as well as the more conventional mononucleotide agonists. A range of dinucleoside polyphosphates, from triphosphates to higher homologs, occurs naturally. Earlier modeling predictions of the P2YRs were not very accurate, but recent findings have provided much detailed structural insight into this receptor family to aid in the rational design of new drugs.

5-OMe-uridine-5'-O-(α-boranodiphosphate), a novel nucleotide derivative highly active at the human P2Y(6) receptor protects against death-receptor mediated glial apoptosis

P2Y receptors are activated by nucleotides and involved in numerous physiological/pathophysiological processes. However, investigations of specific P2Y receptor functions have been hampered by lack of suitable receptor agonists-antagonists. Recently, we identified the nucleotide 5-OMe-UDP as potent and selective agonist for human P2Y6 receptors. We studied a series of derivatives of this analog with a Pα-borano group substituting a non-bridging oxygen and found increased potency and receptor specificity. Rp-5-OMe-UDPαB (Rp-5-OMe-uridine 5'-O-α-boranodiphosphate) was most potent and selective in inducing intracellular calcium signaling in 1321N1 astrocytoma cells expressing the human P2Y6 receptor. Here, we investigated whether Rp-5-OMe-UDPαB evokes cell protection through human P2Y6 receptors. We tested a well-established model, tumor necrosis factor α (TNFα)-induced cell death in 1321N1 astrocytoma cells. Rp-5-OMe-UDPαB inhibited TNFα-induced cell death even stronger than UDP. These first data of a neuro-protective activity of the human P2Y6 receptor emphasize the potential of the stable, selective, and potent Rp-5-OMe-UDPαB analog for exploiting P2Y6 receptor-mediated cellular functions, like cytoprotection in human tissues, with suitability for future neuro-protective drug development.

Highly potent and selective ectonucleotide pyrophosphatase/phosphodiesterase I inhibitors based on an adenosine 5'-(α or γ)-thio-(α,β- or β,γ)-methylenetriphosphate scaffold

Aberrant nucleotide pyrophosphatase/phosphodiesterase-1 (NPP1) activity is associated with chondrocalcinosis, osteoarthritis, and type 2 diabetes. The potential of NPP1 inhibitors as therapeutic agents, and the scarceness of their structure-activity relationship, encouraged us to develop new NPP1 inhibitors. Specifically, we synthesized ATP-α-thio-β,γ-CH2 (1), ATP-α-thio-β,γ-CCl2 (2), ATP-α-CH2-γ-thio (3), and 8-SH-ATP (4) and established their resistance to hydrolysis by NPP1,3 and NTPDase1,2,3,8 (<5% hydrolysis) (NTPDase = ectonucleoside triphosphate diphosphohydrolase). Analogues 1-3 at 100 μM inhibited thymidine 5'-monophosphate p-nitrophenyl ester hydrolysis by NPP1 and NPP3 by >90% and 23-43%, respectively, and only slightly affected (0-40%) hydrolysis of ATP by NTPDase1,2,3,8. Analogue 3 is the most potent NPP1 inhibitor currently known, Ki = 20 nM and IC50 = 0.39 μM. Analogue 2a is a selective NPP1 inhibitor with Ki = 685 nM and IC50 = 0.57 μM. Analogues 1-3 were found mostly to be nonagonists of P2Y1/P2Y2/P2Y11 receptors. Docking analogues 1-3 into the NPP1 model suggested that activity correlates with the number of H-bonds with binding site residues. In conclusion, we propose analogues 2a and 3 as highly promising NPP1 inhibitors.

Alanine-(87)-threonine polymorphism impairs signaling and internalization of the human P2Y11 receptor, when co-expressed with the P2Y1 receptor

The P2Y11 nucleotide receptor detects high extracellular ATP concentrations. Mutations of the human P2RY11 gene can play a role in brain autoimmune responses, and the P2Y11 receptor alanine-87-threonine (A87T) polymorphism has been suggested to affect immune-system functions. We investigated receptor functionality of the P2Y11 A87T mutant using HEK293 and 1321N1 astrocytoma cells. In HEK293 cells, the P2Y11 receptor agonist 3'-O-(4-benzoylbenzoyl)adenosine 5'-triphosphate (BzATP) was completely inactive in evoking intracellular calcium release while the potency of ATP was reduced. ATP was also less potent in triggering cAMP generation. However, 1321N1 astrocytoma cells, which lack any endogenous P2Y1 receptors, did not display a reduction. Only when 1321N1 cells were co-transfected with P2Y11 A87T and P2Y1 receptors, the calcium responses to the P2Y11 receptor-specific agonist BzATP were reduced. It is already known that P2Y1 and P2Y11 receptors interact. We thus conclude that the physiological impact of A87T mutation of the P2Y11 receptor derives from detrimental effects on P2Y1 -P2Y11 receptor interaction. We additionally investigated alanine-87-serine and alanine-87-tyrosine P2Y11 receptor mutants. Both mutations rescue the response to BzATP in HEK293 cells, thus ruling out polarity of amino acid-87 to be the molecular basis for altered receptor characteristics. We further found that the P2Y11 A87T receptor shows complete loss of nucleotide-induced internalization in HEK293 cells. Thus, we demonstrate impaired signaling of the P2Y11 A87T-mutated receptors when co-operating with P2Y1 receptors.

Comparative Modeling of Drug Target Proteins☆

In this perspective, we begin by describing the comparative protein structure modeling technique and the accuracy of the corresponding models. We then discuss the significant role that comparative prediction plays in drug discovery. We focus on virtual ligand screening against comparative models and illustrate the state-of-the-art by a number of specific examples.

Sp-2-propylthio-ATP-α-B and Sp-2-propylthio-ATP-α-B,β-γ-dichloromethylene are novel potent and specific agonists of the human P2Y₁₁ receptor

The human P2Y₁₁ nucleotide receptor mRNA was found in virtually all human tissues, and the receptor serves many physiological roles, such as immune response regulation. The Ala-87-Thr-P2Y₁₁ receptor single nucleotide polymorphism was linked to increased risk for acute myocardial infarction. To facilitate the development of new therapeutic applications involving cells expressing several P2 receptor subtypes, the availability of specific and potent agonists is mandatory. Here, we synthesized a series of novel adenine nucleotide derivatives, based upon the potent P2Y₁₁ receptor agonists AR-C67085. Features of the novel nucleotide derivatives are a propylthio substitution at C2-adenine and a Pα-borano or Pα-thio substitution of non-bridging oxygen atom. The latter substitutions introduce a chiral center at the α-phosphate. Sp-isomers of Pα-borano- and Rp-isomers of Pα-thio-substituted nucleotides are preferred by the P2Y₁₁ receptor. As recently reported by us, diastereoselectivity of the P2Y₁₁ receptor is opposite to that of the P2Y₁ receptor. Therefore, we exploit this characteristic to increase nucleotide selectivity. At the P2Y₁₁ receptor, the Sp-isomers of 2-propylthio-ATP-α-B (2B) and 2-propylthio-ATP-α-B,β-γ-dichloromethylene (4B) were the most potent of the novel nucleotide series, with EC₅₀ values of 0.03 μM for both, being ca. 80-fold more potent than 2-propylthio-ATP and ATP (EC₅₀ = 2.6 μM). We conclude that the borano-substitution at the α-phosphate of 2-propylthio-ATP enhances nucleotide potency at the P2Y₁₁ receptor. The combination with a Pβ-Pγ-dichloromethylene group in 4B results in a nucleotide, which shows higher selectivity for the P2Y₁₁ receptor over the P2Y₁₁ receptor than 2B making it the most promising of the novel P2Y₁₁ receptor agonists.

Highly efficient biocompatible neuroprotectants with dual activity as antioxidants and P2Y receptor agonists

Currently, there is a need for novel, biocompatible, and effective neuroprotectants for the treatment of neurodegenerative diseases and brain injury associated with oxidative damage. Here, we developed nucleotide-based neuroprotectants acting dually as antioxidants and P2Y-R agonists. To improve the potency, selectivity, and metabolic stability of ATP/ADP, we substituted adenine C2-position by Cl and Pα/Pβ position by borano group, 6-9. Nucleotides 6-9 inhibited oxidation in cell-free systems (Fe(II)-H2O2), as detected by ESR (IC50 up to 175 μM), and ABTS assay (IC50 up to 40 μM). They also inhibited FeSO4-induced oxidative stress in PC12 cells (IC50 of 80-200 nM). 2-Cl-ADP(α-BH3), 7a, was found to be the most potent P2Y1-R agonist currently known (EC50 7 nM) and protected primary cortical neurons from FeSO4 insult (EC50 170 nM). In addition, it proved to be metabolically stable in human blood serum (t(1/2) 7 vs 1.5 h for ADP). Hence, we propose 7a as a highly promising neuroprotectant.

Ligand-specific binding and activation of the human adenosine A(2B) receptor

Adenosine A(2B) receptors, which play a role in inflammation and cancer, are of considerable interest as novel drug targets. To gain deeper insights into ligand binding and receptor activation, we exchanged amino acids predicted to be close to the binding pocket. The alanine mutants were stably expressed in CHO cells and characterized by radioligand binding and cAMP assays using three structural classes of ligands: xanthine (antagonist), adenosine, and aminopyridine derivatives (agonists). Asn282(7.45) and His280(7.43) were found to stabilize the binding site by intramolecular hydrogen bond formation as in the related A(2A) receptor subtype. Trp247(6.48), Val250(6.51), and particularly Ser279(7.42) were shown to be important for binding of nucleosidic agonists. Leu81(3.28), Asn186(5.42), and Val250(6.51) were discovered to be crucial for binding of the xanthine-derived antagonist PSB-603. Leu81(3.28), which is not conserved among adenosine receptor subtypes, may be important for the high selectivity of PSB-603. The N186(5.42)A mutant resulted in an increased potency for agonists. The interactions of the non-nucleosidic agonist BAY60-6583 were different from those of the nucleosides: while BAY60-6583 appeared not to interact with Ser279(7.42), its interactions with Trp247(6.48) and Val250(6.51) were significantly weaker compared to those of NECA. Moreover, our results discount the hypothesis of Trp247(6.48) serving as a "toogle switch" because BAY60-6583 was able to activate the corresponding mutant. This study reveals distinct interactions of structurally diverse ligands with the human A(2B) receptor and differences between closely related receptor subtypes (A(2B) and A(2A)). It will contribute to the understanding of G protein-coupled receptor function and advance A(2B) receptor ligand design.

Probing GPCR structure: adenosine and P2Y nucleotide receptors

The adenosine receptors (ARs) provide an example of how to accurately predict ligand recognition, even prior to the availability of a crystallographic structure. Homology modeling has been used to gain structural insight, in conjunction with site-directed mutagenesis, and structure-activity relationships of small molecular ligands. Recent X-ray structures greatly improved the accuracy of knowledge of AR ligand recognition and furthermore characterized conformational changes induced by receptor activation. Now, homology modeling extends these structural insights to related GPCRs and suggests new ligand structures. This strategy is also being applied to the eight subtypes of P2Y receptors for extracellular nucleotides, which lack X-ray structures and are best modeled by homology to the CXCR4 (peptide) receptor. Neoceptors, as studied for three of the four AR subtypes, create a molecular complementarity between a mutant receptor and a chemically tailored agonist ligand to selectively enhance affinity, implying direct physical contact and thus validating docking hypotheses.

Molecular Structure of P2Y Receptors: Mutagenesis, Modeling, and Chemical Probes

There are eight subtypes of P2Y receptors (P2YRs) that are activated, and in some cases inhibited, by a range of extracellular nucleotides. These nucleotides are ubiquitous, but their extracellular concentration can rise dramatically in response to hypoxia, ischemia, or mechanical stress, injury, and release through channels and from vesicles. Two subclasses of P2YRs were defined based on clustering of sequences, second messengers, and receptor sequence analysis. The numbering system for P2YR subtypes is discontinuous; i.e., P2Y(1-14)Rs have been defined, but six of the intermediate-numbered cloned receptor sequences (e.g., P2y(3), P2y(5), P2y(7-10)) are not functional mammalian nucleotide receptors. Of these two clusters, the P2Y(12-14) subtypes couple via Gα(i) to inhibit adenylate cyclase, while the remaining subtypes couple through Gα(q) to activate phospholipase C. Collectively, the P2YRs respond to both purine and pyrimidine nucleotides, in the form of 5'-mono- and dinucleotides and nucleoside-5'-diphosphosugars. In recent years, the medicinal chemistry of P2Y receptors has advanced significantly, to provide selective agonists and antagonists for many but not all of the subtypes. Ligand design has been aided by insights from structural probing using molecular modelling and mutagenesis. Currently, the molecular modelling of the receptors is effectively based on the X-ray structure of the CXCR4 receptor, which is the closest to the P2Y receptors among all the currently crystallized receptors in terms of sequence similarity. It is now a challenge to develop novel and selective P2YR ligands for disease treatment (although antagonists of the P2Y(12)R are already widely used as antithrombotics).

G protein-coupled adenosine (P1) and P2Y receptors: ligand design and receptor interactions

The medicinal chemistry and pharmacology of the four subtypes of adenosine receptors (ARs) and the eight subtypes of P2Y receptors (P2YRs, activated by a range of purine and pyrimidine mono- and dinucleotides) has recently advanced significantly leading to selective ligands. X-ray crystallographic structures of both agonist- and antagonist-bound forms of the A(2A)AR have provided unprecedented three-dimensional detail concerning molecular recognition in the binding site and the conformational changes in receptor activation. It is apparent that this ubiquitous cell signaling system has implications for understanding and treating many diseases. ATP and other nucleotides are readily released from intracellular sources under conditions of injury and organ stress, such as hypoxia, ischemia, or mechanical stress, and through channels and vesicular release. Adenosine may be generated extracellularly or by cellular release. Therefore, depending on pathophysiological factors, in a given tissue, there is often a tonic activation of one or more of the ARs or P2YRs that can be modulated by exogenous agents for a beneficial effect. Thus, this field has provided fertile ground for pharmaceutical development, leading to clinical trials of selective receptor ligands as imaging agents or for conditions including cardiac arrhythmias, ischemia/reperfusion injury, diabetes, pain, thrombosis, Parkinson's disease, rheumatoid arthritis, psoriasis, dry eye disease, pulmonary diseases such as cystic fibrosis, glaucoma, cancer, chronic hepatitis C, and other diseases.

UDP made a highly promising stable, potent, and selective P2Y6-receptor agonist upon introduction of a boranophosphate moiety

P2Y(6) nucleotide receptor (P2Y(6)-R) plays important physiological roles, such as insulin secretion and reduction of intraocular pressure. However, this receptor is still lacking potent and selective agonists to be used as potential drugs. Here, we synthesized uracil nucleotides and dinucleotides, substituted at the C5 and/or P(α) position with methoxy and/or borano groups, 18-22. Compound 18A, R(p) isomer of 5-OMe-UDP(α-B), is the most potent and P2Y(6)-R selective agonist currently known (EC(50) 0.008 μM) being 19-fold more potent than UDP and showing no activity at uridine nucleotide receptors, P2Y(2)- and P2Y(4)-R. Analogue 18A was highly chemically stable under conditions mimicking gastric juice acidity (t(1/2) = 16.9 h). It was more stable to hydrolysis by nucleotide pyrophosphatases (NPP1,3) than UDP (15% and 28% hydrolysis by NPP1 and NPP3, respectively, vs 50% and 51% hydrolysis of UDP) and metabolically stable in blood serum (t(1/2) = 17 vs 2.4, 11.9, and 21 h for UDP, 5-OMe-UDP, and UDP(α-B), respectively). This newly discovered highly potent and physiologically stable P2Y(6)-R agonist may be of future therapeutic potential.

Virtual ligand screening against comparative protein structure models

Virtual ligand screening uses computation to discover new ligands of a protein by screening one or more of its structural models against a database of potential ligands. Comparative protein structure modeling extends the applicability of virtual screening beyond the atomic structures determined by X-ray crystallography or NMR spectroscopy. Here, we describe an integrated modeling and docking protocol, combining comparative modeling by MODELLER and virtual ligand screening by DOCK.

Molecular pharmacology, physiology, and structure of the P2Y receptors

The P2Y receptors are a widely expressed group of eight nucleotide-activated G protein-coupled receptors (GPCRs). The P2Y(1)(ADP), P2Y(2)(ATP/UTP), P2Y(4)(UTP), P2Y(6)(UDP), and P2Y(11)(ATP) receptors activate G(q) and therefore robustly promote inositol lipid signaling responses. The P2Y(12)(ADP), P2Y(13)(ADP), and P2Y(14)(UDP/UDP-glucose) receptors activate G(i) leading to inhibition of adenylyl cyclase and to Gβγ-mediated activation of a range of effector proteins including phosphoinositide 3-kinase-γ, inward rectifying K(+) (GIRK) channels, phospholipase C-β2 and -β3, and G protein-receptor kinases 2 and 3. A broad range of physiological responses occur downstream of activation of these receptors ranging from Cl(-) secretion by epithelia to aggregation of platelets to neurotransmission. Useful structural models of the P2Y receptors have evolved from extensive genetic analyses coupled with molecular modeling based on three-dimensional structures obtained for rhodopsin and several other GPCRs. Selective ligands have been synthesized for most of the P2Y receptors with the most prominent successes attained with highly selective agonist and antagonist molecules for the ADP-activated P2Y(1) and P2Y(12) receptors. The widely prescribed drug, clopidogrel, which results in irreversible blockade of the platelet P2Y(12) receptor, is the most important therapeutic agent that targets a P2Y receptor.

Mining human genome for novel purinergic P2Y receptors: a sequence analysis and molecular modeling approach

The purinergic P2Y receptors are G-protein coupled receptors (GPCRs) that control many physiological processes by mediating cellular responses to purines, pyrimidines and their analogues. They can be used as potential therapeutic targets in a variety of disease conditions. Therefore, it is critical to identify new members of this family of receptors from the human genome and characterize them for their role in health and disease. In the present work, molecular modeling was carried out for the 21 known P2Y receptors. Binding site analysis was done on the basis of docking and site-directed mutagenesis data. Thus, conserved features of P2Y receptors could be formulated. These features can be used to determine the purinergic nature of potential P2Y receptors in the human genome. We applied this knowledge to human genome GPCR sequences found by sensitive sequence search techniques and identified two orphan receptors, namely GPR34 and GP171 that have all the necessary conserved features of P2Y receptors.

Deconstructing nucleotide binding activity of the Mdm2 RING domain

Mdm2, a central negative regulator of the p53 tumor suppressor, possesses a Really Interesting New Gene (RING) domain within its C-terminus. In addition to E3 ubiquitin ligase activity, the Mdm2 RING preferentially binds adenine base nucleotides, and such binding leads to a conformational change in the Mdm2 C-terminus. Here, we present further biochemical analysis of the nucleotide–Mdm2 interaction. We have found that MdmX, an Mdm2 family member with high sequence homology, binds adenine nucleotides with similar affinity and specificity as Mdm2, suggesting that residues involved in nucleotide binding may be conserved between the two proteins and adenosine triphosphate (ATP) binding may have similar functional consequences for both Mdm family members. By generating and testing a series of proteins with deletions and substitution mutations within the Mdm2 RING, we mapped the specific adenine nucleotide binding region of Mdm2 to residues 429–484, encompassing the minimal RING domain. Using a series of ATP derivatives, we demonstrate that phosphate coordination by the Mdm2 P-loop contributes to, but is not primarily responsible for, ATP binding. Additionally, we have identified the 2′ and 3′ hydroxyls of the ribose and the C6 amino group of the adenine base moiety as being essential for binding.

P2Y nucleotide receptors: promise of therapeutic applications

Extracellular nucleotides, such as ATP and UTP, have distinct signaling roles through a class of G-protein-coupled receptors, termed P2Y. The receptor ligands are typically charged molecules of low bioavailability and stability in vivo. Recent progress in the development of selective agonists and antagonists for P2Y receptors and study of knockout mice have led to new drug concepts based on these receptors. The rapidly accelerating progress in this field has already resulted in drug candidates for cystic fibrosis, dry eye disease and thrombosis. On the horizon are novel treatments for cardiovascular diseases, inflammatory diseases and neurodegeneration.

Molecular docking screens using comparative models of proteins

Two orders of magnitude more protein sequences can be modeled by comparative modeling than have been determined by X-ray crystallography and NMR spectroscopy. Investigators have nevertheless been cautious about using comparative models for ligand discovery because of concerns about model errors. We suggest how to exploit comparative models for molecular screens, based on docking against a wide range of crystallographic structures and comparative models with known ligands. To account for the variation in the ligand-binding pocket as it binds different ligands, we calculate "consensus" enrichment by ranking each library compound by its best docking score against all available comparative models and/or modeling templates. For the majority of the targets, the consensus enrichment for multiple models was better than or comparable to that of the holo and apo X-ray structures. Even for single models, the models are significantly more enriching than the template structure if the template is paralogous and shares more than 25% sequence identity with the target.

NF546 [4,4'-(carbonylbis(imino-3,1-phenylene-carbonylimino-3,1-(4-methyl-phenylene)-carbonylimino))-bis(1,3-xylene-alpha,alpha'-diphosphonic acid) tetrasodium salt] is a non-nucleotide P2Y11 agonist and stimulates release of interleukin-8 from human monocyte-derived dendritic cells

The G protein-coupled P2Y(11) receptor is involved in immune system modulation. In-depth physiological evaluation is hampered, however, by a lack of selective and potent ligands. By screening a library of sulfonic and phosphonic acid derivatives at P2Y(11) receptors recombinantly expressed in human 1321N1 astrocytoma cells (calcium and cAMP assays), the selective non-nucleotide P2Y(11) agonist NF546 [4,4'-(carbonylbis(imino-3,1-phenylene-carbonylimino-3,1-(4-methyl-phenylene)carbonylimino))-bis(1,3-xylene-alpha,alpha'-diphosphonic acid) tetrasodium salt] was identified. NF546 had a pEC(50) of 6.27 and is relatively selective for P2Y(11) over P2Y(1), P2Y(2), P2Y(4), P2Y(6), P2Y(12), P2X(1), P2X(2), and P2X(2)-X(3). Adenosine-5'-O-(3-thio)triphosphate (ATPgammaS), a nonhydrolyzable analog of the physiological P2Y(11) agonist ATP, and NF546 use a common binding site as suggested by molecular modeling studies and their competitive behavior toward the nanomolar potency antagonist NF340 [4,4'-(carbonylbis(imino-3,1-(4-methyl-phenylene)carbonylimino))bis(naphthalene-2,6-disulfonic acid) tetrasodium salt] in Schild analysis. The pA(2) of NF340 was 8.02 against ATPgammaS and 8.04 against NF546 (calcium assays). NF546 was further tested for P2Y(11)-mediated effects in monocyte-derived dendritic cells. Similarly to ATPgammaS, NF546 led to thrombospondin-1 secretion and inhibition of lipopolysaccharide-stimulated interleukin-12 release, whereas NF340 inhibited these effects. Further, for the first time, it was shown that ATPgammaS or NF546 stimulation promotes interleukin 8 (IL-8) release from dendritic cells, which could be inhibited by NF340. In conclusion, we have described the first selective, non-nucleotide agonist NF546 for P2Y(11) receptors in both recombinant and physiological expression systems and could show a P2Y(11)-stimulated IL-8 release, further supporting the immunomodulatory role of P2Y(11) receptors.

Interaction of new, very potent non-nucleotide antagonists with Arg256 of the human platelet P2Y12 receptor

The P2Y(12) receptor plays a crucial role in platelet aggregation. In the present study, we analyzed the properties of non-nucleotide antagonists at the recombinant human P2Y(12) receptor and searched for amino acids involved in the molecular interaction. Receptor function was assessed by measuring the cAMP response element (CRE)-directed luciferase expression in Chinese hamster ovary cells. The cellular cAMP production was accelerated by forskolin; 2-methylthio-ADP was used to activate the wild-type P2Y(12) receptor or mutant constructs. 2-Methylthio-ADP inhibited the CRE-dependent luciferase expression with an IC(50) value of approximately 1 nM. The anthraquinone derivative reactive blue 2 used at increasing concentrations shifted the concentration-response curve of 2-methylthio-ADP to the right in a manner compatible with competitive antagonism (pA(2) value, 7.4). Its analog, 1-amino-4-[4-phenylamino-3-sulfophenylamino]-9,10-dioxo-9,10-dihydroanthracene-2-sulfonate (PSB-0739), showed a markedly higher antagonistic potency with a pA(2) value of 9.8. In cells expressing the R256A-mutant receptor, the potencies of both reactive blue 2 (apparent pK(B), 5.9) and PSB-0739 (apparent pK(B), 9.1) were decreased. The same was true for the pure reactive blue 2 meta- and para-isomers and for the ortho-isomer cibacron blue 3GA. In contrast, the analog, 1-amino-4-[4-anilino-phenylamino]-9,10-dioxo-9,10-dihydroanthracene-2-sulfonate, lacking a sulfonic acid residue at ring D (PSB-0826), showed similar pK(B) values at wild-type (8.4) and R256A-mutant receptors (8.3). In summary, the results demonstrate that PSB-0739 is the most potent competitive non-nucleotide antagonist at the human P2Y(12) receptor described so far. The results also indicate that the sulfonic acid residue at ring D is involved in the interaction of antagonists derived from reactive blue 2 with the residue Arg256 of the human P2Y(12) receptor.

Key determinants of nucleotide-activated G protein-coupled P2Y(2) receptor function revealed by chemical and pharmacological experiments, mutagenesis and homology modeling

The P2Y(2) receptor, which is activated by UTP, ATP, and dinucleotides, was studied as a prototypical nucleotide-activated GPCR. A combination of receptor mutagenesis, determination of its effects on potency and efficacy of agonists and antagonists, homology modeling, and chemical experiments was applied. R272 (extracellular loop EL3) was found to play a gatekeeper role, presumably responsible for recognition and orientation of the nucleotides. R272 is also directly involved in binding of dinucleotides, which behaved as partial agonists. Y118A (3.37) mutation led to dramatically reduced efficacy of agonists; it is part of the entry channel as well as the triphosphate binding site. While the Y114A (3.33) mutation did not have any effect on agonist activities, the antagonist Reactive Blue 2 (6) was completely inactive at that mutant. The disulfide bridge Cys25-Cys278 was found to be important for agonist potency but neither for agonist efficacy nor for antagonist potency.

Toward the three-dimensional structure and lysophosphatidic acid binding characteristics of the LPA(4)/p2y(9)/GPR23 receptor: a homology modeling study

Lysophosphatidic acid (LPA) is a naturally occurring phospholipid that initiates a broad array of biological processes, including those involved in cell proliferation, survival and migration via activation of specific G protein-coupled receptors located on the cell surface. To date, at least five receptor subtypes (LPA(1-5)) have been identified. The LPA(1-3) receptors are members of the endothelial cell differentiation gene (Edg) family. LPA(4), a member of the purinergic receptor family, and the recently identified LPA(5) are structurally distant from the canonical Edg LPA(1-3) receptors. LPA(4) and LPA(5) are linked to G(q), G(12/13) and G(s) but not G(i), while LPA(1-3) all couple to G(i) in addition to G(q) and G(12/13). There is also evidence that LPA(4) and LPA(5) are functionally different from the Edg LPA receptors. Computational modeling has provided useful information on the structure-activity relationship (SAR) of the Edg LPA receptors. In this work, we focus on the initial analysis of the structural and ligand-binding properties of LPA(4), a prototype non-Edg LPA receptor. Three homology models of the LPA(4) receptor were developed based on the X-ray crystal structures of the ground state and photoactivated bovine rhodopsin and the recently determined human beta(2)-adrenergic receptor. Docking studies of LPA in the homology models were then conducted, and plausible LPA binding loci were explored. Based on these analyses, LPA is predicted to bind to LPA(4) in an orientation similar to that reported for LPA(1-3), but through a different network of hydrogen bonds. In LPA(1-3), the ligand polar head group is reported to interact with residues at positions 3.28, 3.29 and 7.36, whereas three non-conserved amino acid residues, S114(3.28), T187(EL2) and Y265(6.51), are predicted to interact with the polar head group in the LPA(4) receptor models.

Diastereoselectivity of the P2Y11 nucleotide receptor: mutational analysis

BACKGROUND AND PURPOSE

The P2Y(11) receptor, a member of the group of metabotropic nucleotide receptors, shows a stereospecific ligand recognition of P(alpha)-substituted ATP derivatives (ATP-alpha-S isomers). These compounds are suitable candidates for the development of selective P2Y(11) receptor agonists that might be used as immune modulators. We have analysed the binding mode of ATP at the P2Y(11) receptor by molecular modeling and site-directed mutagenesis. Based on our recent findings, we decided to decipher the molecular determinants of stereoselective recognition at the P2Y(11) receptor.

EXPERIMENTAL APPROACH

Two amino acid residues [Glu186 in the extracellular loop 2 and Arg268 in the transmembrane domain 6 (TM6)], which are part of the nucleotide-binding pocket, were selected and studied by mutational analyses. We expected these residues to be involved in determining the stereospecificity of the P2Y(11) receptor.

KEY RESULTS

After mutation of Arg268 to alanine or glutamine, the stereospecific recognition of the ATP-alpha-S isomers at the P2Y(11) receptor was lost. In contrast, at the Glu186Ala receptor mutant, the stereoselective differentiation between these isomers was increased. On the Arg268Gln/Glu186Ala double mutant we observed no further effect, except for additivity in the decrease in potency of both isomers, as compared with the single-point mutants.

CONCLUSIONS AND IMPLICATIONS

Our results show that the stereospecificity of the P2Y(11) receptor for P(alpha)-substituted ATP derivatives is largely determined by the basic residue Arg268 in TM6. This will allow the design of receptor-subtype selective ligands.

Involvement of basic amino acid residues in transmembrane regions 6 and 7 in agonist and antagonist recognition of the human platelet P2Y(12)-receptor

The P2Y(12)-receptor plays a prominent role in ADP-induced platelet aggregation. In the present study, we searched for amino acid residues involved in ligand recognition of the human P2Y(12)-receptor. Wild-type or mutated receptors were expressed in 1321N1 astrocytoma cells and Chinese hamster ovary (CHO) cells. There were no major differences in cellular expression of the constructs. Cellular cAMP production and cAMP response element (CRE)-dependent luciferase expression was increased by isoproterenol (astrocytoma cells) or forskolin (CHO cells). In cells expressing wild-type receptors, R256K or S101A mutant constructs, 2-methylthio-ADP inhibited the induced cAMP production with IC(50) concentrations of about 0.3nM. In cells expressing R256A constructs, the IC(50) concentration amounted to 25nM. In cells expressing H253A/R256A, Y259D and K280A constructs, 2-methylthio-ADP failed to affect the cellular cAMP production. Moreover, in cells expressing Y259D and K280A constructs, 2-methylthio-ADP did also not change the forskolin-induced CRE-dependent luciferase expression and caused only small increases in the serum response element-dependent luciferase expression. The antagonist cangrelor had similar potencies at wild-type receptors and R256A constructs (apparent pK(B)-value at wild-type receptors: 9.2). In contrast, reactive blue-2 had a lower potency at the R256A construct (apparent pK(B)-value at wild-type receptors: 7.6). In summary, the data indicate the involvement of Arg256, Tyr259 and, possibly, H253 (transmembrane region TM6) as well as Lys280 (TM7) in the function of the human P2Y(12)-receptor. Arg256 appears to play a role in the recognition of nucleotide agonists and the non-nucleotide antagonist reactive blue-2, but no role in the recognition of the nucleotide antagonist cangrelor.

Hetero-oligomerization of the P2Y11 receptor with the P2Y1 receptor controls the internalization and ligand selectivity of the P2Y11 receptor

Nucleotides signal through purinergic receptors such as the P2 receptors, which are subdivided into the ionotropic P2X receptors and the metabotropic P2Y receptors. The diversity of functions within the purinergic receptor family is required for the tissue-specificity of nucleotide signalling. In the present study, hetero-oligomerization between two metabotropic P2Y receptor subtypes is established. These receptors, P2Y1 and P2Y11, were found to associate together when co-expressed in HEK293 cells. This association was detected by co-pull-down, immunoprecipitation and FRET (fluorescence resonance energy transfer) experiments. We found a striking functional consequence of the interaction between the P2Y11 receptor and the P2Y1 receptor where this interaction promotes agonist-induced internalization of the P2Y11 receptor. This is remarkable because the P2Y11 receptor by itself is not able to undergo endocytosis. Co-internalization of these receptors was also seen in 1321N1 astrocytoma cells co-expressing both P2Y11 and P2Y1 receptors, upon stimulation with ATP or the P2Y1 receptor-specific agonist 2-MeS-ADP. 1321N1 astrocytoma cells do not express endogenous P2Y receptors. Moreover, in HEK293 cells, the P2Y11 receptor was found to functionally associate with endogenous P2Y1 receptors. Treatment of HEK293 cells with siRNA (small interfering RNA) directed against the P2Y1 receptor diminished the agonist-induced endocytosis of the heterologously expressed GFP-P2Y11 receptor. Pharmacological characteristics of the P2Y11 receptor expressed in HEK293 cells were determined by recording Ca2+ responses after nucleotide stimulation. This analysis revealed a ligand specificity which was different from the agonist profile established in cells expressing the P2Y11 receptor as the only metabotropic nucleotide receptor. Thus the hetero-oligomerization of the P2Y1 and P2Y11 receptors allows novel functions of the P2Y11 receptor in response to extracellular nucleotides.

A novel nucleotide receptor in Xenopus activates the cAMP second messenger pathway

We describe a Xenopus P2Y receptor that shares only weak homology with members of the mammalian P2Y family, being most similar to human P2Y(11). When activated by nucleotide analogs, it stimulates both calcium and cAMP mobilization pathways, a feature unique, among mammalian P2Y receptors, to P2Y(11). Activity can be blocked by compounds known to act as antagonists of mammalian P2Y(11). Genomic synteny between Xenopus and mammals suggests that the novel gene is a true ortholog of P2Y(11). Xenopus P2Y(11) is transcribed during embryonic development, beginning at gastrulation, and is enriched in the developing nervous system.