A physiologically based toxicokinetic model for the uptake and disposition of waterborne organic chemicals in fish

Adverse effects of Microcystis aeruginosa exudates on the filtration, digestion, and reproduction organs of benthic bivalve Corbicula fluminea

Cyanobacteria bloom and the secondary metabolites released by the microorganism are extremely harmful to aquatic animals, yet study on their adverse effects in zoobenthos is rare. Corbicula fluminea widely distributed in freshwater environment with algal blooms. It is a typical filter feeding zoobenthos that may be affected by the secondary metabolites of cyanobacteria due to its high filtering rate. In this study, C. fluminea was exposed to Microcystis aeruginosa exudates (MaE) for 96 h, which was obtained from 5 × 105 cells/mL and 2.5 × 106 cells/mL exponential stage M. aeruginosa culture solution that represented cyanobacteria cell density needs environmental risk precaution control and emergent control, respectively. The responses of C. fluminea critical organs to MaE were analyzed and evaluated based on histopathological sections, antitoxicity biomarkers, and organ function biomarkers. The results showed that all the organs underwent structural disorders, cell vacuolization, apoptosis, and necrosis, and the damage levels increased as MaE concentration increased. The detoxification and antioxidant defense systems biomarkers in each organ response to MaE exposure differently and the level of reaction improved when MaE concentration increased. The siphon rate and acetylcholinesterase activity showed that the filtration function decreased significantly as the MaE concentration increased. Increased activity of glutathione S-transferase and amylase in the digestive gland indicate that it is the major detoxification organ of C. fluminea. Increased vitellogenin concentration and enlarged oocytes in the gonad indicate that MaE may have an estrogenic effect on C. fluminea. This study demonstrates that cyanobacteria threat benthic bivalves by inducing oxidative stress, inhibiting filtering feeding system, and disturbing digestion system and reproduction potential of C. fluminea.

Interspecies scaling of toxicity reference values in human health versus ecological risk assessments: A critical review

Risk assessments that focus on anthropogenic chemicals in environmental media-whether considering human health or ecological effects-often rely on toxicity data from experimentally studied species to estimate safe exposures for species that lack similar data. Current default extrapolation approaches used in both human health risk assessments and ecological risk assessments (ERAs) account for differences in body weight between the test organisms and the species of interest, but the two default approaches differ in important ways. Human health risk assessments currently employ a default based on body weight raised to the three-quarters power. Ecological risk assessments for wildlife (i.e., mammals and birds) are typically based directly on body weight, as measured in the test organism and receptor species. This review describes differences in the experimental data underlying these default practices and discusses the many factors that affect interspecies variability in chemical exposures. The interplay of these different factors can lead to substantial departures from default expectations. Alternative methodologies for conducting more accurate interspecies extrapolations in ERAs for wildlife are discussed, including tissue-based toxicity reference values, physiologically based toxicokinetic and/or toxicodynamic modeling, chemical read-across, and a system of categorical defaults based on route of exposure and toxic mode of action. Integr Environ Assess Manag 2024;20:749-764. © 2023 SETAC.

Cross-Species Extrapolation of Biological Data to Guide the Environmental Safety Assessment of Pharmaceuticals-The State of the Art and Future Priorities

The extrapolation of biological data across species is a key aspect of biomedical research and drug development. In this context, comparative biology considerations are applied with the goal of understanding human disease and guiding the development of effective and safe medicines. However, the widespread occurrence of pharmaceuticals in the environment and the need to assess the risk posed to wildlife have prompted a renewed interest in the extrapolation of pharmacological and toxicological data across the entire tree of life. To address this challenge, a biological "read-across" approach, based on the use of mammalian data to inform toxicity predictions in wildlife species, has been proposed as an effective way to streamline the environmental safety assessment of pharmaceuticals. Yet, how effective has this approach been, and are we any closer to being able to accurately predict environmental risk based on known human risk? We discuss the main theoretical and experimental advancements achieved in the last 10 years of research in this field. We propose that a better understanding of the functional conservation of drug targets across species and of the quantitative relationship between target modulation and adverse effects should be considered as future research priorities. This pharmacodynamic focus should be complemented with the application of higher-throughput experimental and computational approaches to accelerate the prediction of internal exposure dynamics. The translation of comparative (eco)toxicology research into real-world applications, however, relies on the (limited) availability of experts with the skill set needed to navigate the complexity of the problem; hence, we also call for synergistic multistakeholder efforts to support and strengthen comparative toxicology research and education at a global level. Environ Toxicol Chem 2024;43:513-525. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Fish Physiologically Based Toxicokinetic Modeling Approach for In Vitro-In Vivo and Cross-Species Extrapolation of Endocrine-Disrupting Chemicals in Risk Assessment

High-throughput in vitro assays combined with in vitro-in vivo extrapolation (IVIVE) leverage in vitro responses to predict the corresponding in vivo exposures and thresholds of concern. The integrated approach is also expected to offer the potential for efficient tools to provide estimates of chemical toxicity to various wildlife species instead of animal testing. However, developing fish physiologically based toxicokinetic (PBTK) models for IVIVE in ecological applications is challenging, especially for plausible estimation of an internal effective dose, such as fish equivalent concentration (FEC). Here, a fish PBTK model linked with the IVIVE approach was established, with parameter optimization of chemical unbound fraction, pH-dependent ionization and hepatic clearance, and integration of temperature effect and growth dilution. The fish PBTK-IVIVE approach provides not only a more precise estimation of tissue-specific concentrations but also a reasonable approximation of FEC targeting the estrogenic potency of endocrine-disrupting chemicals. Both predictions were compared with in vivo data and were accurate for most indissociable/dissociable chemicals. Furthermore, the model can help determine cross-species variability and sensitivity among the five fish species. Using the available IVIVE-derived FEC with target pathways is helpful to develop predicted no-effect concentration for chemicals with similar mode of action and support screening-level ecological risk assessment.

Non-negligible vector effect of micro(nano)plastics on tris(1,3-dichloro-2-propyl) phosphate in zebrafish quantified by toxicokinetic model

Micro(nano)plastics (MNPs) inevitably interact with coexisting contaminants and can act as vectors to affect their fate in organisms. However, the quantitative contribution of MNPs in the in vivo bioaccumulation and distribution of their coexisting contaminants remains unclear. Here, by selecting tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) as the typical coexisting contaminant, we quantified the contribution of MNPs to bioaccumulation and distribution of TDCIPP with toxicokinetic models. Results indicated that MNPs differentially facilitated TDCIPP bioaccumulation and distribution, and NPs slowed down TDCIPP depuration more significantly than MPs. Model analysis further revealed increasing contributions of MNPs to whole-fish TDCIPP bioaccumulation over time, with NPs (33-42%) contributing more than MPs (12-32%) at 48 h exposure. NPs contributed more than MPs to TDCIPP distribution in the liver (13-19% for MPs; 36-52% for NPs) and carcass (24-45% for MPs; 57-71% for NPs). The size-dependent vector effect might be attributed to the fact that MNPs promote contaminant transfer by damaging biofilm structure and increasing tissue membrane permeability, with NPs exerting stronger effects. This work demonstrated the effectiveness of using modeling tools to understand the relative importance of MNPs as contaminant vectors in the TK process and highlighted the higher contaminant transfer potential of NPs under combined exposure scenarios.

Multimodal Model to Predict Tissue-to-Blood Partition Coefficients of Chemicals in Mammals and Fish

Tissue-to-blood partition coefficients (Ptb) are key parameters for assessing toxicokinetics of xenobiotics in organisms, yet their experimental data were lacking. Experimental methods for measuring Ptb values are inefficient, underscoring the urgent need for prediction models. However, most existing models failed to fully exploit Ptb data from diverse sources, and their applicability domain (AD) was limited. The current study developed a multimodal model capable of processing and integrating textual (categorical features) and numerical information (molecular descriptors/fingerprints) to simultaneously predict Ptb values across various species, tissues, blood matrices, and measurement methods. Artificial neural network algorithms with embedding layers were used for the multimodal modeling. The corresponding unimodal models were developed for comparison. Results showed that the multimodal model outperformed unimodal models. To enhance the reliability of the model, a method considering categorical features, weighted molecular similarity density, and weighted inconsistency in molecular activities of structure-activity landscapes was used to characterize the AD. The model constrained by the AD exhibited better prediction accuracy for the validation set, with the determination coefficient, root mean-square error, and mean absolute error being 0.843, 0.276, and 0.213 log units, respectively. The multimodal model coupled with the AD characterization can serve as an efficient tool for internal exposure assessment of chemicals.

In silico bioavailability triggers applied to direct and indirect thyroid hormone disruptors

Among endocrine disruption, interference with the thyroid hormone (TH) regulation is of increasing concern. Respective compounds encode through their structural features both the potential for TH disruption, and the bioavailability mitigating the toxicological effect. The aim of this study is to provide a substructure-based screening-level QSAR (quantitative structure-activity relationship) that discriminates bioavailable TH disruptors from not bioavailable counterparts, covering both direct and indirect (retinoid- and AhR-mediated) modes of action. The QSAR has been derived from literature data for 1642 compounds, and takes into account Lipinski's rule-of-five and the brain/blood partition coefficient Kbrain/blood. For its validation, an external test set of 145 substances has been used. For 1787 compounds meeting the model application domain, the model yields only one false negative. The discussion addresses the mechanistic meaning of the bioavailability triggers molecular weight, H-bond donor and acceptor, hydrophobicity (log Kow), and the physicochemical properties underlying log Kbrain/blood. The model may serve as bioavailability-screening step within a decision support system for the predictive assessment of chemicals regarding their potential to disrupt thyroid hormone function in a direct or indirect manner.

Can blood proteome diversity among fish species help explain perfluoroalkyl acid trophodynamics in aquatic food webs?

Per- and polyfluoroalkyl substances (PFAS) are a diverse family of industrially significant synthetic chemicals infamous for extreme environmental persistence and global environmental distribution. Many PFAS are bioaccumulative and biologically active mainly due to their tendency to bind with various proteins. These protein interactions are important in determining the accumulation potential and tissue distribution of individual PFAS. Trophodynamics studies including aquatic food webs present inconsistent evidence for PFAS biomagnification. This study strives to identify whether the observed variability in PFAS bioaccumulation potential among species could correspond with interspecies protein composition differences. Specifically, this work compares the perfluorooctane sulfonate (PFOS) serum protein binding potential and the tissue distribution of ten perfluoroalkyl acids (PFAAs) detected in alewife (Alosa pseudoharengus), deepwater sculpin (Myoxocephalus thompsonii), and lake trout (Salvelinus namaycush) of the Lake Ontario aquatic piscivorous food web. These three fish sera and fetal bovine reference serum all had unique total serum protein concentrations. Serum protein-PFOS binding experiments showed divergent patterns between fetal bovine serum and fish sera, suggesting potentially two different PFOS binding mechanisms. To identify interspecies differences in PFAS-binding serum proteins, fish sera were pre-equilibrated with PFOS, fractionated by serial molecular weight cut-off filter fractionation, followed by liquid chromatography-tandem mass spectrometry analysis of the tryptic protein digests and the PFOS extracts of each fraction. This workflow identified similar serum proteins for all fish species. However, serum albumin was only identified in lake trout, suggesting apolipoproteins are likely the primary PFAA transporters in alewife and deepwater sculpin sera. PFAA tissue distribution analysis provided supporting evidence for interspecies variations in lipid transport and storage, which may also contribute to the varied PFAA accumulation in these species. Proteomics data are available via ProteomeXchange with identifier PXD039145.

An amended in vitro-in vivo extrapolation model that accounts for first pass clearance effects on chemical bioaccumulation in fish

Measured rates of in vitro intrinsic clearance for fish may be extrapolated to the whole animal as a means of estimating a whole-body biotransformation rate constant (k_B; d^-1). This estimate of k_B can then be used as an input to existing bioaccumulation prediction models. Most in vitro-in vivo extrapolation/bioaccumulation (IVIVE/B) modeling efforts to date have focused on predicting the chemical bioconcentration in fish (aqueous only exposure), with less attention paid to dietary exposures. Following dietary uptake, biotransformation in the gut lumen, intestinal epithelia, and liver can reduce chemical accumulation; however, current IVIVE/B models do not consider these first pass clearance effects on dietary uptake. Here we present an amended IVIVE/B model that accounts for first pass clearance. The model is then used to examine how biotransformation in the liver and intestinal epithelia (alone or combined) may impact chemical accumulation that occurs during dietary exposure. First pass clearance by the liver can greatly reduce dietary uptake of contaminants, but these effects are only apparent at rapid rates of in vitro biotransformation (first order depletion rate constant k_DEP ≥ 10 h^-1). The impact of first pass clearance becomes more pronounced when biotransformation in the intestinal epithelia is included in the model. Modelled results suggest that biotransformation in the liver and intestinal epithelia cannot entirely explain reduced dietary uptake reported in several in vivo bioaccumulation tests. This unexplained reduction in dietary uptake is attributed to chemical degradation in the gut lumen. These findings underscore the need for research to directly investigate luminal biotransformation in fish.

An Overview of Physiologically-Based Pharmacokinetic Models for Forensic Science

A physiologically-based pharmacokinetic (PBPK) model represents the structural components of the body with physiologically relevant compartments connected via blood flow rates described by mathematical equations to determine drug disposition. PBPK models are used in the pharmaceutical sector for drug development, precision medicine, and the chemical industry to predict safe levels of exposure during the registration of chemical substances. However, one area of application where PBPK models have been scarcely used is forensic science. In this review, we give an overview of PBPK models successfully developed for several illicit drugs and environmental chemicals that could be applied for forensic interpretation, highlighting the gaps, uncertainties, and limitations.

Dermal Uptake is an Important Pathway for the Bioconcentration of Hydrophobic Organic Compounds by Zebrafish (Danio rerio)

For bioconcentration of hydrophobic organic compounds (HOCs), most of studies assumed that fish absorb HOCs mainly through gills but often ignored the dermal uptake. In this study, deuterated polycyclic aromatic hydrocarbons (PAHs-d₁₀, phenanthrene-d₁₀, and pyrene-d₁₀) and polychlorinated biphenyls (PCB-153) were selected to study whether zebrafish can absorb freely dissolved and dissolved organic matter (DOM)-associated HOCs through dermal uptake. The results showed that the freely dissolved PAHs and PCBs could directly enter the body of zebrafish through its skin. However, PAHs and PCB-153 associated with DOM (~ 10 kDa) could not enter zebrafish through the skin. When gill and dermal exposure coexisted, dermal uptake contributed 2.9 ~ 7.6% and 31.9 ~ 38.4% of PAHs and PCB-153 bioconcentration after exposure for 6 h, respectively. The present study demonstrates that dermal uptake is an important pathway for the bioconcentration of HOCs by fish, which should be considered when studying the toxicodynamics and toxicokinetics of HOCs in organisms.

Enhanced prediction of internal concentrations of phenolic endocrine disrupting chemicals and their metabolites in fish by a physiologically based toxicokinetic incorporating metabolism (PBTK-MT) model

Bisphenol A (BPA), 4-nonylphenol (4-NP), and triclosan (TCS) are phenolic endocrine disrupting chemicals (EDCs), which are widely detected in aquatic environments and further bioaccumulated and metabolized in fish. Physiologically based toxicokinetic (PBTK) models have been used to describe the absorption, distribution, metabolism, and excretion (ADME) of parent compounds in fish, whereas the metabolites are less explored. In this study, a PBTK incorporating metabolism (PBTK-MT) model for BPA, 4-NP, and TCS was established to enhance the performance of the traditional PBTK model. The PBTK-MT model comprised 16 compartments, showing great accuracy in predicting the internal concentrations of three compounds and their glucuronidated and sulfated conjugates in fish. The impact of typical hepatic metabolism on the PBTK-MT model was successfully resolved by optimizing the mechanism for deriving the partition coefficients between the blood and liver. The PBTK-MT model exhibited a potential data gap-filling capacity for unknown parameters through a backward extrapolation approach of parameters. Model sensitivity analysis suggested that only five parameters were sensitive in at least two PBTK-MT models, while most parameters were insensitive. The PBTK-MT model will contribute to a well understanding of the environmental behavior and risks of pollutants in aquatic biota.

Factors Affecting the Binding of Diltiazem to Rainbow Trout Plasma: Implications for the Risk Assessment of Pharmaceuticals in Aquatic Systems

The accumulation of organic toxicants in fish plasma, and how they partition between the bound and unbound fraction once absorbed, are important metrics in models that seek to predict the risk of such contaminants in aquatic settings. Rapid equilibrium dialysis of diltiazem, an ionizable weak base and important human pharmaceutical contaminant of freshwaters, was conducted with rainbow trout (Oncorhynchus mykiss) plasma. The effect of fed state, fish sex, fish strain/size, and dialysis buffer pH on the binding of radiolabeled diltiazem (9 ng ml^-1 ) was assessed. In fed fish, 24.6%-29.5% of diltiazem was free, unbound to plasma proteins. Although starvation of fish resulted in a decrease in plasma protein, the bound fraction of diltiazem remained relatively constant. Consequently, the protein-bound concentration of diltiazem increased with length of starvation. In general, rainbow trout strain was a significant factor affecting plasma binding, although the two strains tested also differed markedly in size. Dialysis buffer pH significantly influenced plasma binding, with a higher unbound diltiazem fraction at pH 6.8 than pH 8.0. These data indicate that empirical measures of plasma binding in fish are important for accurate risk assessment and that the physiological status of a fish is likely to impact its sensitivity to toxicants such as diltiazem. Environ Toxicol Chem 2022;41:3125-3133. © 2022 SETAC.

The toxicokinetics of bisphenol A and its metabolites in fish elucidated by a PBTK model

Bisphenol A (BPA) is a chemical of major concern due to its endocrine disrupting function, high production volume, and persistence in the aquatic environment. Consequently, organisms such as fish are subject to chronic exposure to BPA. However, physiologically-based toxicokinetic (PBTK) models, which are valuable tools to improve the understanding of a chemical's fate in an organism, have never been specifically adapted to model BPA toxicokinetics (TK) in fish. In our work, an existing PBTK developed for four different fish species was modified to model BPA ADME processes (absorption, distribution, metabolization and excretion). The metabolization of BPA into BPA-monoglucuronide (BPA gluc) and BPA-monosulfate (BPA sulf) and their TK in various organs was taking into account in the model. Experiments were performed to generate BPA TK data in a model species commonly used in ecotoxicology, the stickleback. The model structure had to include two sites of metabolization to simulate BPA TK accurately in stickleback organs. Thus, the fish liver may not be the only site of the metabolization of BPA: plasma or gills could also play a role in BPA metabolization. The PBTK model predictive performance evaluated on literature data in zebrafish and rainbow trout concurs with this conclusion. Finally, a calibration mixing data from the three species was compared to the calibration on stickleback data only.

A Generalized Physiologically Based Kinetic Model for Fish for Environmental Risk Assessment of Pharmaceuticals

An increasing number of pharmaceuticals found in the environment potentially impose adverse effects on organisms such as fish. Physiologically based kinetic (PBK) models are essential risk assessment tools, allowing a mechanistic approach to understanding chemical effects within organisms. However, fish PBK models have been restricted to a few species, limiting the overall applicability given the countless species. Moreover, many pharmaceuticals are ionizable, and fish PBK models accounting for ionization are rare. Here, we developed a generalized PBK model, estimating required parameters as functions of fish and chemical properties. We assessed the model performance for five pharmaceuticals (covering neutral and ionic structures). With biotransformation half-lives (HLs) from EPI Suite, 73 and 41% of the time-course estimations were within a 10-fold and a 3-fold difference from measurements, respectively. The performance improved using experimental biotransformation HLs (87 and 59%, respectively). Estimations for ionizable substances were more accurate than any of the existing species-specific PBK models. The present study is the first to develop a generalized fish PBK model focusing on mechanism-based parameterization and explicitly accounting for ionization. Our generalized model facilitates its application across chemicals and species, improving efficiency for environmental risk assessment and supporting an animal-free toxicity testing paradigm.

Screening the baseline fish bioconcentration factor of various types of surfactants using phospholipid binding data

Fish bioconcentration factors (BCFs) are commonly used in chemical hazard and risk assessment. For neutral organic chemicals BCFs are positively correlated with the octanol-water partition ratio (K_OW), but K_OW is not a reliable parameter for surfactants. Membrane lipid-water distribution ratios (D_MLW) can be accurately measured for all kinds of surfactants, using phospholipid-based sorbents. This study first demonstrates that D_MLW values for ionic surfactants are more than 100 000 times higher than the partition ratio to fish-oil, representing neutral storage lipid. A non-ionic alcohol ethoxylate surfactant showed almost equal affinity for both lipid types. Accordingly, a baseline screening BCF value for surfactants (BCF_baseline) can be approximated for ionic surfactants by multiplying D_MLW by the phospholipid fraction in tissue, and for non-ionic surfactants by multiplying D_MLW by the total lipid fraction. We measured D_MLW values for surfactant structures, including linear and branched alkylbenzenesulfonates, an alkylsulfoacetate and an alkylethersulfate, bis(2-ethylhexyl)-surfactants (e.g., docusate), zwitterionic alkylbetaines and alkylamine-oxides, and a polyprotic diamine. Together with sixty previously published D_MLW values for surfactants, structure-activity relationships were derived to elucidate the influence of surfactant specific molecular features on D_MLW. For 23 surfactant types, we established the alkyl chain length at which BCF_baseline would exceed the EU REACH bioaccumulation (B) threshold of 2000 L kg^-1, and would therefore require higher tier assessments to further refine the BCF estimate. Finally, the derived BCF_baseline are compared with measured literature in vivo BCF data where available, suggesting that refinements, most notably reliable estimates of biotransformation rates, are needed for most surfactant types.

Toxicokinetic Models for Bioconcentration of Organic Contaminants in Two Life Stages of White Sturgeon (Acipenser transmontanus)

The white sturgeon (Acipenser transmontanus) is an endangered ancient fish species that is known to be particularly sensitive to certain environmental contaminants, partly because of the uptake and subsequent toxicity of lipophilic pollutants prone to bioconcentration as a result of their high lipid content. To better understand the bioconcentration of organic contaminants in this species, toxicokinetic (TK) models were developed for the embryo-larval and subadult life stages. The embryo-larval model was designed as a one-compartment model and validated using whole-body measurements of benzo[a]pyrene (B[a]P) metabolites from a waterborne exposure to B[a]P. A physiologically based TK (PBTK) model was used for the subadult model. The predictive power of the subadult model was validated with an experimental data set of four chemicals. Results showed that the TK models could accurately predict the bioconcentration of organic contaminants for both life stages of white sturgeon within 1 order of magnitude of measured values. These models provide a tool to better understand the impact of environmental contaminants on the health and the survival of endangered white sturgeon populations.

Tissue-Specific Accumulation, Biotransformation, and Physiologically Based Toxicokinetic Modeling of Benzotriazole Ultraviolet Stabilizers in Zebrafish (Danio rerio)

Benzotriazole ultraviolet stabilizers (BUVSs) are high-production-volume chemicals with ubiquitous occurrence in the aquatic environment. However, little is known about their bioconcentration and biotransformation, and physiologically based toxicokinetic (PBTK) models for BUVSs are lacking. This study selected six BUVSs for which experiments were performed with zebrafish (Danio rerio) exposed to two different levels (0.5 and 10 μg·L^-1). Higher kinetic bioconcentration factors (BCFs) were observed at the lower exposure level with environmental relevance, with BCF of 3.33 × 10³ L·kg^-1 for 2-(2-hydroxy-3,5-di-tert-butylphenyl)-5-chlorobenzotriazole (UV-327). This phenomenon was interpreted by a nonlinear adsorption mechanism, where binding with specific protein sites contributes to bioconcentration. Muscle exhibited the lowest accumulation, in which depuration half-life of UV-327 was 19.5 d. In kidney, muscle, ovary, gill, and skin, logBCF increased with increase in log K_OW of the BUVSs until log K_OW was ca. 6.5, above which logBCF decreased. However, the trend was not observed in the liver and intestine. Six biotransformation products were identified and mainly accumulated in the liver and intestine. Considering the nonlinear adsorption mechanism in the PBTK model, the prediction accuracy of the model was improved, highlighting the binding of xenobiotics with specific protein sites in assessing the bioconcentration of chemicals for their risk assessment.

A Novel Multispecies Toxicokinetic Modeling Approach in Support of Chemical Risk Assessment

Standardized laboratory tests with a limited number of model species are a key component of chemical risk assessments. These surrogate species cannot represent the entire diversity of native species, but there are practical and ethical objections against testing chemicals in a large variety of species. In previous research, we have developed a multispecies toxicokinetic model to extrapolate chemical bioconcentration across species by combining single-species physiologically based toxicokinetic (PBTK) models. This "top-down" approach was limited, however, by the availability of fully parameterized single-species models. Here, we present a "bottom-up" multispecies PBTK model based on available data from 69 freshwater fishes found in Canada. Monte Carlo-like simulations were performed using statistical distributions of model parameters derived from these data to predict steady-state bioconcentration factors (BCFs) for a set of well-studied chemicals. The distributions of predicted BCFs for 1,4-dichlorobenzene and dichlorodiphenyltrichloroethane largely overlapped those of empirical data, although a tendency existed toward overestimation of measured values. When expressed as means, predicted BCFs for 26 of 34 chemicals (82%) deviated by less than 10-fold from measured data, indicating an accuracy similar to that of previously published single-species models. This new model potentially enables more environmentally relevant predictions of bioconcentration in support of chemical risk assessments.

Combined sediment desorption and bioconcentration model to predict levels of dioxin-like chemicals in fish

Flooding and other sediment disturbances can lead to increases in sediment resuspension. In this context, it is of central importance to understand the kinetics of release from these sediments and the uptake of pollutants, such as polychlorinated biphenyls (PCBs) and polychlorinated dioxins and furans (PCDD/Fs), into aquatic organisms. In the present study, we parameterized a sediment desorption model based on experimentally determined rapidly-desorbing fractions of dioxin-like chemicals (DLCs). We coupled this desorption model with a physiologically-based toxicokinetic model for rainbow trout. This combined model was used to predict DLC concentrations in the muscle of exposed fish. The performance of this model was evaluated using a previously published dataset on DLC uptake from sediment suspensions during simulated re-suspension events. Predictions generally differed less than 10-fold from measured values, and the model showed a good global coefficient of determination (R²) of 0.95. The root mean squared error (RMSE) for PCBs was 0.31 log units and 0.53 log units for PCDD/Fs. The results of our study demonstrate that the prediction of bioconcentration and related risk to fish resulting from sediment resuspension can be accurately predicted using coupled desorption and toxicokinetic models.

Development and Application of a Water Temperature Related Physiologically Based Pharmacokinetic Model for Enrofloxacin and Its Metabolite Ciprofloxacin in Rainbow Trout

Enrofloxacin (ENR) has been approved for the treatment of infections in aquaculture, but it may cause tissue residue. This research aimed to develop and validate a water temperature related PBPK model, including both ENR and ciprofloxacin (CIP), in rainbow trout, and to predict further their residue concentrations and the withdrawal periods for ENR at different water temperatures. With the published concentrations data, a flow-limited PBPK model including both ENR and CIP sub-models was developed to predict ENR and CIP concentrations in tissues and plasma/serum after intravenous, oral, or immersion administration. A Monte Carlo simulation including 500 iterations was further incorporated into this model. Based on the model and Monte Carlo analysis, the withdrawal intervals were estimated for different dosage regimens and at different water temperatures, ranging from 80 to 272 degree-days. All of these values were shorter than the labeled withdrawal period (500 degree-days) in fish. This model provided a useful tool for predicting the tissue residues of ENR and CIP in rainbow trout under different dosage regimens and at different water temperatures.

In vitro-in vivo and cross-life stage extrapolation of uptake and biotransformation of benzo[a]pyrene in the fathead minnow (Pimephales promelas)

Understanding internal dose metrics is integral to adequately assess effects environmental contaminants might have on aquatic wildlife, including fish. In silico toxicokinetic (TK) models are a leading approach for quantifying internal exposure metrics for fishes; however, they often do not adequately consider chemicals that are actively biotransformed and have not been validated against early-life stages (ELS) that are often considered the most sensitive to the exposure to contaminants. To address these uncertainties, TK models were parameterized for the rapidly biotransformed chemical benzo[a]pyrene (B[a]P) in embryo-larval and adult life stages of fathead minnows. Biotransformation of B[a]P was determined through measurements of in vitro clearance. Using in vitro-in vivo extrapolation, in vitro clearance was integrated into a multi-compartment TK model for adult fish and a one-compartment model for ELS. Model predictions were validated using measurements of B[a]P metabolites from in vivo flow-through exposures to graded concentrations of water-borne B[a]P. Significantly greater amounts of B[a]P metabolites were observed with exposure to greater concentrations of parent compound in both life stages. However, when assessing biotransformation capacity, no differences in phase I or phase II biotransformation were observed with greater exposures to B[a]P. Results of modelling suggested that biotransformation of B[a]P can be successfully implemented into in silico models to accurately predict life stage-specific abundances of B[a]P metabolites in either whole-body larvae or the bile of adult fish. Models developed increase the scope of applications in which TK models can be used to support environmental risk assessments.

In vitro-in vivo extrapolation of hepatic and gastrointestinal biotransformation rates of hydrophobic chemicals in rainbow trout

Hepatic in vitro biotransformation assays, in combination with in vitro-in vivo extrapolation (IVIVE) and bioaccumulation modeling, can be used to support regulatory bioaccumulation assessments. In most applications, however, these methods ignore the possibility of extrahepatic metabolism. Here we evaluated intestinal biotransformation in rainbow trout using S9 fractions prepared from the upper intestinal (GIT) epithelium. Measured levels of activity determined using standard substrates for phase I and phase II biotransformation enzymes were within 2-fold of activities measured in hepatic S9 fractions. In vitro intrinsic clearance rates for 2-ethylhexyl-4-methoxycinnamate (EHMC; an organic sunscreen agent) and two polycyclic aromatic hydrocarbons (pyrene [PYR] and benzo(a)pyrene [BAP]) were significantly higher in liver S9 fractions than in GIT S9 fractions. For octocrylene (OCT; a second sunscreen agent), however, in vitro intrinsic clearance rates were higher in GIT S9 fractions compared to liver S9 fractions. An existing 'liver only' IVIVE model was expanded to consider biotransformation in both the liver and GIT. Relevant IVIVE scaling factors were developed by morphological, histological, and biochemical evaluation of trout intestines. For chemicals biotransformed at higher rates by hepatic S9 fractions (i.e., BAP, PYR, EHMC), the 'liver & GIT' model yielded whole-body biotransformation rate constants (k_MET) that were within 1.2 to 1.4-fold of those estimated using the 'liver only' model. In contrast to these findings, the mean k_MET for OCT obtained using the 'liver & GIT' model was 3.3 times higher than the mean k_MET derived using the 'liver only' model and was in good agreement with empirical k_MET estimates determined previously for trout (<20 % difference). The results of this study suggest that current 'liver only' IVIVE approaches may underestimate in vivo biotransformation rates for chemicals that undergo substantial biotransformation in the GIT.

Temperature effect on perfluorooctane sulfonate toxicokinetics in rainbow trout (Oncorhynchus mykiss): Exploration via a physiologically based toxicokinetic model

Salmonids are poikilotherms, which means that their internal temperature varies with that of water. Water temperature thus controls many of their lifecycle processes and physiological functions, which could influence the mechanisms of absorption, distribution, metabolism and excretion (ADME) of many substances, including perfluorinated alkyl acids (PFAAs). However, the processes governing the fate of PFAAs are still poorly understood in fish. Here we developed a physiologically-based toxicokinetic (PBTK) model for rainbow trout (Oncorhynchus mykiss) to study changes in physiological functions and PFAA ADME at different temperatures. The model was calibrated using experimental data from dietary exposure to perfluorooctane sulfonate at 7 °C and 19 °C. Predictions of PFOS concentrations were globally satisfactory at both temperatures, when accounting for the influence of temperature on growth, ventilation rate, cardiac output, clearances, and absorption rates. Accounting for the influence of temperature on tissue-plasma partition coefficients significantly improved predicted in-organ PFOS concentrations.

In vitro evaluation of the metabolic stability of nine fragrance chemicals in trout and human hepatocytes

In vitro metabolic stability of nine fragrance chemicals: p-tolyl acetate, cashmeran, ethylene brassylate, celestolide, galaxolide, traseolide, ambretone, tonalide and pentadecanolide, was evaluated in trout and human hepatocytes. The compounds were incubated with trout hepatocytes at 12°C and human hepatocytes at 37°C. Quantification of compound disappearance with time was performed using gas chromatography/mass spectrometry. in vivo hepatic intrinsic clearance values were calculated from the in vitro data. Significant metabolism was observed with trout hepatocytes for five of the nine fragrance chemicals, while all nine were metabolized significantly with human hepatocytes. Previously published models were used to examine expected bioaccumulation and persistence in whole organisms. Calculated half-lives due to metabolism of the nine chemicals are significantly shorter for humans than trout: <1 hour and <1 day, respectively. For all chemicals with demonstrated hepatic metabolism, the models indicate a lack of accumulation. For those where metabolism was demonstrated in trout, calculated bioconcentration factors would not be classified as bioaccumulative under prevailing regulatory systems.

Accumulation of Geosmin and 2-methylisoborneol in European Whitefish Coregonus Lavaretus and Rainbow Trout Oncorhynchus Mykiss in RAS

Geosmin (GSM) and 2-methylisoborneol (MIB)-induced off-flavors can cause serious problems in a recirculating aquaculture system (RAS), such as delayed harvest and increased production costs, but also damage producers’ reputation. Traditionally, off-flavors have been removed by depuration before harvesting. Rainbow trout (Oncorhynchus mykiss) and European whitefish (Coregonus lavaretus) are commercially valuable species produced for consumers, both being suitable for rearing in RAS. In this study, European whitefish and rainbow trout were raised from juvenile up to 240 g (European whitefish) and 660 g (rainbow trout) to monitor the long-term accumulation of off-flavors. The concentrations in fillet of rainbow trout reached 3.6 ng·g−1 (MIB) and 5.6 ng∙g−11 (GSM) with lipid content of 22.5%, while for European whitefish up to 3.2 ng·g−1 (MIB) and 3.9 ng·g−1 (GSM) were found with 14.8% in lipid content. Concentrations up to 58 ng·L−1 (MIB) and 49 ng·L−1 (GSM) were found in the circulating water. Based on the results, the accumulation of MIB proceeds at similar pace for both species. In the case of GSM, the accumulation started similarly for both species but proceeded more quickly for rainbow trout after 140 days of the experiment, with a statistically significant difference (p < 0.05).

A perspective on the role of fugacity and activity for evaluating the PBT properties of organic chemicals and providing a multi-media synoptic indicator of environmental contamination

The hazard and risk from organic chemicals present in the environment are routinely evaluated using P (persistence), B (bioaccumulation) and T (toxicity) criteria. We present a perspective on how models based on the equilibrium criteria of fugacity and chemical activity can contribute to all three evaluations, thus providing a consistent and seamless assessment process. Persistence and its closely related, but fundamentally different property, chemical residence time can be determined from degradation half-lives (typically obtained from monitoring data) and multi-media mass balance models describing chemical fate and transport in real or evaluative environments. Bioaccumulation is driven by equilibrium partitioning processes that can be estimated from fugacity models treating uptake by respiration and diet in single organisms and in food webs, most commonly for aquatic systems for which confirmatory monitoring data can be obtained. Biomagnification is readily evaluated both experimentally and as a fugacity (or activity) increase from prey to predator. Toxicokinetic evaluations of chemical uptake and disposition in a variety of organisms, including humans, can be obtained using fugacity- and physiologically-based pharmacokinetic models. Toxicodynamic processes of chemical interactions with organisms leading to adverse toxic outcomes are less obviously amenable to fugacity assessment with the notable exception of baseline toxicity (narcosis), the most common mode of action of commercial chemicals. It is shown that a range of fugacities can be defined and correlated for specific species thus enabling direct comparison of fugacities estimated or monitored in the environment with fugacities resulting in baseline toxicity. Finally, it is shown that fugacity and activity can serve as multi-media metrics of chemical contamination in ecosystems, thus enabling identification of specific regions and species at greatest risk in a pictorial format that is readily understandable to scientific and regulatory communities and to the general public.

Development of a generic zebrafish embryo PBPK model and application to the developmental toxicity assessment of valproic acid analogs

In order to better explain, predict, or extrapolate to humans the developmental toxicity effects of chemicals to zebrafish (Danio rerio) embryos, we developed a physiologically-based pharmacokinetic (PBPK) model designed to predict organ concentrations of neutral or ionizable chemicals, up to 120 h post-fertilization. Chemicals' distribution is modeled in the cells, lysosomes, and mitochondria of ten organs of the embryo. The model's partition coefficients are calculated with sub-models using physicochemical properties of the chemicals of interest. The model accounts for organ growth and changes in metabolic clearance with time. We compared ab initio model predictions to data obtained on culture medium and embryo concentrations of valproic acid (VPA) and nine analogs during continuous dosing under the OECD test guideline 236. We further improved the predictions by estimating metabolic clearance and partition coefficients from the data by Bayesian calibration. We also assessed the performance of the model at reproducing data published by Brox et al. (2016) on VPA and 16 other chemicals. We finally compared dose-response relationships calculated for mortality and malformations on the basis of predicted whole embryo concentrations versus those based on nominal water concentrations. The use of target organ concentrations substantially shifted the magnitude of dose-response parameters and the relative toxicity ranking of chemicals studied.

Multicompartmental Toxicokinetic Modeling of Discrete Dietary and Continuous Waterborne Uptake of Two Polycyclic Aromatic Hydrocarbons by Zebrafish Danio rerio

In the present study, we developed a multicompartmental toxicokinetic model for two polycyclic aromatic hydrocarbons (phenanthrene and anthracene) in their deuterated form (PAHs-d₁₀) in zebrafish considering continuous waterborne uptake and discrete dietary uptake. We quantified the bioconcentration, bioaccumulation, and depuration of these two PAHs-d₁₀ in zebrafish, and then estimated the kinetic parameters by fitting the model into the experimental data. The experimental and fitting results both showed that there was a peak concentration in each compartment of zebrafish after every dietary uptake, while the peak value depended on the ingestion amount of the PAH-d₁₀ and varied among different compartments. The PAH-d₁₀ amount in the blood reached 20-27% of the total amount bioaccumulated in zebrafish at steady-state, followed by skin (20-26%), and fillet (16-22%). The rank of PAH-d₁₀ steady-state concentrations in each compartment showed inconsistency with its lipid contents, indicating that the distribution of the PAHs-d₁₀ in zebrafish was not merely affected by the lipid content in each compartment, but also affected by their kinetics and biotransformation. This study suggests that discrete dietary uptake caused by intermittent food ingestion significantly affects the bioaccumulation of PAHs in fish. Further studies are needed to investigate such effect on other toxicants that are more resistant to biotransformation than PAHs in fish.

Variability in fish bioconcentration factors: Influences of study design and consequences for regulation

The fish bioconcentration factor (BCF) is an important aspect within bioaccumulation assessments. Several factors have been suggested to influence BCF values - including species, developmental stage, mixture exposure, and calculation method. However, their exact contribution to variance in BCF values is unknown. Within this study we assessed the relative impact of these test characteristics on BCF values and analyzed the reproducibility of aquatic exposure bioconcentration tests. Linear mixed effects analyses were performed on a newly develop database to investigate the relationship between the response variable (i.e. lipid normalized log BCF values) and several test characteristics as fixed effects. Lower BCF values were observed for substances that were simultaneously applied with high molecular weight polycyclic aromatic hydrocarbons compared to single substance exposure (with an average difference of -0.81 log BCF). Also, lower BCFs upon kinetic determination were observed compared to steady-state BCFs (log BCF -0.27), and lower BCFs for species from the Ostariophysi subcohort level (log BCF -0.17 to -0.15). In addition, data analysis showed high variation within BCF values for single substances (average SD = log BCF 0.21), which questions the robustness of the current bioaccumulation assessments. For example, the 95% confidence range of a BCF value of 2500 ranges from 953 ('not-bioaccumulative') to 6561 ('very bioaccumulative'). Our results show that the use of one single BCF leads to a high uncertainty in bioaccumulation assessments. We strongly recommend that within future bioconcentration studies, the used experimental design and test conditions are described in detail and justified to support solid interpretation.

Elucidating the fate of perfluorooctanoate sulfonate using a rainbow trout (Oncorhynchus mykiss) physiologically-based toxicokinetic model

Per- and poly-fluorinated substances (PFAS) are widely found in freshwater ecosystems because of their resistance to degradation. Among them, several long-chain perfluoroalkyl acids bioaccumulate in aquatic vertebrates, but our understanding of the mechanisms of absorption, distribution and elimination is still limited in fish. For this purpose, we developed a 10-compartment physiologically-based toxicokinetic (PBTK) model to elucidate perfluorooctane sulfonate (PFOS) kinetics in adult rainbow trout. This PBTK model included various physiological characteristics: blood perfusion to each organ, plasmatic fraction, PFOS free fraction, and growth of individuals. The parameters were optimized using Bayesian inferences. First, only PFOS absorption by diet was considered in the model as well as its elimination by urine, bile and feces. Then two mechanistic hypotheses, assumed to govern PFOS toxicokinetics in fish, namely the enterohepatic cycle and the absorption and elimination though gills, were tested. Improvement of the model's fit to the data was studied in each organ by comparing predictions with observed data using relative error. The experimental data set was obtained from an exposure experiment, where adult rainbow trout were fed with a PFOS-spiked diet for 42 days, followed by a 35-day depuration period. In all cases, PFOS concentrations were accurately predicted in organs and feces by the model. The results of this PBTK model demonstrated that feces represented the major elimination route for PFOS while urine was a minor route. Also, PFOS branchial uptake can be substantial despite low concentrations of the compound in water, and elimination through gills should not be neglected. Finally, the enterohepatic cycle is likely to play a minor role in PFOS toxicokinetics. Overall, this PBTK model accurately described PFOS distribution in rainbow trout and provides information on the relative contribution of absorption and elimination pathways.

A Petri Net Approach to Physiologically Based Toxicokinetic Modeling

Physiologically based toxicokinetic (PBTK) modeling enables researchers to predict internal tissue concentrations for various species exposed to exogenous compounds through different routes at varying concentrations without having to run in vivo experiments for each scenario. Parameters for the models may be gathered from in vivo or in vitro measurements, cross-species or cross-chemical extrapolations, literature reviews, or other models. The PBTK models, described using ordinary differential equations (ODEs), are then simulated using these parameters for a given compound/exposure/species scenario. Although they are potentially useful for regulatory toxicology, the complexity of ODE programming and simulation remains a barrier for many would-be researchers. Petri nets, a graphical modeling framework, offers a more intuitive approach to PBTK modeling. To demonstrate their utility and ease of use, we present a model of waterborne fluoranthene exposure to rainbow trout (Oncorhynchus mykiss) written and simulated in Snoopy, a graphical Petri net development and simulation software package. We converted an existing ODE PBTK model and evaluated the Petri net model against the ODE model results. The simulated tissue concentrations of the Petri net model closely mirrored the simulated concentrations of the ODE model. To convert the ODE model to a Petri net model, we introduced a new parameter, blood volume (V _BLOOD ). Sensitivity analysis found V _BLOOD to be very robust when varied over an order of magnitude. The resulting Petri net PBTK model has a number of advantages over ODE models, while maintaining equivalent predictive functionality. Environ Toxicol Chem 2019;00:1-10. © 2019 SETAC.

Physiologically based toxicokinetics (PBTK) models for pharmaceuticals and personal care products in wild common carp (Cyprinus carpio)

Pharmaceutical and personal care products (PPCPs) are regarded as "pseudo-persistent" pollutants due to their being continuously loaded into the aquatic environment. Physiologically based toxicokinetics (PBTK) models that can quantitatively describe absorption, distribution, metabolism and excretion processes of chemicals in biota are of importance to predict internal exposure (e.g. doses at specific target tissues/organs) from external exposure concentrations. In this study, PBTK models with up to six compartments including brain, liver, kidney, gills, richly perfused tissues and poorly perfused tissues, were developed for predicting internal distribution of 10 PPCPs in wild common carp (Cyprinus carpio). The PBTK predicted concentrations were close to the measured ones, with deviations less than 1 log unit for most of PPCPs. Sensitivity analysis showed that various partition coefficients of the chemicals exerted significant influence on model outputs.

Investigating the interaction between melamine and cyanuric acid using a Physiologically-Based Toxicokinetic model in rainbow trout

Following outbreaks of feed and food adulterations with a melamine and cyanuric acid mixture in 2007 and melamine in 2008 respectively, the kinetics and toxicodynamics of the mixture have been investigated particularly in sensitive species such as the rainbow trout. Tissue concentrations and intensity of the adverse effect, melamine-cyanurate crystal formation in kidney, were reported in similar experimental conditions. Here, a recent PBTK model for rainbow trout has been applied to model the kinetics of both single compounds based on residue levels in tissues. Both PBTK models for the single compounds were combined and a model of crystal formation for the mixture melamine-cyanuric acid was also added to predict the intensity of crystal formation under the assumptions that crystals formed either in urine or in kidney tissue. Modelling the kinetics of melamine and cyanuric acid provided a better understanding and prediction of intensity of crystal formation in case of sequential exposures with varying intensity or co-exposure. This study demonstrates, for the first time, how fish PBTK models can play a key role in the understanding and prediction of toxicokinetics and toxicodynamics of mixtures. This study also illustrates how adverse effects may potentially occur even when the compounds are not administered together as a mixture.

Evaluation of the bioaccumulation potential of selected alternative brominated flame retardants in marine fish using in vitro metabolic transformation rates

The global consumption of alternative brominated flame retardants (BFRs) has increased with the restriction of the first generation BFRs such as polybrominated diphenyl ethers (PBDEs) and hexabromocyclododecanes (HBCDs). However, many alternative BFRs are suspected to be persistent in the environment and possibly bioaccumulative after their release into the environment because of their chemical properties, which are similar to those of the already banned BFRs. In this study, the bioaccumulation potential of selected alternative BFRs (1,2-bis(2,4,6‑tribromophenoxy)ethane (BTBPE), 1,2,3,4,5,6‑hexabromobenzene (HBB), pentabromoethylbenzene (PBEB), 2,3,4,5,6‑pentabromotoluene (PBT), 2‑ethylhexyl-2,3,4,5-tetrabromobenzoate (TBB), and 2,3,4,5‑tetrabromo-6-chlorotoluene (TBCT)) was evaluated. The in vitro depletion rate constants (k_depl) were measured for the alternative BFRs using liver S9 fractions isolated from five marine fish species (Epinephelus septemfasciatus, Konosirus punctatus, Lateolabrax japonicus, Mugil cephalus, and Sebastes schlegelii) that inhabit the oceans off the Korean coast. The measured k_depl values were converted to in vitro intrinsic clearance rate constants (CL_{in vitro}) to estimate whole-body metabolic rate constants (k_MET) using an in vitro to in vivo extrapolation (IVIVE) model. Finally, the bioconcentration factors (BCF) were determined using a one-compartment model. The transformation kinetics for obtaining k_depl agreed well with first-order chemical kinetics, regardless of initial BFR concentrations. The values of CL_{in vitro} were lower in the selected marine fish species than those in freshwater fish species, implying slower metabolic transformation. The derived BCF values based on the total concentration in water (BCF_TOT) ranged from 16 (TBB in M. cephalus) to 27,000 (HBB in K. punctatus). The BCF values for HBB and PBT were >2000 except for those in M. cephalus suggesting further investigation of BCF values of BFRs whose log K_OW is between 6 and 7.

Generic physiologically-based toxicokinetic modelling for fish: Integration of environmental factors and species variability

One of the goals of environmental risk assessment is to protect the whole ecosystem from adverse effects resulting from exposure to chemicals. Many research efforts have aimed to improve the quantification of dose-response relationships through the integration of toxicokinetics. For this purpose, physiologically-based toxicokinetic (PBTK) models have been developed to estimate internal doses from external doses in a time-dependent manner. In this study, a generic PBTK model was developed and adapted for rainbow trout (Onchorhynchus mykiss), zebrafish (Danio rerio), fathead minnow (Pimephales promelas), and three-spined stickleback (Gasterosteus aculeatus). New mechanistic approaches were proposed for including the effects of growth and temperature in the model. Physiological parameters and their inter-individual variability were estimated based on the results of extensive literature searches or specific experimental data. The PBTK model was implemented for nine environmental contaminants (with log k_ow from -0.9 to 6.8) to predict whole-body concentrations and concentrations in various fish's organs. Sensitivity analyses were performed for a lipophilic and a hydrophilic compound to identify which parameters have most impact on the model's outputs. Model predictions were compared with experimental data according to dataset-specific exposure scenarios and were accurate: 50% of predictions were within a 3-fold factor for six out of nine chemicals and 75% of predictions were within a 3-fold factor for three of the most lipophilic compounds studied. Our model can be used to assess the influence of physiological and environmental factors on the toxicokinetics of chemicals and provide guidance for assessing the effect of those critical factors in environmental risk assessment.

An assessment of applicability of existing approaches to predicting the bioaccumulation of conventional substances in nanomaterials

The experimental determination of bioaccumulation is challenging, and a number of approaches have been developed for its prediction. It is important to assess the applicability of these predictive approaches to nanomaterials (NMs), which have been shown to bioaccumulate. The octanol/water partition coefficient (K_OW ) may not be applicable to some NMs that are not found in either the octanol or water phases but rather are found at the interface. Thus the K_OW values obtained for certain NMs are shown not to correlate well with the experimentally determined bioaccumulation. Implementation of quantitative structure-activity relationships (QSARs) for NMs is also challenging because the bioaccumulation of NMs depends on nano-specific properties such as shape, size, and surface area. Thus there is a need to develop new QSAR models based on these new nanodescriptors; current efforts appear to focus on digital processing of NM images as well as the conversion of surface chemistry parameters into adsorption indices. Water solubility can be used as a screening tool for the exclusion of NMs with short half-lives. Adaptation of fugacity/aquivalence models, which include physicochemical properties, may give some insights into the bioaccumulation potential of NMs, especially with the addition of a biota component. The use of kinetic models, including physiologically based pharmacokinetic models, appears to be the most suitable approach for predicting bioaccumulation of NMs. Furthermore, because bioaccumulation of NMs depends on a number of biotic and abiotic factors, it is important to take these factors into account when one is modeling bioaccumulation and interpreting bioaccumulation results. Environ Toxicol Chem 2018;37:2972-2988. © 2018 SETAC.

Measurement of kinetic parameters for biotransformation of polycyclic aromatic hydrocarbons by trout liver S9 fractions: Implications for bioaccumulation assessment

In vitro substrate depletion methods developed by the pharmaceutical industry are being used with increasing frequency to support chemical bioaccumulation assessments for fish. However, the application of these methods to high log K _ow chemicals poses special challenges. Biotransformation of three polycyclic aromatic hydrocarbons (PAHs) was measured using trout liver S9 fractions. Measured activity declined with incubation time and was reduced by acetone (used as a spiking solvent) at concentrations greater than 0.5%. Addition of alamethicin, a pore-forming peptide used to support UDP-glucuronosyltransferase activity, also reduced activity in a concentration-dependent manner. The substrate concentration dependence of activity was evaluated to estimate K _M and V _max values for each compound. Derived kinetic constants suggested that all three PAHs are transformed by the same reaction pathway and indicated an inverse correlation between K _M and chemical log K _ow. Binding effects on activity were evaluated by measuring unbound chemical concentrations across a range of S9 protein levels. Reaction rates were proportional to the unbound concentration except when these concentrations approached saturating levels, providing a direct demonstration of the free chemical hypothesis. These findings suggest that previous in vitro work with high log K _ow compounds was conducted at inappropriately high substrate concentrations resulting in underestimation of true in vivo activity. Preliminary calculations also indicate that PAH metabolism in fish may approach saturation during standardized in vivo testing efforts, potentially resulting in concentration-dependent accumulation and/or steady-state levels of accumulation greater than those which occur in a natural setting.

Uptake, distribution and elimination of chemicals in fish - Which physiological parameters are the most relevant for toxicokinetics?

Bioconcentration and toxicity studies are regularly conducted for the risk assessment of chemicals. If such tests yield different results for different fish species, this can either be due to differences in toxicokinetics or to differences in toxicodynamics. Here we investigate which physiological parameters could cause major differences in the toxicokinetics in fish. To this end it is important to distinguish physiological parameters that affect the sorption capacity of the fish from those that affect kinetic processes. Variability in the lipid content of a fish is the most influential parameter for the sorption capacity of fish and therefore most relevant for the total concentration in fish under steady-state conditions when metabolism is not relevant. In terms of kinetics, ventilation rate, uptake efficiency from food and metabolism are the most influential factors. While ventilation rate can roughly be estimated from allometric scaling equations, little general information is available on the uptake efficiency from food. The metabolism rate constant appears to be the single most influential toxicokinetic factor. This information cannot be estimated but must be determined experimentally, preferably from in vitro experiments.

Multi-compartmental toxicokinetic modeling of fipronil in tilapia: Accumulation, biotransformation and elimination

Bioaccumulation and biotransformation are critical processes modifying toxicity of easily metabolizable chemicals to aquatic organisms. In this study, tissue-specific accumulation, biotransformation and elimination of a current-use pesticide fipronil in tilapia (Oreochromis niloticus) were quantified by combining in vivo measurements and a newly developed multi-compartmental toxicokinetic model. Waterborne fipronil was taken up via gills and metabolized rapidly and solely to fipronil sulfone. Significant decrease of fipronil residues in liver and intestine during exposure period strongly suggested the induction of metabolism in these two organs. Significant transport of fipronil and fipronil sulfone in the liver-bile-intestine system implied that hepatobiliary excretion and enterohepatic re-absorption played important roles in fipronil metabolism and system circulation of the parent compound and the metabolite. The multi-compartmental model quantitatively described the highly dynamic inter-compartmental transport and rapid branchial clearance of fipronil in fish. Modeling results also suggested that uptake and biotransformation were the stronger driving forces for the inter-compartmental transport of fipronil in fish than the inherent partitioning capacity. Overall, our findings highlight the importance of biotransformation on internal disposition of fipronil in fish, which helps to improve aquatic toxicity assessment of this pesticide.

Screening tools for the bioconcentration potential of monovalent organic ions in fish

Currently the bioaccumulation potential of organic chemicals is assessed in a first tier approach via their octanol-water partition coefficient. This approach has been developed for neutral chemicals and cannot work for ionizable and ionic chemicals because the latter have different sorption-mechanisms and -preferences. Thus, suitable screening tools for the bioconcentration potential of ionic and ionizable chemicals need to be developed because it cannot be expected that these chemicals are non-bioaccumulative per se. Here, we present such screening tools for monovalent ions and ionizable chemicals based on calibrated sorption models for membrane lipids, structural proteins and albumin. The molecular descriptors used for these models arise from quantum chemical calculations and are based on COSMO-RS theory. When we applied our screening tools to 1839 preselected chemicals from the REACH registration data base, we identified 187 chemicals as potentially bioconcentrating (still ignoring any kind of metabolism). Among these were carbon and sulphur based aromatic and aliphatic acids mostly with a rather high molecular surface area. We hope that this outcome will trigger further research on ion specific sorption mechanisms and lead to a re-evaluation of the bioconcentration potential of ionic chemicals.

Equilibrium biopartitioning of organic anions - A case study for humans and fish

In this work we combine partition coefficients between water and membrane lipid, storage lipid, the plasma protein albumin as well as structural protein with the tissue dependent fraction of the respective phases in order to obtain a clearer picture on the relevance of various biological tissues for the bioaccumulation of 31 organic anions. Most of the partition coefficients are based on experimental data, supplemented by some predicted ones. The data suggest that the plasma protein, albumin, will be the major sorption matrix in mammals. Only small fractions of the studied chemicals will occur freely dissolved in an organism. For the investigated acids with pKa <5, partitioning is dominated by the ionic species rather than the corresponding neutral species. Bioconcentration in fish is not expected to occur for many of these acids unless pH in the aqueous environment is low or specific sorption mechanisms are relevant. In contrast, biomagnification in terrestrial mammals would be expected for most organic anions if they are not sufficiently metabolized. We conclude that sorption is important for the toxicokinetics of ionizable organic chemicals and the dominating sorbing matrices are quite different from those for neutral species.

Biotransformation of Benzo[ a]pyrene by Three Rainbow Trout ( Onchorhynchus mykiss) Cell Lines and Extrapolation To Derive a Fish Bioconcentration Factor

Permanent fish cell lines constitute a promising complement or substitute for fish in the environmental risk assessment of chemicals. We demonstrate the potential of a set of cell lines originating from rainbow trout ( Oncorhynchus mykiss) to aid in the prediction of chemical bioaccumulation in fish, using benzo[ a]pyrene (BaP) as a model chemical. We selected three cell lines from different tissues to more fully account for whole-body biotransformation in vivo: the RTL-W1 cell line, representing the liver as major site of biotransformation, and the RTgill-W1 (gill) and RTgutGC (intestine) cell lines, as important environment-organism interfaces, which likely influence chemical uptake. All three cell lines were found to effectively biotransform BaP. However, rates of in vitro clearance differed, with the RTL-W1 cell line being most efficient, followed by RTgutGC. Co-exposures with α-naphthoflavone as potent inhibitor of biotransformation, assessment of CYP1A catalytic activity, and the progression of cellular toxicity upon prolonged BaP exposure revealed that BaP is handled differently in the RTgill-W1 compared to the other two cell lines. Application of the cell-line-derived in vitro clearance rates into a physiology-based toxicokinetic model predicted a BaP bioconcentration factor (BCF) of 909-1057 compared to 920 reported for rainbow trout in vivo.

Toxicokinetics of the neonicotinoid insecticide imidacloprid in rainbow trout (Oncorhynchus mykiss)

Studies were conducted to determine the distribution and elimination of imidacloprid (IMI) in rainbow trout. Animals were injected with a low (47.6 μg/kg), medium (117.5 μg/kg) or high (232.7 μg/kg) dose directly into the bloodstream and allowed to depurate. The fish were then sampled to characterize the loss of IMI from plasma and its appearance in expired water (all dose groups) and urine (medium dose only). In vitro biotransformation of IMI was evaluated using trout liver S9 fractions. Mean total clearance (CL_T) values determined by non-compartmental analysis of plasma time-course data were 21.8, 27.0 and 19.5 mL/h/kg for the low, medium and high dose groups, respectively. Estimated half-lives for the same groups were 67.0, 68.4 and 68.1 h, while fitted values for the steady-state volume of distribution (V_SS) were 1.72, 2.23 and 1.81 L/kg. Branchial elimination rates were much lower than expected, suggesting that IMI is highly bound in blood. Renal clearance rates were greater than measured rates of branchial clearance (60% of CL_T in the medium dose group), possibly indicating a role for renal membrane transporters. There was no evidence for hepatic biotransformation of IMI. Collectively, these findings suggest that IMI would accumulate in trout in continuous waterborne exposures.

The influence of chemical degradation during dietary exposures to fish on biomagnification factors and bioaccumulation factors

The chemical dietary absorption efficiency (E_D) quantifies the amount of chemical absorbed by an organism relative to the amount of chemical an organism is exposed to following ingestion. In particular, E_D can influence the extent of bioaccumulation and biomagnification for hydrophobic chemicals. A new E_D model is developed to quantify chemical process rates in the gastrointestinal tract (GIT). The new model is calibrated with critically evaluated measured E_D values (n = 250) for 80 hydrophobic persistent chemicals. The new E_D model is subsequently used to estimate chemical reaction rate constants (k_R) assumed to occur in the lumen of the GIT from experimental dietary exposure tests (n = 255) for 165 chemicals. The new k_R estimates are corroborated with k_R estimates for the same chemicals from the same data derived previously by other methods. The roles of k_R and the biotransformation rate constant (k_B) on biomagnification factors (BMFs) determined under laboratory test conditions and on BMFs and bioaccumulation factors (BAFs) in the environment are examined with the new model. In this regard, differences in lab and field BMFs are highlighted. Recommendations to address uncertainty in E_D and k_R data are provided.

Evaluating parameter availability for physiologically based pharmacokinetic (PBPK) modeling of perfluorooctanoic acid (PFOA) in zebrafish

Physiologically based pharmacokinetic (PBPK) models are considered useful tools to describe the absorption, distribution, metabolism and excretion of xenobiotics. For accurate predictions, PBPK models require species-specific and compound-specific parameters. Zebrafish are considered an appropriate vertebrate model for investigating the toxicity of a wide variety of compounds. However, no specific mechanistic model exists for the pharmacokinetics of perfluoroalkyl acids (PFAAs) in zebrafish, despite growing concern about this class of ubiquitous environmental contaminants. The purpose of this study was to evaluate the current state of knowledge for the parameters that would be needed to construct such a model for zebrafish. We chose perfluorooctanoic acid (PFOA) as a model PFAA with greater data availability. We have updated a previous PBPK model for rainbow trout to simulate PFOA fate in zebrafish following waterborne exposure. For the first time, the model considers hepatobiliary circulation. In order to evaluate the availability of parameters to implement this model, we performed an extensive literature review to find zebrafish-specific parameters. As in previous approaches, we broadened our search to include mammalian and other fish studies when zebrafish-specific data were lacking. Based on the method used to measure or estimate parameters, or based on their species-specific origin, we scored and ranked the quality of available parameters. These scores were then used in Monte Carlo and partial rank correlation analyses to identify the most critical data gaps. The liver, where fatty acid binding proteins (FABPs) and plasma proteins are considered, represented the best model-data agreement. Lack of agreement in other tissues suggest better parameters are needed. The results of our study highlight the lack of zebrafish-specific parameters. Based on sensitivity and uncertainty analysis, parameters associated with PFAA-protein interactions and passive diffusion need further refinement to enable development of predictive models for these emerging chemicals in zebrafish.

Can poly-parameter linear-free energy relationships (pp-LFERs) improve modelling bioaccumulation in fish?

A wide range of studies have characterized different types of biosorbent, with regard to their interactions with chemicals. This has resulted in the development of poly-parameter linear free energy relationships (pp-LFERs) for the estimation of partitioning of neutral organic compounds to biological phases (e.g., storage lipids, phospholipids and serum albumins). The aims of this study were to explore and evaluate the influence of implementing pp-LFERs both into a one-compartment fish model and a multi-compartment physiologically based toxicokinetic (PBTK) fish model and the associated implications for chemical risk assessment. For this purpose, fish was used as reference biota, due to their important role in aquatic food chains and dietary exposure to humans. The bioconcentration factor (BCF) was utilized as the evaluation metric. Overall, our results indicated that models incorporating pp-LFERs (R² = 0.75) slightly outperformed the single parameter (sp) LFERs approach in the one-compartmental fish model (R² = 0.72). A pronounced enhancement was achieved for compounds with log K_OW between 4 and 5 with increased R² from 0.52 to 0.71. The minimal improvement was caused by the overestimation of lipid contribution and underestimation of protein contribution by the sp-approach, which cancelled each other out. Meanwhile, a greater improvement was observed for multi-compartmental PBTK models with consideration of metabolism, making all predictions fall within a factor of 10 compared with measured data. For screening purposes, the K_OW-based (sp-LFERs) approach should be sufficient to quantify the main partitioning characteristics. Further developments are required for the consideration of ionization and more accurate quantification of biotransformation in biota.

Enhanced uptake of BPA in the presence of nanoplastics can lead to neurotoxic effects in adult zebrafish

Plastic particles have been proven to be abundant in the aquatic environment, raising concerns about their potential toxic effects. In the present study, we determined the bioaccumulation potential of bisphenol A (BPA) in adult zebrafish (Danio rerio) in the absence and presence of nano-sized plastic particles (nanoplastics, NPPs). Results show that BPA can accumulate in the viscera, gill, head and muscle of zebrafish with 85, 43, 20, and 3μg/g ww after 1d exposure. NPPs were also found to accumulate in different tissues of the fish. Relative equilibrium was reached after 1d exposure in different tissues with 39 to 636mg/kg ww. Co-exposure of NPPs and BPA led to a 2.2 and 2.6-fold significant increment of BPA uptake in the head and viscera, if compared with BPA alone treatment after 3d exposure. As such, we further investigated several neurotoxic biomarker alterations in the fish head. It was found that either BPA or NPPs can cause myelin basic protein (MBP)/gene up-regulation in the central nervous system (CNS); meanwhile, both contaminants exhibited significant inhibition of acetylcholinesterase (AChE) activity, which is a well-known representative biomarker for neurotoxicity. Moreover, for the co-exposure treatment, biomarkers of myeline and tubulin protein/gene expressions, dopamine content, and the mRNA expression of mesencephalic astrocyte derived neurotrophic factor (MANF) were all significantly up-regulated, suggesting that an enhanced neurotoxic effects in both CNS and dopaminergic system occurred. However, AChE activity was no more inhibited in the co-exposure treatment, which implies that solely AChE measurement may not be sufficient to identify neurotoxic effects in the cholinergic system. Overall, the present study demonstrates that the presence of NPPs can increase BPA bioavailability and cause neurotoxicity in adult zebrafish.