In silico prediction of mitochondrial toxicity of chemicals using machine learning methods

Developing a GNN-based AI model to predict mitochondrial toxicity using the bagging method

Mitochondrial toxicity has been implicated in the development of various toxicities, including hepatotoxicity. Therefore, mitochondrial toxicity has become a major screening factor in the early discovery phase of drug development. Several models have been developed to predict mitochondrial toxicity based on chemical structures. However, they only provide a binary classification of positive or negative results and do not provide the substructures that contribute to a positive decision. Therefore, we developed an artificial intelligence (AI) model to predict mitochondrial toxicity and visualize structural alerts. To construct the model, we used the open-source software library kMoL, which employs a graph neural network approach that allows learning from chemical structure data. We also utilized the integrated gradient method, which enables the visualization of substructures that contribute to positive results. The dataset used to construct the AI model exhibited a significant imbalance, with significantly more negative than positive data. To address this, we employed the bagging method, which resulted in a model with high predictive performance, as evidenced by an F1 score of 0.839. This model can also be used to visualize substructures that contribute to mitochondrial toxicity using the integrated gradient method. Our AI model predicts mitochondrial toxicity based on chemical structures and may contribute to screening mitochondrial toxicity in the early stages of drug discovery.

EDC-Predictor: A Novel Strategy for Prediction of Endocrine-Disrupting Chemicals by Integrating Pharmacological and Toxicological Profiles

Identification of endocrine-disrupting chemicals (EDCs) is crucial in the reduction of human health risks. However, it is hard to do so because of the complex mechanisms of the EDCs. In this study, we propose a novel strategy named EDC-Predictor to integrate pharmacological and toxicological profiles for the prediction of EDCs. Different from conventional methods that only focus on a few nuclear receptors (NRs), EDC-Predictor considers more targets. It uses computational target profiles from network-based and machine learning-based methods to characterize compounds, including both EDCs and non-EDCs. The best model constructed by these target profiles outperformed those models by molecular fingerprints. In a case study to predict NR-related EDCs, EDC-Predictor showed a wider applicability domain and higher accuracy than four previous tools. Another case study further demonstrated that EDC-Predictor could predict EDCs targeting other proteins rather than NRs. Finally, a free web server was developed to make EDC prediction easier (http://lmmd.ecust.edu.cn/edcpred/). In summary, EDC-Predictor would be a powerful tool in EDC prediction and drug safety assessment.

Predicting GPR40 Agonists with A Deep Learning-Based Ensemble Model

Recent studies have identified G protein-coupled receptor 40 (GPR40) as a promising target for treating type 2 diabetes mellitus, and GPR40 agonists have several superior effects over other hypoglycemic drugs, including cardiovascular protection and suppression of glucagon levels. In this study, we constructed an up-to-date GPR40 ligand dataset for training models and performed a systematic optimization of the ensemble model, resulting in a powerful ensemble model (ROC AUC: 0.9496) for distinguishing GPR40 agonists and non-agonists. The ensemble model is divided into three layers, and the optimization process is carried out in each layer. We believe that these results will prove helpful for both the development of GPR40 agonists and ensemble models. All the data and models are available on GitHub. (https://github.com/Jiamin-Yang/ensemble_model).

Modeling and Insights into the Structural Characteristics of Chemical Mitochondrial Toxicity

Mitochondria are the energy metabolism center of cells and are involved in a number of other processes, such as cell differentiation and apoptosis, signal transduction, and regulation of cell cycle and cell proliferation. It is of great significance to evaluate the mitochondrial toxicity of drugs and other chemicals. In the present study, we aimed to propose easily available artificial intelligence (AI) models for the prediction of chemical mitochondrial toxicity and investigate the structural characteristics with the analysis of molecular properties and structural alerts. The consensus model achieved good predictive results with high total accuracy at 87.21% for validation sets. The models can be accessed freely via https://ochem.eu/article/158582. Besides, several commonly used chemical properties were significantly different between chemicals with and without mitochondrial toxicity. We also detected the structural alerts (SAs) responsible for mitochondrial toxicity and integrated them into the web-server SApredictor (www.sapredictor.cn). The study may provide useful tools for in silico estimation of mitochondrial toxicity and be helpful to understand the mechanisms of mitochondrial toxicity.

In silico modeling-based new alternative methods to predict drug and herb-induced liver injury: A review

New approach methods (NAMs) have been developed to predict a wide range of toxicities through innovative technologies. Liver injury is one of the most extensively studied endpoints due to its severity and frequency, occurring among populations that consume drugs or dietary supplements. In this review, we focus on recent developments of in silico modeling for liver injury prediction using deep learning and in vitro data based on adverse outcome pathways (AOPs). Despite these models being mainly developed using datasets generated from drug-like molecules, they were also applied to the prediction of hepatotoxicity caused by herbal products. As deep learning has achieved great success in many different fields, advanced machine learning algorithms have been actively applied to improve the accuracy of in silico models. Additionally, the development of liver AOPs, combined with big data in toxicology, has been valuable in developing in silico models with enhanced predictive performance and interpretability. Specifically, one approach involves developing structure-based models for predicting molecular initiating events of liver AOPs, while others use in vitro data with structure information as model inputs for making predictions. Even though liver injury remains a difficult endpoint to predict, advancements in machine learning algorithms and the expansion of in vitro databases with relevant biological knowledge have made a huge impact on improving in silico modeling for drug-induced liver injury prediction.

Machine Learning Models to Predict Cytochrome P450 2B6 Inhibitors and Substrates

Cytochrome P450 2B6 (CYP2B6) is responsible for the metabolism of ∼7% of marketed drugs. The in vitro drug interaction studies guidance for industry issued by the FDA stipulates that drug sponsors need to evaluate whether the investigated drugs interact with the major drug-metabolizing P450s including CYP2B6. Therefore, there has been greater attention to the development of predictive models for CYP2B6 inhibitors and substrates. In this study, conventional machine learning and deep learning models were developed to predict CYP2B6 inhibitors and substrates. Our results showed that the best CYP2B6 inhibitor model yielded the AUC values of 0.95 and 0.75 with the 10-fold cross-validation and the test set, respectively, and the best CYP2B6 substrate model produced the AUC values of 0.93 and 0.90 with the 10-fold cross-validation and the test set, respectively. The generalization ability of the CYP2B6 inhibitor and substrate models was assessed by using the external validation sets. Several significant substructural fragments relevant to CYP2B6 inhibitors and substrates were detected via frequency substructure analysis and information gain. In addition, the applicability domain of the models was defined by employing a nonparametric method based on the probability density distribution. We anticipate that our results would be useful for the prediction of potential CYP2B6 inhibitors and substrates in the early stage of drug discovery.

In Silico Prediction of Human Organ Toxicity via Artificial Intelligence Methods

Unpredicted human organ level toxicity remains one of the major reasons for drug clinical failure. There is a critical need for cost-efficient strategies in the early stages of drug development for human toxicity assessment. At present, artificial intelligence methods are popularly regarded as a promising solution in chemical toxicology. Thus, we provided comprehensive in silico prediction models for eight significant human organ level toxicity end points using machine learning, deep learning, and transfer learning algorithms. In this work, our results showed that the graph-based deep learning approach was generally better than the conventional machine learning models, and good performances were observed for most of the human organ level toxicity end points in this study. In addition, we found that the transfer learning algorithm could improve model performance for skin sensitization end point using source domain of in vivo acute toxicity data and in vitro data of the Tox21 project. It can be concluded that our models can provide useful guidance for the rapid identification of the compounds with human organ level toxicity for drug discovery.

XML-CIMT: Explainable Machine Learning (XML) Model for Predicting Chemical-Induced Mitochondrial Toxicity

Organ toxicity caused by chemicals is a serious problem in the creation and usage of chemicals such as medications, insecticides, chemical products, and cosmetics. In recent decades, the initiation and development of chemical-induced organ damage have been related to mitochondrial dysfunction, among several adverse effects. Recently, many drugs, for example, troglitazone, have been removed from the marketplace because of significant mitochondrial toxicity. As a result, it is an urgent requirement to develop in silico models that can reliably anticipate chemical-induced mitochondrial toxicity. In this paper, we have proposed an explainable machine-learning model to classify mitochondrially toxic and non-toxic compounds. After several experiments, the Mordred feature descriptor was shortlisted to be used after feature selection. The selected features used with the CatBoost learning algorithm achieved a prediction accuracy of 85% in 10-fold cross-validation and 87.1% in independent testing. The proposed model has illustrated improved prediction accuracy when compared with the existing state-of-the-art method available in the literature. The proposed tree-based ensemble model, along with the global model explanation, will aid pharmaceutical chemists in better understanding the prediction of mitochondrial toxicity.

Promotion of Bladder Cancer Cell Metastasis by 2-Mercaptobenzothiazole via Its Activation of Aryl Hydrocarbon Receptor Transcription: Molecular Dynamics Simulations, Cell-Based Assays, and Machine Learning-Driven Prediction

2-Mercaptobenzothiazole (MBT) is an industrial chemical widely used for rubber products, corrosion inhibitors, and polymer materials with multiple environmental and exposure pathways. A growing body of evidence suggests its potential bladder cancer (BC) risk as a public health concern; however, the molecular mechanism remains poorly understood. Herein, we demonstrate the activation of the aryl hydrocarbon receptor (AhR) by MBT and reveal key events in carcinogenesis associated with BC. MBT alters conformational changes of AhR ligand binding domain (LBD) as revealed by 500 ns molecular dynamics simulations and activates AhR transcription with upregulation of AhR-target genes CYP1A1 and CYP1B1 to approximately 1.5-fold. MBT upregulates the expression of MMP1, the cancer cell metastasis biomarker, to 3.2-fold and promotes BC cell invasion through an AhR-mediated manner. MBT is further revealed to induce differentially expressed genes (DEGs) most enriched in cancer pathways by transcriptome profiling. The exposure of MBT at environmentally relevant concentrations induces BC risk via AhR signaling disruption, transcriptome aberration, and malignant cell metastasis. A machine learning-based model with an AUC value of 0.881 is constructed to successfully predict 31 MBT analogues. Overall, we provide molecular insight into the BC risk of MBT and develop an effective tool for rapid screening of AhR agonists.

In Silico Prediction of Human and Rat Liver Microsomal Stability via Machine Learning Methods

Liver microsomal stability is an important property considered for the screening of drug candidates in the early stage of drug development. Determination of hepatic metabolic stability can be performed by an in vitro assay, but it requires quite a few resources and time. In recent years, machine learning methods have made much progress. Therefore, development of computational models to predict liver microsomal stability is highly desirable in the drug discovery process. In this study, the in silico classification models for the prediction of the metabolic stability of compounds in rat and human liver microsomes were constructed by the conventional machine learning and deep learning methods. The performance of the models was evaluated using the test and external sets. For the rat liver microsomes (RLM) stability, the best model yielded the AUC values of 0.84 and 0.71 on the test and external validation sets, respectively. For the human liver microsome (HLM) stability, the best model exhibited the AUC values of 0.86 and 0.77 on the test and external validation sets, respectively. In addition, several important substructure fragments were detected using information gain and frequency substructure analysis methods. The applicability domain of the models was defined using the Euclidean distance-based method. We anticipate that our results would be helpful for the prediction of liver microsomal stability of compounds in the early stage of drug discovery.

Integrating cell morphology with gene expression and chemical structure to aid mitochondrial toxicity detection

Mitochondrial toxicity is an important safety endpoint in drug discovery. Models based solely on chemical structure for predicting mitochondrial toxicity are currently limited in accuracy and applicability domain to the chemical space of the training compounds. In this work, we aimed to utilize both -omics and chemical data to push beyond the state-of-the-art. We combined Cell Painting and Gene Expression data with chemical structural information from Morgan fingerprints for 382 chemical perturbants tested in the Tox21 mitochondrial membrane depolarization assay. We observed that mitochondrial toxicants differ from non-toxic compounds in morphological space and identified compound clusters having similar mechanisms of mitochondrial toxicity, thereby indicating that morphological space provides biological insights related to mechanisms of action of this endpoint. We further showed that models combining Cell Painting, Gene Expression features and Morgan fingerprints improved model performance on an external test set of 244 compounds by 60% (in terms of F1 score) and improved extrapolation to new chemical space. The performance of our combined models was comparable with dedicated in vitro assays for mitochondrial toxicity. Our results suggest that combining chemical descriptors with biological readouts enhances the detection of mitochondrial toxicants, with practical implications in drug discovery.

Cell Painting, gene expression, and chemical structural data are used to examine the differences between mitochondrial toxicants and non-toxicants and enhance the detection of mitotoxic compounds for future drug discovery.

Cell Morphological Profiling Enables High-Throughput Screening for PROteolysis TArgeting Chimera (PROTAC) Phenotypic Signature

PROteolysis TArgeting Chimeras (PROTACs) use the ubiquitin-proteasome system to degrade a protein of interest for therapeutic benefit. Advances made in targeted protein degradation technology have been remarkable, with several molecules having moved into clinical studies. However, robust routes to assess and better understand the safety risks of PROTACs need to be identified, which is an essential step toward delivering efficacious and safe compounds to patients. In this work, we used Cell Painting, an unbiased high-content imaging method, to identify phenotypic signatures of PROTACs. Chemical clustering and model prediction allowed the identification of a mitotoxicity signature that could not be expected by screening the individual PROTAC components. The data highlighted the benefit of unbiased phenotypic methods for identifying toxic signatures and the potential to impact drug design.

Mitochondria as the Target of Hepatotoxicity and Drug-Induced Liver Injury: Molecular Mechanisms and Detection Methods

One of the major mechanisms of drug-induced liver injury includes mitochondrial perturbation and dysfunction. This is not a surprise, given that mitochondria are essential organelles in most cells, which are responsible for energy homeostasis and the regulation of cellular metabolism. Drug-induced mitochondrial dysfunction can be influenced by various factors and conditions, such as genetic predisposition, the presence of metabolic disorders and obesity, viral infections, as well as drugs. Despite the fact that many methods have been developed for studying mitochondrial function, there is still a need for advanced and integrative models and approaches more closely resembling liver physiology, which would take into account predisposing factors. This could reduce the costs of drug development by the early prediction of potential mitochondrial toxicity during pre-clinical tests and, especially, prevent serious complications observed in clinical settings.

Multitask CapsNet: An Imbalanced Data Deep Learning Method for Predicting Toxicants

Drug development has a high failure rate, with safety properties constituting a considerable challenge. To reduce risk, in silico tools, including various machine learning methods, have been applied for toxicity prediction. However, these approaches often confront a serious problem: the training data sets are usually biased (imbalanced positive and negative samples), which would result in model training difficulty and unsatisfactory prediction accuracy. Multitask networks obtained significantly better predictive accuracies than single-task methods, and capsule neural networks showed excellent performance in sparse data sets in previous studies. In this study, we developed a new multitask framework based on a capsule neural network (multitask CapsNet) to measure 12 different toxic effects simultaneously. We found that multitask CapsNet excelled in toxicity prediction and outperformed many other computational approaches using the multitask strategy. Only after training on biased data sets did multitask CapsNet achieve significantly improved prediction accuracy on the Tox21 Data Challenge, which gave the largest ratio of highest accuracy (8/12) among compared models. Our model gave a prediction accuracy of 96.6% for the target NR.PPAR.gamma, whose ratio of negative to positive samples was up to 36:1. These results suggested that multitask CapsNet could overcome the bias problems and would provide a novel, accurate, and efficient approach for predicting the toxicities of compounds.