Ligand-based prediction of hERG-mediated cardiotoxicity based on the integration of different machine learning techniques.

Drug-induced cardiotoxicity is a common side effect of drugs in clinical use or under postmarket surveillance and is commonly due to off-target interactions with the cardiac human-ether-a-go-go-related (hERG) potassium channel. Therefore, prioritizing drug candidates based on their hERG blocking potential is a mandatory step in the early preclinical stage of a drug discovery program. Herein, we trained and properly validated 30 ligand-based classifiers of hERG-related cardiotoxicity based on 7,963 curated compounds extracted by the freely accessible repository ChEMBL (version 25). Different machine learning algorithms were tested, namely, random forest, K-nearest neighbors, gradient boosting, extreme gradient boosting, multilayer perceptron, and support vector machine. The application of 1) the best practices for data curation, 2) the feature selection method VSURF, and 3) the synthetic minority oversampling technique (SMOTE) to properly handle the unbalanced data, allowed for the development of highly predictive models (BAMAX = 0.91, AUCMAX = 0.95). Remarkably, the undertaken temporal validation approach not only supported the predictivity of the herein presented classifiers but also suggested their ability to outperform those models commonly used in the literature. From a more methodological point of view, the study put forward a new computational workflow, freely available in the GitHub repository (https://github.com/PDelre93/hERG-QSAR), as valuable for building highly predictive models of hERG-mediated cardiotoxicity.

Monte Carlo Models for Sub-Chronic Repeated-Dose Toxicity: Systemic and Organ-Specific Toxicity.

The risk-characterization of chemicals requires the determination of repeated-dose toxicity (RDT). This depends on two main outcomes: the no-observed-adverse-effect level (NOAEL) and the lowest-observed-adverse-effect level (LOAEL). These endpoints are fundamental requirements in several regulatory frameworks, such as the Registration, Evaluation, Authorization and Restriction of Chemicals (REACH) and the European Regulation of 1223/2009 on cosmetics. The RDT results for the safety evaluation of chemicals are undeniably important; however, the in vivo tests are time-consuming and very expensive. The in silico models can provide useful input to investigate sub-chronic RDT. Considering the complexity of these endpoints, involving variable experimental designs, this non-testing approach is challenging and attractive. Here, we built eight in silico models for the NOAEL and LOAEL predictions, focusing on systemic and organ-specific toxicity, looking into the effects on the liver, kidney and brain. Starting with the NOAEL and LOAEL data for oral sub-chronic toxicity in rats, retrieved from public databases, we developed and validated eight quantitative structure-activity relationship (QSAR) models based on the optimal descriptors calculated by the Monte Carlo method, using the CORAL software. The results obtained with these models represent a good achievement, to exploit them in a safety assessment, considering the importance of organ-related toxicity.

In Silico Methods for Chromosome Damage.

Due to the link with serious adverse health effects, genotoxicity is an important toxicological endpoint in each regulatory setting with respect to human health, including for pharmaceuticals. To this extent, a compound potential to induce gene mutations as well as chromosome damage needs to be addressed. For chromosome damage, i.e., the induction of structural or numerical chromosome aberrations, several in vitro and in vivo test methods are available. In order to rapidly collect toxicological data without the need for test material, several in silico tools for chromosome damage have been developed over the last years. In this chapter, a battery of freely available in silico chromosome damage prediction tools for chromosome damage is applied on a dataset of pharmaceuticals. Examples of the different outcomes obtained with the in silico battery are provided and briefly discussed. Furthermore, results for coumarin are presented in more detail as a case study. Overall, it can be concluded that although they are in general less developed than those for mutagenicity, in silico tools for chromosome damage can provide valuable information, especially when combined in a battery.

QSAR models for soil ecotoxicity: Development and validation of models to predict reproductive toxicity of organic chemicals in the collembola Folsomia candida.

Soil pollution is a critical environmental challenge: the substances released in the soil can adversely affect humans and the ecosystem. Several bioassays were developed to investigate the soil ecotoxicity of chemicals with soil microbes, plants, invertebrates and vertebrates. The 28-day collembolan reproduction test with the springtail Folsomia candida is a recently introduced bioassay described by OECD guideline 232. Although the importance of springtails for maintaining soil quality, toxicity data for Collembola are still limited. We have developed two QSAR models for the prediction of reproductive toxicity induced by organic compounds in Folsomia candida using 28 days NOEC data. We assembled a dataset with the highest number of compounds available so far: 54 compounds were collected from publicly available sources, including plant protection products, reactive intermediates and industrial chemicals, household and cosmetic ingredients, drugs, environmental transformation products and polycyclic aromatic hydrocarbons. The models were developed using partial least squares regression (PLS) and the Monte Carlo technique with respectively the open source tools Small Dataset Modeler and CORAL software. Both QSAR models gave good predictive performance even though based on a small dataset, so they could serve for the ecological risk assessment of chemicals for terrestrial organisms.

Ecotoxicological QSAR modeling of the acute toxicity of organic compounds to the freshwater crustacean Thamnocephalus platyurus.

Growing interest in environmental toxicity assessment using Thamnocephalus platyurus as organism has led to an increased availability of acute toxicity data. Despite this growing interest in tests with this organism, however, to the best of our knowledge there are no computational models to predict the acute toxicity in T. platyurus. In view of the limited number of in silico models for this crustacean, we developed Quantitative Structure-Activity Relationship (QSAR) models for the prediction of acute toxicity towards T. platyurus, reflected by the 24h LC50, using publicly available data according to the ISO 14380:2011 guideline. Two models were developed following the principles of QSAR modeling recommended by the Organization for Economic Cooperation and Development (OECD). We used partial least squares and gradient boosting machine techniques, which gave encouraging statistical quality in our data set.

CATMoS: Collaborative Acute Toxicity Modeling Suite.

BACKGROUND: Humans are exposed to tens of thousands of chemical substances that need to be assessed for their potential toxicity. Acute systemic toxicity testing serves as the basis for regulatory hazard classification, labeling, and risk management. However, it is cost- and time-prohibitive to evaluate all new and existing chemicals using traditional rodent acute toxicity tests. In silico models built using existing data facilitate rapid acute toxicity predictions without using animals. OBJECTIVES: The U.S. Interagency Coordinating Committee on the Validation of Alternative Methods (ICCVAM) Acute Toxicity Workgroup organized an international collaboration to develop in silico models for predicting acute oral toxicity based on five different end points: Lethal Dose 50 (LD50 value, U.S. Environmental Protection Agency hazard (four) categories, Globally Harmonized System for Classification and Labeling hazard (five) categories, very toxic chemicals [LD50 (LD50≤50mg/kg)], and nontoxic chemicals (LD50>2,000mg/kg). METHODS: An acute oral toxicity data inventory for 11,992 chemicals was compiled, split into training and evaluation sets, and made available to 35 participating international research groups that submitted a total of 139 predictive models. Predictions that fell within the applicability domains of the submitted models were evaluated using external validation sets. These were then combined into consensus models to leverage strengths of individual approaches. RESULTS: The resulting consensus predictions, which leverage the collective strengths of each individual model, form the Collaborative Acute Toxicity Modeling Suite (CATMoS). CATMoS demonstrated high performance in terms of accuracy and robustness when compared with in vivo results. DISCUSSION: CATMoS is being evaluated by regulatory agencies for its utility and applicability as a potential replacement for in vivo rat acute oral toxicity studies. CATMoS predictions for more than 800,000 chemicals have been made available via the National Toxicology Program's Integrated Chemical Environment tools and data sets (ice.ntp.niehs.nih.gov). The models are also implemented in a free, standalone, open-source tool, OPERA, which allows predictions of new and untested chemicals to be made. https://doi.org/10.1289/EHP8495.

Integrated In Silico Models for the Prediction of No-Observed-(Adverse)-Effect Levels and Lowest-Observed-(Adverse)-Effect Levels in Rats for Sub-chronic Repeated-Dose Toxicity.

Repeated-dose toxicity (RDT) is a critical endpoint for hazard characterization of chemicals and is assessed to derive safe levels of exposure for human health. Here we present the first attempt to model simultaneously no-observed-(adverse)-effect level (NO(A)EL) and lowest-observed-(adverse)-effect level (LO(A)EL). Classification and regression models were derived based on rat sub-chronic repeated dose toxicity data for 327 compounds from the Fraunhofer RepDose database. Multi-category classification models were built for both NO(A)EL and LO(A)EL though a consensus of statistics- and fragment-based algorithms, while regression models were based on quantitative relationships between the endpoints and SMILES-based attributes. NO(A)EL and LO(A)EL models were integrated, and predictions were compared to exclude inconsistent values. This strategy improved the performance of single models, leading to R2 greater than 0.70, root-mean-square error (RMSE) lower than 0.60 (for regression models), and accuracy of 0.61-0.73 (for classification models) on the validation set, based on the endpoint and the threshold applied for selecting predictions. This study confirms the effectiveness of the modeling strategy presented here for assessing RDT of chemicals using in silico models.

Zebrafish AC50 modelling: (Q)SAR models to predict developmental toxicity in zebrafish embryo.

Developmental toxicity refers to the occurrence of adverse effects on a developing organism as a consequence of exposure to hazardous chemicals. The assessment of developmental toxicity has become relevant to the safety assessment process of chemicals. The zebrafish embryo developmental toxicology assay is an emerging test used to screen the teratogenic potential of chemicals and it is proposed as a promising test to replace teratogenic assays with animals. Supported by the increased availability of data from this test, the developmental toxicity assay with zebrafish has become an interesting endpoint for the in silico modelling. The purpose of this study was to build up quantitative structure-activity relationship (QSAR) models. In this work, new in silico models for the evaluation of developmental toxicity were built using a well-defined set of data from the ToxCastTM Phase I chemical library on the zebrafish embryo. Categorical and continuous QSAR models were built by gradient boosting machine learning and the Monte Carlo technique respectively, in accordance with Organization for Economic Co-operation and Development principles and their statistical quality was satisfactory. The classification model reached balanced accuracy 0.89 and Matthews correlation coefficient 0.77 on the test set. The regression model reached correlation coefficient R2 0.70 in external validation and leave-one-out cross-validated Q2 0.73 in internal validation.

Automated integration of structural, biological and metabolic similarities to improve read-across.

Read-across (RAX) is a popular data-gap filling technique that uses category and analogue approaches to predict toxicological endpoints for a target. Despite its increasing relevance, RAX relies on human expert judgement and lacks a reproducible and automated protocol. It also only relies on structural similarity for identifying the analogues, while other aspects are often neglected. In this paper, we propose an automated procedure for the selection of analogues for data gap-filling. Analogues were identified with a decision algorithm that integrates three similarity metrics, each considering different toxicologically relevant aspects (i.e., structural, biological and metabolic similarity). Structural filters based on the presence of maximum common substructures (MCS) and common functional groups were applied to narrow the chemical space for the analogues search. The procedure has been implemented as a workflow in KNIME and is freely available. The workflow provides informative tabular and graphical outputs to support toxicologists and risk assessors in drawing conclusion based on the RAX approach. The procedure has been validated for its predictive power on two datasets related to high-tier in vivo toxicological endpoints, i.e. human hepatotoxicity and drug-induced liver injury (DILI). The validation results gave good accuracy values (i.e., up to 0.79 for the binary hepatotoxicity classification and up to 0.67 for the three-class DILI classification) that were higher than those returned by RAX based on the sole use of structural similarity. Results confirmed the suitability of the procedure as a source of data to support regulatory decision-making.

SAR and QSAR modeling of a large collection of LD50 rat acute oral toxicity data.

The median lethal dose for rodent oral acute toxicity (LD50) is a standard piece of information required to categorize chemicals in terms of the potential hazard posed to human health after acute exposure. The exclusive use of in vivo testing is limited by the time and costs required for performing experiments and by the need to sacrifice a number of animals. (Quantitative) structure-activity relationships [(Q)SAR] proved a valid alternative to reduce and assist in vivo assays for assessing acute toxicological hazard. In the framework of a new international collaborative project, the NTP Interagency Center for the Evaluation of Alternative Toxicological Methods and the U.S. Environmental Protection Agency's National Center for Computational Toxicology compiled a large database of rat acute oral LD50 data, with the aim of supporting the development of new computational models for predicting five regulatory relevant acute toxicity endpoints. In this article, a series of regression and classification computational models were developed by employing different statistical and knowledge-based methodologies. External validation was performed to demonstrate the real-life predictability of models. Integrated modeling was then applied to improve performance of single models. Statistical results confirmed the relevance of developed models in regulatory frameworks, and confirmed the effectiveness of integrated modeling. The best integrated strategies reached RMSEs lower than 0.50 and the best classification models reached balanced accuracies over 0.70 for multi-class and over 0.80 for binary endpoints. Computed predictions will be hosted on the EPA's Chemistry Dashboard and made freely available to the scientific community.