The molecular determinants of small-molecule ligand binding at P2X receptors

P2X receptors are trimeric eukaryotic ATP-gated cation channels. Extracellular ATP - their physiological ligand - is released as a neurotransmitter and in conditions of cell damage such as inflammation, and substantial evidence implicates P2X receptors in diseases including neuropathic pain, cancer and arthritis. In 2009, the first P2X crystal structure, Danio rerio P2X4 in the apo- state, was published, and this was followed in 2012 by the ATP-bound structure. These structures transformed our understanding of the conformational changes induced by ATP binding and the mechanism of ligand specificity, and enabled homology modelling of mammalian P2X receptors for ligand docking and rational design of receptor modulators. P2X receptors are attractive drug targets, and a wide array of potent, subtype-selective modulators (mostly antagonists) have been developed. In 2016, crystal structures of human P2X3 in complex with the competitive antagonists TNP-ATP and A-317491, and Ailuropoda melanoleuca P2X7 in complex with a series of allosteric antagonists were published, giving fascinating insights into the mechanism of channel antagonism. In this article we not only summarize current understanding of small-molecule modulator binding at P2X receptors, but also use this information in combination with previously published structure-function data and molecular docking experiments, to hypothesise a role for the dorsal fin loop region in differential ATP potency, and describe novel, testable binding conformations for both the semi-selective synthetic P2X7 agonist 2’-(3’)-O-(4-benzoyl)benzoyl ATP (BzATP), and the P2X4-selective positive allosteric modulator ivermectin. We find that the distal benzoyl group of BzATP lies in close proximity to Lys-127, a residue previously implicated in BzATP binding to P2X7, potentially explaining the increased potency of BzATP at rat P2X7 receptors. Furthermore, to our knowledge, we present the first molecular docking of ivermectin to rat P2X4 receptors, illustrating a plausible binding conformation between the first and second transmembrane domains which not only tallies with previous mutagenesis studies, but would also likely have the effect of stabilizing the open channel structure, consistent with the mode of action of this positive allosteric modulator. From our docking simulations and analysis of sequence homology we propose a series of mutations likely to confer ivermectin sensitivity to human P2X1

Kinase-independent phosphoramidate S1P1 receptor agonist benzyl ether derivatives

Previously published S1P receptor modulator benzyl ether derivatives have shown potential as being viable therapeutics for the treatment of neurodegenerative diseases, however, two of the most S1P1-selective compounds are reported as being poorly phosphorylated by kinases in vivo. Phosphoramidates of BED compounds (2a, 2b) were synthesised with the aim of producing kinase-independent S1P receptor modulators. Carboxypeptidase, human serum and cell lysate processing experiments were conducted. ProTide BED analogues were found to have an acceptable level of stability in acidic and basic conditions and in vitro metabolic processing experiments showed that they are processed to the desired pharmacologically active monophosphate. The research describes the development of an entirely novel family of therapeutic agents

Erratum to ‘Combining bioinformatics, cheminformatics, functional genomics and whole organism approaches for identifying epigenetic drug targets in Schistosoma mansoni’:[IJP Drugs Drug Resist. 8 (2018) 559–570]

Schistosomiasis endangers the lives of greater than 200 million people every year and is predominantly controlled by a single class chemotherapy, praziquantel (PZQ). Development of PZQ replacement (to combat the threat of PZQ insensitivity/resistance arising) or combinatorial (to facilitate the killing of PZQ-insensitive juvenile schistosomes) chemotherapies would help sustain this control strategy into the future. Here, we re-categorise two families of druggable epigenetic targets in Schistosoma mansoni, the histone methyltransferases (HMTs) and the histone demethylases (HDMs). Amongst these, a S. mansoni Lysine Specific Demethylase 1 (SmLSD1, Smp_150560) homolog was selected for further analyses. Homology modelling of SmLSD1 and in silico docking of greater than four thousand putative inhibitors identified seven (L1 – L7) showing more favourable binding to the target pocket of SmLSD1 vs Homo sapiens HsLSD1; six of these seven (L1 – L6) plus three structural analogues of L7 (L8 – L10) were subsequently screened against schistosomula using the Roboworm anthelmintic discovery platform. The most active compounds (L10 - pirarubicin > L8 – danunorubicin hydrochloride) were subsequently tested against juvenile (3 wk old) and mature (7 wk old) schistosome stages and found to impede motility, suppress egg production and affect tegumental surfaces. When compared to a surrogate human cell line (HepG2), a moderate window of selectivity was observed for the most active compound L10 (selectivity indices - 11 for schistosomula, 9 for juveniles, 1.5 for adults). Finally, RNA interference of SmLSD1 recapitulated the egg-laying defect of schistosomes co-cultivated in the presence of L10 and L8. These preliminary results suggest that SmLSD1 represents an attractive new target for schistosomiasis; identification of more potent and selective SmLSD1 compounds, however, is essential. Nevertheless, the approaches described herein highlight an interdisciplinary strategy for selecting and screening novel/repositioned anti-schistosomals, which can be applied to any druggable (epigenetic) target encoded by the parasite's genome. Keywords: Anthelmintic drug discovery, Neglected tropical diseases, Schistosoma mansoni, Epigenetics, Lysine specific demethylas

Haptic-driven, interactive drug design: implementing a GPU-based approach to evaluate the induced fit effect

In this paper, we describe a hybrid meta-heuristics method of energy minimization and conformational sampling and its application into our haptic-driven molecular modelling simulator. The proposed method has been designed to suit real-time molecular docking simulations, where the time-lapse between two successive ligand poses is relatively small. In these situations, the energy minimization problem becomes increasingly complex and chaotic. The algorithm is tuned to take advantage of recent advances in GPU computing with asynchronous kernel execution, which has allowed us to include full protein flexibility in the real-time interactive, haptic-driven simulations. Finally, in this paper, we will also discuss the implementation of such high-performance computing approaches in our software, discussing the results of our initial validation studies, highlighting the advantages and limitations of such interactive methodology

Homology Modelling of Human E1 Ubiquitin Activating Enzyme

Human E1 is a key player in protein ubiquitination, however the E1 structure is not available. In this paper, we describe the derivation of a human E1 structure using molecular modelling based on the crystal structure of S. cerevisiae E1 and M. Musculus E1. Key interactions between our E1 model and ubiquitin are also discussed

Rational modifications, synthesis and biological evaluation of new potential antivirals for RSV designed to target the M2-1 protein

Respiratory syncytial virus (RSV) is the main cause of lower respiratory tract diseases in infants and young children, with potentially serious and fatal consequences associated with severe infections. Despite extensive research efforts invested in the identification of therapeutic measures, no vaccine is currently available, while treatment options are limited to ribavirin and palivizumab, which both present significant limitations. While clinical and pre-clinical candidates mainly target the viral fusion protein, the nucleocapsid protein or the viral polymerase, our focus has been the identification of new antiviral compounds targeting the viral M2-1 protein, thanks to the presence of a zinc-ejecting group in their chemical structure. Starting from an anti-RSV hit we had previously identified with an in silico structure-based approach, we have designed, synthesised and evaluated a new series of dithiocarbamate analogues, with which we have explored the antiviral activity of this scaffold. The findings presented in this work may provide the basis for the identification of a new antiviral lead to treat RSV infections

Identification of a novel P2X7 antagonist using structure-based virtual screening

P2X4 and P2X7 receptors are ATP-gated ion channels, which play important roles in neuropathic and inflammatory pain, and as such they are important drug targets in diseases of inflammatory origin. While several compounds targeting P2X4 and P2X7 receptors have been developed using traditional high-throughput screening approaches, relatively few compounds have been developed using structure-based design. We initially set out to develop compounds targeting human P2X4, by performing virtual screening on the orthosteric (ATP-binding) pocket of a molecular model of human P2X4 based on the crystal structure of the Danio rerio receptor. The screening of a library of approximately 300,000 commercially available drug-like compounds led to the initial selection of 17 compounds; however, none of these compounds displayed a significant antagonist effect at P2X4 in a Fluo-4 ATP-induced calcium influx assay. When the same set of compounds was tested against human P2X7 in an ATP-stimulated Yo-Pro1 dye uptake assay, one compound (an indeno(1,2-b)pyridine derivative; GP-25) reduced the response by greater than 50% when applied at a concentration of 30 µM. GP-25 displayed an IC50 value of 8.7 μM at human P2X7 and 24.4 μM at rat P2X7, and was confirmed to be active using whole-cell patch clamp electrophysiology and not cytotoxic. Schild analysis suggested that mode of action of GP-25 was orthosteric. Screening of a further 16 commercially available analogues of GP-25 led to the discovery of five additional compounds with antagonist activity at human P2X7, enabling us to investigate the structure-activity relationship. Finally, docking of the R- and S-enantiomers of GP-25 into the orthosteric pocket of molecular models of human P2X4 and human P2X7 revealed that, while both enantiomers were able to make multiple interactions between their carboxyl moieties and conserved positively charged amino-acids in human P2X7, only the S-enantiomer of GP-25 was able to do this in human P2X4, potentially explaining the lack of activity of GP-25 at this receptor

Identification of the molecular determinants of antagonist potency in the allosteric binding pocket of human P2X4

P2X receptors are a family of ATP-gated cation channels comprising seven subtypes in mammals, which play key roles in nerve transmission, pain sensation and inflammation. The P2X4 receptor in particular has attracted significant interest from pharmaceutical companies due to its physiological roles in neuropathic pain and modulation of vascular tone. A number of potent small-molecule P2X4 receptor antagonists have been developed, including the allosteric P2X4 receptor antagonist BX430, which is approximately 30-fold more potent at human P2X4 compared with the rat isoform. A single amino-acid difference between human and rat P2X4 (I312T), located in an allosteric pocket, has previously been identified as critical for BX430 sensitivity, implying that BX430 binds in this pocket. Using a combination of mutagenesis, functional assay in mammalian cells and in silico docking we confirmed these findings. Induced-fit docking, permitting the sidechains of the amino-acids of P2X4 to move, showed that BX430 could access a deeper portion of the allosteric pocket, and that the sidechain of Lys-298 was important for shaping the cavity. We then performed blind docking of 12 additional P2X4 antagonists into the receptor extracellular domain, finding that many of these compounds favored the same pocket as BX430 from their calculated binding energies. Induced-fit docking of these compounds in the allosteric pocket enabled us to show that antagonists with high potency (IC50 ≤ 100 nM) bind deep in the allosteric pocket, disrupting a network of interacting amino acids including Asp-85, Ala-87, Asp-88, and Ala-297, which are vital for transmitting the conformational change following ATP binding to channel gating. Our work confirms the importance of Ile-312 for BX430 sensitivity, demonstrates that the allosteric pocket where BX430 binds is a plausible binding pocket for a series of P2X4 antagonists, and suggests a mode of action for these allosteric antagonists involving disruption of a key structural motif required for the conformational change induced in P2X4 when ATP binds

A rational design of α-helix-shaped peptides employing the hydrogen-bond surrogate approach: A modulation strategy for Ras-RasGRF1 interaction in neuropsychiatric disorders

In the last two decades, abnormal Ras (rat sarcoma protein)–ERK (extracellular signal-regulated kinase) signalling in the brain has been involved in a variety of neuropsychiatric disorders, including drug addiction, certain forms of intellectual disability, and autism spectrum disorder. Modulation of membrane-receptor-mediated Ras activation has been proposed as a potential target mechanism to attenuate ERK signalling in the brain. Previously, we showed that a cell penetrating peptide, RB3, was able to inhibit downstream signalling by preventing RasGRF1 (Ras guanine nucleotide-releasing factor 1), a neuronal specific GDP/GTP exchange factor, to bind Ras proteins, both in brain slices and in vivo, with an IC50 value in the micromolar range. The aim of this work was to mutate and improve this peptide through computer-aided techniques to increase its inhibitory activity against RasGRF1. The designed peptides were built based on the RB3 peptide structure corresponding to the α-helix of RasGRF1 responsible for Ras binding. For this purpose, the hydrogen-bond surrogate (HBS) approach was exploited to maintain the helical conformation of the designed peptides. Finally, residue scanning, MD simulations, and MM-GBSA calculations were used to identify 18 most promising α-helix-shaped peptides that will be assayed to check their potential activity against Ras-RasGRF1 and prevent downstream molecular events implicated in brain disorders