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

Are we ready to combat the ecotoxicity of COVID-19 pharmaceuticals? An in silico aquatic risk assessment

To fight COVID-19 with uncountable medications and bioproducts throughout the world has taken us to another challenge of ecotoxicity. The indiscriminate usage followed by improper disposal of unused antibacterials, antivirals, antimalarials, immunomodulators, angiotensin II receptor blockers, corticosteroids, anthelmintics, anticoagulants etc. can lead us to an unimaginable ecotoxicity in the long run. A series of studies already identified active pharmaceutical ingredients (APIs) of the mentioned therapeutic classes and their metabolites in aquatic bodies as well as in wastewater treatment plants. Therefore, an initial ecotoxicity assessment of the majorly used pharmaceuticals is utmost requirement of the present time. The present in silico risk assessment study is focused on the aquatic toxicity prediction of 81 pharmaceuticals where 77 are most-used pharmaceuticals for COVID-19 throughout the world based on the literature along with one drug nirmatrelvir [PF-07321332] approved for emergency use by US-FDA and three other molecules under clinical trial. The ecotoxicity of the studied compounds were predicted based on the three aquatic species fish, algae and crustaceans employing the highest quality QSAR models available from the literature as well as using ECOSAR and QSAR Toolbox. To compare the toxicity thresholds, we have also used 4 control pharmaceuticals based on the worldwide occurrence from river, lake, STP, WWTPs, influent and effluent followed by high reported aquatic toxicity over the years as per the literature. Based on the statistical comparison, we have proposed top 3 pharmaceuticals used for the COVID-19 most toxic to the aquatic environment. The study will provide confident predictions of aquatic ecotoxicity data related to abundant use of COVID-19 drugs. The major aim of the study is to fill up the aquatic ecotoxicity data gap of major medications used for COVID-19.

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

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

Prediction reliability of QSAR models: an overview of various validation tools

The reliability of any quantitative structure-activity relationship (QSAR) model depends on multiple aspects such as the accuracy of the input dataset, selection of significant descriptors, the appropriate splitting process of the dataset, statistical tools used, and most notably on the measures of validation. Validation, the most crucial step in QSAR model development, confirms the reliability of the developed QSAR models and the acceptability of each step in the model development. The present review deals with various validation tools that involve multiple techniques that improve the model quality and robustness. The double cross-validation tool helps in building improved quality models using different combinations of the same training set in an inner cross-validation loop. This exhaustive method is also integrated for small datasets (< 40 compounds) in another tool, namely the small dataset modeler tool. The main aim of QSAR researchers is to improve prediction quality by lowering the prediction errors for the query compounds. 'Intelligent' selection of multiple models and consensus predictions integrated in the intelligent consensus predictor tool were found to be more externally predictive than individual models. Furthermore, another tool called Prediction Reliability Indicator was explained to understand the quality of predictions for a true external set. This tool uses a composite scoring technique to identify query compounds as 'good' or 'moderate' or 'bad' predictions. We have also discussed a quantitative read-across tool which predicts a chemical response based on the similarity with structural analogues. The discussed tools are freely available from https://dtclab.webs.com/software-tools or http://teqip.jdvu.ac.in/QSAR_Tools/DTCLab/ and https://sites.google.com/jadavpuruniversity.in/dtc-lab-software/home (for read-across).

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 EC₁₀ and EC₅₀ values. The QSTR models are established using CORAL software by utilizing the target function (TF₂) 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 pEC₁₀ and pEC₅₀ 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 (pEC₁₀ and pEC₅₀). The numerical value of the determination coefficient of the validation set of split 3 for the endpoint pEC₁₀ is 0.7849 and for the endpoint pEC₅₀, it is 0.8150. The structural fractions accountable for the toxicity of chemicals are also extracted. The hydrophilic attributes like 1…n…(… and S…(… Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> … 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.

In the ecotoxicological risk assessment, acute toxicity is one of the most significant criteria.

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.

Examining variable selection methods for the predictive performance of regression models and the proportion of selected variables and selected random variables

The selection of a descriptor, X, is crucial for improving the interpretation and prediction accuracy of a regression model. In this study, the prediction accuracy of models constructed using the selected X was determined and the results of variable selection, according to the number of selected X and number of selected variables that are unrelated to an objective variable, such as activities and properties (y), were investigated to evaluate the variable or feature selection methods. Variable selection methods include least absolute shrinkage and selection operator, genetic algorithm-based partial least squares, genetic algorithm-based support vector regression, and Boruta. Several regression analysis methods were used to test the prediction accuracy of the model constructed using the selected X. The characteristics of each variable selection method were analyzed using eight datasets. The results showed that even when variables unrelated to y were selected by variable selection and the number of unrelated variables was the same as the number of the original variables, a regression model with good accuracy, which ignores the influence of such noise variables, can be constructed by applying various regression analysis methods. Additionally, the variables related to y must not to be deleted. These findings provide a basis for improving the variable selection methods.

Prediction of inhibitory constants of compounds against SARS-CoV 3CLpro enzyme with 2D-QSAR model

Developing broad-spectrum anti-coronavirus drugs is greatly important, since the novel SARS-CoV-2 has rapidly become a threat to the public health and economy worldwide. SARS-CoV 3-chymotrypsin-like protease (3CLpro), as highly conserved in betacoronavirus, is a viable target for anti-SARS drugs. A quantitative structure–activity relationship (QSAR) for inhibitory constants (pKi) of 89 compounds against SARS-CoV 3CLpro enzyme was developed by using support vector machine (SVM) and genetic algorithm. The optimal SVM model (C = 90.2339 and γ = 1.19826 × 10⁻⁵) based on six molecular descriptors has determination coefficients of 0.839 for the training set (65 compounds) and 0.747 for test set (24 compounds), and rms errors of 0.435 and 0.525, respectively. These results are accurate and acceptable compared with that in other models reported, although our SVM model deals with more samples in the dada set. The SVM model could be beneficial for search of novel 3CLpro enzyme inhibitors against SARS-CoV.

A systematic review and ecological risk assessment for organic ultraviolet filters in aquatic environments

Organic ultraviolet filters (OUVFs) are used in a wide range of manufactured products including personal care (e.g. sunscreens) and plastic items. This review summarizes the available data regarding the toxic effects of OUVFs on marine and freshwater organisms and generates the predicted no-effect concentration (PNEC) values necessary for assessing ecological risk. Through a systematic search of the literature, 89 studies were identified and ecotoxicological data extracted. Collectively, these studies described toxicity testing with 39 OUVF from 10 structural classes, with derivatives of benzophenones (49%) and camphors (16%) most studied. There was a bias towards selecting freshwater species (61%), and evaluating single OUVF effects (87%) rather than OUVF mixtures. Short-term (acute) experimentation (58%) was marginally more common than long-term (chronic) testing (42%). Reproductive, developmental, genetic, and neurological toxicity were the most commonly identified effects in aquatic organism, and were associated with molecular interactions with steroid receptors, DNA, or the production of reactive oxygen species. Species sensitivity distribution and/or assessment factors were used to calculate PNECs for 22 OUVFs and the risk quotients for 12 OUVFs. When using maximum concentrations, high risk was observed for six OUVFs in marine environments (4-methylbenzylidene-camphor, octocrylene, padimate-O, benzophenone-1, and oxybenzone, ethylhexyl-4-methoxycinnamate), and for four OUVFs in freshwater environments (ethylhexyl-4-methoxycinnamate, octocrylene, avobenzone and oxybenzone). When using median concentrations, a risk to marine environments was observed for oxybenzone. The results of this review underline that there is limited knowledge of the pathological effects of OUVFs and their metabolites in aquatic environments, and this inhibits the development of informed water-quality guidelines.

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.

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.

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.

Development of a simple, interpretable and easily transferable QSAR model for quick screening antiviral databases in search of novel 3C-like protease (3CLpro) enzyme inhibitors against SARS-CoV diseases

In the context of recently emerged pandemic of COVID-19, we have performed two-dimensional quantitative structure-activity relationship (2D-QSAR) modelling using SARS-CoV-3CLpro enzyme inhibitors for the development of a multiple linear regression (MLR) based model. We have used 2D descriptors with an aim to develop an easily interpretable, transferable and reproducible model which may be used for quick prediction of SAR-CoV-3CLpro inhibitory activity for query compounds in the screening process. Based on the insights obtained from the developed 2D-QSAR model, we have identified the structural features responsible for the enhancement of the inhibitory activity against 3CLpro enzyme. Moreover, we have performed the molecular docking analysis using the most and least active molecules from the dataset to understand the molecular interactions involved in binding, and the results were then correlated with the essential structural features obtained from the 2D-QSAR model. Additionally, we have performed in silico predictions of SARS-CoV 3CLpro enzyme inhibitory activity of a total of 50,437 compounds obtained from two anti-viral drug databases (CAS COVID-19 antiviral candidate compound database and another recently reported list of prioritized compounds from the ZINC15 database) using the developed model and provided prioritized compounds for experimental detection of their performance for SARS-CoV 3CLpro enzyme inhibition.

Application of QSARs in identification of mutagenicity mechanisms of nitro and amino aromatic compounds against Salmonella typhimurium species

In an attempt to describe the underlying causes of mutagenicity mainly due to organic chemicals, quantitative structure-activity relationship (QSAR) models have been developed using two different Salmonella typhimurium mutagenicity endpoints with or without presence of liver metabolic microsomal enzymes (S9) namely TA98-S9 and TA98 + S9. The models were developed using simple 2D variables having definite physicochemical meaning calculated from Dragon, SiRMS, and PaDEL-descriptor software tools. Stepwise regression followed by partial least squares (PLS) regression was used in model development following the strict OECD guidelines for QSAR model development and validation. The models were validated using coefficient of determination R², cross-validation coefficient Q²_LOO (leave one out) while the test set predictions were analyzed using Q²F₁ (coefficient of determination for the test set). Several other internationally accepted validation metrics like MAE_95%train, average r_m(LOO)² and Δr_m(LOO)² (for the training set) were used to check model robustness while predictive efficiency was evaluated using MAE_95%test, average r_m² and Δr_m² (for the test set). The scope of predictions was defined by applicability domain analysis using the DModX approach, a recommended tool for PLS models. The major contributing features related to mutagenicity include lipophilicity, electronegativity, branching and unsaturation, etc. The present manuscript is the first attempt to undertake modeling of two different endpoints (TA98-S9 and TA98 + S9) in order to explore major contributing molecular features linked directly or indirectly to mutagenicity.