Prediction of Unbound Fractions for in Vitro-in Vivo Extrapolation of Biotransformation Data.

For in vitro-in vivo extrapolation of biotransformation data, the different sorptive environments in vitro and in vivo need to be considered. The most common approach for doing so is using the ratio of unbound fractions in vitro and in vivo. In the literature, several algorithms for prediction of these unbound fractions are available. In this study, we present a theoretical evaluation of the most commonly used algorithms for prediction of unbound fractions in S9 assays and blood and compare prediction results with empirical values from the literature. The results of this analysis prove a good performance of "composition-based" algorithms, i.e. algorithms that represent the inhomogeneous composition of in vitro assay and in vivo system and describe sorption to the individual components (lipids, proteins, water) in the same way. For strongly sorbing chemicals, these algorithms yield constant values for the ratio of unbound fractions in vitro and in vivo. This is mechanistically plausible, because in these cases, the chemicals are mostly bound, and the ratio of unbound fractions is determined by the volume ratio of sorbing components in both phases.

Environmental Sorption Behavior of Ionic and Ionizable Organic Chemicals.

Traditionally our tools for environmental risk assessment of organic chemicals have been developed for neutral chemicals. However, many commercial chemicals are ionic or ionizable and require different tools and approaches for their assessment. In recent years this task starts to obtain increasing attention but our understanding for their environmental fate is still far behind that for neutral chemicals. This review first gives an overview on the principles that govern ionic partitioning in environmental systems which are more complex than the simple partition processes of neutral chemicals. Second, a summary of our current knowledge on various topics such as bioaccumulation, sorption in soils, and nonspecific-toxicity reveals that ionic species can actually be quite hydrophobic contrary to commonly held beliefs. Eventually, we discuss existing models for the quantitative prediction of organic ions' sorption in soils and biota. We have to assert that the available model tools are quite restricted in their application range compared to neutral chemicals which is due to the higher complexity of the various ionic sorption processes. In order to further advance our understanding more high-quality sorption data are needed with a focus on multivalent and zwitterionic ions in all partition systems as well as cations in biological matrices.

Predicting Uncoupling Toxicity of Organic Acids Based on Their Molecular Structure Using a Biophysical Model.

We present a purely mechanistic model to predict protonophoric uncoupling activity ECw of organic acids. All required input information can be derived from their chemical structure. This makes it a convenient predictive model to gain valuable information on the toxicity of organic chemicals already at an early stage of development of new commercial chemicals (e.g., in agriculture or pharmaceutical industries). A critical component of the model is the consideration of the possible formation of heterodimers from the neutral and anionic monomer, and its permeation through the membrane. The model was tested against literature data measured in chromatophores, submitochondrial particles, isolated mitochondria, and intact green algae cells with good success. It was also possible to reproduce pH-dependencies in isolated mitochondria and intact cells. Besides the prediction of the ECw, the mechanistic nature of the model allows researchers to draw direct conclusions on the impact of single input factors such as pH- and voltage-gradients across the membrane, the anionic and neutral membrane permeability, and the heterodimerization constant. These insights are of importance in drug design or chemical regulation.

Membrane/Water Partitioning and Permeabilities of Perfluoroalkyl Acids and Four of their Alternatives and the Effects on Toxicokinetic Behavior.

The search for alternatives to bioaccumulative perfluoroalkyl acids (PFAAs) is ongoing. New, still highly fluorinated alternatives are produced in hopes of reducing bioaccumulation. To better estimate this bioaccumulative behavior, we performed dialysis experiments and determined membrane/water partition coefficients, Kmem/w, of six perfluoroalkyl carboxylic acids (PFCAs), three perfluoroalkanesulfonic acids, and four alternatives. We also investigated how passive permeation might influence the uptake kinetics into cells, measuring the passive anionic membrane permeability Pion through planar lipid bilayers for six PFAAs and three alternatives. Experimental Kmem/w and Pion were both predicted well by the COSMO-RS theory (log RMSE 0.61 and 0.46, respectively). Kmem/w values were consistent with the literature data, and alternatives showed similar sorption behavior as PFAAs. Experimental Pion values were high enough to explain observed cellular uptake by passive diffusion with no need to postulate the existence of active uptake processes. However, predicted pKa and neutral permeabilities suggest that also the permeation of the neutral species should be significant in case of PFCAs. This can have direct consequences on the steady-state distribution of PFAAs across cell membranes and thus toxicity. Consequently, we propose a model to predict pH-dependent baseline toxicity based on Kmem/w, which considers the permeation of both neutral and anionic species.

Yolk Sac of Zebrafish Embryos as Backpack for Chemicals?

The zebrafish embryo (Danio rerio) has developed into one of the most important nonsentient animal models for the hazard assessments of chemicals, but the processes governing its toxicokinetics (TK) are poorly understood. This study compares the uptake of seven test compounds into the embryonic body and the yolk sac of the zebrafish embryo using TK experiments, a dialysis approach, thermodynamic calculations, and kinetic modeling. Experimental data show that between 95% (4-iodophenol) and 67% (carbamazepine) of the total internal amount in 26 h post fertilization (hpf) embryos and between 80 and 49% in 74 hpf embryos were found in the yolk. Thus, internal concentrations determined for the whole embryo overestimate the internal concentration in the embryonic body: for the compounds of this study, up to a factor of 5. Partition coefficients for the embryonic body and a one-compartment model with diffusive exchange were calculated for the neutral test compounds and agreed reasonably with the experimental data. For prevalently ionic test compounds at exposure pH (bromoxynil, paroxetine), however, the extent and the speed of uptake were low and could not be modeled adequately. A better understanding of the TK of ionizable test compounds is essential to allow assessment of the validity of this organismic test system for ionic test compounds.

Predicting the Gas/Particle Distribution of SVOCs in the Indoor Environment Using Poly Parameter Linear Free Energy Relationships.

Understanding the partitioning of semi volatile organic compounds (SVOCs) between gas phase and particle phase is essential for exposure analysis and risk assessment in the indoor environment. Numerous attempts have been made to calculate gas/particle partitioning coefficients Kip. Single-parameter adsorption and absorption models, which relate Kip to the vapor pressure Ps or the octanol/air distribution coefficient KOA are usually applied. In this work we use poly parameter Linear Free Energy Relationships (pp-LFER) to describe the partitioning behavior of 14 SVOCs with high relevance for the indoor environment. The pp-LFER concept is based on Abraham descriptors and considers interactions between molecule and particle by separate parameters. van der Waals interactions can be approximated by the logarithm of the hexadecane/air partitioning coefficient (log KHdA = L), which is a key parameter for the 14 polar but nonionizable organic esters being studied here. For many of the examined compounds experimentally determined L-values were not available and had to be measured using gas chromatography. It is shown that the pp-LFER method is a strong alternative to single-parameter approaches and gives reliable coefficients for gas/particle distribution in the indoor environment.

Passive Permeability of Planar Lipid Bilayers to Organic Anions.

The membrane permeability P of organic ions was reported to be governed by the structure of the permeating molecule. Thus far, it is unclear whether the ion structure alters membrane partition or translocation proper across the membrane. Here, we obtained P values for 24 anionic compounds (18 concrete values, 6 upper limits) measuring the current that they carry through folded planar lipid bilayers. The P values range over more than 10 log units. Our measured permeability values correlate well (r = 0.95; logRMSE 0.74) with the hexadecane/water partition coefficients of the respective chemicals predicted by the COSMO-RS theory. Other attempts to predict P from the partition coefficient of the neutral molecule and from the solvation energy (Born energy) that opposes transfer into the membrane once the molecule is charged were unsuccessful. The uncertainties in assigning an effective radius to nonspherical molecules were much too large. The observation underlines that the actual structure of the molecules needs to be considered to predict partition and thus P by the solubility-diffusion model.

In Vitro- in Vivo Extrapolation of Hepatic Metabolism for Different Scenarios - a Toolbox.

The extrapolation of metabolism data from in vitro experiments to in vivo clearances can provide useful information in the fields of pharmacokinetics and toxicokinetics. Depending on the purpose, different toxicokinetic models are used, and these different models require the in vivo metabolic information in different forms. In this study, a comprehensive toolbox for in vitro- in vivo extrapolation (IVIVE) of hepatic metabolism is presented addressing a variety of different extrapolation goals: extrapolation to hepatic blood clearance, extrapolation to organ clearance, extrapolation to whole-body clearance, and extrapolation to clearance at the level of hepatocytes. The use of the extrapolated clearances for calculation of extraction efficiencies and the use in physiologically based pharmacokinetic models are discussed. Furthermore, a sensitivity analysis demonstrates which parameters affect the accuracy of the extrapolation results the most, and the presented extrapolation procedure is evaluated by comparison to experimental data from perfused liver experiments.

Application of Experimental Polystyrene Partition Constants and Diffusion Coefficients to Predict the Sorption of Neutral Organic Chemicals to Multiwell Plates in in Vivo and in Vitro Bioassays.

Sorption to the polystyrene (PS) of multiwell plates can affect the exposure to organic chemicals over time in in vitro and in vivo bioassays. Experimentally determined diffusion coefficients in PS ( DPS) were in a narrow range of 1.25 to 8.0 · 10-16 m2 s-1 and PS-water partition constants ( KPS/w) ranged from 0.04 to 5.10 log-units for 22 neutral organic chemicals. A kinetic model, which explicitly accounts for diffusion in the plastic, was applied to predict the depletion of neutral organic chemicals from different bioassay media by sorption to various multiwell plate formats. For chemicals with log Kow > 3, the medium concentrations decreased rapidly and considerably in the fish embryo toxicity assay but medium concentrations remained relatively constant in the cell-based bioassays with medium containing 10% fetal bovine serum (FBS), emphasizing the ability of the protein- and lipid-rich medium to compensate for losses by multiwell plate sorption. The PS sorption data may serve not only for exposure assessment in bioassays but also to model the contaminant uptake by and release from plastic packaging material and the chemical transport by PS particles in the environment.

The impact of desorption kinetics from albumin on hepatic extraction efficiency and hepatic clearance: a model study.

Until now, the question whether slow desorption of compounds from transport proteins like the plasma protein albumin can affect hepatic uptake and thereby hepatic metabolism of these compounds has not yet been answered conclusively. This work now combines recently published experimental desorption rate constants with a liver model to address this question. For doing so, the used liver model differentiates the bound compound in blood, the unbound compound in blood and the compound within the hepatocytes as three well-stirred compartments. Our calculations show that slow desorption kinetics from albumin can indeed limit hepatic metabolism of a compound by decreasing hepatic extraction efficiency and hepatic clearance. The extent of this decrease, however, depends not only on the value of the desorption rate constant but also on how much of the compound is bound to albumin in blood and how fast intrinsic metabolism of the compound in the hepatocytes is. For strongly sorbing and sufficiently fast metabolized compounds, our calculations revealed a twentyfold lower hepatic extraction efficiency and hepatic clearance for the slowest known desorption rate constant compared to the case when instantaneous equilibrium between bound and unbound compound is assumed. The same desorption rate constant, however, has nearly no effect on hepatic extraction efficiency and hepatic clearance of weakly sorbing and slowly metabolized compounds. This work examines the relevance of desorption kinetics in various example scenarios and provides the general approach needed to quantify the effect of flow limitation, membrane permeability and desorption kinetics on hepatic metabolism at the same time.

Desorption kinetics of organic chemicals from albumin.

When present in blood, most chemicals tend to bind to the plasma protein albumin. For distribution into surrounding tissues, desorption from albumin is necessary, because only the unbound form of a chemical is assumed to be able to cross cell membranes. For metabolism of chemicals, the liver is a particularly important organ. One potentially limiting step for hepatic uptake of the chemicals is desorption from albumin, because blood passes the human liver within seconds. Desorption kinetics from albumin can thus be an important parameter for our pharmacokinetic and toxicokinetic understanding of chemicals. This work presents a dataset of measured desorption rate constants and reveals a possibility for their prediction. Additionally, the obtained extraction profiles directly indicate physiological relevance of desorption kinetics, because desorption of the test chemicals is still incomplete after time frames comparable to the residence time of blood in the liver.

Modeling Exposure in the Tox21 in Vitro Bioassays.

High-throughput in vitro bioassays are becoming increasingly important in the risk characterization of anthropogenic chemicals. Large databases gather nominal effect concentrations (Cnom) for diverse modes of action. However, the biologically effective concentration can substantially deviate due to differences in chemical partitioning. In this study, we modeled freely dissolved (Cfree), cellular (Ccell), and membrane concentrations (Cmem) in the Tox21 GeneBLAzer bioassays for a set of neutral and ionogenic organic chemicals covering a large physicochemical space. Cells and medium constituents were experimentally characterized for their lipid and protein content, and partition constants were either collected from the literature or predicted by mechanistic models. The chemicals exhibited multifaceted partitioning to proteins and lipids with distribution ratios spanning over 8 orders of magnitude. Modeled Cfree deviated over 5 orders of magnitude from Cnom and can be compared to in vivo effect data, environmental concentrations, and the unbound fraction in plasma, which is needed for the in vitro to in vivo extrapolation. Ccell was relatively constant for chemicals with membrane lipid-water distribution ratios of 1000 or higher and proportional to Cnom. Representing a sum parameter for exposure that integrates the entire dose from intracellular partitioning, Ccell is particularly suitable for the effect characterization of chemicals with multiple target sites and the calculation of their relative effect potencies. Effective membrane concentrations indicated that the specific effects of very hydrophobic chemicals in multiple bioassays are occurring at concentrations close to baseline toxicity. The equilibrium partitioning model including all relevant system parameters and a generic bioassay setup is attached as an excel workbook to this paper and can readily be applied to diverse in vitro bioassays.

Partitioning of Organic Ions to Muscle Protein: Experimental Data, Modeling, and Implications for in Vivo Distribution of Organic Ions.

The in vivo partitioning behavior of ionogenic organic chemicals (IOCs) is of paramount importance for their toxicokinetics and bioaccumulation. Among other proteins, structural proteins including muscle proteins could be an important sorption phase for IOCs, because of their high quantity in the human and other animals' body and their polar nature. Binding data for IOCs to structural proteins are, however, severely limited. Therefore, in this study muscle protein-water partition coefficients (KMP/w) of 51 systematically selected organic anions and cations were determined experimentally. A comparison of the measured KMP/w with bovine serum albumin (BSA)-water partition coefficients showed that anionic chemicals sorb more strongly to BSA than to muscle protein (by up to 3.5 orders of magnitude), while cations sorb similarly to both proteins. Sorption isotherms of selected IOCs to muscle protein are linear (i.e., KMP/w is concentration independent), and KMP/w is only marginally influenced by pH value and salt concentration. Using the obtained data set of KMP/w a polyparameter linear free energy relationship (PP-LFER) model was established. The derived equation fits the data well (R(2) = 0.89, RMSE = 0.29). Finally, it was demonstrated that the in vitro measured KMP/w values of this study have the potential to be used to evaluate tissue-plasma partitioning of IOCs in vivo.

Equilibrium Sorption of Structurally Diverse Organic Ions to Bovine Serum Albumin.

Reliable partitioning data are essential for assessing the bioaccumulation potential and the toxicity of chemicals. In contrast to neutral organic chemicals, the partitioning behavior of ionogenic organic chemicals (IOCs) is still a black box for environmental scientists. Partitioning to serum albumin, the major protein in blood plasma, strongly influences the freely dissolved concentration of many chemicals (including IOCs), which affects their transport and distribution in the body. Because consistent data sets for partitioning of IOCs are rarely available, bovine serum albumin-water partition coefficients (KBSA/w) were measured in this study for 45 anionic and 4 cationic organic chemicals, including various substituted benzoic and naphthoic acids, sulfonates and several pesticides and pharmaceuticals. The results of this study suggest that binding to BSA is substantially influenced by the three-dimensional structure of the chemicals and the position of substitutions on the sorbing molecules. For example, we found a difference of >1.5 log units between isomeric chemicals such as 3,4-dichlorobenzoic acid and 2,6-dichlorobenzoic acid, and 1-naphthoic acid and 2-naphthoic acid. Conventional modeling approaches (e.g., based on octanol-water partition coefficients) poorly predict log KBSA/w of organic ions (R(2) ≤ 0.5), partially because they do not capture the observed steric effects. Hence, alternative modeling strategies will be required for accurate prediction of serum albumin-water partition coefficients of organic ions.

Predicting storage-lipid water partitioning of organic solutes from molecular structure.

Partitioning to storage fat is the major process for bioaccumulation of many neutral organic chemicals. In this work, we evaluated the performance of four predictive models, ABSOLV, COSMOtherm, KOWWIN, and SPARC to calculate storage lipid-water partition coefficients. In a first step of the validation, we used over 300 literature data for chemicals with relatively simple molecular structures. For these compounds the overall performance was similar for all models with a root-mean-square error (rmse) between 0.45 and 0.61 log units. Clear differences became visible in the second validation step where a subset with only H-bond-donor compounds was used. Here, COSMOtherm and SPARC performed clearly better with an rmse of 0.35 and 0.42 log units, respectively, compared to ABSOLV and KOWWIN with an rmse of 0.91 and 0.85 log units, respectively. The last step in our validation was a comparison with experimental values for 22 complex, multifunctional chemicals (including pesticides, hormones, mycotoxins) that we measured specifically for this validation purpose. For these chemicals, predictions by all models were less accurate than those for simpler chemicals. COSMOtherm performed the best (rmse 0.71 log units) while the other methods showed considerably poorer results (rmse 1.29 (ABSOLV), 1.25 (SPARC), and 1.62 (KOWWIN) log units).

Calculating Equilibrium Phase Distribution during the Formation of Secondary Organic Aerosol Using COSMOtherm.

Challenges in the parametrization of compound distribution between the gas and particle phase contribute significantly to the uncertainty in the prediction of secondary organic aerosol (SOA) formation and are rooted in the complexity and variability of atmospheric condensed matter, which includes water, salts, and a multitude of organic oxidation products, often in two separated phases. Here, we explore the use of the commercial quantum-chemistry-based software COSMOtherm to predict equilibrium partitioning and Setchenow coefficients of a suite of oxidation products of α-pinene ozonolysis in an aerosol that is assumed to separate into an organic-enriched phase and an electrolyte-enriched aqueous phase. The predicted coefficients are used to estimate the phase distribution of the organic compounds, water and ammonium sulfate, the resulting phase composition, and the SOA yield. Four scenarios that differ in terms of organic loading, liquid water content, and chemical aging are compared. The organic compounds partition preferentially to the organic phase rather than the aqueous phase for the studied aerosol scenarios, partially due to the salting-out effect. Extremely low volatile organic compounds are predicted to be the dominant species in the organic aerosols at low loadings and an important component at higher loadings. The highest concentration of oxidation products in the condensed phase is predicted for a scenario assuming the presence of non-phase-separated cloud droplets. Partitioning into an organic aerosol phase composed of the oxidation products is predicted to be similar to partitioning into a phase composed of a single organic surrogate molecule, suggesting that the calculation procedure can be simplified without major loss of accuracy. COSMOtherm is shown to produce results that are comparable to those obtained using group contribution methods. COSMOtherm is likely to have a much larger application domain than those group contribution methods because it is based on fundamental principles with little calibration.

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.

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.

Effects of Aqueous Boundary Layers and Paracellular Transport on the Efflux Ratio as a Measure of Active Transport Across Cell Layers.

The efflux ratio (ER), determined by Caco-2/MDCK assays, is the standard in vitro metric to establish qualitatively whether a compound is a substrate of an efflux transporter. However, others have also enabled the utilisation of this metric quantitatively by deriving a relationship that expresses the ER as a function of the intrinsic membrane permeability of the membrane (P0) as well as the permeability of carrier-mediated efflux (Ppgp). As of yet, Ppgp cannot be measured directly from transport experiments or otherwise, but the ER relationship provides easy access to this value if P0 is known. However, previous derivations of this relationship failed to consider the influence of additional transport resistances such as the aqueous boundary layers (ABLs) and the filter on which the monolayer is grown. Since single fluxes in either direction can be heavily affected by these experimental artefacts, it is crucial to consider the potential impact on the ER. We present a model that includes these factors and show both mathematically and experimentally that this simple ER relationship also holds for the more realistic scenario that does not neglect the ABLs/filter. Furthermore, we also show mathematically how paracellular transport affects the ER, and we experimentally confirm that paracellular dominance reduces the ER to unity and can mask potential efflux.