Dynamical coupling of PBPK/PD and AUC-based toxicity models for arsenic in tilapia Oreochromis mossambicus from blackfoot disease area in Taiwan

A physiologically based pharmacokinetic and pharmacodynamic (PBPK/PD) models were developed for arsenic (As) in tilapia Oreochromis mossambicus from blackfoot disease area in Taiwan. The PBPK/PD model structure consisted of muscle, gill, gut wall, alimentary canal, and liver, which were interconnected by blood circulation. We integrate the target organ concentrations and dynamic response describing uptake, metabolism, and disposition of As and the associated area-under-curve (AUC)-based toxicological dynamics following an acute exposure. The model validations were compared against the field observations from real tilapia farms and previously published uptake/depuration experimental data, indicating that predicted and measured As concentrations in major organs of tilapia were in good agreement. The model was utilized to reasonably simulate and construct a dose-dependent dynamic response between mortality effect and equilibrium target organ concentrations. Model simulations suggest that tilapia gills may serve as a surrogate sensitive biomarker of short-term exposure to As. This integrated As PBPK/PD/AUC model quantitatively estimates target organ concentration and dynamic response in tilapia and is a strong framework for future waterborne metal model development and for refining a biologically-based risk assessment for exposure of aquatic species to waterborne metals under a variety of scenarios.