A mixture parameterized biologically based dosimetry model to predict body burdens of polycyclic aromatic hydrocarbons in developmental zebrafish toxicity assays.

Polycyclic aromatic hydrocarbons (PAHs) are a group of environmental toxicants found ubiquitously as complex mixtures in human-impacted environments. Developmental zebrafish exposures have been used widely to study PAH toxicity, but most studies report nominal exposure concentrations. Nominal exposure concentrations can be unreliable dose metrics due to differences in toxicant bioavailability resulting from disparate exposure methodologies and chemical properties. Toxicokinetic modeling can predict toxicant tissue doses to facilitate comparison between exposures of different chemicals, methodologies, and biological models. We parameterize a biologically based dosimetry model for developmental zebrafish toxicity assays for 9 PAHs. The model was optimized with measurements from media, tissue, and plastic plate walls throughout a static developmental exposure to a mixture of 10 PAHs of high abundance within the Portland Harbor Superfund Site. Plate binding, volatilization, zebrafish permeability, and tissue-media partitioning coefficients vary widely between PAHs. Model predictions accounted for 83% and 54% of 48 hpf body burdens within a factor of 2 resulting from exposures to mixtures and individual PAHs, respectively. Accounting for solubility significantly improves model performance. Competition for active sites in metabolizing enzymes may change biotransformation kinetics between individual PAH and mixture exposures. Area under the curve estimations of concentrations in zebrafish resulted in altered hazard rankings from nominal exposure concentrations. Future work will be oriented to generalizing the model to other PAHs. This PAH dosimetry model improves the interpretability of developmental zebrafish toxicity assays by providing time-resolved body burdens from nominal exposure concentrations.