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.     

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

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.     

Correlation relationships between physico-chemical properties and gas chromatographic retention index of polychlorinated-dibenzofurans

Correlation relationships between physico-chemical properties including vapor pressures (P), water solubilities (S), Henry's law constants (H(c)), n-octanol-water partition coefficients (K(ow)), sediment-water partition coefficient (K(pw)) and biotic lipid-water partition coefficient (K(bw), bioconcentration factor) of polychlorinated-dibenzofurans (PCDFs) and their gas chromatographic retention indices (GC-RIs) were established. A model equation between GC-RIs (= RI) and these physico-chemical properties (K) of PCDFs was in a form of log K = aRI2 + bRI + c with correlation coefficients (R2) greater than 0.94, except H(c). These equations were derived from six experimental data (five experimental data for log K(bw)) in each physico-chemical properties of PCDFs reported previously. The values of log P, log S, log H(c), log K(ow), log K(pw) and log K(bw) of PCDFs predicted by these equations based on their GC-RIs in the present study derviated from those calculated by the solubility parameters for fate analysis method in a previous study by 0.49, 0.32, 0.11, 0.34, 0.14 and 0.22 log units, respectively.     

The aquatic toxicity of anionic surfactants to Daphnia magna--a comparative QSAR study of linear alkylbenzene sulphonates and ester sulphonates

This paper develops quantitative structure activity relationships (QSARs) for the acute aquatic toxicity of the anionic surfactants linear alkylbenzene sulphonates (LAS) and ester sulphonates (ES) to Daphnia magna, the aim being to investigate the modes of action by comparing the QSARs for the two types of surfactant. The generated data for ES have been used to develop a QSAR correlating toxicity with calculated log P values: log(1/EC50)= 0.78 log P+1.37. This equation has an intercept 1.1 log units lower than a QSAR for linear alkylbenzene sulphonates (LAS). The findings suggest that either ES surfactants act by a different mode of action to LAS and other anionic surfactants or the log P calculation method introduces a systematic overestimate when applied to ES.     

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

Vapour pressures, aqueous solubilities, Henry's Law constants, and octanol/water partition coefficients of a series of mixed halogenated dimethyl bipyrroles

Basic physical-chemical properties of five bromine and chlorine containing mixed halogenated dimethyl bipyrroles (HDBPs) were determined using established methods. Subcooled liquid vapour pressures (P(o)(L,25)), aqueous solubilities (S(w,25)), and octanol/water partition coefficients (K(ow)) were determined using the gas chromatography-retention time, generator column, and slow-stirring methods, respectively. Henry's Law constants (H25) were estimated using experimentally-derived P(o)(L) and S(w,25) data. Values of all four properties were generally similar to those reported for other polyhalogenated aromatic compounds [P(o)(L,25) = (7.55-191) x 10(-6) Pa; S(w,25) = (1.0-1.9) x 10(-5) g/l; log K(ow) = 6.4-6.7; H25 = 0.0020-0.14 Pa m3/mol]. The effect of replacing a chlorine with a bromine atom significantly decreased P(o)(L,25) (log P(o)(L,25) = -0.4197 (# bromine atoms) - 2.643, p<0.01) and H25 (log H25 = -0.508 (# bromine atoms) + 0.394, p<0.02). There were no significant effects of bromine/chlorine substitution on S(w,25) or K(ow). A simple Level I equilibrium partitioning model predicted the environmental behaviour of HDBPs to be similar to a tetrabrominated diphenyl ether. Only slight differences in behaviour amongst HDBP congeners were predicted since substitution of a bromine for a chlorine (Cl/Br substitution) atom had less effect than H/Cl or H/Br substitution on P(o)(L,25), S(w,25), H25, and K(ow).     

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.     

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.