Challenges in the application of quantitative approaches in risk assessment: a case study with di-(2-ethylhexyl)phthalate

A physiologically based pharmacokinetic model of diethyl phthalates in humans

Phthalates are a family of industrial and consumer product chemicals, among which diethyl phthalate (DEP) has been widely used. DEP is metabolized into the active metabolite monoethyl phthalate (MEP) and exposure to DEP may induce male reproductive toxicity, developmental toxicity and hepatotoxicity. To better assess the toxicity of DEP and MEP, it is important to understand and predict their internal concentrations, especially in reproductive organs. Here we present a human physiologically based pharmacokinetic (PBPK) model of DEP. Implemented in R, the PBPK model consists of seven tissue compartments, including blood, gut, liver, fat, skin, gonad, and rest of body (RB). In the blood both DEP and MEP partition into free and bound forms, and tissue distribution is considered as blood flow-limited. DEP is metabolized in the gut and liver into MEP which is further glucuronidated and cleared through the urine. The chemical-specific parameters of the model were predicted in silico or estimated based on published human urinary MEP data after exposure to DEP in the air at 250 or 300 μg/m3 for 3 or 6 h through inhalation and dermal absorption. Sensitivity analysis identified important parameters including partition coefficients of DEP for fat, RB, and skin compartments, and the rate constants for glucuronidation of MEP and urinary excretion, with regard to Cmax, area under the curve (AUC), and clearance half-lives of DEP and MEP. A subset of the sensitive parameters was then included in hierarchical population Bayesian Markov chain Monte Carlo (MCMC) simulations to characterize the uncertainty and variability of these parameters. The model is consistent with the notion that dermal absorption represents a significant route of exposure to DEP in ambient air and clothing can be an effective barrier. The developed human PBPK model can be utilized upon further refinement as a quantitative tool for DEP risk assessment.

Multigenerational Effects of a Complex Human-Relevant Exposure during Folliculogenesis and Preimplantation Embryo Development: The FEDEXPO Study

Animal toxicological studies often fail to mimic the complexity of the human exposome, associating low doses, combined molecules and long-term exposure. Since the reproductive potential of a woman begins in the fetal ovary, the literature regarding the disruption of its reproductive health by environmental toxicants remains limited. Studies draw attention to follicle development, a major determinant for the quality of the oocyte, and the preimplantation embryo, as both of them are targets for epigenetic reprogramming. The "Folliculogenesis and Embryo Development EXPOsure to a mixture of toxicants: evaluation in the rabbit model" (FEDEXPO) project emerged from consideration of these limitations and aims to evaluate in the rabbit model the impacts of an exposure to a mixture of known and suspected endocrine disrupting chemicals (EDCs) during two specific windows, including folliculogenesis and preimplantation embryo development. The mixture combines eight environmental toxicants, namely perfluorooctanesulfonic acid (PFOS), perfluorooctanoic acid (PFOA), dichlorodiphenyldichloroethylene (DDE), hexachlorobenzene (HCB), β-hexachlorocyclohexane (β-HCH), 2,2'4,4'-tetrabromodiphenyl ether (BDE-47), di(2-ethylhexyl) phthalate (DEHP) and bisphenol S (BPS), at relevant exposure levels for reproductive-aged women based on biomonitoring data. The project will be organized in order to assess the consequences of this exposure on the ovarian function of the directly exposed F0 females and monitor the development and health of the F1 offspring from the preimplantation stage. Emphasis will be made on the reproductive health of the offspring. Lastly, this multigenerational study will also tackle potential mechanisms for the inheritance of health disruption via the oocyte or the preimplantation embryo.

Comparison of the exposure assessment of di(2-ethylhexyl) phthalate between the PBPK model-based reverse dosimetry and scenario-based analysis: A Korean general population study

Di (2-ethylhexyl) phthalate (DEHP), classified as a reproductive toxicant, is a ubiquitous pollutant in foodstuffs, dust, and commercial products. In this study, to provide a useful cross-check on the accuracy of the exposure assessment, the estimated daily intake of DEHP was compared using reverse dosimetry with a physiologically-based pharmacokinetic (PBPK) model and a scenario-based probabilistic estimation model for six subpopulations in Korea. For reverse dosimetry analysis, the concentrations of urinary DEHP metabolites, namely mono (2-ethyl-5-hydroxyhexyl) phthalate (MEHHP) and mono (2-ethyl-5-oxohexyl)phthalate (MEOHP), from three human biomonitoring program datasets were used. For the scenario-based model, we evaluated the various exposure sources of DEHP, including diet, air, indoor dust, soil, and personal care products (PCPs), and also determined its levels based on the literature review and measurements of indoor dust. The DEHP exposure doses using both exposure assessment approaches were similar in all cases, except for the 95th percentile exposure doses in toddlers (1-2 years) and young children (3-6 years). The PBPK-reverse dosimetry estimated daily intakes at the 95th percentile ranged between 22.53 and 29.90 μg/kg/day for toddlers and young children. These exceeded the reference dose (RfD) of 20 μg/kg bw/day of the US Environmental Protection Agency (EPA) based on the increased relative liver weight. Although, food was considered the primary source of DEHP, contributing to a total exposure of 50.8-75.1%, the effect of exposure to indoor dust should not be overlooked. The occurrence of high levels of DEHP in indoor dust collected from Korean homes suggests the use of a wide variety of consumer products containing DEHP. Furthermore, more attention should be paid to the high exposure levels of DEHP, especially in young children. Therefore, it is necessary to perform continuous monitoring of the indoor dust, consumer products, and the body burden of children.

Use of Mixture Dosing and Nonlinear Mixed Effect Modeling of Eight Environmental Contaminants in Rabbits to Improve Extrapolation Value of Toxicokinetic Data

BACKGROUND

Although in vivo studies of internal exposure to hazardous substances have been carried out for many years, there is room for progress to improve their informative value while adhering to the four R's: replacement, reduction, refinement, and responsibility rule.

OBJECTIVES

The objective of the study was to illustrate how toxicokinetic (TK) study design and data analysis can be implemented under the 4R rule to plan a chronic dosage regimen for investigating TK/toxicodynamic (TD) relationships.

METHODS

The intravenous (IV) and oral serum concentrations of eight hazardous environmental contaminants including 1,1-Dichloro-2,2-bis(p-chlorophenyl)ethylene (pp'DDE), ß-Hexachlorocyclohexane (β-HCH), hexachlorobenzene (HCB), 2,2'4,4'-tetrabromodiphenyl ether (BDE-47), perfluorooctane sulfonate (PFOS), perfluorooctanoic acid (PFOA), di(2ethylhexyl)phthalate (DEHP), and bisphenol S (BPS) were obtained after mixture dosing in rabbits using a sparse sampling design. Data were comprehensively analyzed using nonlinear mixed effect (NLME) modeling.

RESULTS

The short persistence of BPS and of the DEHP metabolite (mono-2-ethylhexyl phthalate), reflected by their mean residence times (MRT) of a few hours, was due to their efficient clearance (CL, 3.2 and 0.47L/kg/h). The longer MRT of the other compounds (1-48 d) resulted either from their extremely low clearance (lower than 0.01L/kg/h for PFOA and PFOS) or from their very large volume of distribution (VSS) ranging from 33 to 45L/kg. Estimates of CL, VSS, and bioavailability were used to compute the oral loading and daily maintenance doses required to attain a nominal steady-state serum concentration of 1 ng/mL. Simulations with the NLME model were applied to predict the serum concentration profile and to contrast the differential rates of accumulation in the central vs. peripheral compartments.

CONCLUSION

NLME modeling of the IV and oral TK of hazardous environmental contaminants, in rabbits while fulfilling the 4R rule, was able to provide the physiological basis for interspecies extrapolation of exposure rates in a TK/TD approach to risk assessment. https://doi.org/10.1289/EHP8957.

Using quantitative modeling tools to assess pharmacokinetic bias in epidemiological studies showing associations between biomarkers and health outcomes at low exposures

Biomarkers of exposure can be measured at lower and lower levels due to advances in analytical chemistry. Using these sensitive methods, some epidemiology studies report associations between biomarkers and health outcomes at biomarker levels much below those associated with effects in animal studies. While some of these low exposure associations may arise from increased sensitivity of humans compared with animals or from species-specific responses, toxicology studies with drugs, commodity chemicals and consumer products have not generally indicated significantly greater sensitivity of humans compared with test animals for most health outcomes. In some cases, these associations may be indicative of pharmacokinetic (PK) bias, i.e., a situation where a confounding factor or the health outcome itself alters pharmacokinetic processes affecting biomarker levels. Quantitative assessment of PK bias combines PK modeling and statistical methods describing outcomes across large numbers of individuals in simulated populations. Here, we first provide background on the types of PK models that can be used for assessing biomarker levels in human population and then outline a process for considering PK bias in studies intended to assess associations between biomarkers and health outcomes at low levels of exposure. After providing this background, we work through published examples where these PK methods have been applied with several chemicals/chemical classes - polychlorinated biphenyls (PCBs), perfluoroalkyl substances (PFAS), polybrominated biphenyl ethers (PBDE) and phthalates - to assess the possibility of PK bias. Studies of the health effects of low levels of exposure will be improved by developing some confidence that PK bias did not play significant roles in the observed associations.

Development of a physiologically based pharmacokinetic model of diisononyl phthalate (DiNP) in pregnant rat and human

A physiologically based pharmacokinetic (PBPK) model for di-isononyl phthalate (DiNP) was developed by adapting the existing models for di(2-ethylhexyl) phthalate (DEHP) and di-butylphthalate (DBP). Both pregnant rat and human time-course plasma and urine data were used to address the hydrolysis of DiNP in intestinal tract, plasma, and liver as well as hepatic oxidative metabolism and conjugation of the monoester and primary oxidative metabolites. Data in both rats and humans were available to inform the uptake and disposition of mono-isononyl phthalate (MiNP) as well as the three primary oxidative metabolites including hydroxy (7-OH)-, oxo (7-OXO)-, and carboxy (7-COX)-monoisononyl phthalate in plasma and urine. The DiNP model was reliable over a wide range of exposure levels in the pregnant rat as well as the two low exposure levels in humans including capturing the nonlinear behavior in the pregnant rat after repeated 750 mg/kg/day dosing. The presented DiNP PBPK model in pregnant rat and human, based upon an extensive kinetic dataset in both species, may provide a basis for assessing human equivalent exposures based upon either rodent or in vitro points of departure.

Application of a combined aggregate exposure pathway and adverse outcome pathway (AEP-AOP) approach to inform a cumulative risk assessment: A case study with phthalates

Advancements in measurement and modeling capabilities are providing unprecedented access to estimates of chemical exposure and bioactivity. With this influx of new data, there is a need for frameworks that help organize and disseminate information on chemical hazard and exposure in a manner that is accessible and transparent. A case study approach was used to demonstrate integration of the Adverse Outcome Pathway (AOP) and Aggregate Exposure Pathway (AEP) frameworks to support cumulative risk assessment of co-exposure to two phthalate esters that are ubiquitous in the environment and that are associated with disruption of male sexual development in the rat: di(2-ethylhexyl) phthalate (DEHP) and di-n-butyl phthalate (DnBP). A putative AOP was developed to guide selection of an in vitro assay for derivation of bioactivity values for DEHP and DnBP and their metabolites. AEPs for DEHP and DnBP were used to extract key exposure data as inputs for a physiologically based pharmacokinetic (PBPK) model to predict internal metabolite concentrations. These metabolite concentrations were then combined using in vitro-based relative potency factors for comparison with an internal dose metric, resulting in an estimated margin of safety of ~13,000. This case study provides an adaptable workflow for integrating exposure and toxicity data by coupling AEP and AOP frameworks and using in vitro and in silico methodologies for cumulative risk assessment.

Physiology-based toxicokinetic modelling in the frame of the European Human Biomonitoring Initiative

Given the opportunities provided by internal dosimetry modelling in the interpretation of human biomonitoring (HBM) data, the assessment of the links between exposure to chemicals and observed HBM data can be effectively supported by PBTK modelling. This paper gives a comprehensive review of available human PBTK models for compounds selected as a priority by the European Human Biomonitoring Initiative (HBM4EU). We highlight their advantages and deficiencies and suggest steps for advanced internal dose modelling. The review of the available PBTK models highlighted the conceptual differences between older models compared to the ones developed recently, reflecting commensurate differences in research questions. Due to the lack of coordinated strategies for deriving useful biomonitoring data for toxicokinetic properties, significant problems in model parameterisation still remain; these are further increased by the lack of human toxicokinetic data due to ethics issues. Finally, questions arise as well as to the extent they are really representative of interindividual variability. QSARs for toxicokinetic properties is a complementary approach for PBTK model parameterisation, especially for data poor chemicals. This approach could be expanded to model chemico-biological interactions such as intestinal absorption and renal clearance; this could serve the development of more complex generic PBTK models that could be applied to newly derived chemicals. Another gap identified is the framework for mixture interaction terms among compounds that could eventually interact in metabolism. From the review it was concluded that efforts should be shifted toward the development of generic multi-compartmental and multi-route models, supported by targeted biomonitoring coupled with parameterisation by both QSAR approach and experimental (in-vivo and in-vitro) data for newly developed and data poor compounds.

Excretion of Di-2-ethylhexyl phthalate (DEHP) metabolites in urine is related to body mass index because of higher energy intake in the overweight and obese

A higher body mass index (BMI) has been positively associated with the rate of excretion of di-2-ethylhexyl phthalate (DEHP) metabolites in urine in data from the National Health and Nutrition Examination Survey (NHANES), suggesting an association between DEHP exposure and BMI. The association, however, may be due to the association between body mass maintenance and higher energy intake, with higher energy intake being accompanied by a higher intake of DEHP. To examine this hypothesis, we ran a Monte Carlo simulation with a DEHP physiologically-based pharmacokinetic (PBPK) model for adult humans. A realistic exposure sub-model was used, which included the relation of body weight to energy intake and of energy intake to DEHP intake. The model simulation output, when compared with urinary metabolite data from NHANES, supported good model validity. The distribution of BMI in the simulated population closely resembled that in the NHANES population. This indicated that the simulated subjects and DEHP exposure model were closely aligned with the NHANES population of interest. In the simulated population, the ordinary least squares regression coefficient for log(BMI) as a function of log(DEHP nmol/min) was 0.048 (SE 0.001), as compared with the reported value of 0.019 (SE 0.005). In other words, given our model structure, the higher energy intake in the overweight and obese, and the concomitant higher DEHP exposure, describes the reported relationship between BMI and DEHP.

The PPARα-dependent rodent liver tumor response is not relevant to humans: addressing misconceptions

A number of industrial chemicals and therapeutic agents cause liver tumors in rats and mice by activating the nuclear receptor peroxisome proliferator-activated receptor α (PPARα). The molecular and cellular events by which PPARα activators induce rodent hepatocarcinogenesis have been extensively studied elucidating a number of consistent mechanistic changes linked to the increased incidence of liver neoplasms. The weight of evidence relevant to the hypothesized mode of action (MOA) for PPARα activator-induced rodent hepatocarcinogenesis is summarized here. Chemical-specific and mechanistic data support concordance of temporal and dose-response relationships for the key events associated with many PPARα activators. The key events (KE) identified in the MOA are PPARα activation (KE1), alteration in cell growth pathways (KE2), perturbation of hepatocyte growth and survival (KE3), and selective clonal expansion of preneoplastic foci cells (KE4), which leads to the apical event-increases in hepatocellular adenomas and carcinomas (KE5). In addition, a number of concurrent molecular and cellular events have been classified as modulating factors, because they potentially alter the ability of PPARα activators to increase rodent liver cancer while not being key events themselves. These modulating factors include increases in oxidative stress and activation of NF-kB. PPARα activators are unlikely to induce liver tumors in humans due to biological differences in the response of KEs downstream of PPARα activation. This conclusion is based on minimal or no effects observed on cell growth pathways and hepatocellular proliferation in human primary hepatocytes and absence of alteration in growth pathways, hepatocyte proliferation, and tumors in the livers of species (hamsters, guinea pigs and cynomolgus monkeys) that are more appropriate human surrogates than mice and rats at overlapping dose levels. Despite this overwhelming body of evidence and almost universal acceptance of the PPARα MOA and lack of human relevance, several reviews have selectively focused on specific studies that, as discussed, contradict the consensus opinion and suggest uncertainty. In the present review, we systematically address these most germane suggested weaknesses of the PPARα MOA.

Using exposure prediction tools to link exposure and dosimetry for risk-based decisions: A case study with phthalates

A few different exposure prediction tools were evaluated for use in the new in vitro-based safety assessment paradigm using di-2-ethylhexyl phthalate (DEHP) and dibutyl phthalate (DnBP) as case compounds. Daily intake of each phthalate was estimated using both high-throughput (HT) prediction models such as the HT Stochastic Human Exposure and Dose Simulation model (SHEDS-HT) and the ExpoCast heuristic model and non-HT approaches based on chemical specific exposure estimations in the environment in conjunction with human exposure factors. Reverse dosimetry was performed using a published physiologically based pharmacokinetic (PBPK) model for phthalates and their metabolites to provide a comparison point. Daily intakes of DEHP and DnBP were estimated based on the urinary concentrations of their respective monoesters, mono-2-ethylhexyl phthalate (MEHP) and monobutyl phthalate (MnBP), reported in NHANES (2011-2012). The PBPK-reverse dosimetry estimated daily intakes at the 50th and 95th percentiles were 0.68 and 9.58 μg/kg/d and 0.089 and 0.68 μg/kg/d for DEHP and DnBP, respectively. For DEHP, the estimated median from PBPK-reverse dosimetry was about 3.6-fold higher than the ExpoCast estimate (0.68 and 0.18 μg/kg/d, respectively). For DnBP, the estimated median was similar to that predicted by ExpoCast (0.089 and 0.094 μg/kg/d, respectively). The SHEDS-HT prediction of DnBP intake from consumer product pathways alone was higher at 0.67 μg/kg/d. The PBPK-reverse dosimetry-estimated median intake of DEHP and DnBP was comparable to values previously reported for US populations. These comparisons provide insights into establishing criteria for selecting appropriate exposure prediction tools for use in an integrated modeling platform to link exposure to health effects.

A global human health risk assessment for octamethylcyclotetrasiloxane (D4)

Octamethylcyclotetrasiloxane (D₄) is a low-molecular-weight volatile cyclic siloxane, primarily used as an intermediate in the production of some widely-used industrial and consumer silicone based polymers and may be present as a component in a variety of consumer products. A global "harmonized" risk assessment was conducted to meet requirements for substance-specific risk assessments conducted by regulatory agencies such as USEPA's Integrated Risk Information System (IRIS), Health Canada's Chemical Management Program (CMP) and various independent scientific committees of the European Commission (e.g. the Scientific Committee on Consumer Safety (SCCS), the Scientific Committee on Health and Environmental Risks (SCHER)), as well as to provide guidance for chemical safety assessments under REACH in Europe. This risk assessment incorporates global exposure information combined with a Monte Carlo analysis to determine the most significant routes of exposure. Utilization of a multi-species, multi-route physiologically based pharmacokinetic (PBPK) model was included to estimate internal dose metrics, benchmark modeling was used to determine a point of departure (POD), and a margin of safety (MOS) evaluation was used to compare the estimates of intake with the POD. Because of the specific pharmacokinetic behaviors of D₄ including high lipophilicity, high volatility with low blood-to-air partition coefficients and an extensive metabolic clearance that regulates tissue dose after exposure, the use of a PBPK model was essential to provide a comparison of a dose metric that reflects these processes. The characterization of the potential for adverse effects after exposure to D₄ using a MOS approach based on an internal dose metric removes the subjective application of varying uncertainty factors from various regulatory agencies and allows examination of the differences between internal dose metrics associated with exposure and those associated with adverse effects.

Systematic review in chemical risk assessment - A chemical industry perspective

This commentary provides a perspective from the chemicals industry on the potential usefulness of systematic review approaches in chemical risk assessment.

Assessment of urinary metabolite excretion after rat acute exposure to perfluorooctanoic acid and other peroxisomal proliferators

Perfluorooctanoic acid (PFOA) is a persistent environmental contaminant. Activation of the peroxisome proliferator activated receptor alpha (PPARα) resulting from exposure to PFOA has been extensively studied in rodents. However, marked differences in response to peroxisome proliferators prevent extrapolation of rodent PPARα activation to human health risks and additional molecular mechanisms may also be involved in the biological response to PFOA exposure. To further explore the potential involvement of such additional pathways, the effects of PFOA exposure on urinary metabolites were directly compared with those of other well-known PPARα agonists. Male rats were administered PFOA (10, 33, or 100 mg/kg/d), fenofibrate (100 mg/kg/d), or di(2-ethylhexyl) phthalate (100 mg/kg/d) by gavage for 3 consecutive days and allowed to recover for 4 days, and overnight urine was collected. Greater urinary output was observed exclusively in PFOA-treated rats as the total fraction of PFOA excreted in urine increased with the dose administered. Assessment of urinary metabolites (ascorbic acid, quinolinic acid, 8-hydroxy-2'-deoxyguanosine, and malondialdehyde) provided additional information on PFOA's effects on hepatic glucuronic acid and tryptophan-nicotinamide adenine dinucleotide (NAD) pathways and on oxidative stress, whereas increased liver weight and palmitoyl-CoA oxidase activity indicative of PPARα activation and peroxisomal proliferation persisted up to day five after the last exposure.

Urinary environmental chemical concentrations and vitamin D are associated with vision, hearing, and balance disorders in the elderly

BACKGROUND

Acute exposure to environmental chemicals can result in loss of consciousness and upon recovery neurological symptoms, but little evidence exists in large epidemiological human studies. Hence, it was aimed to determine the relationships between urinary environmental chemicals (including heavy metals, environmental bisphenols, pesticides, arsenic, and phthalates) concentrations and vision, hearing, and balance disorders in a national population-based setting.

METHODS

United States National Health and Nutrition Examination Surveys is a national population-based multi-year cross-sectional study. Information on demographics and vision, hearing, and balance disorders was obtained by household interview using questionnaires in the 2003-2004 cohort (aged 50 and above). Urinary environmental chemicals were detected by mass spectrometry in selected but representative people. Analyses involved logistic regression models.

RESULTS

Urinary cadmium, molybdenum, and tungsten concentrations, which are commonly associated with heart disease, were associated with vision disorder. Urinary 2,4,5-trichlorophenol and arsenic acid concentrations and circulating mono-n-butyl phthalate, mono-benzyl phthalate, and mono-(3-carboxypropyl) phthalate metabolites were significantly associated with hearing disorder. Moreover, urinary barium and 4-tert-octyl-phenol concentrations were associated with balance disorder. People who had ears ringing, roaring, or buzzing in the past year tended to have higher urinary barium, 2,4-dichlorophenol, and mono-benzyl phthalate concentrations.

DISCUSSION

Significant correlations were observed in urinary environmental chemicals and neurobehavioural impairment for the first time. However, the causation cannot be established due to its cross-sectional study design. Future studies with a longitudinal aspect and/or in clinical trials are warranted to clearly understand the biological mechanism along the pathway before drawing a firm conclusion on these relationships.