David A. Evans

Crystal structures of the M1 and M4 muscarinic acetylcholine receptors

Muscarinic M1–M5 acetylcholine receptors are G-protein-coupled receptors that regulate many vital functions of the central and peripheral nervous systems. In particular, the M1 and M4 receptor subtypes have emerged as attractive drug targets for treatments of neurological disorders, such as Alzheimer’s disease and schizophrenia, but the high conservation of the acetylcholine-binding pocket has spurred current research into targeting allosteric sites on these receptors. Here we report the crystal structures of the M1 and M4 muscarinic receptors bound to the inverse agonist, tiotropium. Comparison of these structures with each other, as well as with the previously reported M2 and M3 receptor structures, reveals differences in the orthosteric and allosteric binding sites that contribute to a role in drug selectivity at this important receptor family. We also report identification of a cluster of residues that form a network linking the orthosteric and allosteric sites of the M4 receptor, which provides new insight into how allosteric modulation may be transmitted between the two spatially distinct domains.

Extensions to In Silico Bioactivity Predictions Using Pathway Annotations and Differential Pharmacology Analysis: Application to Xenopus laevis Phenotypic Readouts

The simultaneous increase of computational power and the availability of chemical and biological data have contributed to the recent popularity of in silico bioactivity prediction algorithms. Such methods are commonly used to infer the ‘Mechanism of Action’ of small molecules and they can also be employed in cases where full bioactivity profiles have not been established experimentally. However, protein target predictions by themselves do not necessarily capture information about the effect of a compound on a biological system, and hence merging their output with a systems biology approach can help to better understand the complex network modulation which leads to a particular phenotype. In this work, we review approaches and applications of target prediction, as well as their shortcomings, and demonstrate two extensions of this concept which are exemplified using phenotypic readouts from a chemical genetic screen in Xenopus laevis. In particular, the experimental observations are linked to their predicted bioactivity profiles. Predicted targets are annotated with pathways, which lead to further biological insight. Moreover, we subject the prediction to further machine learning algorithms, namely decision trees, to capture the differential pharmacology of ligand‐target interactions in biological systems. Both methodologies hence provide new insight into understanding the Mechanism of Action of compound activities from phenotypic screens.

Extending in silico mechanism-of- action analysis by annotating targets with pathways: application to cellular cytotoxicity readouts

BACKGROUND An in silico mechanism-of-action analysis protocol was developed, comprising molecule bioactivity profiling, annotation of predicted targets with pathways and calculation of enrichment factors to highlight targets and pathways more likely to be implicated in the studied phenotype. RESULTS The method was applied to a cytotoxicity phenotypic endpoint, with enriched targets/pathways found to be statistically significant when compared with 100 random datasets. Application on a smaller apoptotic set (10 molecules) did not allowed to obtain statistically relevant results, suggesting that the protocol requires modification such as analysis of the most frequently predicted targets/annotated pathways. CONCLUSION Pathway annotations improved the mechanism-of-action information gained by target prediction alone, allowing a better interpretation of the predictions and providing better mapping of targets onto pathways.

ChEMBLSpace—a graphical explorer of the chemogenomic space covered by the ChEMBL database

MOTIVATION The ChEMBLSpace graphical explorer enables the identification of compounds from the ChEMBL database, which exhibit a desirable polypharmacology profile. This profile can be predefined or created iteratively, and the tool can be extended to other data sources.

Comparing global and local likelihood score thresholds in multiclass laplacian-modified Naive Bayes protein target prediction.

The increase of publicly available bioactivity data has led to the extensive development and usage of in silico bioactivity prediction algorithms. A particularly popular approach for such analyses is the multiclass Naïve Bayes, whose output is commonly processed by applying empirically-derived likelihood score thresholds. In this work, we describe a systematic way for deriving score cut-offs on a per-protein target basis and compare their performance with global thresholds on a large scale using both 5-fold cross-validation (ChEMBL 14, 189k ligand-protein pairs over 477 protein targets) and external validation (WOMBAT, 63k pairs, 421 targets). The individual protein target cut-offs derived were compared to global cut-offs ranging from -10 to 40 in score bouts of 2.5. The results indicate that individual thresholds had equal or better performance in all comparisons with global thresholds, ranging from 95% of protein targets to 57.96%. It is shown that local thresholds behave differently for particular families of targets (CYPs, GPCRs, Kinases and TFs). Furthermore, we demonstrate the discrepancy in performance when we move away from the training dataset chemical space, using Tanimoto similarity as a metric (from 0 to 1 in steps of 0.2). Finally, the individual protein score cut-offs derived for the in silico bioactivity application used in this work are released, as well as the reproducible and transferable KNIME workflows used to carry out the analysis.

Predicting the Accuracy of Ligand Overlay Methods with Random Forest Models

The accuracy of binding mode prediction using standard molecular overlay methods (ROCS, FlexS, Phase, and FieldCompare) is studied. Previous work has shown that simple decision tree modeling can be used to improve accuracy by selection of the best overlay template. This concept is extended to the use of Random Forest (RF) modeling for template and algorithm selection. An extensive data set of 815 ligand-bound X-ray structures representing 5 gene families was used for generating ca. 70,000 overlays using four programs. RF models, trained using standard measures of ligand and protein similarity and Lipinski-related descriptors, are used for automatically selecting the reference ligand and overlay method maximizing the probability of reproducing the overlay deduced from X-ray structures (i.e., using rmsd < or = 2 A as the criteria for success). RF model scores are highly predictive of overlay accuracy, and their use in template and method selection produces correct overlays in 57% of cases for 349 overlay ligands not used for training RF models. The inclusion in the models of protein sequence similarity enables the use of templates bound to related protein structures, yielding useful results even for proteins having no available X-ray structures.

Discovery of selective RIO2 kinase small molecule ligand.

We report the discovery and initial optimization of diphenpyramide and several of its analogs as hRIO2 kinase ligands. One of these analogs is the most selective hRIO2 ligand reported to date. Diphenpyramide is a Cyclooxygenase 1 and 2 inhibitor that was used as an anti-inflammatory agent. The RIO2 kinase affinity of diphenpyramide was discovered by serendipity while profiling of 13 marketed drugs on a large 456 kinase assay panel. The inhibition values also suggested a relative selectivity of diphenpyramide for RIO2 against the other kinases in the panel. Subsequently three available and eight newly synthesized analogs were assayed, one of which showed a 10 fold increased hRIO2 binding affinity. Additionally, this compound shows significantly better selectivity over assayed kinases, when compared to currently known RIO2 inhibitors. As RIO2 is involved in the biosynthesis of the ribosome and cell cycle regulation, our selective ligand may be useful for the delineation of the biological role of this kinase. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.

Design and synthesis of N-[6-(Substituted Aminoethylideneamino)-2-Hydroxyindan-1-yl]arylamides as selective and potent muscarinic M1 agonists

The observation that cholinergic deafferentation of circuits projecting from forebrain basal nuclei to frontal and hippocampal circuits occurs in Alzheimers disease has led to drug-targeting of muscarinic M1 receptors to alleviate cognitive symptoms. The high homology within the acetylcholine binding domain of this family however has made receptor-selective ligand development challenging. This work presents the synthesis scheme, pharmacokinetic and structure-activity-relationship study findings for M1-selective ligand, LY593093. Pharmacologically the compound acts as an orthosteric ligand. The homology modeling work presented however will illustrate that compound binding spans from the acetylcholine pocket to the extracellular loops of the receptor, a common allosteric vestibule for the muscarinic protein family. Altogether LY593093 represents a growing class of multi-topic ligands which interact with the receptors in both the ortho- and allosteric binding sites, but which exert their activation mechanism as an orthosteric ligand.

Comparative mode-of-action analysis following manual and automated phenotype detection in Xenopus laevis

Given the increasing utilization of phenotypic screens in drug discovery also the subsequent mechanism-of-action analysis gains increased attention. Such analyses frequently use in silico methods, which have become significantly more popular in recent years. However, identifying phenotype-specific mechanisms of action depends heavily on suitable phenotype identification in the first place, many of which rely on human input and are therefore inconsistent. In this work, we aimed at analysing the impact that human phenotype classification has on subsequent in silico mechanism-of-action analysis. To this end, an image analysis application was implemented for the rapid identification of seven high-level phenotypes in Xenopus laevis tadpoles treated with compounds from the National Cancer Institute Diversity Set II. It was found that manual and automated phenotype classifications were in agreement with some of the phenotypes (e.g. 73.9% agreement observed for general morphology abnormality), while this was not the case in others (e.g. melanophore migration with 37.6% agreement between both annotations). Based on both annotations, protein targets of active compounds were predicted in silico, and decision trees were generated to understand mechanisms-of-action behind every phenotype while also taking polypharmacology (combinations of targets) into account. It was found that the automated phenotype categorisation greatly increased the accuracy of the results of the mechanism-of-action model, where it improved the classification accuracy by 9.4%, as well as reducing the tree size by eight nodes and the number of leaves and the depth by three levels. Overall we conclude that consistent phenotype annotations seem to be generally crucial for successful subsequent mechanism-of-action analysis, and this is what we have shown here in Xenopus laevis screens in combination with in silico mechanism-of-action analysis.

Expression-based drug screening of neural progenitor cells from individuals with schizophrenia

A lack of biologically relevant screening models hinders the discovery of better treatments for schizophrenia (SZ) and other neuropsychiatric disorders. Here we compare the transcriptional responses of 8 commonly used cancer cell lines (CCLs) directly with that of human induced pluripotent stem cell (hiPSC)-derived neural progenitor cells (NPCs) from 12 individuals with SZ and 12 controls across 135 drugs, generating 4320 unique drug-response transcriptional signatures. We identify those drugs that reverse post-mortem SZ-associated transcriptomic signatures, several of which also differentially regulate neuropsychiatric disease-associated genes in a cell type (hiPSC NPC vs. CCL) and/or a diagnosis (SZ vs. control)-dependent manner. Overall, we describe a proof-of-concept application of transcriptomic drug screening to hiPSC-based models, demonstrating that the drug-induced gene expression differences observed with patient-derived hiPSC NPCs are enriched for SZ biology, thereby revealing a major advantage of incorporating cell type and patient-specific platforms in drug discovery.Unbiased large scale screening of small molecules for drug discovery in psychiatric disease is technically challenging and financially costly. Here, Readhead and colleagues integrate in silico and in vitro approaches to design and conduct transcriptomic drug screening in schizophrenia patient-derived neural cells, in order to survey novel pathologies and points of intervention.