Studying interaction effects on toxicokinetics in zebrafish combining experimental and modelling approaches

Humans and wildlife are exposed to a complex mixture of anthropogenic chemicals of which only a few have been subjected to regulations. Chemical risk assessment is currently based on evaluating single chemicals, which is costly, time-consuming, and neglect toxicokinetic and toxicodynamic mixture effects. This study focused on interaction effects on the absorption, distribution, metabolism and excretion (ADME) processes of selected chemicals representing potential modulators of these processes. Adult female zebrafish (Danio rerio) were exposed to selected mixture of 11 chemicals and bioconcentration factors (BCFs) on tissue level were determined for 9 of them: bisphenol A (BPA), bisphenol AF (BPAF), bisphenol Z (BPZ), triclosan, tribromophenol, pentachlorophenol, heptafluorobutyric acid (PFBA), perfluorobutanesulfonic acid (PFBS), and perfluorooctanesulfonic acid (PFOS). Comparison of BCFs of bisphenols obtained from single chemical exposure experiments versus the current study revealed no statistically significant differences (p > 0.05), implying no mixture effects on kinetics of bisphenols at investigated concentrations. The same conclusion was reached using two physiologically based kinetic (PBK) models, developed for individual bisphenols and per- and polyfluoroalkyl substances (PFAS), showing good model fit for BPA, BPZ, BPAF, and PFOS. To simulate exposure scenarios where kinetic interaction effects may occur through competitive protein binding in blood, a new PBK model was developed. Simulations where zebrafish were dosed with BPA and BPZ, individually, and combined with varying levels of PFOS, showed that competitive binding to serum proteins alter tissue levels of bisphenols when levels of PFOS exceeded 1 μg/L. This indicates that chemicals acting in concert could perturb ADME but only at higher levels or in complex mixtures.