Alarms about structural alerts

Identifying Substructures That Facilitate Compounds to Penetrate the Blood-Brain Barrier via Passive Transport Using Machine Learning Explainer Models

The local interpretable model-agnostic explanation (LIME) method was used to interpret two machine learning models of compounds penetrating the blood-brain barrier. The classification models, Random Forest, ExtraTrees, and Deep Residual Network, were trained and validated using the blood-brain barrier penetration dataset, which shows the penetrability of compounds in the blood-brain barrier. LIME was able to create explanations for such penetrability, highlighting the most important substructures of molecules that affect drug penetration in the barrier. The simple and intuitive outputs prove the applicability of this explainable model to interpreting the permeability of compounds across the blood-brain barrier in terms of molecular features. LIME explanations were filtered with a weight equal to or greater than 0.1 to obtain only the most relevant explanations. The results showed several structures that are important for blood-brain barrier penetration. In general, it was found that some compounds with nitrogenous substructures are more likely to permeate the blood-brain barrier. The application of these structural explanations may help the pharmaceutical industry and potential drug synthesis research groups to synthesize active molecules more rationally.

Expanding Predictive Capacities in Toxicology: Insights from Hackathon-Enhanced Data and Model Aggregation

In the realm of predictive toxicology for small molecules, the applicability domain of QSAR models is often limited by the coverage of the chemical space in the training set. Consequently, classical models fail to provide reliable predictions for wide classes of molecules. However, the emergence of innovative data collection methods such as intensive hackathons have promise to quickly expand the available chemical space for model construction. Combined with algorithmic refinement methods, these tools can address the challenges of toxicity prediction, enhancing both the robustness and applicability of the corresponding models. This study aimed to investigate the roles of gradient boosting and strategic data aggregation in enhancing the predictivity ability of models for the toxicity of small organic molecules. We focused on evaluating the impact of incorporating fragment features and expanding the chemical space, facilitated by a comprehensive dataset procured in an open hackathon. We used gradient boosting techniques, accounting for critical features such as the structural fragments or functional groups often associated with manifestations of toxicity.

Prediction of Cytochrome P450 Inhibition Using a Deep Learning Approach and Substructure Pattern Recognition

Cytochrome P450 (CYP) is a family of enzymes that are responsible for about 75% of all metabolic reactions. Among them, CYP1A2, CYP2C9, CYP2C19, CYP2D6, and CYP3A4 participate in the metabolism of most drugs and mediate many adverse drug reactions. Therefore, it is necessary to estimate the chemical inhibition of Cytochrome P450 enzymes in drug discovery and the food industry. In the past few decades, many computational models have been reported, and some provided good performance. However, there are still several issues that should be resolved for these models, such as single isoform, models with unbalanced performance, lack of structural characteristics analysis, and poor availability. In the present study, the deep learning models based on python using the Keras framework and TensorFlow were developed for the chemical inhibition of each CYP isoform. These models were established based on a large data set containing 85715 compounds extracted from the PubChem bioassay database. On external validation, the models provided good AUC values with 0.97, 0.94, 0.94, 0.96, and 0.94 for CYP1A2, CYP2C9, CYP2C19, CYP2D6, and CYP3A4, respectively. The models can be freely accessed on the Web server named CYPi-DNNpredictor (cypi.sapredictor.cn), and the codes for the model were made open source in the Supporting Information. In addition, we also analyzed the structural characteristics of chemicals with CYP450 inhibition and detected the structural alerts (SAs), which should be responsible for the inhibition. The SAs were also made available online, named CYPi-SAdetector (cypisa.sapredictor.cn). The models can be used as a powerful tool for the prediction of CYP450 inhibitors, and the SAs should provide useful information for the mechanisms of Cytochrome P450 inhibition.

Pyruvate dehydrogenase is a potential mitochondrial off-target for gentamicin based on in silico predictions and in vitro inhibition studies

During the drug development process, organ toxicity leads to an estimated failure of one-third of novel chemical entities. Drug-induced toxicity is increasingly associated with mitochondrial dysfunction, but identifying the underlying molecular mechanisms remains a challenge. Computational modeling techniques have proven to be a good tool in searching for drug off-targets. Here, we aimed to identify mitochondrial off-targets of the nephrotoxic drugs tenofovir and gentamicin using different in silico approaches (KRIPO, ProBis and PDID). Dihydroorotate dehydrogenase (DHODH) and pyruvate dehydrogenase (PDH) were predicted as potential novel off-target sites for tenofovir and gentamicin, respectively. The predicted targets were evaluated in vitro, using (colorimetric) enzymatic activity measurements. Tenofovir did not inhibit DHODH activity, while gentamicin potently reduced PDH activity. In conclusion, the use of in silico methods appeared a valuable approach in predicting PDH as a mitochondrial off-target of gentamicin. Further research is required to investigate the contribution of PDH inhibition to overall renal toxicity of gentamicin.

PBT assessment of chemicals detected in effluent of wastewater treatment plants by suspected screening analysis

Wastewater treatment plants (WWTPs) are the major sources of contaminants discharged into downstream water bodies. Profiling the contaminants in effluent of WWTPs is crucial to assess the potential eco-risks toward downstream organisms. To this end, this study investigated the contaminants in effluent of 10 WWTPs locating in 10 cities of Yangtze River delta region of China by suspected screening analysis. Further, the persistence, bioaccumulation, toxicity (PBT) and the characteristics sub-structures of PBT-like chemicals were analyzed. Totally, 704 chemicals including 155 chemical products, 31 food additives, 52 natural substances, 112 personal care products, 123 pesticides, 192 pharmaceuticals, 17 hormones and 22 others were found. The results of PBT analysis suggested that 42 chemicals (5.97% among the detected chemicals in WWTPs) were with PBT property. Among them, 31 contaminants were not reported previously. 9 characteristics sub-structures (N-methyleneisobutylamine, 1-naphthaldehyde, 2,3,3-trimethylcyclohexene, cyclohexanol, N-sec-butyl-n-propylamine, (5E)-2,6-dimethylocta-1,5-diene, 2-ethylphenol, pentadecane and 6-methoxyhexane) were found for PBT-like chemicals. The sub-structures of highly linear alkyl partially explained the significantly higher PBT score for personal care products. Present study provides fundamental information on PBT properties of contaminants in effluent of WWTPs, which will benefit to prioritize contaminants with high concerns in effluent of WWTPs.

Lies and Liabilities: Computational Assessment of High-Throughput Screening Hits to Identify Artifact Compounds

Hits from high-throughput screening (HTS) of chemical libraries are often false positives due to their interference with assay detection technology. In response, we generated the largest publicly available library of chemical liabilities and developed "Liability Predictor," a free web tool to predict HTS artifacts. More specifically, we generated, curated, and integrated HTS data sets for thiol reactivity, redox activity, and luciferase (firefly and nano) activity and developed and validated quantitative structure-interference relationship (QSIR) models to predict these nuisance behaviors. The resulting models showed 58-78% external balanced accuracy for 256 external compounds per assay. QSIR models developed and validated herein identify nuisance compounds among experimental hits more reliably than do popular PAINS filters. Both the models and the curated data sets were implemented in "Liability Predictor," publicly available at https://liability.mml.unc.edu/. "Liability Predictor" may be used as part of chemical library design or for triaging HTS hits.

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.

Augmenting Expert Knowledge-Based Toxicity Alerts by Statistically Mined Molecular Fragments

Structural alerts are molecular substructures assumed to be associated with molecular initiating events in various toxic effects and an integral part of in silico toxicology. However, alerts derived using the knowledge of human experts often suffer from a lack of predictivity, specificity, and satisfactory coverage. In this work, we present a method to build hybrid QSAR models by combining expert knowledge-based alerts and statistically mined molecular fragments. Our objective was to find out if the combination is better than the individual systems. Lasso regularization-based variable selection was applied on combined sets of knowledge-based alerts and molecular fragments, but the variable elimination was only allowed to happen on the molecular fragments. We tested the concept on three toxicity end points, i.e., skin sensitization, acute Daphnia toxicity, and Ames mutagenicity, which covered both classification and regression problems. Results showed the predictive performance of such hybrid models is, indeed, better than the models based solely on expert alerts or statistically mined fragments alone. The method also enables the discovery of activating and mitigating/deactivating features for toxicity alerts and the identification of new alerts, thereby reducing false positive and false negative outcomes commonly associated with generic alerts and alerts with poor coverage, respectively.

The Good, The Bad, and The Perplexing: Structural Alerts and Read-Across for Predicting Skin Sensitization Using Human Data

In our earlier work (Golden et al., 2021), we showed 70-80% accuracies for several skin sensitization computational tools using human data. Here, we expanded the data set using the NICEATM human skin sensitization database to create a final data set of 1355 discrete chemicals (largely negative, ∼70%). Using this expanded data set, we analyzed model performance and evaluated mispredictions using Toxtree (v 3.1.0), OECD QSAR Toolbox (v 4.5), VEGA's (1.2.0 BETA) CAESAR (v 2.1.7), and a k-nearest-neighbor (kNN) classification approach. We show that the accuracy on this data set was lower than previous estimates, with balanced accuracies being 63% and 65% for Toxtree and OECD QSAR Toolbox, respectively, 46% for VEGA, and 59% for a kNN approach, with the lower accuracy likely due to the higher percentage of nonsensitizing chemicals. Two hundred eighty seven chemicals were mispredicted by both Toxtree and OECD QSAR Toolbox, which was approximately 20% of the entire data set, and 84% of these were false positives. The absence or presence of metabolic simulation in OECD QSAR Toolbox made no overall difference. While Toxtree is known for overpredicting, 60% of the chemicals in the data set had no alert for skin sensitization, and a substantial number of these chemicals were in fact sensitizers, pointing to sensitization mechanisms not recognized by Toxtree. Interestingly, we observed that chemicals with more than one Toxtree alert were more likely to be nonsensitizers. Finally, a kNN approach tended to mispredict different chemicals than either OECD QSAR Toolbox or Toxtree, suggesting that there was additional information to be garnered from a kNN approach. Overall, the results demonstrate that while there is merit in structural alerts as well as QSAR or read-across approaches (perhaps even more so in their combination), additional improvement will require a more nuanced understanding of mechanisms of skin sensitization.

A high-throughput screening campaign against PFKFB3 identified potential inhibitors with novel scaffolds

The growth of solid tumors depends on tumor vascularization and the endothelial cells (ECs) that line the lumen of blood vessels. ECs generate a large fraction of ATP through glycolysis, and elevation of their glycolytic activity is associated with angiogenic behavior in solid tumors. 6-Phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3 (PFKFB3) positively regulates glycolysis via fructose-2/6-bisphosphate, the product of its kinase activity. Partial inhibition of glycolysis in tumor ECs by targeting PFKFB3 normalizes the otherwise abnormal tumor vessels, thereby reducing metastasis and improving the outcome of chemotherapy. Although a limited number of tool compounds exist, orally available PFKFB3 inhibitors are unavailable. In this study we conducted a high-throughput screening campaign against the kinase activity of PFKFB3, involving 250,240 chemical compounds. A total of 507 initial hits showing >50% inhibition at 20 µM were identified, 66 of them plus 1 analog from a similarity search consistently displayed low IC50 values (<10 µM). In vitro experiments yielded 22 nontoxic hits that suppressed the tube formation of primary human umbilical vein ECs at 10 µM. Of them, 15 exhibited binding affinity to PFKFB3 in surface plasmon resonance assays, including 3 (WNN0403-E003, WNN1352-H007 and WNN1542-F004) that passed the pan-assay interference compounds screening without warning flags. This study provides potential leads to the development of new PFKFB3 inhibitors.

Structure optimization of new tumor-selective Passerini α-acyloxy carboxamides as Caspase-3/7 activators

Selective elimination of tumors has always been the mainstay of oncology research. The on-going research underlying the cellular apoptotic mechanisms reveal caspases activation, especially the key effector caspase-3, as a personalized tumor-selective therapeutic strategy. Our continued research protocol has exploited new optimized Passerini α-acyloxy carboxamides as efficient apoptotic inducers via caspase-3/7 dependent mechanism with highly selective anticancer profiles. The adopted design rationale relied on excluding structural alerts of previous leads, while merging various pharmacophoric motifs of natural and synthetic caspase activators via optimized one-pot Passerini reaction conditions. The prepared compounds resulting from Passerini reaction were screened for their cytotoxic activities against colorectal Caco-2 and liver HepG-2 cancer cells compared to normal fibroblasts utilizing MTT assay. Notably, all compounds exhibited promising low-range submicromolar IC₅₀ against the studied cancer cell lines, with outstanding tumor selectivity (SI values up to 266). Hence, they were superior to 5-fluorouracil. Notably, 7a, 7g, and 7j conferred the highest potencies against Caco-2 and HepG-2 cells and were selected for further mechanistic studies. Caspas-3/7 activation assay of the hit compounds and flow cytometric analysis of the treated apoptotic cancer cells demonstrated their significant caspase activation potential (up to 4.2 folds) and apoptotic induction capacities (up to 58.7%). Further assessment of Bcl2 expression was performed being a physiological caspase-3 substrate. Herein, the three studied Passerini adducts were able to downregulate Bcl2 in the treated Caco-2 cells. Importantly, the mechanistic studies results of the three hits echoed their preliminary MTT antiproliferative potencies data highlighting their caspase-3 dependent apoptotic induction. Finally, the in silico predicted physicochemical and pharmacokinetic profiles, as well as ligand efficiency metrics were drug-like.

MS2Tox Machine Learning Tool for Predicting the Ecotoxicity of Unidentified Chemicals in Water by Nontarget LC-HRMS

To achieve water quality objectives of the zero pollution action plan in Europe, rapid methods are needed to identify the presence of toxic substances in complex water samples. However, only a small fraction of chemicals detected with nontarget high-resolution mass spectrometry can be identified, and fewer have ecotoxicological data available. We hypothesized that ecotoxicological data could be predicted for unknown molecular features in data-rich high-resolution mass spectrometry (HRMS) spectra, thereby circumventing time-consuming steps of molecular identification and rapidly flagging molecules of potentially high toxicity in complex samples. Here, we present MS2Tox, a machine learning method, to predict the toxicity of unidentified chemicals based on high-resolution accurate mass tandem mass spectra (MS²). The MS2Tox model for fish toxicity was trained and tested on 647 lethal concentration (LC₅₀) values from the CompTox database and validated for 219 chemicals and 420 MS² spectra from MassBank. The root mean square error (RMSE) of MS2Tox predictions was below 0.89 log-mM, while the experimental repeatability of LC₅₀ values in CompTox was 0.44 log-mM. MS2Tox allowed accurate prediction of fish LC₅₀ values for 22 chemicals detected in water samples, and empirical evidence suggested the right directionality for another 68 chemicals. Moreover, by incorporating structural information, e.g., the presence of carbonyl-benzene, amide moieties, or hydroxyl groups, MS2Tox outperforms baseline models that use only the exact mass or log K_OW.

PreS/MD: Predictor of Sensitization hazard for chemical substances released from Medical Devices

In the United States, a pre-market regulatory submission for any medical device that comes into contact with either a patient or the clinical practitioner must include an adequate toxicity evaluation of chemical substances that can be released from the device during its intended use. These substances, also referred to as extractables and leachables, must be evaluated for their potential to induce sensitization/allergenicity, which traditionally has been done in animal assays such as the guinea pig maximization test (GPMT). However, advances in basic and applied science are continuously presenting opportunities to employ New Approach Methodologies (NAMs), including computational methods which, when qualified, could replace animal testing methods to support regulatory submissions. Herein, we developed a new computational tool for rapid and accurate prediction of the GPMT outcome that we have named PreS/MD (Predictor of Sensitization for Medical Devices). To enable model development, we (i) collected, curated, and integrated the largest publicly available dataset for GPMT results; (ii) succeeded in developing externally predictive (balanced accuracy of 70-74% as evaluated by both 5-fold external cross-validation and testing of novel compounds) Quantitative Structure-Activity Relationships (QSAR) models for GPMT using machine learning algorithms, including Deep Learning; and (iii) developed a publicly accessible web portal integrating PreS/MD models that can predict GPMT outcomes for any molecule of interest. We expect that PreS/MD will be used by both industry and regulatory scientists in medical device safety assessments and help replace, reduce, or refine the use of animals in toxicity testing. PreS/MD is freely available at https://presmd.mml.unc.edu/.

SApredictor: An Expert System for Screening Chemicals Against Structural Alerts

The rapid and accurate evaluation of chemical toxicity is of great significance for estimation of chemical safety. In the past decades, a great number of excellent computational models have been developed for chemical toxicity prediction. But most machine learning models tend to be “black box”, which bring about poor interpretability. In the present study, we focused on the identification and collection of structural alerts (SAs) responsible for a series of important toxicity endpoints. Then, we carried out effective storage of these structural alerts and developed a web-server named SApredictor (www.sapredictor.cn) for screening chemicals against structural alerts. People can quickly estimate the toxicity of chemicals with SApredictor, and the specific key substructures which cause the chemical toxicity will be intuitively displayed to provide valuable information for the structural optimization by medicinal chemists.

BayeshERG: a robust, reliable and interpretable deep learning model for predicting hERG channel blockers

Unintended inhibition of the human ether-à-go-go-related gene (hERG) ion channel by small molecules leads to severe cardiotoxicity. Thus, hERG channel blockage is a significant concern in the development of new drugs. Several computational models have been developed to predict hERG channel blockage, including deep learning models; however, they lack robustness, reliability and interpretability. Here, we developed a graph-based Bayesian deep learning model for hERG channel blocker prediction, named BayeshERG, which has robust predictive power, high reliability and high resolution of interpretability. First, we applied transfer learning with 300 000 large data in initial pre-training to increase the predictive performance. Second, we implemented a Bayesian neural network with Monte Carlo dropout to calibrate the uncertainty of the prediction. Third, we utilized global multihead attentive pooling to augment the high resolution of structural interpretability for the hERG channel blockers and nonblockers. We conducted both internal and external validations for stringent evaluation; in particular, we benchmarked most of the publicly available hERG channel blocker prediction models. We showed that our proposed model outperformed predictive performance and uncertainty calibration performance. Furthermore, we found that our model learned to focus on the essential substructures of hERG channel blockers via an attention mechanism. Finally, we validated the prediction results of our model by conducting in vitro experiments and confirmed its high validity. In summary, BayeshERG could serve as a versatile tool for discovering hERG channel blockers and helping maximize the possibility of successful drug discovery. The data and source code are available at our GitHub repository (https://github.com/GIST-CSBL/BayeshERG).

A Computational Approach for Respiratory Hazard Identification of Flavor Chemicals in Tobacco Products

Flavor chemicals contribute to the appeal and toxicity of tobacco products, including electronic nicotine delivery systems (ENDS). The assortment of flavor chemicals available for use in tobacco products is extensive. In this study, a chemistry-driven computational approach was used to evaluate flavor chemicals based on intrinsic hazardous structures and reactivity of chemicals. A large library of 3012 unique flavor chemicals was compiled from publicly available information. Next, information was computed and collated based on their (1) physicochemical properties, (2) global harmonization system (GHS) health hazard classification, (3) structural alerts linked to the chemical's reactivity, instability, or toxicity, and (4) common substructure shared with FDA's harmful and potentially harmful constituents (HPHCs) flavor chemicals that are respiratory toxicants. Computational analysis of the constructed flavor library flagged 638 chemicals with GHS classified respiratory health hazards, 1079 chemicals with at least one structural alert, and 2297 chemicals with substructural similarity to FDA's established and proposed list of HPHCs. A subsequent analysis was performed on a subset of 173 chemicals in the flavor library that are respiratory health hazards, contain structural alerts as well as flavor HPHC substructures. Four general toxicophore structures with an increased potential for respiratory toxicity were then identified. In summary, computational methods are efficient tools for hazard identification and understanding structure-toxicity relationship. With appropriate context of use and interpretation, in silico methods may provide scientific evidence to support toxicological evaluations of chemicals in or emitted from tobacco products.

Promoting GAINs (Give Attention to Limitations in Assays) over PAINs Alerts: no PAINS, more GAINs

Many concepts and guidelines in medicinal chemistry have been introduced to aid in successful drug discovery and development. An example is the concept of Pan-Assay Interference Compounds (PAINS) and the elimination of such nuisance compounds from high-throughput screening (HTS) libraries. PAINs, along with other guidelines in medicinal chemistry, are like double-edged swords. If used appropriately, they may be beneficial for drug discovery and development. However, rigid and blind use of such concepts can hinder productivity. In this perspective, we introduce GAINS (give attention to limitations in assays) and highlight its relevance for successful drug discovery.

STopTox: An in Silico Alternative to Animal Testing for Acute Systemic and Topical Toxicity

BACKGROUND

Modern chemical toxicology is facing a growing need to Reduce, Refine, and Replace animal tests (Russell 1959) for hazard identification. The most common type of animal assays for acute toxicity assessment of chemicals used as pesticides, pharmaceuticals, or in cosmetic products is known as a "6-pack" battery of tests, including three topical (skin sensitization, skin irritation and corrosion, and eye irritation and corrosion) and three systemic (acute oral toxicity, acute inhalation toxicity, and acute dermal toxicity) end points.

METHODS

We compiled, curated, and integrated, to the best of our knowledge, the largest publicly available data sets and developed an ensemble of quantitative structure-activity relationship (QSAR) models for all six end points. All models were validated according to the Organisation for Economic Co-operation and Development (OECD) QSAR principles, using data on compounds not included in the training sets.

RESULTS

In addition to high internal accuracy assessed by cross-validation, all models demonstrated an external correct classification rate ranging from 70% to 77%. We established a publicly accessible Systemic and Topical chemical Toxicity (STopTox) web portal (https://stoptox.mml.unc.edu/) integrating all developed models for 6-pack assays.

CONCLUSIONS

We developed STopTox, a comprehensive collection of computational models that can be used as an alternative to in vivo 6-pack tests for predicting the toxicity hazard of small organic molecules. Models were established following the best practices for the development and validation of QSAR models. Scientists and regulators can use the STopTox portal to identify putative toxicants or nontoxicants in chemical libraries of interest. https://doi.org/10.1289/EHP9341.

In Silico Models for Hepatotoxicity

In this chapter, we review the state of the art of predicting human hepatotoxicity using in silico techniques. There has been significant progress in this area over the past 20 years but there are still some challenges ahead. Principally, these challenges are our partial understanding of a very complex biochemical system and our ability to emulate that in a predictive capacity. Here, we provide an overview of the published modeling approaches in this area to date and discuss their design, strengths and weaknesses. It is interesting to note the diversity in modeling approaches, whether they be statistical algorithms or evidenced-based approaches including structural alerts and pharmacophore models. Irrespective of modeling approach, it appears a common theme of access to appropriate, relevant, and high-quality data is a limitation to all and is likely to continue to be the focus of future research.

Screening of 6000 Compounds for Uncoupling Activity: A Comparison Between a Mechanistic Biophysical Model and the Structural Alert Profiler Mitotox

Protonophoric uncoupling of phosphorylation is an important factor when assessing chemicals for their toxicity, and has recently moved into focus in pharmaceutical research with respect to the treatment of diseases such as cancer, diabetes, or obesity. Reliably identifying uncoupling activity is thus a valuable goal. To that end, we screened more than 6000 anionic compounds for in vitro uncoupling activity, using a biophysical model based on ab initio COSMO-RS input parameters with the molecular structure as the only external input. We combined these results with a model for baseline toxicity (narcosis). Our model identified more than 1250 possible uncouplers in the screening dataset, and identified possible new uncoupler classes such as thiophosphoric acids. When tested against 423 known uncouplers and 612 known inactive compounds in the dataset, the model reached a sensitivity of 83% and a specificity of 96%. In a direct comparison, it showed a similar specificity than the structural alert profiler Mitotox (97%), but much higher sensitivity than Mitotox (47%). The biophysical model thus allows for a more accurate screening for uncoupling activity than existing structural alert profilers. We propose to use our model as a complementary tool to screen large datasets for protonophoric uncoupling activity in drug development and toxicity assessment.

Comprehensive strategies of machine-learning-based quantitative structure-activity relationship models

Early quantitative structure-activity relationship (QSAR) technologies have unsatisfactory versatility and accuracy in fields such as drug discovery because they are based on traditional machine learning and interpretive expert features. The development of Big Data and deep learning technologies significantly improve the processing of unstructured data and unleash the great potential of QSAR. Here we discuss the integration of wet experiments (which provide experimental data and reliable verification), molecular dynamics simulation (which provides mechanistic interpretation at the atomic/molecular levels), and machine learning (including deep learning) techniques to improve QSAR models. We first review the history of traditional QSAR and point out its problems. We then propose a better QSAR model characterized by a new iterative framework to integrate machine learning with disparate data input. Finally, we discuss the application of QSAR and machine learning to many practical research fields, including drug development and clinical trials.

Machine Learning for Predicting Risk of Drug-Induced Autoimmune Diseases by Structural Alerts and Daily Dose

An effective approach for assessing a drug’s potential to induce autoimmune diseases (ADs) is needed in drug development. Here, we aim to develop a workflow to examine the association between structural alerts and drugs-induced ADs to improve toxicological prescreening tools. Considering reactive metabolite (RM) formation as a well-documented mechanism for drug-induced ADs, we investigated whether the presence of certain RM-related structural alerts was predictive for the risk of drug-induced AD. We constructed a database containing 171 RM-related structural alerts, generated a dataset of 407 AD- and non-AD-associated drugs, and performed statistical analysis. The nitrogen-containing benzene substituent alerts were found to be significantly associated with the risk of drug-induced ADs (odds ratio = 2.95, p = 0.0036). Furthermore, we developed a machine-learning-based predictive model by using daily dose and nitrogen-containing benzene substituent alerts as the top inputs and achieved the predictive performance of area under curve (AUC) of 70%. Additionally, we confirmed the reactivity of the nitrogen-containing benzene substituent aniline and related metabolites using quantum chemistry analysis and explored the underlying mechanisms. These identified structural alerts could be helpful in identifying drug candidates that carry a potential risk of drug-induced ADs to improve their safety profiles.

Simplex representation of molecular structure as universal QSAR/QSPR tool

We review the development and application of the Simplex approach for the solution of various QSAR/QSPR problems. The general concept of the simplex method and its varieties are described. The advantages of utilizing this methodology, especially for the interpretation of QSAR/QSPR models, are presented in comparison to other fragmentary methods of molecular structure representation. The utility of SiRMS is demonstrated not only in the standard QSAR/QSPR applications, but also for mixtures, polymers, materials, and other complex systems. In addition to many different types of biological activity (antiviral, antimicrobial, antitumor, psychotropic, analgesic, etc.), toxicity and bioavailability, the review examines the simulation of important properties, such as water solubility, lipophilicity, as well as luminescence, and thermodynamic properties (melting and boiling temperatures, critical parameters, etc.). This review focuses on the stereochemical description of molecules within the simplex approach and details the possibilities of universal molecular stereo-analysis and stereochemical configuration description, along with stereo-isomerization mechanism and molecular fragment "topography" identification.

In-silico evaluation of Malawi essential medicines and reactive metabolites for potential drug-induced toxicities

BACKGROUND

Drug-induced toxicity is one of the problems that have negatively impacted on the well-being of populations throughout the world, including Malawi. It results in unnecessary hospitalizations, retarding the development of the country. This study assessed the Malawi Essential Medicines List (MEML) for structural alerts and reactive metabolites with the potential for drug-induced toxicities.

METHODS

This in-silico screening study used StopTox, ToxAlerts and LD-50 values toxicity models to assess the MEML drugs. A total of 296 drugs qualified for the analysis (those that had defined chemical structures) and were screened in each software programme. Each model had its own toxicity endpoints and the models were compared for consensus of their results.

RESULTS

In the StopTox model, 86% of the drugs had potential to cause at least one toxicity including 55% that had the potential of causing eye irritation and corrosion. In ToxAlerts, 90% of the drugs had the potential of causing at least one toxicity and 72% were found to be potentially reactive, unstable and toxic. In LD-50, 70% of the drugs were potentially toxic. Model consensus evaluation results showed that the highest consensus was observed between ToxAlerts and StopTox (80%). The overall consensus amongst the three models was 57% and statistically significant (p < 0.05).

CONCLUSIONS

A large number of drugs had the potential to cause various systemic toxicities. But the results need to be interpreted cautiously since the clinical translation of QSAR-based predictions depends on many factors. In addition, inconsistencies have been reported between screening results amongst different models.

Mechanism-Driven Read-Across of Chemical Hepatotoxicants Based on Chemical Structures and Biological Data

Hepatotoxicity is a leading cause of attrition in the drug development process. Traditional preclinical and clinical studies to evaluate hepatotoxicity liabilities are expensive and time consuming. With the advent of critical advancements in high-throughput screening, there has been a rapid accumulation of in vitro toxicity data available to inform the risk assessment of new pharmaceuticals and chemicals. To this end, we curated and merged all available in vivo hepatotoxicity data obtained from the literature and public resources, which yielded a comprehensive database of 4089 compounds that includes hepatotoxicity classifications. After dividing the original database of chemicals into modeling and test sets, PubChem assay data were automatically extracted using an in-house data mining tool and clustered based on relationships between structural fragments and cellular responses in in vitro assays. The resultant PubChem assay clusters were further investigated. During the cross-validation procedure, the biological data obtained from several assay clusters exhibited high predictivity of hepatotoxicity and these assays were selected to evaluate the test set compounds. The read-across results indicated that if a new compound contained specific identified chemical fragments (ie, Molecular Initiating Event) and showed active responses in the relevant selected PubChem assays, there was potential for the chemical to be hepatotoxic in vivo. Furthermore, several mechanisms that might contribute to toxicity were derived from the modeling results including alterations in nuclear receptor signaling and inhibition of DNA repair. This modeling strategy can be further applied to the investigation of other complex chemical toxicity phenomena (eg, developmental and reproductive toxicities) as well as drug efficacy.

What's in a Name? Drug Nomenclature and Medicinal Chemistry Trends using INN Publications

The World Health Organization assigns international nonproprietary names (INN), also known as common names, to compounds upon request from drug developers. Structures of INNs are publicly available and represent a source, albeit underused, to understand trends in drug research and development. Here, we explain how a common drug name is composed and analyze chemical entities from 2000 to 2021. In the analysis, we describe some changes that intertwine chemical structure, newer therapeutic targets (e.g., kinases), including a significant increase in the use of fluorine and of heterocycles, and some other evolutionary modifications, such as the progressive increase in molecular weight. Alongside these, small signs of change can be spotted, such as the rise in spirocyclic scaffolds and small rings and the emergence of unconventional structural moieties that might forecast the future to come.

Deep Learning Model for Identifying Critical Structural Motifs in Potential Endocrine Disruptors

This paper aims to identify structural motifs within a molecule that contribute the most toward a chemical being an endocrine disruptor. We have developed a deep neural network-based toolkit toward this aim. The trained model can virtually assess a synthetic chemical's potential to be an endocrine disruptor using machine-readable molecular representation, simplified molecular input line entry system (SMILES). Our proposed toolkit is a multilabel or multioutput classification model that combines both convolution and long short-term memory (LSTM) architectures. The toolkit leverages the advantages of an active learning-based framework that combines multiple sources of data. Class activation maps (CAMs) generated from the feature-extraction layers can identify the structural alerts and the chemical environment that determines the specificity of the structural alerts.

PySmash: Python package and individual executable program for representative substructure generation and application

BACKGROUND

Substructure screening is widely applied to evaluate the molecular potency and ADMET properties of compounds in drug discovery pipelines, and it can also be used to interpret QSAR models for the design of new compounds with desirable physicochemical and biological properties. With the continuous accumulation of more experimental data, data-driven computational systems which can derive representative substructures from large chemical libraries attract more attention. Therefore, the development of an integrated and convenient tool to generate and implement representative substructures is urgently needed.

RESULTS

In this study, PySmash, a user-friendly and powerful tool to generate different types of representative substructures, was developed. The current version of PySmash provides both a Python package and an individual executable program, which achieves ease of operation and pipeline integration. Three types of substructure generation algorithms, including circular, path-based and functional group-based algorithms, are provided. Users can conveniently customize their own requirements for substructure size, accuracy and coverage, statistical significance and parallel computation during execution. Besides, PySmash provides the function for external data screening.

CONCLUSION

PySmash, a user-friendly and integrated tool for the automatic generation and implementation of representative substructures, is presented. Three screening examples, including toxicophore derivation, privileged motif detection and the integration of substructures with machine learning (ML) models, are provided to illustrate the utility of PySmash in safety profile evaluation, therapeutic activity exploration and molecular optimization, respectively. Its executable program and Python package are available at https://github.com/kotori-y/pySmash.

Benchmarking the mechanisms of frequent hitters: limitation of PAINS alerts

In 2010, the pan-assay interference compounds (PAINS) rule was proposed to identify false-positive compounds, especially frequent hitters (FHs), in biological screening campaigns, and has rapidly become an essential component in drug design. However, the specific mechanisms remain unknown, and the result validation and follow-up processing schemes are still unclear. In this review, a large benchmark collection of >600,000 compounds sourced from databases and the literature, including six common false-positive mechanisms, was used to evaluate the detection ability of PAINS. In addition, 400 million purchasable molecules from the ZINC database were also applied to PAINS screening. The results indicate that the PAINS rule is not suitable for the screening of all types of false-positive results and needs more improvement.

Modeling the Bioactivation and Subsequent Reactivity of Drugs

Electrophilically reactive drug metabolites are implicated in many adverse drug reactions. In this mechanism-termed bioactivation-metabolic enzymes convert drugs into reactive metabolites that often conjugate to nucleophilic sites within biological macromolecules like proteins. Toxic metabolite-product adducts induce severe immune responses that can cause sometimes fatal disorders, most commonly in the form of liver injury, blood dyscrasia, or the dermatologic conditions toxic epidermal necrolysis and Stevens-Johnson syndrome. This study models four of the most common metabolic transformations that result in bioactivation: quinone formation, epoxidation, thiophene sulfur-oxidation, and nitroaromatic reduction, by synthesizing models of metabolism and reactivity. First, the metabolism models predict the formation probabilities of all possible metabolites among the pathways studied. Second, the exact structures of these metabolites are enumerated. Third, using these structures, the reactivity model predicts the reactivity of each metabolite. Finally, a feedfoward neural network converts the metabolism and reactivity predictions to a bioactivation prediction for each possible metabolite. These bioactivation predictions represent the joint probability that a metabolite forms and that this metabolite subsequently conjugates to protein or glutathione. Among molecules bioactivated by these pathways, we predicted the correct pathway with an AUC accuracy of 89.98%. Furthermore, the model predicts whether molecules will be bioactivated, distinguishing bioactivated and nonbioactivated molecules with 81.06% AUC. We applied this algorithm to withdrawn drugs. The known bioactivation pathways of alclofenac and benzbromarone were identified by the algorithm, and high probability bioactivation pathways not yet confirmed were identified for safrazine, zimelidine, and astemizole. This bioactivation model-the first of its kind that jointly considers both metabolism and reactivity-enables drug candidates to be quickly evaluated for a toxicity risk that often evades detection during preclinical trials. The XenoSite bioactivation model is available at http://swami.wustl.edu/xenosite/p/bioactivation.

Confidence in Inactive and Active Predictions from Structural Alerts

Having a measure of confidence in computational predictions of biological activity from in silico tools is vital when making predictions for new chemicals, for example, in chemical risk assessment. Where predictions of biological activity are used as an indicator of a potential hazard, false-negative predictions are the most concerning prediction; however, assigning confidence in inactive predictions is particularly challenging. How can one confidently identify the absence of activating features? In this study, we present methods for assigning confidence to both active and inactive predictions from structural alerts for protein-binding molecular initiating events (MIEs). Structural alerts were derived through an iterative statistical method. Confidence in the activity predictions is assigned by measuring the Tanimoto similarity between Morgan fingerprints of chemicals in the test set to relevant chemicals in the training set, and suitable cutoff values have been defined to give different confidence categories. To avoid a potential compound series bias in the test set and hence overestimate the performance of the method, we measured the biological activity of 27 compounds with 24 proteins, which gave us an additional 648 experimental measurements; many of the measurements are currently nonexistent in the literature and databases. This data set was complemented with newly measured biological activities published in ChEMBL25 and formed a combined independent validation data set. Applying the confidence categories to the computational predictions for the new data leads to the identification of chemicals for which one should be confident of either an inactive or active prediction, allowing model predictions to be used responsibly.

Mapping Chemical Respiratory Sensitization: How Useful Are Our Current Computational Tools?

Chemical respiratory sensitization is an immunological process that manifests clinically mostly as occupational asthma and is responsible for 1 in 6 cases of adult asthma, although this may be an underestimate of the prevalence, as it is under-diagnosed. Occupational asthma results in unemployment for roughly one-third of those affected due to severe health issues. Despite its high prevalence, chemical respiratory sensitization is difficult to predict, as there are currently no validated models and the mechanisms are not entirely understood, creating a significant challenge for regulatory bodies and industry alike. The Adverse Outcome Pathway (AOP) for respiratory sensitization is currently incomplete. However, some key events have been identified, and there is overlap with the comparatively well-characterized AOP for dermal sensitization. Because of this, and the fact that dermal sensitization is often assessed by in vivo, in chemico, or in silico methods, regulatory bodies are defaulting to the dermal sensitization status of chemicals as a proxy for respiratory sensitization status when evaluating chemical safety. We identified a data set of known human respiratory sensitizers, which we used to investigate the accuracy of a structural alert model, Toxtree, designed for skin sensitization and the Centre for Occupational and Environmental Health (COEH)'s model, a model developed specifically for occupational asthma. While both models had a reasonable level of accuracy, the COEH model achieved the highest balanced accuracy at 76%; when the models agreed, the overall accuracy was 87%. There were important differences between the models: Toxtree had superior performance for some structural alerts and some categories of well-characterized skin sensitizers, while the COEH model had high accuracy in identifying sensitizers that lacked identified skin sensitization reactivity domains. Overall, both models achieved respectable accuracy. However, neither model addresses potency, which, along with data quality, remains a hurdle, and the field must prioritize these issues to move forward.

ChemFLuo: a web-server for structure analysis and identification of fluorescent compounds

BACKGROUND

Fluorescent detection methods are indispensable tools for chemical biology. However, the frequent appearance of potential fluorescent compound has greatly interfered with the recognition of compounds with genuine activity. Such fluorescence interference is especially difficult to identify as it is reproducible and possesses concentration-dependent characteristic. Therefore, the development of a credible screening tool to detect fluorescent compounds from chemical libraries is urgently needed in early stages of drug discovery.

RESULTS

In this study, we developed a webserver ChemFLuo for fluorescent compound detection, based on two large and high-quality training datasets containing 4906 blue and 8632 green fluorescent compounds. These molecules were used to construct a group of prediction models based on the combination of three machine learning algorithms and seven types of molecular representations. The best blue fluorescence prediction model achieved with balanced accuracy (BA) = 0.858 and area under the receiver operating characteristic curve (AUC) = 0.931 for the validation set, and BA = 0.823 and AUC = 0.903 for the test set. The best green fluorescence prediction model achieved the prediction accuracy with BA = 0.810 and AUC = 0.887 for the validation set, and BA = 0.771 and AUC = 0.852 for the test set. Besides prediction model, 22 blue and 16 green representative fluorescent substructures were summarized for the screening of potential fluorescent compounds. The comparison with other fluorescence detection tools and theapplication to external validation sets and large molecule libraries have demonstrated the reliability of prediction model for fluorescent compound detection.

CONCLUSION

ChemFLuo is a public webserver to filter out compounds with undesirable fluorescent properties, which will benefit the design of high-quality chemical libraries for drug discovery. It is freely available at http://admet.scbdd.com/chemfluo/index/.

The Influence of Structural Patterns on Acute Aquatic Toxicity of Organic Compounds

Investigation of the influence of molecular structure of different organic compounds on acute toxicity towards Fathead minnow, Daphnia magna, and Tetrahymena pyriformis has been carried out using 2D simplex representation of molecular structure and two modelling methods: Random Forest (RF) and Gradient Boosting Machine (GBM). Suitable QSAR (Quantitative Structure - Activity Relationships) models were obtained. The study was focused on QSAR models interpretation. The aim of the study was to develop a set of structural fragments that simultaneously consistently increase toxicity toward Fathead minnow, Daphnia magna, Tetrahymena pyriformis. The interpretation allowed to gain more details about known toxicophores and to propose new fragments. The results obtained made it possible to rank the contributions of molecular fragments to various types of toxicity to aquatic organisms. This information can be used for molecular optimization of chemicals. According to the results of structural interpretation, the most significant common mechanisms of the toxic effect of organic compounds on Fathead minnow, Daphnia magna and Tetrahymena pyriformis are reactions of nucleophilic substitution and inhibition of oxidative phosphorylation in mitochondria. In addition acetylcholinesterase and voltage-gated ion channel of Fathead minnow and Daphnia magna are important targets for toxicants. The on-line version of the OCHEM expert system (https://ochem.eu) were used for a comparative QSAR investigation. The proposed QSAR models comply with the OECD principles and can be used to reliably predict acute toxicity of organic compounds towards Fathead minnow, Daphnia magna and Tetrahymena pyriformis with allowance for applicability domain estimation.

Scopy: an integrated negative design python library for desirable HTS/VS database design

BACKGROUND

High-throughput screening (HTS) and virtual screening (VS) have been widely used to identify potential hits from large chemical libraries. However, the frequent occurrence of 'noisy compounds' in the screened libraries, such as compounds with poor drug-likeness, poor selectivity or potential toxicity, has greatly weakened the enrichment capability of HTS and VS campaigns. Therefore, the development of comprehensive and credible tools to detect noisy compounds from chemical libraries is urgently needed in early stages of drug discovery.

RESULTS

In this study, we developed a freely available integrated python library for negative design, called Scopy, which supports the functions of data preparation, calculation of descriptors, scaffolds and screening filters, and data visualization. The current version of Scopy can calculate 39 basic molecular properties, 3 comprehensive molecular evaluation scores, 2 types of molecular scaffolds, 6 types of substructure descriptors and 2 types of fingerprints. A number of important screening rules are also provided by Scopy, including 15 drug-likeness rules (13 drug-likeness rules and 2 building block rules), 8 frequent hitter rules (four assay interference substructure filters and four promiscuous compound substructure filters), and 11 toxicophore filters (five human-related toxicity substructure filters, three environment-related toxicity substructure filters and three comprehensive toxicity substructure filters). Moreover, this library supports four different visualization functions to help users to gain a better understanding of the screened data, including basic feature radar chart, feature-feature-related scatter diagram, functional group marker gram and cloud gram.

CONCLUSION

Scopy provides a comprehensive Python package to filter out compounds with undesirable properties or substructures, which will benefit the design of high-quality chemical libraries for drug design and discovery. It is freely available at https://github.com/kotori-y/Scopy.

Pred-Skin: A Web Portal for Accurate Prediction of Human Skin Sensitizers

Safety assessment is an essential component of the regulatory acceptance of industrial chemicals. Previously, we have developed a model to predict the skin sensitization potential of chemicals for two assays, the human patch test and murine local lymph node assay, and implemented this model in a web portal. Here, we report on the substantially revised and expanded freely available web tool, Pred-Skin version 3.0. This up-to-date version of Pred-Skin incorporates multiple quantitative structure-activity relationship (QSAR) models developed with in vitro, in chemico, and mice and human in vivo data, integrated into a consensus naïve Bayes model that predicts human effects. Individual QSAR models were generated using skin sensitization data derived from human repeat insult patch tests, human maximization tests, and mouse local lymph node assays. In addition, data for three validated alternative methods, the direct peptide reactivity assay, KeratinoSens, and the human cell line activation test, were employed as well. Models were developed using open-source tools and rigorously validated according to the best practices of QSAR modeling. Predictions obtained from these models were then used to build a naïve Bayes model for predicting human skin sensitization with the following external prediction accuracy: correct classification rate (89%), sensitivity (94%), positive predicted value (91%), specificity (84%), and negative predicted value (89%). As an additional assessment of model performance, we identified 11 cosmetic ingredients known to cause skin sensitization but were not included in our training set, and nine of them were accurately predicted as sensitizers by our models. Pred-Skin can be used as a reliable alternative to animal tests for predicting human skin sensitization.

Effect of the structural factors of organic compounds on the acute toxicity toward Daphnia magna

The acute toxicity of organic compounds towards Daphina magna was subjected to QSAR analysis. The two-dimensional simplex representation of molecular structure (2D SiRMS) and the support vector machine (SVM), gradient boosting (GBM) methods were used to develop QSAR models. Adequate regression QSAR models were developed for incubation of 24 h. Their interpretation allowed us to quantitatively describe and rank the well-known toxicophores, to refine their molecular surroundings, and to distinguish the structural derivatives of the fragments that significantly contribute to the acute toxicity (LC₅₀) of organic compounds towards D. magna. Based on the results of the interpretation of the regression models, a molecular design (modification) of highly toxic compounds was performed in order to reduce their hazard. In addition, acceptable classification QSAR models were developed to reliably predict the following mode of action (MOA): specific and non-specific toxicity of organic compounds towards D. magna. When interpreting these models, we were able to determine the structural fragments and the physicochemical characteristics of molecules that are responsible for the manifestation of one of the modes of action. The on-line version of the OCHEM expert system (https://ochem.eu), HYBOT descriptors, and the random forest and SVM methods were used for a comparative QSAR investigation.

QSAR without borders

Prediction of chemical bioactivity and physical properties has been one of the most important applications of statistical and more recently, machine learning and artificial intelligence methods in chemical sciences. This field of research, broadly known as quantitative structure-activity relationships (QSAR) modeling, has developed many important algorithms and has found a broad range of applications in physical organic and medicinal chemistry in the past 55+ years. This Perspective summarizes recent technological advances in QSAR modeling but it also highlights the applicability of algorithms, modeling methods, and validation practices developed in QSAR to a wide range of research areas outside of traditional QSAR boundaries including synthesis planning, nanotechnology, materials science, biomaterials, and clinical informatics. As modern research methods generate rapidly increasing amounts of data, the knowledge of robust data-driven modelling methods professed within the QSAR field can become essential for scientists working both within and outside of chemical research. We hope that this contribution highlighting the generalizable components of QSAR modeling will serve to address this challenge.

Data science-driven analyses of drugs inducing hypertension as an adverse effect

The utilization of approved medication is a requisite to combat certain diseases for health; however, the undesirable adverse effects (AEs) due to medication are generally unavoidable. Hypertension is one of such AEs resulting from approved medication in which blood pressure in the arteries gets elevated and is a risk factor for several diseases including heart and kidney failure. HTs are the approved drugs that can cause hypertension as an AE. Here, the goal of the study is to investigate the structural and functional diversities of HTs. In our quest to unravel the structural parameters of the HTs, a systematic analysis of the HTs having a different number and type of ring systems was conducted. The cellular component, molecular function and biological processes adopted by the gene products were analyzed. Moreover, our systematically done analysis suggests that all the target families are active in a common pathway, that is, nerve transmission. A comparison of the selected structural and functional aspect of HTs with anti-hypertensives suggests that HTs follow certain structural and functional features in spite of many possibilities. Our study provides a promising methodology that considers the influence of structural diversity of AE causing agents on a functional perspective for precursory clinical decision making. This could be extended to explore the structural and functional trends that are adopted by agents causing certain diseases or AEs.

Using Machine Learning Methods and Structural Alerts for Prediction of Mitochondrial Toxicity

Over the last few years more and more organ and idiosyncratic toxicities were linked to mitochondrial toxicity. Despite well-established assays, such as the seahorse and Glucose/Galactose assay, an in silico approach to mitochondrial toxicity is still feasible, particularly when it comes to the assessment of large compound libraries. Therefore, in silico approaches could be very beneficial to indicate hazards early in the drug development pipeline. By combining multiple endpoints, we derived the largest so far published dataset on mitochondrial toxicity. A thorough data analysis shows that molecules causing mitochondrial toxicity can be distinguished by physicochemical properties. Finally, the combination of machine learning and structural alerts highlights the suitability for in silico risk assessment of mitochondrial toxicity.

In silico toxicology: From structure–activity relationships towards deep learning and adverse outcome pathways

In silico toxicology is an emerging field. It gains increasing importance as research is aiming to decrease the use of animal experiments as suggested in the 3R principles by Russell and Burch. In silico toxicology is a means to identify hazards of compounds before synthesis, and thus in very early stages of drug development. For chemical industries, as well as regulatory agencies it can aid in gap‐filling and guide risk minimization strategies. Techniques such as structural alerts, read‐across, quantitative structure–activity relationship, machine learning, and deep learning allow to use in silico toxicology in many cases, some even when data is scarce. Especially the concept of adverse outcome pathways puts all techniques into a broader context and can elucidate predictions by mechanistic insights.

This article is categorized under:

Structure and Mechanism > Computational Biochemistry and Biophysics

Data Science > Chemoinformatics

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.

New in silico models to predict in vitro micronucleus induction as marker of genotoxicity

The evaluation of genotoxicity is a fundamental part of the safety assessment of chemicals due to the relevance of the potential health effects of genotoxicants. Among the testing methods available, the in vitro micronucleus assay with mammalian cells is one of the most used and required by regulations targeting several industrial sectors such as the cosmetic industry and food-related sectors. As an alternative to the testing methods, in recent years, lots in silico methods were developed to predict the genotoxicity of chemicals, including models for the Ames mutagenicity test, the in vitro chromosomal aberrations and the in vivo micronucleus assay. We developed several in silico models for the prediction of genotoxicity as reflected by the in vitro micronucleus assay. The resulting models include both statistical and knowledge-based models. The most promising model is the one based on fragments extracted with the SARpy platform. More than 100 structural alerts were extracted, including also fragments associated with the non-genotoxic activity. The model is characterized by high accuracy and the lowest false negative rate, making this tool suitable for chemical screening according to the regulators' needs. The SARpy model will be implemented on the VEGA platform (https://www.vegahub.eu) and will be freely available.

Towards reproducible computational drug discovery

The reproducibility of experiments has been a long standing impediment for further scientific progress. Computational methods have been instrumental in drug discovery efforts owing to its multifaceted utilization for data collection, pre-processing, analysis and inference. This article provides an in-depth coverage on the reproducibility of computational drug discovery. This review explores the following topics: (1) the current state-of-the-art on reproducible research, (2) research documentation (e.g. electronic laboratory notebook, Jupyter notebook, etc.), (3) science of reproducible research (i.e. comparison and contrast with related concepts as replicability, reusability and reliability), (4) model development in computational drug discovery, (5) computational issues on model development and deployment, (6) use case scenarios for streamlining the computational drug discovery protocol. In computational disciplines, it has become common practice to share data and programming codes used for numerical calculations as to not only facilitate reproducibility, but also to foster collaborations (i.e. to drive the project further by introducing new ideas, growing the data, augmenting the code, etc.). It is therefore inevitable that the field of computational drug design would adopt an open approach towards the collection, curation and sharing of data/code.

Structural Alerts and Random Forest Models in a Consensus Approach for Receptor Binding Molecular Initiating Events

A molecular initiating event (MIE) is the gateway to an adverse outcome pathway (AOP), a sequence of events ending in an adverse effect. In silico predictions of MIEs are a vital tool in a modern, mechanism-focused approach to chemical risk assessment. For 90 biological targets representing important human MIEs, structural alert-based models have been constructed with an automated procedure that uses Bayesian statistics to iteratively select substructures. These models give impressive average performance statistics (an average of 92% correct predictions across targets), significantly improving on previous models. Random Forest models have been constructed from physicochemical features for the same targets, giving similarly impressive performance statistics (93% correct predictions). A key difference between the models is interpretation of predictions-the structural alert models are transparent and easy to interpret, while Random Forest models can only identify the most important physicochemical features for making predictions. The two complementary models have been combined in a consensus model, improving performance compared to each individual model (94% correct predictions) and increasing confidence in predictions. Variation in model performance has been explained by calculating a modelability index (MODI), using Tanimoto coefficient between Morgan fingerprints to identify nearest neighbor chemicals. This work is an important step toward building confidence in the use of in silico tools for assessment of toxicity.

Evaluation of two in silico programs for predicting mutagenicity and carcinogenicity potential for 4-methylimidazole (4-MeI) and known metabolites

4-Methylimidazole (4-MeI) is a nitrogen-containing heterocyclic compound that is used in the manufacture of chemicals, dyes and pharmaceuticals and may be found in a variety of foods following formation during heating. The purpose of this study was to use two different in silico programs, CASE Ultra and Toxtree, to investigate potential structure-activity relationships in 4-MeI and its metabolites for mutagenicity and carcinogenicity, and combine that information with the available literature to draw conclusions regarding the strength of the predictions observed. Neither CASE Ultra nor Toxtree identified any structural alerts that were associated with mutagenic activity. Data for 4-MeI from a single study were used in the development of the CASE Ultra mouse and rat carcinogenicity models, but no additional similar structures were identified in the carcinogenicity model training set. One metabolite, 5-methylhydantoin, was predicted to be positive in the CASE Ultra carcinogenicity male and female mouse models; positive predictivity percentages of 60.9% and 73.7%, respectively. However, low structural similarity between 5-methylhydantoin and the compounds identified in the training set (<25%) decreases confidence in the positive prediction. Three metabolites were predicted to be positive in the CASE Ultra mouse micronucleus model, but again suffered from low structural similarity. Both limited structural similarity and inconsistent responses among the other clastogenicity models suggest that additional structurally similar compounds are needed to assess the predictive capacity of these alerts for biological activity of these compounds.

Human ether-à-go-go-related potassium channel: exploring SAR to improve drug design

hERG is best known as a primary anti-target, the inhibition of which is responsible for serious side effects. A renewed interest in hERG as a desired target, especially in oncology, was sparked because of its role in cellular proliferation and apoptosis. In this study, we survey the most recent advances regarding hERG by focusing on SAR in the attempt to elucidate, at a molecular level, off-target and on-target actions of potential hERG binders, which are highly promiscuous and largely varying in structure. Understanding the rationale behind hERG interactions and the molecular determinants of hERG activity is a real challenge and comprehension of this is of the utmost importance to prioritize compounds in early stages of drug discovery and to minimize cardiotoxicity attrition in preclinical and clinical studies.