Polypharmacology prediction: the long road toward comprehensively anticipating small-molecule selectivity to de-risk drug discovery.

INTRODUCTION: Small molecules often bind to multiple targets, a behavior termed polypharmacology. Anticipating polypharmacology is essential for drug discovery since unknown off-targets can modulate safety and efficacy - profoundly affecting drug discovery success. Unfortunately, experimental methods to assess selectivity present significant limitations and drugs still fail in the clinic due to unanticipated off-targets. Computational methods are a cost-effective, complementary approach to predict polypharmacology. AREAS COVERED: This review aims to provide a comprehensive overview of the state of polypharmacology prediction and discuss its strengths and limitations, covering both classical cheminformatics methods and bioinformatic approaches. The authors review available data sources, paying close attention to their different coverage. The authors then discuss major algorithms grouped by the types of data that they exploit using selected examples. EXPERT OPINION: Polypharmacology prediction has made impressive progress over the last decades and contributed to identify many off-targets. However, data incompleteness currently limits most approaches to comprehensively predict selectivity. Moreover, our limited agreement on model assessment challenges the identification of the best algorithms - which at present show modest performance in prospective real-world applications. Despite these limitations, the exponential increase of multidisciplinary Big Data and AI hold much potential to better polypharmacology prediction and de-risk drug discovery.

Computational repurposing of oncology drugs through off-target drug binding interactions from pharmacological databases.

PURPOSE: Systematic repurposing of approved medicines for another indication may accelerate drug development in oncology. We present a strategy combining biomarker testing with drug repurposing to identify new treatments for patients with advanced cancer. METHODS: Tumours were sequenced with the Illumina TruSight Oncology 500 (TSO-500) platform or the FoundationOne CDx panel. Mutations were screened by two medical oncologists and pathogenic mutations were categorised referencing literature. Variants of unknown significance were classified as potentially pathogenic using plausible mechanisms and computational prediction of pathogenicity. Gain of function (GOF) mutations were evaluated through repurposing databases Probe Miner (PM), Broad Institute Drug Repurposing Hub (Broad Institute DRH) and TOPOGRAPH. GOF mutations were repurposing events if identified in PM, not indexed in TOPOGRAPH and excluding mutations with a known Food and Drug Administration (FDA)-approved biomarker. The computational repurposing approach was validated by evaluating its ability to identify FDA-approved biomarkers. The total repurposable genome was identified by evaluating all possible gene-FDA drug-approved combinations in the PM dataset. RESULTS: The computational repurposing approach was accurate at identifying FDA therapies with known biomarkers (94%). Using next-generation sequencing molecular reports (n = 94), a meaningful percentage of patients (14%) could have an off-label therapeutic identified. The frequency of theoretical drug repurposing events in The Cancer Genome Atlas pan-cancer dataset was 73% of the samples in the cohort. CONCLUSION: A computational drug repurposing approach may assist in identifying novel repurposing events in cancer patients with no access to standard therapies. Further validation is needed to confirm a precision oncology approach using drug repurposing.

Identification of differential biological activity and synergy between the PARP inhibitor rucaparib and its major metabolite.

The (poly)pharmacology of drug metabolites is seldom comprehensively characterized in drug discovery. However, some drug metabolites can reach high plasma concentrations and display in vivo activity. Here, we use computational and experimental methods to comprehensively characterize the kinase polypharmacology of M324, the major metabolite of the PARP1 inhibitor rucaparib. We demonstrate that M324 displays unique PLK2 inhibition at clinical concentrations. This kinase activity could have implications for the efficacy and safety of rucaparib and therefore warrants further clinical investigation. Importantly, we identify synergy between the drug and the metabolite in prostate cancer models and a complete reduction of α-synuclein accumulation in Parkinson's disease models. These activities could be harnessed in the clinic or open new drug discovery opportunities. The study reported here highlights the importance of characterizing the activity of drug metabolites to comprehensively understand drug response in the clinic and exploit our current drug arsenal in precision medicine.

Short-term exposure to environmental levels of nicotine and cotinine impairs visual motor response in zebrafish larvae through a similar mode of action: Exploring the potential role of zebrafish α7 nAChR.

The current view is that environmental levels of nicotine and cotinine, commonly in the ng/L range, are safe for aquatic organisms. In this study, 7 days post-fertilization zebrafish embryos have been exposed for 24 h to a range of environmental concentrations of nicotine (2.0 ng/L-2.5 µg/L) and cotinine (50 pg/L-10 µg/L), as well as to a binary mixture of these emerging pollutants. Nicotine exposure led to hyperactivity, decreased vibrational startle response and increased non-associative learning. However, the more consistent effect found for both nicotine and cotinine was a significant increase in light-off visual motor response (VMR). The effect of both pollutants on this behavior occurred through a similar mode of action, as the joint effects of the binary mixture of both chemicals were consistent with the concentration addition concept predictions. The results from docking studies suggest that the effect of nicotine and cotinine on light-off VMR could be mediated by zebrafish α7 nAChR expressed in retina. The results presented in this study emphasize the need to revisit the environmental risk assessment of chemicals including additional ecologically relevant sublethal endpoints.

Elucidating Compound Mechanism of Action and Polypharmacology with a Large-scale Perturbational Profile Compendium.

The Mechanism of Action (MoA) of a drug is generally represented as a small, non-tissue-specific repertoire of high-affinity binding targets. Yet, drug activity and polypharmacology are increasingly associated with a broad range of off-target and tissue-specific effector proteins. To address this challenge, we have implemented an efficient integrative experimental and computational framework leveraging the systematic generation and analysis of drug perturbational profiles representing >700 FDA-approved and experimental oncology drugs, in cell lines selected as high-fidelity models of 23 aggressive tumor subtypes. Protein activity-based analyses revealed highly reproducible, drug-mediated modulation of tissue-specific targets, leading to generation of a proteome-wide polypharmacology map, characterization of MoA-related drug clusters and off-target effects, and identification and experimental validation of novel, tissue-specific inhibitors of undruggable oncoproteins. The proposed framework, which is easily extended to elucidating the MoA of novel small-molecule libraries, could help support more systematic and quantitative approaches to precision oncology.

A machine learning and live-cell imaging tool kit uncovers small molecules induced phospholipidosis.

Drug-induced phospholipidosis (DIPL), characterized by excessive accumulation of phospholipids in lysosomes, can lead to clinical adverse effects. It may also alter phenotypic responses in functional studies using chemical probes. Therefore, robust methods are needed to predict and quantify phospholipidosis (PL) early in drug discovery and in chemical probe characterization. Here, we present a versatile high-content live-cell imaging approach, which was used to evaluate a chemogenomic and a lysosomal modulation library. We trained and evaluated several machine learning models using the most comprehensive set of publicly available compounds and interpreted the best model using SHapley Additive exPlanations (SHAP). Analysis of high-quality chemical probes extracted from the Chemical Probes Portal using our algorithm revealed that closely related molecules, such as chemical probes and their matched negative controls can differ in their ability to induce PL, highlighting the importance of identifying PL for robust target validation in chemical biology.

Target 2035 - an update on private sector contributions.

Target 2035, an international federation of biomedical scientists from the public and private sectors, is leveraging 'open' principles to develop a pharmacological tool for every human protein. These tools are important reagents for scientists studying human health and disease and will facilitate the development of new medicines. It is therefore not surprising that pharmaceutical companies are joining Target 2035, contributing both knowledge and reagents to study novel proteins. Here, we present a brief progress update on Target 2035 and highlight some of industry's contributions.

canSAR: update to the cancer translational research and drug discovery knowledgebase.

canSAR (https://cansar.ai) is the largest public cancer drug discovery and translational research knowledgebase. Now hosted in its new home at MD Anderson Cancer Center, canSAR integrates billions of experimental measurements from across molecular profiling, pharmacology, chemistry, structural and systems biology. Moreover, canSAR applies a unique suite of machine learning algorithms designed to inform drug discovery. Here, we describe the latest updates to the knowledgebase, including a focus on significant novel data. These include canSAR's ligandability assessment of AlphaFold; mapping of fragment-based screening data; and new chemical bioactivity data for novel targets. We also describe enhancements to the data and interface.

The Chemical Probes Portal: an expert review-based public resource to empower chemical probe assessment, selection and use.

We describe the Chemical Probes Portal (https://www.chemicalprobes.org/), an expert review-based public resource to empower chemical probe assessment, selection and use. Chemical probes are high-quality small-molecule reagents, often inhibitors, that are important for exploring protein function and biological mechanisms, and for validating targets for drug discovery. The publication, dissemination and use of chemical probes provide an important means to accelerate the functional annotation of proteins, the study of proteins in cell biology, physiology, and disease pathology, and to inform and enable subsequent pioneering drug discovery and development efforts. However, the widespread use of small-molecule compounds that are claimed as chemical probes but are lacking sufficient quality, especially being inadequately selective for the desired target or even broadly promiscuous in behaviour, has resulted in many erroneous conclusions in the biomedical literature. The Chemical Probes Portal was established as a public resource to aid the selection and best-practice use of chemical probes in basic and translational biomedical research. We describe the background, principles and content of the Portal and its technical development, as well as examples of its applications and use. The Chemical Probes Portal is a community resource and we therefore describe how researchers can be involved in its content and development.

canSAR chemistry registration and standardization pipeline.

BACKGROUND: Integration of medicinal chemistry data from numerous public resources is an increasingly important part of academic drug discovery and translational research because it can bring a wealth of important knowledge related to compounds in one place. However, different data sources can report the same or related compounds in various forms (e.g., tautomers, racemates, etc.), thus highlighting the need of organising related compounds in hierarchies that alert the user on important bioactivity data that may be relevant. To generate these compound hierarchies, we have developed and implemented canSARchem, a new compound registration and standardization pipeline as part of the canSAR public knowledgebase. canSARchem builds on previously developed ChEMBL and PubChem pipelines and is developed using KNIME. We describe the pipeline which we make publicly available, and we provide examples on the strengths and limitations of the use of hierarchies for bioactivity data exploration. Finally, we identify canonicalization enrichment in FDA-approved drugs, illustrating the benefits of our approach. RESULTS: We created a chemical registration and standardization pipeline in KNIME and made it freely available to the research community. The pipeline consists of five steps to register the compounds and create the compounds' hierarchy: 1. Structure checker, 2. Standardization, 3. Generation of canonical tautomers and representative structures, 4. Salt strip, and 5. Generation of abstract structure to generate the compound hierarchy. Unlike ChEMBL's RDKit pipeline, we carry out compound canonicalization ahead of getting the parent structure, similar to PubChem's OpenEye pipeline. canSARchem has a lower rejection rate compared to both PubChem and ChEMBL. We use our pipeline to assess the impact of grouping the compounds in hierarchies for bioactivity data exploration. We find that FDA-approved drugs show statistically significant sensitivity to canonicalization compared to the majority of bioactive compounds which demonstrates the importance of this step. CONCLUSIONS: We use canSARchem to standardize all the compounds uploaded in canSAR (> 3 million) enabling efficient data integration and the rapid identification of alternative compound forms with useful bioactivity data. Comparison with PubChem and ChEMBL pipelines evidenced comparable performances in compound standardization, but only PubChem and canSAR canonicalize tautomers and canSAR has a slightly lower rejection rate. Our results highlight the importance of compound hierarchies for bioactivity data exploration. We make canSARchem available under a Creative Commons Attribution-ShareAlike 4.0 International License (CC BY-SA 4.0) at https://gitlab.icr.ac.uk/cansar-public/compound-registration-pipeline .

Rare Variation in Drug Metabolism and Long QT Genes and the Genetic Susceptibility to Acquired Long QT Syndrome.

BACKGROUND: Acquired long QT syndrome (aLQTS) is a serious unpredictable adverse drug reaction. Pharmacogenomic markers may predict risk. METHODS: Among 153 aLQTS patients (mean age 58 years [range, 14-88], 98.7% White, 85.6% symptomatic), computational methods identified proteins interacting most significantly with 216 QT-prolonging drugs. All cases underwent sequencing of 31 candidate genes arising from this analysis or associating with congenital LQTS. Variants were filtered using a minor allele frequency <1% and classified for susceptibility for aLQTS. Gene-burden analyses were then performed comparing the primary cohort to control exomes (n=452) and an independent replication aLQTS exome sequencing cohort. RESULTS: In 25.5% of cases, at least one rare variant was identified: 22.2% of cases carried a rare variant in a gene associated with congenital LQTS, and in 4% of cases that variant was known to be pathogenic or likely pathogenic for congenital LQTS; 7.8% cases carried a cytochrome-P450 (CYP) gene variant. Of 12 identified CYP variants, 11 (92%) were in an enzyme known to metabolize at least one culprit drug to which the subject had been exposed. Drug-drug interactions that affected culprit drug metabolism were found in 19% of cases. More than one congenital LQTS variant, CYP gene variant, or drug interaction was present in 7.8% of cases. Gene-burden analyses of the primary cohort compared to control exomes (n=452), and an independent replication aLQTS exome sequencing cohort (n=67) and drug-tolerant controls (n=148) demonstrated an increased burden of rare (minor allele frequency<0.01) variants in CYP genes but not LQTS genes. CONCLUSIONS: Rare susceptibility variants in CYP genes are emerging as potentially important pharmacogenomic risk markers for aLQTS and could form part of personalized medicine approaches in the future.

Identification of different side effects between PARP inhibitors and their polypharmacological multi-target rationale.

AIMS: The aim of this study was to determine the differences and potential mechanistic rationale for observed adverse drug reactions (ADRs) between four approved PARP inhibitors (PARPi). METHODS: The Medicines and Healthcare products Regulatory Authority (MHRA) Yellow Card drug analysis profiles and NHS secondary care medicines database enabled the identification of suspected ADRs associated with the PARPi in the UK from launch to 2020. The polypharmacology of the PARPi were data-mined from several public data sources. RESULTS: The overall ADRs per 100 000 Rx identified across the four PARPi are statistically significant (χ2 test, P < .001). Rucaparib has the greatest relative suspected ADRs, which can be explained by its least clean kinome and physicochemical properties. The suspected gastrointestinal ADRs of rucaparib and niraparib can be ascribed to their kinase polypharmacology. Suspected blood and lymphatic system ADRs of PARPi can be linked to their high volume of distribution (Vd ). The thrombocytopenia rate of niraparib > rucaparib > olaparib tracked with the Vd trend. Hypertension is only associated with niraparib and could be explained by the therapeutically achievable inhibition of DYRK1A and/or transporters. Arrhythmia cases are potentially linked to the structural features of hERG ion-channel inhibition found in rucaparib and niraparib. Enhanced psychiatric/nervous disorders associated with niraparib can be interpreted from the diverse neurotransporter off-targets reported. CONCLUSIONS: Despite their similar mode of action, the differential polypharmacology of PARP inhibitors influences their ADR profile.

AI delivers Michaelis constants as fuel for genome-scale metabolic models.

Michaelis constants (Km) are essential to predict the catalytic rate of enzymes, but are not widely available. A new study in PLOS Biology uses artificial intelligence (AI) to accurately predict Km on a proteome-wide scale, paving the way for dynamic, genome-wide modeling of metabolism.

Evolution of kinase polypharmacology across HSP90 drug discovery.

Most small molecules interact with several target proteins but this polypharmacology is seldom comprehensively investigated or explicitly exploited during drug discovery. Here, we use computational and experimental methods to identify and systematically characterize the kinase cross-pharmacology of representative HSP90 inhibitors. We demonstrate that the resorcinol clinical candidates ganetespib and, to a lesser extent, luminespib, display unique off-target kinase pharmacology as compared with other HSP90 inhibitors. We also demonstrate that polypharmacology evolved during the optimization to discover luminespib and that the hit, leads, and clinical candidate all have different polypharmacological profiles. We therefore recommend the computational and experimental characterization of polypharmacology earlier in drug discovery projects to unlock new multi-target drug design opportunities.

canSAR: update to the cancer translational research and drug discovery knowledgebase.

canSAR (http://cansar.icr.ac.uk) is the largest, public, freely available, integrative translational research and drug discovery knowledgebase for oncology. canSAR integrates vast multidisciplinary data from across genomic, protein, pharmacological, drug and chemical data with structural biology, protein networks and more. It also provides unique data, curation and annotation and crucially, AI-informed target assessment for drug discovery. canSAR is widely used internationally by academia and industry. Here we describe significant developments and enhancements to the data, web interface and infrastructure of canSAR in the form of the new implementation of the system: canSARblack. We demonstrate new functionality in aiding translation hypothesis generation and experimental design, and show how canSAR can be adapted and utilised outside oncology.

Public resources for chemical probes: the journey so far and the road ahead.

High-quality small molecule chemical probes are extremely valuable for biological research and target validation. However, frequent use of flawed small-molecule inhibitors produces misleading results and diminishes the robustness of biomedical research. Several public resources are available to facilitate assessment and selection of better chemical probes for specific protein targets. Here, we review chemical probe resources, discuss their current strengths and limitations, and make recommendations for further improvements. Expert review resources provide in-depth analysis but currently cover only a limited portion of the liganded proteome. Computational resources encompass more proteins and are regularly updated, but have limitations in data availability and curation. We show how biomedical scientists may use these resources to choose the best available chemical probes for their research.

The kinase polypharmacology landscape of clinical PARP inhibitors.

Polypharmacology plays an important role in defining response and adverse effects of drugs. For some mechanisms, experimentally mapping polypharmacology is commonplace, although this is typically done within the same protein class. Four PARP inhibitors have been approved by the FDA as cancer therapeutics, yet a precise mechanistic rationale to guide clinicians on which to choose for a particular patient is lacking. The four drugs have largely similar PARP family inhibition profiles, but several differences at the molecular and clinical level have been reported that remain poorly understood. Here, we report the first comprehensive characterization of the off-target kinase landscape of four FDA-approved PARP drugs. We demonstrate that all four PARP inhibitors have a unique polypharmacological profile across the kinome. Niraparib and rucaparib inhibit DYRK1s, CDK16 and PIM3 at clinically achievable, submicromolar concentrations. These kinases represent the most potently inhibited off-targets of PARP inhibitors identified to date and should be investigated further to clarify their potential implications for efficacy and safety in the clinic. Moreover, broad kinome profiling is recommended for the development of PARP inhibitors as PARP-kinase polypharmacology could potentially be exploited to modulate efficacy and side-effect profiles.

canSAR: update to the cancer translational research and drug discovery knowledgebase.

canSAR (http://cansar.icr.ac.uk) is a public, freely available, integrative translational research and drug discovery knowlegebase. canSAR informs researchers to help solve key bottlenecks in cancer translation and drug discovery. It integrates genomic, protein, pharmacological, drug and chemical data with structural biology, protein networks and unique, comprehensive and orthogonal 'druggability' assessments. canSAR is widely used internationally by academia and industry. Here we describe major enhancements to canSAR including new and expanded data. We also describe the first components of canSARblack-an advanced, responsive, multi-device compatible redesign of canSAR with a question-led interface.