Virtual screening leads to the discovery of novel non-nucleotide P2Y₁ receptor antagonists

Exploring bias in platelet P2Y1 signalling: Host defence versus haemostasis

Platelets are necessary for maintaining haemostasis. Separately, platelets are important for the propagation of inflammation during the host immune response against infection. The activation of platelets also causes inappropriate inflammation in various disease pathologies, often in the absence of changes to haemostasis. The separate functions of platelets during inflammation compared with haemostasis are therefore varied and this will be reflected in distinct pathways of activation. The activation of platelets by the nucleotide adenosine diphosphate (ADP) acting on P2Y1 and P2Y12 receptors is important for the development of platelet thrombi during haemostasis. However, P2Y1 stimulation of platelets is also important during the inflammatory response and paradoxically in scenarios where no changes to haemostasis and platelet aggregation occur. In these events, Rho-GTPase signalling, rather than the canonical phospholipase Cβ (PLCβ) signalling pathway, is necessary. We describe our current understanding of these differences, reflecting on recent advances in knowledge of P2Y1 structure, and the possibility of biased agonism occurring from activation via other endogenous nucleotides compared with ADP. Knowledge arising from these different pathways of P2Y1 stimulation of platelets during inflammation compared with haemostasis may help therapeutic control of platelet function during inflammation or infection, while preserving essential haemostasis. LINKED ARTICLES: This article is part of a themed issue on Platelet purinergic receptor and non-thrombotic disease. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.4/issuetoc.

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.

Update of P2Y receptor pharmacology: IUPHAR Review 27

Eight G protein-coupled P2Y receptor subtypes respond to extracellular adenine and uracil mononucleotides and dinucleotides. P2Y receptors belong to the δ group of rhodopsin-like GPCRs and contain two structurally distinct subfamilies: P2Y₁ , P2Y₂ , P2Y₄ , P2Y₆ , and P2Y₁₁ (principally G_q protein-coupled P2Y₁ -like) and P2Y_12-14 (principally G_i protein-coupled P2Y₁₂ -like) receptors. Brain P2Y receptors occur in neurons, glial cells, and vasculature. Endothelial P2Y₁ , P2Y₂ , P2Y₄ , and P2Y₆ receptors induce vasodilation, while smooth muscle P2Y₂ , P2Y₄ , and P2Y₆ receptor activation leads to vasoconstriction. Pancreatic P2Y₁ and P2Y₆ receptors stimulate while P2Y₁₃ receptors inhibits insulin secretion. Antagonists of P2Y₁₂ receptors, and potentially P2Y₁ receptors, are anti-thrombotic agents, and a P2Y₂ /P2Y₄ receptor agonist treats dry eye syndrome in Asia. P2Y receptor agonists are generally pro-inflammatory, and antagonists may eventually treat inflammatory conditions. This article reviews recent developments in P2Y receptor pharmacology (using synthetic agonists and antagonists), structure and biophysical properties (using X-ray crystallography, mutagenesis and modelling), physiological and pathophysiological roles, and present and potentially future therapeutic targeting.

Pharmacology of P2Y receptors

P2Y receptors are G-protein-coupled receptors (GPCRs) for extracellular nucleotides. There are eight mammalian P2Y receptor subtypes divided into two subgroups (P2Y₁, P2Y₂, P2Y₄, P2Y₆, and P2Y₁₁) and (P2Y₁₂, P2Y₁₃, and P2Y₁₄). The P2Y receptors are expressed in various cell types and play important roles in physiology and pathophysiology including inflammatory responses and neuropathic pain. The antagonism of P2Y₁₂ receptors is used in pharmacotherapy for the prevention and therapy of cardiovascular events. The nucleoside analogue ticagrelor and active metabolites of the thienopyridine compounds ticlopidine, clopidogrel and prasugrel inhibit platelet P2Y₁₂ receptors and reduce thereby platelet aggregation. The P2Y₂ receptor agonist diquafosol is used for the treatment of the dry eye syndrome. The P2Y receptor subtypes differ in their amino acid sequences, their pharmacological profiles and their signaling transduction pathways. Recently, selective receptor ligands have been developed for all subtypes. The published crystal structures of the human P2Y₁ and P2Y₁₂ receptors as well as receptor models will facilitate the development of novel drugs for pharmacotherapy.

Studies towards the development of a PET radiotracer for imaging of the P2Y1 receptors in the brain: synthesis, 18F-labeling and preliminary biological evaluation

Purine nucleotides such as ATP and ADP are important extracellular signaling molecules in almost all tissues activating various subtypes of purinoreceptors. In the brain, the P2Y₁ receptor (P2Y₁R) subtype mediates trophic functions like differentiation and proliferation, and modulates fast synaptic transmission, both suggested to be affected in diseases of the central nervous system. Research on P2Y₁R is limited because suitable brain-penetrating P2Y₁R-selective tracers are not yet available. Here, we describe the first efforts to develop an ¹⁸F-labeled PET tracer based on the structure of the highly affine and selective, non-nucleotidic P2Y₁R allosteric modulator 1-(2-[2-(tert-butyl)phenoxy]pyridin-3-yl)-3-[4-(trifluoromethoxy)phenyl]urea (7). A small series of fluorinated compounds was developed by systematic modification of the p-(trifluoromethoxy)phenyl, the urea and the 2-pyridyl subunits of the lead compound 7. Additionally, the p-(trifluoromethoxy)phenyl subunit was substituted by carborane, a boron-rich cluster with potential applicability in boron neutron capture therapy (BNCT). By functional assays, the new fluorinated derivative 1-{2-[2-(tert-butyl)phenoxy]pyridin-3-yl}-3-[4-(2-fluoroethyl)phenyl]urea (18) was identified with a high P2Y₁R antagonistic potency (IC₅₀ = 10 nM). Compound [¹⁸F]18 was radiosynthesized by using tetra-n-butyl ammonium [¹⁸F]fluoride with high radiochemical purity, radiochemical yield and molar activities. Investigation of brain homogenates using hydrophilic interaction chromatography (HILIC) revealed [¹⁸F]fluoride as major radiometabolite. Although [¹⁸F]18 showed fast in vivo metabolization, the high potency and unique allosteric binding mode makes this class of compounds interesting for further optimizations and investigation of the theranostic potential as PET tracer and BNCT agent.

Demystifying P2Y1 Receptor Ligand Recognition through Docking and Molecular Dynamics Analyses

We performed a molecular modeling analysis of 100 nucleotide-like bisphosphates and 46 non-nucleotide arylurea derivatives previously reported as P2Y₁R binders using the recently solved hP2Y₁R structures. We initially docked the compounds at the X-ray structures and identified the binding modes of representative compounds highlighting key patterns in the structure-activity relationship (SAR). We subsequently subjected receptor complexes with selected key agonists (2MeSADP and MRS2268) and antagonists (MRS2500 and BPTU) to membrane molecular dynamics (MD) simulations (at least 200 ns run in triplicate, simulation time 0.6-1.6 μs per ligand system) while considering alternative protonation states of nucleotides. Comparing the temporal evolution of the ligand-protein interaction patterns with available site-directed mutagenesis (SDM) data and P2Y₁R apo state simulation provided further SAR insights and suggested reasonable explanations for loss/gain of binding affinity as well as the most relevant charged species for nucleotide ligands. The MD analysis also predicted local conformational changes required for the receptor inactive state to accommodate nucleotide agonists.

Molecular Recognition of Agonists and Antagonists by the Nucleotide-Activated G Protein-Coupled P2Y2 Receptor

A homology model of the nucleotide-activated P2Y₂R was created based on the X-ray structures of the P2Y₁ receptor. Docking studies were performed, and receptor mutants were created to probe the identified binding interactions. Mutation of residues predicted to interact with the ribose (Arg110) and the phosphates of the nucleotide agonists (Arg265, Arg292) or that contribute indirectly to binding (Tyr288) abolished activity. The Y114F, R194A, and F261A mutations led to inactivity of diadenosine tetraphosphate and to a reduced response of UTP. Significant reduction in agonist potency was observed for all other receptor mutants (Phe111, His184, Ser193, Phe261, Tyr268, Tyr269) predicted to be involved in agonist recognition. An ionic lock between Asp185 and Arg292 that is probably involved in receptor activation interacts with the phosphate groups. The antagonist AR-C118925 and anthraquinones likely bind to the orthosteric site. The updated homology models will be useful for virtual screening and drug design.

Pharmacology and structure of P2Y receptors

P2Y receptors are G-protein-coupled receptors (GPCRs) for extracellular nucleotides. There are eight mammalian P2Y receptor subtypes (P2Y1, P2Y2, P2Y4, P2Y6, P2Y11, P2Y12, P2Y13, and P2Y14). P2Y receptors are widely expressed and play important roles in physiology and pathophysiology. One important example is the ADP-induced platelet aggregation mediated by P2Y1 and P2Y12 receptors. Active metabolites of the thienopyridine compounds ticlopidine, clopidogrel and prasugrel as well as the nucleoside analogue ticagrelor block P2Y12 receptors and thereby platelet aggregation. These drugs are used for the prevention and therapy of cardiovascular events. Moreover, P2Y receptors play important roles in the nervous system. Adenine nucleotides modulate neuronal activity and neuronal fibre outgrowth by activation of P2Y1 receptors and control migration of microglia by P2Y12 receptors. UDP stimulates microglial phagocytosis through activation of P2Y6 receptors. There is evidence for a role for P2Y2 receptors in Alzheimer's disease pathology. The P2Y receptor subtypes are highly diverse in both their amino acid sequences and their pharmacological profiles. Selective receptor ligands have been developed for the pharmacological characterization of the receptor subtypes. The recently published three-dimensional crystal structures of the human P2Y1 and P2Y12 receptors will facilitate the development of therapeutic agents that selectively target P2Y receptors. This article is part of the Special Issue entitled 'Purines in Neurodegeneration and Neuroregeneration'.

Docking and Virtual Screening Strategies for GPCR Drug Discovery

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

Discovery of 4-aryl-7-hydroxyindoline-based P2Y1 antagonists as novel antiplatelet agents

Adenosine diphosphate (ADP)-mediated platelet aggregation is signaled through two distinct G protein-coupled receptors (GPCR) on the platelet surface: P2Y12 and P2Y1. Blocking P2Y12 receptor is a clinically well-validated strategy for antithrombotic therapy. P2Y1 antagonists have been shown to have the potential to provide equivalent antithrombotic efficacy as P2Y12 inhibitors with reduced bleeding in preclinical animal models. We have previously reported the discovery of a potent and orally bioavailable P2Y1 antagonist, 1. This paper describes further optimization of 1 by introducing 4-aryl groups at the hydroxylindoline in two series. In the neutral series, 10q was identified with excellent potency and desirable pharmacokinetic (PK) profile. It also demonstrated similar antithrombotic efficacy with less bleeding compared with the known P2Y12 antagonist prasugrel in rabbit efficacy/bleeding models. In the basic series, 20c (BMS-884775) was discovered with an improved PK and liability profile over 1. These results support P2Y1 antagonism as a promising new antiplatelet target.

4-Benzothiazole-7-hydroxyindolinyl diaryl ureas are potent P2Y1 antagonists with favorable pharmacokinetics: low clearance and small volume of distribution

Current antithrombotic discovery efforts target compounds that are highly efficacious in thrombus reduction with less bleeding liability than the standard of care. Preclinical data suggest that P2Y1 antagonists may have lower bleeding liabilities than P2Y12 antagonists while providing similar antithrombotic efficacy. This article describes our continuous SAR efforts in a series of 7-hydroxyindolinyl diaryl ureas. When dosed orally, 4-trifluoromethyl-7-hydroxy-3,3-dimethylindolinyl analogue 4 was highly efficacious in a model of arterial thrombosis in rats with limited bleeding. The chemically labile CF3 group in 4 was then transformed to various groups via a novel one-step synthesis, yielding a series of potent P2Y1 antagonists. Among them, the 4-benzothiazole-substituted indolines had desirable PK properties in rats, specifically, low clearance and small volume of distribution. In addition, compound 40 had high i.v. exposure and modest bioavailability, giving it the best overall profile.

Conformationally constrained ortho-anilino diaryl ureas: discovery of 1-(2-(1'-neopentylspiro[indoline-3,4'-piperidine]-1-yl)phenyl)-3-(4-(trifluoromethoxy)phenyl)urea, a potent, selective, and bioavailable P2Y1 antagonist

Preclinical antithrombotic efficacy and bleeding models have demonstrated that P2Y1 antagonists are efficacious as antiplatelet agents and may offer a safety advantage over P2Y12 antagonists in terms of reduced bleeding liabilities. In this article, we describe the structural modification of the tert-butyl phenoxy portion of lead compound 1 and the subsequent discovery of a novel series of conformationally constrained ortho-anilino diaryl ureas. In particular, spiropiperidine indoline-substituted diaryl ureas are described as potent, orally bioavailable small-molecule P2Y1 antagonists with improved activity in functional assays and improved oral bioavailability in rats. Homology modeling and rat PK/PD studies on benchmark compound 3l will also be presented. Compound 3l was our first P2Y1 antagonist to demonstrate a robust oral antithrombotic effect with mild bleeding liability in the rat thrombosis and hemostasis models.

Muscarinic receptors as model targets and antitargets for structure-based ligand discovery

G protein-coupled receptors (GPCRs) regulate virtually all aspects of human physiology and represent an important class of therapeutic drug targets. Many GPCR-targeted drugs resemble endogenous agonists, often resulting in poor selectivity among receptor subtypes and restricted pharmacologic profiles. The muscarinic acetylcholine receptor family exemplifies these problems; thousands of ligands are known, but few are receptor subtype-selective and nearly all are cationic in nature. Using structure-based docking against the M2 and M3 muscarinic receptors, we screened 3.1 million molecules for ligands with new physical properties, chemotypes, and receptor subtype selectivities. Of 19 docking-prioritized molecules tested against the M2 subtype, 11 had substantial activity and 8 represented new chemotypes. Intriguingly, two were uncharged ligands with low micromolar to high nanomolar Ki values, an observation with few precedents among aminergic GPCRs. To exploit a single amino-acid substitution among the binding pockets between the M2 and M3 receptors, we selected molecules predicted by docking to bind to the M3 and but not the M2 receptor. Of 16 molecules tested, 8 bound to the M3 receptor. Whereas selectivity remained modest for most of these, one was a partial agonist at the M3 receptor without measurable M2 agonism. Consistent with this activity, this compound stimulated insulin release from a mouse β-cell line. These results support the ability of structure-based discovery to identify new ligands with unexplored chemotypes and physical properties, leading to new biologic functions, even in an area as heavily explored as muscarinic pharmacology.

Structure-based approaches to ligands for G-protein-coupled adenosine and P2Y receptors, from small molecules to nanoconjugates

Adenosine receptor (ARs) and P2Y receptors (P2YRs) that respond to extracellular nucleosides/nucleotides are associated with new directions for therapeutics. The X-ray structures of the A2AAR complexes with agonists and antagonists are examined in relationship to the G-protein-coupled receptor (GPCR) superfamily and applied to drug discovery. Much of the data on AR ligand structure from early SAR studies now are explainable from the A2AAR X-ray crystallography. The ligand-receptor interactions in related GPCR complexes can be identified by means of modeling approaches, e.g., molecular docking. Thus, molecular recognition in binding and activation processes has been studied effectively using homology modeling and applied to ligand design. Virtual screening has yielded new nonnucleoside AR antagonists, and existing ligands have been improved with knowledge of the receptor interactions. New agonists are being explored for central nervous system and peripheral therapeutics based on in vivo activity, such as chronic neuropathic pain. Ligands for receptors more distantly related to the X-ray template, i.e., P2YRs, have been introduced and are mainly used as pharmacological tools for elucidating the physiological role of extracellular nucleotides. Other ligand tools for drug discovery include fluorescent probes, radioactive probes, multivalent probes, and functionalized nanoparticles.

Discovery of 2-(phenoxypyridine)-3-phenylureas as small molecule P2Y1 antagonists

Two distinct G protein-coupled purinergic receptors, P2Y1 and P2Y12, mediate ADP-driven platelet activation. The clinical effectiveness of P2Y12 blockade is well established. Recent preclinical data suggest that P2Y1 and P2Y12 inhibition provide equivalent antithrombotic efficacy, while targeting P2Y1 has the potential for reduced bleeding liability. In this account, the discovery of a 2-(phenoxypyridine)-3-phenylurea chemotype that inhibited ADP-mediated platelet aggregation in human blood samples is described. Optimization of this series led to the identification of compound 16, 1-(2-(2-tert-butylphenoxy)pyridin-3-yl)-3-4-(trifluoromethoxy)phenylurea, which demonstrated a 68 ± 7% thrombus weight reduction in an established rat arterial thrombosis model (10 mg/kg plus 10 mg/kg/h) while only prolonging cuticle and mesenteric bleeding times by 3.3- and 3.1-fold, respectively, in provoked rat bleeding time models. These results suggest that a P2Y1 antagonist could potentially provide a safe and efficacious antithrombotic profile.

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).