Updating the Biologically Based Dose Response Model for the Nasal Carcinogenicity of Inhaled Formaldehyde in the F344 Rat

Chronic inhalation of formaldehyde by F344 rats causes nasal squamous cell carcinoma (SCC). This outcome is well-characterized: including dose-response and time course data for SCC, mechanistic endpoints, and nasal dosimetry. Conolly et al. (Toxicol. Sci. 75, 432-447, 2003) used these resources to develop a biologically based dose-response (BBDR) model for SCC in F344 rats. This model, scaled up to humans, has informed dose-response conclusions reached by several international regulatory agencies. However, USEPA concluded that uncertainties precluded its use for cancer risk assessment. Here we describe an updated BBDR model that addresses uncertainties through refined dosimetry modeling, revised analysis of labeling index data, and an extended dataset where both inhaled (exogenous) and endogenous formaldehyde (exogF, endoF) form DNA adducts. Further, since Conolly et al. (ibid) was published, it has become clear that, when controls from all F344 inhalation bioassays are considered, accounting for over 4,000 rats, at most one nasal SCC occurred. This low spontaneous incidence constrains possible contribution of endoF to the formation of nasal SCC via DNA reactivity. Further, since both exogF and endoF form DNA adducts, this constraint also applies to exogF. The revised BBDR model therefore drives SCC formation through the cytotoxicity of high concentration exogF. An option for direct mutagenicity associated with DNA adducts is retained to allow estimation of an upper bound on adduct mutagenicity consistent with the lack of a spontaneous SCC incidence. These updates represent an iterative refinement of the 2003 model, incorporating new data and insights to reduce identified model uncertainties.

Incorporation of rapid association/dissociation processes in tissues into the monkey and human physiologically based pharmacokinetic models for manganese

In earlier physiologically based pharmacokinetic (PBPK) models for manganese (Mn), the kinetics of transport of Mn into and out of tissues were primarily driven by slow rates of association and dissociation of Mn with tissue binding sites. However, Mn is known to show rapidly reversible binding in tissues. An updated Mn model for primates, following similar work with rats, was developed that included rapid association/dissociation processes with tissue Mn-binding sites, accumulation of free Mn in tissues after saturation of these Mn-binding sites and rapid rates of entry into tissues. This alternative structure successfully described Mn kinetics in tissues in monkeys exposed to Mn via various routes including oral, inhalation, and intraperitoneal, subcutaneous, or intravenous injection and whole-body kinetics and tissue levels in humans. An important contribution of this effort is showing that the extension of the rate constants for binding and cellular uptake established in the monkey were also able to describe kinetic data from humans. With a consistent model structure for monkeys and humans, there is less need to rely on cadaver data and whole-body tracer studies alone to calibrate a human model. The increased biological relevance of the Mn model structure and parameters provides greater confidence in applying the Mn PBPK models to risk assessment. This model is also well-suited to explicitly incorporate emerging information on the role of transporters in tissue disposition, intestinal uptake, and hepatobiliary excretion of Mn.

Relative contributions of endogenous and exogenous formaldehyde to formation of deoxyguanosine monoadducts and DNA-protein crosslink adducts of DNA in rat nasal mucosa

Understanding the dose-response for formaldehyde-induced nasal cancer in rats is complicated by (1) the uneven distribution of inhaled formaldehyde across the interior surface of the nasal cavity and, (2) the presence of endogenous formaldehyde (endoF) in the nasal mucosa. In this work, we used computational fluid dynamics (CFD) modeling to predict flux of inhaled (exogenous) formaldehyde (exogF) from air into tissue at the specific locations where DNA adducts were measured. Experimental work has identified DNA-protein crosslink (DPX) adducts due to exogF and deoxyguanosine (DG) adducts due to both exogF and endoF. These adducts can be considered biomarkers of exposure for effects of endoF and exogF on DNA that may be part of the mechanism of tumor formation. We describe a computational model linking CFD-predicted flux of formaldehyde from air into tissue, and the intracellular production of endoF, with the formation of DPX and DG adducts. We assumed that, like exogF, endoF can produce DPX. The model accurately reproduces exogDPX, exogDG, and endoDG data after inhalation from 0.7 to 15 ppm. The dose-dependent concentrations of exogDPX and exogDG are predicted to exceed the concentrations of their endogenous counterparts at about 2 and 6 ppm exogF, respectively. At all concentrations examined, the concentrations of endoDPX and exogDPX were predicted to be at least 10-fold higher than that of their DG counterparts. The modeled dose-dependent concentrations of these adducts are suitable to be used together with data on the dose-dependence of cell proliferation to conduct quantitative modeling of formaldehyde-induced rat nasal carcinogenicity.

Time-dependent genomic response in primary human uroepithelial cells exposed to arsenite for up to 60 days

Evidence from both in vivo and in vitro studies suggests that gene expression changes from long-term exposure to arsenite evolve markedly over time, including reversals in the direction of expression change in key regulatory genes. In this study, human uroepithelial cells from the ureter segments of 4 kidney-donors were continuously treated in culture with arsenite at concentrations of 0.1 or 1 μM for 60 days. Gene expression at 10, 20, 30, 40, and 60 days was determined using Affymetrix human genome microarrays and signal pathway analysis was performed using GeneGo Metacore. Arsenic treated cells continued to proliferate for the full 60-day period, whereas untreated cells ceased proliferating after approximately 30 days. A peak in the number of gene changes in the treated cells compared to untreated controls was observed between 30 and 40 days of exposure, with substantially fewer changes at 10 and 60 days, suggesting remodeling of the cells over time. Consistent with this possibility, the direction of expression change for a number of key genes was reversed between 20 and 30 days, including CFOS and MDM2. While the progression of gene changes was different for each subject, a common pattern was observed in arsenic treated cells over time, with early upregulation of oxidative stress responses (HMOX1, NQ01, TXN, TXNRD1) and down-regulation of immune/inflammatory responses (IKKα). At around 30 days, there was a transition to increased inflammatory and proliferative signaling (AKT, CFOS), evidence of epithelial-to-mesenchymal transition (EMT), and alterations in DNA damage responses (MDM2, ATM). A common element in the changing response of cells to arsenite over time appears to involve up-regulation of MDM2 by inflammatory signaling (through AP-1 and NF-κB), leading to inhibition of P53 function.

A Kinetic Analysis of DNA-Deoxy Guanine Adducts in the Nasal Epithelium Produced by Inhaled Formaldehyde in Rats-Assessing Contributions to Adduct Production From Both Endogenous and Exogenous Sources of Formaldehyde

Although formaldehyde is a normal constituent of tissues, lifetime inhalation exposures at 6 h/day, 5 days/week at concentrations ≥6 ppm caused a nonlinear increase in nasal tumors in rats with incidence reaching close to 50% at 15 ppm. Studies with heavy isotope labeled [¹³CD₂]-formaldehyde permit quantification of both the mass-labeled exogenous and endogenous DNA-formaldehyde reaction products. An existing pharmacokinetic model developed initially to describe ¹⁴C-DNA-protein crosslinks (DPX) provided a template for describing the time course of mass-labeled adducts. Published datasets included both DPX and N₂-HO¹³CD₂-dG adducts measured after a single 6-h exposure to 0.7, 2, 6, 9, 10, or 15 ppm formaldehyde, after multi-day exposures to 2 ppm for 6 h/day, 7 days/week with interim sacrifices up to 28 days, and after 28-day exposures for 6 h/day, 7 days/week to 0.3, 0.03, or 0.001 ppm. The existing kinetic model overpredicted endogenous adducts in the nasal epithelium after 1-day [¹³CD₂]-formaldehyde exposure, requiring adjustment of parameters for rates of tissue metabolism and background formaldehyde. After refining tissue formaldehyde parameters, we fit the model to both forms of adducts by varying key parameters and optimizing against all 3 studies. Fitting to all these studies required 2 nonlinear pathways—one for high-exposure saturation of clearance in the nasal epithelial tissues and another for extracellular clearance that restricts uptake into the epithelial tissue for inhaled concentrations below 0.7 ppm. This refined pharmacokinetic model for endogenous and exogenous formaldehyde acetal adducts can assist in updating biologically based dose-response models for formaldehyde carcinogenicity.

Extended Analysis and Evidence Integration of Chloroprene as a Human Carcinogen

β-Chloroprene is used in the production of polychloroprene, a synthetic rubber. In 2010, Environmental Protection Agency (EPA) published the Integrated Risk Information System "Toxicological Review of Chloroprene," concluding that chloroprene was "likely to be carcinogenic to humans." This was based on findings from a 1998 National Toxicology Program (NTP) study showing multiple tumors within and across animal species; results from occupational epidemiological studies; a proposed mutagenic mode of action; and structural similarities with 1,3-butadiene and vinyl chloride. Using mouse data from the NTP study and assuming a mutagenic mode of action, EPA calculated an inhalation unit risk (IUR) for chloroprene of 5 × 10^-4 per µg/m³ . This is among the highest IURs for chemicals classified by IARC or EPA as known or probable human carcinogens and orders of magnitude higher than the IURs for carcinogens such as vinyl chloride, benzene, and 1,3-butadiene. Due to differences in pharmacokinetics, mice appear to be uniquely responsive to chloroprene exposure compared to other animals, including humans, which is consistent with the lack of evidence of carcinogenicity in robust occupational epidemiological studies. We evaluated and integrated all lines of evidence for chloroprene carcinogenicity to assess whether the 2010 EPA IUR could be scientifically substantiated. Due to clear interspecies differences in carcinogenic response to chloroprene, we applied a physiologically based pharmacokinetic model for chloroprene to calculate a species-specific internal dose (amount metabolized/gram of lung tissue) and derived an IUR that is over 100-fold lower than the 2010 EPA IUR. Therefore, we recommend that EPA's IUR be updated.

Incorporation of in vitro metabolism data and physiologically based pharmacokinetic modeling in a risk assessment for chloroprene

Objective: To develop a physiologically based pharmacokinetic (PBPK) model for chloroprene in the mouse, rat and human, relying only on in vitro data to estimate tissue metabolism rates and partitioning, and to apply the model to calculate an inhalation unit risk (IUR) for chloroprene.Materials and methods: Female B6C3F1 mice were the most sensitive species/gender for lung tumors in the 2-year bioassay conducted with chloroprene. The PBPK model included tissue metabolism rate constants for chloroprene estimated from results of in vitro gas uptake studies using liver and lung microsomes. To assess the validity of the PBPK model, a 6-hr, nose-only chloroprene inhalation study was conducted with female B6C3F1 mice in which both chloroprene blood concentrations and ventilation rates were measured. The PBPK model was then used to predict dose measures - amounts of chloroprene metabolized in lungs per unit time - in mice and humans.Results: The mouse PBPK model accurately predicted in vivo pharmacokinetic data from the 6-hr, nose-only chloroprene inhalation study. The PBPK model was used to conduct a cancer risk assessment based on metabolism of chloroprene to reactive epoxides in the lung, the target tissue in mice. The IUR was over100-fold lower than the IUR from the EPA Integrated Risk Information System (IRIS), which was based on inhaled chloroprene concentration. The different result from the PBPK model risk assessment arises from use of the more relevant tissue dose metric, amount metabolized, rather than inhaled concentrationDiscussion and conclusions: The revised chloroprene PBPK model is based on the best available science, including new test animal in vivo validation, updated literature review and a Markov-Chain Monte Carlo analysis to assess parameter uncertainty. Relying on both mouse and human metabolism data also provides an important advancement in the use of quantitative in vitro to in vivo extrapolation (QIVIVE). Inclusion of the best available science is especially important when deriving a toxicity value based on species extrapolation for the potential carcinogenicity of a reactive metabolite.

Assessing children's exposure to manganese in drinking water using a PBPK model

A previously published human PBPK model for manganese (Mn) in infants and children has been updated with Mn in drinking water as an additional exposure source. Built upon the ability to capture differences in Mn source-specific regulation of intestinal uptake in nursing infants who are breast-fed and formula-fed, the updated model now describes the bioavailability of Mn from drinking water in children of ages 0-18. The age-related features, including the recommended age-specific Mn dietary intake, age-specific water consumption rates, and age-specific homeostasis of Mn, are based on the available human data and knowledge of the biology of essential-metal homeostasis. Model simulations suggest that the impact of adding drinking-water exposure to daily Mn exposure via dietary intake and ambient air inhalation in children is not greater than the impacts in adults, even at a drinking-water concentration that is 2 times higher than the USEPA's lifetime health advisory value. This conclusion was also valid for formula-fed infants who are considered at the highest potential exposure to Mn from drinking water compared to all other age groups. Our multi-route, multi-source Mn PBPK model for infants and children provides insights about the potential for Mn-related health effects on growing children and will thereby improve the level of confidence in properly interpreting Mn exposure-health effects relationships in children in human epidemiological studies.

Peak Exposures in Epidemiologic Studies and Cancer Risks: Considerations for Regulatory Risk Assessment

We review approaches for characterizing "peak" exposures in epidemiologic studies and methods for incorporating peak exposure metrics in dose-response assessments that contribute to risk assessment. The focus was on potential etiologic relations between environmental chemical exposures and cancer risks. We searched the epidemiologic literature on environmental chemicals classified as carcinogens in which cancer risks were described in relation to "peak" exposures. These articles were evaluated to identify some of the challenges associated with defining and describing cancer risks in relation to peak exposures. We found that definitions of peak exposure varied considerably across studies. Of nine chemical agents included in our review of peak exposure, six had epidemiologic data used by the U.S. Environmental Protection Agency (US EPA) in dose-response assessments to derive inhalation unit risk values. These were benzene, formaldehyde, styrene, trichloroethylene, acrylonitrile, and ethylene oxide. All derived unit risks relied on cumulative exposure for dose-response estimation and none, to our knowledge, considered peak exposure metrics. This is not surprising, given the historical linear no-threshold default model (generally based on cumulative exposure) used in regulatory risk assessments. With newly proposed US EPA rule language, fuller consideration of alternative exposure and dose-response metrics will be supported. "Peak" exposure has not been consistently defined and rarely has been evaluated in epidemiologic studies of cancer risks. We recommend developing uniform definitions of "peak" exposure to facilitate fuller evaluation of dose response for environmental chemicals and cancer risks, especially where mechanistic understanding indicates that the dose response is unlikely linear and that short-term high-intensity exposures increase risk.

Updating physiologically based pharmacokinetic models for manganese by incorporating rapid association/dissociation processes in tissues

Previously, we developed a series of physiologically based pharmacokinetic (PBPK) models for manganese (Mn) in which saturable tissue binding and dose-dependent increases in biliary excretion captured key aspects of Mn homeostasis biology. These models reproduced the non-linear behavior of Mn kinetics in different tissues, accounting for dose-dependent changes in Mn kinetics. The original model construct had relatively slow association and dissociation rate constants for Mn binding in tissues. In this updated model, both rates of entry into tissue and the interaction of Mn with binding sites are rapid, and the step limiting Mn accumulation is the saturation of tissue binding sites. This binding reflects general cellular requirements for Mn with high affinity but rapid exchange between bound and free forms, which we captured using a dissociation constant (KD) of ~ 0.5 μM across tissues while maintaining different maximum binding capacities in each tissue. Variability in the binding capacities accounted for different background levels of Mn in particular tissues. This alternative structure successfully described Mn kinetics in tissues in adult rats exposed to Mn either in their diet or by inhalation, indicating that both the original and the present models capture the dose-dependent and tissue-specific kinetic behavior of Mn in adult rats. Although the published models that emphasize the role of smaller tissue binding rate constants in non-linear behaviors capture all relevant dose-dependent kinetic behaviors of this metal, increasing biological relevance of the model structure and parameters should provide greater confidence in applying the Mn PBPK models to risk assessment.

Considerations for refining the risk assessment process for formaldehyde: Results from an interdisciplinary workshop

Anticipating the need to evaluate and integrate scientific evidence to inform new risk assessments or to update existing risk assessments, the Formaldehyde Panel of the American Chemistry Council (ACC), in collaboration with the University of North Carolina, convened a workshop: "Understanding Potential Human Health Cancer Risk - From Data Integration to Risk Evaluation" in October 2017. Twenty-four (24) invited-experts participated with expertise in epidemiology, toxicology, science integration and risk evaluation. Including members of the organizing committee, there were 29 participants. The meeting included eleven presentations encompassing an introduction and three sessions: (1) "integrating the formaldehyde science on nasal/nasopharyngeal carcinogenicity and potential for causality"; (2) "integrating the formaldehyde science on lymphohematopoietic cancer and potential for causality; and, (3) "formaldehyde research-data suitable for risk assessment". Here we describe key points from the presentations on epidemiology, toxicology and mechanistic studies that should inform decisions about the potential carcinogenicity of formaldehyde in humans and the discussions about approaches for structuring an integrated, comprehensive risk assessment for formaldehyde. We also note challenges expected when attempting to reconcile divergent results observed from research conducted within and across different scientific disciplines - especially toxicology and epidemiology - and in integrating diverse, multi-disciplinary mechanistic evidence.

Application of the adverse outcome pathway (AOP) approach to inform mode of action (MOA): A case study with inorganic arsenic

The aim of this study was to establish a process for deriving a chemical-specific mode of action (MOA) from chemical-agnostic adverse outcome pathway (AOPs), using inorganic arsenic (iAs) as a case study. The AOP developed for this case study are related to disruption of cellular signaling by chemicals that strongly bind to vicinal dithiols in cellular proteins, leading to disruption of inflammatory and oxidative stress signaling along with inhibition of the DNA damage responses. The proposed MOA for iAs incorporates this AOP, overlaid on a background of increasing oxidative stress and/or co-exposure to mutagenic chemicals or radiation. The most challenging aspect of developing a MOA from AOP is the incorporation of metabolism and dose-response, neither of which may be considered in the development of an AOP. The cellular responses to relatively low concentrations (below 100 parts per billion) of iAs in drinking water appear to be secondary to binding of trivalent arsenite and its trivalent metabolite, monomethyl arsenous acid to key cellular vicinal dithiols in target tissues, resulting in a co-carcinogenic MOA. The proposed AOP may also be applied to non-cancer endpoints, enabling an integrated approach to conducting a risk assessment for iAs.

Quantitative risk assessment of tobacco products: A potentially useful component of substantial equivalence evaluations

Quantitative risk assessment (QRA), a scientific, evidence-based analytical process that combines chemical and biological data to quantify the probability and potential impact of some defined risk, is used by regulatory agencies for decision-making. Thus, in tobacco product regulation, specifically in substantial equivalence (SE) evaluations, QRA can provide a useful, practical, and efficient approach to address questions that might arise regarding human health risk and potential influence on public health. In SE reporting, when differences in product characteristics may necessitate the determination of whether a new product raises different questions of public health, the results from QRA are a valuable metric. An approach for QRA in this context is discussed, which is modeled after the methodology for assessment of constituent mixtures by the US Environmental Protection Agency for environmental Superfund site assessment. Given the intent in both cases is an assessment of the public health impact resulting from the totality of exposure to a mixture of constituents, the application is appropriate. Although some uncertainties in the information incorporated may exist, relying on the most appropriate of the available data increases the confidence and decreases the uncertainty in the risk characterization using this data-driven methodology.

Six years after the NRC review of EPA's Draft IRIS Toxicological Review of Formaldehyde: Regulatory implications of new science in evaluating formaldehyde leukemogenicity

Shortly after the International Agency for Research on Cancer (IARC) determined that formaldehyde causes leukemia, the United States Environmental Protection Agency (EPA) released its Draft IRIS Toxicological Review of Formaldehyde ("Draft IRIS Assessment"), also concluding that formaldehyde causes leukemia. Peer review of the Draft IRIS Assessment by a National Academy of Science committee noted that "causal determinations are not supported by the narrative provided in the draft" (NRC 2011). They offered recommendations for improving the Draft IRIS assessment and identified several important research gaps. Over the six years since the NRC peer review, significant new science has been published. We identify and summarize key recommendations made by NRC and map them to this new science, including extended analysis of epidemiological studies, updates of earlier occupational cohort studies, toxicological experiments using a sensitive mouse strain, mechanistic studies examining the role of exogenous versus endogenous formaldehyde in bone marrow, and several critical reviews. With few exceptions, new findings are consistently negative, and integration of all available evidence challenges the earlier conclusions that formaldehyde causes leukemia. Given formaldehyde's commercial importance, environmental ubiquity and endogenous production, accurate hazard classification and risk evaluation of whether exposure to formaldehyde from occupational, residential and consumer products causes leukemia are critical.

Evaluation of triclosan in Minnesota lakes and rivers: Part II - human health risk assessment

Triclosan, an antimicrobial compound found in consumer products, has been detected in low concentrations in Minnesota municipal wastewater treatment plant (WWTP) effluent. This assessment evaluates potential health risks for exposure of adults and children to triclosan in Minnesota surface water, sediments, and fish. Potential exposures via fish consumption are considered for recreational or subsistence-level consumers. This assessment uses two chronic oral toxicity benchmarks, which bracket other available toxicity values. The first benchmark is a lower bound on a benchmark dose associated with a 10% risk (BMDL₁₀) of 47mg per kilogram per day (mg/kg-day) for kidney effects in hamsters. This value was identified as the most sensitive endpoint and species in a review by Rodricks et al. (2010) and is used herein to derive an estimated reference dose (RfD_(Rodricks)) of 0.47mg/kg-day. The second benchmark is a reference dose (RfD) of 0.047mg/kg-day derived from a no observed adverse effect level (NOAEL) of 10mg/kg-day for hepatic and hematopoietic effects in mice (Minnesota Department of Health [MDH] 2014). Based on conservative assumptions regarding human exposures to triclosan, calculated risk estimates are far below levels of concern. These estimates are likely to overestimate risks for potential receptors, particularly because sample locations were generally biased towards known discharges (i.e., WWTP effluent).

Does occupational exposure to formaldehyde cause hematotoxicity and leukemia-specific chromosome changes in cultured myeloid progenitor cells?

Several cross-sectional studies of a single population of workers exposed to formaldehyde at one of two factories using or producing formaldehyde-melamine resins in China have concluded that formaldehyde exposure induces damage to hematopoietic cells that originate in the bone marrow. Moreover, the investigators interpret observed differences between groups as evidence that formaldehyde induces myeloid leukemias, although the mechanisms for inducing these diseases are not obvious and recently published scientific findings do not support causation. Our objective was to evaluate hematological parameters and aneuploidy in relation to quantitative exposure measures of formaldehyde. We obtained the study data for the original study (Zhang et al. 2010 ) and performed linear regression analyses. Results showed that differences in white blood cell, granulocyte, platelet, and red blood cell counts are not exposure dependent. Among formaldehyde-exposed workers, no association was observed between individual average formaldehyde exposure estimates and frequency of aneuploidy, suggested by the original study authors to be indicators of myeloid leukemia risk.

A tissue dose-based comparative exposure assessment of manganese using physiologically based pharmacokinetic modeling-The importance of homeostatic control for an essential metal

A physiologically-based pharmacokinetic (PBPK) model (Schroeter et al., 2011) was applied to simulate target tissue manganese (Mn) concentrations following occupational and environmental exposures. These estimates of target tissue Mn concentrations were compared to determine margins of safety (MOS) and to evaluate the biological relevance of applying safety factors to derive acceptable Mn air concentrations. Mn blood concentrations measured in occupational studies permitted verification of the human PBPK models, increasing confidence in the resulting estimates. Mn exposure was determined based on measured ambient air Mn concentrations and dietary data in Canada and the United States (US). Incorporating dietary and inhalation exposures into the models indicated that increases in target tissue concentrations above endogenous levels only begin to occur when humans are exposed to levels of Mn in ambient air (i.e. >10μg/m³) that are far higher than those currently measured in Canada or the US. A MOS greater than three orders of magnitude was observed, indicating that current Mn air concentrations are far below concentrations that would be required to produce the target tissue Mn concentrations associated with subclinical neurological effects. This application of PBPK modeling for an essential element clearly demonstrates that the conventional application of default factors to "convert" an occupational exposure to an equivalent continuous environmental exposure, followed by the application of safety factors, is not appropriate in the case of Mn. PBPK modeling demonstrates that the relationship between ambient Mn exposures and dose-to-target tissue is not linear due to normal tissue background levels and homeostatic controls.

Requirements for a Biologically Realistic Cancer Risk Assessment for Inorganic Arsenic

A remarkable feature of the carcinogenicity of inorganic arsenic (As,) is the observation that human exposures to Asi have been strongly associated with increases in skin, lung, and internal cancers, but As, does not typically cause tumors in standard laboratory animal test protocols. Considerable controversy has centered on whether there is epidemiological evidence of a “threshold” for the carcinogenic effects of Asi, or at least of a highly nonlinear dose–response. Saturation of metabolism in the dose-range associated with tumors does not appear to be adequate to produce a major impact on the dose-response for carcinogenicity. If there is a strong nonlinearity, it results from the nature of the carcinogenic mechanism(s) of Asi. However, no single hypothesis for the mechanism of Asi carcinogenicity has widespread support. A biologically realistic cancer risk assessment for Asi would requirea quantitative description of the dose of active arsenic species in target tissues, the interactions between active arsenic and tissue constituents, and the manner in which these interactions result in tumor formation in multiple organs in humans, but not in experimental animals. Although Asi has only infrequently been associated with tumors in animal studies, it has repeatedly been shown to act as a comutagen in vitro and as a cocarcinogen in vivo. Asi is clastogenic, producing chromatid aberrations, but does not produce point mutations at single gene loci. Of particular interest, Asi has been shown to inhibit repair of DNA single-strand breaks, a possible mechanism for its observed comutagenicity and cocarcinogenicity. We propose a cocarcinogenic mode of action in which Asi acts primarily on intermediate cells deficient in cell cycle control at a late stage in a preexisting carcinogenic process. This interaction enhances ge-nomic fragility and accelerates conversion of premalignant lesions to more aggressive, clinically observable tumors. An indirect effect of As, on DNA repair is consistent with the expectation of a nonlinear dose-response rather than the linear dose-response traditionally assumed for mutagenic carcinogens. However, defining the exact nature of this tumor dose-response will require further experimental data on the dose-response for the cellular effects of Asi. Because Asi carcinogenicity is unlikely to be observed in normal experimental animals not exposed to other carcinogens, studies in animals and cell lines deficient in cell cycle control should also be considered. Experimental studies specifically designed to address the key mechanistic and dose-response issues for Asi carcinogenicity are critically needed to support public health policy decisions regarding current environmental exposures to Asi.

Evaluation of gene expression changes in human primary lung epithelial cells following 24-hr exposures to inorganic arsenic and its methylated metabolites and to arsenic trioxide

The concentration response for altered gene expression in primary lung epithelial cells was determined following two treatments with arsenicals: (1) a mixture of trivalent arsenic compounds representative of urinary arsenic concentrations in exposed human populations, and (2) arsenite (As2 O3 ) a common form of inhaled arsenic dust that is frequently used in both in vivo and in vitro experimental exposures. Biochemical assays did not detect any evidence of cytotoxicity at the concentrations used, apart from a concentration-related increase in cellular heme oxygenase that was also indicated by the genomic analysis. Cell signal pathway enrichment analysis indicated similar responses to both treatments, with concentration-related responses in pathways related to cell adhesion, cytoskeleton remodeling, development (morphogenesis), cell cycle control, and to a lesser extent inflammatory responses. These cellular responses to arsenic were consistent with those observed in a previous study with primary uroepithelial cells. Benchmark dose analysis also demonstrated similar potency of the two treatments as well as comparable sensitivity of the two cell types. A number of genes showing similar concentration-dependent expression across individuals in both bladder and lung cells were identified, including heme oxygenase 1, thioredoxin reductase, DNA damage binding protein 2, and thrombomodulin. The data on human primary lung cells from this study, together with the data from human primary uroepithelial cells, support a conclusion that biological responses to arsenic by human cells under study conditions are unlikely to occur at concentrations below 0.1 µM. Environ. Mol. Mutagen. 56:477-490, 2015. © 2015 Wiley Periodicals, Inc.

Risk assessments for chronic exposure of children and prospective parents to ethylbenzene (CAS No. 100-41-4)

Potential chronic health risks for children and prospective parents exposed to ethylbenzene were evaluated in response to the Voluntary Children's Chemical Evaluation Program. Ethylbenzene exposure was found to be predominately via inhalation with recent data demonstrating continuing decreases in releases and both outdoor and indoor concentrations over the past several decades. The proportion of ethylbenzene in ambient air that is attributable to the ethylbenzene/styrene chain of commerce appears to be relatively very small, less than 0.1% based on recent relative emission estimates. Toxicity reference values were derived from the available data, with physiologically based pharmacokinetic models and benchmark dose methods used to assess dose-response relationships. An inhalation non-cancer reference concentration or RfC of 0.3 parts per million (ppm) was derived based on ototoxicity. Similarly, an oral non-cancer reference dose or RfD of 0.5 mg/kg body weight/day was derived based on liver effects. For the cancer assessment, emphasis was placed upon mode of action information. Three of four rodent tumor types were determined not to be relevant to human health. A cancer reference value of 0.48 ppm was derived based on mouse lung tumors. The risk characterization for ethylbenzene indicated that even the most highly exposed children and prospective parents are not at risk for non-cancer or cancer effects of ethylbenzene.

Potential occupational risk of amines in carbon capture for power generation

PURPOSE

While CO2 capture and storage (CCS) technology has been well studied in terms of its efficacy and cost of implementation, there is limited available data concerning the potential for occupational exposure to amines, mixtures of amines, or degradation of by-products from the CCS process. This paper is a critical review of the available data concerning the potential effects of amines and CCS-degradation by-products.

METHODS

A comprehensive review of the occupational health and safety issues associated with exposure to amines and amine by-products at CCS facilities was performed, along with a review of the regulatory status and guidelines of amines, by-products, and CCS process vapor mixtures.

RESULTS

There are no specific guidelines or regulations regarding permissible levels of exposure via air for amines and degradation products that could form atmospheric oxidation of amines released from post-combustion CO2 capture plants. While there has been a worldwide effort to develop legal and regulatory frameworks for CCS, none are directly related to occupational exposures.

CONCLUSIONS

By-products of alkanolamine degradation may pose the most significant health hazard to workers in CCS facilities, with several aldehydes, amides, nitramines, and nitrosamines classified as either known or potential/possible human carcinogens. The absence of large-scale CCS facilities; absence and unreliability of reported data in the literature from pilot facilities; and proprietary amine blends make it difficult to estimate potential amine exposures and predict formation and exposure to degradation products.

Use of mode of action data to inform a dose-response assessment for bladder cancer following exposure to inorganic arsenic

In the recent National Research Council report on conducting a dose-response assessment for inorganic arsenic, the committee remarked that mode of action data should be used, to the extent possible, to extrapolate below the observed range for epidemiological studies to inform the shape of the dose-response curve. Recent in vitro mode of action studies focused on understanding the development of bladder cancer following exposure to inorganic arsenic provide data to inform the dose-response curve. These in vitro data, combined with results of bladder cancer epidemiology studies, inform the dose-response curve in the low-dose region, and include values for both pharmacokinetic and pharmacodynamic variability. Integration of these data provides evidence of a range of concentrations of arsenic for which no effect on the bladder would be expected. Specifically, integration of these results suggest that arsenic exposures in the range of 7-43 ppb in drinking water are exceedingly unlikely to elicit changes leading to key events in the development of cancer or noncancer effects in bladder tissue. These findings are consistent with the lack of evidence for bladder cancer following chronic ingestion of arsenic water concentrations <100 ppb in epidemiological studies.

Formaldehyde exposure and leukemia: critical review and reevaluation of the results from a study that is the focus for evidence of biological plausibility

A recent study (Zhang et al., 2010) has provided results attributed to aneuploidy in circulating stem cells that has been characterized as providing potential support for proposed mechanisms for formaldehyde to impact bone marrow. A critical review of the study, as well as a reanalysis of the underlying data, was performed and the results of this reanalysis suggested factors other than formaldehyde exposure may have contributed to the effects reported. In addition, although the authors stated in their paper that "all scorable metaphase spreads on each slide were analyzed, and a minimum of 150 cells per subject was scored," this protocol was not followed. In fact, the protocol to evaluate the presence of monosomy 7 or trisomy 8 was followed for three or less samples in exposed workers and six or less samples in non-exposed workers. In addition, the assays used (CFU-GM) do not actually measure the proposed events in primitive cells involved in the development of acute myeloid leukemia. Evaluation of these data indicates that the aneuploidy measured could not have arisen in vivo, but rather arose during in vitro culture. The results of our critical review and reanalysis of the data, in combination with recent toxicological and mechanistic studies, do not support a mechanism for a causal association between formaldehyde exposure and myeloid or lymphoid malignancies.

Physiologically based pharmacokinetic/toxicokinetic modeling

Physiologically based pharmacokinetic (PBPK) models differ from conventional compartmental pharmacokinetic models in that they are based to a large extent on the actual physiology of the organism. The application of pharmacokinetics to toxicology or risk assessment requires that the toxic effects in a particular tissue are related in some way to the concentration time course of an active form of the substance in that tissue. The motivation for applying pharmacokinetics is the expectation that the observed effects of a chemical will be more simply and directly related to a measure of target tissue exposure than to a measure of administered dose. The goal of this work is to provide the reader with an understanding of PBPK modeling and its utility as well as the procedures used in the development and implementation of a model to chemical safety assessment using the styrene PBPK model as an example.

Analysis of genomic dose-response information on arsenic to inform key events in a mode of action for carcinogenicity

A comprehensive literature search was conducted to identify information on gene expression changes following exposures to inorganic arsenic compounds. This information was organized by compound, exposure, dose/concentration, species, tissue, and cell type. A concentration-related hierarchy of responses was observed, beginning with changes in gene/protein expression associated with adaptive responses (e.g., preinflammatory responses, delay of apoptosis). Between 0.1 and 10 microM, additional gene/protein expression changes related to oxidative stress, proteotoxicity, inflammation, and proliferative signaling occur along with those related to DNA repair, cell cycle G2/M checkpoint control, and induction of apoptosis. At higher concentrations (10-100 microM), changes in apoptotic genes dominate. Comparisons of primary cell results with those obtained from immortalized or tumor-derived cell lines were also evaluated to determine the extent to which similar responses are observed across cell lines. Although immortalized cells appear to respond similarly to primary cells, caution must be exercised in using gene expression data from tumor-derived cell lines, where inactivation or overexpression of key genes (e.g., p53, Bcl-2) may lead to altered genomic responses. Data from acute in vivo exposures are of limited value for evaluating the dose-response for gene expression, because of the transient, variable, and uncertain nature of tissue exposure in these studies. The available in vitro gene expression data, together with information on the metabolism and protein binding of arsenic compounds, provide evidence of a mode of action for inorganic arsenic carcinogenicity involving interactions with critical proteins, such as those involved in DNA repair, overlaid against a background of chemical stress, including proteotoxicity and depletion of nonprotein sulfhydryls. The inhibition of DNA repair under conditions of toxicity and proliferative pressure may compromise the ability of cells to maintain the integrity of their DNA.

Evaluation of Physiologically Based Pharmacokinetic Models in Risk Assessment: An Example with Perchloroethylene

One of the more problematic aspects of the application of physiologically based pharmacokinetic (PBPK) models in risk assessment is the question of whether the model has been adequately validated to provide confidence in the dose metrics calculated with it. A number of PBPK models have been developed for perchloroethylene (PCE), differing primarily in the parameters estimated for metabolism. All of the models provide reasonably accurate simulations of selected kinetic data for PCE in mice and humans and could thus be considered to be "validated" to some extent. However, quantitative estimates of PCE cancer risk are critically dependent on the prediction of the rate of metabolism at low environmental exposures. Recent data on the urinary excretion of trichloroacetic acid (TCA), the major metabolite of PCE, for human subjects exposed to lower concentrations than those used in previous studies, make it possible to compare the high- to low-dose extrapolation capability of the various published human models. The model of Gearhart et al., which is the only model to include a description of TCA kinetics, provided the closest predictions of the urinary excretion observed in these low-concentration exposures. Other models overestimated metabolite excretion in this study by 5- to 15-fold. A systematic discrepancy between model predictions and experimental data for the time course of the urinary excretion of TCA suggested a contribution from TCA formed by metabolism of PCE in the kidney and excreted directly into the urine. A modification of the model of Gearhart et al. to include metabolism of PCE to TCA in the kidney at 10% of the capacity of the liver, with direct excretion of the TCA formed in the kidney into the urine, markedly improved agreement with the experimental time-course data, without altering predictions of liver metabolism. This case study with PCE demonstrates the danger of relying on parent chemical kinetic data to validate a model that will be used for the prediction of metabolism.

Research toward the development of a biologically based dose response assessment for inorganic arsenic carcinogenicity: a progress report

Cancer risk assessments for inorganic arsenic have been based on human epidemiological data, assuming a linear dose response below the range of observation of tumors. Part of the reason for the continued use of the linear approach in arsenic risk assessments is the lack of an adequate biologically based dose response (BBDR) model that could provide a quantitative basis for an alternative nonlinear approach. This paper describes elements of an ongoing collaborative research effort between the CIIT Centers for Health Research, the U.S. Environmental Protection Agency, ENVIRON International, and EPRI to develop BBDR modeling approaches that could be used to inform a nonlinear cancer dose response assessment for inorganic arsenic. These efforts are focused on: (1) the refinement of physiologically based pharmacokinetic (PBPK) models of the kinetics of inorganic arsenic and its metabolites in the mouse and human; (2) the investigation of mathematical solutions for multi-stage cancer models involving multiple pathways of cell transformation; (3) the review and evaluation of the literature on the dose response for the genomic effects of arsenic; and (4) the collection of data on the dose response for genomic changes in the urinary bladder (a human target tissue for arsenic carcinogenesis) associated with in vivo drinking water exposures in the mouse as well as in vitro exposures of both mouse and human cells. An approach is proposed for conducting a biologically based margin of exposure risk assessment for inorganic arsenic using the in vitro dose response for the expression of genes associated with the obligatory precursor events for arsenic tumorigenesis.

The use of Markov chain Monte Carlo uncertainty analysis to support a Public Health Goal for perchloroethylene

The current Public Health Goal (PHG) for perchloroethylene (PCE) was derived using upper-bound estimates of fractional PCE metabolism in humans. These estimates were in part obtained from a published evaluation of the uncertainty and variability in human PCE metabolism conducted using a physiologically-based pharmacokinetic (PBPK) model in a Markov chain Monte Carlo (MCMC) analysis; however, the data used in that analysis were limited to post-exposure PCE blood and exhaled air concentrations from a single study. A more recent study [Volkel, W., Friedewald, M., Lederer, E., Pahler, A., Parker, J., Dekant, W., 1998. Biotransformation of perchloroethene: dose-dependent excretion of trichloroacetic acid, dichloroacetic acid, and N-acetyl-S-(trichlorovinyl)-l-cysteine in rats and humans after inhalation. Toxicol. Appl. Pharmacol. 153(1), 20-27.] provides data on blood concentrations of PCE and its major metabolite, trichloroacetic acid (TCA), and urinary excretion of TCA following exposure of human subjects to lower concentrations of PCE (10-40ppm) than in previous studies. In the present effort, a new MCMC analysis was performed that focused on data from this study along with two others [Fernandez, J., Guberan, E., Caperos, J., 1976. Experimental human exposures to tetrachloroethylene vapor and elimination in breath after inhalation. Am. Ind. Hyg. Assoc. J. 37, 143-150; Monster, A., Boersma, G., Steenweg, H., 1979. Kinetics of tetrachloroethylene in volunteers; influence of exposure concentration and work load. Int. Arch. Occup. Environ. Health 42, 303-309.] providing data on PCE blood concentrations and urinary excretion of TCA. To provide an accurate prediction of TCA kinetics, the PBPK model used here includes a description of the metabolism of PCE to TCA in both the liver and kidney. The resulting upper 95th percentile estimates of fraction of PCE metabolized by inhalation and oral routes were 2.1 and 5.2%, respectively, compared to 58 and 79% used in the derivation of the PHG.

Revised assessment of cancer risk to dichloromethane II. Application of probabilistic methods to cancer risk determinations

An updated PBPK model of methylene chloride (DCM, dichloromethane) carcinogenicity in mice was recently published using Bayesian statistical methods (Marino et al., 2006). In this work, this model was applied to humans, as recommended by Sweeney et al.(2004). Physiological parameters for input into the MCMC analysis were selected from multiple sources reflecting, in each case, the source that was considered to represent the most current scientific evidence for each parameter. Metabolic data for individual subjects from five human studies were combined into a single data set and population values derived using MCSim. These population values were used for calibration of the human model. The PBPK model using the calibrated metabolic parameters was used to perform a cancer risk assessment for DCM, using the same tumor incidence and exposure concentration data relied upon in the current IRIS entry. Unit risks, i.e., the risk of cancer from exposure to 1 microg/m3 over a lifetime, for DCM were estimated using the calibrated human model. The results indicate skewed distributions for liver and lung tumor risks, alone or in combination, with a mean unit risk (per microg/m3) of 1.05 x 10(-9), considering both liver and lung tumors. Adding the distribution of genetic polymorphisms for metabolism to the ultimate carcinogen, the unit risks range from 0 (which is expected given that approximately 20% of the US population is estimated to be nonconjugators) up to a unit risk of 2.70 x 10(-9) at the 95th percentile. The median, or 50th percentile, is 9.33 x 10(-10), which is approximately a factor of 500 lower than the current EPA unit risk of 4.7 x 10(-7) using a previous PBPK model. These values represent the best estimates to date for DCM cancer risk because all available human data sets were used, and a probabilistic methodology was followed.

Revised assessment of cancer risk to dichloromethane: Part I Bayesian PBPK and dose–response modeling in mice

The current USEPA cancer risk assessment for dichloromethane (DCM) is based on deterministic physiologically based pharmacokinetic (PBPK) modeling involving comparative metabolism of DCM by the GST pathway in the lung and liver of humans and mice. Recent advances in PBPK modeling include probabilistic methods and, in particular, Bayesian inference to quantitatively address variability and uncertainty separately. Although Bayesian analysis of human PBPK models has been published, no such efforts have been reported specifically addressing the mouse, apart from results included in the OSHA final rule on DCM. Certain aspects of the OSHA model, however, are not consistent with current approaches or with the USEPA's current DCM cancer risk assessment. Therefore, Bayesian analysis of the mouse PBPK model and dose-response modeling was undertaken to support development of an improved cancer risk assessment for DCM. A hierarchical population model was developed and prior parameter distributions were selected to reflect parameter values that were considered the most appropriate and best available. Bayesian modeling was conducted using MCSim, a publicly available software program for Markov Chain Monte Carlo analysis. Mean posterior values from the calibrated model were used to develop internal dose metrics, i.e., mg DCM metabolized by the GST pathway/L tissue/day in the lung and liver using exposure concentrations and results from the NTP mouse bioassay, consistent with the approach used by the USEPA for its current DCM cancer risk assessment. Internal dose metrics were 3- to 4-fold higher than those that support the current USEPA IRIS assessment. A decrease of similar magnitude was also noted in dose-response modeling results. These results show that the Bayesian PBPK model in the mouse provides an improved basis for a cancer risk assessment of DCM.

Interspecies dose extrapolation for inhaled dimethyl sulfate: a PBPK model-based analysis using nasal cavity N7-methylguanine adducts

Dimethyl sulfate (DMS) is a volatile sulfuric acid ester used principally as a methylating agent in a wide variety of industrial applications. DMS reacts with organic macromolecules by a SN2 mechanism. The weight of experimental evidence suggests that DMS possesses genotoxic and carcinogenic potential. Inhalation studies have shown that repeated exposure to DMS leads to tumors in the nasal cavity and lower respiratory tract in both rats and mice. Here we present a quantitative assessment for cross-species dose extrapolation for inhaled DMS using a physiologically based pharmacokinetic (PBPK) model. The model is designed to simulate N7-methylguanine (N7 mG) DNA adduct levels in the nasal mucosa following DMS exposure in rats and humans. This model was parameterized and predictions were tested by comparison against experimentally measured N7 mG DNA adduct levels in rat nasal mucosa following inhalation exposure to DMS. The model-based interspecies dose comparison, using N7 mG adduct levels in the nasal respiratory tissue as the appropriate dose metrics, predicts a dose rate seven times higher in rats compared to humans.

Comparison of tissue dosimetry in the mouse following chronic exposure to arsenic compounds

Several chronic bioassays have been conducted in multiple strains of mice in which various concentrations of arsenate or arsenite were administered in the drinking water without a tumorigenic effect. However, one study (Ng et al., 1999) reported a significant increase in tumor incidence in C57Bl/6J mice exposed to arsenic in their drinking water throughout their lifetime, with no tumors reported in controls. A physiologically based pharmacokinetic model for arsenic in the mouse has previously been developed (Gentry et al., 2004) to investigate potential differences in tissue dosimetry of arsenic species across various strains of mice. Initial results indicated no significant differences in blood, liver, or urine dosimetry in B6C3F1 and C57Bl/6 mice for acute or subchronic exposure. The current work was conducted to compare model-predicted estimates of tissue dosimetry to additional kinetic information from the (C57Bl/6 xCBA)F1 and TgAc mouse. The results from the current modeling indicate that the pharmacokinetic parameters derived based on information in the B6C3F1 mouse adequately describe the measured concentrations in the blood/plasma, liver, and urine of both the (C57Bl/6 x CBA)F1 and TgAc mouse, providing further support that the differences in response observed in the chronic bioassays are not related to strain-specific differences in pharmacokinetics. One significant finding was that no increases in skin or lung concentrations of arsenic species in the (C57Bl/6 x CBA)F1 strain were observed following administration of low concentrations (0.2 or 2 mg/U of arsenate in the drinking water, even though differences in response in the skin were reported. These data suggest that pharmacodynamic changes may be observed following exposure to arsenic compounds without an observable change in tissue dosimetry. These results provided further indirect support for the existence of inducible arsenic efflux in these tissues.

Evaluation of the Potential Impact of Age- and Gender-Specific Pharmacokinetic Differences on Tissue Dosimetry 2Current address: Novartis Pharmaceuticals, East Hanover, NJ 07936

The physiological and biochemical processes that determine the tissue concentration time courses (pharmacokinetics) of xenobiotics vary, in some cases significantly, with age and gender. While it is known that age- and gender-specific differences have the potential to affect tissue concentrations and, hence, individual risk, the relative importance of the contributing processes and the quantitative impact of these differences for various life stages are not well characterized. The objective of this study was to identify age- and gender-specific differences in physiological and biochemical processes that affect tissue dosimetry and integrate them into a predictive physiologically based pharmacokinetic (PBPK) life-stage model. The life-stage model was exercised for several environmental chemicals with a variety of physicochemical, biochemical, and mode-of-action properties. In general, predictions of average pharmacokinetic dose metrics for a chemical across life stages were within a factor of two, although larger transient variations were predicted, particularly during the neonatal period. The most important age-dependent pharmacokinetic factor appears to be the potential for decreased clearance of a toxic chemical in the perinatal period due to the immaturity of many metabolic enzyme systems, although this same factor may also reduce the production of a reactive metabolite. Given the potential for age-dependent pharmacodynamic factors during early life, there may be chemicals and health outcomes for which decreased clearance over a relatively brief period could have a substantial impact on risk.

Physiologically based pharmacokinetic modeling of arsenic in the mouse

A remarkable feature of the carcinogenicity of inorganic arsenic is that while human exposures to high concentrations of inorganic arsenic in drinking water are associated with increases in skin, lung, and bladder cancer, inorganic arsenic has not typically caused tumors in standard laboratory animal test protocols. Inorganic arsenic administered for periods of up to 2 yr to various strains of laboratory mice, including the Swiss CD-1, Swiss CR:NIH(S), C57Bl/6p53(+/-), and C57Bl/6p53(+/+), has not resulted in significant increases in tumor incidence. However, Ng et al. (1999) have reported a 40% tumor incidence in C57Bl/6J mice exposed to arsenic in their drinking water throughout their lifetime, with no tumors reported in controls. In order to investigate the potential role of tissue dosimetry in differential susceptibility to arsenic carcinogenicity, a physiologically based pharmacokinetic (PBPK) model for inorganic arsenic in the rat, hamster, monkey, and human (Mann et al., 1996a, 1996b) was extended to describe the kinetics in the mouse. The PBPK model was parameterized in the mouse using published data from acute exposures of B6C3F1 mice to arsenate, arsenite, monomethylarsonic acid (MMA), and dimethylarsinic acid (DMA) and validated using data from acute exposures of C57Black mice. Predictions of the acute model were then compared with data from chronic exposures. There was no evidence of changes in the apparent volume of distribution or in the tissue-plasma concentration ratios between acute and chronic exposure that might support the possibility of inducible arsenite efflux. The PBPK model was also used to project tissue dosimetry in the C57Bl/6J study, in comparison with tissue levels in studies having shorter duration but higher arsenic treatment concentrations. The model evaluation indicates that pharmacokinetic factors do not provide an explanation for the difference in outcomes across the various mouse bioassays. Other possible explanations may relate to strain-specific differences, or to the different durations of dosing in each of the mouse studies, given the evidence that inorganic arsenic is likely to be active in the later stages of the carcinogenic process.

Application of a physiologically based pharmacokinetic model for reference dose and reference concentration estimation for acetone

Recent health risk assessments to propose a Reference Dose (RfD) for acetone (Forsyth, 2001; U.S. EPA, 2001) have been based on the results of an oral subchronic study conducted in rats and mice (Dietz et al., 1991; NTP, 1991). These assessments have utilized the traditional concept of establishing the RfD by determining the lowest experimentally determined No-Observed-Adverse-Effect Level (NOAEL) and applying various Uncertainty Factors (UFs) (U.S. EPA, 1988). This article describes a risk assessment for acetone based on the systemic toxicity observed in subchronic and developmental toxicity studies to estimate an RfD and an inhalation reference concentration (RfC) for acetone. Specifically, this approach examined the subchronic study by Dietz et al. (1991), as well as an inhalation developmental toxicity study on acetone (Mast et al., 1988) and several toxicology studies of isopropanol (IPA). This was accomplished by applying a physiologically based pharmacokinetic (PBPK) model developed previously for IPA and its metabolite acetone (Clewell et al., 2001). The incorporation of the PBPK model into the derivation of an RfD and RfC for acetone allowed for a tissue-based approach rather than an external exposure-based approach, making it possible to derive an oral RfD from an inhalation study. In addition, the use of the PBPK model to analyze data from chronic and reproductive/developmental studies conducted with IPA enabled an assessment of the potential for acetone to produce any of the effects observed in the IPA studies. This analysis provided sufficient information to reduce the need for UFs in the adjustment of the NOAEL from the oral subchronic study for the determination of an RfD. Using the PBPK model in the acetone risk assessment supports a composite UF of 60 for the subchronic study, compared to composite factors of 300 to 3000 in the other recent risk assessments. This difference resulted in an RfD of 16 mg/kg/d, compared to the values of 0.3 to 3 that have previously been estimated (Forsyth, 2001; U.S. EPA, 2001). Considering the results from the inhalation developmental study (Mast et al., 1988) resulted in an RfD of 8.7 mg/kg/d. Using this study also fills a data gap for acetone that exists if only the oral database for acetone is considered for RfD derivation. An RfC of 29 ppm was also estimated for acetone using the Mast et al. (1988) study results in combination with the PBPK model. The potential impact of endogenous acetone on a risk assessment for acetone is also discussed.

Evaluation of the potential impact of age- and gender-specific lung morphology and ventilation rate on the dosimetry of vapors

In recent years, there have been growing concerns that due to differences, both pharmacokinetic and pharmacodynamic, between children and adults, children could be at greater risk of adverse effects following chemical exposure. The specific goal of this study was to demonstrate an approach for using physiologically based pharmacokinetic (PBPK) modeling to compare inhalation dosimetry in the adult and the child of both males and females. Three categories of gases were considered: rapidly and irreversibly reactive in the respiratory tract (ozone), relatively water-soluble and nonreactive (isopropanol), and relatively water-insoluble and nonreactive (styrene, vinyl chloride, and perchloroethylene). The nonreactive chemicals were also selected because they are metabolized in the respiratory tract. The age-related changes observed for the estimated dose metrics were a function of the physiochemical properties of the inhaled vapor and their interactions in the body. Blood concentrations estimated for all vapors, either poorly metabolized (e.g., PERC), moderately metabolized (e.g., ST), or highly metabolized vapors (e.g., IPA and VC), varied less than a factor of two between infants and adults. These changes, moreover, were confined to the first year after birth, a relatively short window compared to the total lifespan of the individual. In contrast, circulating metabolite concentrations estimated in the blood, as well as amounts metabolized in the liver and lung, appeared to be a strong function of age, due to their dependence on the maturity of the pertinent metabolic enzyme systems.

Evaluation of the potential impact of pharmacokinetic differences on tissue dosimetry in offspring during pregnancy and lactation

In recent years efforts have increased to develop a framework for assessing differences, both pharmacokinetic and pharmacodynamic, between children and adults for purposes of assessing risk of adverse effects following chemical exposure. The specific goal of this study was to demonstrate an approach for using PBPK modeling to compare maternal and fetal/neonatal blood and tissue dose metrics during pregnancy and lactation. Six chemical classes were targeted to provide a variety of physicochemical properties (volatility, lipophilicity, water solubility), and surrogate chemicals were selected to represent each class (isopropanol, vinyl chloride, methylene chloride, perchloroethylene, nicotine, and TCDD), based on the availability of pharmacokinetic information. These chemicals were also selected to provide different pharmacokinetic characteristics, including metabolic production of stable or reactive intermediates in the liver and competing pathways for metabolism. Changes in dosimetry during pregnancy predicted by the modeling were mainly attributable to the development of enzymatic pathways in the fetus or to changes in tissue composition in the mother and fetus during pregnancy. In general, blood concentrations were lower in the neonate during the lactation period than in the fetus during gestation. This postnatal decrease varied from only a slight change (for TCDD) to approximately four orders of magnitude (for vinyl chloride). As compared to maternal exposure, fetal/neonatal exposures ranged from approximately twice as great (for TCDD) to several orders of magnitude lower (for isopropanol). The results of this study are in general agreement with the analyses of data on pharmaceutical chemicals, which have suggested that the largest difference in pharmacokinetics observed between children and adults is for the perinatal period. The most important factor appears to be the potential for decreased clearance of toxic chemicals in the perinatal period due to immature metabolic enzyme systems, although this same factor can also reduce the risk from reactive metabolites during the same period.

An approach for the quantitative consideration of genetic polymorphism data in chemical risk assessment: examples with warfarin and parathion

In recent years, a great deal of research has been conducted to identify genetic polymorphisms. One focus has been to characterize variability in metabolic enzyme systems that could impact internal doses of pharmaceuticals or environmental pollutants. Methods are needed for using this metabolic information to estimate the resulting variability in tissue doses associated with chemical exposure. We demonstrate here the use of physiologically based pharmacokinetic (PBPK) modeling in combination with Monte Carlo analysis to incorporate information on polymorphisms into the analysis of toxicokinetic variability. Warfarin and parathion were used as case studies to demonstrate this approach. Our results suggest that polymorphisms in the PON1 gene, that give rise to allelic variants of paraoxonase, which is involved in the metabolism of paraoxon (a metabolite of parathion), make only a minor contribution to the overall variability in paraoxon tissue dose, while polymorphisms in the CYP2C9 gene, which gives rise to allelic variants of the major metabolic enzyme for warfarin, account for a significant portion of the overall variability in (S)-warfarin tissue dose. These analyses were used to estimate chemical-specific adjustment factors (CSAFs) for the human variability in toxicokinetics for both parathion and warfarin. Implications of alternatives in the calculation of CSAFs are explored. Key decision points for applying the PBPK-Monte Carlo approach to evaluate toxicokinetic variability for other chemicals are also discussed.

Genetic polymorphisms in assessing interindividual variability in delivered dose

Increasing sophistication in methods used to account for human variability in susceptibility to toxicants has been one of the success stories in the continuing evolution of risk assessment science. Genetic polymorphisms have been suggested as an important contributor to overall human variability. Recently, data on polymorphisms in metabolic enzymes have been integrated with physiologically based pharmacokinetic (PBPK) modeling as an approach to determining the resulting overall variability. We present an analysis of the potential contribution of polymorphisms in enzymes modulating the disposition of four diverse compounds: methylene chloride, warfarin, parathion, and dichloroacetic acid. Through these case studies, we identify key uncertainties likely to be encountered in the use of polymorphism data and highlight potential simplifying assumptions that might be required to test the hypothesis that genetic factors are a substantive source of human variability in susceptibility to environmental toxicants. These uncertainties include (1) the relative contribution of multiple enzyme systems, (2) the extent of induction/inhibition through coexposure, (3) allelic frequencies of major ethnic groups, (4) the absence of chemical-specific data on the kinetic parameters for the different allelic forms of key enzymes, (5) large numbers of low-frequency alleles, and (6) uncertainty regarding differences between in vitro and in vivo kinetic data. Our effort sets the stage for the acquisition of critical data and further integration of polymorphism data with PBPK modeling as a means to quantitate population variability.

Development of a physiologically based pharmacokinetic model of isopropanol and its metabolite acetone

A physiologically based pharmacokinetic (PBPK) model for isopropanol (IPA) and its major metabolite, acetone, is described. The structure of the parent chemical model, which can be used for either IPA or acetone by choosing the appropriate chemical-specific parameters, is similar to previously published models of volatile organic chemicals such as styrene. However, in order to properly simulate data on the exhalation of IPA and acetone during inhalation exposures, it was necessary to expand the description of the lung compartment to include a subcompartment for the upper respiratory tract mucus layer. This elaboration is consistent with published PBPK models of other water-soluble vapors in which the mucus layer serves to absorb the chemical during inhalation and then release it during exhalation. In the case of IPA exposure, a similar PBPK structure is used to describe the kinetics of the acetone produced from the metabolism of IPA. The resulting model is able to provide a coherent description of IPA and acetone kinetics in the rat and human for exposures to IPA by several routes: intravenous, intraperitoneal, oral, inhalation, and dermal. It is also able to consistently reproduce kinetic data for exposures of rats or humans to acetone. Thus, the model provides a validated framework for performing chemical-specific route-to-route extrapolation and cross-species dosimetry, which can be used in place of generic default calculations in support of risk assessments for IPA and acetone.

Comparison of cancer risk estimates for vinyl chloride using animal and human data with a PBPK model

Vinyl chloride (VC) is a trans-species carcinogen, producing tumors in a variety of tissues, from both inhalation and oral exposures, across a number of species. In particular, exposure to VC has been associated with a rare tumor, liver angiosarcoma, in a large number of studies in mice, rats, and humans. The mode of action for the carcinogenicity of VC appears to be a relatively straightforward example of DNA adduct formation by a reactive metabolite, leading to mutation, mistranscription, and neoplasia. The objective of the present analysis was to investigate the comparative potency of a classic genotoxic carcinogen across species, by performing a quantitative comparison of the carcinogenic potency of VC using data from inhalation and oral rodent bioassays as well as from human epidemiological studies. A physiologically-based pharmacokinetic (PBPK) model for VC was developed to support the target tissue dosimetry for the cancer risk assessment. Unlike previous models, the initial metabolism of VC was described as occurring via two saturable pathways, one representing low capacity-high affinity oxidation by CYP2E1 and the other (in the rodent) representing higher capacity-lower affinity oxidation by other isozymes of P450, producing in both cases chloroethylene oxide (CEO) and chloroacetaldehyde (CAA) as intermediate reactive products. Depletion of glutathione by reaction with CEO and CAA was also described. Animal-based risk estimates for human inhalation exposure to VC using total metabolism estimates from the PBPK model were consistent with risk estimates based on human epidemiological data, and were lower than those currently used in environmental decision-making by a factor of 80.

A hybrid computational fluid dynamics and physiologically based pharmacokinetic model for comparison of predicted tissue concentrations of acrylic acid and other vapors in the rat and human nasal cavities following inhalation exposure

To assist in interspecies dosimetry comparisons for risk assessment of the nasal effects of organic acids, a hybrid computational fluid dynamics (CFD) and physiologically based pharmacokinetic (PBPK) dosimetry model was constructed to estimate the regional tissue dose of inhaled vapors in the rat and human nasal cavity. Application to a specific vapor would involve the incorporation of the chemical-specific reactivity, metabolism, partition coefficients, and diffusivity (in both air and tissue phases) of the vapor. This report describes the structure of the CFD-PBPK model and its application to a representative acidic vapor, acrylic acid, for interspecies tissue concentration comparisons to assist in risk assessment. By using the results from a series of short-term in vivo studies combined with computer modeling, regional nasal tissue dose estimates were developed and comparisons of tissue doses between species were conducted. To make these comparisons, the assumption was made that the susceptibilities of human and rat olfactory epithelium to the cytotoxic effects of organic acids were similar, based on similar histological structure and common mode of action considerations. Interspecies differences in response were therefore assumed to be driven primarily by differences in nasal tissue concentrations that result from regional differences in nasal air flow patterns relative to the species-specific distribution of olfactory epithelium in the nasal cavity. The results of simulations with the seven-compartment CFD-PBPK model suggested that the olfactory epithelium of the human nasal cavity would be exposed to tissue concentrations of acrylic acid similar to that of the rat nasal cavity when the exposure conditions are the same. Similar analysis of CFD data and CFD-PBPK model simulations with a simpler one-compartment model of the whole nasal cavities of rats and humans provides comparable results to averaging over the compartments of the seven-compartment model. These results indicate that the general structure of the hybrid CFD-PBPK model applied in this assessment would be useful for target tissue dosimetry and interspecies dose comparisons for a wide variety of vapors. Because of its flexibility, this CFD-PBPK model is envisioned to be a platform for the construction of case-specific inhalation dosimetry models to simulate in vivo exposures that do not involve significant histopathological damage to the nasal cavity.

Determination of a site-specific reference dose for methylmercury for fish-eating populations

Environmental risk-management decisions in the U.S. involving potential exposures to methylmercury currently use a reference dose (RfD) developed by the U.S. Environmental Protection Agency (USEPA). This RfD is based on retrospective studies of an acute poisoning incident in Iraq in which grain contaminated with a methylmercury fungicide was inadvertently used in the baking of bread. The exposures, which were relatively high but lasted only a few months, were associated with neurological effects in both adults (primarily paresthesia) and infants (late walking, late talking, etc.). It is generally believed that the developing fetus represents a particularly sensitive subpopulation for the neurological effects of methylmercury. The USEPA derived an RfD of 0.1 microg/kg/day based on benchmark dose (BMD) modeling of the combined neurological endpoints reported for children exposed in utero. This RfD included an uncertainty factor of 10 to consider human pharmacokinetic variability and database limitations (lack of data on multigeneration effects or possible long-term sequelae of perinatal exposure). Alcoa signed an Administrative Order of Consent for the conduct of a remedial investigation/feasibility study (RI/FS) at their Point Comfort Operations and the adjacent Lavaca Bay in Texas to address the effects of historical discharges of mercury-containing wastewater. In cooperation with the Texas Natural Resource Conservation Commission and USEPA Region VI, Alcoa conducted a baseline risk assessment to assess potential risk to human health and the environment. As a part of this assessment. Alcoa pursued the development of a site-specific RfD for methylmercury to specifically address the potential human health effects associated with the ingestion of contaminated finfish and shellfish from Lavaca Bay. Application of the published USEPA RfD to this site is problematic; while the study underlying the RfD represented acute exposure to relatively high concentrations of methylmercury, the exposures of concern for the Point Comfort site are from the chronic consumption of relatively low concentrations of methylmercury in fish. Since the publication of the USEPA RfD, several analyses of chronic exposure to methylmercury in fish-eating populations have been reported. The purpose of the analysis reported here was to evaluate the possibility of deriving an RfD for methylmercury, specifically for the case of fish ingestion, on the basis of these new studies. In order to better support the risk-management decisions associated with developing a remediation approach for the site in question, the analysis was designed to provide information on the distribution of acceptable ingestion rates across a population, which could reasonably be expected to be consistent with the results of the epidemiological studies of other fish-eating populations. Based on a review of the available literature on the effects of methylmercury, a study conducted with a population in the Seychelles Islands was selected as the critical study for this analysis. The exposures to methylmercury in this population result from chronic, multigenerational ingestion of contaminated fish. This prospective study was carefully conducted and analyzed, included a large cohort of mother-infant pairs, and was relatively free of confounding factors. The results of this study are essentially negative, and a no-observed-adverse-effect level (NOAEL) derived from the estimated exposures has recently been used by the Agency for Toxic Substances and Disease Registry (ATSDR) as the basis for a chronic oral minimal risk level (MRL) for methylmercury. In spite of the fact that no statistically significant effects were observed in this study, the data as reported are suitable for dose-response analysis using the BMD method. Evaluation of the BMD method used in this analysis, as well as in the current USEPA RfD, has demonstrated that the resulting 95% lower bound on the 10% benchmark dose (BMDL) represents a conservative estimate of the traditional NOAEL, and that it is superior to the use of "average" or "grouped" exposure estimates when dose-response information is available, as is the case for the Seychelles study. A more recent study in the Faroe Islands, which did report statistically significant associations between methylmercury exposure and neurological effects, could not be used for dose-response modeling due to inadequate reporting of the data and confounding from co-exposure to polychlorinated biphenyls (PCBs). BMD modeling over the wide range of neurological endpoints reported in the Seychelles study yielded a lowest BMDL for methylmercury in maternal hair of 21 ppm. This BMDL was then converted to an expected distribution of daily ingestion rates across a population using Monte Carlo analysis with a physiologically based pharmacokinetic (PBPK) model to evaluate the impact of interindividual variability. The resulting distribution of ingestion rates at the BMDL had a geometric mean of 1.60 microg/kg/day with a geometric standard deviation of 1.33; the 1st, 5th, and 10th percentiles of the distribution were 0.86, 1.04, and 1.15 microg/kg/day. In place of the use of an uncertainty factor of 3 for pharmacokinetic variability, as is done in the current RfD, one of these lower percentiles of the daily ingestion rate distribution provides a scientifically based, conservative basis for taking into consideration the impact of pharmacokinetic variability across the population. On the other hand, it was felt that an uncertainty factor of 3 for database limitations should be used in the current analysis. Although there can be high confidence in the benchmark-estimated NOAEL of 21 ppm in the Seychelles study, some results in the New Zealand and Faroe Islands studies could be construed to suggest the possibility of effects at maternal hair concentrations below 10 ppm. In addition, while concerns regarding the possibility of chronic sequelae are not supported by the available data, neither can they be absolutely ruled out. The use of an uncertainty factor of 3 is equivalent to using a NOAEL of 7 ppm in maternal hair, which provides additional protection against the possibility that effects could occur at lower concentrations in some populations. Based on the analysis described above, the distribution of acceptable daily ingestion rates (RfDs) recommended to serve as the basis for site-specific risk-management decisions at Alcoa's Point Comfort Operations ranges from approximately 0.3 to 1.1 microg/kg/day, with a population median (50th percentile) of 0.5 microg/kg/day. By analogy with USEPA guidelines for the use of percentiles in applications of distributions in exposure assessments, the 10th percentile provides a reasonably conservative measure. On this basis, a site-specific RfD of 0.4 microg/kg/day is recommended.