Quantitative structure-activity relationships for the inhibition toxicity to root elongation of Cucumis sativus of selected phenols and interspecies correlation with Tetrahymena pyriformis

The comparative toxicities of selected phenols to higher plants Cucumis sativus were measured and the negative logarithm molar concentration of the root elongation median inhibition (IRC50) were derived. Quantitative structure-activity relationships (QSARs) were developed to explore the toxicity influencing factors and for predictive purpose. The toxicity data, fell into two classes: polar narcosis and bio-reactive. For polar narcotic phenols, a highly significant two-parameter QSAR based on 1-octanol/water partition coefficient (logKow) and energy of the lowest unoccupied orbital (E(lumo)) was derived (IRC50 = 0.77 log Kow - 0.39E(lumo) + 2.36 n = 22 r2 = 0.89). The five bio-reactive chemicals proved to show elevated toxicity due to their typical substructure involved diverse reactive mechanisms. In an effort to model all chemicals, a robust multiple-variable QSAR combining logKow, E(lumo) and Qmax, the most negative net atomic charge, was developed (IRC50 = 0.65 logKow - 0.72E(lumo) + 0.23Qmax + 2.81 n = 27 r2 = 0.94), indicating that hydrophobicity, electrophilicity and hydrogen bond interaction contribute mainly to the phytotoxicity. The toxicological data was compared with Tetrahymena pyriformis 2-d population growth inhibition toxicity (IGC50) and excellent interspecies correlations were observed both for the polar narcotics and for five reactive chemicals (for polar narcotics: IRC50 = 0.95IGC50 + 1.07 n = 16 r2 = 0.89; for bio-reactive chemicals: IRC50 = 0.98IGC50 + 2.19 n = 5 r2 = 0.97; and for all: IRC50 = 0.93IGC50 + 1.63 n = 21 r2 = 0.87). This suggested that T pyriformis toxicity could serve as a surrogate of C. sativus toxicity for phenols and interspecies correlation also could be established for reactive chemicals.     

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.     

An evaluation of global QSAR models for the prediction of the toxicity of phenols to Tetrahymena pyriformis

This study presents an analysis of the ability of a two-parameter response surface, a multiple linear regression and a neural network model to produce global quantitative structure-activity relationships (QSARs) to predict the toxic potency of phenols to Tetrahymena pyriformis. The phenolic toxicity data set analysed is characterised by multiple mechanisms of toxic action. The study aimed to evaluate the confidence that can be applied to the modelling of the differing mechanisms of action. Assessment of confidence was decided in terms of whether the statistics for the global models reflect the ability of the QSARs to model the individual mechanisms of toxic action present in the data set. The results showed that the global statistics only reflected the ability of models to predict the two non-covalent mechanisms (polar narcosis and respiratory uncoupling), with the metabolically transformed and electrophilic mechanism (pre-electrophiles and soft electrophiles) being modelled poorly by all three model building methods. The results confirm the difficulty in modelling electrophilic mechanisms of toxic action. The results also highlight the fact that this poor predictivity is often 'hidden' in good statistical fit of some global models. In particular these results emphasise that for practical predictive purposes the mechanistic applicability domain is required to give confidence to estimated toxicity values.     

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.     

Toxicity of 58 substituted anilines and phenols to algae Pseudokirchneriella subcapitata and bacteria Vibrio fischeri: comparison with published data and QSARs

A congeneric set of 58 substituted anilines and phenols was tested using the 72-h algal growth inhibition assay with Pseudokirchneriella subcapitata and 15-min Vibrio fischeri luminescence inhibition assay. The set contained molecules substituted with one, two or three groups chosen from -chloro, -methyl or -ethyl. For 48 compounds there was no REACH-compatible algal toxicity data available before. The experimentally obtained EC50 values (mg L⁻¹) for algae ranged from 1.43 (3,4,5-trichloroaniline) to 197 (phenol) and for V. fischeri from 0.37 (2,3,5-trichlorophenol) to 491 (aniline). Only five of the tested 58 chemicals showed inhibitory effect to algae at concentrations >100 mg L⁻¹, i.e. could be classified as “not harmful”, 32 chemicals as “harmful” (10-100 mg L⁻¹) and 21 as “toxic” (1-10 mg L⁻¹). The occupied para-position tended to increase toxicity whereas most of the ortho-substituted congeners were the least toxic. As a rule, the higher the number of substituents the higher the hydrophobicity and toxicity. However, in case of both assays, the compounds of similar hydrophobicity showed up to 30-fold different toxicities. There were also assay/organism dependent tendencies: phenols were more toxic than anilines in the V. fischeri assay but not in the algal test. The comparison of the experimental toxicity data to the data available from the literature as well as to QSAR predictions showed that toxicity of phenols to algae can be modeled based on hydrophobicity, whereas the toxicity of anilines to algae as well as toxicity of both anilines and phenols to V. fischeri depended on other characteristics in addition to logK_ow.     

Embryonic development assay with Daphnia magna: application to toxicity of aniline derivatives.

An assay system using Daphnia magna embryos was applied to investigate the adverse effects of aniline derivatives. The data were compared with our previous data for chlorophenols. This new assay provides useful information to evaluate the toxicity of chemicals and the differences in sensitivity between the life stages. The effects of 15 aniline derivatives on embryonic development of D. magna embryos were determined. At the start of exposure, 2-6-h old eggs (between stages 1 and 2, round in shape, diameter approx. 400 microm), were used. In control and solvent control groups, embryonic development from an egg to a free-swimming animal proceeded completely within 3 days with more than 90% hatchability. Median effective concentrations (EC50s) to reduce the numbers hatched were determined and gross morphological abnormalities of hatched animals recorded. Anilines induced no obvious morphological abnormalities and no developmental delay although premature deaths occurred. However, they affected the number of embryos hatching in a dose-dependent manner. In addition, this embryo assay was more sensitive to aniline derivatives (except for aniline) than acute juveniles immobilization assay. Ratios of 48-h EC50 (juvenile)/3-day EC50 (embryo) for eight anilines were greater than 5.0. Particularly, the ratios of 4-methyl-, 4-ethyl- and 3-methylaniline were 77, 23 and 11, respectively. EC50s for embryos and juveniles were poorly correlated (r = 0.41). This indicated that the sensitivities of the two life stages were different to the effects of anilines. EC50s were poorly correlated (r = -0.097) with the log Kow (1-octanol/water partition coefficient). These results were compared with previous results for phenols.     

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.     

Validation of germination rate and root elongation as indicator to assess phytotoxicity with Cucumis sativus

Germination rate and root elongation, as a rapid phytotoxicity test method, possess several advantages, such as sensitivity, simplicity, low cost and suitability for unstable chemicals or samples. These advantages made them suitable for developing a large-scale phytotoxicity database and especially applicable for developing quantitative structure–activity relationship (QSAR) to study mechanisms of phytotoxicity. In this paper, the comparative inhibition of germination rate and root elongation of Cucumis sativus by selected halogen-substituted phenols and anilines were determined. The suitability of germination rate and root elongation as phytotoxicity endpoints was evaluated. Excellent reproducibility and stability of germination rate and root elongation in the control test, relatively greater sensitivity and similar dose–response relations for all tested compounds were observed. These results together with those of a 2-day test were used to demonstrate the suitability of this phytotoxicity test method. A QSAR was developed for the phytotoxicity mode of action of the tested compounds to C. sativus seeds. Models that combined the logarithm of 1-octanol/water partition coefficient (log K_ow) and the energy of the lowest unoccupied molecular orbital (E_lumo) were developed for both germination rate inhibition and root elongation inhibition. The results of these studies indicate that phytotoxicity of substituted phenols and anilines to C. sativus seeds could be explained by a polar narcosis mechanism. This paper will promote the application of germination rate and root elongation method and the development of large-scale phytotoxicity database, which will provide the fundamental data for QSAR and ecological risk assessment of organic pollutants.     

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.     

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

Comparative use of bacterial, algal and protozoan tests to study toxicity of azo- and anthraquinone dyes

Toxicity of two azo dyes (Reactive Orange 16 (RO16); Congo Red (CR)) and two anthraquinone dyes (Remazol Brilliant Blue R (RBBR); Disperse Blue 3 (DB3)) were compared using bacterium Vibrio fischeri, microalga Selenastrum capricornutum and ciliate Tetrahymena pyriformis. The following respective endpoints were involved: acute toxicity measured as bacterial luminescence inhibition, algal growth inhibition, and the effects on the protozoa including viability, growth inhibition, grazing effect and morphometric effects. In addition, mutagenicity of the dyes was determined using Ames test with bacterium Salmonella typhimurium His(-). DB3 dye was the most toxic of all dyes in the bacterial, algal and protozoan tests. In contrast to other dyes, DB3 exhibited mutagenic effects after metabolic activation in vitro in all S. typhimurium strains used. Of the methods applied, the algal test was the most sensitive to evaluate toxicity of the dyes tested.     

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

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.