Comparative Quantitative Structure–Activity–Activity Relationships for Toxicity to Tetrahymena pyriformis and Pimephales promelas

Modelling quantitative structure activity-activity relationships (QSAARs): auto-pass-pass, a new approach to fill data gaps in environmental risk assessment under the REACH regulation

Reviewing the toxicological literature for over the past decades, the key elements of QSAR modelling have been the mechanisms of toxic action and chemical classes. As a result, it is often hard to design an acceptable single model for a particular endpoint without clustering compounds. The main aim here was to develop a Pass-Pass Quantitative Structure-Activity-Activity Relationship (PP QSAAR) model for direct prediction of the toxicity of a larger set of compounds, combing the application of an already predicted model for another species, and molecular descriptors. We investigated a large acute toxicity data set of five aquatic organisms, fish, Daphnia magna, and algae from the VEGA-Hub, as well as Tetrahymena pyriformis and Vibrio fischeri. The statistical quality of the models encouraged us to consider this alternative for the prediction of toxicity using interspecies extrapolation QSAAR models without regard to the toxicity mechanism or chemical class. In the case of algae, the use of activity values from a second species did not improve the models. This can be attributed to the weak interspecies relationships, due to different aquatic toxicity mechanisms in species.

'Ideal correlations' for the predictive toxicity to Tetrahymena pyriformis

OBJECTIVES

Predictive models for toxicity to Tetrahymena pyriformis are an important component of natural sciences. The present study aims to build up a predictive model for the endpoint using the so-called index of ideality of correlation (IIC). Besides, the comparison of the predictive potential of these models with the predictive potential of models suggested in the literature is the task of the present study.

METHODS

The Monte Carlo technique is a tool to build up the predictive model applied in this study. The molecular structure is represented via a simplified molecular input-line entry system (SMILES). The IIC is a statistical characteristic sensitive to both the correlation coefficient and mean absolute error. Applying of the IIC to build up quantitative structure-activity relationships (QSARs) for the toxicity to Tetrahymena pyriformis improves the predictive potential of those models for random splits into the training set and the validation set. The calculation was carried out with CORAL software (http://www.insilico.eu/coral).

RESULTS

The statistical quality of the suggested models is incredibly good for the external validation set, but the statistical quality of the models for the training set is modest. This is the paradox of ideal correlation, which is obtained with applying the IIC.

CONCLUSIONS

The Monte Carlo technique is a convenient and reliable way to build up a predictive model for toxicity to Tetrahymena pyriformis. The IIC is a useful statistical criterion for building up predictive models as well as for the assessment of their statistical quality.

Development of QSAAR and QAAR models for predicting fish early-life stage toxicity with a focus on industrial chemicals

We developed models for predicting fish early-life stage (ELS) toxicities oriented to industrial chemicals. The training set was constructed without data from the Office of Pesticide Programs Pesticide Ecotoxicity Database, the main source for the pesticide-biased training set used in our previous work (SAR QSAR Environ. Res. 29:9, 725-742). In addition to the descriptors from the previous study, we also used water solubility to develop the new models, which were evaluated against the test set used in our previous study so that we could focus on the effects of the different training set and the additional descriptor. The statistics for the new models were hardly better than those for the previous models, which suggests, contrary to our expectations, that pesticide-biased data can successfully be used to develop models for predicting the fish ELS toxicities oriented to industrial chemicals. Acute Daphnia magna toxicity was important for the predictive QSAARs in both studies. A distance-based method for defining the applicability domains indicated that water solubility was a key indicator for detecting underestimated chemicals. The comparison of fish ELS toxicities for chemicals presented in different literatures revealed the uncertainty of the experimental data, which may lead to the low predictivity.

Development of models to predict fish early-life stage toxicity from acute Daphnia magna toxicity$

Herein, we propose models for predicting fish early-life stage (ELS) toxicity from acute Daphnia magna toxicity and various molecular descriptors. Specifically, eight models were developed with fathead minnow (Pimephales promelas) data and were validated against Japanese medaka (Oryzias latipes) data because the quantity of available Japanese medaka data is much smaller than the quantity of fathead minnow data. The training data set for the models consisted of ELS fathead minnow toxicity data for 77 chemicals; data for 67 of the 77 chemicals originated from the OPP Pesticide Ecotoxicity Database of the US Environmental Protection Agency. The training data were biased toward pesticides. A simple quantitative activity-activity relationship (QAAR) model based on the correlation between fish ELS and acute Daphnia magna toxicities showed good predictivity for the chemicals in the external validation data set relative to the predictivities of the other models in this study. However, goodness-of-fit and robustness were better for quantitative structure-activity-activity relationship (QSAAR) models that included molecular descriptors (such as pesticide-related atoms and substructures as well as molecular weight and three-dimensional-structure-based parameters). A battery approach involving the use of both the QAAR and the QSAARs might enhance the reliability of the estimated values and prevent underestimates.

Aquatic ecotoxicological models and their applicability in Arctic regions

Dose-response modeling is one of the most important steps of ecological risk assessment. It requires concentration-effects relationships for the species under consideration. There are very limited studies and experimental data available for the Arctic aquatic species. Lack of toxicity data hinders obtaining dose-response relationships for lethal (LC50 values), sub-lethal and carcinogenic effects. Gaps in toxicity data could be filled using a variety of in-silico ecotoxicological methods. This paper reviews the suitability of such methods for the Arctic scenario. Mechanistic approaches like toxicokinetic and toxicodynamic analysis are found to be better suited for interspecies extrapolation than statistical methods, such as Quantitative Structure-Activity Relationships/Quantitative Structure Activity-Activity Relationship, ICE, and other empirical models, such as Haber's law and Ostwald's equation. A novel approach is proposed where the effects of the toxicant exposure are quantified based on the probability of cellular damage and metabolites interactions. This approach recommends modeling cellular damage using a toxicodynamic model and physiology or metabolites interactions using a toxicokinetic model. Together, these models provide more reliable estimates of toxicity in the Arctic aquatic species, which will assist in conducting ecological risk assessment of Arctic environment.

(Q)SARs to predict environmental toxicities: current status and future needs

The current state of the art of (Quantitative) Structure-Activity Relationships ((Q)SARs) to predict environmental toxicity is assessed along with recommendations to develop these models further. The acute toxicity of compounds acting by the non-polar narcotic mechanism of action can be well predicted, however other approaches, including read-across, may be required for compounds acting by specific mechanisms of action. The chronic toxicity of compounds to environmental species is more difficult to predict from (Q)SARs, with robust data sets and more mechanistic information required. In addition, the toxicity of mixtures is little addressed by (Q)SAR approaches. Developments in environmental toxicology including Adverse Outcome Pathways (AOPs) and omics responses should be utilised to develop better, more mechanistically relevant, (Q)SAR models.

Ecotoxicity interspecies QAAR models from Daphnia toxicity of pharmaceuticals and personal care products

Pharmaceutical and Personal Care Products (PPCPs) became a class of contaminants of emerging concern because are ubiquitously detected in surface water and soil, where they can affect wildlife. Ecotoxicological data are only available for a few PPCPs, thus modelling approaches are essential tools to maximize the information contained in the existing data. In silico methods may be helpful in filling data gaps for the toxicity of PPCPs towards various ecological indicator organisms. The good correlation between toxicity toward Daphnia magna and those on two fish species (Pimephales promelas and Oncorhynchus mykiss), improved by the addition of one theoretical molecular descriptor, allowed us to develop predictive models to investigate the relationship between toxicities in different species. The aim of this work is to propose quantitative activity-activity relationship (QAAR) models, developed in QSARINS and validated for their external predictivity. Such models can be used to predict the toxicity of PPCPs to a particular species using available experimental toxicity data from a different species, thus reducing the tests on organisms of higher trophic level. Similarly, good QAAR models, implemented by molecular descriptors to improve the quality, are proposed here for fish interspecies. We also comment on the relevance of autocorrelation descriptors in improving all studied interspecies correlations.

Interspecies quantitative structure-activity relationships (QSARs) for eco-toxicity screening of chemicals: the role of physicochemical properties

In addition to molecular structure profiles, descriptors based on physicochemical properties are useful for explaining the eco-toxicities of chemicals. In a previous study we reported that a criterion based on the difference between the partition coefficient (log POW) and distribution coefficient (log D) values of chemicals enabled us to identify aromatic amines and phenols for which interspecies relationships with strong correlations could be developed for fish-daphnid and algal-daphnid toxicities. The chemicals that met the log D-based criterion were expected to have similar toxicity mechanisms (related to membrane penetration). Here, we investigated the applicability of log D-based criteria to the eco-toxicity of other kinds of chemicals, including aliphatic compounds. At pH 10, use of a log POW - log D > 0 criterion and omission of outliers resulted in the selection of more than 100 chemicals whose acute fish toxicities or algal growth inhibition toxicities were almost equal to their acute daphnid toxicities. The advantage of log D-based criteria is that they allow for simple, rapid screening and prioritizing of chemicals. However, inorganic molecules and chemicals containing certain structural elements cannot be evaluated, because calculated log D values are unavailable.

Interspecies quantitative structure-activity-activity relationships (QSAARs) for prediction of acute aquatic toxicity of aromatic amines and phenols

We propose interspecies quantitative structure-activity-activity relationships (QSAARs), that is, QSARs with descriptors, to estimate species-specific acute aquatic toxicity. Using training datasets consisting of more than 100 aromatic amines and phenols, we found that the descriptors that predicted acute toxicities to fish (Oryzias latipes) and algae were daphnia toxicity, molecular weight (an indicator of molecular size and uptake) and selected indicator variables that discriminated between the absence or presence of various substructures. Molecular weight and the selected indicator variables improved the goodness-of-fit of the fish and algae toxicity prediction models. External validations of the QSAARs proved that algae toxicity could be predicted within 1.0 log unit and revealed structural profiles of outlier chemicals with respect to fish toxicity. In addition, applicability domains based on leverage values provided structural alerts for the predicted fish toxicity of chemicals with more than one hydroxyl or amino group attached to an aromatic ring, but not for fluoroanilines, which were not included in the training dataset. Although these simple QSAARs have limitations, their applicability is defined so clearly that they may be practical for screening chemicals with molecular weights of ≤364.9.

Predicting aquatic toxicities of benzene derivatives in multiple test species using local, global and interspecies QSTR modeling approaches

A flow diagram showing QSTR modeling strategy for aquatic toxicity prediction of benzene derivatives in multiple test species.

Determining chemical reactivity driving biological activity from SMILES transformations: the bonding mechanism of anti-HIV pyrimidines

Assessing the molecular mechanism of a chemical-biological interaction and bonding stands as the ultimate goal of any modern quantitative structure-activity relationship (QSAR) study. To this end the present work employs the main chemical reactivity structural descriptors (electronegativity, chemical hardness, chemical power, electrophilicity) to unfold the variational QSAR though their min-max correspondence principles as applied to the Simplified Molecular Input Line Entry System (SMILES) transformation of selected uracil derivatives with anti-HIV potential with the aim of establishing the main stages whereby the given compounds may inhibit HIV infection. The bonding can be completely described by explicitly considering by means of basic indices and chemical reactivity principles two forms of SMILES structures of the pyrimidines, the Longest SMILES Molecular Chain (LoSMoC) and the Branching SMILES (BraS), respectively, as the effective forms involved in the anti-HIV activity mechanism and according to the present work, also necessary intermediates in molecular pathways targeting/docking biological sites of interest.

Quantitative relationship between rate constants and molecular structure descriptors for the gas phase hydrogen abstraction reactions

The abstraction of hydrogen by general radicals has a wide role in environmental and also in technological processes because it results in reactive free radicals that play a vital role in atmospheric chemistry and also in biochemical processes. In addition to experimental studies, the theoretical modelling of this elementary reaction has been important for understanding and predicting respective rate constants. In this paper, molecular descriptors in the context of a QSAR approach are used to codify the relationship between molecular structure and rate constants. Unique experimental data is collected from the literature for the reaction R(i)• + R(j)H → R(i)H + R(j)•, where R(i)• = H• and R(j)• are diverse radicals. The four-parameter QSAR model (n = 34, r(2) = 0.81, r(2)(CV) = 0.74, r(2)(scr) = 0.12, s(2) = 0.19) is presented for the bimolecular rate constants, accompanied with model diagnostics and analysis of descriptors in the model.

Prediction of acute aquatic toxicity in Tetrahymena pyriformis--'Eco-Derek', a knowledge-based system approach

A 'proof-of-concept' version of a software tool for making transparent predictions of acute aquatic toxicity has been developed. It is primarily limited to semi-quantitative predictions in one species, the ciliated protozoan, Tetrahymena pyriformis. A freely available system, 'Eco-Derek', was derived by adapting a well-established, knowledge-based structure-activity and reasoning platform (Derek for Windows, Lhasa Limited). The Derek reasoning code was modified to express potency rather than confidence. Structure-activity relationship (SAR) development utilised a curated version of a published dataset, supplemented with the CADASTER Challenge datasets. Forty-five structural alerts were produced. The dependence on log P was examined for each alert and entered into the system as qualitative reasoning rules specifying the predicted potency as Very Low, Low, Moderate, High or Very High. Evaluation studies showed: (a) moderate accuracy for the training set but low accuracy for an external test set; (b) non-linearity in the toxicity-log P relationship for chemicals without identified structural alerts; (c) insufficient differentiation of substituent effects in some of the reactivity-based structural alerts resulting in too few chemicals predicted with Very High toxicity; and (d) the need for additional structural alerts covering polar narcosis and less common reactive or metabolically activated chemical functionality.

Biological Read-Across: Mechanistically-Based Species–Species and Endpoint–Endpoint Extrapolations

This chapter describes the development and use of relationships between toxicity data, for the same chemicals, derived from different species. These relationships provide simple models to predict toxicity of, potentially, a higher species from that of a lower species. Approaches to the formation of these models are described, notably inter-species relationships, quantitative structure-activity-activity relationships and prediction models. It is noted that the best extrapolations are for closely related species i.e. within taxa. In addition, forming groups or categories of compounds according to common mechanisms of toxic action improves the correlation for extrapolations from lower to higher species. A freely available software package, Web-ICE, is introduced as a suitable tool to apply these methods.

Estimating the toxicities of organic chemicals in activated sludge process

The experimental logEC50 toxicity values of 104 compounds causing bioluminescent repression of the bacterium strain Pseudomonas isolated from an industrial wastewater were studied. Using the Best Multilinear Regression method implemented in CODESSA PRO, models with up to 8 theoretical descriptors were obtained. Utilizing a rigorous descriptor selection and validation procedure a reliable QSAR model with four parameters was selected as best. The proposed model emphasizes the importance of the halogen atoms presented in each compound, the possibility of H-bond formation and the flexibility and degree of branching of the molecules. As pointed out by many researchers, the contribution of the octanol-water partition coefficient to the explanation of the toxicity effect was also found to be significant. In addition, the model currently proposed was compared to those reported earlier and its advantages were discussed in detail.

QSAR with quantum topological molecular similarity indices: toxicity of aromatic aldehydes to Tetrahymena pyriformis

Extensive production and utilization of aromatic aldehydes and their derivatives without proper certification is alarming with regard to environmental safety. This concern motivated our construction of predictive quantitative structure-activity relationship (QSAR) models for the toxicity of aldehydes to the ecologically important species Tetrahymena pyriformis. Quantum topological molecular similarity (QTMS) descriptors, along with the lipid-water partition coefficient (log K(o/w)), were used as predictor variables. The QTMS descriptors were calculated at different levels of theory including AM1, HF/3-21G(d), HF/6-31G(d), B3LYP/6-31 + G(d,p), B3LYP/6-311 + G(2d,p) and MP2/6-311+G(2d,p). The data set of 77 aromatic aldehydes was divided into a training set (n = 58) and a test (n = 19) set, and 58 models were developed using partial least squares (PLS) and genetic partial least squares (G/PLS). We evaluated the overall predictive capacity of the models based on leave-one-out predictions for the training set compounds and model derived predictions for the test set compounds. For both PLS and G/PLS, the models built at the HF/6-31G(d) level show better predictivity (based on overall prediction) than the models developed at any of the other five levels. Further validation was also performed utilizing (process and model) randomization tests. We show that improved predictive QSAR models for aldehydic toxicity to Tetrahymena pyriformis can be generated using QTMS descriptors along with log K(o/w).

Quantitative structure-activity relationship (QSAR) modeling of EC50 of aquatic toxicities for Daphnia magna

The experimental EC(50) toxicities toward Daphnia magna for a series of 130 benzoic acids, benzaldehydes, phenylsulfonyl acetates, cycloalkane-carboxylates, benzanilides, and other esters were studied using the Best multilinear regression algorithm (BMLR) implemented in CODESSA. A modified quantitative structure-activity relationships (QSAR) procedure was applied guaranteeing the stability and reproducibility of the results. Separating the initial data set into training and test subsets generated three independent models with an average R(2) of .735. A five-descriptor general model including all 130 compounds, constructed using the descriptors found effective for the independent subsets, was characterized by the following statistical parameters: R(2) = .712; R(2)(cv) = .676; F = 61.331; s(2) = 0.6. The removal of two extreme outliers improved significantly the statistical parameters: R(2) = .759; R(2)(cv) = .728; F = 77.032; s(2) = 0.499. The sensitivity of the general model to chance correlations was estimated by applying a scrambling procedure involving 20 randomizations of the original property values. The resulting R(2) = .192 demonstrated the high robustness of the model proposed. The descriptors appearing in the obtained models are related to the biochemical nature of the adverse effects. An additional study of the EC(50)/LC(50) relationship for a series of 28 compounds (part of our general data set) revealed that these endpoints correlated with R(2) = .98.

Definition of the structural domain of the baseline non-polar narcosis model for Tetrahymena pyriformis

The aim of this work was to develop a high-quality 1-octanol/water partition coefficient-dependent (log P) baseline quantitative structure-activity relationship (QSAR) for the toxicity (log IGC(50)(-1)) of classic non-polar narcotics to Tetrahymena pyriformis, and subsequently use this model to define the domain of applicability for baseline narcosis. The toxicities to T. pyriformis of 514 possible non-polar narcotics were assessed. A QSAR to predict toxicity was created from a training set of 87 classic non-polar narcotics (the saturated alcohols and ketones): log IGC(50)(-1) = 0.78 log P-2.01 (n = 87, r(2) = 0.96). This model was then used to predict the toxicity of the remaining chemicals. The chemicals from the large dataset which were poorly predicted by the model (i.e. the prediction was > +/-0.5 log units from the experimental value) were used to aid the definition of structural categories of chemicals which are not non-polar narcotics. Doing so has enabled the domain for non-polar narcosis to be defined in terms of structural categories. Defining domains of applicability for QSAR models is important if they are to be considered for making predictions of toxicity for regulatory purposes.