Ecotoxicological modelling of cosmetics for aquatic organisms: A QSTR approach

QSTR Models in Dioxins and Dioxin-like Compounds Provide Insights into Gene Expression Dysregulation

Polychlorinated dibenzo-p-dioxins and polychlorinated dibenzo-p-furans (PCDD/Fs) are a group of organic chemicals containing three-ring structures that can be substituted with one to eight chlorine atoms, leading to 75 dioxin and 135 furan congeners. As endocrine-disrupting chemicals (EDCs), they can alter physiological processes causing a number of disorders. In this study, quantitative structure-toxicity relationship (QSTR) studies were used to determine the correlations between the PCDD/Fs' molecular structures and various toxicity endpoints. Strong QSTR models, with the coefficients of determination (r2) values greater than 0.95 and ANOVA p-values less than 0.0001 were established between molecular descriptors and the endpoints of bioconcentration, fathead minnow LC50, and Daphnia magna LC50. The ability of PCDD/Fs to bind to several nuclear receptors was investigated via molecular docking studies. The results show comparable, and in some instances better, binding affinities of PCDD/Fs toward the receptors relative to their natural agonistic and antagonistic ligands, signifying possible interference with the receptors' natural biological activities. These studies were accompanied by the molecular dynamics simulations of the top-binding PCDD/Fs to show changes in the receptor-ligand complexes during binding and provide insights into these compounds' ability to interfere with transcription and thereby modify gene expression. This introspection of PCDD/Fs at the molecular level provides a deeper understanding of these compounds' toxicity and opens avenues for future studies.

Development of an effective QSAR-based hazard threshold prediction model for the ecological risk assessment of aromatic hydrocarbon compounds

The insufficient hazard thresholds of specific individual aromatic hydrocarbon compounds (AHCs) with diverse structures limit their ecological risk assessment. Thus, herein, quantitative structure-activity relationship (QSAR) models for estimating the hazard threshold of AHCs were developed based on the hazardous concentration for 5% of species (HC5) determined using the optimal species sensitivity distribution models and on the molecular descriptors calculated via the PADEL software and ORCA software. Results revealed that the optimal QSAR model, which involved eight descriptors, namely, Zagreb, GATS2m, VR3_Dzs, AATSC2s, GATS2c, ATSC2i, ω, and Vm, displayed excellent performance, as reflected by an optimal goodness of fit (R2adj = 0.918), robustness (Q2LOO = 0.869), and external prediction ability (Q2F1 = 0.760, Q2F2 = 0.782, and Q2F3 = 0.774). The hazard thresholds estimated using the optimal QSAR model were approximately close to the published water quality criteria developed by different countries and regions. The quantitative structure-toxicity relationship demonstrated that the molecular descriptors associated with electrophilicity and topological and electrotopological properties were important factors that affected the risks of AHCs. A new and reliable approach to estimate the hazard threshold of ecological risk assessment for various aromatic hydrocarbon pollutants was provided in this study, which can be widely popularised to similar contaminants with diverse structures.

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.

Toxicity prediction of 1,2,4-triazoles compounds by QSTR and interspecies QSTTR models

1,2,4-triazole derivatives exhibit various biological activities, including antibacterial and antifungal properties. On the other hand, these chemicals may have unique cumulative and harmful effects on living organisms. The goal of this work is to use quantitative structure-toxicity relationship (QSTR) and interspecies quantitative toxicity-toxicity relationship (iQSTTR) models to predict the acute toxicity of 1,2,4-triazole derivatives. The QSTR models were generated by multiple linear regression (MLR) following the OECD recommendations for QSAR model development and validation. The iQSTTR models were constructed using data on acute oral toxicity in rats and mice, as well as the 2D descriptor. The application domain (AD) analysis was used to identify model outliers and determine if the forecast was credible. Six QSTR models were successfully constructed in rats and mice using various delivery methods, and the scatter plots demonstrated excellent consistency across training and test sets. According to external and internal validation criteria, all six QSTR models may be broadly accepted; however, the orally administered mice model was the optimum one among the six species. Several chemicals with leverage values above the requirements were identified as response or structural outliers in the training sets for six QSTR and two iQSTTR models. All outliers, however, fell slightly outside the threshold or had low prediction errors, which may have had little impact on the capacity to forecast and were therefore preserved in the final models. In fact, neither the QSTR nor the iQSTTR test sets contained any response outliers. Additionally, all external and internal validation results for the iQSTTR models were approved, with the iQSTTR models outperforming the comparable QSTR models, which are deemed more dependable. The QSTR and iQSTTR models performed well in predicting toxicity using test sets, which would be beneficial in evaluating and synthesizing newly discovered 1,2,4-triazoles derivatives with low toxicity and environmental hazard.

Ecotoxicological prediction of organic chemicals toward Pseudokirchneriella subcapitata by Monte Carlo approach

In the ecotoxicological risk assessment, acute toxicity is one of the most significant criteria. Green alga Pseudokirchneriella subcapitata has been used for ecotoxicological studies to assess the toxicity of different toxic chemicals in freshwater. Quantitative Structure Activity Relationships (QSAR) are mathematical models to relate chemical structure and activity/physicochemical properties of chemicals quantitatively. Herein, Quantitative Structure Toxicity Relationship (QSTR) modeling is applied to assess the toxicity of a data set of 334 different chemicals on Pseudokirchneriella subcapitata, in terms of EC10 and EC50 values. The QSTR models are established using CORAL software by utilizing the target function (TF2) with the index of ideality of correlation (IIC). A hybrid optimal descriptor computed from SMILES and molecular hydrogen-suppressed graphs (HSG) is employed to construct QSTR models. The results of various statistical parameters of the QSTR model developed for pEC10 and pEC50 range from excellent to good and are in line with the standard parameters. The models prepared with IIC for Split 3 are chosen as the best model for both endpoints (pEC10 and pEC50). The numerical value of the determination coefficient of the validation set of split 3 for the endpoint pEC10 is 0.7849 and for the endpoint pEC50, it is 0.8150. The structural fractions accountable for the toxicity of chemicals are also extracted. The hydrophilic attributes like 1…n…(… and S…(…[double bond, length as m-dash]… exert positive contributions to controlling the aquatic toxicity and reducing algal toxicity, whereas attributes such as c…c…c…, C…C…C… enhance lipophilicity of the molecules and consequently enhance algal toxicity.

Application of QSAR for the identification of key molecular fragments and reliable predictions of effects of textile dyes on growth rate and biomass values of Raphidocelis subcapitata

The current quantitative structure-activity relationship (QSAR) study seeks to explore the underlying causes of fluctuations in growth rate and biomass of microalgae mainly due to textile dyes. The derived QSAR models cover two endpoints: ErC50 (growth rate) and EbC50 (biomass) of Raphidocelis subcapitata. In order to extract the structural features involved, multiple PLS (partial least squares) models have been developed with easy to interpret and uncomplicated 2D descriptors having proper physico-chemical meaning. These descriptors were calculated from Dragon, SiRMS, and PaDEL-descriptor software. Then, the models were developed initially using stepwise regression followed by partial least squares (PLS) regression, and the model development procedure for both the endpoints (ErC50 and EbC50) followed the stringent Organization for Economic Cooperation and Development (OECD) rules. Later on, the model validation was carried out with statistically significant and internationally accepted metrics (both internally and externally) in both the cases. Next, we have used the "Intelligent Consensus Predictor" tool (available from http://teqip.jdvu.ac.in/QSAR_Tools/DTCLab/) to test the prediction quality with an "intelligent" approach to select multiple models. The estimated prediction quality for the appropriate test sets reveals that the consensus models (CM) surpass the quality shown by individual models (IM) for both the endpoints (ErC50 and EbC50). Finally, the developed models were able to identify the major contributing features (hydrophobic units, unsaturation, saturation, electronegativity, branched atoms and charged fragments) related to aquatic toxicity of textile dyes.

In silico modelling of acute toxicity of 1, 2, 4-triazole antifungal agents towards zebrafish (Danio rerio) embryos: Application of the Small Dataset Modeller tool

Nowadays, there is a widespread use of triazole antifungal agents to kill broad classes of fungi in farming lands and to protect herbs, fruits and grains. These agents further deposit into the aquatic systems causing toxicity to the living aquatic creatures, which can then affect human beings. Considering this issue, risk assessment of these toxic chemicals is a very essential task. Due to the inadequate experimental data on acute toxicity of antifungal agents containing the 1, 2, 4-triazole ring, higher testing costs along with the regulatory restrictions and the international regulations to lessen animal testing emphasize on in silico techniques such as quantitative structure-activity relationship (QSAR) studies. The application of QSAR modelling has created an easier avenue to predict activity/property/toxicity of newly synthesized compounds. In the present study, we have used 23 antifungal agents containing the 1, 2, 4-triazole ring to develop 2D-QSAR models and explored their structural attributes crucial for acute toxicity towards embryonic phase of zebrafish (Danio rerio). Here, we have employed simple 2D descriptors to develop the QSAR models. The models were evolved by executing the Small Dataset Modeller tool (https://dtclab.webs.com/software-tools), and the validation of the models was achieved by employing different precise validation principles. The statistical validation metrics confirm that built models are robust, useful and well predictive to forecast the acute toxicity of new compounds.

Robust classification-based molecular modelling of diverse chemical entities as potential SARS-CoV-2 3CLpro inhibitors: theoretical justification in light of experimental evidences

COVID-19 is the most unanticipated incidence of 2020 affecting the human population worldwide. Currently, it is utmost important to produce novel small molecule anti-SARS-CoV-2 drugs urgently that can save human lives globally. Based on the earlier SARS-CoV and MERS-CoV infection along with the general characters of coronaviral replication, a number of drug molecules have been proposed. However, one of the major limitations is the lack of experimental observations with different drug molecules. In this article, 70 diverse chemicals having experimental SARS-CoV-2 3CL^proinhibitory activity were accounted for robust classification-based QSAR analysis statistically validated with 4 different methodologies to recognize the crucial structural features responsible for imparting the activity. Results obtained from all these methodologies supported and validated each other. Important observations obtained from these analyses were also justified with the ligand-bound crystal structure of SARS-CoV-2 3CL^pro enzyme. Our results suggest that molecules should contain a 2-oxopyrrolidine scaffold as well as a methylene (hydroxy) sulphonic acid warhead in proper orientation to achieve higher inhibitory potency against SARS-CoV-2 3CL^pro. Outcomes of our study may be able to design and discover highly effective SARS-CoV-2 3CL^pro inhibitors as potential anticoronaviral therapy to crusade against COVID-19.

A critical review on the detection, occurrence, fate, toxicity, and removal of cannabinoids in the water system and the environment

Cannabinoids are a group of organic compounds found in cannabis. Δ⁹-tetrahydrocannabinol (THC) and cannabidiol (CBD), the two major constituents of cannabinoids, and their metabolites are contaminants of emerging concern due to the limited information on their environmental impacts. As well, their releases to the water systems and environment are expected to increase due to recent legalization. Solid-phase extraction is the most common technique for the extraction and pre-concentration of cannabinoids in water samples as well as a clean-up step after the extraction of cannabinoids from solid samples. Liquid chromatography coupled with mass spectrometry is the most common technique used for the analysis of cannabinoids. THC and its metabolites have been detected in wastewater, surface water, and drinking water. In particular, 11-nor-9-carboxy-Δ⁹-tetrahydrocannabinol (THC-COOH) has been detected at concentrations up to 2590 and 169 ng L^-1 in untreated and treated wastewater, respectively, 79.9 ng L^-1 in surface water, and 1 ng L^-1 in drinking water. High removal of cannabinoids has been observed in wastewater treatment plants; this is likely a result of adsorption due to the low aqueous solubility of cannabinoids. Based on the estrogenicity and cytotoxicity studies and modelling, it has been predicted that THC and THC-COOH pose moderate risk for adverse impact on the environment. While chlorination and photo-oxidation have been shown to be effective in the removal of THC-COOH, they also produce by-products that are potentially more toxic than regulated disinfection by-products. The potential of indirect exposure to cannabinoids and their metabolites through recreational water is of great interest. As cannabinoids and especially their by-products may have adverse impacts on the environment and public health, more studies on their occurrence in various types of water and environmental systems, as well as on their environmental toxicity, would be required to accurately assess their impact on the environment and public health.

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.

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.

Quantitative structure-toxicity relationships of organic chemicals against Pseudokirchneriella subcapitata

Predicting the toxicity of organic toxicants to aquatic life through chemometric approach is challenging area. In this paper, a six-descriptor quantitative structure-activity/toxicity relationship (QSAR/QSTR) model was successfully developed for the toxicity pEC₁₀ of organic chemicals against Pseudokirchneriella subcapitata, by applying support vector machine (SVM) together with genetic algorithm. A sufficiently large data set consisting of 334 organic chemicals was randomly divided into a training set (167 compounds) and a test set (167 compounds) with a ratio of 1:1. The optimal SVM model possesses coefficient of determination R² of 0.76 and mean absolute error (MAE) of 0.60 for the training set and R² of 0.75 and MAE of 0.61 for the test set. Compared with other models reported in the literature, our SVM model for the toxicity pEC₁₀ shows significant statistical quality and satisfactory predictive ability, although it has fewer molecular descriptors and more samples in the test set. A QSTR model for pEC₅₀ of organic chemicals against Pseudokirchneriella subcapitata was also developed with the same subsets and molecular descriptors.

Reliable and representative in silico predictions of freshwater ecotoxicological hazardous concentrations

A reliable quantification of the potential effects of chemicals on freshwater ecosystems requires ecotoxicological response data for a large set of species which is typically not available in practice. In this study, we propose a method to estimate hazardous concentrations (HCs) of chemicals on freshwater ecosystems by combining two in silico approaches: quantitative structure activity relationships (QSARs) and interspecies correlation estimation (ICE) models. We illustrate the principle of our QSAR-ICE method by quantifying the HCs of 51 chemicals at which 50% and 5% of all species are exposed above the concentration causing acute effects. We assessed the bias of the HCs, defined as the ratio of the HC based on measured ecotoxicity data and the HC based on in silico data, as well as the statistical uncertainty, defined as the ratio of the 95th and 5th percentile of the HC. Our QSAR-ICE method resulted in a bias that was comparable to the use of measured data for three species, as commonly used in effect assessments: the average bias of the QSAR-ICE HC50 was 1.2 and of the HC5 2.3 compared to 1.2 when measured data for three species were used for both HCs. We also found that extreme statistical uncertainties (>10⁵) are commonly avoided in the HCs derived with the QSAR-ICE method compared to the use of three measurements with statistical uncertainties up to 10¹². We demonstrated the applicability of our QSAR-ICE approach by deriving HC50s for 1,223 out of the 3,077 organic chemicals of the USEtox database. We conclude that our QSAR-ICE method can be used to determine HCs without the need for additional in vivo testing to help prioritise which chemicals with no or few ecotoxicity data require more thorough assessment.

QSAR and molecular docking studies of lethal factor protease inhibitors against Bacillus anthracis

Bacillus anthracis is considered as a biological warfare agent because it is the causative agent of the serious infectious anthrax disease. Delay in treatment leads to lethal factor-mediated toxaemia which is very critical due to lack of therapeutic options. Consequently, attempts have been made to discover potent lethal factor (LF) protease inhibitors such as small-molecule synthetic 2-thio-1,3-thiazolidine-4-one (rhodanine) compounds. But computed descriptor-based quantitative structure-activity relationship (QSAR) and drug design studies on such aspect are poorly represented. Therefore, an attempt was made for developing QSAR models using structural descriptors for 1,3-thiazolidine-4-one compounds. The models were developed on a series of 49 LF protease inhibitors using the combination of constitutional, functional group, atom-centred fragment and molecular property descriptors. The best QSAR model included four variables, namely, C-040, nR05, GVWAI-80 and ALOGP that correlated well with the anti-LF protease activity with a good correlation coefficient (r = 0.870) of good statistical significance (F_{4, 29} = 14.09 (α = 0.001) F_{4, 29} = 6.19). This model was also validated and explained 58.1% of variances of the Bacillus anthracis inhibitory activities of the studied compounds with r²_pred = 0.710 which denotes external predictability. Finally, molecular docking was carried out to predict the mode of binding of some highly active congeneric compounds. It was shown that VAL 1403 is an important residue for phenyl ring. TYR 1456 and HIS 1418 are responsible for interaction with the rhodanine nucleus. Therefore, these residues are considered responsible for the inhibition of LF protease anthrax and can predict significant dimension of essential structural features of these inhibitors to evaluate, screen and help priorities of the synthesis of the candidates against anthrax bioterrorism.

Ecotoxicological QSAR modelling of organic chemicals against Pseudokirchneriella subcapitata using consensus predictions approach

The present study provides robust consensus quantitative structure-activity relationship (QSAR) models developed from 334 organic chemicals covering a wide chemical domain for the prediction of effective concentrations of chemicals for 50% and 10% inhibition of algal growth. Only 2D descriptors with definite physicochemical meaning were employed for QSAR model building, whereas development, validation and interpretation were achieved following the strict Organization for Economic Co-operation and Development (OECD) recommended guidelines. Genetic algorithm along with stepwise approach was used in feature selection while the final QSAR models were derived using partial least squares regression technique. The applicability domain of the developed models was also checked. The obtained consensus models were then used to predict 64 organic chemicals having no definite observed responses while the confidence of predictions was checked by the 'prediction reliability indicator' tool. The developed models should be applicable for data gap filling in case of new or untested organic chemicals provided they fall within the domain of the model and can also be implemented to design safer alternatives to the environment.

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

In the recent years, ecotoxicological hazard potential of biocidal products has been receiving increasing attention in the industries and regulatory agencies. Biocides/pesticides are currently one of the most studied groups of compounds, and their registration cannot be done without the empirical toxicity information. In view of limited experimental data available for these compounds, we have developed Quantitative Structure-Activity Relationship (QSAR) models for the toxicity of biocides to fish and Daphnia magna following principles of QSAR modeling recommended by the OECD (Organization for Economic Cooperation and Development). The models were developed using simple and interpretable 2D descriptors and validated using stringent tests. Both models showed encouraging statistical quality in terms of determination coefficient R² (0.800 and 0.648), cross-validated leave-one-out Q² (0.760 and 0.602) and predictive R²_pred or Q²_ext (0.875 and 0.817) for fish (n_Training = 66, n_Test = 22) and Daphnia magna (n_Training = 100, n_Test = 33) toxicity datasets, respectively. These models should be applicable for data gap filling in case of new or untested biocidal compounds falling within the applicability domain of the models. In general, the models indicate that the toxicity increases with lipophilicity and decreases with polarity, branching and unsaturation. We have also developed interspecies toxicity models for biocides using the daphnia and fish toxicity data and used the models for data gap filling.

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

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

Chemometric modeling of Daphnia magna toxicity of agrochemicals

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

Ecotoxicological QSAR modeling of endocrine disruptor chemicals

This study reports highly robust externally predictive quantitative structure-toxicity relationship (QSTR) and interspecies quantitative structure-toxicity-toxicity (i-QSTTR) models developed using toxicity data of endocrine disruptor chemicals (EDCs) towards 14 different species falling in four different trophic levels. Genetic algorithm followed by Partial Least Squares (PLS) regression was used in model development following the strict OECD guidelines. The models were developed using 2D descriptors having definite physicochemical meaning and validated by several internationally accepted validation metrics. The scope of predictions was defined by estimating applicability domain of the models. Presence of halogens, sulfur and phosphorus in the molecules greatly influenced the toxicity of EDCs as suggested by continuous repetition of 2D atom pair descriptors. Lipophilic contributions as calculated by logP terms (mainly ALOGP2 and XlogP) were the second most important feature controlling the EDC hazards. Hydrophilic moiety such as functionalities like esters, aliphatic ethers, branching and higher oxygen content reduced the EDC toxicity. Interspecies models were employed in data gap filling following the hierarchy of different species. The reliability of predictions was calculated by the "prediction reliability indicator" tool.

Consensus QSAR modeling of toxicity of pharmaceuticals to different aquatic organisms: Ranking and prioritization of the DrugBank database compounds

In the present work, quantitative structure-activity relationship (QSAR) models have been developed for ecotoxicity of pharmaceuticals on four different aquatic species namely Pseudokirchneriella subcapitata, Daphnia magna, Oncorhynchus mykiss and Pimephales promelas using genetic algorithm (GA) for feature selection followed by Partial Least Squares regression technique according to the Organization for Economic Co-operation and Development (OECD) guidelines. Double cross-validation methodology was employed for selecting suitable models. Only 2D descriptors were used for capturing chemical information and model building, whereas validation of the models was performed by considering various stringent internal and external validation metrics. Interestingly, models could be developed even without using any LogP terms in contrary to the usual dependence of toxicity on lipophilicity. However, the current manuscript proposes highly robust and more predictive models employing computed logP descriptors. The applicability domain study was performed in order to set a predefined chemical zone of applicability for the obtained QSAR models, and the test compounds falling outside the domain were not taken for further analysis while making a prioritized list. An additional comparison was made with ECOSAR, an online expert system for toxicity prediction of organic pollutants, in order to prove predictability of the obtained models. The obtained robust consensus models were utilized to predict the toxicity of a large dataset of approximately 9300 drug-like molecules in order to prioritize the existing drug-like substances in accordance to their acute predicted aquatic toxicities following a scaling technique. Finally, prioritized lists of 500 most toxic chemicals obtained by respective consensus models and those predicted from ECOSAR tool have been reported.