Assessing the impact of the duration and intensity of inhalation exposure on the magnitude of the variability of internal dose metrics in children and adults

Exposure limits for indoor volatile substances concerning the general population: The role of population-based differences in sensory irritation of the eyes and airways for assessment factors

Assessment factors (AFs) are essential in the derivation of occupational exposure limits (OELs) and indoor air quality guidelines. The factors shall accommodate differences in sensitivity between subgroups, i.e., workers, healthy and sick people, and occupational exposure versus life-long exposure for the general population. Derivation of AFs itself is based on empirical knowledge from human and animal exposure studies with immanent uncertainty in the empirical evidence due to knowledge gaps and experimental reliability. Sensory irritation in the eyes and airways constitute about 30-40% of OELs and is an abundant symptom in non-industrial buildings characterizing the indoor air quality and general health. Intraspecies differences between subgroups of the general population should be quantified for the proposal of more 'empirical' based AFs. In this review, we focus on sensitivity differences in sensory irritation about gender, age, health status, and vulnerability in people, based solely on human exposure studies. Females are more sensitive to sensory irritation than males for few volatile substances. Older people appear less sensitive than younger ones. However, impaired defense mechanisms may increase vulnerability in the long term. Empirical evidence of sensory irritation in children is rare and limited to children down to the age of six years. Studies of the nervous system in children compared to adults suggest a higher sensitivity in children; however, some defense mechanisms are more efficient in children than in adults. Usually, exposure studies are performed with healthy subjects. Exposure studies with sick people are not representative due to the deselection of subjects with moderate or severe eye or airway diseases, which likely underestimates the sensitivity of the group of people with diseases. Psychological characterization like personality factors shows that concentrations of volatile substances far below their sensory irritation thresholds may influence the sensitivity, in part biased by odor perception. Thus, the protection of people with extreme personality traits is not feasible by an AF and other mitigation strategies are required. The available empirical evidence comprising age, lifestyle, and health supports an AF of not greater than up to 2 for sensory irritation. Further, general AFs are discouraged for derivation, rather substance-specific derivation of AFs is recommended based on the risk assessment of empirical data, deposition in the airways depending on the substance's water solubility and compensating for knowledge and experimental gaps. Modeling of sensory irritation would be a better 'empirical' starting point for derivation of AFs for children, older, and sick people, as human exposure studies are not possible (due to ethical reasons) or not generalizable (due to self-selection). Dedicated AFs may be derived for environments where dry air, high room temperature, and visually demanding tasks aggravate the eyes or airways than for places in which the workload is balanced, while indoor playgrounds might need other AFs due to physical workload and affected groups of the general population.

Health Risk Assessment in Children Occupationally and Para-Occupationally Exposed to Benzene Using a Reverse-Translation PBPK Model

Benzene is a known human carcinogen and one of the ten chemicals of major public health concern identified by the World Health Organization. Our objective was to evaluate benzene's carcinogenic and non-carcinogenic health risks (current and projected) in highly exposed children in Yucatan, Mexico. Benzene exposure was estimated through a reverse-translation, four-compartment, physiologically based pharmacokinetic model (PBPK) based on previously performed urine trans, trans-muconic acid (benzene metabolite) determinations. Using a risk assessment methodology, the carcinogenic and non-carcinogenic risks of benzene were estimated for 6-12-year-old children from a family of shoemakers. The children's hazard quotients for decreased lymphocyte count were 27 and 53 for 4 and 8 h/day exposure, respectively, and 37 for the projected 8 h/day exposure in adults. The risks of developing leukemia were 2-6 cases in 1000 children exposed 4 h/day; 4-10 cases in 1000 children exposed 8 h/day, and 2-9 cases in 1000 adults with an 8 h/day lifetime exposure. Children in Yucatan working in shoe-manufacturing workshops, or living next to them, are exposed to benzene concentrations above the reference concentration and have unacceptably high risks of presenting with non-carcinogenic and carcinogenic hematologic symptoms, now and in the future. Interventions to prevent further exposure and mitigate health risks are necessary.

Evaluating Phenotypic and Transcriptomic Responses Induced by Low-Level VOCs in Zebrafish: Benzene as an Example

Urban environments are plagued by complex mixtures of anthropogenic volatile organic compounds (VOCs), such as mixtures of benzene, toluene, ethylene, and xylene (BTEX). Sources of BTEX that drive human exposure include vehicle exhaust, industrial emissions, off-gassing of building material, as well as oil spillage and leakage. Among the BTEX mixture, benzene is the most volatile compound and has been linked to numerous adverse health outcomes. However, few studies have focused on the effects of low-level benzene on exposure during early development, which is a susceptible window when hematological, immune, metabolic, and detoxification systems are immature. In this study, we used zebrafish to conduct a VOC exposure model and evaluated phenotypic and transcriptomic responses following 0.1 and 1 ppm benzene exposure during the first five days of embryogenesis (n = 740 per treatment). The benzene body burden was 2 mg/kg in 1 ppm-exposed larval zebrafish pools and under the detection limit in 0.1 ppm-exposed fish. No observable phenotypic changes were found in both larvae except for significant skeletal deformities in 0.1 ppm-exposed fish (p = 0.01) compared with unexposed fish. Based on transcriptomic responses, 1 ppm benzene dysregulated genes that were implicated with the development of hematological system, and the regulation of oxidative stress response, fatty acid metabolism, immune system, and inflammatory response, including apob, nfkbiaa, serpinf1, foxa1, cyp2k6, and cyp2n13 from the cytochrome P450 gene family. Key genes including pik3c2b, pltp, and chia.2 were differentially expressed in both 1 and 0.1 ppm exposures. However, fewer transcriptomic changes were induced by 0.1 ppm compared with 1 ppm. Future studies are needed to determine if these transcriptomic responses during embryogenesis have long-term consequences at levels equal to or lower than 1 ppm.

Estimation of toluene exposure in air from BMA (S-benzylmercapturic acid) urinary measures using a reverse dosimetry approach based on physiologically pharmacokinetic modeling

This study aimed to use a reverse dosimetry PBPK modeling approach to estimate toluene atmospheric exposure from urinary measurements of S-benzylmercapturic acid (BMA) in a small group of individuals and to evaluate the uncertainty associated to urinary spot-sampling compared to 24-h collected urine samples. Each exposure assessment technique was developed namely to estimate toluene air exposure from BMA measurements in 24-h urine samples (24-h-BMA) and from distributions of daily urinary BMA spot measurements (DUBSM). Model physiological parameters were described based upon age, weight, size and sex. Monte Carlo simulations with the PBPK model allowed converting DUBSM distribution (and 24-h-BMA) into toluene air levels. For the approach relying on DUBSM distribution, the ratio between the 95% probability of predicted toluene concentration and its 50% probability in each individual varied between 1.2 and 1.4, while that based on 24-h-BMA varied between 1.0 and 1.1. This suggests more variability in estimated exposure from spot measurements. Thus, estimating toluene exposure based on DUBSM distribution generated about 20% more uncertainty. Toluene levels estimated (0.0078-0.0138 ppm) are well below Health Canada's maximum chronic air guidelines. PBPK modeling and reverse dosimetry may be combined to interpret urinary metabolites data of VOCs and assess related uncertainties.

A Review of Volatile Organic Compound Contamination in Post-Industrial Urban Centers: Reproductive Health Implications Using a Detroit Lens

Volatile organic compounds (VOCs) are a group of aromatic or chlorinated organic chemicals commonly found in manufactured products that have high vapor pressure, and thus vaporize readily at room temperature. While airshed VOCs are well studied and have provided insights into public health issues, we suggest that belowground VOCs and the related vapor intrusion process could be equally or even more relevant to public health. The persistence, movement, remediation, and human health implications of subsurface VOCs in urban landscapes remain relatively understudied despite evidence of widespread contamination. This review explores the state of the science of subsurface movement and remediation of VOCs through groundwater and soils, the linkages between these poorly understood contaminant exposure pathways and health outcomes based on research in various animal models, and describes the role of these contaminants in human health, focusing on birth outcomes, notably low birth weight and preterm birth. Finally, this review provides recommendations for future research to address knowledge gaps that are essential for not only tackling health disparities and environmental injustice in post-industrial cities, but also protecting and preserving critical freshwater resources.

A generic PBTK model implemented in the MCRA platform: Predictive performance and uses in risk assessment of chemicals

Physiologically-based toxicokinetic (PBTK) models are important tools for in vitro to in vivo or inter-species extrapolations in health risk assessment of foodborne and non-foodborne chemicals. Here we present a generic PBTK model implemented in the EuroMix toolbox, MCRA 9 and predict internal kinetics of nine chemicals (three endocrine disrupters, three liver steatosis inducers, and three developmental toxicants), in data-rich and data-poor conditions, when increasingly complex levels of parametrization are applied. At the first stage, only QSAR models were used to determine substance-specific parameters, then some parameter values were refined by estimates from substance-specific or high-throughput in vitro experiments. At the last stage, elimination or absorption parameters were calibrated based on available in vivo kinetic data. The results illustrate that parametrization plays a capital role in the output of the PBTK model, as it can change how chemicals are prioritized based on internal concentration factors. In data-poor situations, estimates can be far from observed values. In many cases of chronic exposure, the PBTK model can be summarized by an external to internal dose factor, and interspecies concentration factors can be used to perform interspecies extrapolation. We finally discuss the implementation and use of the model in the MCRA risk assessment platform.

Model qualification of the PK-Sim® pediatric module for pediatric exposure assessment of CYP450 metabolized compounds

Pediatric physiologically based pharmacokinetic (PBPK) models facilitate the estimation of pharmacokinetic (PK) parameters in children under specific exposure conditions. In human health risk assessment, PBPK modeling has been used to determine a chemical-specific human kinetic adjustment factor (HKAF). Due to increased demands in regulatory assessment, model evaluation and qualification have gained growing attention. The aim of this study was to undertake model qualification of pediatric PBPK models for compounds that are primarily metabolized by cytochrome P450 (CYP) enzymes. The objectives were to determine the appropriateness of the virtual individual creating algorithm in PK-Sim® in predicting PK parameters and their variability in children and identify critical system-specific inputs. PBPK models in adults were constructed for several pharmaceuticals (grouped by major clearance process such as CYP3A4). Several age groups of virtual individuals were created to represent children in pediatric clinical studies. The mean and variance of clearance (CL) from virtual populations were compared to observed values. Sensitivity analysis on area under the curve (AUC) was performed. System-specific parameters of virtual children that contribute to inter-individual PK properties were assessed. Eighty-one percent of the comparisons between simulated and observed clearance values were within twofold error. The mean fold errors were 1.1, 1, 0.7 and 1.8 in adolescents, children, infants and neonates, respectively. CL variability was reasonably predicted for 70% of the comparisons with comparable coefficients of variation between observed and predicted. The sensitivity analysis revealed that fraction unbound in plasma, parameters related to CYP enzyme-mediated metabolism and liver volumewere most important in the estimation of pediatric exposure. A comparison of variabilities in weight, height and liver volume in virtual children showed reliable agreement with observed data. The presented results of predictive performance and properties of virtual populations provide confidence in the use of PK-Sim for pediatric PBPK modeling in toxicological applications including PBPK-based-HKAF derivation.

Evaluation of potential health effects associated with occupational and environmental exposure to styrene - an update

The potential chronic health risks of occupational and environmental exposure to styrene were evaluated to update health hazard and exposure information developed since the Harvard Center for Risk Analysis risk assessment for styrene was performed in 2002. The updated hazard assessment of styrene's health effects indicates human cancers and ototoxicity remain potential concerns. However, mechanistic research on mouse lung tumors demonstrates these tumors are mouse-specific and of low relevance to human cancer risk. The updated toxicity database supports toxicity reference levels of 20 ppm (equates to 400 mg urinary metabolites mandelic acid + phenylglyoxylic acid/g creatinine) for worker inhalation exposure and 3.7 ppm and 2.5 mg/kg bw/day, respectively, for general population inhalation and oral exposure. No cancer risk value estimates are proposed given the established lack of relevance of mouse lung tumors and inconsistent epidemiology evidence. The updated exposure assessment supports inhalation and ingestion routes as important. The updated risk assessment found estimated risks within acceptable ranges for all age groups of the general population and workers with occupational exposures in non-fiber-reinforced polymer composites industries and fiber-reinforced polymer composites (FRP) workers using closed-mold operations or open-mold operations with respiratory protection. Only FRP workers using open-mold operations not using respiratory protection have risk exceedances for styrene and should be considered for risk management measures. In addition, given the reported interaction of styrene exposure with noise, noise reduction to sustain levels below 85 dB(A) needs be in place.

Reverse dosimetry modeling of toluene exposure concentrations based on biomonitoring levels from the Canadian health measures survey

Biomonitoring might provide useful estimates of population exposure to environmental chemicals. However, data uncertainties stemming from interindividual variability are common in large population biomonitoring surveys. Physiologically based pharmacokinetic (PBPK) models might be used to account for age- and gender-related variability in internal dose. The objective of this study was to reconstruct air concentrations consistent with blood toluene measures reported in the third Canadian Health Measures Survey using reverse dosimetry PBPK modeling techniques. Population distributions of model's physiological parameters were described based upon age, weight, and size for four subpopulations (12-19, 20-39, 40-59, and 60-79 years old). Monte Carlo simulations applied to PBPK modeling allowed converting the distributions of venous blood measures of toluene obtained from CHMS into related air levels. Based upon blood levels observed at the 50^th, 90^th and 95^th percentiles, corresponding air toluene concentrations were estimated for teenagers aged 12-19 years as being, respectively, 0.009, 0.04 and 0.06 ppm. Similarly, values were computed for adults aged 20-39 years (0.007, 0.036, and 0.06 ppm), 40-59 years (0.007, 0.036 and 0.06 ppm) and 60-79 years (0.006, 0.022 and 0.04 ppm). These estimations are well below Health Canada's maximum recommended chronic air guidelines for toluene. In conclusion, PBPK modeling and reverse dosimetry may be combined to help interpret biomonitoring data for chemical exposure in large population surveys and estimate the associated toxicological health risk.

Assessing human variability in kinetics for exposures to multiple environmental chemicals: a physiologically based pharmacokinetic modeling case study with dichloromethane, benzene, toluene, ethylbenzene, and m-xylene

The objective of this study was to compare the magnitude of interindividual variability in internal dose for inhalation exposure to single versus multiple chemicals. Physiologically based pharmacokinetic models for adults (AD), neonates (NEO), toddlers (TODD), and pregnant women (PW) were used to simulate inhalation exposure to "low" (RfC-like) or "high" (AEGL-like) air concentrations of benzene (Bz) or dichloromethane (DCM), along with various levels of toluene alone or toluene with ethylbenzene and xylene. Monte Carlo simulations were performed and distributions of relevant internal dose metrics of either Bz or DCM were computed. Area under the blood concentration of parent compound versus time curve (AUC)-based variability in AD, TODD, and PW rose for Bz when concomitant "low" exposure to mixtures of increasing complexities occurred (coefficient of variation (CV) = 16-24%, vs. 12-15% for Bz alone), but remained unchanged considering DCM. Conversely, AUC-based CV in NEO fell (15 to 5% for Bz; 12 to 6% for DCM). Comparable trends were observed considering production of metabolites (AMET), except for NEO's CYP2E1-mediated metabolites of Bz, where an increased CV was observed (20 to 71%). For "high" exposure scenarios, Cmax-based variability of Bz and DCM remained unchanged in AD and PW, but decreased in NEO (CV= 11-16% to 2-6%) and TODD (CV= 12-13% to 7-9%). Conversely, AMET-based variability for both substrates rose in every subpopulation. This study analyzed for the first time the impact of multiple exposures on interindividual variability in toxicokinetics. Evidence indicates that this impact depends upon chemical concentrations and biochemical properties, as well as the subpopulation and internal dose metrics considered.

Characterization of the human kinetic adjustment factor for the health risk assessment of environmental contaminants

A default uncertainty factor of 3.16 (√10) is applied to account for interindividual variability in toxicokinetics when performing non-cancer risk assessments. Using relevant human data for specific chemicals, as WHO/IPCS suggests, it is possible to evaluate, and replace when appropriate, this default factor by quantifying chemical-specific adjustment factors for interindividual variability in toxicokinetics (also referred to as the human kinetic adjustment factor, HKAF). The HKAF has been determined based on the distributions of pharmacokinetic parameters (e.g., half-life, area under the curve, maximum blood concentration) in relevant populations. This article focuses on the current state of knowledge of the use of physiologically based algorithms and models in characterizing the HKAF for environmental contaminants. The recent modeling efforts on the computation of HKAF as a function of the characteristics of the population, chemical and its mode of action (dose metrics), as well as exposure scenario of relevance to the assessment are reviewed here. The results of these studies, taken together, suggest the HKAF varies as a function of the sensitive subpopulation and dose metrics of interest, exposure conditions considered (route, duration, and intensity), metabolic pathways involved and theoretical model underlying its computation. The HKAF seldom exceeded the default value of 3.16, except in very young children (i.e., <≈ 3 months) and when the parent compound is the toxic moiety. Overall, from a public health perspective, the current state of knowledge generally suggest that the default uncertainty factor is sufficient to account for human variability in non-cancer risk assessments of environmental contaminants.

Assessing the impact of child/adult differences in hepatic first-pass effect on the human kinetic adjustment factor for ingested toxicants

The objective of this study was to evaluate the impact of interindividual differences in hepatic first-pass effect (FPE) on the magnitude of the human kinetic adjustment factor (HKAF) for ingested toxicants. This factor aims at replacing a default value of 3.2 used in non-cancer risk assessment. Coupled with Monte Carlo simulations, steady-state equations that account for FPE were used to obtain distributions of arterial blood concentrations (CAss) and rates of metabolism in adults, neonates, infants and toddlers continuously exposed to an oral dose of 1 μg/kg/d of theoretical CYP2E1 and CYP1A2 substrates. For such substrates exhibiting a range of blood:air partition coefficients (Pb: 1-10,000) and hepatic extraction ratios in an average adult (E(ad): 0.01-0.99), HKAFs were computed as the ratio of the 95th percentile of dose metrics for each subpopulation over the 50th percentile value in adults. The reduced hepatic enzyme content in neonates as compared to adults resulted in correspondingly diminished FPE. Consequently, HKAFs greater than 3.2 could be observed, based on CAss only, in the following cases: for some CYP2E1 substrates with E(ad) ≤ 0.3, in neonates (max.: 6.3); and for some CYP1A2 substrates with E(ad) ≤ 0.1 and 0.7, in, respectively, neonates and infants (max.: 28.3). Overall, this study pointed out the importance of accounting for child/adult differences in FPE when determining the HKAF for oral exposure.