A Physiologically Based Pharmacokinetic (PBPK) Modeling Framework for Mixtures of Dioxin-like Compounds

Humans are exposed to persistent organic pollutants, such as dioxin-like compounds (DLCs), as mixtures. Understanding and predicting the toxicokinetics and thus internal burden of major constituents of a DLC mixture is important for assessing their contributions to health risks. PBPK models, including dioxin models, traditionally focus on one or a small number of compounds; developing new or extending existing models for mixtures often requires tedious, error-prone coding work. This lack of efficiency to scale up for multi-compound exposures is a major technical barrier toward large-scale mixture PBPK simulations. Congeners in the DLC family, including 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), share similar albeit quantitatively different toxicokinetic and toxicodynamic properties. Taking advantage of these similarities, here we reported the development of a human PBPK modeling framework for DLC mixtures that can flexibly accommodate an arbitrary number of congeners. Adapted from existing TCDD models, our mixture model contains the blood and three diffusion-limited compartments—liver, fat, and rest of the body. Depending on the number of congeners in a mixture, varying-length vectors of ordinary differential equations (ODEs) are automatically generated to track the tissue concentrations of the congeners. Shared ODEs are used to account for common variables, including the aryl hydrocarbon receptor (AHR) and CYP1A2, to which the congeners compete for binding. Binary and multi-congener mixture simulations showed that the AHR-mediated cross-induction of CYP1A2 accelerates the sequestration and metabolism of DLC congeners, resulting in consistently lower tissue burdens than in single exposure, except for the liver. Using dietary intake data to simulate lifetime exposures to DLC mixtures, the model demonstrated that the relative contributions of individual congeners to blood or tissue toxic equivalency (TEQ) values are markedly different than those to intake TEQ. In summary, we developed a mixture PBPK modeling framework for DLCs that may be utilized upon further improvement as a quantitative tool to estimate tissue dosimetry and health risks of DLC mixtures.