QSAR modeling of Daphnia magna and fish toxicities of biocides using 2D descriptors

Safer and greener chemicals for the aquatic ecosystem: Chemometric modeling of the prolonged and chronic aquatic toxicity of chemicals on Oryzias latipes

In the modern era, chemicals and their products have been used everywhere like agriculture, healthcare, food, cosmetics, pharmaceuticals, household products, clothing industry, etc. These chemicals find their way to reach the aquatic ecosystem (directly/indirectly) and cause severe chronic and prolonged toxic effects to aquatic species which is also then translated to human beings. Prolonged and chronic toxicity data of many chemicals that are used daily is not available due to high experimentation testing costs, time investment, and the requirement of a large number of animal sacrifices. Thus, in silico approaches (e.g., QSAR (quantitative structure-activity relationship)) are the best alternative for chronic and prolonged toxicity predictions. The present work offers multi-endpoint (five endpoints: chronic_LOEC, prolonged_14D_LC50, prolonged_14D_NOEC, prolonged_21D_LC50, prolonged_21D_NOEC) QSAR models for addressing the prolonged and chronic aquatic toxicity of chemicals toward fish (O. latipes). The statistical results (R2 =0.738-0.869, QLOO2 =0.712-0.831, Q(F1)2 =0.618-0.731) of the developed models show that they were robust, reliable, reproducible, accurate, and predictive. Some of the features that are responsible for prolonged and chronic toxicity of chemicals towards O. latipes are as follows: the presence of substituted benzene, hydrophobicity, unsaturation, electronegativity, the presence of long-chain fragments, the presence of a greater number of atoms at conjugation, and the presence of halogen atoms. On the other hand, hydrophilicity and graph density descriptors retard the aquatic chronic and prolonged toxicity of chemicals toward O. latipes. The PPDB (pesticide properties database) and experimental and investigational classes of drugs from the DrugBank database were also screened using the developed model. Thus, these multi-endpoint models will be helpful for data-gap filling and provide a broad range of applicability. Therefore, this research will aid in the in silico QSAR (quantitative structure-activity relationship) prediction (non-animal testing) of the prolonged and chronic toxicity of untested and new toxic chemicals/drugs/pesticides, design and development of eco-friendly, novel, and safer chemicals, and help to protect the aquatic ecosystem from exposure to toxic and hazardous chemicals.

The VERA software: Implementation of the acute fish toxicity endpoint and its application to pharmaceutical compounds

The Virtual Extensive Read-Across software (VERA) is a new tool for read-across using a global similarity score, molecular groups, and structural alerts to find clusters of similar substances; these clusters are then used to identify suitable similar substances and make an assessment for the target substance. A beta version of VERA GUI is free and available at vegahub.eu; the source code of the VERA algorithm is available on GitHub. In the past we described its use to assess carcinogenicity, a classification endpoint. The aim here is to extend the automated read-across approach to assess continuous endpoints as well. We addressed acute fish toxicity. VERA evaluation on the acute fish toxicity endpoint was done on a dataset containing general substances (pesticides, industrial products, biocides, etc.), obtaining an overall R2 of 0.68. We employed the VERA algorithm also on active pharmaceutical ingredients (APIs). We included a portion of the APIs in the training dataset to predict APIs, successfully achieving an overall R2 of 0.63. VERA evaluates the assessment's reliability, and we reached an R2 of 0.78 and Root Mean Square Error (RMSE) of 0.44 for predictions with high reliability.

Quantitative structure-activity relationship predicting toxicity of pesticides towards Daphnia magna

Global pesticide usage reaching 2.7 million metric tons annually, brings a grave threat to non-target organisms, especially aquatic organisms, resulting in serious concerns. Predicting aquatic toxicity of pesticides towards Daphnia magna is significant. In this work, random forest (RF) algorithm, together with ten Dragon molecular descriptors, was successfully utilized to develop a quantitative structure-activity/toxicity relationship (QSAR/QSTR) model for the toxicity pEC50 of 745 pesticides towards Daphnia magna. The optimal QSTR model (RF Model I) based on the RF parameters of ntree = 50, mtry = 3 and nodesize = 5, yielded R2 = 0.877, MAE = 0.570, rms = 0.739 (training set of 596 pEC50), R2 = 0.807, MAE = 0.732, rms = 0.902 (test set of 149 pEC50), and R2 = 0.863, MAE = 0.602, rms = 0.774 (total set of 745 pEC50), which are accurate and satisfactory. The optimal RF model is comparable to other published QSTR models for Daphnia magna, although the optimal RF model possessed a small descriptor subset and dealt with a large dataset of pesticide toxicity pEC50. Thus, the investigation in this work provides a reliable, applicable QSTR model for predicting the toxicity pEC50 of pesticides towards Daphnia magna.

First report on chemometrics-driven multilayered lead prioritization in addressing oxysterol-mediated overexpression of G protein-coupled receptor 183

Contemporary research has convincingly demonstrated that upregulation of G protein-coupled receptor 183 (GPR183), orchestrated by its endogenous agonist, 7α,25-dihydroxyxcholesterol (7α,25-OHC), leads to the development of cancer, diabetes, multiple sclerosis, infectious, and inflammatory diseases. A recent study unveiled the cryo-EM structure of 7α,25-OHC bound GPR183 complex, presenting an untapped opportunity for computational exploration of potential GPR183 inhibitors, which served as our inspiration for the current work. A predictive and validated two-dimensional QSAR model using genetic algorithm (GA) and multiple linear regression (MLR) on experimental GPR183 inhibition data was developed. QSAR study highlighted that structural features like dissimilar electronegative atoms, quaternary carbon atoms, and CH2RX fragment (X: heteroatoms) influence positively, while the existence of oxygen atoms with a topological separation of 3, negatively affects GPR183 inhibitory activity. Post assessment of true external set prediction capability, the MLR model was deployed to screen 12,449 DrugBank compounds, followed by a screening pipeline involving molecular docking, druglikeness, ADMET, protein-ligand stability assessment using deep learning algorithm, molecular dynamics, and molecular mechanics. The current findings strongly evidenced DB05790 as a potential lead for prospective interference of oxysterol-mediated GPR183 overexpression, warranting further in vitro and in vivo validation.

First report on pesticide sub-chronic and chronic toxicities against dogs using QSAR and chemical read-across

Excessive use of chemicals is the outcome of the industrialization of agricultural sectors which leads to disturbance of ecological balance. Various agrochemicals are widely used in agricultural fields, urban green areas, and to protect from various pest-associated diseases. Due to their long-term health and environmental hazards, chronic toxicity assessment is crucial. Since in vivo and in vitro toxicity assessments are costly, lengthy, and require a large number of animal experiments, in silico toxicity approaches are better alternatives to save time, cost, and animal experimentation. We have developed the first regression-based 2D-QSAR models using different sub-chronic and chronic toxicity data of pesticides against dogs employing 2D descriptors. From the statistical results (n train = 53 - 62 ,   r 2  = 0.614 to 0.754, Q L O O 2  = 0.501 to 0.703 and   Q F 1 2  = 0.531 to 0.718, Q F 2 2 = 0.523 - 0.713 ), it was concluded that the models are robust, reliable, interpretable, and predictive. Similarity-based read-across algorithm was also used to improve the predictivity (Q F 1 2 = 0.595 - 0.813 , Q F 2 2 = 0.573 - 0.809 ) of the models. 5132 chemicals obtained from the CPDat and 1694 pesticides obtained from the PPDB database were also screened using the developed models, and their predictivity and reliability were checked. Thus, these models will be helpful for eco-toxicological data-gap filling, toxicity prediction of untested pesticides, and development of novel, safer & eco-friendly pesticides.

Unveiling first report on in silico modeling of aquatic toxicity of organic chemicals to Labeo rohita (Rohu) employing QSAR and q-RASAR

Labeo rohita, a fish species within the Carp family, holds significant dietary and aquacultural importance in South Asian countries. However, the habitats of L. rohita often face exposure to various harmful pesticides and organic compounds originating from industrial and agricultural runoff. It is challenging to individually investigate the effects of each potentially harmful compound. In such cases, in silico techniques like Quantitative Structure-Activity Relationship (QSAR) and quantitative Read-Across Structure-Activity Relationship (q-RASAR) can be employed to construct algorithmic models capable of simultaneously assessing the toxicity of numerous compounds. We utilized the US EPA's ToxValDB database to curate data regarding acute median lethal concentration (LC50) toxicity for L. rohita. The experimental variables included study type (mortality), study duration (ranging from 0.25 h to 4 h), exposure route (static, flowthrough, and renewal), exposure method (drinking water), and types of chemicals (industrial chemicals and pharmaceuticals). Using this dataset, we developed regression-based QSAR and q-RASAR models to predict chemical toxicity to L. rohita based on chemical descriptors. The key descriptors for predicting the toxicity of L. rohita in the regression-based QSAR model include F05[S-Cl], SpMax_EA(ri), s4_relPathLength_2, and SpDiam_AEA(ed). These descriptors can be employed to estimate the toxicity of untested compounds and aid in the development of compounds with lower toxicity based on the presence or absence of these descriptors. Both the QSAR and q-RASAR models serve as valuable tools for understanding the chemicals' structural features responsible for toxicity and for filling gaps in aquatic toxicity data by predicting the toxicity of newly untested compounds in relation to L. rohita. Finally, the developed best model was employed to predict 297 external chemicals, the most toxic substances to L. rohita were identified as cyhalothrin, isobornyl thiocyanatoacetate, and paclobutrzol, while the least toxic ones included ethyl acetate, ethylthiourea, and n-butyric acid.

First report on chemometric modeling of tilapia fish aquatic toxicity to organic chemicals: Toxicity data gap filling

The Toxic Substances Control Act (TSCA) mandates the Environmental Protection Agency (EPA) to document chemicals entering the US. Due to the vast range of toxicity endpoints, experimental toxicological study for all chemicals is impossible to conduct. To address this, in silico methods like QSAR and read-across are strategically used to prioritize testing for chemicals lacking ecotoxicity data. Aquatic toxicity is one of the most critical endpoints directly related to aquatic species, mainly fish, followed by direct to indirect effects on humans through drinking water and fish as food, respectively. Therefore, we have employed the ToxValDB database to curate acute LC50 toxicity data for three Tilapia species covering two different genera, an ideal species for aquatic toxicity testing. Employing the curated dataset, we have developed multiple robust and predictive QSAR and quantitative read-across structure-activity relationship (q-RASAR) models for Tilapia zillii, Oreochromis niloticus, and Oreochromis mossambicus which helped to understand the toxicological mode of action (MoA) of the modeled chemicals and predict the aquatic toxicity of new untested chemicals followed by toxicity data gap filling. The best three QSAR models showed encouraging statistical quality in terms of determination coefficient R2 (0.94, 0.74, and 0.77), cross-validated leave-one-out Q2 (0.90, 0.67 and 0.70), and predictive capability in terms of R2pred (0.95, 0.77, and 0.74) for T. zillii, O. niloticus, and O. mossambicus datasets, respectively. The developed best mathematical models were used for the prediction of aquatic toxicity in terms of pLC50 for 297 untested organic chemicals across three major Tilapia species ranging from 1.841 to 8.561 M in terms of environmental risk assessment.

q-RASTR modelling for prediction of diverse toxic chemicals towards T. pyriformis

A series of diverse organic compounds impose serious detrimental effects on the health of living organisms and the environment. Determination of the structural aspects of compounds that impart toxicity and evaluation of the same is crucial before public usage. The present study aims to determine the structural characteristics of compounds for Tetrahymena pyriformis toxicity using the q-RASTR (Quantitative Read Across Structure-Toxicity Relationship) model. It was developed using RASTR and 2-D descriptors for a dataset of 1792 compounds with defined endpoint (pIGC50) against a model organism, T. pyriformis. For the current study, the whole dataset was divided based on activity/property into the training and test sets, and the q-RASTR model was developed employing six descriptors (three latent variables) having r2, Q2F1 and Q2 values of 0.739, 0.767, and 0.735, respectively. The generated model was thoroughly validated using internationally recognized internal and external validation criteria to assess the model's dependability and predictability. It was highlighted that high molecular weight, aromatic hydroxyls, nitrogen, double bonds, and hydrophobicity increase the toxicity of organic compounds. The current study demonstrates the applicability of the RASTR algorithm in QSTR model development for the prediction of toxic chemicals (pIGC50) towards T. pyriformis.

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.

Global classification models for predicting acute toxicity of chemicals towards Daphnia magna

Molecular descriptors reflecting structural information on hydrophobicity, reactivity, polarizability, hydrogen bond and charged groups, were used to predict the toxicity (pLC50) of chemicals towards Daphnia magna with global quantitative structure-activity/toxicity relationship (QSAR/QSTR) models. A sufficiently large dataset including 1517 chemical toxicity to Daphnia magna was divided into a training set (758 pLC50) and a test set (759 pLC50). By applying random forest algorithm, two classification models, Class Model A and Class Model B were developed, having prediction accuracy, sensitivity and specificity above 85% for Class 1 (with pLC50 ≤ 4.48) and Class 2 (with pLC50 > 4.48). The Class Model A was based on nine molecular descriptors and RF parameters of nodesize = 1, ntree = 80 and mtry = 2, and yielded accuracy of 92.3% (training set), 85.6% (test set) and 88.9% (total data set). Class Model B was based on ten descriptors and parameters, nodesize = 1, ntree = 90 and mtry = 2, produced accuracy of 88.3% (training set), 86.8% (test set) and 87.5% (total data set). The two classification models were satisfactory compared with other classification model reported in the literature, although classification models in this work dealt with more samples. Thus, the two classification models with a larger applicability domain provided efficient tools for assessing chemical aquatic toxicity towards Daphnia magna.

Quantitative Read-across structure-activity relationship (q-RASAR): A new approach methodology to model aquatic toxicity of organic pesticides against different fish species

We have developed quantitative toxicity prediction models for organic pesticides of agricultural importance considering different fish species using a novel quantitative Read-across structure-activity relationship (q-RASAR) approach. The current study uses experimental (Log 1/LC50) data of organic pesticides to various fish species, including Rainbow trout (RT: Oncorhynchus mykiss: 715 data points), Lepomis (LP: Lepomis macrochirus: 136 data points), and Miscellaneous (Pimephales promelas, Brachydanio rerio: 226 data points). This study has also discussed the validation of the developed models and the analysis of structural features that are important for aquatic toxicity towards fishes. The read-across-derived similarity, error, and concordance measures (RASAR descriptors) have been extracted from the preliminary 0D-2D descriptors; the combined pool of RASAR and selected 0D-2D descriptors have been used to develop the final models by employing partial least squares algorithm. All the q-RASAR models are acceptable in terms of goodness of fit, robustness, and external predictivity, superseding the quality of the respective QSAR models, as seen from the computed validation metrics. The q-RASAR is an effective approach that has the potential to be used as a good alternative way to enhance external predictivity, interpretability, and transferability for aquatic toxicity prediction as well as ecotoxicity potential identification.

QSAR assessment of aquatic toxicity potential of diverse agrochemicals

The fast-increasing number of commercially produced chemicals challenges the experimental ecotoxicity assessment methods, which are costly, time-consuming, and dependent on the sacrifice of animals. In this regard, Quantitative Structure-Property/Activity Relationships (QSPR/QSAR) have led the way in developing ecotoxicity assessment models. In this study, QSAR models have been developed using the pEC50 values of 82 diverse agrochemicals or agro-molecules against a planktonic crustacean Daphnia magna with easily interpretable 2D descriptors. Moreover, a link among octanol-water partition coefficient (KOW), bio-concentration factor (BCF), and critical body residue (CBR) has been addressed, and their imputation for the prediction of the toxicity endpoint (EC50) has been done with an objective of the advanced exploration of several ecotoxicological parameters for toxic chemicals. The developed partial least squares (PLS) models were validated rigorously and proved to be robust, sound, and immensely well-predictive. The final Daphnia toxicity model derived from experimental derived properties along with computed descriptors emerged better in statistical quality and predictivity than those obtained solely from computed descriptors. Additionally, the pEC50 and other important properties (log KOW, log BCF, and log CBR) for a set of external agro-molecules, not employed in model development, were predicted to show the predictive ability of the models.

Virtual screening and invitro evaluation of cyclooxygenase inhibitors from Tinospora cordifolia using the machine learning tool

Tinospora cordifolia has a variety of compounds, and some of these compounds may have anti-inflammatory and antioxidant properties. In the present study, we identified the compounds in the leaf extract of T. cordifolia through Gas Chromatography-Mass Spectrometry (GC-MS) analysis and found the various metabolites. The compounds are screened virtually using a machine learning model, followed by molecular docking and simulation study to identify top-hit compounds as cyclooxygenase (COX) inhibitors. The molecular docking revealed that the compound 7,9-Di-tert-butyl-1-oxaspiro (4,5) deca-6,9-diene-2,8-dione (CID:545303) exhibited the lowest binding energies of -7.1 and -6.8 kcal/mol against COX 1 and COX 2 respectively. The interactions are favored by hydrogen bonding and hydrophobic interaction inside the binding pocket. The 100 ns MD simulation study for these compounds was performed to know the stability and found the RMSD around 2 Å and around 1.0 Å with minimal fluctuations indicating a stable complex throughout the simulation of 100 ns. Based on these findings, we proposed 7,9-Di-tertbutyl- 1-oxaspiro (4,5) deca-6,9-diene-2,8-dione could be used as a dual inhibitor of COX enzymes and a drug-like molecule for treating inflammation after evaluation of their biological properties. The methanolic extract of T. cordifolia was subjected to in vitro DPPH, ABTS, nitric oxide, anti-microbial, COX, and LOX inhibition activity. The results exhibited possible positive effects against the above activities.Communicated by Ramaswamy H. Sarma.

The system of self-consistent models for pesticide toxicity to Daphnia magna

Quantitative structure-property/activity relationships (QSPRs/QSARs) are a tool of modern theoretical and computational chemistry. The self-consistent model system is both a method to build up a group of QSPR/QSAR models and an approach to checking the reliability of these models. Here, a group of models of pesticide toxicity toward Daphnia magna for different distributions into training and test sub-sets is compared. This comparison is the basis for formulating the system of self-consistent models. The so-called index of the ideality of correlation (IIC) has been used to improve the above models' predictive potential of pesticide toxicity. The predictive potential of the suggested models should be classified as high since the average value of the determination coefficient for the validation sets is 0.841, and the dispersion is 0.033 (on all five models). The best model (number 4) has an average determination coefficient of 0.89 for the external validation sets (related to all five splits).

Predictive in silico models for aquatic toxicity of cosmetic and personal care additive mixtures

As emerging environmental contaminants, cosmetic and personal care additives (CPCAs) may have less oversight than other consumer products. Their continuous release and pseudopersistence could cause long-term harm to the aquatic environment. Since CPCAs generally exist in the form of mixtures in the environment, prediction and analysis of their mixture toxicity are crucial for ecological risk assessment. In this study, the acute toxicity of five typical CPCA mixtures to Daphnia magna was tested. The combined toxicity of binary mixtures was examined with the traditional concentration addition (CA) and independent action (IA) model. Overall, the synergistic effect of the five CPCAs may be caused mainly by methylparaben. In addition, reliable approaches for quantitative structure-activity relationship (QSAR) model development were explored. Specifically, 18 QSAR models were developed by three dataset partitioning techniques (Kennard-Stone's algorithm division, Euclidean distance based division, and sorted activity based division), two descriptor filtering methods (genetic algorithm and stepwise multiple linear regression) and three regression methods (multiple linear regression, partial least squares and support vector machine). Sixteen equations were applied for the calculation of the mixture descriptors to screen the functional expression of the mixture descriptors with the largest contribution to the mixture toxicity. A new comprehensive parameter that integrates internal and external validation was proposed for QSAR models evaluation. The mixture toxicity is mainly related the 3D distribution of atomic masses and the spatial distribution of the molecule electronic properties. Rigorously validated and externally predictive QSAR models were developed for predicting the toxicity of binary CPCAs mixtures with any ratio, in the applicability domain. The best possible work frame for construction and validation of QSAR models to provide reliable predictions on the mixture toxicity was proposed.

In silico prediction of acute chemical toxicity of biocides in marine crustaceans using machine learning

Biocides are a heterogeneous group of chemical substances intended to control the growth or kill undesired organisms. Due to their extensive use, they enter marine ecosystems via non-point sources and may pose a threat to ecologically important non-target organisms. Consequently, industries and regulatory agencies have recognized the ecotoxicological hazard potential of biocides. However, the prediction of biocide chemical toxicity on marine crustaceans has not been previously evaluated. This study aims to provide in silico models capable of classifying structurally diverse biocidal chemicals into different toxicity categories and predict acute chemical toxicity (LC₅₀) in marine crustaceans using a set of calculated 2D molecular descriptors. The models were built following the guidelines recommended by the OECD (Organization for Economic Cooperation and Development) and validated through stringent processes (internal and external validation). Six machine learning (ML) models were built and compared (linear regression: LR; support vector machine: SVM; random forest: RF; feed-forward backpropagation-based artificial neural network: ANN; decision trees: DT and naïve Bayes: NB) for regression and classification analysis to predict toxicities. All the models displayed encouraging results with high generalisability: the feed-forward-based backpropagation method showed the best results with determination coefficient R² values of 0.82 and 0.94, respectively, for training set (TS) and validation set (VS). For classification-based modelling, the DT model performed the best with an accuracy (ACC) of 100 % and an area under curve (AUC) value of 1 for both TS and VS. These models showed the potential to replace animal testing for the chemical hazard assessment of untested biocides if they fall within the applicability domain of the proposed models. In general, the models are highly interpretable and robust, with good predictive performance. The models also displayed a trend indicating that toxicity is largely influenced by factors such as lipophilicity, branching, non-polar bonding and saturation of molecules.

Ecotoxicological QSAR modelling of the acute toxicity of fused and non-fused polycyclic aromatic hydrocarbons (FNFPAHs) against two aquatic organisms: Consensus modelling and comparison with ECOSAR

Fused and non-fused polycyclic aromatic hydrocarbons (FNFPAHs) are a type of organic compounds widely occurring in the environment that pose a potential hazard to ecosystem and public health, and thus receive extensive attention from various regulatory agencies. Here, quantitative structure-activity relationship (QSAR) models were constructed to model the ecotoxicity of FNFPAHs against two aquatic species, Daphnia magna and Oncorhynchus mykiss. According to the stringent OECD guidelines, we used genetic algorithm (GA) plus multiple linear regression (MLR) approach to establish QSAR models of the two aquatic toxicity endpoints: D. magna (48 h LC₅₀) and O. mykiss (96 h LC₅₀). The models were established using simple 2D descriptors with explicit physicochemical significance and evaluated using various internal/external validation metrics. The results clearly show that both models are statistically robust (Q_LOO² = 0.7834 for D. magna and Q_LOO² = 0.8162 for O. mykiss), have good internal fitness (R² = 0.8159 for D. magna and R² = 0.8626 for O. mykiss and external predictive ability (D. magna: R_test² = 0.8259, Q_Fn² = 0.7640∼0.8140, CCC_test = 0.8972; O. mykiss:R_test² = 0.8077, Q_Fn² = 0.7615∼0.7722, CCC_test = 0.8910). To prove the predictive performance of the developed models, an additional comparison with the standard ECOSAR tool obviously shows that our models have lower RMSE values. Subsequently, we utilized the best models to predict the true external set compounds collected from the PPDB database to further fill the toxicity data gap. In addition, consensus models (CMs) that integrate all validated individual models (IMs) were more externally predictive than IMs, of which CM2 has the best prediction performance towards the two aquatic species. Overall, the models presented here could be used to evaluate unknown FNFPAHs inside the domain of applicability (AD), thus being very important for environmental risk assessment under current regulatory frameworks.

First report on soil ecotoxicity prediction against Folsomia candida using intelligent consensus predictions and chemical read-across

Soil invertebrates serve as an outstanding biological indicator of the terrestrial ecosystem and overall soil quality, considering their high sensitivity when compared to other indicators of soil quality. In this study, the available soil ecotoxicity data (pEC50) against the soil invertebrate Folsomia candida (C. name: Springtail) (n = 45) were collated from the database of ECOTOX (cfpub.epa.gov/ecotox) and subjected to QSAR analysis using 2D descriptors. Four partial least squares (PLS) models were built based on the features selected through genertic algorithm followed by the best subset selection. These four models were then used as inputs for Intelligent Consensus Predictor version 1.2 (PLS version) to get the final consensus predictions, using the best selection of predictions (compound-wise) from four "qualified" individual models. Both internal and external validations metrics of the consensus predictions are well- balanced and within the acceptable range as per the OECD criteria. The consensus model was found to be better than the previous developed models for this endpoint. Predictions were also made using the Chemical Read-across approach, which showed even better external validation metric values than the consensus predictions. From the selected features in the QSAR models, it has been found out that molecular weight and presence of a di-thiophosphate group, electron donor groups, and polyhalogen substitutions have a significant impact on the soil ecotoxicity. The soil ecotoxicological risk assessment of organic chemicals can therefore be prioritized by these features. The models developed from diverse structural organic compounds can be applied to any new query compound for data gap filling.

Versatile in silico modelling of microplastics adsorption capacity in aqueous environment based on molecular descriptor and machine learning

To comprehensively evaluate the hazards of microplastics and their coexisting organic pollutants, the sorption capacity of microplastics is a major issue that is quantified through the microplastic-aqueous sorption coefficient (K_d). Almost all quantitative structure-property relationship (QSPR) models that describe K_d apply only to narrow, relatively homogeneous groups of reactants. Herein, non-hybrid QSPR-based models were developed to predict PE-water (K_PE-w), PE-seawater (K_PE-sw), PVC-water (K_PVC-w) and PP-seawater (K_PP-sw) sorption coefficients at different temperatures, with eight machine learning algorithms. Moreover, novel hybrid intelligent models for predicting K_d more accurately were innovatively developed by applying GA, PSO and AdaBoost algorithms to optimize MLP and ELM models. The results indicated that all three optimization algorithms could improve the robustness and predictability of the standalone MLP and ELM models. In all models trained with K_PE-w, K_PE-sw, K_PVC-w and K_PP-sw data sets, GBDT-1 and XGBoost-1 models, MLP-GA-2 and MLP-PSO-2 models, MLR-3 and MLR-4 models performed better in terms of goodness of fit (R_adj²: 0.907-0.999), robustness (Q_BOOT²: 0.900-0.937) and predictability (R_ext²: 0.889-0.970), respectively. Analyzing the descriptors revealed that temperature, lipophilicity, ionization potential and molecular size were correlated closely with the adsorption capacity of microplastics to organic pollutants. The proposed QSPR models may assist in initial environmental exposure assessments without imposing heavy costs in the early experimental phase.

Prediction of acute toxicity of pesticides for Americamysis bahia using linear and nonlinear QSTR modelling approaches

Globally, pesticides are toxic substances with wide applications. However, the widespread use of pesticides has received increasing attention from regulatory agencies due to their various acute and chronic effects on multiple organisms. In this study, Quantitative Structure-Toxicity Relationship (QSTR) models were established using Multiple Linear Regression (MLR) and five Machine Learning (ML) algorithms to predict pesticide toxicity in Americamysis bahia. The most influential descriptors included in the MLR model are RBF, JGI2, nCbH, nRCOOR, nRSR, nPO4 and 'Cl-090', with positive contributions to the dependent variable (negative decimal logarithm of median lethal concentration at 96-h). The Random Forest (RF) regression model was superior amongst the five ML models. We observed higher values of R² (0.812) and lower values of RMSE (0.595) and MAE (0.462) in the cross-validation training set and external validation set. Similarly, this study had a high level of fitness and was internally robust and externally predictive compared to models presented in similar studies. The results suggest that the developed QSTR models are suitable for reliably predicting the aquatic toxicity of structurally diverse pesticides and can be used for screening, prioritising new pesticides, filling data gaps and overcoming the limitations of in vivo and in vitro tests.

Predicting the Surface Tension of Deep Eutectic Solvents: A Step Forward in the Use of Greener Solvents

Deep eutectic solvents (DES) are an important class of green solvents that have been developed as an alternative to toxic solvents. However, the large-scale industrial application of DESs requires fine-tuning their physicochemical properties. Among others, surface tension is one of such properties that have to be considered while designing novel DESs. In this work, we present the results of a detailed evaluation of Quantitative Structure-Property Relationships (QSPR) modeling efforts designed to predict the surface tension of DESs, following the Organization for Economic Co-operation and Development (OECD) guidelines. The data set used comprises a large number of structurally diverse binary DESs and the models were built systematically through rigorous validation methods, including ‘mixtures-out’- and ‘compounds-out’-based data splitting. The most predictive individual QSPR model found is shown to be statistically robust, besides providing valuable information about the structural and physicochemical features responsible for the surface tension of DESs. Furthermore, the intelligent consensus prediction strategy applied to multiple predictive models led to consensus models with similar statistical robustness to the individual QSPR model. The benefits of the present work stand out also from its reproducibility since it relies on fully specified computational procedures and on publicly available tools. Finally, our results not only guide the future design and screening of novel DESs with a desirable surface tension but also lays out strategies for efficiently setting up silico-based models for binary mixtures.

A critical review of effect modeling for ecological risk assessment of plant protection products

A wide diversity of plant protection products (PPP) is used for crop protection leading to the contamination of soil, water, and air, which can have ecotoxicological impacts on living organisms. It is inconceivable to study the effects of each compound on each species from each compartment, experimental studies being time consuming and cost prohibitive, and animal testing having to be avoided. Therefore, numerous models are developed to assess PPP ecotoxicological effects. Our objective was to provide an overview of the modeling approaches enabling the assessment of PPP effects (including biopesticides) on the biota. Six categories of models were inventoried: (Q)SAR, DR and TKTD, population, multi-species, landscape, and mixture models. They were developed for various species (terrestrial and aquatic vertebrates and invertebrates, primary producers, micro-organisms) belonging to diverse environmental compartments, to address different goals (e.g., species sensitivity or PPP bioaccumulation assessment, ecosystem services protection). Among them, mechanistic models are increasingly recognized by EFSA for PPP regulatory risk assessment but, to date, remain not considered in notified guidance documents. The strengths and limits of the reviewed models are discussed together with improvement avenues (multigenerational effects, multiple biotic and abiotic stressors). This review also underlines a lack of model testing by means of field data and of sensitivity and uncertainty analyses. Accurate and robust modeling of PPP effects and other stressors on living organisms, from their application in the field to their functional consequences on the ecosystems at different scales of time and space, would help going toward a more sustainable management of the environment. Graphical Abstract Combination of the keyword lists composing the first bibliographic query. Columns were joined together with the logical operator AND. All keyword lists are available in Supplementary Information at https://doi.org/10.5281/zenodo.5775038 (Larras et al. 2021).

A QSTR model for toxicity prediction of pesticides towards Daphnia magna

Because of the large amount of pesticides discharged into rivers, adverse effects could be induced to aquatic organisms. Daphnia magna is often used as an indicator organism to evaluate the toxicity of pesticides. In this study, a quantitative structure-toxicity relationship (QSTR) model was established based on norm descriptors for predicting the acute toxicity of pesticides to Daphnia magna. The model results showed the good predictability (R_training² = 0.8045, R_testing² = 0.8224). The validation results of internal validation, external validation, Y-randomization test and application domain analysis demonstrated the model's stability, reliability and robustness. Therefore, the above results indicate that norm descriptors might be universal for describing the relationship between the toxicity and structures of pesticides compounds. Moreover, some pesticides' toxicities without experimental data were also predicted by this model.

Surface functional groups determine adsorption of pharmaceuticals and personal care products on polypropylene microplastics

The pervasiveness of microplastics (MPs), which can absorb pharmaceuticals and personal care products (PPCPs), has a certain impact on pollutant migration in natural waters. The adsorption behaviors of PPCPs on the aged polypropylene (PP) followed the pseudo-second-order kinetics and Langmuir isotherm, and the adsorption capacity (q_e) on the aged PP was much higher than that on the fresh PP. The Weber-Morris and Boyd models confirmed that the liquid-film and intra-particle diffusion affected the adsorption of PPCPs on the aged PP while the surface diffusion was a rate-limiting step for the fresh PP. The analysis of SEM-EDS, BET, FT-IR, and XPS further showed that changes in the type and content of the surface functional groups of PP led to differences in adsorption capacity and adsorption interactions. The Dragon-descriptor-based LFER and the quantum-chemical-descriptor-based QSAR models reflected the difference in adsorption interaction mechanisms. The examined models showed that the adsorption of the fresh PP toward PPCPs relied on hydrophobic and hydrogen bonding interaction, while for the aged PP electrostatic interaction and hydrogen bonding controlled the adsorption. The findings clarified interactions between PPCPs and MPs and provided a theoretical basis for the assessment of environmental behavior and ecological risk when MPs and PPCPs coexist.

Chemometric modeling of plant protection products (PPPs) for the prediction of acute contact toxicity against honey bees (A. mellifera): A 2D-QSAR approach

Honey bees (Apis mellifera) are vital for economic, viable agriculture and for food safety. Although Plant Protection Products (PPPs) are of undeniable importance in the global agricultural system, these have become potential threats for non-target organisms like pollinators (e.g., honey bees etc.), resulting in the disruption of the ecological balance. In the current work, we have used the 113 PPP analogs to develop a 2D-QSAR model and explored the structural features modulating the toxic effects on honey bees, following the Organization for Economic Co-operation and Development (OECD) guidelines. The extensive validation of the developed model has been performed using internal and external validation metrics to make sure that the model is statistically sound and interpretable enough to be acceptable. The obtained results (R² = 0.666, Q² = 0.594, Q²_F1 = 0.647 and Q²_F2 = 0.646) determine the predictability and reliability of the developed model. This model should be useful for the predictions (acute contact toxicity (LD₅₀)) of the new and untested compounds located inside the applicability domain of the developed model. Moreover, we have performed the in-silico prediction of toxicity against honey bees of a total of 709 compounds obtained from the pesticide properties database (PPDB) using the developed model.

Systematic QSAR and iQCCR modelling of fused/non-fused aromatic hydrocarbons (FNFAHs) carcinogenicity to rodents: reducing unnecessary chemical synthesis and animal testing

The prediction of new or untested FNFAHs will reduce unnecessary chemical synthesis and animal testing, and contribute to the design of safer chemicals for production activities.

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.

Chemometric QSAR modeling of acute oral toxicity of Polycyclic Aromatic Hydrocarbons (PAHs) to rat using simple 2D descriptors and interspecies toxicity modeling with mouse

The information of the acute oral toxicity for most polycyclic aromatic hydrocarbons (PAHs) in mammals are lacking due to limited experimental resources, leading to a need to develop reliable in silico methods to evaluate the toxicity endpoint. In this study, we developed the quantitative structure-activity relationship (QSAR) models by genetic algorithm (GA) and multiple linear regression (MLR) for the rat acute oral toxicity (LD₅₀) of PAHs following the strict validation principles of QSAR modeling recommended by OECD. The best QSAR model comprised eight simple 2D descriptors with definite physicochemical meaning, which showed that maximum atom-type electrotopological state, van der Waals surface area, mean atomic van der Waals volume, and total number of bonds are main influencing factors for the toxicity endpoint. A true external set (554 compounds) without rat acute oral toxicity values, and 22 limit test compounds, were firstly predicted along with reliability assessment. We also compared our proposed model with the OPERA predictions and recently published literature to prove the prediction reliability. Furthermore, the interspecies toxicity (iST) models of PAHs between rat and mouse were also established, validated and employed for filling data gap. Overall, our developed models should be applicable to new or untested or not yet synthesized PAHs falling within the applicability domain (AD) of the models for rapid acute oral toxicity prediction, thus being important for environmental or personal exposure risk assessment under regulatory frameworks.

Identification of concealed structural alerts using QSTR modeling for Pseudokirchneriella subcapitata

In the present work, QSTR modeling was conducted for microalga Pseudokirchneriella subcapitata using a data set of 271 molecules belonging to different types of chemical classes for the prediction of EC₅₀ for 72 hr based assays. The balanced QSTR model encompasses seven easily interpretable molecular descriptors and possesses statistical robustness with high predictive ability. This Genetic Algorithm Multi-linear regression (GA-MLR) model was subjected to internal validation, Y-randomization test, applicability domain analysis, and external validation as per the recommended OECD guidelines. The newly developed model fulfilled the threshold values for more than 20 recommended validation parameters including R² = 0.72, Q²_LOO = 0.70, etc. The developed QSTR model was successful in identifying the type of hybridization or specific type of atoms of previously reported and newer structural alerts. Thus, the model could be useful for data gap filling and expanding mechanistic interpretation of toxicity for different chemicals.

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 LC_50, 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.

Density of Deep Eutectic Solvents: The Path Forward Cheminformatics-Driven Reliable Predictions for Mixtures

Deep eutectic solvents (DES) are often regarded as greener sustainable alternative solvents and are currently employed in many industrial applications on a large scale. Bearing in mind the industrial importance of DES—and because the vast majority of DES has yet to be synthesized—the development of cheminformatic models and tools efficiently profiling their density becomes essential. In this work, after rigorous validation, quantitative structure-property relationship (QSPR) models were proposed for use in estimating the density of a wide variety of DES. These models were based on a modelling dataset previously employed for constructing thermodynamic models for the same endpoint. The best QSPR models were robust and sound, performing well on an external validation set (set up with recently reported experimental density data of DES). Furthermore, the results revealed structural features that could play crucial roles in ruling DES density. Then, intelligent consensus prediction was employed to develop a consensus model with improved predictive accuracy. All models were derived using publicly available tools to facilitate easy reproducibility of the proposed methodology. Future work may involve setting up reliable, interpretable cheminformatic models for other thermodynamic properties of DES and guiding the design of these solvents for applications.

Development of blood brain barrier permeation prediction models for organic and inorganic biocidal active substances

Biocidal products are broadly used in homes and industries. However, the safety of biocidal active substances (BASs) is not yet fully understood. In particular, the neurotoxic action of BASs needs to be studied as diverse epidemiological studies have reported associations between exposure to BASs and neural diseases. In this study, we developed in silico models to predict the blood-brain barrier (BBB) permeation of organic and inorganic BASs. Due to a lack of BBB data for BASs, the chemical space of BASs and BBB dataset were compared in order to select BBB data that were structurally similar to BASs. In silico models to predict log-scaled BBB penetration were developed using support vector regression for organic BASs and multiple linear regression for inorganic BASs. The model for organic BASs was developed with 231 compounds (training set: 153 and test set: 78) and achieved good prediction accuracy on an external test set (R² = 0.64), and the model outperformed the model for pharmaceuticals. The model for inorganic BASs was developed with 11 compounds (R² = 0.51). Applicability domain (AD) analysis of the models clarified molecular structures reliably predicted by the models. Therefore, the models developed in this study can be used for predicting BBB permeable BASs in human. These models were developed according to the Quantitative Structure-Activity Relationship validation principles proposed by the Organization for Economic Cooperation and Development.

Support vector machine-based model for toxicity of organic compounds against fish

Predicting the toxicity of chemicals to various fish species through chemometric approach is crucial for ecotoxicological assessment of existing as well as not yet synthesized chemicals. This paper reports a quantitative structure-activity/toxicity relationship (QSAR/QSTR) model for the toxicity pLC₅₀ of organic chemicals against various fish species. Only six descriptors were used to develop the QSTR model, by applying support vector machine (SVM) together with genetic algorithm. The QSTR model was trained and established on a sufficiently large data set of 840 organic compounds and evaluated with a test set (281 compounds). Compared with other QSTRs reported in the literature, the optimal SVM model for fish toxicity produces better statistical results with determination coefficients R² above 0.70 for both the training set and test set, although the QSTR model in this work possesses fewer molecular descriptors. Applying SVM and genetic algorithm to develop the QSTR model for pLC₅₀ of organic compounds against various fish species is successful.

Predictive classification-based QSTR models for toxicity study of diverse pesticides on multiple avian species

Protection and restoration of different endangered bird species from pesticide exposure is crucial from the point of safety assessment of ecosystem. Toxicity predictions or risk assessment of pesticides by chemometric tools is one of the challenging fields in recent era. In the present study, classification-based quantitative structure toxicity relationship (QSTR) models were developed for a large dataset (516) of diverse pesticides on multiple avian species mallard duck, bobwhite quail, and zebra finch according to the Organization for Economic Co-operation and Development guidelines. The QSTR models were developed by linear discriminant analysis method with genetic algorithm for feature selection from 2D descriptors using QSAR-Co software. Different statistical metrics assured the reliability and robustness of the developed models. External compound prediction highlighted predictive nature of the models. The mechanistic interpretation suggested that presence of phosphate, halogens (Cl, Br), ether linkage, and NCOO influence the avian toxicity. Furthermore, model reliability was checked by the application of the standardization approach of the applicability domain (AD). Finally, the developed models provided a priori toxic and non-toxic classification for unknown pesticides (inside AD), with particular emphasis on organophosphate pesticides. The interspecies toxicity correlation and predictions encouraged for their further applicability for the fulfilment of data gaps in vital missing species.

Defining the Human-Biota Thresholds of Toxicological Concern for Organic Chemicals in Freshwater: The Proposed Strategy of the LIFE VERMEER Project Using VEGA Tools

Several tons of chemicals are released every year into the environment and it is essential to assess the risk of adverse effects on human health and ecosystems. Risk assessment is expensive and time-consuming and only partial information is available for many compounds. A consolidated approach to overcome this limitation is the Threshold of Toxicological Concern (TTC) for assessment of the potential health impact and, more recently, eco-TTCs for the ecological aspect. The aim is to allow a safe assessment of substances with poor toxicological characterization. Only limited attempts have been made to integrate the human and ecological risk assessment procedures in a "One Health" perspective. We are proposing a strategy to define the Human-Biota TTCs (HB-TTCs) as concentrations of organic chemicals in freshwater preserving both humans and ecological receptors at the same time. Two sets of thresholds were derived: general HB-TTCs as preliminary screening levels for compounds with no eco- and toxicological information, and compound-specific HB-TTCs for chemicals with known hazard assessment, in terms of Predicted No effect Concentration (PNEC) values for freshwater ecosystems and acceptable doses for human health. The proposed strategy is based on freely available public data and tools to characterize and group chemicals according to their toxicological profiles. Five generic HB-TTCs were defined, based on the ecotoxicological profiles reflected by the Verhaar classes, and compound-specific thresholds for more than 400 organic chemicals with complete eco- and toxicological profiles. To complete the strategy, the use of in silico models is proposed to predict the required toxicological properties and suitable models already available on the VEGAHUB platform are listed.

In vivo toxicity of nitroaromatic compounds to rats: QSTR modelling and interspecies toxicity relationship with mouse

Nitroaromatic compounds (NACs) in the environment can cause serious public health and environmental problems due to their potential toxicity. This study established quantitative structure-toxicity relationship (QSTR) models for the acute oral toxicity of NACs towards rats following the stringent OECD principles for QSTR modelling. All models were assessed by various internationally accepted validation metrics and the OECD criteria. The best QSTR model contains seven simple and interpretable 2D descriptors with defined physicochemical meaning. Mechanistic interpretation indicated that van der Waals surface area, presence of C-F at topological distance 6, heteroatom content and frequency of C-N at topological distance 9 are main factors responsible for the toxicity of NACs. This proposed model was successfully applied to a true external set (295 compounds), and prediction reliability was analysed and discussed. Moreover, the rat-mouse and mouse-rat interspecies quantitative toxicity-toxicity relationship (iQTTR) models were also constructed, validated and employed in toxicity prediction for true external sets consisting of 67 and 265 compounds, respectively. These models showed good external predictivity that can be used to rapidly predict the rat oral acute toxicity of new or untested NACs falling within the applicability domain of the models, thus being beneficial in environmental risk assessment and regulatory purposes.

Ecosystem ecology: Models for acute toxicity of pesticides towards Daphnia magna

Quantitative structure - activity relationships (QSARs) which are obtained with a representation of the molecular architecture via simplified molecular input-line entry system (SMILES) are applied to build up predictive models of acute toxicity of pesticides towards Daphnia magna. The acute toxicity towards Daphnia magna is an adequate measure of the ecological impact of various substances. The Monte Carlo technique is the basis to build up the above QSAR models. The statistical quality of suggested models is good: the best model is characterized by n = 103, R² = 0.76, RMSE = 0.91 (training set); n = 53, R² = 0.82, RMSE = 0.87 (validation set). The approach provides the mechanistic interpretation (e.g. aromaticity and branching of carbon skeleton are promoters of increase for toxicity towards Daphnia magna in the case of the examined set of pesticides). The approach is attractive to build up predictive models since instead of a large number of different molecular descriptors the corresponding model is based on solely one optimal descriptor calculated with SMILES and all necessary calculations can be done using the CORAL software available on the Internet (http://ww.insilico.eu/coral).

QSAR modeling the toxicity of pesticides against Americamysis bahia

The widespread use of pesticides has received increasing attention in regulatory agencies because their extensive overuse and various adverse effects on all living organisms. Organizations such as EPA and ECHA have published laws that pesticides should be fully evaluated before bring them to market. In the present study, we evaluated the pesticides toxicity using the Quantitative Structural-Activity Relationship (QSAR) method. The models for the single class pesticides (herbicides, insecticides and fungicides) as well as the general class pesticides (the combined dataset plus some microbicides, molluscicides, etc.) were developed using the Genetic Algorithm and Multiple Linear Regression method. The internal and external validation results suggested that all the obtained models were stable and predictive. According to the modeling descriptors, the lipophilic descriptors contributed positively while all the electrotopological state descriptors showed a negative contribution, their presences in every model verified the conspicuous influence of molecular lipophilicity and hydrophilicity on the pesticides toxicity. However, the influence of topological structure descriptors was different and varies with the physiochemical information they encode. Finally, the models presented in this paper would help assess the pesticides toxicity against Americamysis bahia, shorten test time, and reduce the cost of pesticides risk assessment.

Evaluation of modes of action of pesticides to Daphnia magna based on QSAR, excess toxicity and critical body residues

Agricultural pesticides serve as effective controls of unwanted weeds and pests. However, these same chemicals can exert toxic effects in non-target organisms. To determine chemical modes of action, the toxicity ratio (TR) and critical body residues (CBRs) of 57 pesticides were calculated for Daphnia magna. Results showed that the CBR values of inert compounds were close to a constant while the CBR values of pesticides varied over a wider range. Although herbicides are categorized as specifically-acting compounds to plants, herbicides did not exhibit excess toxicity to Daphnia magna and were categorized as inert compounds with an average logTR = 0.41, which was less than a threshold of one. Conversely, fungicides and insecticides exhibited strong potential for toxic effects to Daphnia magna with an average logTR >2. Many of these chemicals act via disruption of the nervous, respiratory, or reproductive system, with high ligand-receptor binding activity which leads to higher toxicity for Daphnia magna. Molecular docking using acetylcholinesterase revealed that fungicides and insecticides bind more easily with the biological macromolecule when compared with inert compounds. Quantitative structure-activity relationship (QSAR) analysis revealed that the toxicity of fungicides was mainly dependent upon the heat of formation and polar surface area, while the toxicity of insecticides was more related to hydrogen-bond properties. This comprehensive analysis reveals that there are specific differences in toxic mechanisms between fungicides and insecticides. These results are useful for determining relative risk associated with pesticide exposure to aquatic crustaceans, such as Daphnia magna.

QSAR models for the acute toxicity of 1,2,4-triazole fungicides to zebrafish (Danio rerio) embryos

In recent decades, the 1,2,4-triazole fungicides are widely used for crop diseases control, and their toxicity to wild lives and pollution to ecosystem have attracted more and more attention. However, how to quickly and efficiently evaluate the toxicity of these compounds to environmental organisms is still a challenge. In silico method, such like Quantitative Structure-Activity Relationship (QSAR), provides a good alternative to evaluate the environmental toxicity of a large number of chemicals. At the present study, the acute toxicity of 23 1,2,4-triazole fungicides to zebrafish (Danio rerio) embryos was firstly tested, and the LC₅₀ (median lethal concentration) values were used as the bio-activity endpoint to conduct QSAR modelling for these triazoles. After the comparative study of several QSAR models, the 2D-QSAR model was finally constructed using the stepwise multiple linear regression algorithm combining with two physicochemical parameters (logD and μ), an electronic parameter (Q_N1) and a topological parameter (X^vPC₄). The optimal model could be mathematically described as following: pLC₅₀ = -7.24-0.30X^vPC₄ + 0.76logD - 26.15Q_N1 - 0.08μ. The internal validation by leave-one-out (LOO) cross-validation showed that the R²_adj (adjusted noncross-validation squared correlation coefficient), Q² (cross-validation correlation coefficient) and RMSD (root-mean-square error) was 0.88, 0.84 and 0.17, respectively. The external validation indicated the model had a robust predictability with the q² (predictive squared correlation coefficient) of 0.90 when eliminated tricyclazole. The present study provided a potential tool for predicting the acute toxicity of new 1,2,4-triazole fungicides which contained an independent triazole ring group in their molecules to zebrafish embryos, and also provided a reference for the development of more environmentally-friendly 1,2,4-triazole pesticides in the future.

QSAR model for pesticides toxicity to Rainbow Trout based on "ideal correlations"

Pesticides have an impact on the aquatic environment, with ecological effects. The regulation of this impact is of key importance. One of the components of the planning of agricultural and industrial activities is the development of databases and models in order to identify substances that may cause damage. In this study, a quantitative structure-activity relationship (QSAR) approach was established for the prediction of acute toxicity toward rainbow trout of various pesticides. The so-called index of ideality of correlation is the main component of this approach. The validation of this approach has been carried out with three random splits into the training and validation sets. The range of statistical quality of models obtained here for the validation set is R² = [0.81-0.86] and RMSE = [0.55-0.65].

Consensus QSAR models estimating acute toxicity to aquatic organisms from different trophic levels: algae, Daphnia and fish

We report new consensus models estimating acute toxicity for algae, Daphnia and fish endpoints. We assembled a large collection of 3680 public unique compounds annotated by, at least, one experimental value for the given endpoint. Support Vector Machine models were internally and externally validated following the OECD principles. Reasonable predictive performances were achieved (RMSE_ext = 0.56-0.78) which are in line with those of state-of-the-art models. The known structural alerts are compared with analysis of the atomic contributions to these models obtained using the ISIDA/ColorAtom utility. A benchmarking against existing tools has been carried out on a set of compounds considered more representative and relevant for the chemical space of the current chemical industry. Our model scored one of the best accuracy and data coverage. Nevertheless, industrial data performances were noticeably lower than those on public data, indicating that existing models fail to meet the industrial needs. Thus, final models were updated with the inclusion of new industrial compounds, extending the applicability domain and relevance for application in an industrial context. Generated models and collected public data are made freely available.

A Multi-layered Variable Selection Strategy for QSAR Modeling of Butyrylcholinesterase Inhibitors

Background:

Alzheimer’s disease (AD), a neurological disorder, is the most common cause of senile dementia. Butyrylcholinesterase (BuChE) enzyme plays a vital role in regulating the brain acetylcholine (ACh) neurotransmitter, but in the case of Alzheimer’s disease (AD), BuChE activity gradually increases in patients with a decrease in the acetylcholine (ACh) concentration via hydrolysis. ACh plays an essential role in regulating learning and memory as the cortex originates from the basal forebrain, and thus, is involved in memory consolidation in these sites.

Methods:

In this work, we have developed a partial least squares (PLS)-regression based two dimensional quantitative structure-activity relationship (2D-QSAR) model using 1130 diverse chemical classes of compounds with defined activity against the BuChE enzyme. Keeping in mind the strict Organization for Economic Co-operation and Development (OECD) guidelines, we have tried to select significant descriptors from the large initial pool of descriptors using multi-layered variable selection strategy using stepwise regression followed by genetic algorithm (GA) followed by again stepwise regression technique and at the end best subset selection prior to development of final model thus reducing noise in the input. Partial least squares (PLS) regression technique was employed for the development of the final model while model validation was performed using various stringent validation criteria.

Results:

The results obtained from the QSAR model suggested that the quality of the model is acceptable in terms of both internal (R2= 0.664, Q2= 0.650) and external (R2 Pred= 0.657) validation parameters. The QSAR studies were analyzed, and the structural features (hydrophobic, ring aromatic and hydrogen bond acceptor/donor) responsible for enhancement of the activity were identified. The developed model further suggests that the presence of hydrophobic features like long carbon chain would increase the BuChE inhibitory activity and presence of amino group and hydrazine fragment promoting the hydrogen bond interactions would be important for increasing the inhibitory activity against BuChE enzyme.

Conclusion:

Furthermore, molecular docking studies have been carried out to understand the molecular interactions between the ligand and receptor, and the results are then correlated with the structural features obtained from the QSAR models. The information obtained from the QSAR models are well corroborated with the results of the docking study.

Exploring QSAR models for assessment of acute fish toxicity of environmental transformation products of pesticides (ETPPs)

Environmental transformation products of pesticides (ETPPs) have a great deal of ecological impact owing to their ability to cause toxicity to the aquatic organisms, which can then be translated to the humans. The limited experimental data on biochemical and toxic effects of ETPPs, the high test costs together with regulatory limitations and the international push to reduce animal testing encourage greater dependence on predictive in silico techniques like quantitative structure-activity relationship (QSAR) models. The aim of the present work was to explore the key structural features, which regulate the toxicity towards fishes, for 85 ETPPs using a partial least squares (PLS) regression based chemometric model developed according to Organisation for Economic Co-operation and Development (OECD) guidelines. The model was extensively validated using both internal and external validation metrics, and the results so obtained justify the reliability and usefulness of the developed model (Q² = 0.648, R²_pred or Q²_F1 = 0.734 and Q²_F2 = 0.733). From the developed model, we can conclude that lipophilicity, polarity, presence of branching and the functional form of O-atom in the transformed structures of pesticides are the important features that are to be considered during ecotoxicity assessment of ETPPs. The information obtained from the descriptors of the developed model could be utilized in the future for assessing ETPPs with the benefit of providing an early warning of their potentially detrimental effect on fishes for regulatory purposes.

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.

Inclusion of molecular descriptors in predictive models improves pesticide soil-air partitioning estimates

The soil-air exchange of pesticides is one potential fate and exposure pathways, and this process is generally thought to be governed by soil properties and environmental conditions. The experimental determination of soil-air partitioning coefficient (Ksa) is laborious and costly and typically, Ksa's are predicted from a semiempirical or a simple linear regression approach with soil and environmental variables. Here we developed a model that combined linear regression of soil, environmental and molecular parameters with the quantitative structural-property relationship (QSPR) to predict Ksa for pesticides. The values of theoretical descriptors of pesticides were calculated and the best descriptors selected using the Boruta Algorithm. Seventy-six experimental logKsa values for 17 pesticides were used in model development. Multiple linear regression (MLR) with a soil (organic carbon fraction), physicochemical (octanol-air partitioning coefficient), environmental (temperature and humidity) and molecular descriptor (Gmin, a 2D E-state molecular parameter), called as MLR-QSPR combined model exhibited better predictability (adj. r² = 0.95) of logKsa compared to MLR (adj. r² = 0.87) or QSPR (adj. r² = 0.82) itself. MLR-QSPR also showed the best performance in five-fold cross-validation (adj. r² = 0.94) and test set verification (adj. r² = 0.96). The developed model was validated and characterized by the applicability domain. Results showed that the proposed MLR-QSPR approach is highly predictive and statistically robust with >95% of predictions within ±0.5 log unit of the measured Ksa. Therefore, this approach can be used in estimating the soil-air partitioning of pesticides to better predict it's fate and transport in environments.

Prediction Model of Aryl Hydrocarbon Receptor Activation by a Novel QSAR Approach, DeepSnap-Deep Learning

The aryl hydrocarbon receptor (AhR) is a ligand-dependent transcription factor that senses environmental exogenous and endogenous ligands or xenobiotic chemicals. In particular, exposure of the liver to environmental metabolism-disrupting chemicals contributes to the development and propagation of steatosis and hepatotoxicity. However, the mechanisms for AhR-induced hepatotoxicity and tumor propagation in the liver remain to be revealed, due to the wide variety of AhR ligands. Recently, quantitative structure–activity relationship (QSAR) analysis using deep neural network (DNN) has shown superior performance for the prediction of chemical compounds. Therefore, this study proposes a novel QSAR analysis using deep learning (DL), called the DeepSnap–DL method, to construct prediction models of chemical activation of AhR. Compared with conventional machine learning (ML) techniques, such as the random forest, XGBoost, LightGBM, and CatBoost, the proposed method achieves high-performance prediction of AhR activation. Thus, the DeepSnap–DL method may be considered a useful tool for achieving high-throughput in silico evaluation of AhR-induced hepatotoxicity.

QSPR Studies on the Octane Number of Toluene Primary Reference Fuel Based on the Electrotopological State Index

The quantitative structure-property relationship (QSPR) models for predicting the octane number (ON) of toluene primary reference fuel (TPRF; blends of n-heptane, isooctane, and toluene) was investigated. The electrotopological state (E-state) index of TPRF components was computed and weight-summed to generate the quantitative descriptor of TPRF samples. The partial least squares (PLS) technique was used to build up the regression model between the ON and weight-summed E-state index of the investigated samples. The QSPR models for the research octane number (RON) and motor octane number (MON) of TPRF were built. The prediction performance of the obtained PLS models was assessed by the external test set validation and leave-one-out cross-validation. The validation results demonstrate that the proposed PLS models are feasible for predicting the ON, both RON and MON, of TPRF. In addition, several other QSPR models for the ON of TPRF were developed by employing the stepwise regression and Scheffé polynomials methods, and the prediction performance of these models were compared with that of the PLS models. The comparison result shows that the proposed PLS models are slightly better than multiple linear regression models and Scheffé models. It is demonstrated that the combination of the E-state index and PLS is an easy-to-use and promising method for studying and forecasting the ON of TPRF.