Molecular dynamics simulation of the P2Y14 receptor. Ligand docking and identification of a putative binding site of the distal hexose moiety

Kinetic Mechanism of Surfactant-Based Molecular Recognition: Selective Permeability across an Oil-Water Interface Regulated by Supramolecular Aggregates

Lipids are known to play a vital role in the molecular organization of all cellular life. Molecular recognition is another fundamental biological process that is generally attributed to biological polymers, such as proteins and nucleic acids. However, there is evidence that aggregates of lipids and lipid-like molecules are also capable of selectively binding to or regulating the partitioning of other molecules. We previously demonstrated that a model two-phase octanol/water system can selectively partition Red 40 and Blue 1 dyes added to an aqueous phase, with the selectivity depending on the surfactant (e.g., cetyltrimethylammonium bromide) dissolved in the organic phase. Here, we elucidate the mechanism of molecular recognition in this system by using quantitative partitioning experiments and molecular dynamics (MD) simulations. Our results indicate that the selectivity for the red dye is thermodynamically favored at all surfactant concentrations, while selectivity for the blue dye is kinetically favored at high surfactant concentrations. The kinetic selectivity for the blue dye correlates with the presence of molecular aggregation at the oil-water interface. Coarse-grained MD simulations elucidate nanoscale supramolecular structures that can preferentially bind one small molecule rather than another at an interface, providing a selectively permeable barrier in the absence of proteins. The results suggest a new supramolecular mechanism for molecular recognition with potential applications in drug delivery, drug discovery, and biosensing.

Role of the Purinergic P2Y2 Receptor in Pulmonary Hypertension

Pulmonary arterial hypertension (PAH), group 1 pulmonary hypertension (PH), is a fatal disease that is characterized by vasoconstriction, increased pressure in the pulmonary arteries, and right heart failure. PAH can be described by abnormal vascular remodeling, hyperproliferation in the vasculature, endothelial cell dysfunction, and vascular tone dysregulation. The disease pathomechanisms, however, are as yet not fully understood at the molecular level. Purinergic receptors P2Y within the G-protein-coupled receptor family play a major role in fluid shear stress transduction, proliferation, migration, and vascular tone regulation in systemic circulation, but less is known about their contribution in PAH. Hence, studies that focus on purinergic signaling are of great importance for the identification of new therapeutic targets in PAH. Interestingly, the role of P2Y2 receptors has not yet been sufficiently studied in PAH, whereas the relevance of other P2Ys as drug targets for PAH was shown using specific agonists or antagonists. In this review, we will shed light on P2Y receptors and focus more on the P2Y2 receptor as a potential novel player in PAH and as a new therapeutic target for disease management.

Tools and drugs for uracil nucleotide-activated P2Y receptors

P2Y receptors (P2YRs) are a family of G protein-coupled receptors activated by extracellular nucleotides. Physiological P2YR agonists include purine and pyrimidine nucleoside di- and triphosphates, such as ATP, ADP, UTP, UDP, nucleotide sugars, and dinucleotides. Eight subtypes exist, P2Y₁, P2Y₂, P2Y₄, P2Y₆, P2Y₁₁, P2Y₁₂, P2Y₁₃, and P2Y₁₄, which represent current or potential future drug targets. Here we provide a comprehensive overview of ligands for the subgroup of the P2YR family that is activated by uracil nucleotides: P2Y₂ (UTP, also ATP and dinucleotides), P2Y₄ (UTP), P2Y₆ (UDP), and P2Y₁₄ (UDP, UDP-glucose, UDP-galactose). The physiological agonists are metabolically unstable due to their fast hydrolysis by ectonucleotidases. A number of agonists with increased potency, subtype-selectivity and/or enzymatic stability have been developed in recent years. Useful P2Y₂R agonists include MRS2698 (6-01, highly selective) and PSB-1114 (6-05, increased metabolic stability). A potent and selective P2Y₂R antagonist is AR-C118925 (10-01). For studies of the P2Y₄R, MRS4062 (3-15) may be used as a selective agonist, while PSB-16133 (10-06) is a selective antagonist. Several potent P2Y₆R agonists have been developed including 5-methoxyuridine 5'-O-((R_p)α-boranodiphosphate) (6-12), PSB-0474 (3-11), and MRS2693 (3-26). The isocyanate MRS2578 (10-08) is used as a selective P2Y₆R antagonist, although its reactivity and low water-solubility are limiting. With MRS2905 (6-08), a potent and metabolically stable P2Y₁₄R agonist is available, while PPTN (10-14) represents a potent and selective P2Y₁₄R antagonist. The radioligand [³H]UDP can be used to label P2Y₁₄Rs. In addition, several fluorescent probes have been developed. Uracil nucleotide-activated P2YRs show great potential as drug targets, especially in inflammation, cancer, cardiovascular and neurodegenerative diseases.

Molecular modeling of the human P2Y14 receptor: A template for structure-based design of selective agonist ligands

The P2Y14 receptor (P2Y14R) is a Gi protein-coupled receptor that is activated by uracil nucleotides UDP and UDP-glucose. The P2Y14R structure has yet to be solved through X-ray crystallography, but the recent agonist-bound crystal structure of the P2Y12R provides a potentially suitable template for its homology modeling for rational structure-based design of selective and high-affinity ligands. In this study, we applied ligand docking and molecular dynamics refinement to a P2Y14R homology model to qualitatively explain structure-activity relationships of previously published synthetic nucleotide analogues and to probe the quality of P2Y14R homology modeling as a template for structure-based design. The P2Y14R model supports the hypothesis of a conserved binding mode of nucleotides in the three P2Y12-like receptors involving functionally conserved residues. We predict phosphate group interactions with R253(6.55), K277(7.35), Y256(6.58) and Q260(6.62), nucleobase (anti-conformation) π-π stacking with Y102(3.33) and the role of F191(5.42) as a means for selectivity among P2Y12-like receptors. The glucose moiety of UDP-glucose docked in a secondary subpocket at the P2Y14R homology model. Thus, P2Y14R homology modeling may allow detailed prediction of interactions to facilitate the design of high affinity, selective agonists as pharmacological tools to study the P2Y14R.

Structure-based drug design using GPCR homology modeling: toward the discovery of novel selective CysLT2 antagonists

3D structure of CysLT2 receptor was constructed by using homology modeling and molecular simulations. The binding pocket of CysLT2 receptor and the proposition of the interaction mode between CysLT2 and HAMI3379 were identified. A series of dicarboxylated chalcones was then virtually evaluated through molecular docking studies. A total of six compounds 13a-f with preferable scores was further synthesized and tested for CysLT2 antagonistic activities by determination of the cytosolic free Ca(2+) levels in HEK293 cells. Compounds 13e and 13f exhibited potent and selective CysLT2 antagonistic activities with IC50 values being 7.5 and 0.25 μM, respectively.

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.

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.

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.

Forced unbinding of GPR17 ligands from wild type and R255I mutant receptor models through a computational approach

Background

GPR17 is a hybrid G-protein-coupled receptor (GPCR) activated by two unrelated ligand families, extracellular nucleotides and cysteinyl-leukotrienes (cysteinyl-LTs), and involved in brain damage and repair. Its exploitment as a target for novel neuro-reparative strategies depends on the elucidation of the molecular determinants driving binding of purinergic and leukotrienic ligands. Here, we applied docking and molecular dynamics simulations (MD) to analyse the binding and the forced unbinding of two GPR17 ligands (the endogenous purinergic agonist UDP and the leukotriene receptor antagonist pranlukast from both the wild-type (WT) receptor and a mutant model, where a basic residue hypothesized to be crucial for nucleotide binding had been mutated (R255I) to Ile.

Results

MD suggested that GPR17 nucleotide binding pocket is enclosed between the helical bundle and extracellular loop (EL) 2. The driving interaction involves R255 and the UDP phosphate moiety. To support this hypothesis, steered MD experiments showed that the energy required to unbind UDP is higher for the WT receptor than for R255I. Three potential binding sites for pranlukast where instead found and analysed. In one of its preferential docking conformations, pranlukast tetrazole group is close to R255 and phenyl rings are placed into a subpocket highly conserved among GPCRs. Pulling forces developed to break polar and aromatic interactions of pranlukast were comparable. No differences between the WT receptor and the R255I receptor were found for the unbinding of pranlukast.

Conclusions

These data thus suggest that, in contrast to which has been hypothesized for nucleotides, the lack of the R255 residue doesn't affect the binding of pranlukast a crucial role for R255 in binding of nucleotides to GPR17. Aromatic interactions are instead likely to play a predominant role in the recognition of pranlukast, suggesting that two different binding subsites are present on GPR17.

Human P2Y(14) receptor agonists: truncation of the hexose moiety of uridine-5'-diphosphoglucose and its replacement with alkyl and aryl groups

Uridine-5'-diphosphoglucose (UDPG) activates the P2Y(14) receptor, a neuroimmune system GPCR. P2Y(14) receptor tolerates glucose substitution with small alkyl or aryl groups or its truncation to uridine 5'-diphosphate (UDP), a full agonist at the human P2Y(14) receptor expressed in HEK-293 cells. 2-Thiouracil derivatives displayed selectivity for activation of the human P2Y(14) vs the P2Y(6) receptor, such as 2-thio-UDP 4 (EC(50) = 1.92 nM at P2Y(14), 224-fold selectivity vs P2Y(6)) and its beta-propyloxy ester 18. EC(50) values of the beta-methyl ester of UDP and its 2-thio analogue were 2730 and 56 nM, respectively. beta-tert-Butyl ester of 4 was 11-fold more potent than UDPG, but beta-aryloxy or larger, branched beta-alkyl esters, such as cyclohexyl, were less potent. Ribose replacement of UDP with a rigid North or South methanocarba (bicyclo[3.1.0]hexane) group abolished P2Y(14) receptor agonist activity. alpha,beta-Methylene and difluoromethylene groups were well tolerated at the P2Y(14) receptor and are expected to provide enhanced stability in biological systems. alpha,beta-Methylene-2-thio-UDP 11 (EC(50) = 0.92 nM) was 2160-fold selective versus P2Y(6). Thus, these nucleotides and their congeners may serve as important pharmacological probes for the detection and characterization of the P2Y(14) receptor.

Enhanced potency of nucleotide-dendrimer conjugates as agonists of the P2Y14 receptor: multivalent effect in G protein-coupled receptor recognition

The P2Y(14) receptor is a G protein-coupled receptor activated by uridine-5'-diphosphoglucose and other nucleotide sugars that modulates immune function. Covalent conjugation of P2Y(14) receptor agonists to PAMAM (polyamidoamine) dendrimers enhanced pharmacological activity. Uridine-5'-diphosphoglucuronic acid (UDPGA) and its ethylenediamine adduct were suitable functionalized congeners for coupling to several generations (G2.5-6) of dendrimers (both terminal carboxy and amino). Prosthetic groups, including biotin for avidin complexation, a chelating group for metal complexation (and eventual magnetic resonance imaging), and a fluorescent moiety, also were attached with the eventual goals of molecular detection and characterization of the P2Y(14) receptor. The activities of conjugates were assayed in HEK293 cells stably expressing the human P2Y(14) receptor. A G3 PAMAM conjugate containing 20 bound nucleotide moieties (UDPGA) was 100-fold more potent (EC(50) 2.4 nM) than the native agonist uridine-5'-diphosphoglucose. A molecular model of this conjugate docked in the human P2Y(14) receptor showed that the nucleotide-substituted branches could extend far beyond the dimensions of the receptor and be available for multivalent docking to receptor aggregates. Larger dendrimer carriers and greater loading favored higher potency. A similar conjugate of G6 with 147 out of 256 amino groups substituted with UDPGA displayed an EC(50) value of 0.8 nM. Thus, biological activity was either retained or dramatically enhanced in the multivalent dendrimer conjugates in comparison with monomeric P2Y(14) receptor agonists, depending on size, degree of substitution, terminal functionality, and attached prosthetic groups.

Molecular recognition in the P2Y(14) receptor: Probing the structurally permissive terminal sugar moiety of uridine-5'-diphosphoglucose

The P2Y(14) receptor, a nucleotide signaling protein, is activated by uridine-5'-diphosphoglucose 1 and other uracil nucleotides. We have determined that the glucose moiety of 1 is the most structurally permissive region for designing analogues of this P2Y(14) agonist. For example, the carboxylate group of uridine-5'-diphosphoglucuronic acid proved to be suitable for flexible substitution by chain extension through an amide linkage. Functionalized congeners containing terminal 2-acylaminoethylamides prepared by this strategy retained P2Y(14) activity, and molecular modeling predicted close proximity of this chain to the second extracellular loop of the receptor. In addition, replacement of glucose with other sugars did not diminish P2Y(14) potency. For example, the [5'']ribose derivative had an EC(50) of 0.24muM. Selective monofluorination of the glucose moiety indicated a role for the 2''- and 6''-hydroxyl groups of 1 in receptor recognition. The beta-glucoside was twofold less potent than the native alpha-isomer, but methylene replacement of the 1''-oxygen abolished activity. Replacement of the ribose ring system with cyclopentyl or rigid bicyclo[3.1.0]hexane groups abolished activity. Uridine-5'-diphosphoglucose also activates the P2Y(2) receptor, but the 2-thio analogue and several of the potent modified-glucose analogues were P2Y(14)-selective.

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.

Development of selective high affinity antagonists, agonists, and radioligands for the P2Y1 receptor

The P2Y(1) receptor is a member of the P2Y family of nucleotide-activated G protein-coupled receptors, and it is an important therapeutic target based on its broad tissue distribution and essential role in platelet aggregation. We have designed a set of highly selective and diverse pharmacological tools for studying the P2Y(1) receptor using a rational approach to ligand design. Based on the discovery that bisphosphate analogues of the P2Y(1) receptor agonist, ADP, are partial agonists/competitive antagonists of this receptor, an iterative approach was used to develop competitive antagonists with enhanced affinity and selectivity. Halogen substitutions of the 2-position of the adenine ring provided increased affinity while an N(6) methyl substitution eliminated partial agonist activity. Furthermore, various replacements of the ribose ring with symmetrically branched, phosphorylated acyclic structures revealed that the ribose is not necessary for recognition at the P2Y(1) receptor. Finally, replacement of the ribose ring with a five member methanocarba ring constrained in the Northern conformation conferred dramatic increases in affinity to both P2Y(1) receptor antagonists as well as agonists. These combined structural modifications have resulted in a series of selective high affinity antagonists of the P2Y(1) receptor, two broadly applicable radioligands, and a high affinity agonist capable of selectively activating the P2Y(1) receptor in human platelets. Complementary receptor modeling and computational ligand docking have provided a putative structural framework for the drug-receptor interactions. A similar rational approach is being applied to develop selective ligands for other subtypes of P2Y receptors.

Development of selective agonists and antagonists of P2Y receptors

Although elucidation of the medicinal chemistry of agonists and antagonists of the P2Y receptors has lagged behind that of many other members of group A G protein-coupled receptors, detailed qualitative and quantitative structure–activity relationships (SARs) were recently constructed for several of the subtypes. Agonists selective for P2Y₁, P2Y₂, and P2Y₆ receptors and nucleotide antagonists selective for P2Y₁ and P2Y₁₂ receptors are now known. Selective nonnucleotide antagonists were reported for P2Y₁, P2Y₂, P2Y₆, P2Y₁₁, P2Y₁₂, and P2Y₁₃ receptors. At the P2Y₁ and P2Y₁₂ receptors, nucleotide agonists (5′-diphosphate derivatives) were converted into antagonists of nanomolar affinity by altering the phosphate moieties, with a focus particularly on the ribose conformation and substitution pattern. Nucleotide analogues with conformationally constrained ribose-like rings were introduced as selective receptor probes for P2Y₁ and P2Y₆ receptors. Screening chemically diverse compound libraries has begun to yield new lead compounds for the development of P2Y receptor antagonists, such as competitive P2Y₁₂ receptor antagonists with antithrombotic activity. Selective agonists for the P2Y₄, P2Y₁₁, and P2Y₁₃ receptors and selective antagonists for P2Y₄ and P2Y₁₄ receptors have not yet been identified. The P2Y₁₄ receptor appears to be the most restrictive of the class with respect to modification of the nucleobase, ribose, and phosphate moieties. The continuing process of ligand design for the P2Y receptors will aid in the identification of new clinical targets.

Mapping the UDP-glucuronic acid binding site in UDP-glucuronosyltransferase-1A10 by homology-based modeling: confirmation with biochemical evidence

The UDP-glucuronosyltransferase (UGT) isozyme system is critical for protecting the body against endogenous and exogenous chemicals by linking glucuronic acid donated by UDP-glucuronic acid to a lipophilic acceptor substrate. UGTs convert metabolites, dietary constituents, and environmental toxicants to highly excretable glucuronides. Because of difficulties associated with purifying endoplasmic reticulum-bound UGTs for structural studies, we carried out homology-based computer modeling to aid analysis. The search found structural homology in Escherichia coli UDP-galactose 4-epimerase. Consistent with predicted similarities involving the common UDP moiety in substrate/inhibitor, UDP-glucose and UDP-hexanol amine caused competitive inhibition by Lineweaver-Burk plots. Among predicted binding sites N292, K314, K315, and K404 in UGT1A10, two informative sets of mutants K314R/Q/A/E/G and K404R/E had null activities or 2.7-fold higher/50% less activity, respectively. Scatchard analysis of binding data of the affinity ligand, 5-azidouridine-[beta- (32)P]diphosphoglucuronic acid, to purified UGT1A10-His or UGT1A7-His revealed high- and low-affinity binding sites. 2-Nitro-5-thiocyanobenzoic acid-digested UGT1A10-His bound with the radiolabeled affinity ligand revealed an 11.3 and 14.3 kDa peptide associated with K314 and K404, respectively, in a discontinuous SDS-PAGE system. Similar treatment of 1A10His-K314A bound with the ligand lacked both peptides; 1A10-HisK404R- and 1A10-HisK404E showed 1.3-fold greater and 50% less label in the 14.3 kDa peptide, respectively, compared to 1A10-His without affecting the 11.3 kDa peptide. Scatchard analysis of binding data of the affinity ligand to 1A10His-K404R and -K404E showed a 6-fold reduction and a large increase in K d, respectively. Our results indicate that K314 and K404 are required UDP-glcA binding sites in 1A10, that K404 controls activity and high-affinity sites, and that K314 and K404 are strictly conserved in 70 aligned UGTs, except for S321, equivalent to K314, in UGT2B15 and 2B17 and I321 in the inactive UGT8, which suggests UGT2B15 and 2B17 contain suboptimal activity. Hence our data strongly support UDP-glcA binding to K314 and K404 in UGT1A10.

GPR17: molecular modeling and dynamics studies of the 3-D structure and purinergic ligand binding features in comparison with P2Y receptors

Background

GPR17 is a G-protein-coupled receptor located at intermediate phylogenetic position between two distinct receptor families: the P2Y and CysLT receptors for extracellular nucleotides and cysteinyl-LTs, respectively. We previously showed that GPR17 can indeed respond to both classes of endogenous ligands and to synthetic compounds active at the above receptor families, thus representing the first fully characterized non-peptide "hybrid" GPCR. In a rat brain focal ischemia model, the selective in vivo knock down of GPR17 by anti-sense technology or P2Y/CysLT antagonists reduced progression of ischemic damage, thus highlighting GPR17 as a novel therapeutic target for stroke. Elucidation of the structure of GPR17 and of ligand binding mechanisms are the necessary steps to obtain selective and potent drugs for this new potential target. On this basis, a 3-D molecular model of GPR17 embedded in a solvated phospholipid bilayer and refined by molecular dynamics simulations has been the first aim of this study. To explore the binding mode of the "purinergic" component of the receptor, the endogenous agonist UDP and two P2Y receptor antagonists demonstrated to be active on GPR17 (MRS2179 and cangrelor) were then modeled on the receptor.

Results

Molecular dynamics simulations suggest that GPR17 nucleotide binding pocket is similar to that described for the other P2Y receptors, although only one of the three basic residues that have been typically involved in ligand recognition is conserved (Arg255). The binding pocket is enclosed between the helical bundle and covered at the top by EL2. Driving interactions are H-bonds and salt bridges between the 6.55 and 6.52 residues and the phosphate moieties of the ligands. An "accessory" binding site in a region formed by the EL2, EL3 and the Nt was also found.

Conclusion

Nucleotide binding to GPR17 occurs on the same receptor regions identified for already known P2Y receptors. Agonist/antagonist binding mode are similar, but not identical. An accessory external binding site could guide small ligands to the deeper principal binding site in a multi-step mechanism of activation. The nucleotide binding pocket appears to be unable to allocate the leukotrienic type ligands in the same effective way.

Structure-activity relationship of uridine 5'-diphosphoglucose analogues as agonists of the human P2Y14 receptor

UDP-glucose (UDPG) and derivatives are naturally occurring agonists of the Gi protein-coupled P2Y14 receptor, which occurs in the immune system. We synthesized and characterized pharmacologically novel analogues of UDPG modified on the nucleobase, ribose, and glucose moieties, as the basis for designing novel ligands in conjunction with modeling. The recombinant human P2Y14 receptor expressed in COS-7 cells was coupled to phospholipase C through an engineered Galpha-q/i protein. Most modifications of the uracil or ribose moieties abolished activity; this is among the least permissive P2Y receptors. However, a 2-thiouracil modification in 15 (EC50 49 +/- 2 nM) enhanced the potency of UDPG (but not UDP-glucuronic acid) by 7-fold. 4-Thio analogue 13 was equipotent to UDPG, but S-alkylation was detrimental. Compound 15 was docked in a rhodposin-based receptor homology model, which correctly predicted potent agonism of UDP-fructose, UDP-mannose, and UDP-inositol. The hexose moiety of UDPG interacts with multiple H-bonding and charged residues and provides a fertile region for agonist modification.