Prediction of environmental partition coefficients and the Henry's law constants for 135 congeners of chlorodibenzothiophene

Occurrence, physicochemical properties and environmental behavior of polychlorinated dibenzothiophenes: A comprehensive review and future perspectives

Polychlorinated dibenzothiophenes (PCDTs) are a form of emerging pollutant that has attracted great attention due to their structural resemblance to dioxins, which cast detrimental influence on the ecosystem and human health. This review shows the current status of research on PCDTs, focusing on their environmental occurrence, physicochemical properties, environmental behavior, and toxicity. Studies have suggested that the steps leading to the formation of PCDTs resemble those generating polychlorinated dibenzo-p-dioxin/dibenzofurans (PCDD/Fs), indicating their probable origin from the same sources. Furthermore, they may undergo a dechlorination process as a result of their photodegradation in the environment and metabolic reaction occurring within organisms, which could result in the conversion of these substances into additional pollutants like dibenzothiophene. PCDTs exist widely in the environmental media and have high logKOW values (>4.0), indicating their tendency to bioaccumulate. Moreover, the prediction results of EPI (Estimation Program Interface) Suite demonstrated a strong accumulation capacity for tetra-CDTs in fish compared to other chlorinated PCDTs. The biotransformation half-life of PCDTs would prolong with an increasing number of substituted Cl atoms in fish. A limited number of studies have also suggested that PCDTs can cause damage to the liver and immune system in living organisms, and the toxicity of PCDTs depends on the number and position of substituted Cl atoms. Future studies should be conducted on processes causing PCDT toxicity as well as their behavior and fate in actual environments.

An Insight into the Sorption Behavior of 2,3,7,8-Tetrachlorodibenzothiophene on the Sediments and Paddy Soil from Chaohu Lake Basin

Considering the frequent detection of polychlorinated dibenzothiophenes (PCDTs) in various environmental matrices and the potential ecological health risks, the environmental behavior of such compounds needs to be elucidated further. In this work, the sorption behavior of 2,3,7,8-tetrachlorodibenzothiophene (2,3,7,8-TCDT) onto three sediments and paddy soil from Chaohu Lake were investigated via batch equilibration experiments. From the perspective of sorption kinetics and isotherms, the sorption characteristics and mechanism of 2,3,7,8-TCDT on the above four carriers were compared, and the relationship between their structural characteristics and soil sorption capacity was discussed. Results suggested that rapid sorption played the primary role during the sorption process of 2,3,7,8-TCDT and the corresponding sorption isotherms were well fitted using the Freundlich logarithmic model. Moreover, the effects of pH and dissolved organic matter (DOM) on the sorption of 2,3,7,8-TCDT were investigated. The maximum sorption capacity of 2,3,7,8-TCDT on sediment was under acidic pH condition (pH = 4.0). Meanwhile, DOM at a low level promoted the sorption capacity of sediment toward 2,3,7,8-TCDT, while the high concentration of DOM inhibited this effect. In addition, the values of logK_oc were obtained using high-performance liquid chromatography (HPLC) and did not show any significant correlation with organic carbon (OC) contents, thereby indicating that the partition effect was the dominating influencing factor for the sorption of 2,3,7,8-TCDT both on sediments and soil. This work provides useful data to understand the sorption behavior of 2,3,7,8-TCDT on sediments and soil and assess its potential environmental risk.

Effect of Water Molecule on the Complete Series Reactions of Chlorothiobenzenes with H/·OH: A Theoretical Study

The chlorothiobenzenes (CTBs) are the principal precursors for the formation of polychlorinated thianthrene/dibenzothiophenes (PCTA/DTs), which have high toxicity and wide distribution in the environment. Under the pyrolysis or combustion conditions, CTBs can react with H/·OH radicals to form the chlorothiobenzyl radicals (CTBRs) through abstraction of the chlorothiobenzyl-hydrogen. The water molecule can play an important role in this process. The coupling of CTBRs is the essential first step in forming PCTA/DTs. In this paper, quantum chemical calculations were carried out to investigate the formation of CTBRs from the complete series reactions of 19 chlorothiobenzene (CTB) congeners with H/·OH radicals in the presence of the water molecule. Using the MPWB1K/6-311 + G(3df,2p)//MPWB1K/6-31 + G(d,p) energy level, schematic energy profiles were constructed with the water molecule and then compared with the non-hydrated case. The present study shows that structural parameters and thermal data, as well as CTBRs formation potential from CTBs, are strongly dominated by the chlorine substitution at the ortho-position of CTBs. Meanwhile, the water molecule can promote the CTBR formation from CTBs abstracted by H/·OH, which has a stronger catalysis effect on the H abstraction from CTBs by OH than from CTBs by H. This study may provide reference parameters for future experimental research, which would enhance measures to reduce dioxin emission and establish dioxin control strategies.

Occurrence, distribution and partitioning of polychlorinated dibenzothiophenes (PCDTs) in Chaohu Lake, Southeast China

Polychlorinated dibenzothiophenes (PCDTs) are a class of compounds structurally similar to dioxins that possess various toxicological impacts on living organisms. Unfortunately, information on the levels of PCDTs in freshwater lakes in China is still scarce. In this work, the occurrence of 14 congeners of PCDTs in different matrices (i.e., sediment, suspended particulate matter (SPM), and water) of Chaohu Lake was investigated. It was determined that the concentrations of 14 PCDTs (Σ₁₄PCDTs) in the sediment, SPM, and surface water were 0.40-3.55 ng g^-1 (dry weight, d.w.), 0.38-2.95 ng·g^-1 d.w., and 0.34-2.61 ng L^-1, respectively. The dominant congener found in sediments was 1,2,3,4,7-penta-CDT (19.54%), and 1,3,9-tri-CDT was the predominant congener in SPM (19.13%) and water (20.08%). Medium- and high-chlorinated PCDTs were detected as the major compounds in sediments and SPM. The low-chlorinated PCDTs (e.g., mono-CDTs) have higher relative percentages in the water than those detected in the sediment samples. The annual Σ₁₄PCDT input of the eight main tributaries to Chaohu Lake was 19.90 kg. A strong linear correlation between the Σ₁₄PCDT levels and the organic carbon (OC) content demonstrated that OC had an important influence on the PCDT redistribution in Chaohu Lake. In addition, the organic carbon normalized partitioning coefficient (logK_OC) of PCDTs in the SPM-water system in Chaohu Lake was 1.95-2.49 mL g^-1, and correlations between logK_OC and other typical environment-related properties of PCDTs were established. This study provided useful data on the evaluation of ecological risks of PCDTs in Chaohu Lake.

Mechanism and kinetic properties for the complete series reactions of chloro(thio)phenols with O(3P) under high temperature conditions

Polychlorinated dibenzo-p-dioxins/dibenzofurans (PCDD/Fs) and polychlorinated dibenzothiophenes/thianthrenes (PCDT/TAs) are two groups of dioxin-like compounds with oxygen and sulfur substitution, respectively. Chlorophenols (CPs) and chlorothiophenols (CTPs) are direct precursors in PCDD/F and PCDT/TA formation. The formation of chlorophenoxy radicals (CPRs) and chlorothiophenoxy radicals (CTPRs) from chlorophenols (CPs) and chlorothiophenols (CTPs) with O(3P) is an important initial step for the formation of PCDD/Fs and PCDT/TAs, respectively. In this paper, the formation of CPRs/CTPRs from the complete series reactions of 19 CP/CTP congeners with O(3P) was studied using the density functional theory (DFT) method. The rate constants of each reaction were calculated using canonical variational transition state (CVT) theory along with a small-curvature tunneling (SCT) contribution over a wide temperature range of 600-1200 K. The effect of the chlorine substitution pattern on the structural parameters, thermochemical properties and rate constants in both CPs and CTPs was discussed. This study shows that the reactions between CPs and O(3P) can be affected by the chlorine substitution at the para-position, and the reactions between CTPs and O(3P) are mostly influenced by both ortho-substitutions. The thiophenoxyl-hydrogen abstraction from CTPs by O(3P) is more likely to occur than the phenoxyl-hydrogen abstraction from CPs by O(3P). Comparison of the reactivity of CP/CTPs with O(3P) with our previous work on CP/CTPs with H and OH shows that the order for phenoxyl-hydrogen abstraction potential is CP + OH > CP + O(3P) > CP + H, and the order for thiophenoxyl-hydrogen abstraction potential is CTP + O(3P) > CTP + H > CTP + OH.

The homogeneous gas-phase formation mechanisms of PCPTs/PCDTs/PCDFs from the radical/radical cross-condensation of 2-CPR and 2-CTPR: a theoretical, mechanistic and kinetics study

Polychlorinated phenoxathiins (PCPTs) are one group of dioxin-like compounds, which can be considered to be one-oxygen-substituted polychlorinated thianthrene (PCTA) compounds or one-sulfur-substituted polychlorinated dibenzo-p-dioxin (PCDD) compounds. Owing to their high toxicity and wide distribution, clarifying the formation and emission of PCPTs due to combustion and thermal processes can deepen our understanding of the dioxin formation mechanism and allow reduced-emission and dioxin-control strategies to be established. Chlorophenols (CPs) and chlorothiophenols (CTPs) are direct precursors in PCPT formation. In this paper, the homogeneous gas-phase formation mechanisms of PCPTs, as well as polychlorinated dibenzofurans (PCDFs) and polychlorinated dibenzothiophenes (PCDTs), from the cross-condensation of 2-chlorophenoxy radicals (2-CPRs) and 2-chlorothiophenoxy radicals (2-CTPRs) under thermal and combustion conditions were investigated theoretically using a density functional theory (DFT) method. The reaction priorities and effects of water molecules on the formation mechanisms were discussed. The rate constants of crucial elementary steps were calculated from 600-1200 K. The acute and chronic toxicities of the main products were predicted at three trophic levels. This study shows that routes starting with oxygen-carbon condensation are favored over those starting with sulfur-carbon condensation for PCPT formation, and routes ending with Cl loss can occur more easily than those ending with H loss. Water molecules have a negative catalytic effect on CH-S H-transfer steps but a positive catalytic effect on CH-O H-transfer steps.

Quantum Chemical and Kinetic Study on Radical/Molecule Formation Mechanism of Pre-Intermediates for PCTA/PT/DT/DFs from 2-Chlorothiophenol and 2-Chlorophenol Precursors

Polychlorinated phenoxathiins (PCPTs), polychlorinated dibenzothiophenes (PCDTs), and polychlorinated thianthrenes (PCTAs) are sulfur analogues of polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans (PCDD/DFs). Chlorothiophenols (CTPs) and chlorophenols (CPs) are key precursors for the formation of PCTA/PT/DTs, which can react with H or OH to form chloro(thio)phenoxy radical, sulfydryl/hydroxyl-substituted phenyl radicals, and (thio)phenoxyl diradicals. However, previous radical/radical PCTA/DT formation mechanisms in the literature failed to explain the higher concentration of PCDTs than that of PCTAs under the pyrolysis or combustion conditions. In this work, a detailed thermodynamics and kinetic calculations were carried out to investigate the pre-intermediate formation for PCTA/PT/DTs from radical/molecule coupling of the 2-C(T)P with their key radical species. Our study showed that the radical/molecule coupling mechanism explains the gas-phase formation of PCTA/PT/DTs in both thermodynamic and kinetic perspectives. The S/C coupling modes to form thioether-(thio)enol intermediates are preferable over the O/C coupling modes to form ether-(thio)enol intermediates. Thus, although the radical/molecule coupling of chlorophenoxy radical with 2-C(T)P has no effect on the PCDD/PT formation, the radical/molecule coupling of chlorothiophenoxy radical with 2-C(T)P plays an important role in the PCTA/PT formation. Most importantly, the pre-PCDT intermediates formation pathways from the couplings of sulfydryl/hydroxyl-substituted phenyl radical with 2-C(T)P and (thio)phenoxyl diradicals with 2-C(T)P are more favorable than pre-PCTA/PT intermediates formation pathways from the coupling of chlorothiophenoxy radical with 2-C(T)P, which provides reasonable explanation for the high PCDT-to-PCTA ratio in the environment.

Modeling adsorption of organic pollutants onto single-walled carbon nanotubes with theoretical molecular descriptors using MLR and SVM algorithms

Prediction of adsorption equilibrium coefficients (K) of organic compounds onto single walled carbon nanotubes (SWNTs) from in silico molecular descriptors is of importance for probing potential applications of SWNTs as well as for evaluating environmental behavior and ecological risks of organic pollutants and SWNTs. In this study, two models for predicting logK were developed with multiple linear regression (MLR) and support vector machine (SVM) algorithms. The two models have satisfactory goodness-of-fit, robustness and predictive ability, and the SVM model performs slightly better than the MLR model. The two models are based on the up-to-date experimental dataset consisting of 61 logK values, and the applicability domains cover diverse organic compounds with functional groups > CC<, CC, C₆H₅, >CO, COOH, C(O)O, OH, O, F, Cl, Br, NH₂, NH, >N, >NN<, NO₂, >NC(O)NH₂, >NC(O)NH, S and S(O)(O). The adsorption of organic compounds toward SWNTs is mainly determined by van der Waals forces and hydrophobic interactions. Since only in silico molecular descriptors were employed for the modeling, the developed models are beneficial for prediction purposes.

Metal Oxide Nanomaterial QNAR Models: Available Structural Descriptors and Understanding of Toxicity Mechanisms

Metal oxide nanomaterials are widely used in various areas; however, the divergent published toxicology data makes it difficult to determine whether there is a risk associated with exposure to metal oxide nanomaterials. The application of quantitative structure activity relationship (QSAR) modeling in metal oxide nanomaterials toxicity studies can reduce the need for time-consuming and resource-intensive nanotoxicity tests. The nanostructure and inorganic composition of metal oxide nanomaterials makes this approach different from classical QSAR study; this review lists and classifies some structural descriptors, such as size, cation charge, and band gap energy, in recent metal oxide nanomaterials quantitative nanostructure activity relationship (QNAR) studies and discusses the mechanism of metal oxide nanomaterials toxicity based on these descriptors and traditional nanotoxicity tests.

Mechanistic and Kinetic Studies on the Homogeneous Gas-Phase Formation of PCTA/DTs from 2,4-Dichlorothiophenol and 2,4,6-Trichlorothiophenol

Polychlorinated thianthrene/dibenzothiophenes (PCTA/DTs) are sulfur analogues compounds to polychlorinated dibenzo-p-dioxin/dibenzofurans (PCDD/Fs). Chlorothiophenols (CTPs) are key precursors to form PCTA/DTs. 2,4-DCTP has the minimum number of Cl atoms to form 2,4,6,8-tetrachlorinated dibenzothiophenes (2,4,6,8-TeCDT), which is the most important and widely detected of the PCDTs. In this paper, quantum chemical calculations were carried out to investigate the homogeneous gas-phase formation of PCTA/DTs from 2,4-DCTP and 2,4,6-TCTP precursors at the MPWB1K/6-311+G(3df,2p)//MPWB1K/6-31+G(d,p) level. Several energetically feasible pathways were revealed to compare the formation potential of PCTA/DT products. The rate constants of the crucial elementary reactions were evaluated by the canonical variational transition-state (CVT) theory with the small curvature tunneling (SCT) correction over a wide temperature range of 600–1200 K. This study shows that pathways that ended with elimination of Cl step were dominant over pathways ended with elimination of the H step. The water molecule has a negative catalytic effect on the H-shift step and hinders the formation of PCDTs from 2,4-DCTP. This study, together with works already published from our group, clearly illustrates an increased propensity for the dioxin formation from CTPs over the analogous CPs.

Formation of Chlorotriophenoxy Radicals from Complete Series Reactions of Chlorotriophenols with H and OH Radicals

The chlorothiophenoxy radicals (CTPRs) are key intermediate species in the formation of polychlorinated dibenzothiophenes/thianthrenes (PCDT/TAs). In this work, the formation of CTPRs from the complete series reactions of 19 chlorothiophenol (CTP) congeners with H and OH radicals were investigated theoretically by using the density functional theory (DFT) method. The profiles of the potential energy surface were constructed at the MPWB1K/6-311+G(3df,2p)//MPWB1K/6-31+G(d,p) level. The rate constants were evaluated by the canonical variational transition-state (CVT) theory with the small curvature tunneling (SCT) contribution at 600–1200 K. The present study indicates that the structural parameters, thermal data, and rate constants as well as the formation potential of CTPRs from CTPs are strongly dominated by the chlorine substitution at the ortho-position of CTPs. Comparison with the study of formation of chlorophenoxy radicals (CPRs) from chlorophenols (CPs) clearly shows that the thiophenoxyl-hydrogen abstraction from CTPs by H is more efficient than the phenoxyl-hydrogen abstraction from CPs by H, whereas the thiophenoxyl-hydrogen abstraction from CTPs by OH is less impactful than the phenoxyl-hydrogen abstraction from CPs by OH. Reactions of CTPs with H can occur more readily than that of CTPs with OH, which is opposite to the reactivity comparison of CPs with H and OH.

QSPR modeling of n-octanol/air partition coefficients and liquid vapor pressures of polychlorinated dibenzo-p-dioxins

The molecular geometries of 75 polychlorinated dibenzo-p-dioxins (PCDDs) were optimized using B3LYP/6-31G(*) method. The calculated structural parameters were taken as theoretical descriptors to establish two new novel QSPR models for n-octanol/air partition coefficients (log K(OA)) and subcooled liquid vapor pressure (log P(L)) of PCDDs. The R(2) values of the two models are 0.983 and 0.985 respectively. Their standard deviations of prediction in modeling (SD) are 0.174 and 0.230 respectively. The results of leave-one-out (LOO) cross-validation for training set show that the two models exhibited optimum stability and good predictive power.

Advancing risk assessment of engineered nanomaterials: application of computational approaches

Nanotechnology that develops novel materials at size of 100nm or less has become one of the most promising areas of human endeavor. Because of their intrinsic properties, nanoparticles are commonly employed in electronics, photovoltaic, catalysis, environmental and space engineering, cosmetic industry and - finally - in medicine and pharmacy. In that sense, nanotechnology creates great opportunities for the progress of modern medicine. However, recent studies have shown evident toxicity of some nanoparticles to living organisms (toxicity), and their potentially negative impact on environmental ecosystems (ecotoxicity). Lack of available data and low adequacy of experimental protocols prevent comprehensive risk assessment. The purpose of this review is to present the current state of knowledge related to the risks of the engineered nanoparticles and to assess the potential of efficient expansion and development of new approaches, which are offered by application of theoretical and computational methods, applicable for evaluation of nanomaterials.

Mixed poly-brominated/chlorinated biphenyls (PXBs): widespread food and environmental contaminants

Mixed poly-brominated/chlorinated biphenyls (PXBs) are a new class of emerging contaminants. Their environmental occurrence is confirmed by the reported occurrence in foods and human tissues, and the patterns of occurrence suggest that different sources may be contributing to food occurrence in different parts of the world. Important sources of PXBs are thought to include emissions from the combustion of bromine and chlorine containing waste and consumer products and inadvertent contamination in industrial chemicals. Of specific interest are the dioxin-like-PXBs. These are the 62 non-ortho-substituted PXBs with eight congeners of 3,3',4,4'-substitution type (#77A-H), ten of 3,4,4',5-type (#81A-J), twenty four of 3,4,5,3',4'-type (#126A-Z) and twenty of 3,4,5,3',4',5'-type (#169A-V) and many mono-ortho substituted compounds. The toxicological clarification on these contaminants continues, with reports confirming dioxin-like effects at low concentrations, suggesting a greater toxicological significance than PCBs.

QSPR models of n-octanol/water partition coefficients and aqueous solubility of halogenated methyl-phenyl ethers by DFT method

The possible molecular geometries of 134 halogenated methyl-phenyl ethers were optimized at B3LYP/6-31G(*) level with Gaussian 98 program. The calculated structural parameters were taken as theoretical descriptors to establish two new novel QSPR models for predicting aqueous solubility (-lgS(w,l)) and n-octanol/water partition coefficient (lgK(ow)) of halogenated methyl-phenyl ethers. The two models achieved in this work both contain three variables: energy of the lowest unoccupied molecular orbital (E(LUMO)), most positive atomic partial charge in molecule (q(+)), and quadrupole moment (Q(yy) or Q(zz)), of which R values are 0.992 and 0.970 respectively, their standard errors of estimate in modeling (SD) are 0.132 and 0.178, respectively. The results of leave-one-out (LOO) cross-validation for training set and validation with external test sets both show that the models obtained exhibited optimum stability and good predictive power. We suggests that two QSPR models derived here can be used to predict S(w,l) and K(ow) accurately for non-tested halogenated methyl-phenyl ethers congeners.

The aqueous solubility of some herbicidal by-side toxic impurities: predicted data of the 399 chlorinated trans-azoxybenzene congeners

The quantitative structure - property relationship (QSPR) and the artificial neural networks (ANNs) methods were used to estimate aqueous solubility (log S and μg/L) of polychlorinated trans-azoxybenzenes (PCt-ABs). These QSPR and ANN models are based on geometry optimalization and quantum-chemical structural descriptors, which were computed on the level of density functional theory (DFT) using B3LYP functional and 6-311++G** basis set in Gaussian 03 software and the semi-empirical quantum chemistry method for property parameterization (RM1) in the molecular orbital package (MOPAC) software. The predicted solubility of PCt-AOBs by RM1 and DFT models and depending on a congener varied within a homologue class between 47-19498 and 371-1738 μg/L for Mono-; 33-11481 and 7.9-3630 μg/L for Di-; 6.1-4786 and 4.7-12882 μg/L for Tri-; 1.3-1174 and 0.3-14791 μg/L for Tetra-; 0.4-646 and 0.1-38904 μg/L for Penta-; 0.1-155 and 0.2-63096 μg/L for Hexa-; 0.2-27 and 0.1-646 μg/L for Hepta-; < 0.1-6.2 and 0.8-282 μg/L for Octa-; 0.6-2.6 and 0.8-12 μg/L for NonaCt-AOBs; and 1.2 and 0.5 μg/L for DecaCt-AOB, respectively. Both computational models used were characterized by good predictive abilities and small errors, while calculations by RM1 method were highly competitive compared to a much more time-consuming and expensive method by DFT.

Estimation of K(OA) values of 209 polychlorinated trans-azobenzenes by PM6 and DFT methods

The octanol-air partition coefficients (K(OA)) of all 209 PCt-ABs were determined computationally to fill gaps on their environmentally relevant physical and chemical properties. These properties have been determined using two computational approaches: the semi-empirical quantum chemistry method for property parameterization (PM6) of the molecular orbital package (MOPAC) and density functional theory (DFT) method using B3LYP functional and 6-311++G** basis set in Gaussian 03 software and artificial neural network (ANN) predicting abilities. Both computational methods enabled estimation of log K(OA) partition coefficients of PCt-ABs with a similar accuracy and precision. The PM6 method compared to DFT was highly superior because it requires much less time, manpower and cost of hardware. The determined log K(OA) values of the investigated PCt-ABs for standard condition (25 °C) varied between 8.30 and 8.75 for Mono-; 8.71 and 9.92 for Di-; 9.58 and 10.72 for Tri-; 10.11 and 11.34 for Tetra-, 10.83 and 11.85 for Penta-; 11.24 and 12.36 for Hexa-; 11.87 and 12.66 for Hepta-; 12.31 and 12.97 for Octa-; 12.89 and 13.21 for Nona-Ct-ABs; and 13.17- and 13.49 for Deca-Ct-AB. PCt-ABs, in view of these log K(OA) values, can be classified as compounds of relatively low (Mono-, Di- and some of Tri- Ct-ABs with values of log K(OA) around 8 to 10) environmental mobility (most of Tri- to Nona-Ct-ABs and Deca-Ct-AB homologues with values of log K(OA) >10), and with a potential to be adsorbed by soil particles.

Use of quantitative-structure property relationship (QSPR) and artificial neural network (ANN) based approaches for estimating the octanol-water partition coefficients of the 209 chlorinated trans-azobenzene congeners

Polychlorinated azobenzenes (PCABs) can be found as contaminant by products in 3,4-dichloroaniline and its derivatives and in the herbicides Diuron, Linuron, Methazole, Neburon, Propanil and SWEP. Trans congeners of PCABs are physically and chemically more stable and so are environmentally relevant, when compared to unstable cis congeners. In this study, to fulfill gaps on environmentally relevant partitioning properties of PCABs, the values of n-octanol/water partition coefficients (log K(OW)) have been determined for 209 congeners of chloro-trans-azobenzene (Ct-AB) by means of quantitative structure-property relationship (QSPR) approach and artificial neural networks (ANN) predictive ability. The QSPR methods used based on geometry optimalization and quantum-chemical structural descriptors, which were computed on the level of density functional theory (DFT) using B3LYP functional and 6-311++G basis set in Gaussian 03 and of the semi-empirical quantum chemistry method (PM6) of the molecular orbital package (MOPAC). Polychlorinated dibenzo-p-dioxins (PCDDs), -furans (PCDFs) and -biphenyls (PCBs), to which PCABs are related, were reference compounds in this study. An experimentally obtained data on physical and chemical properties of PCDD/Fs and PCBs were reference data for ANN predictions of log K(OW) values of Ct-ABs in this study. Both calculation methods gave similar results in term of absolute log K(OW) values, while the models generated by PM6 are considered highly efficient in time spent, when compared to these by DFT. The estimated log K(OW) values of 209 Ct-ABs varied between 5.22-5.57 and 5.45-5.60 for Mono-, 5.56-6.00 and 5.59-6.07 for Di-, 5.89-6.56 and 5.91-6.46 for Tri-, 6.10-7.05 and 6.13-6.80 for Tetra-, 6.43-7.39 and 6.48-7.14 for Penta-, 6.61-7.78 and 6.98-7.42 for Hexa-, 7.41-7.94 and 7.34-7.86 for Hepta-, 7.99-8.17 and 7.72-8.20 for Octa-, 8.35-8.42 and 8.10-8.62 for NonaCt-ABs, and 8.52-8.60 and 8.81-8.83 for DecaCt-AB. These log K(OW) values shows that Ct-ABs are compounds of relatively low environmental mobility (log K(OW) > 4.5) and of significant bioaccumulation potential.

QSPR for prediction of subcooled vapor pressures (log PL) of polychlorinated trans-azobenzenes

In this study the values of subcooled vapor pressures (log P(L)) were estimated for 209 trans chloroazobenzenes (Ct-ABs) that fill some gaps in analytical and experimental data on these compounds. There are 209 chloro derivatives of trans azobenzenes that are relatively stable and more environmentally relevant than 209 chloro cis congeners. The calculations models were based on the Quantitative Structure-Property Relationship (QSPR) scheme using the semi-empirical method (PM6) in molecular package (MOPAC) software and density functional theory (DFT) method using B3LYP functional and 6-311++G** basis set) in Gaussian 03 software method and the artificial neural networks (ANNs) prediction. The values of log P(L) predicted by models used varied between -3.94 to -2.66 for Mono-; -4.85 to -2.97 for Di-; -5.18 to -3.17 for Tri-; -6.02 to -3.77 for Tetra-; -6.64 to -4.64 for Penta-; -7.36 to -4.76 for Hexa-; -7.54 to -5.79 for Hepta-; -7.75 to -6.64 for Octa-; -7.89 to -7.44 for Nona-Ct-Abs; and -8.09 and -8.13 for Deca-Ct-AB. Based on these values Ct-ABs can be grouped localized among relatively low (log P(L) -4 to -2) and low (log P(L) < -4) mobile Persistent Organic Pollutants (POPs). Both the calculation methods employed were characterized by similar prediction ability of subcooled vapor pressure values of Ct-ABs, while those of PM6 are much more efficient due to a cheaper hardware used and around 300-fold less time spent on calculations.

Prediction of subcooled vapor pressures (log PL) of 399 polychlorinated trans-azoxybenzenes by using the QSPR and ANN approach

Environmentally relevant partitioning properties such as the sub-cooled vapor pressures (log PL) have been predicted for 399 congeners of chloro-trans-azoxybenzene (C-t-AOBs) by two computational methods. The quantitative structure-property relationship (QSPR), an approach which is based on geometry optimalization and quantum-chemical structural descriptors in RM1 and DFT methods and artificial neural networks (ANNs), an approach that predicts abilities that give similar results of estimated log P(L) and the accuracy of the methods was also similar. The RM1 method was less time consuming and less costly compared to calculations by the DFT method. Estimated from the RM1 and DFT methods of log P(L) values of 399 Ct-AOBs varied between -1.98 to -0.93 and -1.83 to -0.79 for Mono-, 3.12 to -1.46 and -3.00 to -1.46 for Di-, -4.03 to -1.39 and -3.53 to -1.67 for Tri-, -4.75 to -2.33 and -4.59 to -1.91 for Tetra-, -5.37 to -2.59 and -5.42 to -2.09 for Penta-, -5.82 to -2.88 and -5.66 to -2.58 for Hexa-, -5.88 to -3.24 and -5.60 to -2.93 for Hepta-, -6.28 to -4.33 and -5.60 to -4.29 for Octa-, -6.54 to -5.28 and -5.66 to -4.93 for NonaCt-AOBs, and -6.59 and -5.61 for DecaCt-AOB. According to a common classification of environmental contaminants and by sub-cooled vapor pressure values, MonoCt-AOBs and a few of the Di- and TriCt-AOBs (log P(L)from -2 to 0) fall into the group of compounds that are relatively well mobile in the ambient environment, while most of the Di- to HeptaCt-AOBs (log P(L) < -4 to -2) mobility is relatively weak. Octa- and NonaCt-AOBs and DecaCt-AOB (log P(L) < -4) are also weak mobile contaminants.

QSPR models for prediction of the soil sorption coefficient (log KOC) values of 209 polychlorinated trans-azobenzenes (PCt-ABs)

The values of the soil sorption coefficient (K(OC)) have been computed for 209 environmentally relevant trans polychlorinated azobenzenes (PCABs) lacking experimental partitioning data. The quantitative structure-property relationship (QSPR) approach and artificial neural networks (ANN) predictive ability used in models based on geometry optimalization and quantum-chemical structural descriptors, which were computed on the level of density functional theory (DFT) using B3LYP functional and 6-311++G** basis set and of the semi-empirical quantum chemistry method for property parameterization (PM6) of the molecular orbital package (MOPAC). An experimentally available data on physical and chemical properties of PCDD/Fs and PCBs were used as reference data for the QSPR models and ANNs predictions in this study. Both calculation methods gave similar results in term of absolute log K(OC) values, while the PM6 model generated in the MOPAC was a much more efficient compared to the DFT model in GAUSSIAN. The estimated values of log K(OC) varied between 4.93 and 5.62 for mono-, 5.27 and 7.46 for di-, 6.46 and 8.09 for tri-, 6.65 and 9.11 for tetra-, 6.75 and 9.68 for penta-, 6.44 and 10.24 for hexa-, 7.00 and 10.36 for hepta-, 7.09 and 9.82 octa-, 8.94 and 9.71 for nona-Ct-ABs, and 9.26 and 9.34 for deca-Ct-AB. Because of high log K(OC) values PCt-ABs could be classified as compounds with high affinity to the particles of soil, sediments and organic matter.

QSAR and ANN for the estimation of water solubility of 209 polychlorinated trans-azobenzenes

Polychlorinated trans-azobenzenes (PCt-ABs) are less studied a highly toxic impurity in 2,3-dichloroaniline (2,3-D) and some herbicides and are compounds of environmental relevance lacking experimental physical and chemical properties data. In this study, to fulfill gaps on environmentally relevant partitioning properties of PCABs, the values of water solubility (μg/L and log S) have been determined for 209 congeners of chloro-trans-azobenzene (Ct-AB) by means of quantitative structure - property relationship (QSPR) approach and artificial neural networks (ANN) predictive ability. The quantitative structure - property relationship (QSPR) approach used based on geometry optimalization and quantum-chemical structural descriptors, which were computed on the level of density functional theory (DFT) using B3LYP functional and 6-311++G** basis set in Gaussian 03 and the semi-empirical quantum chemistry method for property parameterization (PM6) in the molecular orbital package (MOPAC) software. The predicted solubility of PCt-ABs by PM6 and DFT models and depending on a congener within a homologue class varied between 1995-11481 and 5370-15135 μg/L for mono-; 170-5495 and 138-9332 μg/L for di-; 36-1950 and 209-5248 μg/L for tri-; 15-794 and 41-3715 μg/L for tetra-; 5.5-209 and 39-1259 μg/L for penta-; 1.8-98 and 3.5-1096 μg/L for hexa-; 1.5-34 and 4.7-214 μg/L for hepta-; 0.71-6.2 and 0.76-26 μg/L for octa-; 0.83-1.7 and 0.69-1.2 μg/L for nonaCt-ABs; and between 0.36 and 0.04 μg/L for decaCt-AB, respectively. The calculations by PM6 were highly efficient and inexpensive compared to these by DFT, while both models gave data of similar accuracy.

N-octanol-water partition coefficients (log K(OW)) of 399 congeners of polychlorinated azoxybenzenes (PCAOBs) determined by QSPR- and ANN-based approach

Polychlorinated azoxybenzenes (PCAOBs) theoretically consist of 798 congeners with 399 in cis and 399 in trans configuration. PCAOBs in trans configuration are largerly planar compounds and some are highly toxic and environmentally relevant compared to cis congeners. Trans-PCAOBs can be found as by-side products in 3,4-dichloroaniline and some herbicides. To fulfill gaps in physical and chemical properties of PCAOBs, the values of log K(OW) were determined for 399 congeners of t-CAOB using a computational approach. We used the semi-empirical RM1 in MOPAC and DFT B3LYP in Gaussian 03 methods, artificial neural net (ANN) predictions, and the standardized variables with and without the normal varimax rotation. The models created predicted the values of log K(OW) of all 399 chlorinated derivatives of trans-azoxybenzenes (C-t-AOBs). The values of log K(OW) of C-t-AOBs varied between 5.08 and 5.42 for Mono-, 5.16 and 5.96 for Di-, 5.79 and 6.73 for Tri-, 6.26 and 7.18 for Tetra-, 6.65 and 7.54 for Penta-, 7.13 and 7.94 for Hexa-, 7.20 and 8.20 for Hepta-, 7.96 and 8.32 for Octa-, 8.32 and 8.43 for Nonachloro-t-AOBs and 8.55 and 8.97 for Decachloro-t-AOB. These log K(OW) values were similar per chloro-t-AOB congener and independent of the calculation method. C-t-AOBs have log K(OW) values above 4.5, and what relates to contaminants of low or very low environmentally mobility but a high predilection to the soil and sediment particles and with potential for bioaccumulation. The models that used the standardized variables had smallest errors and higher correlation coefficients compared to the models that based on the normal varimax rotation of standardized structural descriptors. In light of these data, the semi-empirical RM1 calculations in MOPAC software and followed by ANN were a much less time consuming and less expensive compared to the DFT B3LYP method.

Modeling the overall persistence and environmental mobility of sulfur-containing polychlorinated organic compounds

BACKGROUND, AIM, AND SCOPE

Experimental data on partition coefficients and environmental half-lives of sulfur analogs of polychlorinated organic compounds are scarce. Consequently, little is known about their overall persistence and long-range transport potential, which are the most vital measures in the environmental exposure assessment. We performed Multimedia Modeling of environmental fate and transport to complement this paucity of scientific data. The main aim of our study was to investigate whether the sulfur analogs of polychlorinated dibenzo-p-dioxins, -dibenzofurans, and -diphenylethers are as environmentally persistent and/or mobile as their oxygen counterparts and to propose the environmental exposure-related classification of the examined sulfur compounds.

MATERIALS AND METHODS

Our study included all possible congeners of the sulfur analogs generated in a combinatorial approach. We predicted (1) lacking data on partition coefficients (log K OW, log K OA and log K AW) for oxygen and sulfur analogs using Quantitative Structure-Property Relationship (QSPR) modeling and (2) their half-lives in air, water, and soil using US EPA tool 'The PBT Profiler, v. 1.203 2006'. Subsequently, we introduced these results into multimedia mass balance model 'The OECD POV and LRTP Screening Tool, v. 2.2'.

RESULTS

Our study revealed that log K OW and log K OA are increasing by constant values of 0.60 and 1.07, respectively, and the values of log K AW are decreasing by 0.90, whenever one oxygen atom in the carbon skeleton is replaced by sulfur. The persistence ranking performed by the PBT Profiler showed that PCDDs, PCDFs, PCDEs, and their sulfur analogs belong to one half-life class.

DISCUSSION

The Multimedia Modeling by the means of 'The OECD POV and LRTP Screening Tool, v. 2.2' suggested that the long-range transport potential depends on the presence/absence of oxygen/sulfur atoms in particular molecules, their substitution pattern and the parent carbon skeleton. Sulfur analogs are generally less mobile than their oxygen analogs, but have similar overall persistence and much higher bioaccumulation potential. Thus, according to the classification of chemicals proposed by Klasmeier et al. (Environ Sci Technol 40:53-60, 2006), some of them show POP-like POV and LRTP characteristics while the rest shows POP-like P OV characteristics.

CONCLUSIONS

The sulfur analogs of PCDDs, PCDFs, or PCDEs bring environmental mobility comparable with the risk related to the oxygen ones; they belong to the pollutants of 'highest' or 'intermediate' priority.

RECOMMENDATIONS AND PERSPECTIVES

Further studies that would verify the necessity to include the studied sulfur molecules in the international lists of high-priority environmental pollutants are recommended.

Toward the development of "nano-QSARs": advances and challenges

The most significant achievements and challenges relating to an application of quantitative structure-activity relationship (QSAR) approach in the risk assessment of nanometer-sized materials are highlighted. Recent advances are discussed in the context of "classical" QSAR methodology. The possible ways for the structural characterization of compounds existing at the nanoscale (at least one dimension of 100 nm or less) are briefly reviewed. The applicability of the existing toxicological data for developing QSAR models is evaluated. Finally, the existing models are presented. The need to develop new interpretative descriptors for the nanosystems is also highlighted. It is suggested that, due to high variability in the molecular structures and different mechanisms of toxicity, individual classes of nanoparticles should be modeled separately.

Prediction of supercooled liquid vapor pressures and n-octanol/air partition coefficients for polybrominated diphenyl ethers by means of molecular descriptors from DFT method

The molecular geometries of 209 polybrominated diphenyl ethers (PBDEs) were optimized at the B3LYP/6-31G level with Gaussian 98 program. The calculated structural parameters were taken as theoretical descriptors to establish two novel QSPR models for predicting supercooled liquid vapor pressures (P(L)) and octanol/air partition coefficients (K(OA)) of PBDEs based on the theoretical linear solvation energy relationship (TLSER) model, respectively. The two models achieved in this work both contain three variables: most negative atomic partial charge in molecule (q(-)), dipole moment of the molecules (mu) and mean molecular polarizability (alpha), of which R(2) values are both as high as 0.997, their root-mean-square errors in modeling (RSMEE) are 0.069 and 0.062 respectively. In addition, the F-value of two models are both evidently larger than critical values F(0.05) and the variation inflation factors (VIF) of variables herein are all less than 5.0, suggesting obvious statistic significance of the P(L) and K(OA) predicting models. The results of Leave-One-Out (LOO) cross-validation for training set and validation with external test set both show that the two models obtained exhibited optimum stability and good predictive power. We suggest that the QSPRs derived here can be used to predict accurately P(L) and K(OA) for non-tested PBDE congeners from Mono-BDEs to Hepta-BDEs and from Mono-BDEs to Hexa-BDEs, respectively.

Application and comparison of different chemometric approaches in QSPR modelling of supercooled liquid vapour pressures for chloronaphthalenes

Molecular descriptors from calculations at the level of Density Functional Theory (B3LYP/6-311++G**) were effectively applied in QSPR estimation of supercooled liquid vapour pressures (P(L)) for individual chloronaphthalene congeners. The estimated values of log P(L) varied from 1.05 Pa to 5.6 x 10(-5) Pa, depending on the number of chlorine substituents present in the molecule and the substitution pattern. Comparison of the five chemometrical methods of modelling (approaches) led to the final conclusion, that the use of relatively simple PLS combined with one of the variable pre-selection algorithms (UVE or GA) seems to be the optimal choice in such computational studies for persistent organic pollutants. The best GA-PLS model was characterized by the value of root mean square error of prediction RMSEP = 0.108 logarithmic Pascal units.