Benchmark Dose Risk Assessment for Formaldehyde Using Airflow Modeling and a Single-Compartment, DNA-Protein Cross-Link Dosimetry Model to Estimate Human Equivalent Doses

Formaldehyde induced squamous-cell carcinomas in the nasal passages of F344 rats in two inhalation bioassays at exposure levels of 6 ppm and above. Increases in rates of cell proliferation were measured by T. M. Monticello and colleagues at exposure levels of 0.7 ppm and above in the same tissues from which tumors arose. A risk assessment for formaldehyde was conducted at the CIIT Centers for Health Research, in collaboration with investigators from Toxicological Excellence in Risk Assessment (TERA) and the U.S. Environmental Protection Agency (U.S. EPA) in 1999. Two methods for dose-response assessment were used: a full biologically based modeling approach and a statistically oriented analysis by benchmark dose (BMD) method. This article presents the later approach, the purpose of which is to combine BMD and pharmacokinetic modeling to estimate human cancer risks from formaldehyde exposure. BMD analysis was used to identify points of departure (exposure levels) for low-dose extrapolation in rats for both tumor and the cell proliferation endpoints. The benchmark concentrations for induced cell proliferation were lower than for tumors. These concentrations were extrapolated to humans using two mechanistic models. One model used computational fluid dynamics (CFD) alone to determine rates of delivery of inhaled formaldehyde to the nasal lining. The second model combined the CFD method with a pharmacokinetic model to predict tissue dose with formaldehyde-induced DNA-protein cross-links (DPX) as a dose metric. Both extrapolation methods gave similar results, and the predicted cancer risk in humans at low exposure levels was found to be similar to that from a risk assessment conducted by the U.S. EPA in 1991. Use of the mechanistically based extrapolation models lends greater certainty to these risk estimates than previous approaches and also identifies the uncertainty in the measured dose-response relationship for cell proliferation at low exposure levels, the dose-response relationship for DPX in monkeys, and the choice between linear and nonlinear methods of extrapolation as key remaining sources of uncertainty.     

Bayesian Estimation of Pharmacokinetic and Pharmacodynamic Parameters in a Mode-of-Action-Based Cancer Risk Assessment for Chloroform

Chloroform is a carcinogen in rodents and its carcinogenicity is secondary to events associated with cytotoxicity and regenerative cell proliferation. In this study, a physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) model that links the processes of chloroform metabolism, reparable cell damage, cell death, and regenerative cellular proliferation was developed to support a new cancer dose-response assessment for chloroform. Model parameters were estimated using Markov Chain Monte Carlo (MCMC) analysis in a two-step approach: (1) metabolism parameters for male and female mice and rats were estimated against available closed chamber gas uptake data; and (2) PD parameters for each of the four rodent groups were estimated from hepatic and renal labeling index data following inhalation exposures. Subsequently, the resulting rodent PD parameters together with literature values for human age-dependent physiological and metabolism parameters were used to scale up the rodent model to a human model. The human model was used to predict exposure conditions under which chloroform-mediated cytolethality is expected to occur in liver and kidney of adults and children. Using the human model, inhalation Reference Concentrations (RfCs) and oral Reference Doses (RfDs) were derived using an uncertainty factor of 10. Based on liver and kidney dose metrics, the respective RfCs were 0.9 and 0.09 ppm; and the respective RfDs were 0.4 and 3 mg/kg/day.     

An Improved Approximation to the Exact Solution of the Two-Stage Clonal Growth Model of Cancer

Multistage clonal growth models are of interest for cancer risk assessment because they can explicitly incorporate data on cell replication. Both approximate and exact formulations of the two stage growth model have been described. The exact solution considers the conditional probability of tumors arising in previously tumor-free animals; the approximate solution estimates total probability of tumor formation. The exact solution is much more computationally intensive when time-dependent cell growth parameters are included. The approximate solution deviates from the exact solution at high incidences and probabilities of tumor. This report describes a computationally tractable,'improved approximation'to the exact solution. Our improved approximation includes a correction term to adjust the unconditional expectation of intermediate cells based on the time history of formation of intermediate cells by mutation of normal cells (recruitment) or by cell division in the intermediate cell population (expansion). The improved approximation provided a much better match to the exact solution than the approximate solution for a wide range of parameter values. The correction term also appears to provide insight into the biological factors that contribute to the variance of the expectation for the number of intermediate cells over time.