New Models to Predict the Acute and Chronic Toxicities of Representative Species of the Main Trophic Levels of Aquatic Environments.

To assess the impact of chemicals on an aquatic environment, toxicological data for three trophic levels are needed to address the chronic and acute toxicities. The use of non-testing methods, such as predictive computational models, was proposed to avoid or reduce the need for animal models and speed up the process when there are many substances to be tested. We developed predictive models for Raphidocelis subcapitata, Daphnia magna, and fish for acute and chronic toxicities. The random forest machine learning approach gave the best results. The models gave good statistical quality for all endpoints. These models are freely available for use as individual models in the VEGA platform and for prioritization in JANUS software.

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.

Towards an Understanding of the Mode of Action of Human Aromatase Activity for Azoles through Quantum Chemical Descriptors-Based Regression and Structure Activity Relationship Modeling Analysis.

Aromatase is an enzyme member of the cytochrome P450 superfamily coded by the CYP19A1 gene. Its main action is the conversion of androgens into estrogens, transforming androstenedione into estrone and testosterone into estradiol. This enzyme is present in several tissues and it has a key role in the maintenance of the balance of androgens and estrogens, and therefore in the regulation of the endocrine system. With regard to chemical safety and human health, azoles, which are used as agrochemicals and pharmaceuticals, are potential endocrine disruptors due to their agonist or antagonist interactions with the human aromatase enzyme. This theoretical study investigated the active agonist and antagonist properties of "chemical classes of azoles" to determine the relationships of azole interaction with CYP19A1, using stereochemical and electronic properties of the molecules through classification and multilinear regression (MLR) modeling. The antagonist activities for the same substituent on diazoles and triazoles vary with its chemical composition and its position and both heterocyclic systems require aromatic substituents. The triazoles require the spherical shape and diazoles have to be in proper proportion of the branching index and the number of ring systems for the inhibition. Considering the electronic aspects, triazole antagonist activity depends on the electrophilicity index that originates from interelectronic exchange interaction (ωHF) and the LUMO energy ( E LUMO PM 7 ), and the diazole antagonist activity originates from the penultimate orbital ( E HOMONL PM 7 ) of diazoles. The regression models for agonist activity show that it is opposed by the static charges but favored by the delocalized charges on the diazoles and thiazoles. This study proposes that the electron penetration of azoles toward heme group decides the binding behavior and stereochemistry requirement for antagonist activity against CYP19A1 enzyme.

QSAR Models for Human Carcinogenicity: An Assessment Based on Oral and Inhalation Slope Factors.

Carcinogenicity is a crucial endpoint for the safety assessment of chemicals and products. During the last few decades, the development of quantitative structure-activity relationship ((Q)SAR) models has gained importance for regulatory use, in combination with in vitro testing or expert-based reasoning. Several classification models can now predict both human and rat carcinogenicity, but there are few models to quantitatively assess carcinogenicity in humans. To our knowledge, slope factor (SF), a parameter describing carcinogenicity potential used especially for human risk assessment of contaminated sites, has never been modeled for both inhalation and oral exposures. In this study, we developed classification and regression models for inhalation and oral SFs using data from the Risk Assessment Information System (RAIS) and different machine learning approaches. The models performed well in classification, with accuracies for the external set of 0.76 and 0.74 for oral and inhalation exposure, respectively, and r2 values of 0.57 and 0.65 in the regression models for oral and inhalation SFs in external validation. These models might therefore support regulators in (de)prioritizing substances for regulatory action and in weighing evidence in the context of chemical safety assessments. Moreover, these models are implemented on the VEGA platform and are now freely downloadable online.

Development of a neurotoxicity assay that is tuned to detect mitochondrial toxicants.

Many neurotoxicants affect energy metabolism in man, but currently available test methods may still fail to predict mito- and neurotoxicity. We addressed this issue using LUHMES cells, i.e., human neuronal precursors that easily differentiate into mature neurons. Within the NeuriTox assay, they have been used to screen for neurotoxicants. Our new approach is based on culturing the cells in either glucose or galactose (Glc-Gal-NeuriTox) as the main carbohydrate source during toxicity testing. Using this Glc-Gal-NeuriTox assay, 52 mitochondrial and non-mitochondrial toxicants were tested. The panel of chemicals comprised 11 inhibitors of mitochondrial respiratory chain complex I (cI), 4 inhibitors of cII, 8 of cIII, and 2 of cIV; 8 toxicants were included as they are assumed to be mitochondrial uncouplers. In galactose, cells became more dependent on mitochondrial function, which made them 2-3 orders of magnitude more sensitive to various mitotoxicants. Moreover, galactose enhanced the specific neurotoxicity (destruction of neurites) compared to a general cytotoxicity (plasma membrane lysis) of the toxicants. The Glc-Gal-NeuriTox assay worked particularly well for inhibitors of cI and cIII, while the toxicity of uncouplers and non-mitochondrial toxicants did not differ significantly upon glucose ↔ galactose exchange. As a secondary assay, we developed a method to quantify the inhibition of all mitochondrial respiratory chain functions/complexes in LUHMES cells. The combination of the Glc-Gal-NeuriTox neurotoxicity screening assay with the mechanistic follow up of target site identification allowed both, a more sensitive detection of neurotoxicants and a sharper definition of the mode of action of mitochondrial toxicants.

Impact of REACH legislation on the production and importation of CMR (carcinogen, mutagen and reproductive) and explosive chemicals in Italy from 2011 to 2015.

On 1 June 2007, the European Commission issued the Registration, Evaluation, Authorization and Restriction of Chemicals (REACH) to protect both the environment and human health. We analyzed the impact of REACH in the Italian market considering the presence of chemicals, their diversity, importation and production during the period 2011-2015, with particular attention to products with toxic or explosive properties. There was a reduction of the chemicals on the market, in terms of tons but also the absolute numbers of types of compounds. The production reduction was particularly noticeable for explosive chemicals: -14.7%. CMR products did not show any statistically significant reduction in term of tons: -2.3%.

QSAR Development for Plasma Protein Binding: Influence of the Ionization State.

PURPOSE: This study explored several strategies to improve the performance of literature QSAR models for plasma protein binding (PPB), such as a suitable endpoint transformation, a correct representation of chemicals, more consistency in the dataset, and a reliable definition of the applicability domain. METHODS: We retrieved human fraction unbound (Fu) data for 670 compounds from the literature and carefully checked them for consistency. Descriptors were calculated taking account of the ionization state of molecules at physiological pH (7.4), in order to better estimate the affinity of molecules to blood proteins. We used different algorithms and chemical descriptors to explore the most suitable strategy for modeling the endpoint. SMILES (simplified molecular input line entry system)-based string descriptors were also tested with the CORAL software (CORelation And Logic). We did an outlier analysis to establish the models to use (or not to use) in case of well recognized families. RESULTS: Internal validation of the selected models returned Q2 values close to 0.60. External validation also gave r2 values always greater than 0.60. The CORAL descriptor based model for √fu was the best, with r2 0.74 in external validation. CONCLUSIONS: Performance in prediction confirmed the robustness of all the derived models and their suitability for real-life purposes, i.e. screening chemicals for their ADMET profiling. Optimization of descriptors can be useful in order to obtain the correct results with a ionized molecule.

QSAR Modeling of ToxCast Assays Relevant to the Molecular Initiating Events of AOPs Leading to Hepatic Steatosis.

Nonalcoholic hepatic steatosis is a worldwide epidemiological concern since it is among the most prominent hepatic diseases. Indeed, research in toxicology and epidemiology has gathered evidence that exposure to endocrine disruptors can perturb cellular homeostasis and cause this disease. Therefore, assessing the likelihood of a chemical to trigger hepatic steatosis is a matter of the utmost importance. However, systematic in vivo testing of all the chemicals humans are exposed to is not feasible for ethical and economical reasons. In this context, predicting the molecular initiating events (MIE) leading to hepatic steatosis by QSAR modeling is an issue of practical relevance in modern toxicology. In this article, we present QSAR models based on random forest classifiers and DRAGON molecular descriptors for the prediction of in vitro assays that are relevant to MIEs leading to hepatic steatosis. These assays were provided by the ToxCast program and proved to be predictive for the detection of chemical-induced steatosis. During the modeling process, special attention was paid to chemical and toxicological data curation. We adopted two modeling strategies (undersampling and balanced random forests) to develop robust QSAR models from unbalanced data sets. The two modeling approaches gave similar results in terms of predictivity, and most of the models satisfy a minimum percentage of correctly predicted chemicals equal to 75%. Finally, and most importantly, the developed models proved to be useful as an effective in silico screening test for hepatic steatosis.

Skin sensitization quantitative QSAR models based on mechanistic structural alerts.

Allergic contact dermatitis is increasingly of interest for the hazard characterization of chemicals. in vivo animal testing is usually adopted but in silico approaches are becoming the new frontier due to their swiftness and economic efficiency. Indeed, in silico models can rationalise the experimental outcomes besides having predictive ability. The aim of the present work was to explore the electrophilic chemical behaviour responsible for allergic contact dermatitis using quantitative QSAR regression models. Eight models were proposed, using an experimental LLNA dataset of 366 chemicals. Each model is unique to encode a type of electrophilic reactivity domain. The models were obtained using autocorrelation, electro-topological and atom centered fragment based on two-dimensional descriptors, which incorporated the electronic and stereochemical features of substances interacting with skin proteins to induce skin cell proliferation. Finally, simple steps were proposed to integrate the eight models for the application on the test chemicals.

Prediction of the Partition Coefficient between Adipose Tissue and Blood for Environmental Chemicals: From Single QSAR Models to an Integrated Approach.

The adipose tissue:blood partition coefficient is a key-endpoint to predict the pharmacokinetics of chemicals in humans and animals, since other organ:blood affinities can be estimated as a function of this parameter. We performed a search in the literature to select all the available rat in vivo data. This approach resulted into two improvements to existing models: a homogeneous definition of the endpoint and an expanded data collection. The resulting dataset was used to develop QSAR models as a function of linear and non-linear algorithms. Several applicability domain definitions were assessed and the definition corresponding to a good balance between performance and coverage was retained. We assessed the pertinence of combining single models into integrated approaches to increase the accuracy in predictions. The best integrated model outperformed the single models and it was characterized by an external mean absolute error (MAE) equal to 0.26, while preserving an adequate coverage (84 %). This performance is comparable to experimental variability and it highlights the pertinence of the integrated model.

SpheraCosmolife: a new tool for the risk assessment of cosmetic products.

A new, freely available software for cosmetic products has been designed that considers the regulatory framework for cosmetics. The software allows an overall toxicological evaluation of cosmetic ingredients without the need for additional testing and, depending on the product type, it applies defined exposure scenarios to derive risk for consumers. It takes regulatory thresholds into account and uses either experimental values, if available, or predictions. Based on the exper­imental or predicted no observed adverse effect level (NOAEL), the software can define a point of departure (POD), which is used to calculate the margin of safety (MoS) of the query chemicals. The software also provides other toxico­logical properties, such as mutagenicity, skin sensitization, and the threshold of toxicological concern (TTC) to provide an overall evaluation of the potential chemical hazard. Predictions are calculated using in silico models implemented within the VEGA software. The full list of ingredients of a cosmetic product can be processed at the same time, at the effective concentrations in the product as given by the user. SpheraCosmolife is designed as a support tool for safety assessors of cosmetic products and can be used to prioritize the cosmetic ingredients or formulations according to their potential risk to consumers. The major novelty of the tool is that it wraps a series of models (some of them new) into a single, user-friendly software system.

Integrating QSAR models predicting acute contact toxicity and mode of action profiling in honey bees (A. mellifera): Data curation using open source databases, performance testing and validation.

Honey bees (Apis mellifera) provide key ecosystem services as pollinators bridging agriculture, the food chain and ecological communities, thereby ensuring food production and security. Ecological risk assessment of single Plant Protection Products (PPPs) requires an understanding of the exposure and toxicity. In silico tools such as QSAR models can play a major role for the prediction of structural, physico-chemical and pharmacokinetic properties of chemicals as well as toxicity of single and multiple chemicals. Here, the first integrative honey bee QSAR model has been developed for PPPs using EFSA's OpenFoodTox, US-EPA ECOTOX and Pesticide Properties DataBase i) to predict acute contact toxicity (LD50) and ii) to profile the Mode of Action (MoA) of pesticides active substances. Three different classification-based and four regression-based models were developed and tested for their performance, thus identifying two models providing the most reliable predictions based on k-NN algorithm. The two-category QSAR model (toxic/non-toxic; n = 411) was validated using sensitivity (=0.93), specificity (=0.85), balanced accuracy (=0.90), and Matthews correlation coefficient (MCC = 0.78) as statistical parameters. The regression-based model (n = 113) was validated for its reliability and robustness (R2 = 0.74; MAE = 0.52). Current study proposes the MoA profiling for 113 pesticides active substances and the first harmonised MoA classification scheme for acute contact toxicity in honey bees, including LD50s data points from three different databases. The classification allows to further define MoAs and the target site of PPPs active substances, thus enabling regulators and scientists to refine chemical grouping and toxicity extrapolations for single chemicals and component-based mixture risk assessment of multiple chemicals. Relevant future perspectives are briefly addressed to integrate MoA, adverse outcome pathways (AOPs) and toxicokinetic information for the refinement of single-chemical/combined toxicity predictions and risk estimates at different levels of biological organization in the bee health context.

Correction to: A new semi-automated workflow for chemical data retrieval and quality checking for modeling applications.

It was highlighted that the original article [1] contained an error in the Funding section. This Correction article states the correct and incorrect versions of the Funding section.

Ecotoxicological QSAR modeling of organic compounds against fish: Application of fragment based descriptors in feature analysis.

Organic compounds (OCs) constitute an enormously large class of highly persistent and toxic chemicals widely used for various purposes throughout the world. Their increased detection in water bodies, mainly sewage treatment plants via industrial discharge, has rendered them to become a cause for ecological concern. The limited availability of experimental toxicological data has necessitated development of models that can help us identify the most hazardous and potentially toxic compounds thus prioritizing the experiments on the selected chemicals. Computational tools such as quantitative structure-activity relationship (QSAR) can be used to predict the missing data and classify the chemicals based on their acute predicted responses for existing as well as not yet synthesized chemicals. In the current study, novel, externally validated, highly robust local QSAR models for different chemical classes and moderately robust global QSAR models were developed using partial least squares (PLS) regression technique using a large dataset of 1121 OCs for the fish mortality endpoint. For feature selection, genetic algorithm along with stepwise regression was used while model validation was performed using various stringent validation criteria following the strict rules of OECD guidelines of QSAR validation. The variables included in the models were obtained from simplex representation of molecular structures (SiRMS) (Version 4.1.2.270), Dragon (Version 7.0) and PaDEL-descriptor software (Version 2.20). The final developed models were robust, externally predictive and characterized by a large chemical as well as biological domain. The predictive efficiency of the developed models was then compared with the ECOSAR tool in order to justify the applicability of the developed models in ecotoxicological predictions for organic chemicals. Better predictive efficiency of the developed QSAR models compared to the ECOSAR derived predictions signifies their applicability in early risk assessment of known as well as untested chemicals in order to design safer alternatives to the environment.

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.

Investigating combined toxicity of binary mixtures in bees: Meta-analysis of laboratory tests, modelling, mechanistic basis and implications for risk assessment.

Bees are exposed to a wide range of multiple chemicals "chemical mixtures" from anthropogenic (e.g. plant protection products or veterinary products) or natural origin (e.g. mycotoxins, plant toxins). Quantifying the relative impact of multiple chemicals on bee health compared with other environmental stressors (e.g. varroa, viruses, and nutrition) has been identified as a priority to support the development of holistic risk assessment methods. Here, extensive literature searches and data collection of available laboratory studies on combined toxicity data for binary mixtures of pesticides and non-chemical stressors has been performed for honey bees (Apis mellifera), wild bees (Bombus spp.) and solitary bee species (Osmia spp.). From 957 screened publications, 14 publications provided 218 binary mixture toxicity data mostly for acute mortality (lethal dose: LD50) after contact exposure (61%), with fewer studies reporting chronic oral toxicity (20%) and acute oral LC50 values (19%). From the data collection, available dose response data for 92 binary mixtures were modelled using a Toxic Unit (TU) approach and the MIXTOX modelling tool to test assumptions of combined toxicity i.e. concentration addition (CA), and interactions (i.e. synergism, antagonism). The magnitude of interactions was quantified as the Model Deviation Ratio (MDR). The CA model applied to 17% of cases while synergism and antagonism were observed for 72% (MDR > 1.25) and 11% (MDR < 0.83) respectively. Most synergistic effects (55%) were observed as interactions between sterol-biosynthesis-inhibiting (SBI) fungicides and insecticide/acaricide. The mechanisms behind such synergistic effects of binary mixtures in bees are known to involve direct cytochrome P450 (CYP) inhibition, resulting in an increase in internal dose and toxicity of the binary mixture. Moreover, bees are known to have the lowest number of CYP copies and other detoxification enzymes in the insect kingdom. In the light of these findings, occurrence of these binary mixtures in relevant crops (frequency and concentrations) would need to be investigated. Addressing this exposure dimension remains critical to characterise the likelihood and plausibility of such interactions to occur under field realistic conditions. Finally, data gaps and further work for the development of risk assessment methods to assess multiple stressors in bees including chemicals and non-chemical stressors in bees are discussed.

A new semi-automated workflow for chemical data retrieval and quality checking for modeling applications.

The quality of data used for QSAR model derivation is extremely important as it strongly affects the final robustness and predictive power of the model. Ambiguous or wrong structures need to be carefully checked, because they lead to errors in calculation of descriptors, hence leading to meaningless results. The increasing amounts of data, however, have often made it hard to check of very large databases manually. In the light of this, we designed and implemented a semi-automated workflow integrating structural data retrieval from several web-based databases, automated comparison of these data, chemical structure cleaning, selection and standardization of data into a consistent, ready-to-use format that can be employed for modeling. The workflow integrates best practices for data curation that have been suggested in the recent literature. The workflow has been implemented with the freely available KNIME software and is freely available to the cheminformatics community for improvement and application to a broad range of chemical datasets.