Mechanism of concentration addition toxicity: they are different for nonpolar narcotic chemicals, polar narcotic chemicals and reactive chemicals

According to the toxicity mechanism of the individual chemicals, the concentration addition toxicity mechanism is revealed for nonpolar-narcotic-chemical mixtures, polar-narcotic-chemical mixtures and reactive-chemical mixtures, respectively. For nonpolar-narcotic-chemical mixtures, the partitioning of individual chemicals from water to biophase was determined, and the result shows that their concentration additive effect results from no competitive partitioning among individual chemicals. For polar-narcotic-chemical mixtures, their toxicity are contributed by two factors (the total baseline toxicity and the hydrogen bond donor activity of individual chemicals), and it is the concentration additive effect for either of these two factors that leads to their concentration addition toxicity. In addition, the interactions between the reactive chemicals and the biological macromolecules are discussed thoroughly. The results suggest that the net effect of these interactions is zero, and it is this zero net effect that leads to the concentration addition toxicity mechanism for reactive-chemical mixtures.     

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.     

3D QSAR studies of dioxins and dioxin-like compounds using CoMFA and CoMSIA

In the present study we have performed comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) on structurally diverse ligands of Ah (dioxin) receptor to explore the physico-chemical requirements for binding. All CoMFA and CoMSIA models have given q(2) value of more than 0.5 and r(2) value of more than 0.84. The predictive ability of the models was validated by an external test set, which gave satisfactory predictive r(2) values. Best predictions were obtained with CoMFA model of combined modified training set (q(2) = 0.631, r(2) = 0.900), giving predictive residual value = 0.02 log unit for the test compound. Addition of CoMSIA study has elucidated the role of hydrophobicity and hydrogen bonding along with the effect of steric and electrostatic properties revealed by CoMFA. We have suggested a model comprises of four structurally different compounds, which offers a good predictability for various ligands. Our QSAR model is consistent with all previously established QSAR models with less structurally diverse ligands.     

Quantitative structure-activity relationship (QSAR) models for polycyclic aromatic hydrocarbons (PAHs) dissipation in rhizosphere based on molecular structure and effect size

Rhizoremediation is a significant form of bioremediation for polycyclic aromatic hydrocarbons (PAHs). This study examined the role of molecular structure in determining the rhizosphere effect on PAHs dissipation. Effect size in meta-analysis was employed as activity dataset for building quantitative structure-activity relationship (QSAR) models and accumulative effect sizes of 16 PAHs were used for validation of these models. Based on the genetic algorithm combined with partial least square regression, models for comprehensive dataset, Poaceae dataset, and Fabaceae dataset were built. The results showed that information indices, calculated as information content of molecules based on the calculation of equivalence classes from the molecular graph, were the most important molecular structural indices for QSAR models of rhizosphere effect on PAHs dissipation. The QSAR model, based on the molecular structure indices and effect size, has potential to be used in studying and predicting the rhizosphere effect of PAHs dissipation.     

The effect of nitroaromatics' composition on their toxicity in vivo: novel, efficient non-additive 1D QSAR analysis

Novel 1D QSAR approach that allows analysis of non-additive effects of molecular fragments on toxicity has been proposed. Twenty-eight nitroaromatic compounds including some well-known explosives have been chosen for this study. The 50% lethal dose concentration for rats (LD50) was used as the estimation of toxicity in vivo to develop 1D QSAR models on the framework of Simplex representation of molecular structure. The results of 1D QSAR analysis show that even the information about the composition of molecules provides the main trends of toxicity changes. The necessity of consideration of substituents' mutual impacts for the development of adequate QSAR models of nitroaromatics' toxicity was demonstrated. Statistic characteristics for all the developed partial least squares QSAR models, except the additive ones are quite satisfactory (R2=0.81-0.92; Q2=0.64-0.83; R2 test=0.84-0.87). A successful performance of such models is due to their non-additivity i.e. possibility of taking into account the mutual influence of substituents in benzene ring which plays the governing role for toxicity change and could be mediated through the different C-H fragments of the ring. The correspondence between observed and predicted by these models toxicity values is good. This allowing combine advantages of such approaches and develop adequate consensus model that can be used as a toxicity virtual screening tool.     

The importance of outlier detection and training set selection for reliable environmental QSAR predictions

Empirical QSAR models are only valid in the domain they were trained and validated. Application of the model to substances outside the domain of the model can lead to grossly erroneous predictions. Partial least squares (PLS) regression provides tools for prediction diagnostics that can be used to decide whether or not a substance is within the model domain, i.e. if the model prediction can be trusted. QSAR models for four different environmental end-points are used to demonstrate the importance of appropriate training set selection and how the reliability of QSAR predictions can be increased by outlier diagnostics. All models showed consistent results; test set prediction errors were very similar in magnitude to training set estimation errors when prediction outlier diagnostics were used to detect and remove outliers in the prediction data. Test set prediction errors for substances classified as outliers were much larger. The difference in the number of outliers between models with a randomly and systematically selected training illustrates well the need of representative training data.     

Holographic QSAR of selected esters

The HQSAR (Holographic QSAR) method, which has been recently developed, can offer the ability to rapidly and easily generate QSAR models of high statistical quality and predictive value. HQSAR analysis requires selecting values for parameters that specify the size of the hologram that is to be used, and the size and type of fragment substructures that are to be encoded. The color coding is provided by HQSAR to reflect which molecular fragments may be important contributors to the biological activity. In this work, we studied the quantitative structure activity relationship of selected esters using the HQSAR method. A robust HQSAR model with r² (non-cross-validated regression coefficient) of 0.981 and q² (cross-validated regression coefficient) of 0.912, was developed after optimizing the fragment size and the hologram length. The color coding analysis, which has rarely been reported before, was done here to explain the outlier successfully.     

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.     

A comparison of DEMETRA individual QSARs with an index for evaluation of uncertainty

Quantitative structure-activity relationships (QSARs) urgently need to be applied in regulatory programs. Many QSAR models can predict the effect of a wide range of substances to different endpoints, particularly in the case of ecotoxicity, but it is difficult to choose the most appropriate model on the basis of the requirements of the application. During the EC-funded project DEMETRA (www.demetra-tox.net) a huge number of QSAR models have been developed for the prediction of different ecotoxicological endpoints. DEMETRA individual models on rainbow trout LC50 after 96 h, water flea LC50 after 48 h and honey bee LD50 after 48 h have been used as a QSAR database to test the advantages of a new index for evaluating model uncertainty. This index takes into consideration the number of outliers (weighted on the total number of compounds) and their root mean square error. Application on the DEMETRA QSAR database indicated that the index can identify the models with the best performance with regard to outliers, and can be used, together with other classical statistical measures (e.g., the squared correlation coefficient), to support the evaluation of QSAR models.     

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.     

A new hybrid system of QSAR models for predicting bioconcentration factors (BCF)

The aim was to develop a reliable and practical quantitative structure-activity relationship (QSAR) model validated by strict conditions for predicting bioconcentration factors (BCF). We built up several QSAR models starting from a large data set of 473 heterogeneous chemicals, based on multiple linear regression (MLR), radial basis function neural network (RBFNN) and support vector machine (SVM) methods. To improve the results, we also applied a hybrid model, which gave better prediction than single models. All models were statistically analysed using strict criteria, including an external test set. The outliers were also examined to understand better in which cases large errors were to be expected and to improve the predictive models. The models offer more robust tools for regulatory purposes, on the basis of the statistical results and the quality check on the input data.     

Toxicology of benzyl alcohols: a QSAR analysis

There is an evidence that benzyl alcohols may exhibit toxicity via a radical mechanism. To test this possibility, we studied the toxicity of para substituted benzyl alcohols on rapidly dividing cancer cells (L1210 leukemia). This system has previously found utility in studying the apparent radical toxicity of a variety of phenols. However, no evidence could be found for an electronic effect and the cellular toxicity was associated primarily with hydrophobicity. Comparison of this quantitative structure-activity relationships (QSAR) with others for the reactions of benzyl alcohols in diverse systems provides insight into mechanisms of action. A QSAR for the interaction of benzyl alcohols with protozoa yields an equation that is dependent on both hydrophobicity and acidity of the OH group versus a mixture of bacteria and fungi, the critical dependence on hydrophobicity prevails with a small dependence on a resonance-stabilized, radical mediated electronic effect. The chloramphenicols provide an instructive example, where the radical mediated electronic effect overshadows the hydrophobic contribution to bacterial toxicity. These various QSAR for benzyl alcohols indicate that mechanisms of growth inhibition in vitro vary depending on cell/organism type, the strength of the bond and lability of the hydrogen, and the strength of the initiating radical reagent.     

Computer-based QSARs for predicting mixture toxicity of benzene and its derivatives

During the past decades, the Quantitative structure-activity relationships (QSARs) have been proven to be reliable tools when little or no empirical data are available in medicinal chemistry, biochemistry, toxicology, and environmental sciences. However, only few studies that quantitatively predict mixtures toxicity have been reported. In this study, the QASR models for the binary mixtures toxicity of 12 benzene and its derivatives, including eight non-polar-narcotic compounds and four polar narcotic compounds were developed, without reference to exact toxicity mechanisms of single compounds. All parameters for the QSAR studies were defined on the basis of quantum mechanical calculations and these parameters were selected by the stepwise procedure. The results of this study provided a simple means of predicting the binary mixtures toxicity from the chemical structure.