Partitioning of polychlorinated biphenyls into human cells and adipose tissues: evaluation of octanol, triolein, and liposomes as surrogates

Polyparameter Linear Free Energy Relationships for Partitioning of Neutral Organic Compounds to Storage Lipids

Storage lipids are an important compartment in the bioaccumulation of neutral organic compounds. Reliable models for predicting storage lipid-water and storage lipid-air partition coefficients (Kislip/w and Kislip/a), as well as their temperature dependence, are considered useful. Polyparameter linear free energy relationships (PP-LFERs) are accurate, general, and mechanistically clear models for predicting partitioning-related physicochemical quantities. About a decade ago, PP-LFERs were calibrated for Kislip/w at the physiological temperature of 37 °C. However, to date, a comprehensive collection and sufficiently reliable PP-LFERs for Kislip/w and Kislip/a at the most common standard temperature of 25 °C are still lacking. In this study, experimentally based Kislip/w and/or Kislip/a values at 25 °C for 278 compounds were extensively collected or converted from the literature. Subsequently, PP-LFERs were calibrated for Kislip/w and Kislip/a at 25 °C, performing well over 10 orders of magnitude with root-mean-square errors of 0.17-0.21 log units for compounds with reliable descriptors. Furthermore, standard internal energy changes of transfer from water or air to storage lipids for 158 compounds were derived and used to calibrate PP-LFERs for estimating the temperature dependence of Kislip/w or Kislip/a. Additionally, using PP-LFERs, low-density polyethylene was confirmed to be a better storage lipid analogue than silicone and polyoxymethylene in the equilibrium passive sampling of nonpolar and H-bond acceptor polar compounds.

Miniaturized method for the quantification of persistent organic pollutants and their metabolites in HepG2 cells: assessment of their biotransformation

Biotransformation can greatly influence the accumulation and, subsequently, toxicity of substances in living beings. Although traditionally these studies to quantify metabolization of a compound have been carried out with in vivo species, currently, in vitro test methods with very different cell lines are being developed for their evaluation. However, this is still a very limited field due to multiple variables of a very diverse nature. So, an increasing number of analytical chemists are working with cells or other similar biological samples of very small size. This makes it necessary to address the development of analytical methods that allow determining their concentration both inside the cells and in their exposure medium. The aim of this study is to develop a set of analytical methodologies for the quantification of polycyclic aromatic hydrocarbons, PAHs (phenanthrene, PHE), and polybrominated diphenyl ethers, PBDEs (2,2',4,4'-tetrabromodiphenyl ether, BDE-47), and their major metabolites in cells and their exposure medium. Analytical methodologies, based on miniaturized ultrasound probe-assisted extraction, gas chromatography-mass spectrometry-microelectron capture detector (GC-MS-µECD), and liquid chromatography-fluorescence detector (LC-FL) determination techniques, have been optimized and then applied to a biotransformation study in HepG2 at 48 h of exposure. Significant concentrations of the major metabolites of PHE (1-OH, 2-OH, 3-OH, 4-OH-, and 9-OH-PHE) and BDE-47 (5-MeO-, 5-OH-, and 3-OH-BDE-47) were detected and quantified inside the cells and in the exposure medium. These results provide a new method for determination and improve information on the metabolization ratios for a better knowledge of the metabolic pathways and their toxicity.

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.

Adsorption of neutral organic compounds on polar and nonpolar microplastics: Prediction and insight into mechanisms based on pp-LFERs

Adsorption of 18 neutral organic compounds (OCs) on polar (polybutylene succinate (PBS) and polycaprolactone (PCL)) and nonpolar (low-density polyethylene (LDPE) and polystyrene (PS)) microplastics (MPs) were investigated. The adsorption coefficients (K_d) varied with ranges of 130-42,002, 124-27,768, 6.40-10,713, and 1.52-10,332 L kg^-1 for adsorption on PCL, PBS, LDPE, and PS MPs, respectively. The polar MPs showed greater adsorption capacities than nonpolar MPs. Non-specific interaction, i.e. hydrophobic partition played a crucial role in the adsorption of OCs on all MPs, while polar interactions also contributed significantly to the greater adsorption on polar MPs. Poly-parameter linear free energy relationships (pp-LFERs) with multiple linear regression (MLR) and feedforward network (FN) were then employed to model the adsorption of OCs on MPs so as to obtain deep insights into adsorption mechanisms. The MLR models achieved R_adj² of 0.90-0.97 and root mean square error (RMSE) of 0.13-0.38 log units, while the FN models achieved R_adj² of 0.85-0.90 and RMSE of 0.21-0.60 log units. The MLR models are more accurate under selected equilibrium concentrations while FN models are capable of making predictions under varying equilibrium concentrations. Lastly, both MLR and FN models showed good prediction on literature adsorption data on nonpolar MPs.

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.

Prediction of Partition Coefficients of Environmental Toxins Using Computational Chemistry Methods

The partitioning of compounds between aqueous and other phases is important for predicting toxicity. Although thousands of octanol-water partition coefficients have been measured, these represent only a small fraction of the anthropogenic compounds present in the environment. The octanol phase is often taken to be a mimic of the inner parts of phospholipid membranes. However, the core of such membranes is typically more hydrophobic than octanol, and other partition coefficients with other compounds may give complementary information. Although a number of (cheap) empirical methods exist to compute octanol-water (log k _OW) and hexadecane-water (log k _HW) partition coefficients, it would be interesting to know whether physics-based models can predict these crucial values more accurately. Here, we have computed log k _OW and log k _HW for 133 compounds from seven different pollutant categories as well as a control group using the solvation model based on electronic density (SMD) protocol based on Hartree-Fock (HF) or density functional theory (DFT) and the COSMO-RS method. For comparison, XlogP3 (log k _OW) values were retrieved from the PubChem database, and KowWin log k _OW values were determined as well. For 24 of these compounds, log k _OW was computed using potential of mean force (PMF) calculations based on classical molecular dynamics simulations. A comparison of the accuracy of the methods shows that COSMO-RS, KowWin, and XlogP3 all have a root-mean-square deviation (rmsd) from the experimental data of ≈0.4 log units, whereas the SMD protocol has an rmsd of 1.0 log units using HF and 0.9 using DFT. PMF calculations yield the poorest accuracy (rmsd = 1.1 log units). Thirty-six out of 133 calculations are for compounds without known log k _OW, and for these, we provide what we consider a robust prediction, in the sense that there are few outliers, by averaging over the methods. The results supplied may be instrumental when developing new methods in computational ecotoxicity. The log k _HW values are found to be strongly correlated to log k _OW for most compounds.

SSP silicone-, lipid- and SPMD-water partition coefficients of seventy hydrophobic organic contaminants and evaluation of the water concentration calculator for SPMD

Passive sampling is increasingly applied for monitoring neutral hydrophobic compounds (HOC) in various environmental media like water, sediment, air and also soft biota tissue. Passive samplers for HOC are often constructed from permeable polymers like silicone and polyethylene (PE), while also SPMD are often applied. Their HOC uptake can be converted to freely dissolved or equivalent lipid-based concentrations using appropriate partition coefficients with or without the use of kinetic uptake models to adjust for non-equilibrium. To facilitate such conversions for seventy HOC partition coefficients are derived by combining polymer-water for Altesil™ silicone and PE, with new and earlier published polymer-polymer, polymer-lipid partition coefficients. Derived SSP silicone-water, lipid-water (K_lip/w), and SPMD-water (K_spmd/w) partition coefficients demonstrate good agreement with literature data, except for K_spmd/w. For SPMD, this work demonstrates a linear K_spmd/w - K_ow relationship (R² = 0.99) in contrast to the parabolic K_spmd/w - K_ow relationship utilized in the USGS "SPMD Water Concentrations Calculator". Following a thorough evaluation of this Calculator it is recommended that in combination with revised K_spmd/w, a radical different model approach should be used for obtaining accurate water concentrations from passive sampling with SPMD.

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.

Occurrence of artificial sweeteners in human liver and paired blood and urine samples from adults in Tianjin, China and their implications for human exposure

In this study, acesulfame (ACE), saccharin (SAC) and cyclamate (CYC) were found in all paired urine and blood samples collected from healthy adults, with mean values of 4070, 918 and 628 ng mL(-1), respectively, in urine and 9.03, 20.4 and 0.72 ng mL(-1), respectively, in blood. SAC (mean: 84.4 ng g(-1)) and CYC (4.29 ng g(-1)) were detectable in all liver samples collected from liver cancer patients, while ACE was less frequently detected. Aspartame (ASP) was not found in any analyzed human sample, which can be explained by the fact that this chemical metabolized rapidly in the human body. Among all adults, significantly positive correlations between SAC and CYC levels were observed (p < 0.001), regardless of human matrices. Nevertheless, no significant correlations between concentrations of SAC (or CYC) and ACE were found in any of the human matrices. Our results suggest that human exposure to SAC and CYC is related, whereas ACE originates from a discrete source. Females (or young adults) were exposed to higher levels of SAC and CYC than males (or elderly). The mean renal clearance of SAC was 730 mL per day per kg in adults, which was significantly (p < 0.001) lower than those for CYC (10 800 mL per day per kg) and ACE (10 300 mL per day per kg). The average total daily intake of SAC and ACE was 9.27 and 33.8 μg per kg bw per day, respectively.

Equations for water-triolein partition coefficients for neutral species; comparison with other water-solvent partitions, and environmental and toxicological processes

Linear free energy relationships, LFERs, have been constructed for water-triolein partition coefficients for neutral species. It is shown that separate equations are required for wet and dry triolein. From a comparison of the equation coefficients for water-wet triolein with those for 52 other water-solvent systems it is shown that there is little correspondence between triolein and any of the 52 other solvents - only the water-isopropyl myristate system is close to the water-wet triolen system. A comparison of equation coefficients for the water-wet triolein system with LFER coefficients of 16 environmentally important processes shows that wet triolein is not a suitable model for any of the processes, although a number of other water-solvent systems are possible models for some of the environmental processes. A comparison of LFER coefficients with those of 17 aqueous toxicological processes reveals that most of the water-solvent systems, including water-wet triolein, will be poor models for any of the toxicological systems, but the water-lower alcohol systems show promise as models for a number of the toxicological systems. Our method of comparison of coefficients for LFERs that have exactly the same independent variables can be extended to various other types of system.

Applications of polyparameter linear free energy relationships in environmental chemistry

Partitioning behavior of organic chemicals has tremendous influences on their environmental distribution, reaction rates, bioaccumulation, and toxic effects. Polyparameter linear free energy relationships (PP-LFERs) have been proven to be useful to characterize the equilibrium partitioning of organic chemicals in various environmental and technical partitioning systems and predict the respective partition coefficients. Over the past decade, PP-LFER solute descriptors for numerous environmentally relevant organic chemicals and system parameters for environmentally important partitioning systems have been determined, extending substantially the applicability of the PP-LFER approaches. However, the information needed for the use of PP-LFERs including descriptors and parameters is scattered over a large number of publications. In this work, we review the state of the art of the PP-LFER approaches in environmental chemical applications. The solute descriptors and system parameters reported in the literature and the availability of their database are summarized, and their calibration and prediction methods are overviewed. We also describe tips and pitfalls associated with the use of the PP-LFER approaches and identify research needs to improve further the usefulness of PP-LFERs for environmental chemistry.