Addressing the Metabolic Stability of Antituberculars through Machine Learning

Generation of new inhibitors of selected cytochrome P450 subtypes- In silico study

Physicochemical and pharmacokinetic compound profile has crucial impact on compound potency to become a future drug. Ligands with desired activity profile cannot be used for treatment if they are characterized by unfavourable physicochemical or ADMET properties. In the study, we consider metabolic stability and focus on selected subtypes of cytochrome P450 - proteins, which take part in the first phase of compound transformations in the organism. We develop a protocol for generation of new potential inhibitors of selected cytochrome isoforms. Its subsequent stages are composed of generation and assessment of new derivatives of known cytochrome inhibitors, docking and evaluation of the compound possible inhibition on the basis of the obtained ligand-protein complexes. Besides the library of new potential agents inhibiting particular cytochrome subtypes, we also prepare a graph neural network that predicts the change in activity for all modifications of the starting molecule. In addition, we perform a systematic statistical study on the influence of particular substitutions on the potential inhibition properties of generated compounds (both mono- and di-substitutions are considered), provide explanations of the inhibitory predictions and prepare an on-line visualization platform enabling manual inspection of the results. The developed methodology can greatly support the design of new cytochrome P450 inhibitors with the overarching goal of generation of new metabolically stable compounds. It enables instant evaluation of possible compound-cytochrome interactions and selection of ligands with the highest potential of possessing desired biological activity.

Random Forest Model Predictions Afford Dual-Stage Antimalarial Agents

The need for novel antimalarials is apparent given the continuing disease burden worldwide, despite significant drug discovery advances from the bench to the bedside. In particular, small-molecule agents with potent efficacy against both the liver and blood stages of Plasmodium parasite infection are critical for clinical settings as they would simultaneously prevent and treat malaria with a reduced selection pressure for resistance. While experimental screens for such dual-stage inhibitors have been conducted, the time and cost of these efforts limit their scope. Here, we have focused on leveraging machine learning approaches to discover novel antimalarials with such properties. A random forest modeling approach was taken to predict small molecules with in vitro efficacy versus liver-stage Plasmodium berghei parasites and a lack of human liver cell cytotoxicity. Empirical validation of the model was achieved with the realization of hits with liver-stage efficacy after prospective scoring of a commercial diversity library and consideration of structural diversity. A subset of these hits also demonstrated promising blood-stage Plasmodium falciparum efficacy. These 18 validated dual-stage antimalarials represent novel starting points for drug discovery and mechanism of action studies with significant potential for seeding a new generation of therapies.

How can SHAP values help to shape metabolic stability of chemical compounds?

Background

Computational methods support nowadays each stage of drug design campaigns. They assist not only in the process of identification of new active compounds towards particular biological target, but also help in the evaluation and optimization of their physicochemical and pharmacokinetic properties. Such features are not less important in terms of the possible turn of a compound into a future drug than its desired affinity profile towards considered proteins. In the study, we focus on metabolic stability, which determines the time that the compound can act in the organism and play its role as a drug. Due to great complexity of xenobiotic transformation pathways in the living organisms, evaluation and optimization of metabolic stability remains a big challenge.

Results

Here, we present a novel methodology for the evaluation and analysis of structural features influencing metabolic stability. To this end, we use a well-established explainability method called SHAP. We built several predictive models and analyse their predictions with the SHAP values to reveal how particular compound substructures influence the model’s prediction. The method can be widely applied by users thanks to the web service, which accompanies the article. It allows a detailed analysis of SHAP values obtained for compounds from the ChEMBL database, as well as their determination and analysis for any compound submitted by a user. Moreover, the service enables manual analysis of the possible structural modifications via the provision of analogous analysis for the most similar compound from the ChEMBL dataset.

Conclusions

To our knowledge, this is the first attempt to employ SHAP to reveal which substructural features are utilized by machine learning models when evaluating compound metabolic stability. The accompanying web service for metabolic stability evaluation can be of great help for medicinal chemists. Its significant usefulness is related not only to the possibility of assessing compound stability, but also to the provision of information about substructures influencing this parameter. It can assist in the design of new ligands with improved metabolic stability, helping in the detection of privileged and unfavourable chemical moieties during stability optimization. The tool is available at https://metstab-shap.matinf.uj.edu.pl/.

Bayesian Modeling and Intrabacterial Drug Metabolism Applied to Drug-Resistant Staphylococcus aureus

We present the application of Bayesian modeling to identify chemical tools and/or drug discovery entities pertinent to drug-resistant Staphylococcus aureus infections. The quinoline JSF-3151 was predicted by modeling and then empirically demonstrated to be active against in vitro cultured clinical methicillin- and vancomycin-resistant strains while also exhibiting efficacy in a mouse peritonitis model of methicillin-resistant S. aureus infection. We highlight the utility of an intrabacterial drug metabolism (IBDM) approach to probe the mechanism by which JSF-3151 is transformed within the bacteria. We also identify and then validate two mechanisms of resistance in S. aureus: one mechanism involves increased expression of a lipocalin protein, and the other arises from the loss of function of an azoreductase. The computational and experimental approaches, discovery of an antibacterial agent, and elucidated resistance mechanisms collectively hold promise to advance our understanding of therapeutic regimens for drug-resistant S. aureus.

Random Forest Model Prediction of Compound Oral Exposure in the Mouse

An early hurdle in the optimization of small-molecule chemical probes and drug discovery entities is the attainment of sufficient exposure in the mouse via oral administration of the compound. While computational approaches have attempted to predict molecular properties related to the mouse pharmacokinetic (PK) profile, we present herein a machine learning approach to specifically predict the oral exposure of a compound as measured in the mouse snapshot PK assay. A random forest workflow was found to produce the best cross-validation and external test set statistics after processing of the input data set and optimization of model features. The modeling approach should be useful to the chemical biology and drug discovery communities to predict this key molecular property and afford chemical entities of translational significance.

Pruned Machine Learning Models to Predict Aqueous Solubility

Solubility is a key metric for therapeutic compounds. Conversely, insoluble compounds cloud the accuracy of assays at all stages of chemical biology and drug discovery. Herein, we disclose naïve Bayesian classifier models to predict aqueous solubility. Publicly accessible aqueous solubility data were used to create two full, or nonpruned, training sets. These two sets were also combined to create a full fused set, and a training set comprised of a literature collation of solubility data was also considered as a reference. We tested different extents of data pruning on the training sets and constructed machine learning models that were evaluated with two independent, external test sets that contained compounds that were different from the training sets. The best pruned and fused model was significantly more accurate, in comparison to either the full model or the full fused model, with the prediction of these external test sets. By carefully removing data from the training set, less information can be used to create more accurate machine learning models for aqueous solubility. This knowledge and the curated training sets should prove useful to future machine learning approaches.

Machine Learning Platform to Discover Novel Growth Inhibitors of Neisseria gonorrhoeae

PURPOSE

To advance fundamental biological and translational research with the bacterium Neisseria gonorrhoeae through the prediction of novel small molecule growth inhibitors via naïve Bayesian modeling methodology.

METHODS

Inspection and curation of data from the publicly available ChEMBL web site for small molecule growth inhibition data of the bacterium Neisseria gonorrhoeae resulted in a training set for the construction of machine learning models. A naïve Bayesian model for bacterial growth inhibition was utilized in a workflow to predict novel antibacterial agents against this bacterium of global health relevance from a commercial library of >10⁵ drug-like small molecules. Follow-up efforts involved empirical assessment of the predictions and validation of the hits.

RESULTS

Specifically, two small molecules were found that exhibited promising activity profiles and represent novel chemotypes for agents against N. gonorrrhoeae.

CONCLUSIONS

This represents, to the best of our knowledge, the first machine learning approach to successfully predict novel growth inhibitors of this bacterium. To assist the chemical tool and drug discovery fields, we have made our curated training set available as part of the Supplementary Material and the Bayesian model is accessible via the web. Graphical Abstract.

Current status and future directions of high-throughput ADME screening in drug discovery

During the last decade high-throughput in vitro absorption, distribution, metabolism and excretion (HT-ADME) screening has become an essential part of any drug discovery effort of synthetic molecules. The conduct of HT-ADME screening has been "industrialized" due to the extensive development of software and automation tools in cell culture, assay incubation, sample analysis and data analysis. The HT-ADME assay portfolio continues to expand in emerging areas such as drug-transporter interactions, early soft spot identification, and ADME screening of peptide drug candidates. Additionally, thanks to the very large and high-quality HT-ADME data sets available in many biopharma companies, in silico prediction of ADME properties using machine learning has also gained much momentum in recent years. In this review, we discuss the current state-of-the-art practices in HT-ADME screening including assay portfolio, assay automation, sample analysis, data processing, and prediction model building. In addition, we also offer perspectives in future development of this exciting field.

Exploiting machine learning for end-to-end drug discovery and development

A variety of machine learning methods such as naive Bayesian, support vector machines and more recently deep neural networks are demonstrating their utility for drug discovery and development. These leverage the generally bigger datasets created from high-throughput screening data and allow prediction of bioactivities for targets and molecular properties with increased levels of accuracy. We have only just begun to exploit the potential of these techniques but they may already be fundamentally changing the research process for identifying new molecules and/or repurposing old drugs. The integrated application of such machine learning models for end-to-end (E2E) application is broadly relevant and has considerable implications for developing future therapies and their targeting.

Naïve Bayesian Models for Vero Cell Cytotoxicity

PURPOSE

To advance translational research of potential therapeutic small molecules against infectious microbes, the compounds must display a relative lack of mammalian cell cytotoxicity. Vero cell cytotoxicity (CC₅₀) is a common initial assay for this metric. We explored the development of naïve Bayesian models that can enhance the probability of identifying non-cytotoxic compounds.

METHODS

Vero cell cytotoxicity assays were identified in PubChem, reformatted, and curated to create a training set with 8741 unique small molecules. These data were used to develop Bayesian classifiers, which were assessed with internal cross-validation, external tests with a set of 193 compounds from our laboratory, and independent validation with an additional diverse set of 1609 unique compounds from PubChem.

RESULTS

Evaluation with independent, external test and validation sets indicated that cytotoxicity Bayesian models constructed with the ECFP_6 descriptor were more accurate than those that used FCFP_6 fingerprints. The best cytotoxicity Bayesian model displayed predictive power in external evaluations, according to conventional and chance-corrected statistics, as well as enrichment factors.

CONCLUSIONS

The results from external tests demonstrate that our novel cytotoxicity Bayesian model displays sufficient predictive power to help guide translational research. To assist the chemical tool and drug discovery communities, our curated training set is being distributed as part of the Supplementary Material. Graphical Abstract Naive Bayesian models have been trained with publically available data and offer a useful tool for chemical biology and drug discovery to select for small molecules with a high probability of exhibiting acceptably low Vero cell cytotoxicity.

Recent advances for identification of new scaffolds and drug targets for Mycobacterium tuberculosis

Tuberculosis (TB) is a leading cause of mortality and morbidity with an estimated 1.7 billion people latently infected with the pathogen worldwide. Clinically, TB infection presents itself as an asymptomatic infection, which gradually manifests to life threatening disease. The emergence of various drug resistant strains of Mycobacterium tuberculosis and lengthy duration of chemotherapy are major challenges in the field of TB drug development. Hence, there is an urgent need to develop scaffolds that possess a novel mechanism of action, can shorten the duration of therapy, and are active against both drug resistant and susceptible strains. In this review, we will discuss recent progress made in the field of TB drug development with emphasis on screening methods and drug targets from M. tuberculosis. The current review provides insights into mechanism of action of new scaffolds that are being evaluated in various stages of clinical trials. © 2018 IUBMB Life, 70(9):905-916, 2018.