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

Fate and transport of monoterpenes through soils. Part I. Prediction of temperature dependent soil fate model input-parameters

Monoterpenes are C10H(n)O(n') compounds of natural origin and are potentially environmentally safe substitutes for traditional pesticides. Still, an assessment of their environmental behaviour is required. As a first step in a theoretical study focussing on monoterpenes applied as pesticides to terrestrial environments, soil fate model input-parameters were determined for 20 monoterpenes with widely different structural characteristics. Input-parameters are the water solubility (S(W)), vapour pressure (P), n-octanol-water partition coefficient (K(OW)), atmospheric air and bulk water diffusion coefficients (D(A)air and D(W)water), first order biodegradation rate constants (k), and their temperature dependence. Values for these parameters were estimated or taken from previous experimental work. The quality of the estimations was discussed by focussing on their statistics and by comparison with available experimental data. From these properties, the air-water partition coefficient (K(AW), Henry's Law constant), the interface-water partition coefficient (K(IW)) and the organic matter-water partition coefficient (K(OM)) could be estimated with varying levels of accuracy. In general, little experimental data turned out to be available on biodegradation rate constants and on the temperature dependence of physico-chemical parameters.     

Vapour pressures, aqueous solubilities, Henry's law constants, partition coefficients between gas/water (Kgw), n-octanol/water (Kow) and gas/n-octanol (Kgo) of 106 polychlorinated diphenyl ethers (PCDE)

Modelling the environmental fate of persistent organic pollutants like polychlorinated diphenyl ethers (PCDE) requires the knowledge of a number of fundamental physico-chemical properties of these compounds. We report here the physico-chemical properties of 106 PCDEs, which are over 50% of all possible congeners. Vapour pressures P(OL), water solubilities S(H2O), and n-octanol/water partition coefficients K(OW) were determined with chromatographic methods. With these experimental data the Henry's law constants H, gas/water K(GW) and gas/n-octanol K(GO) partition coefficients were calculated. Vapour pressures and water solubilities and n-octanol/water partition coefficients of the PCDEs are close to those of similar groups of organochlorine compounds like polychlorinated biphenyls (PCBs) and dibenzofurans (PCDFs). A similar environmental fate can be predicted and was partially already been observed.     

Estimation of selected physicochemical properties for methylated naphthalene compounds

Liquid aqueous solubility (S(w,L)), octanol/water partition coefficients (K(ow)), liquid vapor pressure (P(v,L)), and Henry's law constants (H(c)) were estimated for 20 methylated naphthalenes ranging from monomethyl to tetramethylnaphthalenes. Chromatographic methods were used for the estimation. Chromatographic retention measurements were conducted for 11 reference compounds and regressions were fit between the retention indices and the physicochemical properties. HPLC octadecylsilyl column with acetonitrile/water eluent was used for the estimation of S(w,L) and K(ow). Two GC columns, HP5-MS and a more hydrophobic HP-1, were tested for the estimation of P(v,L). Measured retention indices for the methylated naphthalenes were entered to the regression equations to calculate the physicochemical properties for these compounds. Literature values, where available, were used to validate the calculated values. The method accurately estimated the physicochemical properties. Estimated S(w,L) and P(v,L) decreased with the number of methyl groups. K(ow) increased with the number of methyl groups. There was no obvious relation between H(c) and the number of methyl groups. Log S(w,L) ranged from 0.885 for 1,2,5,6-tetramethylnaphthalene to 2.269 for 1-methylnaphthalene (mmol/m(3)). Log K(ow) varied from 3.89 for 1-methylnaphthalene to 4.95 for 1,2,5,6-tetramethylnaphthalene. Log P(v,L) ranged from -0.983 for 1,2,5,6-tetramethylnaphthalene to 0.789 for 2-methylnaphthalene (Pa). Log H(c) varied from 1.03 for 1,4,5-trimethylnaphthalene to 1.73 for 2,6-dimethylnaphthalene (Pa m(3)/mol). There were no apparent effects of GC column hydrophobicity on the accuracy of the results. Estimation of S(w,L) and K(ow) based on GC retention indices was not as accurate as with HPLC. Comparison of the estimated values with values predicted by EPIWIN indicated that EPIWIN is useful in giving order-of-magnitude prediction of physicochemical properties.     

Physicochemical properties of selected polybrominated diphenyl ethers and extension of the UNIFAC model to brominated aromatic compounds

The aqueous solubilities (S(w)) at various temperatures from 283 K to 308 K and 1-octanol/water partition coefficients (K(ow)) for four polybrominated diphenyl ethers (PBDEs: 4,4'-dibromodiphenyl ether (BDE-15), 2,2',4,4'-tetrabromodiphenyl ether (BDE-47), 2,2',4,4',5-pentabromodiphenyl ether (BDE-99), and 2,2',4,4',5,5'-hexabromodiphenyl ether (BDE-153)) were measured by the generator column method. The S(w) and K(ow) data revealed the effect of bromine substitution and basic structure on S(w) and K(ow). To estimate the infinite dilution activity coefficients (gamma(i)(w,infinity)) of the PBDEs in water from the S(w) data, enthalpies of fusion and melting points for those compounds were measured with a differential scanning calorimeter. Henry's Law constants (H(w)) of the PBDEs were derived from the determined gamma(i)(w,infinity) and literature vapor pressure data. Some physicochemical characteristics of PBDEs were also suggested by comparing the present property data with that of polychlorinated dibenzo-p-dioxins, brominated phenols and brominated benzenes in past studies. Furthermore, in order to represent different phase equilibria including solubility and partition equilibrium for other brominated aromatic compounds using the UNIFAC model, a pair of UNIFAC group interaction parameters between the bromine and water group were determined from the S(w) and K(ow) data of PBDEs and brominated benzenes. The ability of the determined parameters to represent both properties of brominated aromatics was evaluated.     

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

Temperature dependent properties of environmentally important synthetic musks

Environmental fate determining physical properties including their temperature dependence for five nitro musks and for seven polycyclic musks were estimated. The properties evaluated were vapor pressure in a solid and liquid state (PS and PL), solubility in water (S), Henry's law coefficient (H = PL/S) and log octanol-water partition coefficient (log KOW). Gas chromatography for starting values of vapor pressure estimation and HPLC experiments at 5-20 degrees C for comparison of the theoretical versus experimental solubilities in water were performed. The values of temperature (T) dependence coefficients (Ai and Bi) in equations: log (Property)i = Ai - Bi/T were determined. Values of properties were compared with literature-based data, and an example of their use in environmental hazard estimation by fate modeling was given.     

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.     

Predicting organic contaminant concentrations in sediment porewater using solid-phase microextraction

Because of its cost and time saving features, solid-phase microextraction (SPME) is a leading candidate as a biomimic technique in assessing the bioavailable fraction of hydrophobic organic contaminants (HOCs) in sediment porewater. However, no predictive modeling framework in which to systematically address the effect of key parameters on SPME performance for this application exists. In this study, we derived two governing equations to predict (1) the minimum sediment volume (V(s)min) required to achieve non-depletive conditions, and (2) dissolved phase HOC porewater concentrations (C(pw)) as functions of HOC- and sediment specific characteristics in a conceptual three compartment system. The resulting model predicted that V(s)min was independent of HOC concentrations both in sediment and porewater, but did vary with hydrophobicity (characterized by logK(ow)), the fraction of sediment porewater (f(pw)), and the volume (V(f)) of the SPME sorbent phase. Moreover, the effects of these parameters were minimized (i.e., V(s)min reached plateaus) as logK(ow) approached 4-5. Model predictions of C(pw), a surrogate for SPME-based detection limits in porewater, decreased with increasing sediment volume (V(s)) at low V(s) values, but rapidly leveled off as V(s) increased. A third result suggested that the sediment HOC concentration required for SPME is completely independent of K(ow). These results suggest that relatively small sediment volumes participate in exchange equilibria among sediment, porewater and the SPME fiber, and that large sediment HOC reservoirs are not needed to improve the detection sensitivity of SPME-based porewater samplers. The ultimate utility of this modeling framework will be to assist future experimental designs and help predict in situ bioavailability of sediment-associated HOCs.     

Temperature dependencies of Henry's law constants and octanol/water partition coefficients for key plant volatile monoterpenoids

To model the emission dynamics and changes in fractional composition of monoterpenoids from plant leaves, temperature dependencies of equilibrium coefficients must be known. Henry's law constants (H(pc), Pa m3 mol(-1) and octanol/water partition coefficients (K(OW), mol mol(-1)) were determined for 10 important plant monoterpenes at physiological temperature ranges (25-50 degrees C for H(pc) and 20-50 degrees C for K(OW)). A standard EPICS procedure was established to determine H(pc) and a shake flask method was used for the measurements of K(OW). The enthalpy of volatilization (deltaH(vol)) varied from 18.0 to 44.3 kJ mol(-1) among the monoterpenes, corresponding to a range of temperature-dependent increase in H(pc) between 1.3- and 1.8-fold per 10 degrees C rise in temperature. The enthalpy of water-octanol phase change varied from -11.0 to -23.8 kJ mol(-1), corresponding to a decrease of K(OW) between 1.15- and 1.32-fold per 10 degrees C increase in temperature. Correlations among physico-chemical characteristics of a wide range of monoterpenes were analyzed to seek the ways of derivation of H(pc) and K(OW) values from other monoterpene physico-chemical characteristics. H(pc) was strongly correlated with monoterpene saturated vapor pressure (P(v)), and for lipophilic monoterpenes, deltaH(vol) scaled positively with the enthalpy of vaporization that characterizes the temperature dependence of P(v) Thus, P(v) versus temperature relations may be employed to derive the temperature relations of H(pc) for these monoterpenes. These data collectively indicate that monoterpene differences in H(pc) and K(OW) temperature relations can importantly modify monoterpene emissions from and deposition on plant leaves.     

Estimating K(oc) for persistent organic pollutants: limitations of correlations with K(ow)

The n-octanol/water partition coefficient (K(ow)) is commonly used to predict the soil or aquatic particle water partition coefficient normalized to organic carbon (K(oc)). Many correlations are available covering several chemical classes and ranges of hydrophobicity. This work indicates the K(ow) may not be a strong predictor for persistent organic pollutants (POPs) which are defined here as chemicals with logK(ow) > 5.0. In addition, the correlation developed in this work for POPs will still result in a predicted value which is of by a factor of 15. Accordingly, care must be taken when applying K(oc) estimations using K(ow) for POPs until more suitable correlations are developed.     

Correlating the physical-chemical properties of phthalate esters using the 'three solubility' approach

A quantitative structure-property relationship (QSPR) method for the correlation of physical-chemical properties and partition coefficients, namely the 'three solubility' approach, is described and applied to a group of 22 phthalate esters. The solubilities or 'apparent-solubilities' of these substances in the liquid state are compiled and correlated against Le Bas molar volume in the three primary media of air, water and octanol. From these solubilities the air-water (K(AW)), octanol-water (K(OW)) and octanol-air (K(OA)) partition coefficients are deduced. Estimated solubilities in water and octanol-water partition coefficients are shown to compare favourably with more recent accurate measurements. A set of selected values is presented, with error limits, which is recommended for use in modelling and assessment studies. Some environmental implications are discussed of the large range in property values for this series.     

Soot-water distribution coefficients for polychlorinated dibenzo-p-dioxins,polychlorinated dibenzofurans and polybrominated diphenylethers determined with the soot cosolvency-column method

For many types of hydrophobic compounds, sorption non-linearity and solid-water distributions in the field well above expectations from organic matter partitioning models have lead to the proposition that strong adsorption to soot surfaces may not be limited to polycyclic aromatic hydrocarbons but may extend as a significant process for many aromatic compound classes. Here, the soot-water distribution coefficients (Ksc) were determined with the soot cosolvency-column method for homolog series of five polychlorinated dibenzo-p-dioxins (PCDDs), five polychlorinated dibenzofurans (PCDFs) and for two polybrominated diphenylethers (PBDEs). All compounds exhibited significantly stronger association with soot carbon than expected from estimates of their bulk organic-carbon normalized partition coefficients (Koc). The Ksc/Koc ratios (at aqueous concentrations of around 0.1-1 microg/l) were for PCDDs (up to tetrachlorination) 19-130 (median 25), for PCDFs (also up to tetrachlorination) 150-490 (median 300), and for both the tetra- and pentabrominated PBDEs a factor of 60. The particularly strong soot sorption for the PCDFs is of similar enhancement factors as previously elucidated for polycyclic aromatic hydrocarbons. Compound-class specific correlations between log Ksc and octanol-water partition coefficients (log Kow) were significant for both PCDDs and PCDFs (and with R2 > 98%). These may prove useful for anticipating variable fractions of dissolved exposures between different environmental regimes and putative remediation objects.     

Estimation of the volatilization of organic compounds from soil surfaces

Several simple models for the estimation of the half-life (t(1/2)) for the depletion of an organic chemical from a soil surface to air were examined. For moist surfaces, two models are proposed: the first requires knowledge of the soil/organic carbon partition coefficient (K(oc)) and the Henry's law constant (H) and the second the vapor pressure (P(s)) of the chemical involved. Due to uncertainties in the experimental K(oc) values those ones predicted by the group-contribution model of Meylan et al. [Environ. Sci. Technol. 26 (1992) 1560]-and proposed by the U.S. Environmental Protection Agency (EPA)-should be used. If reliable experimental P(s) values are not available, the first model is proposed, where in cases when H values are not available, predicted ones by the Bond-Contribution method of Meylan and Howard [Environ. Toxicol. Chem. 10 (1991) 1283]-and also proposed by EPA-can be used. In general, the agreement of the predicted t(1/2) values with the measured ones is within a factor of 3-5. Similar expressions, but with somewhat poorer results, are presented for dry field soils. In all cases, the obtained results represent a substantial improvement over those obtained with the currently used Dow method: t(1/2) = 1.58 x 10(-8)((K(oc) x S)/P(S)), where S is the solubility of the compound in water.     

Evaluating coal tar-water partitioning coefficient estimation methods and solute-solvent molecular interactions in tar phase

Equilibrium partitioning coefficients between an industrial coal tar sample and water (K(CT/w)) were determined for 41 polar and nonpolar solutes in batch systems. Together with literature values, 69 K(CT/w) data were analyzed using the following model approaches: Raoult's law, the single parameter linear free energy relationship (SPLFER) with octanol-water partitioning coefficients (K(ow)), the linear solvation energy relationships (LSERs), SPARC and COSMOtherm. Estimations by Raoult's law and the SPLFER agreed well with the experimental log K(CT/w) values for the investigated coal tar, with root mean square errors (RMSE) of 0.31 and 0.33, respectively. LSER resulted in as good estimations (RMSE=0.29) as the previous two. The LSER analysis revealed significant hydrogen (H)-bond acceptor properties of the studied coal tar phase. Using naphthalene as a surrogate solvent for the coal tar phase, SPARC and COSMOtherm provided fairly good predictions (RMSE of 0.63 and 0.65, respectively) of log K(CT/w), without any additional empirical parameter. Further calculations using SPARC and COSMOtherm for partitioning between water and other tar-components (e.g., benzofuran, phenol and quinoline) suggested that minor components in coal tar do not significantly influence K(CT/w) of nonpolar solutes, and that Raoult's law and the SPLFER thus may be generally applied to these types of solutes, e.g., polycyclic aromatic hydrocarbons and alkylbenzenes, regardless of coal tar compositions. In contrast, partitioning of H-bonding solutes (e.g., phenols) can significantly vary depending on the amount of polar tar-components such as N-heterocyclic aromatic compounds. Therefore, the presented successful applications of Raoult's law and SPLFER to the studied coal tar could be a special case, and these simple approaches may not provide reasonable estimations for partitioning of H-bonding solutes from compositionally different coal tars.     

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.     

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