Dietary and Inhalation Intake of Lead and Estimation of Blood Lead Levels in Adults and Children in Kanpur, India

This research was initiated to study lead levels in various food items in the city of Kanpur, India, to assess the dietary intake of lead and to estimate blood lead (PbB) levels, a biomarker of lead toxicity. For this purpose, sampling of food products, laboratory analysis, and computational exercises were undertaken. Specifically, six food groups (leafy vegetables, nonleafy vegetables, fruits, pulses, cereals, and milk), drinking water, and lead air concentration were considered for estimating lead intake. Results indicated highest lead content in leafy vegetables followed by pulses. Fruits showed low lead content and drinking water lead levels were always within tolerable limits. It was estimated that average daily lead intake through diet was about 114 microg/day for adults and 50 microg/day in children; tolerable limit is 250 microg/day for adults and 90 microg/day for children. The estimated lead intakes were translated into the resultant PbB concentrations for children and adults using a physiologically-based pharmacokinetic (PBPK) model. Monte Carlo simulation of PbB level variations for adults showed that probability of exceeding the tolerable limit of PbB (i.e.,10 microg/dL) was 0.062 for the pre-unleaded and 0.000328 for the post-unleaded gasoline period. The probability of exceeding tolerable limits in PbB level was reduced by a factor of 189 in the post-unleaded scenario. The study also suggested that in spite of the introduction of unleaded gasoline, children continue to be at a high risk (probability of exceeding 10 microg/dL = 0.39) because of a high intake of lead per unit body weight.     

Uncertainties in Pharmacokinetic Modeling for Perchloroethylene: II. Comparison of Model Predictions with Data for a Variety of Different Parameters

In this paper we compare expectations derived from 10 different human physiologically based pharmacokinetic models for perchloroethylene with data on absorption via inhalation, and concentrations in alveolar air and venous blood. Our most interesting finding is that essentially all of the models show a time pattern of departures of predictions of air and blood levels relative to experimental data that might be corrected by more sophisticated model structures incorporating either (a) heterogeneity of the fat compartment (with respect to either perfusion or partition coefficients or both) or (b) intertissue diffusion of perchloroethylene between the fat and muscle/VRG groups. Similar types of corrections have recently been proposed to reduce analogous anomalies in the fits of pharmacokinetic models to the data for several volatile anesthetics.^(17-20) A second finding is that models incorporating resting values for alveolar ventilation in the region of 5.4 L/min seemed to be most compatible with the most reliable set of perchloroethylene uptake data.     

Use of Markov Chain Monte Carlo Analysis with a Physiologically-Based Pharmacokinetic Model of Methylmercury to Estimate Exposures in U.S. Women of Childbearing Age

A Bayesian approach, implemented using Markov Chain Monte Carlo (MCMC) analysis, was applied with a physiologically‐based pharmacokinetic (PBPK) model of methylmercury (MeHg) to evaluate the variability of MeHg exposure in women of childbearing age in the U.S. population. The analysis made use of the newly available National Health and Nutrition Survey (NHANES) blood and hair mercury concentration data for women of age 16–49 years (sample size, 1,582). Bayesian analysis was performed to estimate the population variability in MeHg exposure (daily ingestion rate) implied by the variation in blood and hair concentrations of mercury in the NHANES database. The measured variability in the NHANES blood and hair data represents the result of a process that includes interindividual variation in exposure to MeHg and interindividual variation in the pharmacokinetics (distribution, clearance) of MeHg. The PBPK model includes a number of pharmacokinetic parameters (e.g., tissue volumes, partition coefficients, rate constants for metabolism and elimination) that can vary from individual to individual within the subpopulation of interest. Using MCMC analysis, it was possible to combine prior distributions of the PBPK model parameters with the NHANES blood and hair data, as well as with kinetic data from controlled human exposures to MeHg, to derive posterior distributions that refine the estimates of both the population exposure distribution and the pharmacokinetic parameters. In general, based on the populations surveyed by NHANES, the results of the MCMC analysis indicate that a small fraction, less than 1%, of the U.S. population of women of childbearing age may have mercury exposures greater than the EPA RfD for MeHg of 0.1 μg/kgg/day, and that there are few, if any, exposures greater than the ATSDR MRL of 0.3 μgg/kgg/day. The analysis also indicates that typical exposures may be greater than previously estimated from food consumption surveys, but that the variability in exposure within the population of U.S. women of childbearing age may be less than previously assumed.     

Development of a Screening Approach to Interpret Human Biomonitoring Data on Volatile Organic Compounds: Reverse Dosimetry on Biomonitoring Data for Trichloroethylene

A screening approach is developed for volatile organic compounds (VOCs) to estimate exposures that correspond to levels measured in fluids and/or tissues in human biomonitoring studies. The approach makes use of a generic physiologically-based pharmacokinetic (PBPK) model coupled with exposure pattern characterization, Monte Carlo analysis, and quantitative structure property relationships (QSPRs). QSPRs are used for VOCs with minimal data to develop chemical-specific parameters needed for the PBPK model. The PBPK model is capable of simulating VOC kinetics following multiple routes of exposure, such as oral exposure via water ingestion and inhalation exposure during shower events. Using published human biomonitoring data of trichloroethylene (TCE), the generic model is evaluated to determine how well it estimates TCE concentrations in blood based on the known drinking water concentrations. In addition, Monte Carlo analysis is conducted to characterize the impact of the following factors: (1) uncertainties in the QSPR-estimated chemical-specific parameters; (2) variability in physiological parameters; and (3) variability in exposure patterns. The results indicate that uncertainty in chemical-specific parameters makes only a minor contribution to the overall variability and uncertainty in the predicted TCE concentrations in blood. The model is used in a reverse dosimetry approach to derive estimates of TCE concentrations in drinking water based on given measurements of TCE in blood, for comparison to the U.S. EPA's Maximum Contaminant Level in drinking water. This example demonstrates how a reverse dosimetry approach can be used to facilitate interpretation of human biomonitoring data in a health risk context by deriving external exposures that are consistent with a biomonitoring data set, thereby permitting comparison with health-based exposure guidelines.     

A Trivariate Distribution for the Height,Weight, and Fat of Adult Men*

Using exploratory data analysis, probability plots, scatterplots, and computer animations to rotate and visualize the data, we fit a trivariate Normal distribution to data for the height, the natural logarithm of body weight, and the body fat for 646 men between the ages of 50 and 80 years as reported by the medical staff of the U.S. Veterans Administration's “Normative Aging Study” in Boston, MA. Although these data do not include any children, women, or young men, the measurements represent the best data that we could find through a 4-year search. We believe that these data are well measured and reliable for men in the specified age range and that these data reveal an interesting statistical pattern for use in probabilistic PBPK models.     

Book Reviews

Books reviewed:
Office International Des Epizooties (World Organization for Animal Health) World Animal Health in 1998, Part 1: Reports on the Animal Health Status and Disease Control Methods and Tables on Incidence of List A Diseases and Part 2: Tables on the Animal Health Status and Disease Control Methods
Peter Sedlmeier. Lawrence Erlbaum Associates, Mahwah Improving Statistical Reasoning: Theoretical Models and Practical Implications
Andrew Ford. Island Press Modeling the Environment. An Introduction to System Dynamics Modeling of Environmental Systems
Jacob I. Bregman. Lewis Publishers, Boca Raton Environmental Impact Statements, Second Edition     

Improving Cancer Dose–Response Characterization by Using Physiologically Based Pharmacokinetic Modeling: An Analysis of Pooled Data for Acrylonitrile-Induced Brain Tumors to Assess Cancer Potency in the Rat

Historically, U.S. regulators have derived cancer slope factors by using applied dose and tumor response data from a single key bioassay or by averaging the cancer slope factors of several key bioassays. Recent changes in U.S. Environmental Protection Agency (EPA) guidelines for cancer risk assessment have acknowledged the value of better use of mechanistic data and better dose–response characterization. However, agency guidelines may benefit from additional considerations presented in this paper. An exploratory study was conducted by using rat brain tumor data for acrylonitrile (AN) to investigate the use of physiologically based pharmacokinetic (PBPK) modeling along with pooling of dose–response data across routes of exposure as a means for improving carcinogen risk assessment methods. In this study, two contrasting assessments were conducted for AN-induced brain tumors in the rat on the basis of (1) the EPA's approach, the dose–response relationship was characterized by using administered dose/concentration for each of the key studies assessed individually; and (2) an analysis of the pooled data, the dose–response relationship was characterized by using PBPK-derived internal dose measures for a combined database of ten bioassays. The cancer potencies predicted for AN by the contrasting assessments are remarkably different (i.e., risk-specific doses differ by as much as two to four orders of magnitude), with the pooled data assessments yielding lower values. This result suggests that current carcinogen risk assessment practices overestimate AN cancer potency. This methodology should be equally applicable to other data-rich chemicals in identifying (1) a useful dose measure, (2) an appropriate dose–response model, (3) an acceptable point of departure, and (4) an appropriate method of extrapolation from the range of observation to the range of prediction when a chemical's mode of action remains uncertain.     

An Adjustment Factor for Mode-of-Action Uncertainty with Dual-Mode Carcinogens: The Case of Naphthalene-Induced Nasal Tumors in Rats

The U.S. Environmental Protection Agency (USEPA) guidelines for cancer risk assessment recognize that some chemical carcinogens may have a site-specific mode of action (MOA) involving mutation and cell-killing-induced hyperplasia. The guidelines recommend that for such dual MOA (DMOA) carcinogens, judgment should be used to compare and assess results using separate "linear" (genotoxic) versus "nonlinear" (nongenotoxic) approaches to low-level risk extrapolation. Because the guidelines allow this only when evidence supports reliable risk extrapolation using a validated mechanistic model, they effectively prevent addressing MOA uncertainty when data do not fully validate such a model but otherwise clearly support a DMOA. An adjustment-factor approach is proposed to address this gap, analogous to reference-dose procedures used for classic toxicity endpoints. By this method, even when a "nonlinear" toxicokinetic model cannot be fully validated, the effect of DMOA uncertainty on low-dose risk can be addressed. Application of the proposed approach was illustrated for the case of risk extrapolation from bioassay data on rat nasal tumors induced by chronic lifetime exposure to naphthalene. Bioassay data, toxicokinetic data, and pharmacokinetic analyses were determined to indicate that naphthalene is almost certainly a DMOA carcinogen. Plausibility bounds on rat-tumor-type-specific DMOA-related uncertainty were obtained using a mechanistic two-stage cancer risk model adapted to reflect the empirical link between genotoxic and cytotoxic effects of the most potent identified genotoxic naphthalene metabolites, 1,2- and 1,4-naphthoquinone. Bound-specific adjustment factors were then used to reduce naphthalene risk estimated by linear extrapolation (under the default genotoxic MOA assumption), to account for the DMOA exhibited by this compound.     

A Simple Toxicokinetic Model Exhibiting Complex Dynamics and Nonlinear Exposure Response

Uncertainty in model predictions of exposure response at low exposures is a problem for risk assessment. A particular interest is the internal concentration of an agent in biological systems as a function of external exposure concentrations. Physiologically based pharmacokinetic (PBPK) models permit estimation of internal exposure concentrations in target tissues but most assume that model parameters are either fixed or instantaneously dose-dependent. Taking into account response times for biological regulatory mechanisms introduces new dynamic behaviors that have implications for low-dose exposure response in chronic exposure. A simple one-compartment simulation model is described in which internal concentrations summed over time exhibit significant nonlinearity and nonmonotonicity in relation to external concentrations due to delayed up- or downregulation of a metabolic pathway. These behaviors could be the mechanistic basis for homeostasis and for some apparent hormetic effects.     

Stochasticity in Physiologically Based Kinetics Models: Implications for Cancer Risk Assessment

In case of low-dose exposure to a substance, its concentration in cells is likely to be stochastic. Assessing the consequences of this stochasticity in toxicological risk assessment requires the coupling of macroscopic dynamics models describing whole-body kinetics with microscopic tools designed to simulate stochasticity. In this article, we propose an approach to approximate stochastic cell concentration of butadiene in the cells of diverse organs. We adapted the dynamics equations of a physiologically based pharmacokinetic (PBPK) model and used a stochastic simulator for the system of equations that we derived. We then coupled kinetics simulations with a deterministic hockey stick model of carcinogenicity. Stochasticity induced substantial modifications relative to dose-response curve, compared with the deterministic situation. In particular, there was nonlinearity in the response and the stochastic apparent threshold was lower than the deterministic one. The approach that we developed could easily be extended to other biological studies to assess the influence of stochasticity at macroscopic scale for compound dynamics at the cell level.