Target prediction utilising negative bioactivity data covering large chemical space

A general prediction model for compound-protein interactions based on deep learning

Background

The identification of compound-protein interactions (CPIs) is crucial for drug discovery and understanding mechanisms of action. Accurate CPI prediction can elucidate drug-target-disease interactions, aiding in the discovery of candidate compounds and effective synergistic drugs, particularly from traditional Chinese medicine (TCM). Existing in silico methods face challenges in prediction accuracy and generalization due to compound and target diversity and the lack of largescale interaction datasets and negative datasets for model learning.

Methods

To address these issues, we developed a computational model for CPI prediction by integrating the constructed large-scale bioactivity benchmark dataset with a deep learning (DL) algorithm. To verify the accuracy of our CPI model, we applied it to predict the targets of compounds in TCM. An herb pair of Astragalus membranaceus and Hedyotis diffusaas was used as a model, and the active compounds in this herb pair were collected from various public databases and the literature. The complete targets of these active compounds were predicted by the CPI model, resulting in an expanded target dataset. This dataset was next used for the prediction of synergistic antitumor compound combinations. The predicted multi-compound combinations were subsequently examined through in vitro cellular experiments.

Results

Our CPI model demonstrated superior performance over other machine learning models, achieving an area under the Receiver Operating Characteristic curve (AUROC) of 0.98, an area under the precision-recall curve (AUPR) of 0.98, and an accuracy (ACC) of 93.31% on the test set. The model's generalization capability and applicability were further confirmed using external databases. Utilizing this model, we predicted the targets of compounds in the herb pair of Astragalus membranaceus and Hedyotis diffusaas, yielding an expanded target dataset. Then, we integrated this expanded target dataset to predict effective drug combinations using our drug synergy prediction model DeepMDS. Experimental assay on breast cancer cell line MDA-MB-231 proved the efficacy of the best predicted multi-compound combinations: Combination I (Epicatechin, Ursolic acid, Quercetin, Aesculetin and Astragaloside IV) exhibited a half-maximal inhibitory concentration (IC50) value of 19.41 μM, and a combination index (CI) value of 0.682; and Combination II (Epicatechin, Ursolic acid, Quercetin, Vanillic acid and Astragaloside IV) displayed a IC50 value of 23.83 μM and a CI value of 0.805. These results validated the ability of our model to make accurate predictions for novel CPI data outside the training dataset and evaluated the reliability of the predictions, showing good applicability potential in drug discovery and in the elucidation of the bioactive compounds in TCM.

Conclusion

Our CPI prediction model can serve as a useful tool for accurately identifying potential CPI for a wide range of proteins, and is expected to facilitate drug research, repurposing and support the understanding of TCM.

MolScore: a scoring, evaluation and benchmarking framework for generative models in de novo drug design

Generative models are undergoing rapid research and application to de novo drug design. To facilitate their application and evaluation, we present MolScore. MolScore already contains many drug-design-relevant scoring functions commonly used in benchmarks such as, molecular similarity, molecular docking, predictive models, synthesizability, and more. In addition, providing performance metrics to evaluate generative model performance based on the chemistry generated. With this unification of functionality, MolScore re-implements commonly used benchmarks in the field (such as GuacaMol, MOSES, and MolOpt). Moreover, new benchmarks can be created trivially. We demonstrate this by testing a chemical language model with reinforcement learning on three new tasks of increasing complexity related to the design of 5-HT2a ligands that utilise either molecular descriptors, 266 pre-trained QSAR models, or dual molecular docking. Lastly, MolScore can be integrated into an existing Python script with just three lines of code. This framework is a step towards unifying generative model application and evaluation as applied to drug design for both practitioners and researchers. The framework can be found on GitHub and downloaded directly from the Python Package Index.Scientific ContributionMolScore is an open-source platform to facilitate generative molecular design and evaluation thereof for application in drug design. This platform takes important steps towards unifying existing benchmarks, providing a platform to share new benchmarks, and improves customisation, flexibility and usability for practitioners over existing solutions.

MAVEN: compound mechanism of action analysis and visualisation using transcriptomics and compound structure data in R/Shiny

BACKGROUND

Understanding the Mechanism of Action (MoA) of a compound is an often challenging but equally crucial aspect of drug discovery that can help improve both its efficacy and safety. Computational methods to aid MoA elucidation usually either aim to predict direct drug targets, or attempt to understand modulated downstream pathways or signalling proteins. Such methods usually require extensive coding experience and results are often optimised for further computational processing, making them difficult for wet-lab scientists to perform, interpret and draw hypotheses from.

RESULTS

To address this issue, we in this work present MAVEN (Mechanism of Action Visualisation and Enrichment), an R/Shiny app which allows for GUI-based prediction of drug targets based on chemical structure, combined with causal reasoning based on causal protein-protein interactions and transcriptomic perturbation signatures. The app computes a systems-level view of the mechanism of action of the input compound. This is visualised as a sub-network linking predicted or known targets to modulated transcription factors via inferred signalling proteins. The tool includes a selection of MSigDB gene set collections to perform pathway enrichment on the resulting network, and also allows for custom gene sets to be uploaded by the researcher. MAVEN is hence a user-friendly, flexible tool for researchers without extensive bioinformatics or cheminformatics knowledge to generate interpretable hypotheses of compound Mechanism of Action.

CONCLUSIONS

MAVEN is available as a fully open-source tool at https://github.com/laylagerami/MAVEN with options to install in a Docker or Singularity container. Full documentation, including a tutorial on example data, is available at https://laylagerami.github.io/MAVEN .

The Coming of Age of AI/ML in Drug Discovery, Development, Clinical Testing, and Manufacturing: The FDA Perspectives

Artificial intelligence (AI) and machine learning (ML) represent significant advancements in computing, building on technologies that humanity has developed over millions of years-from the abacus to quantum computers. These tools have reached a pivotal moment in their development. In 2021 alone, the U.S. Food and Drug Administration (FDA) received over 100 product registration submissions that heavily relied on AI/ML for applications such as monitoring and improving human performance in compiling dossiers. To ensure the safe and effective use of AI/ML in drug discovery and manufacturing, the FDA and numerous other U.S. federal agencies have issued continuously updated, stringent guidelines. Intriguingly, these guidelines are often generated or updated with the aid of AI/ML tools themselves. The overarching goal is to expedite drug discovery, enhance the safety profiles of existing drugs, introduce novel treatment modalities, and improve manufacturing compliance and robustness. Recent FDA publications offer an encouraging outlook on the potential of these tools, emphasizing the need for their careful deployment. This has expanded market opportunities for retraining personnel handling these technologies and enabled innovative applications in emerging therapies such as gene editing, CRISPR-Cas9, CAR-T cells, mRNA-based treatments, and personalized medicine. In summary, the maturation of AI/ML technologies is a testament to human ingenuity. Far from being autonomous entities, these are tools created by and for humans designed to solve complex problems now and in the future. This paper aims to present the status of these technologies, along with examples of their present and future applications.

Design and Synthesis of New Anthranyl Phenylhydrazides: Antileishmanial Activity and Structure-Activity Relationship

Leishmaniasis is a neglected tropical disease affecting millions of people worldwide. A centenary approach to antimonial-based drugs was first initiated with the synthesis of urea stibamine by Upendranath Brahmachari in 1922. The need for new drug development led to resistance toward antimoniates. New drug development to treat leishmaniasis is urgently needed. In this way, searching for new substances with antileishmanial activity, we synthesized ten anthranyl phenylhydrazide and three quinazolinone derivatives and evaluated them against promastigotes and the intracellular amastigotes of Leishmania amazonensis. Three compounds showed good activity against promastigotes 1b, 1d, and 1g, with IC50 between 1 and 5 μM. These new phenylhydrazides were tested against Leishmania arginase, but they all failed to inhibit this parasite enzyme, as we have shown in a previous study. To explain the possible mechanism of action, we proposed the enzyme PTR1 as a new target for these compounds based on in silico analysis. In conclusion, the new anthranyl hydrazide derivatives can be a promising scaffold for developing new substances against the protozoa parasite.

ProfhEX: AI-based platform for small molecules liability profiling

Off-target drug interactions are a major reason for candidate failure in the drug discovery process. Anticipating potential drug's adverse effects in the early stages is necessary to minimize health risks to patients, animal testing, and economical costs. With the constantly increasing size of virtual screening libraries, AI-driven methods can be exploited as first-tier screening tools to provide liability estimation for drug candidates. In this work we present ProfhEX, an AI-driven suite of 46 OECD-compliant machine learning models that can profile small molecules on 7 relevant liability groups: cardiovascular, central nervous system, gastrointestinal, endocrine, renal, pulmonary and immune system toxicities. Experimental affinity data was collected from public and commercial data sources. The entire chemical space comprised 289'202 activity data for a total of 210'116 unique compounds, spanning over 46 targets with dataset sizes ranging from 819 to 18896. Gradient boosting and random forest algorithms were initially employed and ensembled for the selection of a champion model. Models were validated according to the OECD principles, including robust internal (cross validation, bootstrap, y-scrambling) and external validation. Champion models achieved an average Pearson correlation coefficient of 0.84 (SD of 0.05), an R² determination coefficient of 0.68 (SD = 0.1) and a root mean squared error of 0.69 (SD of 0.08). All liability groups showed good hit-detection power with an average enrichment factor at 5% of 13.1 (SD of 4.5) and AUC of 0.92 (SD of 0.05). Benchmarking against already existing tools demonstrated the predictive power of ProfhEX models for large-scale liability profiling. This platform will be further expanded with the inclusion of new targets and through complementary modelling approaches, such as structure and pharmacophore-based models. ProfhEX is freely accessible at the following address: https://profhex.exscalate.eu/ .

CXCR4 expression is elevated in TNBC patient derived samples and Z-guggulsterone abrogates tumor progression by targeting CXCL12/CXCR4 signaling axis in mice model

Environmental factors such as exposure to ionizing radiations, certain environmental pollutants, and toxic chemicals are considered as risk factors in the development of breast cancer. Triple-negative breast cancer (TNBC) is a molecular variant of breast cancer that lacks therapeutic targets such as progesterone receptor, estrogen receptor, and human epidermal growth factor receptor-2 which makes the targeted therapy ineffective in TNBC patients. Therefore, identification of new therapeutic targets for the treatment of TNBC and the discovery of new therapeutic agents is the need of the hour. In this study, CXCR4 was found to be highly expressed in majority of breast cancer tissues and metastatic lymph nodes derived from TNBC patients. CXCR4 expression is positively correlated with breast cancer metastasis and poor prognosis of TNBC patients suggesting that suppression of CXCR4 expression could be a good strategy in the treatment of TNBC patients. Therefore, the effect of Z-guggulsterone (ZGA) on the expression of CXCR4 in TNBC cells was examined. ZGA downregulated protein and mRNA expression of CXCR4 in TNBC cells and proteasome inhibition or lysosomal stabilization had no effect on the ZGA-induced CXCR4 reduction. CXCR4 is under the transcriptional control of NF-κB, whereas ZGA was found to downregulate transcriptional activity NF-κB. Functionally, ZGA downmodulated the CXCL12-driven migration/invasion in TNBC cells. Additionally, the effect of ZGA on growth of tumor was investigated in the orthotopic TNBC mice model and ZGA presented good inhibition of tumor growth and liver/lung metastasis in this model. Western blotting and immunohistochemical analysis indicated a reduction of CXCR4, NF-κB, and Ki67 in tumor tissues. Computational analysis suggested PXR agonism and FXR antagonism as targets of ZGA. In conclusion, CXCR4 was found to be overexpressed in majority of patient-derived TNBC tissues and ZGA abrogated the growth of TNBC tumors by partly targeting the CXCL12/CXCR4 signaling axis.

Mechanism of action deconvolution of the small-molecule pathological tau aggregation inhibitor Anle138b

BACKGROUND

A key histopathological hallmark of Alzheimer's disease (AD) is the presence of neurofibrillary tangles of aggregated microtubule-associated protein tau in neurons. Anle138b is a small molecule which has previously shown efficacy in mice in reducing tau aggregates and rescuing AD disease phenotypes.

METHODS

In this work, we employed bioinformatics analysis-including pathway enrichment and causal reasoning-of an in vitro tauopathy model. The model consisted of cultured rat cortical neurons either unseeded or seeded with tau aggregates derived from human AD patients, both of which were treated with Anle138b to generate hypotheses for its mode of action. In parallel, we used a collection of human target prediction models to predict direct targets of Anle138b based on its chemical structure.

RESULTS

Combining the different approaches, we found evidence supporting the hypothesis that the action of Anle138b involves several processes which are key to AD progression, including cholesterol homeostasis and neuroinflammation. On the pathway level, we found significantly enriched pathways related to these two processes including those entitled "Superpathway of cholesterol biosynthesis" and "Granulocyte adhesion and diapedesis". With causal reasoning, we inferred differential activity of SREBF1/2 (involved in cholesterol regulation) and mediators of the inflammatory response such as NFKB1 and RELA. Notably, our findings were also observed in Anle138b-treated unseeded neurons, meaning that the inferred processes are independent of tau pathology and thus represent the direct action of the compound in the cellular system. Through structure-based ligand-target prediction, we predicted the intracellular cholesterol carrier NPC1 as well as NF-κB subunits as potential targets of Anle138b, with structurally similar compounds in the model training set known to target the same proteins.

CONCLUSIONS

This study has generated feasible hypotheses for the potential mechanism of action of Anle138b, which will enable the development of future molecular interventions aiming to reduce tau pathology in AD patients.

Computer-aided Evaluation of Polyvalent Medications' Pharmacological Potential. Multiphytoadaptogen as a Case Study

Many human diseases including cancer, degenerative and autoimmune disorders, diabetes and others are multifactorial. Pharmaceutical agents acting on a single target do not provide their efficient curation. Multitargeted drugs exhibiting pleiotropic pharmacological effects have certain advantages due to the normalization of the complex pathological processes of different etiology. Extracts of medicinal plants (EMP) containing multiple phytocomponents are widely used in traditional medicines for multifactorial disorders' treatment. Experimental studies of pharmacological potential for multicomponent compositions are quite expensive and time-consuming. In silico evaluation of EMP the pharmacological potential may provide the basis for selecting the most promising directions of testing and for identifying potential additive/synergistic effects. Multiphytoadaptogen (MPhA) containing 70 major phytocomponents of different chemical classes from 40 medicinal plant extracts has been studied in vitro, in vivo and in clinical researches. Antiproliferative and anti-tumor activities have been shown against some tumors as well as evidence-based therapeutic effects against age-related pathologies. In addition, the neuroprotective, antioxidant, antimutagenic, radioprotective, and immunomodulatory effects of MPhA were confirmed. Analysis of the PASS profiles of the biological activity of MPhA phytocomponents showed that most of the predicted anti-tumor and anti-metastatic effects were consistent with the results of laboratory and clinical studies. Antimutagenic, immunomodulatory, radioprotective, neuroprotective and anti-Parkinsonian effects were also predicted for most of the phytocomponents. Effects associated with positive effects on the male and female reproductive systems have been identified too. Thus, PASS and PharmaExpert can be used to evaluate the pharmacological potential of complex pharmaceutical compositions containing natural products.

Systematic Analysis and Prediction of the Target Space of Bioactive Food Compounds: Filling the Chemobiological Gaps

Food compounds and their molecular interactions are crucial for health and provide new chemotypes and targets for drug and nutraceutic design. Here, we retrieve and analyze the complete set of published interactions of food compounds with human proteins using the FooDB as a compound set and ChEMBL as a source of interactions. The data are analyzed in terms of 19 target classes and 19 compound classes, showing a small fraction of target assignment for the compounds (1.6%) and unraveling multiple gaps in the chemobiological space for these molecules. By using well-established cheminformatic approaches [similarity ensemble approach (SEA) combined with the maximum Tanimoto coefficient to the nearest bioactive, "SEA + TC"], we achieve a much enhanced target assignment (64.2%), filling many of the gaps with target hypothesis for fast focused testing. By publishing these data sets and analyses, we expect to provide a set of resources to speed up the full clarification of the chemobiological space of food compounds, opening new opportunities for drug and nutraceutic design.

On modeling and utilizing chemical compound information with deep learning technologies: A task-oriented approach

A large number of chemical compounds are available in databases such as PubChem and ZINC. However, currently known compounds, though large, represent only a fraction of possible compounds, which is known as chemical space. Many of these compounds in the databases are annotated with properties and assay data that can be used for drug discovery efforts. For this goal, a number of machine learning algorithms have been developed and recent deep learning technologies can be effectively used to navigate chemical space, especially for unknown chemical compounds, in terms of drug-related tasks. In this article, we survey how deep learning technologies can model and utilize chemical compound information in a task-oriented way by exploiting annotated properties and assay data in the chemical compounds databases. We first compile what kind of tasks are trying to be accomplished by machine learning methods. Then, we survey deep learning technologies to show their modeling power and current applications for accomplishing drug related tasks. Next, we survey deep learning techniques to address the insufficiency issue of annotated data for more effective navigation of chemical space. Chemical compound information alone may not be powerful enough for drug related tasks, thus we survey what kind of information, such as assay and gene expression data, can be used to improve the prediction power of deep learning models. Finally, we conclude this survey with four important newly developed technologies that are yet to be fully incorporated into computational analysis of chemical information.

Harnessing Protein-Ligand Interaction Fingerprints to Predict New Scaffolds of RIPK1 Inhibitors

Necroptosis has emerged as an exciting target in oncological, inflammatory, neurodegenerative, and autoimmune diseases, in addition to acute ischemic injuries. It is known to play a role in innate immune response, as well as in antiviral cellular response. Here we devised a concerted in silico and experimental framework to identify novel RIPK1 inhibitors, a key necroptosis factor. We propose the first in silico model for the prediction of new RIPK1 inhibitor scaffolds by combining docking and machine learning methodologies. Through the data analysis of patterns in docking results, we derived two rules, where rule #1 consisted of a four-residue signature filter, and rule #2 consisted of a six-residue similarity filter based on docking calculations. These were used in consensus with a machine learning QSAR model from data collated from ChEMBL, the literature, in patents, and from PubChem data. The models allowed for good prediction of actives of >90, 92, and 96.4% precision, respectively. As a proof-of-concept, we selected 50 compounds from the ChemBridge database, using a consensus of both molecular docking and machine learning methods, and tested them in a phenotypic necroptosis assay and a biochemical RIPK1 inhibition assay. A total of 7 of the 47 tested compounds demonstrated around 20−25% inhibition of RIPK1’s kinase activity but, more importantly, these compounds were discovered to occupy new areas of chemical space. Although no strong actives were found, they could be candidates for further optimization, particularly because they have new scaffolds. In conclusion, this screening method may prove valuable for future screening efforts as it allows for the exploration of new areas of the chemical space in a very fast and inexpensive manner, therefore providing efficient starting points amenable to further hit-optimization campaigns.

Predicting clinical trial outcomes using drug bioactivities through graph database integration and machine learning

The ability to estimate the probability of a drug to receive approval in clinical trials provides natural advantages to optimizing pharmaceutical research workflows. Success rates of clinical trials have deep implications for costs, duration of development, and under pressure due to stringent regulatory approval processes. We propose a machine learning approach that can predict the outcome of the trial with reliable accuracies, using biological activities, physicochemical properties of the compounds, target-related features, and NLP-based compound representation. In the above list, biological activities have never been used as an independent variable towards the prediction of clinical trial outcomes. We have extracted the drug-disease pair from clinical trials and mapped target(s) to that pair using multiple data sources. Empirical results demonstrate that ensemble learning outperforms independently trained, small-data ML models. We report results and inferences derived from a Random forest classifier with an average accuracy of 93%, and an F1 score of 0.96 for the "Pass" class. "Pass" refers to one of the two classes (Pass/Fail) of all clinical trials, and the model performed well in predicting the "Pass" category. Through the analysis of feature contributions to predictive capability, we have demonstrated that bioactivity plays a statistically significant role in predicting clinical trial outcome. A significant effort has gone into the production of the dataset that, for the first time, integrates clinical trial information with protein targets. Cleaned, organized, integrated data and code to map these entities, created as a part of this work, are available open-source. This reproducibility and the freely available code ensure that researchers with access to deep curated and proprietary clinical trial databases (we only use open-source data in this study) can further expand the scope of the results.

Yes SIR! On the structure-inactivity relationships in drug discovery

In analogy with structure-activity relationships (SARs), which are at the core of medicinal chemistry, studying structure-inactivity relationships (SIRs) is essential to understanding and predicting biological activity. Current computational methods should predict or distinguish 'activity' and 'inactivity' with the same confidence because both concepts are complementary. However, the lack of inactivity data, in particular in the public domain, limits the development of predictive models and its broad application. In this review, we encourage the scientific community to disclose and analyze high-confidence activity data considering both the labeled 'active' and 'inactive' compounds.

Prediction of In Vivo Pharmacokinetic Parameters and Time-Exposure Curves in Rats Using Machine Learning from the Chemical Structure

Animal pharmacokinetic (PK) data as well as human and animal in vitro systems are utilized in drug discovery to define the rate and route of drug elimination. Accurate prediction and mechanistic understanding of drug clearance and disposition in animals provide a degree of confidence for extrapolation to humans. In addition, prediction of in vivo properties can be used to improve design during drug discovery, help select compounds with better properties, and reduce the number of in vivo experiments. In this study, we generated machine learning models able to predict rat in vivo PK parameters and concentration-time PK profiles based on the molecular chemical structure and either measured or predicted in vitro parameters. The models were trained on internal in vivo rat PK data for over 3000 diverse compounds from multiple projects and therapeutic areas, and the predicted endpoints include clearance and oral bioavailability. We compared the performance of various traditional machine learning algorithms and deep learning approaches, including graph convolutional neural networks. The best models for PK parameters achieved R² = 0.63 [root mean squared error (RMSE) = 0.26] for clearance and R² = 0.55 (RMSE = 0.46) for bioavailability. The models provide a fast and cost-efficient way to guide the design of molecules with optimal PK profiles, to enable the prediction of virtual compounds at the point of design, and to drive prioritization of compounds for in vivo assays.

Machine learning-based predictive models for identifying high active compounds against HIV-1 integrase

HIV-integrase is an important drug target because it catalyzes chromosomal integration of proviral DNA towards establishing latent infection. Computer-aided drug design has immensely contributed to identifying and developing novel antiviral drugs. We have developed various machine learning-based predictive models for identifying high activity compounds against HIV-integrase. Multiclass models were built using support vector machine with reasonable accuracy on the test and evaluation sets. The developed models were evaluated by rigorous validation approaches and the best features were selected by Boruta method. As compared to the model developed from all descriptors set, a slight improvement was observed among the selected descriptors. Validated models were further used for virtual screening of potential compounds from ChemBridge library. Of the six high active compounds predicted from selected models, compounds 9103124, 6642917 and 9082952 showed the most reasonable binding-affinity and stable-interaction with HIV-integrase active-site residues Asp64, Glu152 and Asn155. This was in agreement with previous reports on the essentiality of these residues against a wide range of inhibitors. We therefore highlight the rigorosity of validated classification models for accurate prediction and ranking of high active lead drugs against HIV-integrase.

Computational analyses of mechanism of action (MoA): data, methods and integration

The elucidation of a compound's Mechanism of Action (MoA) is a challenging task in the drug discovery process, but it is important in order to rationalise phenotypic findings and to anticipate potential side-effects. Bioinformatic approaches, advances in machine learning techniques and the increasing deposition of high-throughput data in public databases have significantly contributed to recent advances in the field, but it is not straightforward to decide which data and methods are most suitable to use in a given case. In this review, we focus on these methods and data and their applications in generating MoA hypotheses for subsequent experimental validation. We discuss compound-specific data such as -omics, cell morphology and bioactivity data, as well as commonly used supplementary prior knowledge such as network and pathway data, and provide information on databases where this data can be accessed. In terms of methodologies, we discuss both well-established methods (connectivity mapping, pathway enrichment) as well as more developing methods (neural networks and multi-omics integration). Finally, we review case studies where the MoA of a compound was successfully suggested from computational analysis by incorporating multiple data modalities and/or methodologies. Our aim for this review is to provide researchers with insights into the benefits and drawbacks of both the data and methods in terms of level of understanding, biases and interpretation - and to highlight future avenues of investigation which we foresee will improve the field of MoA elucidation, including greater public access to -omics data and methodologies which are capable of data integration.

"Efficacy-Nature-Structure" Relationship of Traditional Chinese Medicine Based on Chemical Structural Data and Bioinformatics Analysis

Traditional Chinese medicines (TCMs) have wide pharmacological activities, and the ingredients in individual TCMs determine their efficacies. To understand the "efficacy-nature-structure" relationship of TCM, compounds from 2444 kinds of herbs were collected, and the associations between family, structure, nature, and biological activities were mined and analyzed. Bernoulli Naïve Bayes profiling and a data analysis method were used to predict the targets of compounds. The results show that genetic material determined the representation of ingredients from herbs and the nature of TCMs and that the superior scaffolds of compounds of cold nature were 2-phenylochrotinone, anthraquinone, and coumarin, while the compounds of hot nature were cyclohexene. The results of the similarity analysis and distribution for molecular descriptors of compounds show that compounds associated with the same nature were similar and compounds associated with different natures occurred as a transition in part. As for integral compounds from 2-phenylochrotinone, anthraquinone, coumarin, and cyclohexene, the value of the shape index increased, indicating the transition of scaffolds from a spherical structure to a linear structure, with various molecular descriptors decreasing. Three medicines and three recipes prescribed based on "efficacy-nature-structure" had a higher survival rate in the clinic and provided powerful evidence for TCM principles. The research improves the understanding of the "efficacy-nature-structure" relationship and extends TCM applications.

(R)-(+)-Lasiodiplodin isolated from the endophytic fungus Sordaria tamaensis exhibits potent antimycobacterial and anti-inflammatory activities in vitro and in vivo: a dual approach for the treatment of severe pulmonary tuberculosis

OBJECTIVES

This study aimed to evaluate endophytic fungi isolated from Tocoyena bullata and Humiria balsamifera plant species for their antimycobacterial and anti-inflammatory activities, focusing on severe pulmonary tuberculosis cases which are often associated with exacerbated inflammation.

METHODS

Mycobacterium suspensions were incubated with the samples for 5 days. RAW 264.7 macrophages stimulated with LPS were also incubated with them for 24 h to assess the inhibition of inflammatory mediator production and cytotoxicity. C57BL/6 mice were infected with Mtb M299 and treated for 15 days with lasiodiplodin (Lasio).

KEY FINDINGS

Endophytic fungus Sordaria tamaensis, obtained from T. bullata, was the most promising. Its ethanolic extract impaired mycobacterial growth with MIC50 (µg/ml): 1.5 ± 0.6 (BCG), 66.8 ± 0.1 (H37Rv) and 80.0 ± 0.1 (M299). (R)-(+)-Lasio showed MIC50 92.2 ± 1.8 µg/ml (M299). In addition, Lasio was able to inhibit NO, IL-1β and TNF-α production and was not cytotoxic for macrophages. M. tuberculosis-infected C57BL/6 animals treated by Lasio reduced the number of acid-fast bacilli, lung pathology, leucocyte influx and proinflammatory cytokine production in the lungs. The class IIa fructose 1,6-bisphosphate aldolase was the predicted hypothetical target of Lasio.

CONCLUSIONS

(R)-(+)-Lasio stood out as a promising anti-TB compound, exhibiting anti-inflammatory and antimycobacterial effects, as well as low cytotoxicity.

Natural products from Brazilian biodiversity identified as potential inhibitors of PknA and PknB of M. tuberculosis using molecular modeling tools

Mycobacterium tuberculosis was discovered in 1882 by Robert Koch but, since its discovery, the tuberculosis (TB) epidemic has endured, being one of the top 10 causes of death worldwide. Drug-resistant TB continues to be a public health threat and bioactive compounds with a new mode of action (MoA) are needed to overcome this. Since natural products are described as important sources for the development of new drugs, the objective of this work was to identify potential ligands from Brazilian natural products (NPs) for M. tuberculosis targets using molecular modeling tools. Using chemogenomics we identified the Serine/Threonine Protein Kinase PknB as a putative target for 13 NPs from a database from Brazilian biodiversity (NuBBE). Literature data supported further investigation of NuBBE105, NuBBE598, NuBBE936, NuBBE964, NuBBE1045, and NuBBE1180 by molecular docking and dynamics. Key interactions were observed with PknB and simulations confirmed stability and favorable binding energies. Considering structural similarity with PknB, we further explored binding of the NPs to PknA, critical for M. tuberculosis survival, and all of them resembled important interactions with the enzyme, showing stable and favorable binding energies, whilst van der Waals interactions seem to play a key role for binding to PknA and PknB. NuBBE936 and NuBBE1180 have already had their antimycobacterial activity reported and our results can provide a basis for their MoA. Finally, the other NPs which have not been tested against M. tuberculosis deserve further investigation, aiming at the discovery of antimycobacterial drug candidates with innovative MoA.

Development of Machine Learning Models for Accurately Predicting and Ranking the Activity of Lead Molecules to Inhibit PRC2 Dependent Cancer

Disruption of epigenetic processes to eradicate tumor cells is among the most promising interventions for cancer control. EZH2 (Enhancer of zeste homolog 2), a catalytic component of polycomb repressive complex 2 (PRC2), methylates lysine 27 of histone H3 to promote transcriptional silencing and is an important drug target for controlling cancer via epigenetic processes. In the present study, we have developed various predictive models for modeling the inhibitory activity of EZH2. Binary and multiclass models were built using SVM, random forest and XGBoost methods. Rigorous validation approaches including predictiveness curve, Y-randomization and applicability domain (AD) were employed for evaluation of the developed models. Eighteen descriptors selected from Boruta methods have been used for modeling. For binary classification, random forest and XGBoost achieved an accuracy of 0.80 and 0.82, respectively, on external test set. Contrastingly, for multiclass models, random forest and XGBoost achieved an accuracy of 0.73 and 0.75, respectively. 500 Y-randomization runs demonstrate that the models were robust and the correlations were not by chance. Evaluation metrics from predictiveness curve show that the selected eighteen descriptors predict active compounds with total gain (TG) of 0.79 and 0.59 for XGBoost and random forest, respectively. Validated models were further used for virtual screening and molecular docking in search of potential hits. A total of 221 compounds were commonly predicted as active with above the set probability threshold and also under the AD of training set. Molecular docking revealed that three compounds have reasonable binding energy and favorable interactions with critical residues in the active site of EZH2. In conclusion, we highlighted the potential of rigorously validated models for accurately predicting and ranking the activities of lead molecules against cancer epigenetic targets. The models presented in this study represent the platform for development of EZH2 inhibitors.

Sulfated and Sulfur-Containing Steroids and Their Pharmacological Profile

The review focuses on sulfated steroids that have been isolated from seaweeds, marine sponges, soft corals, ascidians, starfish, and other marine invertebrates. Sulfur-containing steroids and triterpenoids are sourced from sedentary marine coelenterates, plants, marine sediments, crude oil, and other geological deposits. The review presents the pharmacological profile of sulfated steroids, sulfur-containing steroids, and triterpenoids, which is based on data obtained using the PASS program. In addition, several semi-synthetic and synthetic epithio steroids, which represent a rare group of bioactive lipids that have not yet been found in nature, but possess a high level of antitumor activity, were included in this review for the comparative pharmacological characterization of this class of compounds. About 140 steroids and triterpenoids are presented in this review, which demonstrate a wide range of biological activities. Therefore, out of 71 sulfated steroids, thirteen show strong antitumor activity with a confidence level of more than 90%, out of 50 sulfur-containing steroids, only four show strong antitumor activity with a confidence level of more than 93%, and out of eighteen epithio steroids, thirteen steroids show strong antitumor activity with a confidence level of 91% to 97.4%.

Comparison of Chemical Structure and Cell Morphology Information for Multitask Bioactivity Predictions

The understanding of the mechanism-of-action (MoA) of compounds and the prediction of potential drug targets play an important role in small-molecule drug discovery. The aim of this work was to compare chemical and cell morphology information for bioactivity prediction. The comparison was performed using bioactivity data from the ExCAPE database, image data (in the form of CellProfiler features) from the Cell Painting data set (the largest publicly available data set of cell images with ∼30,000 compound perturbations), and extended connectivity fingerprints (ECFPs) using the multitask Bayesian matrix factorization (BMF) approach Macau. We found that the BMF Macau and random forest (RF) performance were overall similar when ECFPs were used as compound descriptors. However, BMF Macau outperformed RF in 159 out of 224 targets (71%) when image data were used as compound information. Using BMF Macau, 100 (corresponding to about 45%) and 90 (about 40%) of the 224 targets were predicted with high predictive performance (AUC > 0.8) with ECFP data and image data as side information, respectively. There were targets better predicted by image data as side information, such as β-catenin, and others better predicted by fingerprint-based side information, such as proteins belonging to the G-protein-Coupled Receptor 1 family, which could be rationalized from the underlying data distributions in each descriptor domain. In conclusion, both cell morphology changes and chemical structure information contain information about compound bioactivity, which is also partially complementary, and can hence contribute to in silico MoA analysis.

Antiprotozoal and Antitumor Activity of Natural Polycyclic Endoperoxides: Origin, Structures and Biological Activity

Polycyclic endoperoxides are rare natural metabolites found and isolated in plants, fungi, and marine invertebrates. The purpose of this review is a comparative analysis of the pharmacological potential of these natural products. According to PASS (Prediction of Activity Spectra for Substances) estimates, they are more likely to exhibit antiprotozoal and antitumor properties. Some of them are now widely used in clinical medicine. All polycyclic endoperoxides presented in this article demonstrate antiprotozoal activity and can be divided into three groups. The third group includes endoperoxides, which show weak antiprotozoal activity with a reliability of up to 70%, and this group includes only 1.1% of metabolites. The second group includes the largest number of endoperoxides, which are 65% and show average antiprotozoal activity with a confidence level of 70 to 90%. Lastly, the third group includes endoperoxides, which are 33.9% and show strong antiprotozoal activity with a confidence level of 90 to 99.6%. Interestingly, artemisinin and its analogs show strong antiprotozoal activity with 79 to 99.6% confidence against obligate intracellular parasites which belong to the genera Plasmodium, Toxoplasma, Leishmania, and Coccidia. In addition to antiprotozoal activities, polycyclic endoperoxides show antitumor activity in the proportion: 4.6% show weak activity with a reliability of up to 70%, 65.6% show an average activity with a reliability of 70 to 90%, and 29.8% show strong activity with a reliability of 90 to 98.3%. It should also be noted that some polycyclic endoperoxides, in addition to antiprotozoal and antitumor properties, show other strong activities with a confidence level of 90 to 97%. These include antifungal activity against the genera Aspergillus, Candida, and Cryptococcus, as well as anti-inflammatory activity. This review provides insights on further utilization of polycyclic endoperoxides by medicinal chemists, pharmacologists, and the pharmaceutical industry.

Comparative study of Danshen and Siwu decoction based on the molecular structures of the components and predicted targets

BACKGROUND

The sentence of "Danshen (Salvia Miltiorrhizae Radix et Rhizoma) and Siwu decoction are similar in function" was first recorded in an ancient Chinese medical book "Fu Ren Ming Li Lun". This theory has far-reaching influence on the clinical practice of Chinese medicine and is highly respected by Chinese medical doctors. However, the theory has limitations and controversial part for there is no in-depth and system comparative study.

METHODS

We collected the molecular structures of 129 compounds of Danshen and 81 compounds of Siwu decoction from the literatures. MACCS fingerprints and Tanimoto similarity were calculated based on the molecular structures for comparing the structural feature. Molecular descriptors which represent physical and chemical properties were calculated by Discovery Studio. Principal component analysis (PCA) of was performed based on the descriptors. The ADMET properties were predicted by FAF-Drugs4. The effect targets for the compounds with good ADMET properties were confirmed from experimental data and predicted using the algorithm comprising Bernoulli Naive Bayes profiling.

RESULTS

Based on the molecular structures, the presented study compared the structural feature, physical and chemical properties, ADMET properties, and effect targets of compounds of Danshen and Siwu decoction. It is found that Danshen and Siwu decoction do not have the same main active components. Moreover, the 2D structure of compounds from Danshen and Siwu decoction is not similar. Some of the compounds of Danshen and Siwu decoction are similar in 3D structure. The compounds with good ADMET properties of Danshen and Siwu decoction have same predicted targets, but some have different targets.

CONCLUSIONS

It can be inferred from the result that Danshen and Siwu decoction have some similarities, but also present differences from each other in the structure of the compounds and predicted targets. This may be the material basis of the similar and different traditional efficacy of Danshen and Siwu decoction. The setence of " Danshen and Siwu decoction are similar in function. " which is used in clinical has its material basis and target connotation to some extent. However, the traditional effects of Danshen and Siwu decoction are not exactly the same.

Prediction and mechanistic analysis of drug-induced liver injury (DILI) based on chemical structure

BACKGROUND

Drug-induced liver injury (DILI) is a major safety concern characterized by a complex and diverse pathogenesis. In order to identify DILI early in drug development, a better understanding of the injury and models with better predictivity are urgently needed. One approach in this regard are in silico models which aim at predicting the risk of DILI based on the compound structure. However, these models do not yet show sufficient predictive performance or interpretability to be useful for decision making by themselves, the former partially stemming from the underlying problem of labeling the in vivo DILI risk of compounds in a meaningful way for generating machine learning models.

RESULTS

As part of the Critical Assessment of Massive Data Analysis (CAMDA) "CMap Drug Safety Challenge" 2019 ( http://camda2019.bioinf.jku.at ), chemical structure-based models were generated using the binarized DILIrank annotations. Support Vector Machine (SVM) and Random Forest (RF) classifiers showed comparable performance to previously published models with a mean balanced accuracy over models generated using 5-fold LOCO-CV inside a 10-fold training scheme of 0.759 ± 0.027 when predicting an external test set. In the models which used predicted protein targets as compound descriptors, we identified the most information-rich proteins which agreed with the mechanisms of action and toxicity of nonsteroidal anti-inflammatory drugs (NSAIDs), one of the most important drug classes causing DILI, stress response via TP53 and biotransformation. In addition, we identified multiple proteins involved in xenobiotic metabolism which could be novel DILI-related off-targets, such as CLK1 and DYRK2. Moreover, we derived potential structural alerts for DILI with high precision, including furan and hydrazine derivatives; however, all derived alerts were present in approved drugs and were over specific indicating the need to consider quantitative variables such as dose.

CONCLUSION

Using chemical structure-based descriptors such as structural fingerprints and predicted protein targets, DILI prediction models were built with a predictive performance comparable to previous literature. In addition, we derived insights on proteins and pathways statistically (and potentially causally) linked to DILI from these models and inferred new structural alerts related to this adverse endpoint.

Prediction of pharmacological activities from chemical structures with graph convolutional neural networks

Many therapeutic drugs are compounds that can be represented by simple chemical structures, which contain important determinants of affinity at the site of action. Recently, graph convolutional neural network (GCN) models have exhibited excellent results in classifying the activity of such compounds. For models that make quantitative predictions of activity, more complex information has been utilized, such as the three-dimensional structures of compounds and the amino acid sequences of their respective target proteins. As another approach, we hypothesized that if sufficient experimental data were available and there were enough nodes in hidden layers, a simple compound representation would quantitatively predict activity with satisfactory accuracy. In this study, we report that GCN models constructed solely from the two-dimensional structural information of compounds demonstrated a high degree of activity predictability against 127 diverse targets from the ChEMBL database. Using the information entropy as a metric, we also show that the structural diversity had less effect on the prediction performance. Finally, we report that virtual screening using the constructed model identified a new serotonin transporter inhibitor with activity comparable to that of a marketed drug in vitro and exhibited antidepressant effects in behavioural studies.

Artificial intelligence in drug discovery: what is realistic, what are illusions? Part 1: Ways to make an impact, and why we are not there yet

Although artificial intelligence (AI) has had a profound impact on areas such as image recognition, comparable advances in drug discovery are rare. This article quantifies the stages of drug discovery in which improvements in the time taken, success rate or affordability will have the most profound overall impact on bringing new drugs to market. Changes in clinical success rates will have the most profound impact on improving success in drug discovery; in other words, the quality of decisions regarding which compound to take forward (and how to conduct clinical trials) are more important than speed or cost. Although current advances in AI focus on how to make a given compound, the question of which compound to make, using clinical efficacy and safety-related end points, has received significantly less attention. As a consequence, current proxy measures and available data cannot fully utilize the potential of AI in drug discovery, in particular when it comes to drug efficacy and safety in vivo. Thus, addressing the questions of which data to generate and which end points to model will be key to improving clinically relevant decision-making in the future.

Adding Stochastic Negative Examples into Machine Learning Improves Molecular Bioactivity Prediction

Multitask deep neural networks learn to predict ligand-target binding by example, yet public pharmacological data sets are sparse, imbalanced, and approximate. We constructed two hold-out benchmarks to approximate temporal and drug-screening test scenarios, whose characteristics differ from a random split of conventional training data sets. We developed a pharmacological data set augmentation procedure, Stochastic Negative Addition (SNA), which randomly assigns untested molecule-target pairs as transient negative examples during training. Under the SNA procedure, drug-screening benchmark performance increases from R² = 0.1926 ± 0.0186 to 0.4269 ± 0.0272 (122%). This gain was accompanied by a modest decrease in the temporal benchmark (13%). SNA increases in drug-screening performance were consistent for classification and regression tasks and outperformed y-randomized controls. Our results highlight where data and feature uncertainty may be problematic and how leveraging uncertainty into training improves predictions of drug-target relationships.

Uncertainty quantification in drug design

Machine learning and artificial intelligence are increasingly being applied to the drug-design process as a result of the development of novel algorithms, growing access, the falling cost of computation and the development of novel technologies for generating chemically and biologically relevant data. There has been recent progress in fields such as molecular de novo generation, synthetic route prediction and, to some extent, property predictions. Despite this, most research in these fields has focused on improving the accuracy of the technologies, rather than on quantifying the uncertainty in the predictions. Uncertainty quantification will become a key component in autonomous decision making and will be crucial for integrating machine learning and chemistry automation to create an autonomous design-make-test-analyse cycle. This review covers the empirical, frequentist and Bayesian approaches to uncertainty quantification, and outlines how they can be used for drug design. We also outline the impact of uncertainty quantification on decision making.

Enhancing Chemogenomics with Predictive Pharmacology

One of the grand challenges in contemporary chemical biology is the generation of a probe for every member of the human proteome. Probe selection and optimization strategies typically rely on experimental bioactivity data to determine the potency and selectivity of candidate molecules. However, this approach is profoundly limited by the sparsity of the known data, the annotation bias often found in the literature, and the cost of physical screening. Recent advancements in predictive pharmacology, such as the application of multitask and transfer learning, as well as the use of biologically motivated, structure-agnostic features to characterize molecules, should serve to mitigate these issues. Computational modeling likely offers the only cost-effective approach to substantially increasing the bioactivity annotation density both on the local and global scale and thus, we argue, will need to make a substantial contribution if the ambitious goals of probing the human proteome are to be realized in the foreseeable future.

Computational Drug Repositioning: Current Progress and Challenges

Novel drug discovery is time-consuming, costly, and a high-investment process due to the high attrition rate. Therefore, many trials are conducted to reuse existing drugs to treat pressing conditions and diseases, since their safety profiles and pharmacokinetics are already available. Drug repositioning is a strategy to identify a new indication of existing or already approved drugs, beyond the scope of their original use. Various computational and experimental approaches to incorporate available resources have been suggested for gaining a better understanding of disease mechanisms and the identification of repurposed drug candidates for personalized pharmacotherapy. In this review, we introduce publicly available databases for drug repositioning and summarize the approaches taken for drug repositioning. We also highlight and compare their characteristics and challenges, which should be addressed for the future realization of drug repositioning.

DTiGEMS+: drug-target interaction prediction using graph embedding, graph mining, and similarity-based techniques

In silico prediction of drug–target interactions is a critical phase in the sustainable drug development process, especially when the research focus is to capitalize on the repositioning of existing drugs. However, developing such computational methods is not an easy task, but is much needed, as current methods that predict potential drug–target interactions suffer from high false-positive rates. Here we introduce DTiGEMS+, a computational method that predicts Drug–Target interactions using Graph Embedding, graph Mining, and Similarity-based techniques. DTiGEMS+ combines similarity-based as well as feature-based approaches, and models the identification of novel drug–target interactions as a link prediction problem in a heterogeneous network. DTiGEMS+ constructs the heterogeneous network by augmenting the known drug–target interactions graph with two other complementary graphs namely: drug–drug similarity, target–target similarity. DTiGEMS+ combines different computational techniques to provide the final drug target prediction, these techniques include graph embeddings, graph mining, and machine learning. DTiGEMS+ integrates multiple drug–drug similarities and target–target similarities into the final heterogeneous graph construction after applying a similarity selection procedure as well as a similarity fusion algorithm. Using four benchmark datasets, we show DTiGEMS+ substantially improves prediction performance compared to other state-of-the-art in silico methods developed to predict of drug-target interactions by achieving the highest average AUPR across all datasets (0.92), which reduces the error rate by 33.3% relative to the second-best performing model in the state-of-the-art methods comparison.

DNMT1: A key drug target in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer. Altered epigenetics regulation including DNA hypermethylation by DNA methyltransferase 1 (DNMT1) has been implicated as one of the causes of TNBC tumorigenesis. In this review, the oncogenic functions rendered by DNMT1 in TNBCs, and DNMT1 inhibitors targeting TNBC cells are presented and discussed. In summary, DNMT1 expression is associated with poor breast cancer survival, and it is overexpressed in TNBC subtype. The oncogenic roles of DNMT1 in TNBCs include: (1) Repression of estrogen receptor (ER) expression; (2) Promotion of epithelial-mesenchymal transition (EMT) required for metastasis; (3) Induces cellular autophagy and; (4) Promotes the growth of cancer stem cells in TNBCs. DNMT1 confers these phenotypes by hypermethylating the promoter regions of ER, multiple tumor suppressor genes, microRNAs and epithelial markers involved in suppressing EMT. DNMT1 inhibitors exert anti-tumorigenic effects against TNBC cells. This includes the hypomethylating agents azacitidine, decitabine and guadecitabine that might sensitize TNBC patients to immune checkpoint blockade therapy. DNMT1 represents an epigenetic target for TNBC cells destruction as well as to derail their metastatic and aggressive phenotypes.

Using Predicted Bioactivity Profiles to Improve Predictive Modeling

Predictive modeling is a cornerstone in early drug development. Using information for multiple domains or across prediction tasks has the potential to improve the performance of predictive modeling. However, aggregating data often leads to incomplete data matrices that might be limiting for modeling. In line with previous studies, we show that by generating predicted bioactivity profiles, and using these as additional features, prediction accuracy of biological endpoints can be improved. Using conformal prediction, a type of confidence predictor, we present a robust framework for the calculation of these profiles and the evaluation of their impact. We report on the outcomes from several approaches to generate the predicted profiles on 16 datasets in cytotoxicity and bioactivity and show that efficiency is improved the most when including the p-values from conformal prediction as bioactivity profiles.

Industry-scale application and evaluation of deep learning for drug target prediction

Artificial intelligence (AI) is undergoing a revolution thanks to the breakthroughs of machine learning algorithms in computer vision, speech recognition, natural language processing and generative modelling. Recent works on publicly available pharmaceutical data showed that AI methods are highly promising for Drug Target prediction. However, the quality of public data might be different than that of industry data due to different labs reporting measurements, different measurement techniques, fewer samples and less diverse and specialized assays. As part of a European funded project (ExCAPE), that brought together expertise from pharmaceutical industry, machine learning, and high-performance computing, we investigated how well machine learning models obtained from public data can be transferred to internal pharmaceutical industry data. Our results show that machine learning models trained on public data can indeed maintain their predictive power to a large degree when applied to industry data. Moreover, we observed that deep learning derived machine learning models outperformed comparable models, which were trained by other machine learning algorithms, when applied to internal pharmaceutical company datasets. To our knowledge, this is the first large-scale study evaluating the potential of machine learning and especially deep learning directly at the level of industry-scale settings and moreover investigating the transferability of publicly learned target prediction models towards industrial bioactivity prediction pipelines.

Elucidating Compound Mechanism of Action and Predicting Cytotoxicity Using Machine Learning Approaches, Taking Prediction Confidence into Account

The modes of action (MoAs) of drugs frequently are unknown, because many are small molecules initially identified from phenotypic screens, giving rise to the need to elucidate their MoAs. In addition, the high attrition rate for candidate drugs in preclinical studies due to intolerable toxicity has motivated the development of computational approaches to predict drug candidate (cyto)toxicity as early as possible in the drug-discovery process. Here, we provide detailed instructions for capitalizing on bioactivity predictions to elucidate the MoAs of small molecules and infer their underlying phenotypic effects. We illustrate how these predictions can be used to infer the underlying antidepressive effects of marketed drugs. We also provide the necessary functionalities to model cytotoxicity data using single and ensemble machine-learning algorithms. Finally, we give detailed instructions on how to calculate confidence intervals for individual predictions using the conformal prediction framework. © 2019 by John Wiley & Sons, Inc.

Applications of Deep-Learning in Exploiting Large-Scale and Heterogeneous Compound Data in Industrial Pharmaceutical Research

In recent years, the development of high-throughput screening (HTS) technologies and their establishment in an industrialized environment have given scientists the possibility to test millions of molecules and profile them against a multitude of biological targets in a short period of time, generating data in a much faster pace and with a higher quality than before. Besides the structure activity data from traditional bioassays, more complex assays such as transcriptomics profiling or imaging have also been established as routine profiling experiments thanks to the advancement of Next Generation Sequencing or automated microscopy technologies. In industrial pharmaceutical research, these technologies are typically established in conjunction with automated platforms in order to enable efficient handling of screening collections of thousands to millions of compounds. To exploit the ever-growing amount of data that are generated by these approaches, computational techniques are constantly evolving. In this regard, artificial intelligence technologies such as deep learning and machine learning methods play a key role in cheminformatics and bio-image analytics fields to address activity prediction, scaffold hopping, de novo molecule design, reaction/retrosynthesis predictions, or high content screening analysis. Herein we summarize the current state of analyzing large-scale compound data in industrial pharmaceutical research and describe the impact it has had on the drug discovery process over the last two decades, with a specific focus on deep-learning technologies.

An Integrated System Biology Approach Yields Drug Repositioning Candidates for the Treatment of Heart Failure

Identifying effective pharmacological treatments for heart failure (HF) patients remains critically important. Given that the development of drugs de novo is expensive and time consuming, drug repositioning has become an increasingly important branch. In the present study, we propose a two-step drug repositioning pipeline and investigate the novel therapeutic effects of existing drugs approved by the US Food and Drug Administration to discover potential therapeutic drugs for HF. In the first step, we compared the gene expression pattern of HF patients with drug-induced gene expression profiles to obtain preliminary candidates. In the second step, we performed a systems biology approach based on the known protein–protein interaction information and targets of drugs to narrow down preliminary candidates to obtain final candidates. Drug set enrichment analysis and literature search were applied to assess the performance of our repositioning approach. We also constructed a mode of action network for each candidate and performed pathway analysis for each candidate using genes contained in their mode of action network to uncover pathways that potentially reflect the mechanisms of candidates’ therapeutic efficacy to HF. We discovered numerous preliminary candidates, some of which are used in clinical practice and supported by the literature. The final candidates contained nearly all of the preliminary candidates supported by previous studies. Drug set enrichment analysis and literature search support the validity of our repositioning approach. The mode of action network for each candidate not only displayed the underlying mechanisms of drug efficacy but also uncovered potential biomarkers and therapeutic targets for HF. Our two-step drug repositioning approach is efficient to find candidates with potential therapeutic efficiency to HF supported by the literature and might be of particular use in the discovery of novel effective pharmacological therapies for HF.

Prediction Model with High-Performance Constitutive Androstane Receptor (CAR) Using DeepSnap-Deep Learning Approach from the Tox21 10K Compound Library

The constitutive androstane receptor (CAR) plays pivotal roles in drug-induced liver injury through the transcriptional regulation of drug-metabolizing enzymes and transporters. Thus, identifying regulatory factors for CAR activation is important for understanding its mechanisms. Numerous studies conducted previously on CAR activation and its toxicity focused on in vivo or in vitro analyses, which are expensive, time consuming, and require many animals. We developed a computational model that predicts agonists for the CAR using the Toxicology in the 21st Century 10k library. Additionally, we evaluate the prediction performance of novel deep learning (DL)-based quantitative structure-activity relationship analysis called the DeepSnap-DL approach, which is a procedure of generating an omnidirectional snapshot portraying three-dimensional (3D) structures of chemical compounds. The CAR prediction model, which applies a 3D structure generator tool, called CORINA-generated and -optimized chemical structures, in the DeepSnap-DL demonstrated better performance than the existing methods using molecular descriptors. These results indicate that high performance in the prediction model using the DeepSnap-DL approach may be important to prepare suitable 3D chemical structures as input data and to enable the identification of modulators of the CAR.

Recent applications of deep learning and machine intelligence on in silico drug discovery: methods, tools and databases

The identification of interactions between drugs/compounds and their targets is crucial for the development of new drugs. In vitro screening experiments (i.e. bioassays) are frequently used for this purpose; however, experimental approaches are insufficient to explore novel drug-target interactions, mainly because of feasibility problems, as they are labour intensive, costly and time consuming. A computational field known as 'virtual screening' (VS) has emerged in the past decades to aid experimental drug discovery studies by statistically estimating unknown bio-interactions between compounds and biological targets. These methods use the physico-chemical and structural properties of compounds and/or target proteins along with the experimentally verified bio-interaction information to generate predictive models. Lately, sophisticated machine learning techniques are applied in VS to elevate the predictive performance. The objective of this study is to examine and discuss the recent applications of machine learning techniques in VS, including deep learning, which became highly popular after giving rise to epochal developments in the fields of computer vision and natural language processing. The past 3 years have witnessed an unprecedented amount of research studies considering the application of deep learning in biomedicine, including computational drug discovery. In this review, we first describe the main instruments of VS methods, including compound and protein features (i.e. representations and descriptors), frequently used libraries and toolkits for VS, bioactivity databases and gold-standard data sets for system training and benchmarking. We subsequently review recent VS studies with a strong emphasis on deep learning applications. Finally, we discuss the present state of the field, including the current challenges and suggest future directions. We believe that this survey will provide insight to the researchers working in the field of computational drug discovery in terms of comprehending and developing novel bio-prediction methods.

Prediction and Experimental Confirmation of Novel Peripheral Cannabinoid-1 Receptor Antagonists

Small molecules targeting peripheral CB1 receptors have therapeutic potential in a variety of disorders including obesity-related, hormonal, and metabolic abnormalities, while avoiding the psychoactive effects in the central nervous system. We applied our in-house algorithm, iterative stochastic elimination, to produce a ligand-based model that distinguishes between CB1R antagonists and random molecules by physicochemical properties only. We screened ∼2 million commercially available molecules and found that about 500 of them are potential candidates to antagonize the CB1R. We applied a few criteria for peripheral activity and narrowed that set down to 30 molecules, out of which 15 could be purchased. Ten out of those 15 showed good affinity to the CB1R and two of them with nanomolar affinities (K_i of ∼400 nM). The eight molecules with top affinities were tested for activity: two compounds were pure antagonists, and five others were inverse agonists. These molecules are now being examined in vivo for their peripheral versus central distribution and subsequently will be tested for their effects on obesity in small animals.

Cardioprotective Melatonin: Translating from Proof-of-Concept Studies to Therapeutic Use

In this review we summarized the actual clinical data for a cardioprotective therapeutic role of melatonin, listed melatonin and its agonists in different stages of development, and evaluated the melatonin cardiovascular target tractability and prediction using machine learning on ChEMBL. To date, most clinical trials investigating a cardioprotective therapeutic role of melatonin are in phase 2a. Selective melatonin receptor agonists Tasimelteon, Ramelteon, and combined melatonergic-serotonin Agomelatine, and other agonists with registered structures in CHEMBL were not yet investigated as cardioprotective or cardiovascular drugs. As drug-able for these therapeutic targets, melatonin receptor agonists have the benefit over melatonin of well-characterized pharmacologic profiles and extensive safety data. Recent reports of the X-ray crystal structures of MT1 and MT2 receptors shall lead to the development of highly selective melatonin receptor agonists. Predictive models using machine learning could help to identify cardiovascular targets for melatonin. Selecting ChEMBL scores > 4.5 in cardiovascular assays, and melatonin scores > 4, we obtained 284 records from 162 cardiovascular assays carried out with 80 molecules with predicted or measured melatonin activity. Melatonin activities (agonistic or antagonistic) found in these experimental cardiovascular assays and models include arrhythmias, coronary and large vessel contractility, and hypertension. Preclinical proof-of-concept and early clinical studies (phase 2a) suggest a cardioprotective benefit from melatonin in various heart diseases. However, larger phase 3 randomized interventional studies are necessary to establish melatonin and its agonists’ actions as cardioprotective therapeutic agents.

Pioneering topological methods for network-based drug-target prediction by exploiting a brain-network self-organization theory

The bipartite network representation of the drug–target interactions (DTIs) in a biosystem enhances understanding of the drugs’ multifaceted action modes, suggests therapeutic switching for approved drugs and unveils possible side effects. As experimental testing of DTIs is costly and time-consuming, computational predictors are of great aid. Here, for the first time, state-of-the-art DTI supervised predictors custom-made in network biology were compared—using standard and innovative validation frameworks—with unsupervised pure topological-based models designed for general-purpose link prediction in bipartite networks. Surprisingly, our results show that the bipartite topology alone, if adequately exploited by means of the recently proposed local-community-paradigm (LCP) theory—initially detected in brain-network topological self-organization and afterwards generalized to any complex network—is able to suggest highly reliable predictions, with comparable performance with the state-of-the-art-supervised methods that exploit additional (non-topological, for instance biochemical) DTI knowledge. Furthermore, a detailed analysis of the novel predictions revealed that each class of methods prioritizes distinct true interactions; hence, combining methodologies based on diverse principles represents a promising strategy to improve drug–target discovery. To conclude, this study promotes the power of bio-inspired computing, demonstrating that simple unsupervised rules inspired by principles of topological self-organization and adaptiveness arising during learning in living intelligent systems (like the brain) can efficiently equal perform complicated algorithms based on advanced, supervised and knowledge-based engineering.