Physiologically based pharmacokinetic modeling of the temperature-dependent dermal absorption of chloroform by humans following bath water exposures

Determination of chloroform concentration and human exposure assessment in the swimming pool

This cross-sectional study aimed to examine the concentration of the by-products of chlorination in the swimming pool and estimate human health risk for the swimmers of Shiraz University of Medical Sciences. In this study, the chloroform concentrations of 16 samples were measured using Gas Chromatography (GC). All the measured concentrations were less than the allowed amount announced by the World Health Organization (WHO). The results of the cancer risk (CR) and hazard index (HI) showed that the major exposure routes were found to be dermal during swimming and the 95 percentile of estimated CR and HI for the male group were 1.38 × 10^-10 and 1.82 × 10^-5 respectively, which is higher than the values of 5.48 × 10^-10 and 2.25 × 10^-5 respectively, for the women group. Sensitivity analyses indicated that the swimming exposure time (ET), and chloroform concentration were the most relevant variables in the health risk model. Therefore, knowledge about the sources of micro-pollutants in swimming pools might help promote the health methods of the pool environment.

Health risk assessment of trihalomethanes mixtures from daily water-related activities via multi-pathway exposure based on PBPK model

In this study, the concentrations of trihalomethanes (THMs) in tap water and direct drinking water were analyzed, and based on the human behavior patterns and building parameters, the concentrations of THMs in indoor air were simulated with the water-air concentration conversion model. In addition, concentrations of THMs in human tissues were predicted based on physiologically based pharmacokinetic (PBPK) model, and the health risk of THMs for participants were estimated. Furthermore, the carcinogenic risk of mixtures according to the method proposed by USEPA and PBPK model based method were calculated and compared. The concentrations of chloroform, bromodichloromethane, dibromochloromethane and bromoform in tap water were 11.28-16.21, 4.83-6.28, 0.81-1.32 and 0.08-0.21 μg/L, and those in direct drinking water were 3.29-6.88, 0.35-0.47, 0.03-0.08 and 0.04-0.08 μg/L, respectively. The results of water-air concentration conversion model demonstrated that pollutants in air had a strong correlation with water-related activities. Multi-pathway PBPK model showed that THMs concentrations in liver, kidney and richly perfused tissue were higher than those in other tissues. The results of risk assessment showed that the mean risk levels of mixtures were 1.69 × 10^-4 and 1.72 × 10^-4 calculated by the USEPA recommended method and PBPK based method, which seriously exceeded the acceptable level. TCM and BDCM were the major risk factors, and inhalation was the main exposure route of THMs.

Monte-Carlo and multi-exposure assessment for the derivation of criteria for disinfection byproducts and volatile organic compounds in drinking water: Allocation factors and liter-equivalents per day

The probability distributions of total potential doses of disinfection byproducts and volatile organic compounds via ingestion, inhalation, and dermal exposure were estimated with Monte Carlo simulations, after conducting physiologically based pharmacokinetic model simulations to takes into account the differences in availability between the three exposures. If the criterion that the 95th percentile estimate equals the TDI (tolerable daily intake) is regarded as protecting the majority of a population, the drinking water criteria would be 140 (trichloromethane), 66 (bromodichloromethane), 157 (dibromochloromethane), 203 (tribromomethane), 140 (dichloroacetic acid), 78 (trichloroacetic acid), 6.55 (trichloroethylene, TCE), and 22 μg/L (perchloroethylene). The TCE criterion was lower than the Japanese Drinking Water Quality Standard (10 μg/L). The latter would allow the intake of 20% of the population to exceed the TDI. Indirect inhalation via evaporation from water, especially in bathrooms, was the major route of exposure to compounds other than haloacetic acids (HAAs) and accounted for 1.2-9 liter-equivalents/day for the median-exposure subpopulation. The ingestion of food was a major indirect route of exposure to HAAs. Contributions of direct water intake were not very different for trihalomethanes (30-45% of TDIs) and HAAs (45-52% of TDIs).

Malathion dermal permeability in relation to dermal load: Assessment by physiologically based pharmacokinetic modeling of in vivo human data

Estimates of dermal permeability (K_p), obtained by fitting an updated human PBPK model for malathion to previously reported data on excreted urinary metabolites after 29 volunteers were dermally exposed to measured values of [¹⁴C]malathion dermal load (L), were used to examine the empirical relationship between K_p and L. The PBPK model was adapted from previously reported human biokinetic and PBPK models for malathion, fit to previously reported urinary excretion data after oral [¹⁴C]malathion intake by volunteers, and then augmented to incorporate a standard K_p approach to modeling dermal-uptake kinetics. Good to excellent PBPK-model fits were obtained to all of 29 sets of cumulative urinary metabolite-excretion data (ave. [±1 SD] R² = 0.953 [±0.064]). Contrary to the assumption that K_p and L are independent typically applied for dermally administered liquids or solutions, the 29 PBPK-based estimates of K_p obtained for malathion exhibit a strong positive association with the 2/3rds power of L (log-log Pearson correlation = 0.925, p = ∼0). Possible explanations of this observation involving physico-chemical characteristics and/or in vivo cutaneous effects of malathion are discussed. The PBPK model presented, and our observation that K_p estimates obtained by fitting this model to human experimental urinary-excretion data correlate well with L^2/3, allow more realistic assessments of absorbed and metabolized dose during or after a variety of scenarios involving actual or potential dermal or multi-route malathion exposures, including for pesticide workers or farmers who apply malathion to crops.

Sex Differences in Human and Animal Toxicology

Sex, the states of being female or male, potentially interacts with all xenobiotic exposures, both inadvertent and deliberate, and influences their toxicokinetics (TK), toxicodynamics, and outcomes. Sex differences occur in behavior, exposure, anatomy, physiology, biochemistry, and genetics, accounting for female-male differences in responses to environmental chemicals, diet, and pharmaceuticals, including adverse drug reactions (ADRs). Often viewed as an annoying confounder, researchers have studied only one sex, adjusted for sex, or ignored it. Occupational epidemiology, the basis for understanding many toxic effects in humans, usually excluded women. Likewise, Food and Drug Administration rules excluded women of childbearing age from drug studies for many years. Aside from sex-specific organs, sex differences and sex × age interactions occur for a wide range of disease states as well as hormone-influenced conditions and drug distribution. Women have more ADRs than men; the classic sex hormone paradigm (gonadectomy and replacement) reveals significant interaction of sex and TK including absorption, distribution, metabolisms, and elimination. Studies should be designed to detect sex differences, describe the mechanisms, and interpret these in a broad social, clinical, and evolutionary context with phenomena that do not differ. Sex matters, but how much of a difference is needed to matter remains challenging.

Using physiologically based pharmacokinetic models to estimate the health risk of mixtures of trihalomethanes from reclaimed water

To estimate the health risk of mixture of trihalomethanes (THMs) from reclaimed water during toilet flushing, the interaction-based Hazard Index (HI(interaction-based)) and the mixture carcinogenic risk (CRM) according to tissue dose were conducted through the integrated use of both the exposure concentrations model and the physiologically based pharmacokinetic (PBPK) model of THMs. Monte Carlo simulations were employed to implement the probabilistic risk analysis and sensitivity analysis. Nine samples were analyzed, which were collected from J Water Reclamation Plant (JWRP) in Tianjin of China. The results indicated that the mean HI(interaction-based) (=0.85) was lower than the acceptable risk level (=1). The probability that the HI(interaction-based) exceeded the acceptable risk level is 22.97%. For carcinogenic risk, the CRM ranges from 9.41×10(-7) to 3.54×10(-5), with a mean of 5.49×10(-6). Moreover, the probability of exceeding the acceptable risk level (1×10(-6)) is near 100%. And the values of HI(interaction-based) from sample no. 1, 5, and 7 exceeded 1, while the values of CRM for all samples exceeded 1×10(-6). Consequently, the reclaimed water used for flushing toilets should be paid more attention, though non-carcinogenic effect is relatively small. Furthermore, the concentrations of DBCM had greater impact on both the carcinogenic and non-carcinogenic risk based on sensitivity analysis.

Estimation of in vivo dose of dermally applied chemicals leading to estrogen/androgen receptor-mediated toxicity from in vitro data--Illustration with four reproductive toxicants

We present a quantitative in vitro-in vivo extrapolation framework enabling the estimation of the external dermal exposure dose from in vitro experimental data relevant to a toxicity pathway of interest. The framework adapts elements of the biological pathway altering dose (BPAD) method [Judson et al. Chem Res Toxicol 2011;24:451] to the case of dermal exposure. Dermal doses of four toxicants equivalent to concentrations characterizing their effect on estrogen receptor α or androgen receptor activity in chemical-activated luciferase expression (CALUX) assays are estimated. The analysis shows that dermal BPADs, calculated from one in vitro concentration, can differ by up to a factor of 55, due to the impact applied dose and dermal exposure scenarios can have on skin permeation kinetics. These features should therefore be taken into account in risk assessment of dermally applied chemicals.

Passive exposures of children to volatile trihalomethanes during domestic cleaning activities of their parents

Domestic cleaning has been proposed as a determinant of trihalomethanes (THMs) exposure in adult females. We hypothesized that parental housekeeping activities could influence children's passive exposures to THMs from their mere physical presence during domestic cleaning. In a recent cross-sectional study (n = 382) in Cyprus [41 children (< 18 y) and 341 adults (≥ 18 y)], we identified 29 children who met the study's inclusion criteria. Linear regression models were applied to understand the association between children sociodemographic variables, their individual practices influencing ingestion and noningestion exposures to ΣTHMs, and their urinary THMs levels. Among the children-specific variables, age alone showed a statistically significant inverse association with their creatinine-adjusted urinary ΣTHMs (rS = -0.59, p < 0.001). A positive correlation was observed between urinary ΣTHMs (ng g(-1)) of children and matched-mothers (rS = 0.52, p = 0.014), but this was not the case for their matched-fathers (rS = 0.39, p = 0.112). Time spent daily by the matched-mothers for domestic mopping, toilet and other cleaning activities using chlorine-based cleaning products was associated with their children's urinary THMs levels (rS = 0.56, p = 0.007). This trend was not observed between children and their matched-fathers urinary ΣTHMs levels, because of minimum amount of time spent by the latter in performing domestic cleaning. The proportion of variance of creatinine-unadjusted and adjusted urinary ΣTHMs levels in children that was explained by the matched-mothers covariates was 76% and 74% (p < 0.001), respectively. A physiologically-based pharmacokinetic model adequately predicted urinary chloroform excretion estimates, being consistent with the corresponding measured levels. Our findings highlighted the influence of mothers' domestic cleaning activities towards enhancing passive THMs exposures of their children. The duration of such activities could be further tested as a valid indicator of children's THMs body burden.

Dermal uptake of 18 dilute aqueous chemicals: in vivo disappearance-method measures greatly exceed in vitro-based predictions

Average rates of total dermal uptake (Kup ) from short-term (e.g., bathing) contact with dilute aqueous organic chemicals (DAOCs) are typically estimated from steady-state in vitro diffusion-cell measures of chemical permeability (Kp ) through skin into receptor solution. Widely used ("PCR-vitro") methods estimate Kup by applying diffusion theory to increase Kp predictions made by a physico-chemical regression (PCR) model that was fit to a large set of Kp measures. Here, Kup predictions for 18 DAOCs made by three PCR-vitro models (EPA, NIOSH, and MH) were compared to previous in vivo measures obtained by methods unlikely to underestimate Kup . A new PCR model fit to all 18 measures is accurate to within approximately threefold (r = 0.91, p < 10(-5) ), but the PCR-vitro predictions (r > 0.63) all tend to underestimate the Kup measures by mean factors (UF, and p value for testing UF = 1) of 10 (EPA, p < 10(-6) ), 11 (NIOSH, p < 10(-8) ), and 6.2 (MH, p = 0.018). For all three PCR-vitro models, log(UF) correlates negatively with molecular weight (r(2) = 0.31 to 0.84, p = 0.017 to < 10(-6) ) but not with log(vapor pressure) as an additional predictor (p > 0.05), so vapor pressure appears not to explain the significant in vivo/PCR-vitro discrepancy. Until this discrepancy is explained, careful in vivo measures of Kup should be obtained for more chemicals, the expanded in vivo database should be compared to in vitro-based predictions, and in vivo data should be considered in assessing aqueous dermal exposure and its uncertainty.

Studying permethrin exposure in flight attendants using a physiologically based pharmacokinetic model

Assessment of potential health risks to flight attendants from exposure to pyrethroid insecticides, used for aircraft disinsection, is limited because of (a) lack of information on exposures to these insecticides, and (b) lack of tools for linking these exposures to biomarker data. We developed and evaluated a physiologically based pharmacokinetic (PBPK) model to assess the exposure of flight attendants to the pyrethroid insecticide permethrin attributable to aircraft disinsection. The permethrin PBPK model was developed by adapting previous models for pyrethroids, and was parameterized using currently available metabolic parameters for permethrin. The human permethrin model was first evaluated with data from published human studies. Then, it was used to estimate urinary metabolite concentrations of permethrin in flight attendants who worked in aircrafts, which underwent residual and pre-flight spray treatments. The human model was also applied to analyze the toxicokinetics following permethrin exposures attributable to other aircraft disinsection scenarios. Predicted levels of urinary 3-phenoxybenzoic acid (3-PBA), a metabolite of permethrin, following residual disinsection treatment were comparable to the measurements made for flight attendants. Simulations showed that the median contributions of the dermal, oral and inhalation routes to permethrin exposure in flight attendants were 83.5%, 16.1% and 0.4% under residual treatment scenario, respectively, and were 5.3%, 5.0% and 89.7% under pre-flight spray scenario, respectively. The PBPK model provides the capability to simulate the toxicokinetic profiles of permethrin, and can be used in the studies on human exposure to permethrin.

The influence of physicochemical properties on the internal dose of trihalomethanes in humans following a controlled showering exposure

Although disinfection of domestic water supply is crucial for protecting public health from waterborne diseases, this process forms potentially harmful by-products, such as trihalomethanes (THMs). We evaluated the influence of physicochemical properties of four THMs (chloroform, bromodichloromethane, dibromochloromethane, and bromoform) on the internal dose after showering. One hundred volunteers showered for 10 min in a controlled setting with fixed water flow, air flow, and temperature. We measured THMs in shower water, shower air, bathroom air, and blood samples collected at various time intervals. The geometric mean (GM) for total THM concentration in shower water was 96.2 μg/l. The GM of total THM in air increased from 5.8 μg/m(3) pre shower to 351 μg/m(3) during showering. Similarly, the GM of total-blood THM concentration increased from 16.5 ng/l pre shower to 299 ng/l at 10 min post shower. THM levels were significantly correlated between different matrices (e.g. dibromochloromethane levels) in water and air (r=0.941); blood and water (r=0.845); and blood and air (r=0.831). The slopes of best-fit lines for THM levels in water vs air and blood vs air increased with increasing partition coefficient of water/air and blood/air. The slope of the correlation plot of THM levels in water vs air decreased in a linear (r=0.995) fashion with increasing Henry's law constant. The physicochemical properties (volatility, partition coefficients, and Henry's law constant) are useful parameters for predicting THM movement between matrices and understanding THM exposure during showering.

Effects of temperature, surfactants and skin location on the dermal penetration of haloacetonitriles and chloral hydrate

Dermal exposure has been recognized as an important contributor to the total internal dose to disinfection-by-products (DBPs) in water. However, the effect of the use of surfactants, water temperature and area of the body exposed to DBPs on their dermal flux has not been characterized and was the focus of the present study using an in-vitro system. The dermal flux of mg/l concentrations of haloacetonitriles and chloral hydrate (CH), important cytotoxic DBPs, increased by approximately 50% to 170% with increasing temperature from 25 °C to 40 °C. The fluxes for the torso and dorsum of the hand were much higher than that of palm and scalp skin. An increase in flux was observed for chloroacetonitrite and dichloroacetonitrile, two less lipophilic HANs, but not for trichloroacetonitrile or CH, with the addition of 2% sodium lauryl sulfate or 2% sodium laureth sulfate, two surfactants commonly used in soaps and shampoos used in showering and bathing. Thus, factors such as temperature, surfactants and skin location affect dermal penetration and should be considered when evaluating dermal absorption.

Application of physiologically based pharmacokinetic models in chemical risk assessment

Post-exposure risk assessment of chemical and environmental stressors is a public health challenge. Linking exposure to health outcomes is a 4-step process: exposure assessment, hazard identification, dose response assessment, and risk characterization. This process is increasingly adopting "in silico" tools such as physiologically based pharmacokinetic (PBPK) models to fine-tune exposure assessments and determine internal doses in target organs/tissues. Many excellent PBPK models have been developed. But most, because of their scientific sophistication, have found limited field application-health assessors rarely use them. Over the years, government agencies, stakeholders/partners, and the scientific community have attempted to use these models or their underlying principles in combination with other practical procedures. During the past two decades, through cooperative agreements and contracts at several research and higher education institutions, ATSDR funded translational research has encouraged the use of various types of models. Such collaborative efforts have led to the development and use of transparent and user-friendly models. The "human PBPK model toolkit" is one such project. While not necessarily state of the art, this toolkit is sufficiently accurate for screening purposes. Highlighted in this paper are some selected examples of environmental and occupational exposure assessments of chemicals and their mixtures.

Translational research to develop a human PBPK models tool kit-volatile organic compounds (VOCs)

Toxicity and exposure evaluations remain the two of the key components of human health assessment. While improvement in exposure assessment relies on a better understanding of human behavior patterns, toxicity assessment still relies to a great extent on animal toxicity testing and human epidemiological studies. Recent advances in computer modeling of the dose-response relationship and distribution of xenobiotics in humans to important target tissues have advanced our abilities to assess toxicity. In particular, physiologically based pharmacokinetic (PBPK) models are among the tools than can enhance toxicity assessment accuracy. Many PBPK models are available to the health assessor, but most are so difficult to use that health assessors rarely use them. To encourage their use these models need to have transparent and user-friendly formats. To this end the Agency for Toxic Substances and Disease Registry (ATSDR) is using translational research to increase PBPK model accessibility, understandability, and use in the site-specific health assessment arena. The agency has initiated development of a human PBPK tool-kit for certain high priority pollutants. The tool kit comprises a series of suitable models. The models are recoded in a single computer simulation language and evaluated for use by health assessors. While not necessarily being state-of-the-art code for each chemical, the models will be sufficiently accurate to use for screening purposes. This article presents a generic, seven-compartment PBPK model for six priority volatile organic compounds (VOCs): benzene (BEN), carbon tetrachloride (CCl(4)), dichloromethane (DCM), perchloroethylene (PCE), trichloroethylene (TCE), and vinyl chloride (VC). Limited comparisons of the generic and original model predictions to published kinetic data were conducted. A goodness of fit was determined by calculating the means of the sum of the squared differences (MSSDs) for simulation vs. experimental kinetic data using the generic and original models. Using simplified solvent exposure assumptions for oral ingestion and inhalation, steady-state blood concentrations of each solvent were simulated for exposures equivalent to the ATSDR Minimal Risk Levels (MRLs). The predicted blood levels were then compared to those reported in the National Health and Nutrition Examination Survey (NHANES). With the notable exception of BEN, simulations of combined oral and inhalation MRLs using our generic VOC model yielded blood concentrations well above those reported for the 95th percentile blood concentrations for the U.S. populations, suggesting no health concerns. When the PBPK tool kit is fully developed, risk assessors will have a readily accessible tool for evaluating human exposure to a variety of environmental pollutants.

A generic, cross-chemical predictive PBTK model with multiple entry routes running as application in MS Excel; design of the model and comparison of predictions with experimental results

AIM

Physiologically based toxicokinetic (PBTK) models are computational tools, which simulate the absorption, distribution, metabolism, and excretion of chemicals. The purpose of this study was to develop a physiologically based pharmacokinetic (PBPK) model with a high level of transparency. The model should be able to predict blood and urine concentrations of environmental chemicals and metabolites, given a certain environmental or occupational exposure scenario.

MODEL

The model refers to a reference human of 70 kg. The partition coefficients of the parent compound and its metabolites (blood:air and tissue:blood partition coefficients of 11 organs) are estimated by means of quantitative structure-property relationship, in which five easily available physicochemical properties of the compound are the independent parameters. The model gives a prediction of the fate of the compound, based on easily available chemical properties; therefore, it can be applied as a generic model applicable to multiple compounds. Three routes of uptake are considered (inhalation, dermal, and/or oral) as well as two built-in exercise levels (at rest and at light work). Dermal uptake is estimated by the use of a dermal diffusion-based module that considers dermal deposition rate and duration of deposition. Moreover, evaporation during skin contact is fully accounted for and related to the volatility of the substance. Saturable metabolism according to Michaelis-Menten kinetics can be modelled in any of 11 organs/tissues or in liver only. Renal tubular resorption is based on a built-in algorithm, dependent on the (log) octanol:water partition coefficient. Enterohepatic circulation is optional at a user-defined rate. The generic PBTK model is available as a spreadsheet application in MS Excel. The differential equations of the model are programmed in Visual Basic. Output is presented as numerical listing over time in tabular form and in graphs. The MS Excel application of the PBTK model is available as freeware.

EXPERIMENTAL

The accuracy of the model prediction is illustrated by simulating experimental observations. Published experimental inhalation and dermal exposure studies on a series of different chemicals (pyrene, N-methyl-pyrrolidone, methyl-tert-butylether, heptane, 2-butoxyethanol, and ethanol) were selected to compare the observed data with the model-simulated data. The examples show that the model-predicted concentrations in blood and/or urine after inhalation and/or transdermal uptake have an accuracy of within an order of magnitude.

CONCLUSIONS

It is advocated that this PBTK model, called IndusChemFate, is suitable for 'first tier assessments' and for early explorations of the fate of chemicals and/or metabolites in the human body. The availability of a simple model with a minimum burden of input information on the parent compound and its metabolites might be a stimulation to apply PBTK modelling more often in the field of biomonitoring and exposure science.

Development of a human Physiologically Based Pharmacokinetic (PBPK) Toolkit for environmental pollutants

Physiologically Based Pharmacokinetic (PBPK) models can be used to determine the internal dose and strengthen exposure assessment. Many PBPK models are available, but they are not easily accessible for field use. The Agency for Toxic Substances and Disease Registry (ATSDR) has conducted translational research to develop a human PBPK model toolkit by recoding published PBPK models. This toolkit, when fully developed, will provide a platform that consists of a series of priority PBPK models of environmental pollutants. Presented here is work on recoded PBPK models for volatile organic compounds (VOCs) and metals. Good agreement was generally obtained between the original and the recoded models. This toolkit will be available for ATSDR scientists and public health assessors to perform simulations of exposures from contaminated environmental media at sites of concern and to help interpret biomonitoring data. It can be used as screening tools that can provide useful information for the protection of the public.

A Bayesian population PBPK model for multiroute chloroform exposure

A Bayesian hierarchical model was developed to estimate the parameters in a physiologically based pharmacokinetic (PBPK) model for chloroform using prior information and biomarker data from different exposure pathways. In particular, the model provides a quantitative description of the changes in physiological parameters associated with hot-water bath and showering scenarios. Through Bayesian inference, uncertainties in the PBPK parameters were reduced from the prior distributions. Prediction of biomarker data with the calibrated PBPK model was improved by the calibration. The posterior results indicate that blood flow rates varied under two different exposure scenarios, with a two-fold increase of the skin's blood flow rate predicted in the hot-bath scenario. This result highlights the importance of considering scenario-specific parameters in PBPK modeling. To demonstrate the application of a probability approach in toxicological assessment, results from the posterior distributions from this calibrated model were used to predict target tissue dose based on the rate of chloroform metabolized in liver. This study demonstrates the use of the Bayesian approach to optimize PBPK model parameters for typical household exposure scenarios.

Approaches for applications of physiologically based pharmacokinetic models in risk assessment

Physiologically based pharmacokinetic (PBPK) models are particularly useful for simulating exposures to environmental toxicants for which, unlike pharmaceuticals, there is often little or no human data available to estimate the internal dose of a putative toxic moiety in a target tissue or an appropriate surrogate. This article reviews the current state of knowledge and approaches for application of PBPK models in the process of deriving reference dose, reference concentration, and cancer risk estimates. Examples drawn from previous U.S. Environmental Protection Agency (EPA) risk assessments and human health risk assessments in peer-reviewed literature illustrate the ways and means of using PBPK models to quantify the pharmacokinetic component of the interspecies and intraspecies uncertainty factors as well as to conduct route to route, high dose to low dose and duration extrapolations. The choice of the appropriate dose metric is key to the use of the PBPK models for the various applications in risk assessment. Issues related to whether uncertainty factors are most appropriately applied before or after derivation of human equivalent dose (or concentration) continue to be explored. Scientific progress in the understanding of life stage and genetic differences in dosimetry and their impacts on variability in susceptibility, as well as ongoing development of analytical methods to characterize uncertainty in PBPK models, will make their use in risk assessment increasingly likely. As such, it is anticipated that when PBPK models are used to express adverse tissue responses in terms of the internal target tissue dose of the toxic moiety rather than the external concentration, the scientific basis of, and confidence in, risk assessments will be enhanced.

Computational toxicology of chloroform: reverse dosimetry using Bayesian inference, Markov chain Monte Carlo simulation, and human biomonitoring data

BACKGROUND

One problem of interpreting population-based biomonitoring data is the reconstruction of corresponding external exposure in cases where no such data are available.

OBJECTIVES

We demonstrate the use of a computational framework that integrates physiologically based pharmacokinetic (PBPK) modeling, Bayesian inference, and Markov chain Monte Carlo simulation to obtain a population estimate of environmental chloroform source concentrations consistent with human biomonitoring data. The biomonitoring data consist of chloroform blood concentrations measured as part of the Third National Health and Nutrition Examination Survey (NHANES III), and for which no corresponding exposure data were collected.

METHODS

We used a combined PBPK and shower exposure model to consider several routes and sources of exposure: ingestion of tap water, inhalation of ambient household air, and inhalation and dermal absorption while showering. We determined posterior distributions for chloroform concentration in tap water and ambient household air using U.S. Environmental Protection Agency Total Exposure Assessment Methodology (TEAM) data as prior distributions for the Bayesian analysis.

RESULTS

Posterior distributions for exposure indicate that 95% of the population represented by the NHANES III data had likely chloroform exposures < or = 67 microg/L [corrected] in tap water and < or = 0.02 microg/L in ambient household air.

CONCLUSIONS

Our results demonstrate the application of computer simulation to aid in the interpretation of human biomonitoring data in the context of the exposure-health evaluation-risk assessment continuum. These results should be considered as a demonstration of the method and can be improved with the addition of more detailed data.

Approaches for evaluating the relevance of multiroute exposures in establishing guideline values for drinking water contaminants

In establishing the guideline values for chemical contaminants in drinking water, the contribution of inhalation and dermal routes associated with showering/bathing needs to be evaluated. The present article reviews the current approaches available for evaluating the importance of inhalation and dermal routes of exposure to drinking water contaminants (DWCs) and integrates them within a 2-tier approach. Accordingly, tier 1 would evaluate whether the dermal or inhalation route is likely to contribute to at least 10% of the dose received from ingestion of drinking water (i.e., 0.15 L-equivalent per day based on the daily water intake rate of 1.5 L/day typically used in Health Canada assessments). Based on the route-specific exposure parameters (i.e., area of skin exposed, effective skin permeability coefficient [K(p)], and air to water concentration ratio during use conditions [F(air-water)], breathing rate, duration of contact, and fraction absorbed), it was determined that for DWCs with K(p) less than 0.024 cm/hr and F(air - water) less than 0.0063, the dermal and inhalation routes during showering or bathing are unlikely to contribute significantly to the total dose. For DWCs with K(p) value equal to or greater than 0.025 cm/hr, dermal notation is implied, and as such, tier 2 calculation of L-equivalent associated with dermal exposure needs to be performed. Similarly, for DWCs with F(air-water) greater than 0.00063, inhalation notation is implied, and detailed evaluation of the L-equivalent associated with inhalation exposure (i.e., tier 2) is suggested. In general, data from human volunteer studies, observational measurements, and targeted modeling studies are useful for deriving L-equivalents, reflective of the magnitude of dose received via dermal and inhalation routes relative to the oral route. However, in resource-limited situations, these approaches can be integrated within a 2-tier approach for prioritizing and providing quantitative evaluations of the relevance of dermal and inhalation routes for developing exposure guidelines for DWCs.

Reverse dosimetry: interpreting trihalomethanes biomonitoring data using physiologically based pharmacokinetic modeling

Biomonitoring data provide evidence of exposure of environmental chemicals but are not, by themselves, direct measures of exposure. To use biomonitoring data in understanding exposure, physiologically based pharmacokinetic (PBPK) modeling can be used in a reverse dosimetry approach to assess a distribution of exposures possibly associated with specific blood or urine levels of compounds. Reverse dosimetry integrates PBPK modeling with exposure pattern characterization, Monte Carlo analysis, and statistical tools to estimate a distribution of exposures that are consistent with biomonitoring data in a population. The present study used an existing PBPK model for chloroform as a generic framework to develop PBPK models for other trihalomethanes (THMs). Using Monte Carlo sampling techniques, probabilistic information about pharmacokinetics and exposure patterns was included to estimate distributions of THMs concentrations in blood in relation to various exposure patterns in a diverse population. In addition, the possibility of inhibition of hepatic metabolism among THMs was evaluated under the scenarios of household exposure. These studies demonstrated how PBPK modeling can be used as a tool to estimate a population distribution of exposures that could have resulted in particular biomonitoring results. When toxicity level is known, this tool can also be used to estimate proportion of population above levels associated with health risk.

Assessment of exposure of workers and swimmers to trihalomethanes in an indoor swimming pool

A simultaneous study on workers' and swimmers' exposure to trihalomethanes (THMs) in an indoor swimming pool has been carried out by analyzing urine samples using the headspace and gas chromatography-mass spectrometry technique. The subjects of this study were male and female workers of an indoor swimming pool as well as swimmers regularly attending the pool. The results reported show that only chloroform and bromodichloromethane were detected in the urine of those people exposed, which can be used as a specific index of exposure to these compounds. THM uptake of swimmers after 1 h of swimming was higher than that of workers after a 4 h work shift since THM levels in the workers' urine were associated only with inhalation, while levels in swimmers' urine were mainly associated with dermal absorption, apart from inhalation and occasional ingestion, as well as increased uptake due to the physical stress (swimming). The kinetics of THM excretion in the urine of the participants exposed has been calculated after termination of the exposure to select the sampling time and determine the elimination process. An interval of 15 min after exposure was selected as the sampling time, and the absorbed dosage was eliminated by 2 h after exposure. A good correlation between THM concentrations found in the swimming pool water and the urinary THM concentrations of the people affected after exposure has also been obtained.

Development of physiologically based toxicokinetic models for improving the human indoor exposure assessment to water contaminants: trichloroethylene and trihalomethanes

Generally, ingestion is the only route of exposure that is considered in the risk assessment of drinking water contaminants. However, it is well known that a number of these contaminants are volatile and lipophilic and therefore highly susceptible to being absorbed through other routes, mainly inhalation and dermal. The objective of this study was to develop physiologically based human toxicokinetic (PBTK) models for trihalomethanes (THM) and trichloroethylene (TCE) that will facilitate (1) the estimation of internal exposure to these chemicals for various multimedia indoor exposure scenarios, and (2) consideration of the impact of biological variability in the estimation of internal doses. Five PBTK models describing absorption through ingestion, inhalation and skin were developed for these contaminants. Their concentrations in ambient air were estimated from their respective tap water concentrations and their physicochemical characteristics. Algebraic descriptions of the physiological parameters, varying as a function of age, gender and diverse anthropometric parameters, allow the prediction of the influence of interindividual variations on absorbed dose and internal dosimetry. Simulations for various scenarios were done for a typical human (i.e., 70 kg, 1.7 m) as well as for humans of both genders varying in age from 1 to 90 years. Simulations show that ingestion contributes to less than 50% of the total absorbed dose or metabolized dose for all chemicals. This contribution to internal dosimetry, such as maximal venous blood concentrations (Cmax) and the area under the venous blood concentration time curve (AUC), decreases markedly (e.g., as low as 0.9% of Cmax for bromodichloromethane). The importance of this contribution varies mainly as a function of shower duration. Moreover, model simulations indicate that multimedia exposure is more elevated in children than adults (i.e., up to 200% of the adult internal dose). The models developed in this study allow characterization of the influence of the different routes of exposure and an improved estimation of the realistic multimedia exposure to volatile organic chemicals present in drinking water. Hence, such models will greatly improve health risk assessment for these chemicals.

Spatial distribution of cutaneous microvasculature and local drug clearance after drug application on the skin

We analysed quantitatively blood microvessels distribution in normal skin. We conclude that the segmental area of blood vessels peaks approximately 0.1 mm below the skin surface, where the upper cutaneous blood plexus resides. Total blood vessels area then decreases quasi-exponentially to a depth of approx. -0.75 mm, with a decay length of approximately 0.1 mm, which is site and skin condition dependent, but at greater depths the decrease is approx. 6-times less steep. The corresponding permeability sink exhibits a similar, but superficially steeper, depth-profile. The lateral localisation of superficial blood vessels is such that ensures maximum diffusion from and into the capillaries, which affects transdermal drug delivery: each hairpin-like loop is in the centre of a papilla that corresponds to a cluster of corneocytes surrounded by main diffusion pathways. The aggregate area of blood vessels in the skin is >or=2.5-fold greater than total organ surface area under normal physiological conditions. The molecules diffusing through the skin barrier are thus largely cleared in outermost 20% of the organ, which may create a drug concentration maximum in the dermis, if clearance increases significantly with time. Skin microdialysis data are therefore extremely sensitive to cutaneous blood flow (distribution) and sampling. Skin microvasculature and its distribution must consequently be considered in all topical or transdermal drug transport studies, for example, by including suitably formulated clearance term into generalised diffusion equation.

Use of a physiologically based pharmacokinetic model to identify exposures consistent with human biomonitoring data for chloroform

Biomonitoring data provide evidence of human exposure to environmental chemicals by quantifying the chemical or its metabolite in a biological matrix. To better understand the correlation between biomonitoring data and environmental exposure, physiologically based pharmacokinetic (PBPK) modeling can be of use. The objective of this study was to use a combined PBPK model with an exposure model for showering to estimate the intake concentrations of chloroform based on measured blood and exhaled breath concentrations of chloroform. First, the predictive ability of the combined model was evaluated with three published studies describing exhaled breath and blood concentrations in people exposed to chloroform under controlled showering events. Following that, a plausible exposure regimen was defined combining inhalation, ingestion, and dermal exposures associated with residential use of water containing typical concentrations of chloroform to simulate blood and exhaled breath concentrations of chloroform. Simulation results showed that inhalation and dermal exposure could contribute substantially to total chloroform exposure. Next, sensitivity analysis and Monte Carlo analysis were performed to investigate the sources of variability in model output. The variability in exposure conditions (e.g., shower duration) was shown to contribute more than the variability in pharmacokinetics (e.g., body weight) to the predicted variability in blood and exhaled breath concentrations of chloroform. Lastly, the model was used in a reverse dosimetry approach to estimate distributions of exposure consistent with concentrations of chloroform measured in human blood and exhaled breath.

Dermal uptake of chloroform and haloketones during bathing

Dermal contact with some organic disinfection by-products (DBPs) such as trihalomethanes in chlorinated drinking water has been established to be an important exposure route. We evaluated dermal absorption of two haloketones (1,1-dichloropropanone and 1,1,1-trichloropropanone) and chloroform while bathing, by collecting and analyzing time profiles of expired breath samples of six human subjects during and following a 30-min bath. The DBP concentrations in breath increased towards a maximum concentration during bathing. The maximum haloketone breath concentration during dermal exposure ranged from 0.1 to 0.9 microg / m(3), which was approximately two orders of magnitude lower than the maximum chloroform breath concentration during exposure. Based on a one-compartment model, the in vivo permeability of chloroform, 1,1-dichloropropanone, and 1,1,1-trichloropropanone were approximated to be 0.015, 7.5 x 10(- 4), and 4.5 x 10(- 4) cm / h, respectively. Thus, haloketones are much less permeable across human skin under normal bathing conditions than is chloroform. These findings will be useful for future assessment of total human exposure and consequent health risk of these DBPs.

Risk assessment of chemicals and pharmaceuticals in the pediatric population: a workshop report

ATSDR and RIVM organized an Expert Panel Workshop on the Differences Between Children and Adults and Their Relevance to Risk Assessment. The workshop was held in June 2003, in Brussels, Belgium. The purpose of the workshop was to identify data gaps in current scientific knowledge related to children's health and to recognize areas of mutual interest that would serve as the basis for upcoming ATSDR/RIVM cooperative projects. The aim for both agencies is a better understanding of the issues related to children's health, and the improvement of scientifically based (chemical) risk assessment in children. Topics discussed included clinical trials/toxicity studies, testing in juvenile animals, PBPK modeling in children, and children's risk assessment.

Benzene and naphthalene in air and breath as indicators of exposure to jet fuel

AIMS

To estimate exposures to benzene and naphthalene among military personnel working with jet fuel (JP-8) and to determine whether naphthalene might serve as a surrogate for JP-8 in studies of health effects.

METHODS

Benzene and naphthalene were measured in air and breath of 326 personnel in the US Air Force, who had been assigned a priori into low, moderate, and high exposure categories for JP-8.

RESULTS

Median air concentrations for persons in the low, moderate, and high exposure categories were 3.1, 7.4, and 252 microg benzene/m3 air, 4.6, 9.0, and 11.4 microg benzene/m3 breath, 1.9, 10.3, and 485 microg naphthalene/m3 air, and 0.73, 0.93, and 1.83 microg naphthalene/m3 breath, respectively. In the moderate and high exposure categories, 5% and 15% of the benzene air concentrations, respectively, were above the 2002 threshold limit value (TLV) of 1.6 mg/m3. Multiple regression analyses of air and breath levels revealed prominent background sources of benzene exposure, including cigarette smoke. However, naphthalene exposure was not unduly influenced by sources other than JP-8. Among heavily exposed workers, dermal contact with JP-8 contributed to air and breath concentrations along with several physical and environmental factors.

CONCLUSIONS

Personnel having regular contact with JP-8 are occasionally exposed to benzene at levels above the current TLV. Among heavily exposed workers, uptake of JP-8 components occurs via both inhalation and dermal contact. Naphthalene in air and breath can serve as useful measures of exposure to JP-8 and uptake of fuel components in the body.

Whole body pharmacokinetic models

The aim of the current review is to summarise the present status of physiologically based pharmacokinetic (PBPK) modelling and its applications in drug research, and thus serve as a reference point to people interested in the methodology. The review is structured into three major sections. The first discusses the existing methodologies and techniques of PBPK model development. The second describes some of the most interesting PBPK model implementations published. The final section is devoted to a discussion of the current limitations and the possible future developments of the PBPK modelling approach. The current review is focused on papers dealing with the pharmacokinetics and/or toxicokinetics of medicinal compounds; references discussing PBPK models of environmental compounds are mentioned only if they represent considerable methodological developments or reveal interesting interpretations and/or applications.The major conclusion of the review is that, despite its significant potential, PBPK modelling has not seen the development and implementation it deserves, especially in the drug discovery, research and development processes. The main reason for this is that the successful development and implementation of a PBPK model is seen to require the investment of significant experience, effort, time and resources. Yet, a substantial body of PBPK-related research has been accumulated that can facilitate the PBPK modelling and implementation process. What is probably lagging behind is the expertise component, where the demand for appropriately qualified staff far outreaches availability.

Adduction of the chloroform metabolite phosgene to lysine residues of human histone H2B

Human exposure to trihalomethanes such as chloroform has been associated with both cancer and reproductive toxicity. While there is little evidence for chloroform mutagenicity or DNA adduct formation, in vivo studies in rats have demonstrated adduction to histones and other nuclear proteins. Histones play a key role in controlling DNA expression particularly through the acetylation of lysine residues in their N-termini. Therefore, we studied the reaction of phosgene, the major active metabolite of chloroform, with the N-terminus of human histone H2B (Hpep, Pro-Glu-Pro-Ala-Lys-Ser-Ala-Pro-Ala-Pro-Lys-Lys-Gly-Ser-Lys-Lys-Ala-Val-Thr-Lys-Ala-Gln-Lys) in a model chemical system. The aim of this study was to assess whether phosgene is able to form irreversible adducts with this peptide and to investigate which residues are most susceptible. Hpep was reacted with a range of phosgene concentrations (0.03-36 mM) at 37 degrees C, pH 7.4. The products of these reactions, analyzed by matrix-assisted laser desorption ionization MS, showed that up to three CO moieties could be adducted to the peptide. The singly and doubly adducted peptides were purified by HPLC and then hydrolyzed with trypsin to produce a series of fragments that were analyzed by HPLC-MS. The tryptic products showed that adduction occurred principally at lysine residues, and that all seven lysine residues of the peptide were subject to adduction. Collision-induced dissociation analysis using ion trap MS-MS of the tryptic fragment [Pro-Glu-Pro-Ala-Lys-Ser-Ala-Pro-Ala-Pro-Lys + CO] and of the full-length singly adducted peptide supported the role of lysine residues in adduction; the data also indicated that the N-terminal proline and the serine residues are susceptible. Addition of glutathione to the reaction mixture only partially attenuated adduct formation and allowed production of another adducted species, i.e., Hpep-CO-glutathione. The occurrence of such reactions to the N-termini of histones, if confirmed by in vivo studies, could help to explain the mechanism of chloroform carcinogenicity.

Evaluation of the dermal absorption of aqueous toluene in F344 rats using real-time breath analysis and physiologically based pharmacokinetic modeling

Toluene is a ubiquitous chemical that is commonly used for its solvent properties in industry and manufacturing, and is a component of many paint products. Because of its widespread use, there is potential for both occupational and nonoccupational dermal exposure to toluene. To understand the significance of these exposures, the dermal bioavailability of toluene was assessed in F344 male rats using a combination of real-time exhaled breath analysis and physiologically based pharmacokinetic (PBPK) modeling. Animals were exposed to toluene at 0.5 or 0.2 mg/ml aqueous concentration (0.05% or 0.02%) using a 2.5-cm-diameter occluded glass patch system attached to a clipper-shaved area on the back of the rat. Immediately following exposure, individual animals were placed in a glass off-gassing chamber and exhaled breath was monitored as chamber concentration in real time using an ion-trap mass spectrometer (MS/MS). The real-time exhaled breath profile clearly demonstrated the rapid absorption of toluene, with peak chamber concentrations observed within 1 h from the start of exposure. The PBPK model describing the exposure and off-gassing chamber was used to estimate a dermal permeability coefficient (K(p)) to describe each set of exhaled breath data. Regardless of exposure level, a single K(p) value of 0.074 +/- 0.005 cm/h provided a good fit to all data sets. These rat studies using aqueous toluene will form the basis for comparing the dermal bioavailability of toluene in various paint products and may ultimately aid in understanding human health risk under a variety of exposure scenarios.

Route-of-entry and brain tissue partition coefficients for common superfund contaminants

Various organic solvents may be encountered in contaminated water supplies at U.S. Environmental Protection Agency-designated Superfund sites. Human exposure to these environmental contaminants may occur by oral, dermal, or inhalation routes. The estimation of human health risk associated with exposure to these solvents can be improved through the use of physiologically based pharmacokinetic (PBPK) models to describe the absorption, tissue distribution, metabolism, and elimination of the compounds following any route of exposure. However, development of these PBPK models requires information on the relative solubility, or partition coefficient, of each compound in blood and various tissues. A number of investigators have provided partition coefficient information on different tissues in various species; however, the data for route of entry organs (i.e., skin, lung, stomach) and brain tissue are not complete. Therefore, the objective of this work was to replicate partition coefficient studies for several commonly encountered environmental contaminants using an in vitro gas-phase vial equilibration technique and to include tissues to evaluate brain, lung, stomach, and skin. A comparison of the partition coefficient values determined here with values reported in the literature, where available, showed good agreement in nearly all cases. An additional study was conducted to compare the liver-to-air partition coefficient values for toluene, benzene, and o-xylene introduced as single chemicals to partition coefficient values determined with the chemicals introduced as a mixture of all three compounds. The similarities of the resulting values suggest that both labor and laboratory resources may be reduced when partition coefficients are determined as chemical mixtures.

Methods for assessing risks of dermal exposures in the workplace

The skin as a route of entry for toxic chemicals has caused increasing concern over the last decade. The assessment of systemic hazards from dermal exposures has evolved over time, often limited by the amount of experimental data available. The result is that there are many methods being used to assess safety of chemicals in the workplace. The process of assessing hazards of skin contact includes estimating the amount of substance that may end up on the skin and estimating the amount that might reach internal organs. Most times, toxicology studies by the dermal route are not available and extrapolations from other exposure routes are necessary. The hazards of particular chemicals can be expressed as "skin notations", actual exposure levels, or safe exposure times. Characterizing the risk of a specific procedure in the workplace involves determining the ratio of exposure standards to an expected exposure. The purpose of this review is to address each of the steps in the process and describe the assumptions that are part of the process. Methods are compared by describing their strengths and weaknesses. Recommendations for research in this area are also included.

PBPK modeling of the percutaneous absorption of perchloroethylene from a soil matrix in rats and humans

Perchloroethylene (PCE) is a widely used volatile organic chemical. Exposures to PCE are primarily through inhalation and dermal contact. The dermal absorption of PCE from a soil matrix was compared in rats and humans using real-time MS/MS exhaled breath technology and physiologically based pharmacokinetic (PBPK) modeling. Studies with rats were performed to compare the effects of loading volume, concentration, and occlusion. In rats, the percutaneous permeability coefficient (K(P)) for PCE was 0.102 +/- 0.017, and was independent of loading volume, concentration, or occlusion. Exhaled breath concentrations peaked within 1 h in nonoccluded exposures, but were maintained over the 5 h exposure period when the system was occluded. Three human volunteers submerged a hand in a container of PCE-laden soil for 2 h and their exhaled breath was continually monitored during and for 2.5 h following exposure. The absorption and elimination kinetics of PCE were slower in these subjects than initially predicted based upon the PBPK model developed from rat dermal kinetic data. The resulting K(P) for humans was over 100-fold lower than for the rat utilizing a single, well-stirred dermal compartment. Therefore, two additional PBPK skin compartment models were evaluated: a parallel model to simulate follicular uptake and a layered model to portray a stratum corneum barrier. The parallel dual dermal compartment model was not capable of describing the exhaled breath kinetics, whereas the layered model substantially improved the fit of the model to the complex kinetics of dermal absorption through the hand. In real-world situations, percutaneous absorption of PCE is likely to be minimal.

Skin absorption and human risk assessment

A common practice is to assume that percutaneous absorption does not significantly contribute to total bioavailability and therefore, absorption through other routes is more important to human risk assessment. The skin can represent a significant barrier to absorption, but some substances are absorbed to a significant extent. Since there is a potential for percutaneous penetration that is not consistent between species or substances, the assessment of the potential contribution of total body burden from dermal exposures should be considered. This review briefly discusses some theories, practices, and factors that affect percutaneous absorption with an emphasis on how percutaneous absorption evaluations apply to human risk assessment.

A real-time in-vivo method for studying the percutaneous absorption of volatile chemicals

Realistic estimates of percutaneous absorption following exposures to solvents in the workplace, or through contaminated soil and water, are critical to understanding human health risks. A method was developed to determine dermal uptake of solvents under non-steady-state conditions using real-time breath analysis in rats, monkeys, and humans. The exhaled breath was analyzed using an ion-trap mass spectrometer, which can quantitate chemicals in the exhaled breath stream in the 1-5 ppb range. The resulting data were evaluated using physiologically-based pharmacokinetic (PBPK) models to estimate dermal permeability constants (Kp) under various exposure conditions. The effects of exposure matrix (soil versus water), occlusion versus non-occlusion, and species differences on the absorption of methyl chloroform, trichloroethylene, and benzene were compared. Exposure concentrations were analyzed before and at 0.5-hour intervals throughout the exposures. The percentage of each chemical absorbed and the corresponding Kp were estimated by optimization of the PBPK model to the medium concentration and the exhaled-breath data. The method was found to be sufficiently sensitive for animal and human dermal studies at low exposure concentrations over small body surface areas, for short periods, using non-steady-state exposure conditions.