Liposome and protein-water partitioning of polybrominated diphenyl ethers (PBDEs)

Mitigation of PBDE net discharge in hazardous waste thermal treatment system through reintroducion of sludge and fly ash into GASMILD operations

The presence of polybrominated diphenyl ethers (PBDEs) in consumer products, waste treatment processes, and treated ashes poses a significant environmental threat. Due to the lack of research on the removal of PBDEs during waste incineration, this study investigated the effectiveness of a Hazardous Waste Thermal Treatment System (HAWTTS) utilizing reburning of sludge and fly ash (SFA) with gasification-moderate or intense low-oxygen dilution (GASMILD) combustion for PBDE removal. The closed-loop treatment of sludge and ash within the HAWTTS provides a potential pathway for near-zero PBDE emissions. The GASMILD combustion addresses potential combustion issues associated with fly ash recirculation. The system achieved an impressive overall removal efficiency of 98.4% for PBDEs, with minimal stack emissions (2.45 ng/Nm³) and a negative net discharge rate (-1.02 μg/h). GASMILD combustion played a crucial role (92.7%-97.6% destruction) in addressing challenges associated with high-moisture feedstocks and SFA residues. Debromination of highly brominated PBDEs occurred within the incinerator, resulting in an increased proportion of lower brominated PBDEs in the bottom slag compared to the feedstock. Air Pollution Control Devices (APCDs) achieved a total PBDE removal efficiency of 74.4%. However, the hydrophobic nature of PBDEs limited removal efficiency in scrubbers (36.0%) and cyclonic demisters (37.86%). This study demonstrates that reintroducing SFA into the GASMILD combustion process offers an effective and environmentally sustainable strategy for reducing net PBDE levels in hazardous waste. This approach also provides additional benefits such as energy conservation, reduced carbon emissions, and lower operating costs associated with secondary treatment of thermally treated byproducts.

How error-prone bioaccumulation experiments affect the risk assessment of hydrophobic chemicals and what could be improved

Bioaccumulation is one of the three criteria for the PBT assessment of chemicals, where P stands for persistence, B for bioaccumulation, and T for toxicity, which is a cornerstone for the "Registration, Evaluation, Authorization, and Restriction of Chemicals" (REACH) in the EU. Registrants are required by REACH to submit data on bioaccumulation if the chemical is manufactured in and/or imported to the European Economic Area at more than 100 t/year. Most of the experimental bioaccumulation studies submitted were on the bioconcentration factor (BCF) and were conducted prior to 2012, before the OECD Test Guideline 305 on Bioaccumulation in Fish was updated. An analysis of the submitted data revealed that many of the experimental data, but also the data from QSARs and other calculation methods, underestimate the actual bioaccumulation potential of hydrophobic substances considerably. One of the main reasons in the nonexperimental studies is that the BCF is related there to the total concentration of the chemical in water and not to the dissolved chemical concentration. There is therefore an urgent need to reassess the bioaccumulation potential of the hydrophobic substances registered under REACH. Based on the model calculations in the present study, between 332 and 584 substances that are registered under REACH are likely to bioaccumulate in the aquatic environment-many more than have so far been identified in the B assessment. Integr Environ Assess Manag 2023;19:792-803. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Chemometers: an integrative tool for chemical assessment in multimedia environments

We propose novel chemometers - passive equilibrium samplers of, e.g., silicone - as an integrative tool for the assessment of hydrophobic organic compounds in multimedia environments. The traditional way of assessing levels of organic pollutants across different environmental compartments is to compare the chemical concentration normalized to the major sorptive phase in two or more media. These sorptive phases for hydrophobic organic compounds differ between compartments, e.g., lipids in biota and organic carbon in sediments. Hence, comparability across media can suffer due to differences in sorptive capacities, but also extraction protocols and bioavailability. Chemometers overcome these drawbacks; they are a common, universal and well-defined polymer reference phase for sampling of a large range of nonpolar organic pollutants in different matrices like biota, sediment and water. When bringing the chemometer into direct contact with the sample, the chemicals partition between the sample and the polymer until thermodynamic equilibrium partitioning is established. At equilibrium, the chemical concentrations in the chemometers can be determined and directly compared between media, e.g., between organisms of different trophic levels or inhabiting different areas, between organs within an organism or between biotic and abiotic compartments, amongst others. Chemometers hence allow expressing the data on a common basis, as the equilibrium partitioning concentrations in the polymer, circumventing normalizations. The approach is based on chemical activity rather than total concentrations, and as such, gives a measure of the "effective concentration" of a compound or a mixture. Furthermore, chemical activity is the main driver for partitioning, biouptake and toxicity. As an additional benefit, the extracts of the chemometers only require limited cleanup efforts, avoiding introduction of a bias between chemicals of different persistence, and can be submitted to both chemical analysis and/or bioanalytical profiling.

Estimating the Bioaccumulation Potential of Hydrophobic Ultraviolet Stabilizers Using Experimental Partitioning Properties

Although hydrophobic ultraviolet (UV) stabilizers are an emerging environmental concern because of their widespread occurrence, persistence, and bioaccumulation potential, experimental values of their partitioning properties required for risk assessment are scarce. In this study, n-octanol-water partition (K_ow) and lipid–water partition constants (K_lipw), which are key parameters for environmental risk assessment, were experimentally determined for five selected hydrophobic UV stabilizers (UV326, UV327, UV328, UV329, and UV531) based on third-phase partitioning among polydimethylsiloxane (PDMS), water, and n-octanol/lipid. The partition constants between PDMS and water (K_PDMSw), obtained using the dynamic permeation method were used to derive K_ow and K_lipw. The obtained log K_ow and log K_lipw values were in the ranges of 7.08–7.94 and 7.50–8.34, respectively, indicating that the UV stabilizers exhibited a high bioaccumulation potential in aquatic environments. The experimental K_ow and K_lipw values obtained in this study provide valuable information for the evaluation of the fate, distribution, bioavailability, and toxicity of the UV stabilizers in aquatic environments.

Leaching of PBDEs from microplastics under simulated gut conditions: Chemical diffusion and bioaccumulation

Considerable efforts on exposure assessment of microplastics (MPs) as an agent in transport of toxic contaminants have been performed in organisms. However, chemical diffusion of inherent hydrophobic organic contaminants from MPs under simulated gut conditions is poorly examined. The present study examined the transfer kinetics of polybrominated diphenyl ethers (PBDEs) from polystyrene (PS), acrylonitrile butadiene styrene (ABS), and polypropylene (PP) MPs under gut surfactants (sodium taurocholate) at two relevant body temperatures of marine organisms, and evaluated the importance of MP ingestion in bioaccumulation of PBDEs in lugworm by a biodynamic model. Diffusion coefficients of PBDEs range from 5.82 × 10^-23 to 7.96 × 10^-20 m² s^-1 in PS, 5.49 × 10^-23 to 3.45 × 10^-20 m² s^-1 in ABS, and 5.58 × 10^-21 to 5.79 × 10^-17 m² s^-1 in PP, with apparent activation energies in the range of 33-148 kJ mol^-1. The biota-plastic accumulation factors of PBDEs leached from these plastics range from 1.44 × 10^-8 to 7.15 × 10^-5. Although ingestion of MPs with the common size (>0.5 mm) showed the negligible contribution to bioaccumulation of PBDEs in lugworm, their contribution in PBDEs transfer can be increased with gradual breakdown of MPs.

Critical Membrane Concentration and Mass-Balance Model to Identify Baseline Cytotoxicity of Hydrophobic and Ionizable Organic Chemicals in Mammalian Cell Lines

All chemicals can interfere with cellular membranes and this leads to baseline toxicity, which is the minimal toxicity any chemical elicits. The critical membrane burden is constant for all chemicals; that is, the dosing concentrations to trigger baseline toxicity decrease with increasing hydrophobicity of the chemicals. Quantitative structure-activity relationships, based on hydrophobicity of chemicals, have been established to predict nominal concentrations causing baseline toxicity in human and mammalian cell lines. However, their applicability is limited to hydrophilic neutral compounds. To develop a prediction model that includes more hydrophobic and charged organic chemicals, a mass balance model was applied for mammalian cells (AREc32, AhR-CALUX, PPARγ-BLA, and SH-SY5Y) considering different bioassay conditions. The critical membrane burden for baseline toxicity was converted into nominal concentration causing 10% cytotoxicity by baseline toxicity (IC_10,baseline) using a mass balance model whose main chemical input parameter was the liposome-water partition constants (K_lip/w) for neutral chemicals or the speciation-corrected D_lip/w(pH 7.4) for ionizable chemicals plus the bioassay-specific protein, lipid, and water contents of cells and media. In these bioassay-specific models, log(1/IC_10,baseline) increased with increasing hydrophobicity, and the relationship started to level off at log D_lip/w around 2. The bioassay-specific models were applied to 392 chemicals covering a broad range of hydrophobicity and speciation. Comparing the predicted IC_10,baseline and experimental cytotoxicity IC₁₀, known baseline toxicants and many additional chemicals were identified as baseline toxicants, while the others were classified based on specificity of their modes of action in the four cell lines, confirming excess toxicity of some fungicides, antibiotics, and uncouplers. Given the similarity of the bioassay-specific models, we propose a generalized baseline-model for adherent human cell lines: log[1/IC_10,baseline (M)] = 1.23 + 4.97 × (1 - e^{-0.236 log D_lip/w}). The derived models for baseline toxicity may serve for specificity analysis in reporter gene and neurotoxicity assays as well as for planning the dosing for cell-based assays.

Kinetics of Equilibrium Passive Sampling of Organic Chemicals with Polymers in Diverse Mammalian Tissues

Equilibrium passive sampling employing polydimethylsiloxane (PDMS) as a sampling phase can be used for the extraction of complex mixtures of organic chemicals from lipid-rich biota. We extended the method to lean tissues and more hydrophilic chemicals by implementing a mass-balance model for partitioning between lipids, proteins, and water in tissues and by accelerating uptake kinetics with a custom-built stirrer that effectively decreased time to equilibrium to less than 8 days even for a homogenized liver tissue with an only 4% lipid content. The partition constants log K_lipid/PDMS between tissues and PDMS were derived from measured concentration in PDMS and the mass-balance model and were very similar for 40 neutral chemicals with octanol-water partition constants 1.4 < log K_ow < 8.7, that is, log K_lipid/PDMS of 1.26 (95% CI, 1.13-1.39) for the adipose tissue, 1.16 (1.00-1.33) for the liver, and 0.58 (0.42-0.73) for the brain. This conversion factor can be applied to interpret chemical analysis and in vitro bioassays after additionally accounting for a small fraction of coextracted lipids of <0.7% of the PDMS weight. PDMS is more widely applicable for passive sampling of mammalian tissues than previously thought, both, in terms of diversity of chemicals and the range of lipid contents of tissues and, therefore, an ideal method for human biomonitoring to be combined with in vitro bioassays.

Development and evaluation of two-parameter linear free energy models for the prediction of human skin permeability coefficient of neutral organic chemicals

The experimental values of skin permeability coefficients, required for dermal exposure assessment, are not readily available for many chemicals. The existing estimation approaches are either less accurate or require many parameters that are not readily available. Furthermore, current estimation methods are not easy to apply to complex environmental mixtures. We present two models to estimate the skin permeability coefficients of neutral organic chemicals. The first model, referred to here as the 2-parameter partitioning model (PPM), exploits a linear free energy relationship (LFER) of skin permeability coefficient with a linear combination of partition coefficients for octanol–water and air–water systems. The second model is based on the retention time information of nonpolar analytes on comprehensive two-dimensional gas chromatography (GC × GC). The PPM successfully explained variability in the skin permeability data (n = 175) with R² = 0.82 and root mean square error (RMSE) = 0.47 log unit. In comparison, the US-EPA’s model DERMWIN™ exhibited an RMSE of 0.78 log unit. The Zhang model—a 5-parameter LFER equation based on experimental Abraham solute descriptors (ASDs)—performed slightly better with an RMSE value of 0.44 log unit. However, the Zhang model is limited by the scarcity of experimental ASDs. The GC × GC model successfully explained the variance in skin permeability data of nonpolar chemicals (n = 79) with R² = 0.90 and RMSE = 0.23 log unit. The PPM can easily be implemented in US-EPA’s Estimation Program Interface Suite (EPI Suite™). The GC × GC model can be applied to the complex mixtures of nonpolar chemicals.

Leaching of polybrominated diphenyl ethers from microplastics in fish oil: Kinetics and bioaccumulation

Microplastics (MPs) contain high levels of polybrominated diphenyl ethers (PBDEs), which can leach to organism tissues upon ingestion, thereby leading to increased chemical exposure. However, leaching kinetics of inherent contaminants from ingested MPs are poorly studied. The present study characterized the leaching kinetics of PBDEs from sub-millimeter sized MPs in fish oil at relevant body temperatures for marine organisms and assessed exposure risk of MPs for cod fish by a biodynamic model. Diffusion coefficients (D_p) of PBDEs are in the ranges of 1.98 × 10^-19-2.35 × 10^-16 m²·s^-1 in polystyrene, 1.89 × 10^-20-2.07 × 10^-18 m²·s^-1 in acrylonitrile butadiene styrene, and 4.26 × 10^-18-1.72 × 10^-15 m²·s^-1 in polypropylene. A linear function obtained between log D_p of BDE-209 and glass-transition temperature of MPs allows estimation of D_p of BDE-209 contained in other common types of MPs present in the gastrointestinal lipid. The biota-plastic accumulation factors of PBDEs for three plastics were in the range of 4.77 × 10^-14-4.03 × 10^-7. Although bioaccumulation of MPs-affiliated PBDEs is accelerated by oil in the gastrointestinal tract, the modeled steady-state concentrations of PBDEs in cod tissue lipid through ingestion of MPs under the most likely conditions were below the lower end of the global PBDE concentration, implicating that ingestion of MPs by organisms remains a negligible pathway in general.

In vitro methods for predicting the bioconcentration of xenobiotics in aquatic organisms

The accumulation of anthropogenic chemical substances in aquatic organisms is an immensely important issue from the point of view of environmental protection. In the context of the increasing number and variety of compounds that may potentially enter the environment, there is a need for efficient and reliable solutions to assess the risks. However, the classic approach of testing with fish or other animals is not sufficient. Due to very high costs, significant time and labour intensity, as well as ethical concerns, in vivo methods need to be replaced by new laboratory-based tools. So far, many models have been developed to estimate the bioconcentration potential of chemicals. However, most of them are not sufficiently reliable and their predictions are based on limited input data, often obtained with doubtful quality. The octanol-water partition coefficient is still often used as the main laboratory tool for estimating bioconcentration. However, according to current knowledge, this method can lead to very unreliable results, both for neutral species and, above all, for ionic compounds. It is therefore essential to start using new, more advanced and credible solutions on a large scale. Over the last years, many in vitro methods have been newly developed or improved, allowing for a much more adequate estimation of the bioconcentration potential. Therefore, the aim of this work was to review the most recent laboratory methods for assessing the bioconcentration potential and to evaluate their applicability in further research.

Development of liposome/water partition coefficients predictive models for neutral and ionogenic organic chemicals

Membrane/water partition coefficient (K_m/w) is a vital parameter used to characterize the membrane permeability of compounds. Considering the K_m/w value is difficult to observe experimentally for real biological membranes, liposome/water partition coefficient (K_lip/w) is employed to approximate K_m/w. Here, quantitative structure property relationship (QSPR) models for logK_lip/w of the neutral organic chemicals and the neutral form of ionogenic organic chemicals (IOCs) (logK_{lip/w-neutral}), ionic form of IOCs (logK_lip/w-ionic), the speciation-corrected liposome-water distribution ratios at a pH = 7.40 (logD_{lip/w-(pH=7.40)}) were developed. In the modeling, two modeling methods (multiple linear regressions (MLR) and k-nearest neighbor (kNN)) were used. The predictive variables employed here could be calculated from the molecular structure directly. For logK_{lip/w-neutral} and logD_{lip/w-(pH=7.40)}, the logK_OW and logD_OW-based, non-logK_OW and non-logD_OW-based kNN-QSPR and MLR-QSPR models were developed, respectively. The evaluation results implied that the predictive performance of kNN-QSPR models is better than that of MLR-QSPR models. For logK_lip/w-ionic, only one acceptable MLR-QSPR model was developed for cation and anion, respectively. The model quality of the derived models was evaluated following the OECD QSPR models validation guideline. The determination coefficient (R²), leave-one-out cross validation Q² (Q²_LOO) and bootstrapping coefficient (Q²_BOOT), the external validation coefficient (Q²_EXT) of all the models met the acceptable criteria (Q² > 0.600, R² > 0.700); while the root-mean-square error (RMSE) range from 0.351 to 0.857. All the results implied that the models had good goodness-of-fit, robustness and predictive ability. Therefore, the developed models could be used to fill the data gap for substances within the applicability domain on their missing logK_{lip/w-neutral}, logK_lip/w-ionic, logD_{lip/w-(pH=7.40)} values.

Magnetic metal-organic frameworks nanocomposites for negligible-depletion solid-phase extraction of freely dissolved polyaromatic hydrocarbons

The bioavailability of a pollutant is usually evaluated based on its freely dissolved concentration (C_free), which can be measured by negligible-depletion equilibrium extraction that is commonly suffered from long equilibration time. Herein, metal-organic framework (MOF) composites (Fe₃O₄@MIL-101), consists of a magnetic Fe₃O₄ core and a MIL-101 (Cr) MOF shell, is developed as sorbents for negligible-depletion magnetic solid-phase extraction (nd-MSPE) of freely dissolved polyaromatic hydrocarbons (PAHs) in environmental waters. The freely dissolved PAHs in 1000 mL water samples are extracted with 1.5 mg MOF composites, and desorbed with 0.9 mL of acetonitrile under sonication for 5 min. The MOF composites exclude the extraction of dissolved organic matter (DOM) and DOM-associated PAHs by size exclusion. Additionally, the combined interactions (hydrophobic, π-π and π-complexation) between PAHs and composites markedly reduced the extraction equilibration time to < 60 min for all the studied PAHs with logK_OW up to 5.74. Moreover, the porous coordination polymers property of the MOFs makes the proposed nd-MSPE based on the partitioning of PAHs and thus excludes the competitive adsorption of coexisting substances. The developed nd-MSPE approach provides low detection limits (0.08-0.82 ng L^-1), wide linear range (1-1000 ng L^-1) and high precision (relative standard deviations (RSDs) (3.3-4.8%) in determining C_free of PAHs. The measured C_free of PAHs in environmental waters are in good agreement with that of verified method. Given the large diversity in structure and pore size of MOFs, various magnetic MOFs can be fabricated for task-specific nd-MSPE of analytes, presenting a prospective strategy for high-efficiency measuring C_free of contaminants in environments.

Bioavailability of Pyrene Associated with Different Types of Protein Compounds: Direct Evidence for Its Uptake by Daphnia magna

The protein-like dissolved organic matter (DOM) is ubiquitous in aquatic environments. However, the bioavailability of protein-like DOM-associated hydrophobic organic compounds (HOCs) is not well-understood, and in particular, the direct evidence of their uptake by organisms is scarce. In the present work, tryptone (2000 Da), bovine serum albumin (BSA; 66 000 Da), and phycocyanin (120 000 Da) were chosen as model protein-like DOM, which were labeled by commercial fluorescein (cy5) to investigate the uptake mechanisms of protein compound-associated pyrene (a typical HOC) by Daphnia magna. The pyrene concentration in the tissues except the gut and immobilization of D. magna were detected to calculate the bioavailable fraction of protein compound-associated pyrene when the freely dissolved pyrene concentration was controlled through passive dosing devices. The results demonstrated that the tryptone could permeate cellular membrane and directly enter the tissues of  D. magna from the exposure solutions, whereas BSA and phycocyanin might indirectly enter the tissues from the gut. A part of pyrene associated with protein compounds was bioavailable to D. magna; the order of their bioavailable fractions was trypone (54.6-58.1%) > phycocyanin (21.6-32.8%) > BSA (17.7-26.8%). The difference was principally related to the uptake mechanisms of pyrene associated with different types of protein. This work suggests that the protein compound-associated HOCs should be considered to evaluate the bioavailability and eco-environmental hazard of HOCs in natural waters.

Development of chicken and fish muscle protein - Water partition coefficients predictive models for ionogenic and neutral organic chemicals

Muscle protein was one of critical accumulation protein for anthropogenic chemicals. However, few predictive models were constructed for muscle protein up to now. In addition, some ionizable chemicals classes e.g. sulfonates were not successfully modeled in previously models, indicating considerable work would be needed. The major objective of this study was to develop quantitative structure-activity relationship (QSAR) models for predicting the muscle protein-water partition coefficient (logK_MP/w) of chicken and fish. In the modeling, the n-octanol/water distribution coefficient (logD), functional groups, atom-centred fragments and chemical form adjusted descriptors were used to construct the models. The application domain of the derived models was defined by the Euclidean distance-based method and Williams plot. The modeling results indicated that the determination coefficient (R²), leave-one out cross validation Q² (Q²_LOO) and bootstrapping coefficient (Q²_BOOT) of the QSAR models for chicken and fish were 0.882 and 0.929, 0.844 and 0.906, 0.779 and 0.792, respectively, implying the models had good goodness-of-fit and robustness. The coefficient determination (R²_EXT) and external validation coefficient (Q²_EXT) of the validation set for the two models were 0.874 and 0.937, 0.869 and 0.915, respectively, indicating the models had good predictive ability. The predictor variables selected to construct the logK_MP/w models of chicken and fish included logD, the function groups, and the fraction of the ionized species (δ_I). Considering the molecular descriptors used here can be calculated from their molecular structures directly, the developed models could be easily used to fill the logK_MP/w data gap for other chemicals within the applicability domain.

Diagnosis of complex mixture toxicity in sediments: Application of toxicity identification evaluation (TIE) and effect-directed analysis (EDA)

Determining causality of sediment toxicity is of great importance in aquatic risk assessment, but there are tremendous challenges due to joint toxicity of trace pollutants in complex sediment matrices. Two approaches, namely toxicity identification evaluation (TIE) and effect-directed analysis (EDA) have been developed. Conventional sediment TIEs take the advantage of environmental relevance by using whole organism bioassays; however, they suffer from lack of effective methods for specifically identifying major contributors as it typically only evaluates contaminant class rather than specific contaminants. Alternatively, EDA is a powerful tool in identifying causes of sediment toxicity with sophisticated fractionation and chemical analysis of targeted and non-targeted non-polar organic toxicants, but it is not always environmentally relevant due to the use of in-vitro bioassays and exhaustive solvent extraction. An integrated TIE and EDA method would provide an environmentally relevant and toxicant specific approach to effectively determine causality of sediment toxicity by combining the merits of the two methods. Bioavailability-based extraction and dosing techniques are recommended to be incorporated into the integrated method to improve the accuracy of toxicity diagnosis. Besides considering bioavailability in the integrated TIE and EDA approach, the premise of adverse outcome pathways should also be considered. Generally speaking, both TIE and EDA have focused on adverse effects at cellular and organism levels. The addition of trait-based approaches in screening multiple toxicological endpoints helps to extend effects on cellular and organism levels to population level, and provides a better understanding of potential impacts to the community and ecosystem. The outcome pathway underlies the critical role of determining causality in interpreting impacts of complex mixtures to benthic community and aquatic ecosystem.

Application of Fe3O4@MIL-100 (Fe) core-shell magnetic microspheres for evaluating the sorption of organophosphate esters to dissolved organic matter (DOM)

Organophosphate esters (OPEs) are widely used as flame retardants and plasticizers in many products and materials. Because of the potential biologic toxicity on human beings, OPEs are regarded as a class of emerging pollutants. Dissolved organic matters (DOM) have significant effects on the bioavailability and toxicity of the pollutants in the environment. Negligible-depletion solid-phase microextraction (nd-SPME) is an efficient way for measuring the freely dissolved pollutants but suffers from long equilibrium time. Metal-organic frameworks (MOFs) are a class of porous crystalline materials with unique properties such as high pore volume, regular porosity, and tunable pore size, being widely used for the extraction of various organic compounds. Here we developed a novel method for quick determination the sorption coefficients of OPEs to DOM in aquatic phase using Fe₃O₄@MIL-100 (Fe) core-shell magnetic microspheres. The mesoporous structures of the as-synthesized microspheres hindered the extraction of OPEs which associated with humic acid due to the volume exclusion effect. However, the freely dissolved OPEs can access into the mesoporous and then were extracted by MIL-100 (Fe). Due to the small pore size (4.81 nm), large surface area (141 m² g^-1), high pore volume (0.17 g³ g^-1), and ultra-thin MOFs layers, Fe₃O₄@MIL-100 (Fe) core-shell magnetic microspheres have large contact area for the analytes in aqueous phase and therefore the diffusion distance was largely shortened. Besides, the microspheres can be collected conveniently after the extraction process by applying a magnetic field. Compared to the nd-SPME method with 35 h equilibration time (t_90%), the proposed method for these studied OPEs only need 24 min to achieve equilibration. The sorption coefficients (logK_DOC) of the OPEs to humic acid were ranged from 3.84-5.28, which were highly consistent with the results by using polyacrylate-coated fiber and polydimethylsiloxane-coated fiber with nd-SPME.

Improving the accuracy of effect-directed analysis: the role of bioavailability

Aquatic ecosystems have been suffering from contamination by multiple stressors. Traditional chemical-based risk assessment usually fails to explain the toxicity contributions from contaminants that are not regularly monitored or that have an unknown identity. Diagnosing the causes of noted adverse outcomes in the environment is of great importance in ecological risk assessment and in this regard effect-directed analysis (EDA) has been designed to fulfill this purpose. The EDA approach is now increasingly used in aquatic risk assessment owing to its specialty in achieving effect-directed nontarget analysis; however, a lack of environmental relevance makes conventional EDA less favorable. In particular, ignoring the bioavailability in EDA may cause a biased and even erroneous identification of causative toxicants in a mixture. Taking bioavailability into consideration is therefore of great importance to improve the accuracy of EDA diagnosis. The present article reviews the current status and applications of EDA practices that incorporate bioavailability. The use of biological samples is the most obvious way to include bioavailability into EDA applications, but its development is limited due to the small sample size and lack of evidence for metabolizable compounds. Bioavailability/bioaccessibility-based extraction (bioaccessibility-directed and partitioning-based extraction) and passive-dosing techniques are recommended to be used to integrate bioavailability into EDA diagnosis in abiotic samples. Lastly, the future perspectives of expanding and standardizing the use of biological samples and bioavailability-based techniques in EDA are discussed.

Determining equilibrium partition coefficients between lipid/protein and polydimethylsiloxane for highly hydrophobic organic contaminants using preloaded disks

Bioaccumulation of hydrophobic organic contaminants is of great concern and understanding their partitioning to biological phases is crucial for estimating their bioaccumulation potential. The estimation, however, was of large uncertainty for highly hydrophobic organic contaminants (HHOCs) with log K_OW>9 due to the challenge of quantifying their water concentrations. In the present study, partition coefficients between polydimethylsiloxane (PDMS) and storage lipid (K_SL,PDMS), membrane lipid (K_ML,PDMS) and protein (K_pro,PDMS) were measured for 21 polychlorinated biphenyls (PCBs), 14 polybrominated diphenyl ethers (PBDEs), dechlorane plus (DP) and decabromodiphenyl ethane (DBDPE), covering log K_OW from 5.07 to 11.6, using a preloaded PDMS depletion method. The values of K_SL,PDMS, K_ML,PDMS and K_pro,PDMS were in the ranges of 5.36-52.5, 0.286-11.8 and 0.067-2.62g/g, respectively, being relatively constant although their K_OW values extend more than six orders of magnitude. The relative sorption capacity of the biological phases showed storage lipid was the dominant sorption phase in biota, followed by membrane lipid and protein was the lowest. The K_PDMS,pro values of the compounds with log K_OW<9 were similar (0.382-14.9g/g) regardless of the thickness of preloaded PDMS disks (58-209μm). For HHOCs, however, K_PDMS,pro values dropped when thinner PDMS disks were used, as a result of slow diffusion of HHOCs in PDMS. The K_PDMS,pro values of HHOCs measured by 58-μm PDMS disks ranged from 1.78 to 6.85g/g, which was consistent with compounds with log K_OW<9. This validated that partition coefficients between PDMS and biological phases were independent of chemical hydrophobicity, showing the advantage of using PDMS-based methods to directly estimate bioaccumulation potential of HHOCs.

Development and evaluation of predictive model for bovine serum albumin-water partition coefficients of neutral organic chemicals

The binding of organic chemicals to serum albumin can significantly reduce their unbound concentration in blood and affect their biological reactions. In this study, we developed a new QSAR model for bovine serum albumin (BSA) - water partition coefficients (K_BSA/W) of neutral organic chemicals with large structural variance, logK_BSA/W values covering 3.5 orders of magnitude (1.19-4.76). All chemical geometries were optimized by semi-empirical PM6 algorithm. Several quantum chemical parameters that reflect various intermolecular interactions as well as hydrophobicity were selected to develop QSAR model. The result indicates the regression model derived from logK_ow, the most positive net atomic charges on an atom, Connolly solvent excluded volume, polarizability, and Abraham acidity could explain the partitioning mechanism of organic chemicals between BSA and water. The simulated external validation and cross validation verifies the developed model has good statistical robustness and predictive ability, thus can be used to estimate the logK_BSA/W values for chemicals in application domain, accordingly to provide basic data for the toxicity assessment of the chemicals.

Predicting Partitioning and Diffusion Properties of Nonpolar Chemicals in Biotic Media and Passive Sampler Phases by GC × GC

The chemical parameters needed to explain and predict bioavailability, biodynamics, and baseline toxicity are not readily available for most nonpolar chemicals detected in the environment. Here, we demonstrate that comprehensive two-dimensional gas chromatography (GC × GC) retention times can be used to predict 26 relevant properties for nonpolar chemicals, specifically: partition coefficients for diverse biotic media and passive sampler phases; aquatic baseline toxicity; and relevant diffusion coefficients. The considered biotic and passive sampler phases include membrane and storage lipids, serum and muscle proteins, carbohydrates, algae, mussels, polydimethylsiloxane, polyethylene, polyoxymethylene, polyacrylate, polyurethane, and semipermeable membrane devices. GC × GC-based chemical property predictions are validated with a compilation of 1038 experimental property data collected from the literature. As an example application, we overlay a map of baseline toxicity to fathead minnows onto the separated analyte signal of a polychlorinated alkanes (chlorinated paraffins) technical mixture that contains 7820 congeners. In a second application, GC × GC-estimated properties are used to parametrize multiphase partitioning models for mammalian tissues and organs. In a third example, we estimate chemical depuration kinetics for mussels. Finally, we illustrate an approach to screen the GC × GC chromatogram for nonpolar chemicals of potentially high concern, defined based on their GC × GC-estimated biopartitioning properties, diffusion properties, and baseline toxicity.

Solubility enhancement of dioxins and PCBs by surfactant monomers and micelles quantified with polymer depletion techniques

Partitioning of super-hydrophobic organic contaminants (SHOCs) to dissolved or colloidal materials such as surfactants can alter their behaviour by enhancing apparent aqueous solubility. Relevant partition constants are, however, challenging to quantify with reasonable accuracy. Partition constants to colloidal surfactants can be measured by introducing a polymer (PDMS) as third phase with known PDMS-water partition constant in combination with the mass balance approach. We quantified partition constants of PCBs and PCDDs (log KOW 5.8-8.3) between water and sodium dodecyl sulphate monomers (KMO) and micelles (KMI). A refined, recently introduced swelling-based polymer loading technique allowed highly precise (4.5-10% RSD) and fast (<24 h) loading of SHOCs into PDMS, and due to the miniaturisation of batch systems equilibrium was reached in <5 days for KMI and <3 weeks for KMO. SHOC losses to experimental surfaces were substantial (8-26%) in monomer solutions, but had a low impact on KMO (0.10-0.16 log units). Log KMO for PCDDs (4.0-5.2) were approximately 2.6 log units lower than respective log KMI, which ranged from 5.2 to 7.0 for PCDDs and 6.6-7.5 for PCBs. The linear relationship between log KMI and log KOW was consistent with more polar and moderately hydrophobic compounds. Apparent solubility increased with increasing hydrophobicity and was highest in micelle solutions. However, this solubility enhancement was also considerable in monomer solutions, up to 200 times for OCDD. Given the pervasive presence of surfactant monomers in typical field scenarios, these data suggest that low surfactant concentrations may be effective long-term facilitators for subsurface transport of SHOCs.

Fate of polychlorinated biphenyls in a contaminated lake ecosystem: combining equilibrium passive sampling of sediment and water with total concentration measurements of biota

Equilibrium sampling devices can be applied to study and monitor the exposure and fate of hydrophobic organic chemicals on a thermodynamic basis. They can be used to determine freely dissolved concentrations and chemical activity ratios and to predict equilibrium partitioning concentrations of hydrophobic organic chemicals in biota lipids. The authors' aim was to assess the equilibrium status of polychlorinated biphenyls (PCBs) in a contaminated lake ecosystem and along its discharge course using equilibrium sampling devices for measurements in sediment and water and by also analyzing biota. The authors used equilibrium sampling devices (silicone rubber and polyethylene [PE]) to determine freely dissolved concentrations and chemical activities of PCBs in the water column and sediment porewater and calculated for both phases the corresponding equilibrium concentrations and chemical activities in model lipids. Overall, the studied ecosystem appeared to be in disequilibrium for the studied phases: sediment, water, and biota. Chemical activities of PCBs were higher in sediment than in water, which implies that the sediment functioned as a partitioning source of PCBs and that net diffusion occurred from the sediment to the water column. Measured lipid-normalized PCB concentrations in biota were generally below equilibrium lipid concentrations relative to the sediment (CLip ⇌Sed ) or water (CLip ⇌W ), indicating that PCB levels in the organisms were below the maximum partitioning levels. The present study shows the application versatility of equilibrium sampling devices in the field and facilitates a thermodynamic understanding of exposure and fate of PCBs in a contaminated lake and its discharge course.

Biotransformation of dichlorodiphenyltrichloroethane in the benthic polychaete, Nereis succinea: quantitative estimation by analyzing the partitioning of chemicals between gut fluid and lipid

Biotransformation plays an important role in the bioaccumulation and toxicity of a chemical in biota. Dichlorodiphenyltrichloroethane (DDT) commonly co-occurs with its metabolites (dichlorodiphenyldichloroethane [DDD] and dichlorodiphenyldichloroethylene [DDE]), in the environment; thus it is a challenge to accurately quantify the biotransformation rates of DDT and distinguish the sources of the accumulated metabolites in an organism. The present study describes a method developed to quantitatively analyze the biotransformation of p,p'-DDT in the benthic polychaete, Nereis succinea. The lugworms were exposed to sediments spiked with DDT at various concentrations for 28 d. Degradation of DDT to DDD and DDE occurred in sediments during the aging period, and approximately two-thirds of the DDT remained in the sediment. To calculate the biotransformation rates, residues of individual compounds measured in the bioaccumulation testing (after biotransformation) were compared with residues predicted by analyzing the partitioning of the parent and metabolite compounds between gut fluid and tissue lipid (before biotransformation). The results suggest that sediment ingestion rates decreased when DDT concentrations in sediment increased. Extensive biotransformation of DDT occurred in N. succinea, with 86% of DDT being metabolized to DDD and <2% being transformed to DDE. Of the DDD that accumulated in the lugworms, approximately 70% was the result of DDT biotransformation, and the remaining 30% was from direct uptake of sediment-associated DDD. In addition, the biotransformation was not dependent on bulk sediment concentrations, but rather on bioaccessible concentrations of the chemicals in sediment, which were quantified by gut fluid extraction. The newly established method improved the accuracy of prediction of the bioaccumulation and toxicity of DDTs.

Sorptive capacity of membrane lipids, storage lipids, and proteins: a preliminary study of partitioning of organochlorines in lean fish from a PCB-contaminated freshwater lake

Knowledge on the internal distribution of halogenated organic chemicals (HOCs) would improve our understanding of dose-effect relationships and subsequently improve risk assessment of contaminated sites. Herein, we determine the concentrations of HOCs based on equilibrium partitioning in storage lipids, membrane lipids, and proteins in field-contaminated fish using equilibrium sampling devices. The study shows the importance of protein as a sorptive phase in lean fish. Our results provide a basis for using species-specific equilibrium partitioning coefficients between sorptive tissues and fish internal water as a substitute for K(ow) in, for example, upgrading models that simulate food-chain accumulation of the chemical.

Bioaccumulation of highly hydrophobic organohalogen flame retardants from sediments: application of toxicokinetics and passive sampling techniques

Highly hydrophobic organohalogen flame retardants (HHOFRs) are found ubiquitously in the environment; therefore, a better understanding of their bioavailability is needed. In the current study, bioaccumulation testing using the oligochaete, Lumbriculus variegatus, and passive sampling (solid-phase microextraction (SPME)) were performed to study the bioaccumulation potential of HHOFRs, including decabromodiphenyl ether (deca-BDE), decabromodiphenyl ethane (DBDPE), and dechlorane plus (DP), in laboratory-spiked and field-collected sediments. The HHOFRs were bioavailable to L. variegatus even though their biota-sediment accumulation factors were low (0.016 ± 0.002 to 0.48 ± 0.082 g organic carbon/g lipid, syn-DP > anti-DP > deca-BDE > DBDPE). Hydrophobicity and stereoisomerism affected HHOFR bioavailability. Meanwhile, HHOFR concentrations on the SPME fibers (Cf) correlated with those in biota (Cb), suggesting the potential application of SPME in bioavailability prediction for those compounds. The log Cf to log Cb correlation for deca-BDE and DP had a greater intercept than that for polychlorinated biphenyls (data obtained from the literature) although the slopes were similar, while data for DBDPE fell on the regression line for PCBs, implying some uncertainty in application of SPMEs across chemical classes. The increasing sorptive ability of proteins for HHOFRs in comparison to the less-brominated BDEs suggested that protein-binding should be considered when estimating bioaccumulation potential of HHOFRs in benthic invertebrates.

Predicting partition coefficients of Polyfluorinated and organosilicon compounds using polyparameter linear free energy relationships (PP-LFERs)

The environmental behavior, fate, and effects of polyfluorinated compounds (PFCs) and organosilicon compounds (OSCs) have received increasing attention in recent years. In this study, polyparameter linear free energy relationships (PP-LFERs) were evaluated for predicting partition coefficients of neutral PFCs and OSCs, using experimental data for fluorotelomer alcohols (FTOHs) and cyclic volatile methylsiloxanes (cVMS) reported in the literature and measured newly for this work. It was found that the recently proposed PP-LFER model that uses the McGowan characteristic volume (V), the logarithmic hexadecane-air partition coefficient (L), and three polar interaction descriptors can accurately describe partition coefficients of PFCs and OSCs. The prediction errors were <1 log unit when literature descriptors were used, and the errors were reduced to <0.2 log units on average by further optimization of the descriptors. Surprisingly, the conventional forms of PP-LFERs that include the excess molar refraction (E) sometimes led to substantial errors (>1 log unit) even with optimized parameters. The system parameters for octanol-water, air-water, octanol-air, oil-water, liposome-water, and organic carbon-water partition coefficients as well as the solute descriptors for FTOHs and cVMS were recalibrated in this work, which should provide even more reliable predictions of partition coefficients. The results also confirm the consistency of the published experimental partition coefficients for FTOHs and cVMS.

Desorption kinetics of hydrophobic organic contaminants from marine plastic pellets

This study investigated the desorption behavior of polychlorinated biphenyls (PCBs) from marine plastic pellets. Long-term desorption experiments were conducted using field-collected polyethylene (PE) pellets. The results indicate that the desorption kinetics highly depends on the PE-water partition coefficients of PCB congeners. After 128 d of the experiment, the smallest congener considered (CB 8) had desorbed nearly completely (98%), whereas major fractions (90-99%) of highly chlorinated congeners remained in the pellets. An intraparticle diffusion model mostly failed to reproduce the desorption kinetics, whereas an aqueous boundary layer (ABL) diffusion model well approximated the data. The desorption half-lives are estimated to 14 d to 210 years for CB 8 to CB 209 in an actively stirred solution (ABL thickness: 30 μm). Addition of methanol to water enhanced the desorption to a large extent. A need for further work to explore roles of organic matter in facilitating solute transfer is suggested.

Applicability of passive sampling to bioanalytical screening of bioaccumulative chemicals in marine wildlife

Quantification of bioaccumulative contaminants in biota is time and cost-intensive and the required extensive cleanup steps make it selective toward targeted chemical groups. Therefore tissue extracts prepared for chemical analysis are not amenable to assess the combined effects of unresolved complex mixtures. Passive equilibrium sampling with polydimethylsiloxane (PDMS) has the potential for unbiased sampling of mixtures, and the PDMS extracts can be directly dosed into cell-based bioassays. The passive sampling approach was tested by exposing PDMS to lipid-rich tissue (dugong blubber; 85% lipid) spiked with a known mixture of hydrophobic contaminants (five congeners of tetra- to octachloro-dibenzo-p-dioxins). The equilibrium was attained within 24 h. Lipid-PDMS partition coefficients (Klip-PDMS) ranged from 20 to 38, were independent of hydrophobicity, and within the range of those previously measured for organochlorine compounds. To test if passive sampling can be combined with bioanalysis without the need for chemical cleanup, spiked blubber-PDMS extracts were dosed into the CAFLUX bioassay, which specifically targets dioxin-like chemicals. Small quantities of lipids coextracted by the PDMS were found to affect the kinetics in the regularly applied 24-h bioassay; however, this effect was eliminated by a longer exposure period (72 h). The validated method was applied to 11 unspiked dugong blubber samples with known (native) dioxin concentrations. These results provide the first proof of concept for linking passive sampling of lipid-rich tissue with cell-based bioassays, and could be further extended to other lipid rich species and a wider range of bioanalytical end points.

Determination of polyparameter linear free energy relationship (pp-LFER) substance descriptors for established and alternative flame retardants

Polyparameter linear free energy relationships (pp-LFERs) can predict partition coefficients for a multitude of environmental and biological phases with high accuracy. In this work, the pp-LFER substance descriptors of 40 established and alternative flame retardants (e.g., polybrominated diphenyl ethers, hexabromocyclododecane, bromobenzenes, trialkyl phosphates) were determined experimentally. In total, 251 data for gas-chromatographic (GC) retention times and liquid/liquid partition coefficients (K) were measured and used to calibrate the pp-LFER substance descriptors. Substance descriptors were validated through a comparison between predicted and experimental log K for the systems octanol/water (K(ow)), water/air (K(wa)), organic carbon/water (K(oc)) and liposome/water (K(lipw)), revealing a high reliability of pp-LFER predictions based on our descriptors. For instance, the difference between predicted and experimental log K(ow) was <0.3 log units for 17 out of 21 compounds for which experimental values were available. Moreover, we found an indication that the H-bond acceptor value (B) depends on the solvent for some compounds. Thus, for predicting environmentally relevant partition coefficients it is important to determine B values using measurements in aqueous systems. The pp-LFER descriptors calibrated in this study can be used to predict partition coefficients for which experimental data are unavailable, and the predicted values can serve as references for further experimental measurements.