Molecular orbital studies on brominated diphenyl ethers. Part II--reactivity and quantitative structure-activity (property) relationships

Polybrominated diphenyl ethers (PBDEs) are widely used as flame retardants and are increasingly turning up in the environment. Their structural similarities to polychlorinated biphenyls and thyroid hormones suggest they may be a risk to human health. The present study examines the reactivity of brominated diphenyl ethers (BDEs) on the basis of the electronic structures as calculated by semiempirical AM1 self-consistent field molecular orbital (SCF-MO) method. Frontier orbital energies were used to elucidate the reactivity of BDEs in electrophilic, nucleophilic and photolytic reactions. From an examination of the frontier electron densities, the regioselectivity, or orientation, of metabolic reactions of BDEs was predicted. Furthermore, satisfactory quantitative structure-activity (property) relationship (QSAR and QSPR) models were derived to calculate gas chromatographic and ultraviolet spectral properties and luciferase induction activities from the AM1-computed electronic parameters.     

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.     

Quantitative structure-activity relationships of nitroaromatics toxicity to the algae (Scenedesmus obliguus)

The DFT-B3LYP method, with the basis set 6-311G( * *), was employed to calculate the molecular geometries and electronic structures of 25 nitroaromatics. The acute toxicity (-lgEC(50)) of these compounds to the algae (Scenedesmus obliguus) along with hydrophobicity described by logK(OW), and two quantum chemical parameters-energy of the lowest unoccupied molecular orbital, E(LUMO), and the charge of the nitro group, [ForQ(NO2), were used to establish the quantitative structure-activity relationships (QSARs). For 18 mononitro derivatives, the hydrophobicity parameter logK(OW) could interpret the toxic mechanism successfully. Dinitro aromatic compounds were susceptible to be reduced to aniline for their electrophilic nature. Their toxicity was controlled mainly by electronic factors instead of hydrophobicity. The electronic parameters, E(LUMO) and Q(NO2), were used to yield the following model: -lg EC(50) = 3.746 - 25.053 E(LUMO) + 6.481 Q(NO2) (n=22, R=0.926, SE=0.206, F=56.854, P<0.001). The predicted toxic values using the above equation are in good agreement with the experimental values.     

Quantitative structure-activity relationships for toxicity and genotoxicity of halogenated aliphatic compounds: wing spot test of Drosophila melanogaster

Halogenated aliphatic compounds were evaluated for toxic and genotoxic effects in the somatic mutation and recombination test employing Drosophila melanogaster. The tested chemicals included chlorinated, brominated and iodinated; mono-, di- and tri-substituted; saturated and unsaturated alkanes: 1,2-dibromoethane, 1-bromo-2-chloroethane, 1-iodopropane, 2,3-dichloropropene, 3-bromo-1-propene, epibromohydrin, 2-iodobutane, 3-chloro-2-methylpropene, 1,2,3-trichloropropane, 1,2-dichloroethane, 1,2-dichlorobutane, 1-chloro-2-methylpropane, 1,3-dichloropropane, 1,2-dichloropropane, 2-chloroethymethylether, 1-bromo-2-methylpropane and 1-chloropentane. N-methyl-N-nitrosourea served as the positive and distilled water as the negative control. The set of chemicals for the toxicological testing was selected by the use of statistical experiment design. Group of unsaturated aliphatic hydrocarbons were generally more toxic than saturated analogues. The genotoxic effect was observed with 14 compounds in the wing spot test, while 3 substances did not show any genotoxicity by using the wing spot test at 50% lethal concentration. The highest number of wing spots was observed in genotoxicity assay with 1-bromo-2-chloroethane, 1,2-dichloroethane, 1,2-dibromoethane and 1-iodopropane. Nucleophilic superdelocalizability calculated by quantum mechanics appears to be a good parameter for prediction of both toxicity and genotoxicity effects of halogenated aliphatic compounds.     

Polychlorinated dibenzo-p-dioxins and dibenzofurans: Correlation between in vitro and in vivo structure-activity relationships (SARs)

Polychlorinated dibenzofurans (PCDFs) and dibenzo-p-dioxins (PCDDs) elicit a number of common biologic and toxic responses which are triggered by their initial binding to a cytosolic receptor protein. These effects include the induction of several cytochrome P-448 dependent monooxygenases (eg, aryl hydrocarbon hydroxylase, AHH), body weight loss and thymic atrophy. The dose-response effects of selected PCDFs on AHH induction in rat hepatoma H-4-II E cells and cytosolic receptor binding affinities have been determined. The results of these and studies demonstrate the remarkable effects of structure on the activity of PCDFs. A systematic study of each of the four different position for chlorine substitution in the dibenzofuran ring system showed that the toxic and biologic potencies of these compounds varied with respect to differential chlorine substitution at all four position, i.e. C-3(7) > C-2(8) >C-4(6) > C-1(9). SARs for PCDDs confirmed the importance of the lateral CI substituents and also showed that 1,2(or 6,7-) substituted PCDDs were more active than the corresponding 1,3-dichloro analogs. In addition, there were significant decreases in activity with increasing non-lateral CI substitution. The SARs for PCDFs were different from the PCDDs and this was directly related to the asymmetric structure of the former group of compounds.     

Studies of 3D-quantitative structure-activity relationships on a set of nitroaromatic compounds: CoMFA,advanced CoMFA and CoMSIA

By using comparative molecular field analysis (CoMFA), advanced CoMFA and comparative molecular similarity index analysis (CoMSIA) methods, the 3D relationships between the structures of 35 nitroaromatic compounds and their toxicities have been investigated to yield statistically reliable models of considerable predictive power. In contrast to CoMFA, CoMSIA produces better results for the correlation. Moreover, the obtained CoMSIA contour maps that interpret the correlations in terms of field contributions allow physicochemical properties relevant for binding to be easily mapped back onto molecular structures, and thus elucidate structural features among ligands that are responsible for toxicities. Besides, most of the highlighted regions in CoMSIA and CoMFA contour maps are mirrored by features in the surrounding environment. Thereby, CoMFA and CoMSIA both help to give explanations of the toxic mechanism of tested compounds.     

Quantitative structure-activity relationships for polycyclic aromatic hydrocarbons: Correlation between molecular connectivity,physico-chemical properties,bioconcentration and toxicity in Daphnia pulex

High correlation coefficients (r = 0.9444 – 0.9996) were calculated for the linear one-parameter correlation between molecular weight, molecular connectivity indices, HPLC and TLC retention indices, partition coefficients between n-octanol or sediment and water, and bioconcentration factors or lethal concentrations in Daphnia Pulex of unsubstituted polycyclic aromatic hydrocarbons.     

Molecular hologram derived quantitative structure-property relationships to predict physico-chemical properties of polychlorinated biphenyls

Polychlorinated biphenyls (PCBs) congeners with various degrees of chlorination and substitution patterns are among the most widespread and persistent man-made organic pollutants. They are toxic, lipophilic and tend to be bioaccumulated. The knowledge of the physico-chemical properties is very useful to explain the environmental behavior of PCBs and to perform an exposure assessment. In this paper, we have used a new molecular representation, the molecular hologram, to generate quantitative structure-property relationship models to predict the physico-chemical properties of biphenyl and all of its chlorinated congeners. The investigated properties include 1-octanol/water partition coefficient (logK(ow)), aqueous solubility (-logS(w)), aqueous activity coefficient (-logY(w)), Total molecular surface area, Henry's law constant (logH). The results show that this new quantitative structure-activity relationship approach presents highly predictive models for important physico-chemical properties of PCBs.     

Reverse quantitative structure-activity relationship for modelling the sorption of esfenvalerate to dissolved organic matter. A multivariate approach

The sorption of the pyrethroid, esfenvalerate, to the dissolved and/or dispersed fraction of eight different natural humic compounds has been investigated. The dissolved organic matters (DOMs) included in this study originate from ground water, soil pore water, and surface waters. Sorption was modelled at DOM concentration levels where equilibrium partitioning of esfenvalerate between DOM and the aqueous bulk phase prevails. The inherent characteristics of the eight different humic materials, quantified in the preceding paper by Thomsen et al. (2002, this issue (PII: S0045-6535(02)00335-1)), have been used as explanatory variables for modelling this equilibrium partitioning. Using a reverse QSAR approach based on by projection-into-latent-structure regression (PLS-R) inherent sorbent properties determining for the sorption affinity of esfenvalerate to DOM were analysed. For all humic substances a decrease in the DOM-normalised equilibrium-partitioning coefficient, K_DOM, with increasing concentration of DOM was observed. Significant variations in K_DOM values, as function of the inherent characteristics of the individual humic substances, were found at DOM concentrations of 75 and 100 ppm, respectively. The latter is a strong indication of variations in sorption mechanisms of esfenvalerate to DOM of varying inherent properties. Groupings in the principal property space quantifying DOMs may indicate that separate models are needed for quantifying the equilibrium partitioning to different classes of DOM.     

Three-dimensional quantitative structure activity relationships of quorum-sensing and biofilm inhibitors in gram-negative bacteria

Thirty N-acyl homoserine lactone (AHL) analogs with variable antibacterial activity and displaying inhibition of biofilm formation were selected to develop models for establishing three-dimensional quantitative structure-activity relationships (3D-QSAR). Comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA) were carried out to determine the optimum structural requirements for selectivity and potency of quorum-sensing and bacterial biofilm inhibition. The best CoMFA model predicted a q² value of 0.519 and an r² value of 0.984 and revealed that electrostatic and steric properties play a significant role in potency and selectivity. The CoMSIA model predicted a q² value of 0.411 and an r² value of 0.938 based on a combination of steric, electrostatic, and hydrophobic effects. The analysis of the contour maps from each model provide insight into the structural requirements for increasing the activity of a compound. Consequently, manipulating the chemical and physical properties of substituted acyl groups on the homoserine lactone moiety can provide important information toward enhancing the antibacterial properties of the target chemical compound.     

Three-dimensional quantitative structure activity relationships of quorum-sensing and biofilm inhibitors in gram-negative bacteria

Thirty N-acyl homoserine lactone (AHL) analogs with variable antibacterial activity and displaying inhibition of biofilm formation were selected to develop models for establishing three-dimensional quantitative structure-activity relationships (3D-QSAR). Comparative molecular field analysis (CoMFA) and comparative similarity indices analysis (CoMSIA) were carried out to determine the optimum structural requirements for selectivity and potency of quorum-sensing and bacterial biofilm inhibition. The best CoMFA model predicted a q2 value of 0.519 and an r2 value of 0.984 and revealed that electrostatic and steric properties play a significant role in potency and selectivity. The CoMSIA model predicted a q2 value of 0.411 and an r2 value of 0.938 based on a combination of steric, electrostatic, and hydrophobic effects. The analysis of the contour maps from each model provide insight into the structural requirements for increasing the activity of a compound. Consequently, manipulating the chemical and physical properties of substituted acyl groups on the homoserine lactone moiety can provide important information toward enhancing the antibacterial properties of the target chemical 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)).     

多元逐步回归对苯胺类化合物结构与毒性模型研究

采用Chemoffice6．0中MOPAC-AMl量子化学法计算了24种苯胺类化合物的6种量子化学结构参数．其中取17个化合物作为样本集对-lgEC50进行多元逐步回归分析．得到最佳方程．经自由度校正的回归系数R=0．985。应用所建立的QSAR模型验证了苯胺类化合物的EC50值．并通过“Jackknife”中的逐一抽取法进行模型检验,得出该模型具有很好的稳定性．平均残差仅为0．05个对数单位．小于文献值。经过7个预测样本对该模型进行验证．结果表明．该模型具有很好的预测能力。同时分析了苯胺类化合物的毒性机理。     

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.