In vitro to in vivo extrapolation of biotransformation rates for assessing bioaccumulation of hydrophobic organic chemicals in mammals

Methods for assessing the bioaccumulation of hydrocarbons and related substances in terrestrial organisms: A critical review

This study investigates and reviews methods for the assessment of the terrestrial bioaccumulation potential of hydrocarbons and related organic substances. The study concludes that the unitless biomagnification factor (BMF) and/or the trophic magnification factor (TMF) are appropriate, practical, and thermodynamically meaningful metrics for identifying bioaccumulative substances in terrestrial food chains. The study shows that various methods, including physical-chemical properties like the KOA and KOW , in vitro biotransformation assays, quantitative structure-activity relationships, in vivo pharmacokinetic and dietary bioaccumulation tests, and field-based trophic magnification studies, can inform on whether a substance has the potential to biomagnify in a terrestrial food chain as defined by a unitless BMF exceeding 1. The study further illustrates how these methods can be arranged in a four-tier evaluation scheme for the purpose of screening assessments that aim to minimize effort and costs and expediate bioaccumulation assessment of the vast numbers of organic substances in commerce, identifies knowledge gaps, and provides recommendations for further research to improve bioaccumulation assessment. Integr Environ Assess Manag 2023;19:1433-1456. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Models Used to Predict Chemical Bioaccumulation in Fish from in Vitro Biotransformation Rates Require Accurate Estimates of Blood-Water Partitioning and Chemical Volume of Distribution

Methods for extrapolating measured in vitro intrinsic clearance to a whole-body biotransformation rate constant (kB ) have been developed to support modeled bioaccumulation assessments for fish. The inclusion of extrapolated kB values into existing bioaccumulation models improves the prediction of chemical bioconcentration factors (BCFs), but there remains a tendency for these methods to overestimate BCFs relative to measured values. Therefore, a need exists to evaluate the extrapolation procedure to assess potential sources of error in predicted kB values. We examined how three different approaches (empirically based, composition based, and polyparameter linear free energy relationships [ppLFERs]) used to predict chemical partitioning in vitro (liver S9 system; KS9W ), in blood (KBW ), and in whole fish tissues (KFW ) impact the prediction of a chemical's hepatic clearance binding term (fU ) and apparent volume of distribution (VD ), both of which factor into the calculation of kB and the BCF. Each approach yielded different KS9W , KBW , and KFW values, but resulted in fU values that were of similar magnitude and remained relatively constant at log octanol-water partition ratios (KOW ) greater than 4. This is because KBW and KS9W values predicted by any given approach exhibit a similar dependence on log KOW (i.e., regression slope), which results in a cancelation of "errors" when fU is calculated. In contrast, differences in KBW values predicted by the three approaches translate to differences in VD , and by extension kB and the BCF, which become most apparent at log KOW greater than 6. There is a need to collect KBW and VD data for hydrophobic chemicals in fish that can be used to evaluate and improve existing partitioning prediction approaches in extrapolation models for fish. Environ Toxicol Chem 2023;42:33-45. © 2022 SETAC.

Developing an internal threshold of toxicological concern (iTTC)

BACKGROUND

Threshold of Toxicological Concern (TTC) approaches are used for chemical safety assessment and risk-based priority setting for data poor chemicals. TTCs are derived from in vivo No Observed Effect Level (NOEL) datasets involving an external administered dose from a single exposure route, e.g., oral intake rate. Thus, a route-specific TTC can only be compared to a route-specific exposure estimate and such TTCs cannot be used for other exposure scenarios such as aggregate exposures.

OBJECTIVE

Develop and apply a method for deriving internal TTCs (iTTCs) that can be used in chemical assessments for multiple route-specific exposures (e.g., oral, inhalation or dermal) or aggregate exposures.

METHODS

Chemical-specific toxicokinetics (TK) data and models are applied to calculate internal concentrations (whole-body and blood) from the reported administered oral dose NOELs used to derive the Munro TTCs. The new iTTCs are calculated from the 5th percentile of cumulative distributions of internal NOELs and the commonly applied uncertainty factor of 100 to extrapolate animal testing data for applications in human health assessment.

RESULTS

The new iTTCs for whole-body and blood are 0.5 nmol/kg and 0.1 nmol/L, respectively. Because the iTTCs are expressed on a molar basis they are readily converted to chemical mass iTTCs using the molar mass of the chemical of interest. For example, the median molar mass in the dataset is 220 g/mol corresponding to an iTTC of 22 ng/L-blood (22 pg/mL-blood). The iTTCs are considered broadly applicable for many organic chemicals except those that are genotoxic or acetylcholinesterase inhibitors. The new iTTCs can be compared with measured or estimated whole-body or blood exposure concentrations for chemical safety screening and priority-setting.

SIGNIFICANCE

Existing Threshold of Toxicological Concern (TTC) approaches are limited in their applications for route-specific exposure scenarios only and are not suitable for chemical risk and safety assessments under conditions of aggregate exposure. New internal Threshold of Toxicological Concern (iTTC) values are developed to address data gaps in chemical safety estimation for multi-route and aggregate exposures.

Bioaccumulation Screening of Neutral Hydrophobic Organic Chemicals in Air-Breathing Organisms Using In Vitro Rat Liver S9 Biotransformation Assays

To advance methods for bioaccumulation assessment of organic substances in air-breathing organisms, the present study developed an in vitro approach for screening neutral hydrophobic organic substances for their bioaccumulation potential in air-breathing organisms consisting of (1) depletion assays for chemicals in rat liver S9 subcellular fractions, (2) in vitro-in vivo extrapolation, and (3) whole-organism bioaccumulation modeling to assess the biomagnification potential of neutral organic substances in the rat. Testing of the in vitro method on 14 test chemicals of potentially biomagnifying substances showed that the bioassays could be conducted with a high level of reproducibility and that in vitro-derived elimination rate constants were in good agreement with in vivo-determined elimination rate constants in the rat. Exploring the potential of the in vitro approach for screening organic chemicals for bioaccumulation in air-breathing organisms indicated that chemical substances that exhibit a depletion rate constant in the S9 in vitro bioassay ≥0.3 h^-1 are not expected to biomagnify in rats independent of their octanol-water partitioning coefficient (K_OW ) or octanol-air partitioning coefficient (K_OA ). The high level of reproducibility achieved in the test, combined with the good agreement between in vitro-derived and in vivo-determined depuration rates, suggests that the in vitro approach in combination with a K_OA - and K_OW -based screening approach has good potential for screening chemicals in commerce for their bioaccumulation potential in air-breathing organisms in a cost-effective and expedient manner, especially if the bioassay can be automated. Environ Toxicol Chem 2022;41:2565-2579. © 2022 SETAC.

Relevance of desorption kinetics and permeability for in vitro-based predictions of hepatic clearance in fish

The impact of desorption kinetics and permeation kinetics on in vitro-based predictions of in vivo hepatic blood clearances is investigated in the present study. Most commonly, possible limitations due to slow desorption of chemicals from albumin or slow permeation of chemicals through cellular membranes are not considered when in vivo clearances are predicted from in vitro biotransformation rate constants. To evaluate whether the most commonly used extrapolation models might thus overlook important kinetic limitations, we compare predictions of in vivo clearance that explicitly consider desorption and permeation kinetics with predictions of in vivo clearance that neglect these aspects. Our results show that strong limitations due to slow permeation kinetics are possible depending on the assumed permeability value. While permeability values estimated with a mechanistic approach are fast enough to avoid significant limitations, other experimentally derived permeability values lead to dramatically decreased in vivo clearance predictions. These latter values lead to unrealistically low in vivo biotransformation estimates. Furthermore, we also evaluated the implications of desorption kinetics using experimentally determined desorption rate constants. These evaluations show that slow desorption kinetics are unlikely to limit in vivo clearance.

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.

Comparison of a simple and a complex model for BCF prediction using in vitro biotransformation data

A promising approach for bioaccumulation assessment with reduced animal use is the prediction of bioconcentration factors (BCFs) using in vitro biotransformation data. However, it has been recognized that the BCFs predicted using current models often are in poor agreement with experimental BCFs. Furthermore, extrahepatic biotransformation (e.g. in gill or GIT) is usually not accounted for. Here, we compare two BCF prediction models: a simple one-compartment and a more advanced multi-compartment model. Both models are implemented in a two-in-one calculation tool for the prediction of BCFs using in vitro data. Furthermore, both models were set up in a way that in vitro data for extrahepatic biotransformation can be easily considered, if desired. The models differ in their complexity: the one-compartment model is attractive because its simplicity, while the multi-compartment model is characterized by its refined closeness to reality. A comparison of the results shows that both models yield almost identical results for the presently evaluated cases with plausible physiological data. For regulatory purposes, there is thus no reason not to use the simple one-compartment model. However, if it is desired to represent special in vivo characteristics, e.g. first-pass effects or the direct GIT-to-liver blood flow, the multi-compartment model should be used.

Comparison of Alternative Methods for Bioaccumulation Assessment: Scope and Limitations of In Vitro Depletion Assays with Rainbow Trout and Bioconcentration Tests in the Freshwater Amphipod Hyalella azteca

Bioaccumulation assessment predominantly relies on the bioconcentration factor (BCF) as the sole decisive metric. The test guideline 305 by the Organisation for Economic Co-operation and Development (OECD) provides the standard procedure for deriving this in vivo fish BCF, which is not only expensive and labor-intensive, but also requires many animals. Accordingly, there is a great need for and interest in alternative methods that can help to reduce, replace, and refine vertebrate tests, as described in the 3R principles. Two alternative approaches have been developed: the bioconcentration test with the freshwater amphipod Hyalella azteca and the OECD test guideline 319 which provides a method to determine experimentally derived in vitro metabolism rates that can then be incorporated into in silico prediction models for rainbow trout BCF calculation. In the present study both alternative methods were applied to 5 substances of different physicochemical characteristics. The results were compared with literature values of fish in vivo BCFs and additional BCFs obtained with the alternative methods, if available. Potential differences between the results of the test methods are discussed utilizing information such as in vivo metabolism rates. The currently available data set suggests that these 2 alternative methods pose promising alternatives to predict bioaccumulation in fish, although defined applicability domains have yet to be determined. Environ Toxicol Chem 2020;39:1813-1825. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Examining Uncertainty in In Vitro-In Vivo Extrapolation Applied in Fish Bioconcentration Models

In vitro biotransformation rates were determined for 30 chemicals, mostly fragrance ingredients, using trout liver S9 fractions (RT-S9) and incorporated into in vitro-in vivo extrapolation (IVIVE) models to predict bioconcentration factors (BCFs). Predicted BCFs were compared against empirical BCFs to explore potential major uncertainties involved in the in vitro methods and IVIVE models: (i) in vitro chemical test concentrations; (ii) different gill uptake rate constant calculations (k₁); (iii) protein binding (different calculations and measurement of the fraction of unbound chemical, f_U); (iv) species differences; and (v) extrahepatic biotransformation. Predicted BCFs were within 0.5 log units for 44% of the chemicals compared to empirical BCFs, whereas 56% were overpredicted by >0.5 log units. This trend of overprediction was reduced by alternative k₁ calculations to 32% of chemicals being overpredicted. Moreover, hepatic in vitro rates scaled to whole body biotransformation rates (k_B) were compared against in vivo k_B estimates. In vivo k_B was underestimated for 79% of the chemicals. Neither lowering the test concentration, nor incorporation of new measured f_U values, nor species matching avoided the tendency to overpredict BCFs indicating that further improvements to the IVIVE models are needed or extrahepatic biotransformation plays an underestimated role.

Biotransformation of Polycyclic Aromatic Hydrocarbons by Trout Liver S9 Fractions: Evaluation of Competitive Inhibition Using a Substrate Depletion Approach

Environmental contaminants frequently occur as part of a chemical mixture, potentially resulting in competitive inhibition among multiple substrates metabolized by the same enzyme. Trout liver S9 fractions were used to evaluate the biotransformation of 3 polycyclic aromatic hydrocarbons (PAHs): phenanthrene, pyrene, and benzo[a]pyrene, tested as binary mixtures. Initial rates of biotransformation were determined using a substrate-depletion approach. The resulting data were then fitted by simultaneous nonlinear regression to a competitive inhibition model. In each case, the PAH possessing the lower Michaelis-Menten affinity constant (K_M ) competitively inhibited biotransformation of the other compound. Inhibition constants determined for the lower-K_M compound were generally close to previously determined K_M values, consistent with the suggestion that phase I biotransformation of PAHs is largely catalyzed by one or a small number of cytochrome P450 enzymes. The use of a substrate-depletion approach to perform enzyme-inhibition studies imposes practical limitations on experimental design and complicates the interpretation of derived kinetic constants. Nevertheless, the resulting information may have utility for chemical hazard assessments as well as the design and interpretation of controlled laboratory studies. Depletion experiments informed by measured chemical concentrations in tissues may also provide a means of determining whether enzyme inhibition occurs under relevant environmental conditions. Environ Toxicol Chem 2019;38:2729-2739. Published 2019 Wiley Periodicals, Inc. on behalf of SETAC. This article is a US government work, and as such, is in the public domain in the United States of America.

Concentration dependence of in vitro biotransformation rates of hydrophobic organic sunscreen agents in rainbow trout S9 fractions: Implications for bioaccumulation assessment

In vitro biotransformation studies were performed to support the bioaccumulation assessment of 3 hydrophobic organic ultraviolet filters (UVFs), 4-methylbenzylidene camphor (4-MBC), 2-ethylhexyl-4-methoxycinnamate (EHMC), and octocrylene. In vitro depletion rate constants (kdep ) were determined for each UVF using rainbow trout liver S9 fractions. Incubations performed with and without added cofactors showed complete (4-MBC) or partial (EHMC and octocrylene) dependence of kdep on addition of the reduced form of nicotinamide adenine dinucleotide phosphate (NADPH), suggesting that hydrolysis of EHMC and octocrylene by NADPH-independent enzymes (e.g., carboxylesterases) is an important metabolic route. The concentration dependence of kdep was then evaluated to estimate Michaelis-Menten parameters (KM and Vmax ) for each UVF. Measured kdep values were then extrapolated to apparent whole-body biotransformation rate constants using an in vitro-in vivo extrapolation (IVIVE) model. Bioconcentration factors (BCFs) calculated from kdep values measured at concentrations well below KM were closer to empirical BCFs than those calculated from kdep measured at higher test concentrations. Modeled BCFs were sensitive to in vitro binding assumptions employed in the IVIVE model, highlighting the need for further characterization of chemical binding effects on hepatic clearance. The results suggest that the tested UVFs are unlikely to accumulate to levels exceeding the European Union Registration, Evaluation, Authorisation, and Restriction regulation criterion for bioaccumulative substances (BCF > 2000 L kg-1 ). However, consideration of appropriate in vitro test concentrations and binding correction factors are important when IVIVE methods are used to refine modeled BCFs. Environ Toxicol Chem 2019;38:548-560. © 2018 SETAC.

Development of three-dimensional (3D) spheroid cultures of the continuous rainbow trout liver cell line RTL-W1

In vitro experimental systems based on continuous piscine cell lines can be used as an alternative to animal tests for obtaining qualitative and quantitative information on the possible fate and effect of chemicals in fish. However, their capability to reproduce complex metabolic processes and toxic responses as they occur in vivo is limited due to the lack of organ-specific tissue architecture and functions. Here we introduce a three-dimensional (3D) in vitro experimental system based on spheroidal aggregate cultures (spheroids) of the continuous rainbow trout liver cell line RTL-W1 and provide a first description of their structural and functional properties including growth, viability/longevity, metabolic activity, ultrastructure and cytochrome P450 1A (CYP1A) expression determined by bright-field, multi-photon fluorescence and transmission electron microscopy as well as RT-qPCR analysis. Our results show that RTL-W1 cells in 3D spheroids (ø ~ 150 µm) (including those in the interior) were viable, metabolically active and had higher basal and β-naphthoflavone-induced CYP1A expression levels than conventional 2D cell cultures. Furthermore, they displayed ultrastructural characteristics similar to differentiated hepatocytes. The available evidence suggests that 3D RTL-W1 spheroids may have enhanced hepatotypic functions and be a superior in vitro model to assess hepatic biotransformation, bioaccumulation and chronic toxicity compared to conventional cell monolayer cultures.