QSAR models for biocides: The example of the prediction of Daphnia magna acute toxicity

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.

In silico environmental risk assessment improves efficiency for pesticide safety management

Pesticides are indispensable to maintain crop quality and food production worldwide, but their use also poses environmental risks. Pesticide risk assessment involves a series of complex, expensive and time-consuming toxicity tests. To improve the efficiency and accuracy for assessing the environmental impact of pesticides, numerous computational tools have been developed. However, there is a notable deficiency in critical analysis or a systematic summary of environmental risk assessment tools and their applicable contexts. Here, many of the current approaches and tools for assessing environmental risks posed by pesticides are reviewed, and the question of whether these tools are fit for use on complex multicomponent scenarios is discussed. We analyze the adaptations of these tools to aquatic and terrestrial ecosystems, followed by the provision of resources for predicting pesticide concentrations in environmental medias, including air, soil and water. The successful application of computational tools for risk assessment and interpretation of predicted results will also be discussed. This assessment serves as a valuable resource, enabling scientists to utilize suitable models to enhance the robustness of pesticides risk assessments.

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.

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).

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.

In silico prediction of mosquito repellents for clothing application

Use of protective clothing is a simple and efficient way to reduce the contacts with mosquitoes and consequently the probability of transmission of diseases spread by them. This mechanical barrier can be enhanced by the application of repellents. Unfortunately the number of available repellents is limited. As a result, there is a crucial need to find new active and safer molecules repelling mosquitoes. In this context, a structure-activity relationship (SAR) model was proposed for the design of repellents active on clothing. It was computed from a dataset of 2027 chemicals for which repellent activity on clothing was measured against Aedes aegypti. Molecules were described by means of 20 molecular descriptors encoding physicochemical properties, topological information and structural features. A three-layer perceptron was used as statistical tool. An accuracy of 87% was obtained for both the training and test sets. Most of the wrong predictions can be explained. Avenues for increasing the performances of the model have been proposed.

QSAR modelling of organic dyes for their acute toxicity in Daphnia magna using 2D-descriptors

The present study reports quantitative structure-activity relationship (QSAR) models for 22 organic dyes spanning a broad chemical domain to predict their toxicity in Daphnia magna [log (1/EC₅₀)]. Only two-dimensional descriptors with clear physicochemical meaning were used to construct the QSAR models. The process of development, validation, and interpretation of models adheres to the stringent recommendations of the Organization for Economic Cooperation and Development (OECD) guidelines. In this study, the multi-layered stepwise regression method and linear discriminant analysis (LDA) method were employed for the deployment of regression - and classification-based models respectively; however, the final regression-based QSAR models were obtained through the partial least squares (PLS) regression. Additionally, the applicability domain of the developed models was verified. The constructed models should be applicable in the absence of toxicity data of new or untested dye structures, particularly when the compounds fall within the developed models' scope, and also implementable to develop more environmentally friendly alternatives.

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.

The Monte Carlo method to build up models of the hydrolysis half-lives of organic compounds

The hydrolysis of organic chemicals such as pesticides, pollutants, or drugs can affect the fate and behaviour of environmental contaminants, so it is of interest to evaluate the stability of substances in water for various purposes. For the registration of organic compounds in Europe, information on hydrolysis must be presented. However, the experimental measurements of all chemicals would require enormous resources, and computational models may become attractive. Applying the CORAL software (http://www.insilico.eu/coral) quantitative structure-property relationships (QSPRs) were built up to model hydrolysis. The 2D-optimal descriptor is calculated with so-called correlation weights for attributes of simplified molecular input-line entry systems (SMILES). The correlation weights are obtained as results of the special Monte Carlo optimization. The nature of (five- and six-member) rings is an important component of this approach. Another important component is the atom pair proportions for nitrogen, oxygen, and sulphur. The statistical quality of the best model is: n = 44, r² = 0.74 (training set); n = 14, r² = 0.75 (calibration set); and n = 12, r² = 0.80 (validation set).

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.

A Monte Carlo method based QSPR model for prediction of reaction rate constants of hydrated electrons with organic contaminants

The Monte Carlo algorithm was applied to formulate a robust quantitative structure-property relationship (QSPR) model to compute the reactions rate constants of hydrated electron values for a data set of 309 water contaminants containing 125 aliphatic and 184 phenyl-based chemicals. The QSPR models were computed with the hybrid optimal descriptors which were procured by combining the SMILES and hydrogen-suppressed molecular graph for both classes of compounds. Approximately 75% of the total experimental data set was randomly divided into training and invisible training sets, while approximately 25% was divided into calibration and validation sets. The authenticity and robustness of the developed QSPR models were also judged by the Index of Ideality of Correlation. In QSPR modelling of aliphatic compounds, the numerical values of r T r a i n i n g 2 , r V a l i d a t i o n 2 , Q T r a i n i n g 2 and Q V a l i d a t i o n 2 were in the range of 0.852-0.905, 0.815-0.894, 0.839-0.897 and 0.737-0.867, respectively. Whereas, in the QSPR modelling of phenyl-based compounds, the numerical values of r T r a i n i n g 2 , r V a l i d a t i o n 2 , Q T r a i n i n g 2 and Q V a l i d a t i o n 2 were in the range of 0.867-0.896, 0.852-0.865, 0.816-0.850 and 0.760-0.762, respectively. The structural attributes, which are promoters of l o g K e a q - increase/decrease are also extracted from the SMILES notation for mechanistic interpretation. These QSPR models can also be applied to compute the reaction rate constants of organic contaminants.

In silico enhancement of azo dye adsorption affinity for cellulose fibre through mechanistic interpretation under guidance of QSPR models using Monte Carlo method with index of ideality correlation

Azo dyes are a group of chemical moieties joined by azo (-N=N-) group with potential usefulness in different industrial applications. But these dyes are not devoid of hazardous consequence because of poor affinity for the fibre and discharge into the water stream. The chemical aspects of 72 azo dyes towards cellulose fibre in terms of their affinity by QSPR have been explored in the present work. We have employed two approaches, namely balance of correlation without IIC (TF₁) and balance of correlation with IIC (TF₂), to generate 16 QSAR models from 8 splits. The determination coefficient of calibration and validation set was found higher when the QSPR models were developed using the index of ideality correlation (IIC) parameter (TF₂). The model developed with TF₂ for split 3 was considered as a prominent model because the determination coefficient of the validation set was maximum (r ² = 0.9468). The applicability domain (AD) was also analysed based on 'statistical defect', d(A) for a SMILES attribute. The mechanistic interpretation was done by identifying the SMILES attributes responsible for the promoter of endpoint increase and promoter of endpoint decrease. These SMILES attributes were applied to design 15 new dyes with higher affinity for cellulose fibre.