A pharmacokinetic analysis of interspecies extrapolation in dioxin risk assessment

Methods for assessing the bioaccumulation of hydrocarbons and related substances in terrestrial organisms: A critical review

This study investigates and reviews methods for the assessment of the terrestrial bioaccumulation potential of hydrocarbons and related organic substances. The study concludes that the unitless biomagnification factor (BMF) and/or the trophic magnification factor (TMF) are appropriate, practical, and thermodynamically meaningful metrics for identifying bioaccumulative substances in terrestrial food chains. The study shows that various methods, including physical-chemical properties like the KOA and KOW , in vitro biotransformation assays, quantitative structure-activity relationships, in vivo pharmacokinetic and dietary bioaccumulation tests, and field-based trophic magnification studies, can inform on whether a substance has the potential to biomagnify in a terrestrial food chain as defined by a unitless BMF exceeding 1. The study further illustrates how these methods can be arranged in a four-tier evaluation scheme for the purpose of screening assessments that aim to minimize effort and costs and expediate bioaccumulation assessment of the vast numbers of organic substances in commerce, identifies knowledge gaps, and provides recommendations for further research to improve bioaccumulation assessment. Integr Environ Assess Manag 2023;19:1433-1456. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Physiologically based pharmacokinetic toolkit to evaluate environmental exposures: Applications of the dioxin model to study real life exposures

Chlorinated dibenzo-p-dioxins (CDDs) are a series of mono- to octa-chlorinated homologous chemicals commonly referred to as polychlorinated dioxins. One of the most potent, well-known, and persistent member of this family is 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). As part of translational research to make computerized models accessible to health risk assessors, we present a Berkeley Madonna recoded version of the human physiologically based pharmacokinetic (PBPK) model used by the U.S. Environmental Protection Agency (EPA) in the recent dioxin assessment. This model incorporates CYP1A2 induction, which is an important metabolic vector that drives dioxin distribution in the human body, and it uses a variable elimination half-life that is body burden dependent. To evaluate the model accuracy, the recoded model predictions were compared with those of the original published model. The simulations performed with the recoded model matched well with those of the original model. The recoded model was then applied to available data sets of real life exposure studies. The recoded model can describe acute and chronic exposures and can be useful for interpreting human biomonitoring data as part of an overall dioxin and/or dioxin-like compounds risk assessment.

An assessment of dioxin exposure across gestation and lactation using a PBPK model and new data from Seveso

On July 10, 1976, an explosion at a chemical plant in Seveso, Italy, released up to 30kg of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-the most potent dioxin congener. Twenty years later, the Seveso Women's Health Study (SWHS) initiated a follow-up assessment of a cohort of female Seveso residents. Researchers collected serial blood, measured for TCDD levels, and recorded information about the women's medical history after the explosion. The study's aims were to: 1) modify the human PBPK model for TCDD (Emond et al. 2004; Emond et al. 2005; NCEA-USEPA, 2010) to include repetitive gestation and lactation; 2) simulate TCDD blood concentrations during different life stages including pregnancy and lactation, under different exposure scenarios; and 3) use this PBPK model to compare the influence of gestation and lactation on elimination of TCDD. After optimization of the model, it was assessed using data from the SWHS cohort. The 23 women in Subcohort A, were 4-39years old and in Subcohort B, the 18 women were 3-17years old when the explosion occurred. The model accurately predicted the blood concentrations during the 20years post-exposure, including periods of pregnancy and lactation. The model was also used to analyze the contribution of gestation and lactation to the mother's elimination of TCDD. The results suggest that gestation and lactation do not significantly impact TCDD blood elimination. Future efforts will focus on using additional data to evaluate the PBPK model and improving the mathematical descriptions of lactation and multiple gestations.

A human PBPK/PD model to assess arsenic exposure risk through farmed tilapia consumption

The purpose of this study was to develop a biologically based risk assessment model for human health through consumption of arsenic (As) contaminated farmed tilapia (Oreochromis mossambicus) from blackfoot disease (BFD)-endemic area in Taiwan for estimating the consumption advice. We linked a physiologically based pharmacokinetic (PBPK) and a pharmacodynamic (PD) model to account for the exposure and dose-response profiles of As in human. Risk analysis indicates that consumption of farmed tilapia poses no significant threat from As-induced lung and bladder cancers. The predicted risk-based median consumption advice was no more than 5-17 meals month(-1) (or 2-6 g day(-1)).

Risk characterization and exposure assessment in arseniasis-endemic areas of Taiwan

This paper examines the age-specific human health risks exposed to inorganic arsenic through arsenic-contaminated farmed fish/shrimp and groundwater consumptions in arseniasis-endemic areas of blackfoot disease (BFD)-endemic area and Lanyang Plain in Taiwan, based on an probabilistic integrated risk assessment framework. We employ an age-dependent predictive physiologically-based pharmacokinetic model to account for arsenic concentrations in target organs. We reconstruct age-specific dose-response profiles for arsenicosis and arsenic-induced cancers by best fitting a pharmacodynamics-based three-parameter Hill equation model to published epidemiological data from West Bengal and Taiwan. The predicted median arsenic concentrations in age group-specific skin, lung, and bladder ranged from 2.24-5.70, 3.76-9.46, and 5.11-20.71 micro g g(-1) in BFD-endemic area, whereas 4.98-12.04, 8.23-19.92, and 11.07-43.45 micro g g(-1) in Lanyang Plain, respectively. Risk analysis indicates that consumption of arsenic-contaminated farmed fish/shrimp and groundwater in arseniasis-endemic areas may increase threat to prevalence of arsenicosis for all age groups, whereas adults may undergo potential risks of arsenic-induced skin, lung and bladder cancers. We show that peoples in Lanyang Plain are more readily associated with higher morbidities for arsenicosis and skin cancer as well as fatalities for lung and bladder cancers than that of peoples in BFD-endemic area. Here we report the first case in which theoretical human health risks for consuming As-contaminated farmed fish/shrimp and groundwater in the arseniasis-endemic areas are alarming under a conservative condition based on a probabilistic risk assessment framework.

Estimated cancer risk of dioxins to humans using a bioassay and physiologically based pharmacokinetic model

The health risk of dioxins and dioxin-like compounds to humans was analyzed quantitatively using experimental data and mathematical models. To quantify the toxicity of a mixture of three dioxin congeners, we calculated the new relative potencies (REPs) for 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 1,2,3,7,8-pentachlorodibenzo-p-dioxin (PeCDD), and 2,3,4,7,8- pentachlorodibenzofuran (PeCDF), focusing on their tumor promotion activity. We applied a liver foci formation assay to female SD rats after repeated oral administration of dioxins. The REP of dioxin for a rat was determined using dioxin concentration and the number of the foci in rat liver. A physiologically based pharmacokinetic model (PBPK model) was used for interspecies extrapolation targeting on dioxin concentration in liver. Toxic dose for human was determined by back-estimation with a human PBPK model, assuming that the same concentration in the target tissue may cause the same level of effect in rats and humans, and the REP for human was determined by the toxic dose obtained. The calculated REPs for TCDD, PeCDD, and PeCDF were 1.0, 0.34, and 0.05 for rats, respectively, and the REPs for humans were almost the same as those for rats. These values were different from the toxic equivalency factors (TEFs) presented previously (Van den Berg, M., Birnbaum, L., Bosveld, A.T.C., Brunstrom, B., Cook, P., Feeley, M., Giesy, J.P., Hanberg, A., Hasegawa, R., Kennedy, S.W., Kubiak, T., Larsen, J.C., Rolaf van Leeuwen, F.X., Liem, A.K.D., Nolt, C., Peterson, R.E., Poellinger. L., Safe, S., Schrenk, D., Tillitt, D, Tysklind, M., Younes, M., Waern, F., Zacharewski, T., 1998. Toxic equivalency factors (TEFs) for PCBs, PCDDs, PCDFs for humans and wildlife. Environ. Health Perspect. 106, 775-792). The relative risk of excess liver cancer for Japanese people in general was 1.7-6.5 x 10(-7) by TCDD only, and 2.9-11 x 10(-7) by the three dioxins at the present level of contamination.

Magnitude and mechanistic determinants of the interspecies toxicokinetic uncertainty factor for organic chemicals

The interspecies uncertainty factor, UF(AH), of 10 was recently subdivided into two components to account separately for interspecies differences in toxicokinetics and toxicodynamics (UF(AH-TK)=3.16, UF(AH-TD)=3.16). Even though the UF(AH) in its composite or dissociated form is used during the health risk assessment of systemic toxicants, there is no convincing scientific basis to justify the use of the same UF for all chemicals. In this study, we use equations that describe the toxicokinetics of chemicals at steady-state to characterize the magnitude and mechanistic determinants of UF(AH-TK) for several volatile organic chemicals (VOCs). Further, algorithms have been developed to determine the magnitude of the components of UF(AH-TK), namely the UF(AH-TK-ABS) (accounting for interspecies differences in dose absorbed during identical inhalation exposure conditions), UF(AH-TK-MET) (referring to the factor by which the blood concentration of unchanged parent chemical differs from one species to another, due to metabolic clearance, when both species receive identical doses) and UF(AH-TK-DIS) (reflecting the magnitude of difference in chemical concentrations distributed in target tissues of two species when the arterial blood concentration in both species is identical). The results show that the body weight, the rate of alveolar ventilation, the fraction of cardiac output flowing to the liver, partition coefficients (blood:air and tissue:blood), and the hepatic extraction ratio are the only parameters that play a critical role in the extrapolation of tissue and blood concentrations across species. Further, the magnitude of the UF(AH-TK-ABS) (means+/-SD, 0.19+/-0.04), UF(AH-TK-MET) (means+/-SD, 0.24+/-0.05) and UF(AH-TK-DIS) (mean range: 1.76-0.93) obtained in this study for several VOCs compares well with that obtained previously using physiologically based toxicokinetic models.

Dioxin health risk to infants using simulated tissue concentrations

Dioxin concentrations in infant and child were simulated using physiologically based pharmacokinetic (PBPK) models developed for these groups. The infant model was validated by comparing the simulated concentration with the measured concentration from the literature, and they showed good agreement. Simulations with our PBPK model showed temporal patterns in concentrations in various tissues. For risk assessment, estimated concentrations of 29 dioxins in the liver were summed up in a toxic equivalency (TEQ) basis to be compared with actual 2,3,7,8-TCDD concentrations in rat liver associated with toxicity. Maximum liver concentrations in breast-fed and formula-fed infants were 16.8pg TEQ/g and 3.5pg TEQ/g, respectively. The level in breast-fed infant liver was approximately 1/300 of the level associated with hepatocellular carcinoma and 1/5 of the level found in maternal rat liver associated with alterations in reproductive organs in the next generation. Based on our analysis, the present contamination level is not safe enough, but further dose-response data is required for a quantitative risk assessment.

Simulation of dioxin accumulation in human tissues and analysis of reproductive risk

Polychlorinated dibenzo-p-dioxin, polychlorinated dibenzofuran and dioxin-like polychlorinated biphenyl concentrations in human liver, kidney, fat, blood, muscle, richly perfused tissue (brain, lung etc.) and skin were simulated to assess the health risk for Japanese fetuses. A 40-year time course of dioxin accumulation via food ingestion was simulated using a physiologically based pharmacokinetic (PBPK) model. In richly perfused tissue, the concentration estimated by the PBPK model showed better agreement with measured concentrations than that calculated by the one-compartment model. Fetal dioxin concentration was simulated based on the assumption that the fetal concentration was almost equal to the concentration in the mother's richly perfused tissue. To assess the reproductive risk, the estimated concentration in human fetus was compared with that in rat fetus in which reproductive function showed signs of alteration by 2,3,7,8-TCDD in previous reports [Toxicol. Appl. Pharmacol. 114 (1992) 118; 146 (1997) 11; Toxicol. Sci. 53 (2000) 411; 57 (2000) 275]. The present daily intake of 2,3,7,8-TCDD is approximately 1/50 of the amount that leads to possible reproductive toxicity in the next generation. However, when 29 kinds of dioxin congeners are considered, the present level is 1/5 of the hazardous levels. For species extrapolation of dioxin risk, further study on tissue concentration versus toxicity is required.

Issues related to the use of canines in toxicologic pathology--issues with pharmacokinetics and metabolism

The dog is a commonly used animal model by virtue of its size, well-characterized physiology, and ease of handling. For these reasons and others, dogs are also useful in pharmacokinetic and metabolism studies during the development of both human and veterinary pharmaceutical products. In comparison with humans, or with other animals, dogs have some unique physiologic attributes that can affect the disposition of drugs. Species differences in gastrointestinal physiology, metabolism, renal function, and protein binding can affect the correlation of the pharmacokinetics and toxicology of dogs with those of other species. With the use of relevant examples, this article will provide an introduction to characteristics of dog physiology and their impact on pharmacokinetics, metabolism, drug disposition, toxicity, and dose selection.

Possible range of dioxin concentration in human tissues: simulation with a physiologically based model

In risk evaluation of dioxins, monitoring chemical concentrations in human tissues is an important step, and these concentration data can be utilized along with animal toxicity data for extrapolation of human manifestation. However large differences in dioxin concentrations usually exist even among individuals who have never been accidentally exposed to high quantities of dioxin, and this may cause problems in risk analysis. Body size, age, and history of food consumption are factors responsible for these interindividual differences in addition to exposure levels. Using a physiologically based pharmacokinetic (PBPK) model, the influence of differences in body weight, gastrointestinal absorption, and half-life and intake of dioxin were examined on tissue chemical concentration. Dioxin concentrations over a 40-yr time course in human liver, kidneys, fat, blood, muscle and richly perfused tissue were simulated for polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), and coplanar polychlorinated biphenyls (CoPCBs). Model parameters such as tissue-blood partition coefficients for CoPCBs were prepared, and sensitivity analysis was also performed on these parameters. The range of tissue concentrations was approximately 0.17 to 4.1 times the standard concentration, which was calculated using standard model parameters. The simulated ranges included more than 80% of the individual anatomical data for 2,3,7,8-tetrachlorodibenzo-p-dioxin, 1,2,3,7,-pentachlorodibenzo-p-dioxin, and 3,3',4,4',5-pentachlorobiphenyl in liver, fat, and blood. These results suggest that differences in body weight, gastrointestinal absorption, and food intake behavior may partially explain variation in tissue concentrations among individuals, and the possible interindividual uncertainty, which is approximately 24 for the general Japanese population.

Determination of tissue-blood partition coefficients for a physiological model for humans, and estimation of dioxin concentration in tissues

The tissue-blood partition coefficients for a physiologically based pharmacokinetic (PBPK) model were determined, and the concentrations of 17 congeners of polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans (PCDD/Fs) in tissues in Japanese people were estimated using the model. According to the PBPK model established by Lawrence and Gobas [Chemosphere 35 (1997) 427-452], we assumed a steady-state fugacity model for Japanese people in general, and set the route of PCDD/Fs exposure only from food intake. The required partition coefficients for liver, kidney, adipose, muscle, skin, bile, gut and viscera (richly perfused tissue) were calculated using available autopsy data from eight Japanese men and women who were not accidentally exposed to PCDD/Fs. For validation of the partition coefficients, estimated PCDD/F concentrations in liver, kidney, fat, blood and muscle using the model were compared to other two sets of measured concentration data in Japanese tissues. Good agreement was obtained between estimated data and measured data, and most of the measured data were within the simulated concentration range in liver, kidney, blood and muscle. From these results, our model and calculated partition coefficients seem applicable for the estimation of congener-specific concentrations in human tissues.

Disposition of polychlorinated dibenzo-p-dioxins, dibenzofurans, and non-ortho polychlorinated biphenyls in pregnant long evans rats and the transfer to offspring

Pharmacokinetic properties of polychlorinated dibenzo-p-dioxins, dibenzofurans (PCDFs), and non-ortho biphenyls (PCBs) play a critical role in their relative toxicity. The present study examined the transfer of these chemicals to offspring and placenta. Pregnant Long Evans rats received 0.0 (control), 0.05, 0.2, 0.8, and 1.0 microg/kg of dioxin toxic equivalence (TEQ) by oral gavage on the 15th gestational day (GD 15), using a dosing mixture that contained 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), 2,3,7,8-tetrachlorodibenzofuran (TCDF), 1,2,3,7,8-pentachlorodibenzo-p-dioxin (PeCDD), 1,2,3,7,8-pentachlorodibenzofuran (1-PeCDF), 2,3,4,7,8-pentachlorodibenzofuran (4-PeCDF), octachlorodibenzofuran (OCDF), 3,3',4,4'-tetrachlorobiphenyl (PCB 77), 3,3',4,4',5-pentachlorobiphenyl (PCB 126), and 3,3',4,4',5,5'-hexachlorobiphenyl (PCB 169) in ratios approximating that in food. Rats were euthanized on GD 16, GD 21, and postnatal day 4 (PND 4). The chemical concentrations in fetus, pup, placenta, and maternal liver, serum, and adipose tissue were determined using gas chromatography/high-resolution mass spectrometry. A dose-dependent increase in hepatic sequestration was seen with TCDD, PeCDD, 4-PeCDF, OCDF, PCB 126, and PCB 169, and the transfer to offspring was reduced at higher doses. 4-PeCDF, PeCDD and PCB 126 showed higher liver affinity than TCDD. TCDF, 1-PeCDF, and PCB 77 were metabolized rapidly. On GD 16, TCDD and the three PCBs reached equilibration between the fetus and placenta, but this did not occur with PeCDD and 4-PeCDF until GD 21, according to the lipid-based concentrations. Offspring compartments received more of the dosed compounds lactationally than transplacentally (7-28% versus 0.5-3%). The behavior of each congener was dose-dependent; therefore, extrapolation of high-dose experimental data should be used with caution.

A generic model of human lifetime exposure to persistent organic contaminants: development and application to PCB-101

We have developed a model which successfully reconstructs the lifetime polychlorinated biphenyl (PCB)-101 burden of the UK population for individuals born between 1920 and 1980. It not only follows burdens and clearance of persistent organic contaminants throughout a human lifetime--taking changes in age and body composition into account--but also, importantly, incorporates changing environmental concentrations of the compound of interest. Predicted results agree well with available measured lipid concentrations in human tissues. Its unique construction takes into account both changing environmental concentrations of PCBs in principal food groups and changing dietary habits during the time period. Because environmental burdens of persistent organic contaminants have changed over the last 60 years, residues in food will also have mirrored this change. Critically in this respect, the year in which an individual was born determines the shape and magnitude of their exposure profile for a given compound. Observed trends with age represent an historical legacy of exposure and are not simply a function of equal yearly cumulative inputs. We can demonstrate that the release profile of PCB-101 controls levels in the food supply and ultimately the burden of individuals throughout their life. This effect is expected to be similar for other PCB congeners and persistent organic compounds such as polychlorinated dibenzo-p-dioxins/furans (PCDD/Fs). Models of this type have important applications as predictive tools to estimate the likely impact of source-reduction strategies on human tissue concentrations.

Interspecies differences in cancer susceptibility and toxicity

One of the most complex challenges to the toxicologist represents extrapolation from laboratory animals to humans. In this article, we review interspecies differences in metabolism and toxicity of heterocyclic amines, aflatoxin B1, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), and related compounds, endocrine disrupters, polycyclic aromatic hydrocarbons, tamoxifen, and digitoxin. As far as possible, extrapolations to human toxicity and carcinogenicity are performed. Humans may be more susceptible to the carcinogenic effect of heterocyclic amines than monkeys, rats, and mice. Especially, individuals with high CYP1A2 and 3A4 activities and the rapid acetylator phenotype may be expected to have an increased risk. Striking interspecies variation in susceptibility to aflatoxin B1 carcinogenesis is known, with rats representing the most sensitive and mice the most resistant species, refractory to dietary levels three orders of magnitude higher than rats. An efficient conjugation with glutathione, catalyzed by glutathione S-transferase mYc, confers aflatoxin B1 resistance to mice. Extremely large interspecies differences in TCDD-induced toxicity are known. The guinea pig is the most susceptible mammal known, with an LD50 in the range 1-2 micrograms TCDD/kg, whereas the hamster is the most resistant species with an LD50 greater than 3000 micrograms/kg. A number of experts have pointed out to the fact that humans appear to be less sensitive to TCDD than most laboratory animals. Human exposure to background levels of TCDD is not likely to cause an incremental cancer risk. A clear cause--effect relationship has been shown between environmental endocrine-disrupting contaminants and adverse health effects in wildlife, whereas the effects seem to be less critical for humans. Studies on DNA adduct formation and metabolism of the nonsteroidal antiestrogen tamoxifen indicate that rats and mice are orders of magnitude more susceptible than humans.

Dioxin-like and non-dioxin-like toxic effects of polychlorinated biphenyls (PCBs): implications for risk assessment

Polychlorinated biphenyls (PCBs) are persistent, bioaccumulative, and toxic contaminants in the environment. Individual PCB congeners exhibit different physicochemical properties and biological activities that result in different environmental distributions and toxicity profiles. The variable composition of PCB residues in environmental matrices and their different mechanisms of toxicity complicate the development of scientifically based regulations for the risk assessment. In this article various approaches for the assessment of risks of PCBs have been critically examined. Recent developments in the toxic equivalency factor (TEF) approach for the assessment of toxic effects due to dioxin-like PCBs have been examined. PCB exposure studies that describe non-dioxin-like toxic effects, particularly neurobehavioral effects and their effective doses in animals were compiled. A comparative assessment of effective doses for dioxin-like and non-dioxin-like effects by PCBs has been made to evaluate the relative significance of non-ortho-and ortho-substituted PCBs in risk assessment. Using mink as an example, relative merits and implications of using TEF and total PCB approaches for assessing the potential for toxic effects in wildlife was examined. There are several advantages and limitations associated with each method used for PCB risk assessment. Toxic effects due to coplanar PCBs occur at relatively smaller concentrations than those due to non-dioxin-like PCBs and therefore the TEF approach derives the risk assessment of PCBs, in the environment. The need for the refinement of TEF approach for more accurate assessment of risks is discussed.