Comparison of global and mode of action-based models for aquatic toxicity

In Silico Acute Aquatic Hazard Assessment and Prioritization Using a Grouped Target Site Model: A Case Study of Organic Substances Reported in Permian Basin Hydraulic Fracturing Operations

Hydraulic fracturing (HF) is commonly used to enhance onshore recovery of oil and gas during production. This process involves the use of a variety of chemicals to support the physical extraction of oil and gas, maintain appropriate conditions downhole (e.g., redox conditions, pH), and limit microbial growth. The diversity of chemicals used in HF presents a significant challenge for risk assessment. The objective of the present study is to establish a transparent, reproducible procedure for estimating 5th percentile acute aquatic hazard concentrations (e.g., acute hazard concentration 5th percentiles [HC5s]) for these substances and validating against existing toxicity data. A simplified, grouped target site model (gTSM) was developed using a database (n = 1696) of diverse compounds with known mode of action (MoA) information. Statistical significance testing was employed to reduce model complexity by combining 11 discrete MoAs into three general hazard groups. The new model was trained and validated using an 80:20 allocation of the experimental database. The gTSM predicts toxicity using a combination of target site water partition coefficients and hazard group-based critical target site concentrations. Model performance was comparable to the original TSM using 40% fewer parameters. Model predictions were judged to be sufficiently reliable and the gTSM was further used to prioritize a subset of reported Permian Basin HF substances for risk evaluation. The gTSM was applied to predict hazard groups, species acute toxicity, and acute HC5s for 186 organic compounds (neutral and ionic). Toxicity predictions and acute HC5 estimates were validated against measured acute toxicity data compiled for HF substances. This case study supports the gTSM as an efficient, cost-effective computational tool for rapid aquatic hazard assessment of diverse organic chemicals. Environ Toxicol Chem 2024;00:1-12. © 2024 ExxonMobil Petroleum and Chemical BV. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Application and Progress of Machine Learning in Pesticide Hazard and Risk Assessment

Long-term exposure to pesticides is associated with the incidence of cancer. With the exponential increase in the number of new pesticides being synthesized, it becomes more and more important to evaluate the toxicity of pesticides by means of simulated calculations. Based on existing data, machine learning methods can train and model the predictions of the effects of novel pesticides, which have limited available data. Combined with other technologies, this can aid the synthesis of new pesticides with specific active structures, detect pesticide residues, and identify their tolerable exposure levels. This article mainly discusses support vector machines, linear discriminant analysis, decision trees, partial least squares, and algorithms based on feedforward neural networks in machine learning. It is envisaged that this article will provide scientists and users with a better understanding of machine learning and its application prospects in pesticide toxicity assessment.

Construction of an In Silico Structural Profiling Tool Facilitating Mechanistically Grounded Classification of Aquatic Toxicants

The performance of chemical safety assessment within the domain of environmental toxicology is often impeded by a shortfall of appropriate experimental data describing potential hazards across the many compounds in regular industrial use. In silico schemes for assigning aquatic-relevant modes or mechanisms of toxic action to substances, based solely on consideration of chemical structure, have seen widespread employment─including those of Verhaar, Russom, and later Bauer (MechoA). Recently, development of a further system was reported by Sapounidou, which, in common with MechoA, seeks to ground its classifications in understanding and appreciation of molecular initiating events. Until now, this Sapounidou scheme has not seen implementation as a tool for practical screening use. Accordingly, the primary purpose of this study was to create such a resource─in the form of a computational workflow. This exercise was facilitated through the formulation of 183 structural alerts/rules describing molecular features associated with narcosis, chemical reactivity, and specific mechanisms of action. Output was subsequently compared relative to that of the three aforementioned alternative systems to identify strengths and shortcomings as regards coverage of chemical space.

Relationships Between Aquatic Toxicity, Chemical Hydrophobicity, and Mode of Action: Log Kow Revisited

Relationships between toxicity and chemical hydrophobicity have been known for nearly 100 years in mammals and fish, typically using the log of the octanol:water partition coefficient (Kow). The current study reassessed the influence of mode of action (MOA) on acute aquatic toxicity-log Kow relationships using a comprehensive database of 617 organic chemicals with curated and standardized acute toxicity data that did not exceed solubility limits, their consensus log Kow values, and weight of evidence-based MOA classifications (including 6 broad and 26 specific MOAs). A total of 166 significant (p < 0.05) log Kow-toxicity models were developed across six taxa groups that included QSARs for 5 of the broad and 13 of the specific MOAs. In this study, we demonstrate that QSARs based on MOAs can significantly increase LC50 prediction accuracy for specific acting chemicals. Prediction accuracy increases when QSARs are built based on highly specific MOAs, rather than broad MOA classifications. Additionally, we demonstrate that building QSAR models with chemicals in specific MOA groupings, rather than broader MOA groups leads to significantly better estimates. We also evaluated the differences between models developed from mass-based (µg/L) and mole-based (µmol/L) toxicity data and demonstrate that both are suitable for QSAR development with no clear trend in greater model accuracy. Overall, the results reveal that, despite high variance in all taxa and MOA groups, specific MOA-based models can improve the accuracy of aquatic toxicity predictions over more general groupings.Please check and confirm that the authors and their respective affiliations have been correctly identified and amend if necessary.The affiliations are correct.

Prior Knowledge for Predictive Modeling: The Case of Acute Aquatic Toxicity

Early assessment of the potential impact of chemicals on health and the environment requires toxicological properties of the molecules. Predictive modeling is often used to estimate the property values in silico from pre-existing experimental data, which is often scarce and uncertain. One of the ways to advance the predictive modeling procedure might be the use of knowledge existing in the field. Scientific publications contain a vast amount of knowledge. However, the amount of manual work required to process the enormous volumes of information gathered in scientific articles might hinder its utilization. This work explores the opportunity of semiautomated knowledge extraction from scientific papers and investigates a few potential ways of its use for predictive modeling. The knowledge extraction and predictive modeling are applied to the field of acute aquatic toxicity. Acute aquatic toxicity is an important parameter of the safety assessment of chemicals. The extensive amount of diverse information existing in the field makes acute aquatic toxicity an attractive area for investigation of knowledge use for predictive modeling. The work demonstrates that the knowledge collection and classification procedure could be useful in hybrid modeling studies concerning the model and predictor selection, addressing data gaps, and evaluation of models' performance.

Humic acid promoted activation of peroxymonosulfate by Fe3S4 for degradation of 2,4,6-trichlorophenol: An experimental and theoretical study

Chlorophenols are difficult to degrade and biohazardous in the natural environment. This study demonstrated that humic acid (HA) could promote Fe₃S₄ activation of peroxymonosulfate (PMS) to degrade 2,4,6-trichlorophenol (TCP), the degradation efficiency of TCP was increased by 33%. The system of Fe₃S₄-HA/PMS produced more reactive oxygen species, and •OH was the dominant ROS. The genealogy of iron oxides together with S⁰ on the Fe₃S₄ surface inhibited PMS activation leading to the significant reduction of TCP degraded (< 70%). These problems could be solved successfully through introducing HA, which facilitated electron transfer and increased the continuous release of iron ions by 2 times. In accordance with the determined density functional theory (DFT), the degradation pathway was put forward, which indicated that TCP dechlorination and oxidation to 2,6-dichloro-1,4-benzoquinone constituted the main degradation pathway. Furthermore, the intermediates that were produced in the main degradation processes of TCP showed lower toxicity than TCP according to results that were obtained utilizing the calculations of quantitative structure-activity relationship (QSAR) together with Toxicity Estimation Software Tool (TEST). Thus, the Fe₃S₄-HA/PMS system was demonstrated to be an efficient and safe technology for organic pollutant degradation in contaminated groundwater and surface water environments.

Implementation of in silico toxicology protocols within a visual and interactive hazard assessment platform

Mechanistically-driven alternative approaches to hazard assessment invariably require a battery of tests, including both in silico models and experimental data. The decision-making process, from selection of the methods to combining the information based on the weight-of-evidence, is ideally described in published guidelines or protocols. This ensures that the application of such approaches is defendable to reviewers within regulatory agencies and across the industry. Examples include the ICH M7 pharmaceutical impurities guideline and the published in silico toxicology protocols. To support an efficient, transparent, consistent and fully documented implementation of these protocols, a new and novel interactive software solution is described to perform such an integrated hazard assessment based on public and proprietary information.

Implementation of In Silico Toxicology Protocols in Leadscope

In silico toxicology protocols help ensure the results from computational toxicology models are performed and documented in a standardized, consistent, transparent, and accepted manner. In silico toxicology protocols for skin sensitization and genetic toxicology have been previously published and have been implemented within the Leadscope software. The following chapter outlines how such protocols have been deployed in the Leadscope platform including integration with toxicology databases and a battery of computational toxicology models.

Detection and toxicity modeling of anthraquinone dyes and chlorinated side products from a colored smoke pyrotechnic reaction

"Green" pyrotechnics seek to remove known environmental pollutants and health hazards from their formulations. This chemical engineering approach often focuses on maintaining performance effects upon replacement of objectionable ingredients, yet neglects the chemical products formed by the exothermic reaction. In this work, milligram quantities of a lab-scale pyrotechnic red smoke composition were functioned within a thermal probe for product identification by pyrolysis-gas chromatography-mass spectrometry. Thermally decomposed ingredients and new side product derivatives were identified at lower relative abundances to the intact organic dye (as the engineered sublimation product). Side products included chlorination of the organic dye donated by the chlorate oxidizer. Machine learning quantitative structure-activity relationship models computed impacts to health and environmental hazards. High to very high toxicities were predicted for inhalation, mutagenicity, developmental, and endocrine disruption for common military pyrotechnic dyes and their analogous chlorinated side products. These results underscore the need to revise objectives of "green" pyrotechnic engineering.

Generating accurate in silico predictions of acute aquatic toxicity for a range of organic chemicals: Towards similarity-based machine learning methods

There has been an increase in the use of non-animal approaches, such as in silico and/or in vitro methods, for assessing the risks of hazardous chemicals. A number of machine learning algorithms link molecular descriptors that interpret chemical structural properties with their biological activity. These computer-aided methods encounter several challenges, the most significant being the heterogeneity of datasets; more efficient and inclusive computational methods that are able to process large and heterogeneous chemical datasets are needed. In this context, this study verifies the utility of similarity-based machine learning methods in predicting the acute aquatic toxicity of diverse organic chemicals on Daphnia magna and Oryzias latipes. Two similarity-based methods were tested that employ a limited training dataset, most similar to a given fitting point, instead of using the entire dataset that encompasses a wide range of chemicals. The kernel-weighted local polynomial approach had a number of advantages over the distance-weighted k-nearest neighbor (k-NN) algorithm. The results highlight the importance of lipophilicity, electrophilic reactivity, molecular polarizability, and size in determining acute toxicity. The rigorous model validation ensures that this approach is an important tool for estimating toxicity in new or untested chemicals.

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 (LD₅₀) 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 (R² = 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 LD_50s 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.

MaNGA: a novel multi-niche multi-objective genetic algorithm for QSAR modelling

SUMMARY

Quantitative structure-activity relationship (QSAR) modelling is currently used in multiple fields to relate structural properties of compounds to their biological activities. This technique is also used for drug design purposes with the aim of predicting parameters that determine drug behaviour. To this end, a sophisticated process, involving various analytical steps concatenated in series, is employed to identify and fine-tune the optimal set of predictors from a large dataset of molecular descriptors (MDs). The search of the optimal model requires to optimize multiple objectives at the same time, as the aim is to obtain the minimal set of features that maximizes the goodness of fit and the applicability domain (AD). Hence, a multi-objective optimization strategy, improving multiple parameters in parallel, can be applied. Here we propose a new multi-niche multi-objective genetic algorithm that simultaneously enables stable feature selection as well as obtaining robust and validated regression models with maximized AD. We benchmarked our method on two simulated datasets. Moreover, we analyzed an aquatic acute toxicity dataset and compared the performances of single- and multi-objective fitness functions on different regression models. Our results show that our multi-objective algorithm is a valid alternative to classical QSAR modelling strategy, for continuous response values, since it automatically finds the model with the best compromise between statistical robustness, predictive performance, widest AD, and the smallest number of MDs.

AVAILABILITY AND IMPLEMENTATION

The python implementation of MaNGA is available at https://github.com/Greco-Lab/MaNGA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mesoporous Silica Platforms with Potential Applications in Release and Adsorption of Active Agents

In recent years, researchers focused their attention on mesoporous silica nanoparticles (MSNs) owing to the considerable advancements of the characterization methods, especially electron microscopy methods, which allowed for a clear visualization of the pore structure and the materials encapsulated within the pores, along with the X-ray diffraction (small angles) methods and specific surface area determination by Brunauer–Emmett–Teller (BET) technique. Mesoporous silica gained important consideration in biomedical applications thanks to its tunable pore size, high surface area, surface functionalization possibility, chemical stability, and pore nature. Specifically, the nature of the pores allows for the encapsulation and release of anti-cancer drugs into tumor tissues, which makes MSN ideal candidates as drug delivery carriers in cancer treatment. Moreover, the inner and outer surfaces of the MSN provide a platform for further functionalization approaches that could enhance the adsorption of the drug within the silica network and the selective targeting and controlled release to the desired site. Additionally, stimuli-responsive mesoporous silica systems are being used as mediators in cancer therapy, and through the release of the therapeutic agents hosted inside the pores under the action of specific triggering factors, it can selectively deliver them into tumor tissues. Another important application of the mesoporous silica nanomaterials is related to its ability to extract different hazardous species from aqueous media, some of these agents being antibiotics, pesticides, or anti-tumor agents. The purpose of this paper is to analyze the methods of MSN synthesis and related characteristics, the available surface functionalization strategies, and the most important applications of MSN in adsorption as well as release studies. Owing to the increasing antibiotic resistance, the need for developing materials for antibiotic removal from wastewaters is important and mesoporous materials already proved remarkable performances in environmental applications, including removal or even degradation of hazardous agents such as antibiotics and pesticides.

Mode of Action Classifications in the EnviroTox Database: Development and Implementation of a Consensus MOA Classification

Multiple mode of action (MOA) frameworks have been developed in aquatic ecotoxicology, mainly based on fish toxicity. These frameworks provide information on a key determinant of chemical toxicity, but the MOA categories and level of specificity remain unique to each of the classification schemes. The present study aimed to develop a consensus MOA assignment within EnviroTox, a curated in vivo aquatic toxicity database, based on the following MOA classification schemes: Verhaar (modified) framework, Assessment Tool for Evaluating Risk, Toxicity Estimation Software Tool, and OASIS. The MOA classifications from each scheme were first collapsed into one of 3 categories: non-specifically acting (i.e., narcosis), specifically acting, or nonclassifiable. Consensus rules were developed based on the degree of concordance among the 4 individual MOA classifications to attribute a consensus MOA to each chemical. A confidence rank was also assigned to the consensus MOA classification based on the degree of consensus. Overall, 40% of the chemicals were classified as narcotics, 17% as specifically acting, and 43% as unclassified. Sixty percent of chemicals had a medium to high consensus MOA assignment. When compared to empirical acute toxicity data, the general trend of specifically acting chemicals being more toxic is clearly observed for both fish and invertebrates but not for algae. EnviroTox is the first approach to establishing a high-level consensus across 4 computationally and structurally distinct MOA classification schemes. This consensus MOA classification provides both a transparent understanding of the variation between MOA classification schemes and an added certainty of the MOA assignment. In terms of regulatory relevance, a reliable understanding of MOA can provide information that can be useful for the prioritization (ranking) and risk assessment of chemicals. Environ Toxicol Chem 2019;38:2294-2304. © 2019 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals, Inc. on behalf of SETAC.

Chemometric modeling of Daphnia magna toxicity of agrochemicals

Over the past few years, the ecotoxicological hazard potential of agrochemicals has received much attention in the industries and regulatory agencies. In the current work, we have developed quantitative structure-activity relationship (QSAR) models for Daphnia magna toxicities of different classes of agrochemicals (fungicides, herbicides, insecticides and microbiocides) individually as well as for the combined set with the application of Organization for Economic Co-operation and Development (OECD) recommended guidelines. The models for the individual data sets as well as for the combined set were generated employing only simple and interpretable two-dimensional descriptors, and subsequently strictly validated using test set compounds. The validated individual models were used to generate consensus models, with the objective to improve the prediction quality and reduced prediction errors. All the individual models of different classes of agrochemicals as well as the global set of agrochemicals showed encouraging statistical quality and prediction ability. The general observations from the derived models suggest that the toxicity increases with lipophilicity and decreases with polarity. The generated models of different classes of agrochemicals and also for the combined set should be applicable for data gap filling for new or untested agrochemical compounds falling within the applicability domain of the developed models.

Modeling the Sensitivity of Aquatic Macroinvertebrates to Chemicals Using Traits

In this study, a trait-based macroinvertebrate sensitivity modeling tool is presented that provides two main outcomes: (1) it constructs a macroinvertebrate sensitivity ranking and, subsequently, a predictive trait model for each one of a diverse set of predefined Modes of Action (MOAs) and (2) it reveals data gaps and restrictions, helping with the direction of future research. Besides revealing taxonomic patterns of species sensitivity, we find that there was not one genus, family, or class which was most sensitive to all MOAs and that common test taxa were often not the most sensitive at all. Traits like life cycle duration and feeding mode were identified as important in explaining species sensitivity. For 71% of the species, no or incomplete trait data were available, making the lack of trait data the main obstacle in model construction. Research focus should therefore be on completing trait databases and enhancing them with finer morphological traits, focusing on the toxicodynamics of the chemical (e.g., target site distribution). Further improved sensitivity models can help with the creation of ecological scenarios by predicting the sensitivity of untested species. Through this development, our approach can help reduce animal testing and contribute toward a new predictive ecotoxicology framework.

Target site model: Predicting mode of action and aquatic organism acute toxicity using Abraham parameters and feature-weighted k-nearest neighbors classification

A database of 1480 chemicals with 47 associated modes of action compiled from the literature encompasses a wide range of chemical classes (alkanes, polycyclic aromatic hydrocarbons, pesticides, and polar compounds) and includes toxicity data for 79 different aquatic genera. The data were split into a calibration group and a validation group (80/20) to apply k-nearest neighbors (k-NN) methodology to predict the toxic mode of action for the compound. Other approaches were tested (support vector machines and linear discriminant analysis) as well as variations in the k-NN technique (distance weighting, feature weighting). Best-prediction results were found with k = 3, in a voting platform with optimized feature weighting. Using the predicted mode of action, the appropriate polyparameter target site model for that mode of action is applied to calculate the 50% lethal concentration (LC50). Predicted LC50s for the validation database resulted in a root-mean squared error (RMSE) of 0.752. This can be compared to an RMSE of 0.655 for the same validation set using the reference mode of action labels. The complete database resulted in an RMSE of 0.793 for reference mode of action labels. This confirms that the classification model has sufficient accuracy for predicting the mode of action and for determining toxicity using the target site model. Environ Toxicol Chem 2019;38:375-386. © 2018 SETAC.

Target site model: Application of the polyparameter target lipid model to predict aquatic organism acute toxicity for various modes of action

A database of 2049 chemicals with 47 associated modes of action (MoA) was compiled from the literature. The database includes alkanes, polycyclic aromatic hydrocarbons, pesticides, inorganic, and polar compounds. Brief descriptions of some critical MoA classification groups are provided. The MoA from the 14 sources were assigned using a variety of reliable experimental and modeling techniques. Toxicity information, chemical parameters, and solubility limits were combined with the MoA label information to create the data set used for model development. The model database was used to generate linear free energy relationships for each specific MoA using multilinear regression analysis. The model uses chemical-specific Abraham solute parameters estimated from AbSolv to determine MoA-specific solvent parameters. With this procedure, critical target site concentrations are determined for each genus. Statistical analysis showed a wide range in values of the solvent parameters for the significant MoA. Environ Toxicol Chem 2019;38:222-239. © 2018 SETAC.

In silico toxicology protocols

The present publication surveys several applications of in silico (i.e., computational) toxicology approaches across different industries and institutions. It highlights the need to develop standardized protocols when conducting toxicity-related predictions. This contribution articulates the information needed for protocols to support in silico predictions for major toxicological endpoints of concern (e.g., genetic toxicity, carcinogenicity, acute toxicity, reproductive toxicity, developmental toxicity) across several industries and regulatory bodies. Such novel in silico toxicology (IST) protocols, when fully developed and implemented, will ensure in silico toxicological assessments are performed and evaluated in a consistent, reproducible, and well-documented manner across industries and regulatory bodies to support wider uptake and acceptance of the approaches. The development of IST protocols is an initiative developed through a collaboration among an international consortium to reflect the state-of-the-art in in silico toxicology for hazard identification and characterization. A general outline for describing the development of such protocols is included and it is based on in silico predictions and/or available experimental data for a defined series of relevant toxicological effects or mechanisms. The publication presents a novel approach for determining the reliability of in silico predictions alongside experimental data. In addition, we discuss how to determine the level of confidence in the assessment based on the relevance and reliability of the information.

Expert judgment based multicriteria decision models to assess the risk of pesticides on reproduction failures of grey partridge

A suite of models is proposed for estimating the risk of pesticides against the grey partridge (Perdix perdix) and their clutches. Radio-tracked data of females, description and location of the clutches, and data on the pesticide treatments during the laying periods of the partridges were used as basic information. Quantitative structure-activity relationship (QSAR) and quantitative structure-property relationship (QSPR) modelling allowed us to characterize the pesticides by their 1-octanol/water partition coefficient (log P), vapour pressure, primary and ultimate biodegradation potential, acute toxicity (LD₅₀) on P. perdix, and endocrine disruption potential. From these physicochemical and toxicological data, the system of integration of risk with interaction of scores (SIRIS) method was used to design scores of risk for pesticides, alone or in mixture. A program, written in R (version 3.1.1), called Simulation of Toxicity in Perdix perdix (SimToxPP), was designed for estimating the risk of substances, considered alone or in mixture, against the grey partridge during breeding. The software tool is flexible enough to simulate realistic in situ scenarios. Different examples of applications are shown. The advantages and limitations of the approach are briefly discussed.

Global concentration additivity and prediction of mixture toxicities, taking nitrobenzene derivatives as an example

The toxicity of a mixture depends not only on the mixture concentration level but also on the mixture ratio. For a multiple-component mixture (MCM) system with a definite chemical composition, the mixture toxicity can be predicted only if the global concentration additivity (GCA) is validated. The so-called GCA means that the toxicity of any mixture in the MCM system is the concentration additive, regardless of what its mixture ratio and concentration level. However, many mixture toxicity reports have usually employed one mixture ratio (such as the EC₅₀ ratio), the equivalent effect concentration ratio (EECR) design, to specify several mixtures. EECR mixtures cannot simulate the concentration diversity and mixture ratio diversity of mixtures in the real environment, and it is impossible to validate the GCA. Therefore, in this paper, the uniform design ray (UD-Ray) was used to select nine mixture ratios (rays) in the mixture system of five nitrobenzene derivatives (NBDs). The representative UD-Ray mixtures can effectively and rationally describe the diversity in the NBD mixture system. The toxicities of the mixtures to Vibrio qinghaiensis sp.-Q67 were determined by the microplate toxicity analysis (MTA). For each UD-Ray mixture, the concentration addition (CA) model was used to validate whether the mixture toxicity is additive. All of the UD-Ray mixtures of five NBDs are global concentration additive. Afterwards, the CA is employed to predict the toxicities of the external mixtures from three EECR mixture rays with the NOEC, EC₃₀, and EC₇₀ ratios. The predictive toxicities are in good agreement with the experimental toxicities, which testifies to the predictability of the mixture toxicity of the NBDs.

Mode of Action (MOA) Assignment Classifications for Ecotoxicology: An Evaluation of Approaches

The mode of toxic action (MOA) is recognized as a key determinant of chemical toxicity and as an alternative to chemical class-based predictive toxicity modeling. However, MOA classification has never been standardized in ecotoxicology, and a comprehensive comparison of classification tools and approaches has never been reported. Here we critically evaluate three MOA classification methodologies using an aquatic toxicity data set of 3448 chemicals, compare the approaches, and assess utility and limitations in screening and early tier assessments. The comparisons focused on three commonly used tools: Verhaar prediction of toxicity MOA, the U.S. Environmental Protection Agency (EPA) ASsessment Tool for Evaluating Risk (ASTER) QSAR (quantitative structure activity relationship) application, and the EPA Mode of Action and Toxicity (MOAtox) database. Of the 3448 MOAs predicted using the Verhaar scheme, 1165 were classified by ASTER, and 802 were available in MOAtox. Of the subset of 432 chemicals with MOA assignments for each of the three schemes, 42% had complete concordance in MOA classification, and there was no agreement for 7% of the chemicals. The research shows the potential for large differences in MOA classification between the five broad groups of the Verhaar scheme and the more mechanism-based assignments of ASTER and MOAtox. Harmonization of classification schemes is needed to use MOA classification in chemical hazard and risk assessment more broadly.

Prediction of pesticide acute toxicity using two-dimensional chemical descriptors and target species classification

Previous modelling of the median lethal dose (oral rat LD₅₀) has indicated that local class-based models yield better correlations than global models. We evaluated the hypothesis that dividing the dataset by pesticidal mechanisms would improve prediction accuracy. A linear discriminant analysis (LDA) based-approach was utilized to assign indicators such as the pesticide target species, mode of action, or target species - mode of action combination. LDA models were able to predict these indicators with about 87% accuracy. Toxicity is predicted utilizing the QSAR model fit to chemicals with that indicator. Toxicity was also predicted using a global hierarchical clustering (HC) approach which divides data set into clusters based on molecular similarity. At a comparable prediction coverage (~94%), the global HC method yielded slightly higher prediction accuracy (r² = 0.50) than the LDA method (r² ~ 0.47). A single model fit to the entire training set yielded the poorest results (r² = 0.38), indicating that there is an advantage to clustering the dataset to predict acute toxicity. Finally, this study shows that whilst dividing the training set into subsets (i.e. clusters) improves prediction accuracy, it may not matter which method (expert based or purely machine learning) is used to divide the dataset into subsets.

In silico prediction of toxicity of phenols to Tetrahymena pyriformis by using genetic algorithm and decision tree-based modeling approach

Risk assessment of chemicals is an important issue in environmental protection; however, there is a huge lack of experimental data for a large number of end-points. The experimental determination of toxicity of chemicals involves high costs and time-consuming process. In silico tools such as quantitative structure-toxicity relationship (QSTR) models, which are constructed on the basis of computational molecular descriptors, can predict missing data for toxic end-points for existing or even not yet synthesized chemicals. Phenol derivatives are known to be aquatic pollutants. With this background, we aimed to develop an accurate and reliable QSTR model for the prediction of toxicity of 206 phenols to Tetrahymena pyriformis. A multiple linear regression (MLR)-based QSTR was obtained using a powerful descriptor selection tool named Memorized_ACO algorithm. Statistical parameters of the model were 0.72 and 0.68 for R_training² and R_test², respectively. To develop a high-quality QSTR model, classification and regression tree (CART) was employed. Two approaches were considered: (1) phenols were classified into different modes of action using CART and (2) the phenols in the training set were partitioned to several subsets by a tree in such a manner that in each subset, a high-quality MLR could be developed. For the first approach, the statistical parameters of the resultant QSTR model were improved to 0.83 and 0.75 for R_training² and R_test², respectively. Genetic algorithm was employed in the second approach to obtain an optimal tree, and it was shown that the final QSTR model provided excellent prediction accuracy for the training and test sets (R_training² and R_test² were 0.91 and 0.93, respectively). The mean absolute error for the test set was computed as 0.1615.

A Bayesian network model for predicting aquatic toxicity mode of action using two dimensional theoretical molecular descriptors

The mode of toxic action (MoA) has been recognized as a key determinant of chemical toxicity, but development of predictive MoA classification models in aquatic toxicology has been limited. We developed a Bayesian network model to classify aquatic toxicity MoA using a recently published dataset containing over one thousand chemicals with MoA assignments for aquatic animal toxicity. Two dimensional theoretical chemical descriptors were generated for each chemical using the Toxicity Estimation Software Tool. The model was developed through augmented Markov blanket discovery from the dataset of 1098 chemicals with the MoA broad classifications as a target node. From cross validation, the overall precision for the model was 80.2%. The best precision was for the AChEI MoA (93.5%) where 257 chemicals out of 275 were correctly classified. Model precision was poorest for the reactivity MoA (48.5%) where 48 out of 99 reactive chemicals were correctly classified. Narcosis represented the largest class within the MoA dataset and had a precision and reliability of 80.0%, reflecting the global precision across all of the MoAs. False negatives for narcosis most often fell into electron transport inhibition, neurotoxicity or reactivity MoAs. False negatives for all other MoAs were most often narcosis. A probabilistic sensitivity analysis was undertaken for each MoA to examine the sensitivity to individual and multiple descriptor findings. The results show that the Markov blanket of a structurally complex dataset can simplify analysis and interpretation by identifying a subset of the key chemical descriptors associated with broad aquatic toxicity MoAs, and by providing a computational chemistry-based network classification model with reasonable prediction accuracy.

Adverse Outcome Pathway (AOP) Informed Modeling of Aquatic Toxicology: QSARs, Read-Across, and Interspecies Verification of Modes of Action

Alternative approaches have been promoted to reduce the number of vertebrate and invertebrate animals required for the assessment of the potential of compounds to cause harm to the aquatic environment. A key philosophy in the development of alternatives is a greater understanding of the relevant adverse outcome pathway (AOP). One alternative method is the fish embryo toxicity (FET) assay. Although the trends in potency have been shown to be equivalent in embryo and adult assays, a detailed mechanistic analysis of the toxicity data has yet to be performed; such analysis is vital for a full understanding of the AOP. The research presented herein used an updated implementation of the Verhaar scheme to categorize compounds into AOP-informed categories. These were then used in mechanistic (quantitative) structure-activity relationship ((Q)SAR) analysis to show that the descriptors governing the distinct mechanisms of acute fish toxicity are capable of modeling data from the FET assay. The results show that compounds do appear to exhibit the same mechanisms of toxicity across life stages. Thus, this mechanistic analysis supports the argument that the FET assay is a suitable alternative testing strategy for the specified mechanisms and that understanding the AOPs is useful for toxicity prediction across test systems.

QSAR study of the acute toxicity to fathead minnow based on a large dataset

Acute fathead minnow toxicity is an important basis of hazard and risk assessment for compounds in the aquatic environment. In this paper, a large dataset consisting of 963 organic compounds with acute toxicity towards fathead minnow was studied with a QSAR approach. All molecular structures of compounds were optimized by the hybrid density functional theory method. Dragon molecular descriptors and log Kow were selected to describe molecular information. Genetic algorithm and multiple linear regression analysis were combined to develop models. A global prediction model for compounds without known mode of action and two local models for organic compounds that exhibit narcosis toxicity and excess toxicity were developed, respectively. For all developed models, internal validations were performed by cross-validation and external validations were implemented by the setting of validation set. In addition, applicability domains of models were evaluated using a leverage method and outliers were listed and checked using toxicological knowledge.

Investigation of the Verhaar scheme for predicting acute aquatic toxicity: improving predictions obtained from Toxtree ver. 2.6

Assessment of the potential of compounds to cause harm to the aquatic environment is an integral part of the REACH legislation. To reduce the number of vertebrate and invertebrate animals required for this analysis alternative approaches have been promoted. Category formation and read-across have been applied widely to predict toxicity. A key approach to grouping for environmental toxicity is the Verhaar scheme which uses rules to classify compounds into one of four mechanistic categories. These categories provide a mechanistic basis for grouping and any further predictive modelling. A computational implementation of the Verhaar scheme is available in Toxtree v2.6. The work presented herein demonstrates how modifications to the implementation of Verhaar between version 1.5 and 2.6 of Toxtree have improved performance by reducing the number of incorrectly classified compounds. However, for the datasets used in this analysis, version 2.6 classifies more compounds as outside of the domain of the model. Further amendments to the classification rules have been implemented here using a post-processing filter encoded as a KNIME workflow. This results in fewer compounds being classified as outside of the model domain, further improving the predictivity of the scheme. The utility of the modification described herein is demonstrated through building quality, mechanism-specific Quantitative Structure Activity Relationship (QSAR) models for the compounds within specific mechanistic categories.