Mechanism-based quantitative structure-activity relationships for the inhibition of substituted phenols on germination rate of Cucumis sativus

Comparative inhibition activity (GC50) of 42 structurally diverse substituted phenols on seed germination rate of Cucumis sativus was investigated. Quantitative structure-activity relationships (QSARs) were developed by using hydrophobicity (1-octanol/water partition coefficient, logKow) and electrophilicity (the energy of the lowest unoccupied molecule orbital, Eluma) for the toxicity of phenols according to their modes of toxic action. Most phenols elicited their response via a polar narcotic mechanism and a highly significant log Kow-based model was obtained (GC50 = 0.92 log Kow + 1.99, r2 0.84, n = 29). The inclusion of E(lumo) greatly improved the predictive power of the polar narcotic QSAR (GC50 = 0.88 log Kow - 0.30E(lumo) + 1.99, r2 = 0.93, n = 29). pKa proved to be an insignificant influencing factor in this study. Poor correlation with hydrophobicity and strong correlation with electrophilicity were observed for the nine bio-reactive chemicals. Their elevated toxicity was considerably underestimated by the polar narcotic logKow-dependent QSAR. The nine chemicals consist of selected nitro-substituted phenols, hydroquinone, catechol and 2-aminophenol. Their excess toxic potency could be explained by their molecular structure involving in vivo reaction with bio-macromolecules. Strong dissociation of carboxyl group of the four benzoic acid derivatives greatly decreased their observed toxicity. In an effort to model all chemicals including polar narcotics and bio-reactive chemicals, a response-surface analysis with the toxicity, logKow and E(lumo) was performed. This resulted in a highly predictive two-parameter QSAR for most of the chemicals (GC50 = 0. 70 logKow - 0.66E(lumo) + 2.17, r2 = 0.89, n = 36). Catechol and 2,4-dinitrophenol proved to be outliers of this model and their much high toxicity was explained.     

Regression comparisons of Tetrahymena pyriformis and Poecilia reticulata toxicity

The toxicity data of chemicals common to both the Poecilia reticulata mortality assay and the Tetrahymena pyriformis growth impairment assay were evaluated. Two chemicals were not toxic at saturation in the T. pyriformis assay. In addition, due to abiotic transformation, a third chemical was removed from further consideration. Each chemical was a priori assigned a mode of toxic action: neutral non-covalent, polar non-covalent, or electrophilic covalent toxicity. To further investigate comparisons between endpoints, polar and electrophilic chemicals were separated into class-based groups. The polar non-covalent chemicals were separated into phenols and anilines, while the electrophilic chemicals were separated into those reacting via Schiff-base formation (i.e., aldehydes) and those reacting via bimolecular substitution to a nucleophile (i.e., selected nitroaromatics). A comparison of toxic potency as a collective set was statistically described by the relationship; log(LC50(-1)) = 1.05(log(IGC50(-1))) + 0.56, n = 124; r2 = 0.85; s = 0.42; F = 682; Pr > F = 0.0001. The relationship between endpoints was inversely proportional to reactivity associated with the mode of action. While the comparative toxicity for neutral narcotics exhibited an excellent fit (r2 = 0.94), the fits for polar narcotics and electrophiles were poorer, r2 = 0.69 and 0.62, respectively. Investigations into class-based groupings indicated fit of toxic potency data for aldehydes (r2 = 0.85) and phenols (r2 = 0.81) were quite good. However, fits for anilines (r2 = 0.43) and nitroaromatics (r2 = 0.68) revealed that toxicity was not as well related between endpoints for these chemicals.     

Acute toxicity of benzene derivatives to the tadpoles (Rana japonica) and QSAR analyses

Acute lethal toxicity (the negative logarithm of molar concentrations of 12 h acute median lethal, expressed as 12 h-log1/LC50) of 46 benzene derivatives to Rana japonica tadpoles was determined. 1-octanol/water partition coefficient (logKow)-dependent models were developed to study the toxicity of different categories chemicals. In an effort to model all chemicals, response surface analyses and stepwise multiple regression analyses were performed and successful models were obtained. A general and robust QSAR model was achieved with the combined application of variables reflecting hydrophobicity, electric property, and molecular size respectively (12h-log1/LC50 = 0.393logKow - 0.428Elumo + 0.0110Vol. + 1.362 n = 51, r2 = 0.834) using stepwise multiple regression analyses. Because of strong dissociation of carboxyl group greatly decreasing their observed toxicity, using logDow in instead of logKow the quality of the models is greatly improved. The conventional r2 and cross-validation r2(CV) were 0.914 and 0.785, respectively, indicating that QSAR was both internally consistent and highly predictive.     

Structure-activity relationships and response-surface analysis of nitroaromatics toxicity to the yeast (Saccharomyces cerevisiae)

Inhibition of growth of the yeast Saccharomyces cerevisiae (Cmiz, the minimum concentration that produced a clear inhibition zone within 12 h) for 24 nitroaromatic compounds was investigated and a quantitative structure-activity relationship (QSAR) developed based on hydrophobicity expressed as the l-octanol/water partition coefficient in logarithm form, log K(ow), electrophilicity based on the energy of the lowest unoccupied orbital (E(lumo)). All nitrobenzene derivatives exhibited enhanced reactive toxicity than baseline. The toxicities of mono-nitrobenzenes and di-nitrobenzenes were elicited by different mechanisms of toxic action. For mono-nitro-derivatives, both significant log K(ow) based and strong E(lumo)-dependent relationships were observed indicating that their toxicities were affected both by the penetration process and the interaction with target sites of interaction. The toxicities of di-nitrobenzenes were greater than mono-nitrobenzenes and no log K(ow)-dependent but highly significant E(lumo)-based relationship was obtained. This suggests that toxicity of di-nitrobenzenes was highly electrophilic and involved mainly their in vivo electrophilic interaction with biomacromolecules. In an effort to model the elevated toxicity of all nitrobenzenes, a response-surface analysis was performed and this resulted in a highly predictive two-variable QSAR without reference to their exact mechanisms (Cmiz = 0.41 log K(ow) - 0.89 E(lumo) - 0.46, r2 = 0.87, Q2 = 0.86, n = 24).     

Acute toxicity of substituted phenols to Rana japonica tadpoles and mechanism-based quantitative structure-activity relationship (QSAR) study.

Acute 12 h and 24 h lethal toxicity (12 h-LC50 and 24 h-LC50) of 31 substituted phenols to Rana japonica tadpoles was determined. Results indicate that toxicity of phenols to tadpoles varied only slightly with length of exposure and the 12-h test could serve as surrogate of the 24-h test. A mechanism-based quantitative structure-activity relationship (QSAR) method was employed and 1-octanol/water partition coefficient (log K(ow))-dependent models were developed to study different modes of toxic action. Most phenols elicited their response via a polar narcotic mechanism and an excellent logK(ow)-dependent model was obtained. Soft electrophilicity and pro-electrophilicity were observed for some phenols and a good log K(ow)-dependent model was also achieved. Additionally, the significant dissociation of carboxyl on benzoic acid derivatives sharply reduced their toxicity. A statistically robust QSAR model was developed for all studied compounds with the combined application of log K(ow), energy of lowest unoccupied orbital (E(lumo)), heat of formation (HOF) and the first-order path molecular connectivity dices (1chi(p)).     

QSAR study on the toxicity of substituted benzenes to the algae (Scenedesmus obliquus).

50% effective inhibition concentration 48h-EC50 of 40 substituted benzenes to the algae (Scenedesmus obliquus) was determined. The energy of the lowest unoccupied molecular orbital (E(LUMO)) was calculated by the quantum chemical method MOPAC6.0-AM1. By using E(LUMO) and the hydrophobicity parameter log K(OW) the quantitative structure-activity relationship model (QSAR) was developed: log1/EC50=0.272 logK(OW) - 0.659E(LUMO) + 2.54, R2 = 0.793, S.E. = 0.316, F = 71.07, n = 40. A series of equations were obtained about the measured EC50 values of different subclasses of compounds. For those compounds containing double -NO2, their toxicity may be related chiefly to the intracellular reduction of -NO2 obtaining electron, while for anilines and phenols, K(OW) contributes most to the QSAR and E(LUMO) very little.     

Comparison of in vitro and in vivo acute fish toxicity in relation to toxicant mode of action

Acute toxicity to fish hepatoma cell line PLHC-1 and to juvenile rainbow trout was examined for 18 plant protection products. The main objective was to explore whether hepatoma cells could be used to predict acute toxicity in fish taking into account the mode of toxic action and compound properties. Acute fish toxicity was determined using the OECD guideline test 203 and compared to predicted baseline LC50 of acute fish toxicity calculated with a quantitative structure-activity relationship (QSAR) derived for guppy fish. Cytotoxicity was determined through the inhibition of neutral red uptake (NR(50)) into lysosomes and compared to predicted baseline cytotoxicity derived for goldfish GFS cells. In general, NR50 values were higher by a factor ranging from 3 to 3000 than the corresponding acute LC50. A weak correlation between NR50 and LC50 values was found (log/log: r2=0.62). Also the lipophilicity (log K(ow)) was not a good predictor for cytotoxicity (r2=0.43) and lethality (r2=0.57) of these pesticides. The neutral red assay is detecting general baseline toxicity only. Comparing LC50 data to QSAR results, the compounds can be classified to act as narcotics or reactive compounds with a specific mode of toxic action in fish. The results indicate that limitation of the neutral red assay in predicting acute fish toxicity. A promising alternative might be the assessment of toxicity in a set of in vitro systems addressing also cell-specific functions which are related to the mode of toxic action of the compound.     

Prediction of mixture toxicity with its total hydrophobicity

Based on the C18 Empore disk/water partition coefficient of a mixture, quantitative structure-activity relationships (QSARs) are presented, which are used to predict the toxicity of mixed halogenated benzenes to P. phosphoreum. The predicted toxicity of 10 other related mixtures based on the QSAR model, agree well with the observed data with r2 = 0.973, SE = 0.113 and F = 287.785 at a level of significance P < 0.0001. The joint effect of these chemicals is simple similar action and the toxicity of the mixtures can be predicted from total hydrophobicity and is independent of hydrophobicity of the components or the ratio of the individual chemicals.     

Statistical uncertainty in hazardous terrestrial concentrations estimated with aquatic ecotoxicity data

Since chemicals’ ecotoxic effects depend for most soil species on the dissolved concentration in pore water, the equilibrium partitioning (EP) method is generally used to estimate hazardous concentrations (HC50) in the soil from aquatic toxicity tests. The present study analyzes the statistical uncertainty in terrestrial HC50s derived by the EP-method. For 47 organic chemicals, we compared freshwater HC50s derived from standard aquatic ecotoxicity tests with porewater HC50s derived from terrestrial ecotoxicity tests. Statistical uncertainty in the HC50s due to limited species sample size and in organic carbon–water partitioning coefficients due to predictive error was treated with probability distributions propagated by Monte Carlo simulations. Particularly for specifically acting chemicals, it is very important to base the HC50 on a representative sample of species, composed of both target and non-target species. For most chemical groups, porewater HC50 values were approximately a factor of 3 higher than freshwater HC50 values. The ratio of the porewater HC50/freshwater HC50 was typically 3.0 for narcotic chemicals (2.8 for nonpolar and 3.4 for polar narcotics), 0.8 for reactive chemicals, 2.9 for neurotoxic chemicals (4.3 for AChE agents and 0.1 for the cyclodiene type), and 2.5 for herbicides–fungicides. However, the statistical uncertainty associated with this ratio was large (typically 2.3 orders of magnitude). For 81% of the organic chemicals studied, there was no statistical difference between the hazardous concentration of aquatic and terrestrial species. We conclude that possible systematic deviations between the HC50s of aquatic and terrestrial species appear to be less prominent than the overall statistical uncertainty.     

Estrogenicity and acute toxicity of selected anilines using a recombinant yeast assay

Suspected estrogen modulators include industrial organic chemicals (i.e., xenoestrogens), and have been shown to consist of alkylphenols, bisphenols, biphenylols, and some hydroxy-substituted polycyclic aromatic hydrocarbons. The most prominent structural feature identified to be important for estrogenic activity is a polar group capable of donating hydrogen bonds (i.e., hydroxyl) on an aromatic system. The present study was undertaken to explore the estrogenic activity and acute toxicity of chemicals containing a weaker hydrogen bond donor group on aromatic systems, i.e., the amino substituent. There is a great deal of chemical similarity between aromatic amines (anilines) and aromatic alcohols (phenols). The chemicals chosen for the current study contained an amino-substituted benzene ring with hydrophobic constituents varying in size and shape. Thus, 37 substituted aromatic amines were assayed for estrogenic activity EC50 and acute toxicity LC50 using the Saccharomyces cerevisiae recombinant yeast assay. While the EC50 of 17-beta-estradiol occurs at the 10(-10) range, the aniline with the greatest activity had an EC50 of 10(-6) M. Thus, anilines, in general, are capable only of very weak estrogenic activity in this assay. A comparison of estrogenic potency between the present group of anilines and a set of previously tested analogous phenols indicated that anilines are consistently less estrogenic than phenols. A comparison of hazard indices (EC50/LC50) of these chemicals revealed that, for the vast majority of anilines, the EC50 and LC50 were in the same order of magnitude. More specifically, estrogenic activity of para-substituted alkylanilines increases with alkyl group size up to 5 carbons in length, after which the acute toxicity of the larger alkyl-substituents precluded the ability of the compound to induce the estrogenic response.     

First report on development of quantitative interspecies structure-carcinogenicity relationship models and exploring discriminatory features for rodent carcinogenicity of diverse organic chemicals using OECD guidelines

Different regulatory agencies in food and drug administration and environmental protection worldwide are employing quantitative structure-activity relationship (QSAR) models to fill the data gaps related with properties of chemicals affecting the environment and human health. Carcinogenicity is a toxicity endpoint of major concern in recent times. Interspecies toxicity correlations may provide a tool for estimating sensitivity towards toxic chemical exposure with known levels of uncertainty for a diversity of wildlife species. In this background, we have developed quantitative interspecies structure-carcinogenicity correlation models for rat and mouse [rodent species according to the Organization for Economic Cooperation and Development (OECD) guidelines] based on the carcinogenic potential of 166 organic chemicals with wide diversity of molecular structures, spanning a large number of chemical classes and biological mechanisms. All the developed models have been assessed according to the OECD principles for the validation of QSAR models. Consensus predictions for carcinogenicity of the individual compounds are presented here for any one species when the data for the other species are available. Informative illustrations of the contributing structural fragments of chemicals which are responsible for specific carcinogenicity endpoints are identified by the developed models. The models have also been used to predict mouse carcinogenicities of 247 organic chemicals (for which rat carcinogenicities are present) and rat carcinogenicities of 150 chemicals (for which mouse carcinogenicities are present). Discriminatory features for rat and mouse carcinogenicity values have also been explored.     

Quantitative structure-toxicity relationships of organophosphorous pesticides to fish (Cyprinus carpio)

This study was conducted to determine the relationships between 1381 chemical and structural parameters of 43 organophosphorus pesticides (OPs) and their toxicity to fish, Cyprinus carpio, using ChemOffice 8.03 and Dragon 2.1. By multivariate linear regression and intervariable regression analyses, various equations have been derived to calculate the lethal toxicity value, LC(50), for 43 OPs found in fish with different levels of toxicity. Results show that for all selected OPs, especially those of low toxic OPs (LC(50)< 2.5 mM), the equation, LC(50) = 56.259 - 13.071 lg K(ow)+17.510 MATS8P-17.455 Mor24u - 0.085 MW + 1.706 (lg K(ow))(2) + 2.306 (Mor14e)(2) + 6.849 Mor20 m (n = 43, F = 36.815, r = 0.942, r(adj)(2) = 0.862, SE = 2.899, p < 10(-6)), could account for 86.2% of the variability of the toxic effect. The steric and electronic characteristics and the hydrophobicity of OPs, in particular, are among the most important parameters determining the toxicity of OPs to fish. For the OPs with high toxicity, different structural parameters were introduced into the following two equations: LC(50)=3.795-1.195 (H1p)(2)-0.037 U-2.225 MATS3v-19.593 Tcon (n = 16, F = 56.820, r = 0.977, r(adj)(2) = 0.937, SE = 0.143, p < 10(-6)), where LC(50) is less than 2.5 mM, and LC(50) = 0.341-0.561 (HOMA)(2) + 0.231 HOMA (n = 3,r(adj)(2) = 1), where LC(50) is less than 0.3 mM. These results suggest that chemical and structural parameters could be useful in modeling chemical reactivity within homologous series of OP compounds and elucidating possible mechanisms associated with different levels of toxicity to fish.     

Ecotoxicological Characterisation and Classification of Existing Chemicals

GOAL, SCOPE AND BACKGROUND: In 1998, the International Council of Chemical Associations (ICCA) launched a global initiative to investigate more than 1,000 HPV chemicals (High Production Volume, > or = 1,000 t/a) within the refocused OECD HPV Chemicals Programme. Up to the OECD SIDS Initial Assessment Meeting in April 2004 (SIAM 18) 147 ICCA dossiers (ca. 230 CAS-No) have been assessed based on a harmonised data set. The environmental profile and an ecotoxicological characterisation of these chemicals are presented here. Data for acute aquatic toxicity were correlated among each other, as well as data for fish (LC50, LD50) and rodents (LD50). The data for acute aquatic toxicity are compared with other existing chemicals. METHODS: Data of the ICCA HPV chemicals from the OECD SIAM 11-18 are presented for: log Kow (as an indicator for bioaccumulation potential), biodegradation, acute aquatic toxicity and availability of long-term toxicity data. Correlation analysis was performed with log transformed data and a linear regression model was fitted to the data, if a significant correlation was found. Acute toxicity for fish and acute oral toxicity for rodents were correlated on a molar basis. Acute aquatic toxicity of the chemicals is compared with data from BUA reports 1-234 and a random EINECS sample (Knacker et al. 1995). RESULTS AND DISCUSSION: According to the dossier information, 71 of the 147 ICCA chemicals are not 'readily biodegradable', 21 have a log Kow > or = 3, and 44 are 'toxic' (LC/EC50 < or = 10 mg/L) or 'very toxic' (LC/EC50 < or = 1 mg/L) to aquatic organisms. For 77, only the base set (acute fish, Daphnia and algae) is available, for the rest at least one long-term test (fish or Daphnia) is available and three tests for a mere 14 others. Based on the data presented, the SIAM gives recommendations for Environment and Human Health. 22 chemicals have been identified as a 'candidate for further work' for Environment and 16 for Human Health. The highest correlation coefficient was obtained correlating fish and Daphnia (r2 = 0.79). LC50 (fish) is significantly correlated with LD50 (rodent), but data are widely scattered. The correlation is not improved after transforming LC50 (fish) to LD50 (fish), using BCF QSAR. Based on acute aquatic toxicity, 25.1% of the chemicals from the BUA reports 1-234 are classified as 'very toxic' (LC/EC50 < or = 1 mg/L). This proportion is 2.5-fold higher than the ICCA HPV chemicals and 1.4-fold higher than the random EINECS sample. CONCLUSIONS: Correlation coefficients for aquatic toxicity data are rather uniform (0.57-0.79) compared with literature data, but also the best correlation was observed between fish and Daphnia. Because the scatter around the regression lines is still considerable, simple predictions of ecotoxicity between species are not possible. Correlation of LC50 (fish) and LD50 (rodent) indicates that toxicity is different. Surprisingly, the correlation of fish and rodent toxicity is not improved by transforming LC50 values to internal LD50s. The selection of ICCA chemicals by market significance (production volume) leads to a classification of toxicity, which is more comparable to a random sample of EINECS chemicals than to German BUA chemicals. The latter were chosen for concern (for Environment or Human Health). RECOMMENDATIONS AND OUTLOOK: Of 147 dossiers assessed between SIAM 11-18, ca. 75% were sponsored by the three following countries: Germany (42), USA (37) and Japan (33). The current output is about 50 dossiers per year (70-100 CAS-No), but a trend for an increase of output is noticeable. Industry, national authorities, and OECD work on a further development to speed up the output. The number of chemicals with 'low priority for further work' and the work recommended for the 'candidates' (mainly exposure assessment) indicate that the data presented were adequate for an initial hazard assessment according to OECD requirements. From the ICCA HPV list (n = 880, state of 1999) 44% of the chemicals have data available to cover all SIDS endpoints for Environment and only 33% for Human Health (Allanou et al. 1999). This indicates the importance of the Initiative to provide information on existing chemicals. The authors agree with the expectation "...that the scientific information provided by this global initiative will be considered as an internationally accepted and harmonised basis for further steps of chemicals management." (ICCA 2002 b).     

A novel approach to predict aquatic toxicity from molecular structure

The main aim of the study was to develop quantitative structure-activity relationship (QSAR) models for the prediction of aquatic toxicity using atom-based non-stochastic and stochastic linear indices. The used dataset consist of 392 benzene derivatives, separated into training and test sets, for which toxicity data to the ciliate Tetrahymena pyriformis were available. Using multiple linear regression, two statistically significant QSAR models were obtained with non-stochastic (R2=0.791 and s=0.344) and stochastic (R2=0.799 and s=0.343) linear indices. A leave-one-out (LOO) cross-validation procedure was carried out achieving values of q2=0.781 (scv=0.348) and q2=0.786 (scv=0.350), respectively. In addition, a validation through an external test set was performed, which yields significant values of Rpred2 of 0.762 and 0.797. A brief study of the influence of the statistical outliers in QSAR's model development was also carried out. Finally, our method was compared with other approaches implemented in the Dragon software achieving better results. The non-stochastic and stochastic linear indices appear to provide an interesting alternative to costly and time-consuming experiments for determining toxicity.     

Discrimination of excess toxicity from narcotic effect: Comparison of toxicity of class-based organic chemicals to Daphnia magna and Tetrahymena pyriformis

The discrimination of excess toxicity from narcotic effect plays a crucial role in the study of modes of toxic action for organic compounds. In this paper, the toxicity data of 758 chemicals to Daphnia magna and 993 chemicals to Tetrahymena pyriformis were used to investigate the excess toxicity. The result showed that mode of toxic action of chemicals is species dependent. The toxic ratio (TR) calculated from baseline model over the experimentally determined values showed that some classes (e.g. alkanes, alcohols, ethers, aldehydes, esters and benzenes) shared same modes of toxic action to both D. magna and T. pyriformis. However, some classes may share different modes of toxic action to T. pyriformis and D. magna (e.g. anilines and their derivatives). For the interspecies comparison, same reference threshold need to be used between species toxicity. The excess toxicity indicates that toxicity enhancement is driven by reactive or specific toxicity. However, not all the reactive compounds exhibit excess toxicity. In theory, the TR threshold should not be related with the experimental uncertainty. The experimental uncertainty only brings the difficulty for discriminating the toxic category of chemicals. The real threshold of excess toxicity which is used to identify baseline from reactive chemicals should be based on the critical concentration difference inside body, rather than critical concentration outside body (i.e. EC₅₀ or IGC₅₀). The experimental bioconcentration factors can be greatly different from predicted bioconcentration factors, resulting in different toxic ratios and leading to mis-classification of toxic category and outliers.     

Relative sensitivity of one freshwater and two marine acute toxicity tests as determined by testing 30 offshore E & P chemicals

Acute toxicity of 30 offshore E & P (Exploration and Production) chemicals was measured using the three standard test organisms Daphnia magna (freshwater cladoceran), Acartia tonsa (marine copepod) and Skeletonema costatum (marine diatom alga). Test chemicals included 20 water-soluble and 10 (partially) non-soluble products. For 22 out of the 30 chemicals, the difference in sensitivity between the three tests varied within one order of magnitude. A very good correlation was found between the two marine tests (r = 0.96, P < 0.01, n = 30), and a correlation coefficient of r = 0.78 (P < 0.01, n = 30) was found between D. magna and both A. tonsa and S. costatum, individually. When the comparison of D. magna and A. tonsa sensitivity was based only on the water-soluble chemicals, a significantly higher correlation was obtained (r = 0.84, n = 20), indicating that the sample preparation method used for the (partially) non-soluble chemicals (the water accommodated fraction (WAF) method) induces additional variation between tests performed with different test media. (Partially) non-soluble chemicals are characterised by phase separation or precipitation at the concentrations used for testing. In a WAF-based test, each test concentration/exposure level is prepared separately, and following mixing and separation, only the water phase is used for testing. Toxicity is related to the amount of substance originally added to the mixing vessels. For 25 of the 30 chemicals, D. magna was found to be less sensitive than the marine copepod by a factor >2. The generally higher sensitivity of the marine toxicity tests compared to the Daphnia test emphasise the importance of using marine data for environmental hazard classification as well as for environmental risk assessment purposes.