Dosimetric adjustment factors for methyl methacrylate derived from a steady-state analysis of a physiologically based clearance-extraction model

Development of a Gestational and Lactational Physiologically Based Pharmacokinetic (PBPK) Model for Perfluorooctane Sulfonate (PFOS) in Rats and Humans and Its Implications in the Derivation of Health-Based Toxicity Values

BACKGROUND

There is a great concern on potential adverse effects of exposure to perfluorooctane sulfonate (PFOS) in sensitive subpopulations, such as pregnant women, fetuses, and neonates, due to its reported transplacental and lactational transfer and reproductive and developmental toxicities in animals and humans.

OBJECTIVES

This study aimed to develop a gestational and lactational physiologically based pharmacokinetic (PBPK) model in rats and humans for PFOS to aid risk assessment in sensitive human subpopulations.

METHODS

Based upon existing PBPK models for PFOS, the present model addressed a data gap of including a physiologically based description of basolateral and apical membrane transporter-mediated renal reabsorption and excretion in kidneys during gestation and lactation. The model was calibrated with published rat toxicokinetic and human biomonitoring data and was independently evaluated with separate data. Monte Carlo simulation was used to address the interindividual variability.

RESULTS

Model simulations were generally within 2-fold of observed PFOS concentrations in maternal/fetal/neonatal plasma and liver in rats and humans. Estimated fifth percentile human equivalent doses (HEDs) based on selected critical toxicity studies in rats following U.S. Environmental Protection Agency (EPA) guidelines ranged from 0.08 to 0.91 μ g / kg  per day . These values are lower than the HEDs estimated in U.S. EPA guidance (0.51 - 1.6 μ g / kg  per day ) using an empirical toxicokinetic model in adults.

CONCLUSIONS

The results support the importance of renal reabsorption/excretion during pregnancy and lactation in PFOS dosimetry and suggest that the derivation of health-based toxicity values based on developmental toxicity studies should consider gestational/lactational dosimetry estimated from a life stage-appropriate PBPK model. This study provides a quantitative tool to aid risk reevaluation of PFOS, especially in sensitive human subpopulations, and it provides a basis for extrapolating to other per- and polyfluoroalkyl substances (PFAS). All model codes and detailed tutorials are provided in the Supplemental Materials to allow readers to reproduce our results and to use this model. https://doi.org/10.1289/EHP7671.

Bayesian evaluation of a physiologically based pharmacokinetic (PBPK) model for perfluorooctane sulfonate (PFOS) to characterize the interspecies uncertainty between mice, rats, monkeys, and humans: Development and performance verification

A challenge in the risk assessment of perfluorooctane sulfonate (PFOS) is the large interspecies differences in its toxicokinetics that results in substantial uncertainty in the dosimetry and toxicity extrapolation from animals to humans. To address this challenge, the objective of this study was to develop an open-source physiologically based pharmacokinetic (PBPK) model accounting for species-specific toxicokinetic parameters of PFOS. Considering available knowledge about the toxicokinetic properties of PFOS, a PBPK model for PFOS in mice, rats, monkeys, and humans after intravenous and oral administrations was created. Available species-specific toxicokinetic data were used for model calibration and optimization, and independent datasets were used for model evaluation. Bayesian statistical analysis using Markov chain Monte Carlo (MCMC) simulation was performed to optimize the model and to characterize the uncertainty and interspecies variability of chemical-specific parameters. The model predictions well correlated with the majority of datasets for all four species, and the model was validated with independent data in rats, monkeys, and humans. The model was applied to predict human equivalent doses (HEDs) based on reported points of departure in selected critical toxicity studies in rats and monkeys following U.S. EPA's guidelines. The lower bounds of the model-derived HEDs were overall lower than the HEDs estimated by U.S. EPA (e.g., 0.2 vs. 1.3 μg/kg/day based on the rat plasma data). This integrated and comparative analysis provides an important step towards improving interspecies extrapolation and quantitative risk assessment of PFOS, and this open-source model provides a foundation for developing models for other perfluoroalkyl substances.

A hybrid CFD-PBPK model for naphthalene in rat and human with IVIVE for nasal tissue metabolism and cross-species dosimetry

A PBPK model for naphthalene in the rat and human that incorporates a hybrid CFD-PBPK description of the upper respiratory tract was developed to support cross-species dosimetry comparisons of naphthalene concentrations and tissue normalized rate of metabolism in the nasal respiratory and olfactory epithelium, lung and liver. In vitro measurements of metabolic rates from microsomal incubations published for rat and monkey (surrogate for human) were scaled to the specific tissue based on the tissue microsomal content and volume of tissue. The model reproduces time courses for naphthalene blood concentrations from intravenous and inhalation exposures in rats and upper respiratory tract extraction data in both naïve rats and rats pre-treated to inhibit nasal metabolism. This naphthalene model was applied to estimate human equivalent inhalation concentrations (HECs) corresponding to several NOAELs or LOAELs for the non-cancer effects of naphthalene in rats. Two approaches for cross-species extrapolation were compared: (1) equivalence based on tissue naphthalene concentration and (2) equivalence based on amount metabolized per minute (normalized to tissue volume). At the NOAEL of 0.1 ppm, the regional gas dosimetry ratio (RGDR) based on naphthalene concentration was 0.18 for the dorsal olfactory region; however, the RGDR rises to 5.4 when based on the normalized amount metabolized due to the lower of expression of CYP isozymes in the nasal epithelium of primates and humans. The resulting HEC is 0.12 ppm (0.63 mg/m(3)) continuous exposure at the rat NOAEL of 0.1 ppm (6 h/day, 5 days/week).

Sex differences in the acute nasal antioxidant/antielectrophilic response of the rat to inhaled naphthalene

Naphthalene is a nasal carcinogen, inducing respiratory adenomas in male and olfactory neuroblastomas in female rats, respectively. The reasons for the site and sex-specific tumorigenic response are unknown. Naphthalene is bioactivated to electrophilic metabolites; cytotoxicity followed by regenerative cell proliferation is likely involved in the tumorigenic response. To examine sex differences in the acute nasal response to naphthalene, male and female F344 rats were nose-only exposed to 0, 1, 3, 10, or 30 ppm naphthalene vapor for 4 or 6 h. Following exposure, respiratory/transitional mucosa (RTM) and olfactory mucosa (OM) were isolated and analyzed for markers of oxidant/electrophilic stress and/or toxicity, including reduced/oxidized glutathione levels (GSH/GSSG), mRNA levels of electrophile-responsive genes, and epithelial cytoxicity (as measured by membrane permeability to ethidium homodimer-1). Naphthalene caused significant depletion of GSH in RTM and OM with no increase in GSSG. Cytotoxicity was apparent at concentrations of 15 and 30 ppm. No consistent sex differences were observed in these responses. Sex differences were observed in the induction of antielectrophilic genes in OM: glutamyl cysteine ligase (catalytic subunit) (Gclc), NADPH quinone oxidase 1 (Nqo1), and heme oxygenase 1 (Hmox1) were all induced to a greater extent in the male OM compared with the female. No consistent sex differences were observed in the RTM. Although the mechanism of the sex difference in the RTM adenoma response remains enigmatic, sex differences in the induction of antioxidant/electrophile-responsive genes may contribute to the heightened sensitivity of the female OM to the carcinogenic effects of naphthalene.

Hypothesis-based weight-of-evidence evaluation of methyl methacrylate olfactory effects in humans and derivation of an occupational exposure level

Over 40 years of scientific evidence indicates that methyl methacrylate (MMA) causes olfactory effects in rodents that are relevant to humans. More recent scientific studies have focused on understanding the apparent lack of species concordance between the rodent and human studies. Toxicokinetic studies and a physiologically based pharmacokinetic (PBPK) model describing inhalation dosimetry of MMA in the upper respiratory tract (URT) of rats and humans point to differences in nasal morphology and biochemistry that could explain and reconcile these differences as species-specific manifestations of a common toxicological process. We have applied the hypothesis-based weight-of-evidence (HBWoE) approach to evaluate the concordance of the available data and the hypothesis that the observed difference in sensitivity between rats and humans may be the expected result of physiological and biochemical differences. Our WoE analysis indicates that when the several lines of evidence (i.e., animal, human, mode-of-action, and toxicokinetics data) are integrated, they inform interpretation of one another and, overall, support use of the human data for derivation of an MMA occupational exposure level (OEL) of 50 ppm.

Physiologically-based pharmacokinetic modeling for absorption, transport, metabolism and excretion

The seminal paper on the liver physiologically-based pharmacokinetic (PBPK) model by Rowland et al. (J Pharmacokinet Biopharm 1:123-136, 1973) that described the influence of blood flow, intrinsic clearance, and binding on hepatic clearance had inspired further development of PBPK modeling of the liver, kidney and intestine as well as whole body. Shortly thereafter, a series of papers from Pang and Rowland compared the well-stirred and parallel-tube liver models and sparked further development on clearance concepts in the liver, including those described by the dispersion model. From 2005 onwards, several seminal papers by Rodgers and Rowland, in their recognition of the binding of molecules to tissue acidic and neutral phospholipids, improved the methodology in providing estimates of the tissue-to-plasma coefficient and rendering easy calculation of these hard-to-get constants. The improvement has strongly consolidated the basic premise on PBPK modeling and simulations and these basics have allowed scientists to focus on other important variables: membrane barriers, and transporter and enzyme and their heterogeneities that further impact drug disposition. In particular, the PBPK models have delved into sequential metabolism and futile cycling to illustrate how transporters and enzymes could affect the metabolism of drugs and metabolites. PBPK models that are especially pertinent to metabolite kinetics are being utilized in drug studies and risk assessment. These types of PBPK modeling reveal differences in kinetics between the formed vs. preformed metabolite, showing special considerations for membrane barriers, and the influence of competing pathways and competing organs.

Magnetic resonance imaging and computational fluid dynamics (CFD) simulations of rabbit nasal airflows for the development of hybrid CFD/PBPK models

The percentages of total airflows over the nasal respiratory and olfactory epithelium of female rabbits were calculated from computational fluid dynamics (CFD) simulations of steady-state inhalation. These airflow calculations, along with nasal airway geometry determinations, are critical parameters for hybrid CFD/physiologically based pharmacokinetic models that describe the nasal dosimetry of water-soluble or reactive gases and vapors in rabbits. CFD simulations were based upon three-dimensional computational meshes derived from magnetic resonance images of three adult female New Zealand White (NZW) rabbits. In the anterior portion of the nose, the maxillary turbinates of rabbits are considerably more complex than comparable regions in rats, mice, monkeys, or humans. This leads to a greater surface area to volume ratio in this region and thus the potential for increased extraction of water soluble or reactive gases and vapors in the anterior portion of the nose compared to many other species. Although there was considerable interanimal variability in the fine structures of the nasal turbinates and airflows in the anterior portions of the nose, there was remarkable consistency between rabbits in the percentage of total inspired airflows that reached the ethmoid turbinate region (approximately 50%) that is presumably lined with olfactory epithelium. These latter results (airflows reaching the ethmoid turbinate region) were higher than previous published estimates for the male F344 rat (19%) and human (7%). These differences in regional airflows can have significant implications in interspecies extrapolations of nasal dosimetry.

Deposition of inhaled nanoparticles in the rat nasal passages: dose to the olfactory region

In vivo experiments have shown that nanoparticles depositing in the rat olfactory region can translocate to the brain via the olfactory nerve. Quantitative predictions of the dose delivered by inhalation to the olfactory region are needed to clarify this route of exposure and to evaluate the dose-response effects of exposure to toxic nanoparticles. Previous in vivo and in vitro studies quantified the percentage of inhaled nanoparticles that deposit in the rat nasal passages, but olfactory dose was not determined. The dose to specific nasal epithelium types is expected to vary with inhalation rate and particle size. The purpose of this investigation, therefore, was to develop estimates of nanoparticle deposition in the nasal and, more specifically, olfactory regions of the rat. A three-dimensional, anatomically accurate, computational fluid dynamics (CFD) model of the rat nasal passages was employed to simulate inhaled airflow and to calculate nasal deposition efficiency. Particle sizes from 1 to 100 nm and airflow rates of 288, 432, and 576 ml/min (1, 1.5, and 2 times the estimated resting minute volume) were simulated. The simulations predicted that olfactory deposition is maximum at 6-9% of inhaled material for 3- to 4-nm particles. The spatial distribution of deposited particles was predicted to change significantly with particle size, with 3-nm particles depositing mostly in the anterior nose, while 30-nm particles were more uniformly distributed throughout the nasal passages.

Investigations of the use of bioavailability data to adjust occupational exposure limits for active pharmaceutical ingredients

Occupational exposure limits (OELs) for active pharmaceutical ingredients have traditionally been established using no-observed-adverse-effect levels derived from clinical studies employing po and iv routes of administration and by applying default uncertainty factors or chemical-specific adjustment factors. However, exposure by the inhalation or dermal route is more relevant in terms of occupational safety. In this investigation, to explore new methods for route-to-route extrapolation, the bioavailability of MK-0679, a leukotriene D(4) receptor antagonist, was compared following iv, po, intranasal (in), or intratracheal (it) administration. The relative bioavailability of MK-0679 was iv congruent with it > po congruent with in. Bioavailability correction factors (BCFs) of 2.0 and 0.6 were derived from these data to adjust a hypothetical OEL of 0.1 mg/m(3) for MK-0679 with particle sizes of 10 and 50 mum, respectively. These BCFs were used to adjust the OEL established using po clinical data, to reflect the differences in bioavailability following deposition in different regions of the respiratory tract. To further investigate how bioavailability data could be used in setting OELs, a preliminary pharmacokinetic (PK) model was developed to describe the time course of plasma concentrations using the data from the route comparison study. An inhalation study was then performed to test the validity of using either empirical data or modeling approaches to derive BCFs when setting OELs. These investigations demonstrated how the use of route-specific PK data could reduce some of the uncertainties associated with route-to-route extrapolation and allow for improved precision and quantitative adjustments when establishing OELs. Further investigations are needed to better understand the factors responsible for differences in systemic uptake following deposition in different regions of the respiratory tract and how these can be generalized across different classes of soluble compounds.

Derivation of an inhalation reference concentration based upon olfactory neuronal loss in male rats following subchronic acetaldehyde inhalation

Acetaldehyde inhalation induces neoplastic and nonneoplastic responses in the rodent nasal cavity. This experiment further characterizes the dose-response relationship for nasal pathology, nasal epithelial cell proliferation, and DNA-protein cross-link formation in F-344 rats exposed subchronically to acetaldehyde. Animals underwent whole-body exposure to 0, 50, 150, 500, or 1500 ppm acetaldehyde for 6 h/day, 5 days/wk for up to 65 exposure days. Respiratory tract histopathology was evaluated after 4, 9, 14, 30, and 65 exposure days. Acetaldehyde exposure was not associated with reduced body weight gain or other evidence of systemic toxicity. Histologic evaluation of the nasal cavity showed an increased incidence of olfactory neuronal loss (ONL) following acute to subchronic exposure to > or = 150 ppm acetaldehyde and increased olfactory epithelial cell proliferation following exposure to 1500 ppm acetaldehyde. The severity of the ONL demonstrated dose- and temporal-dependent behaviors, with minimal effects noted at 150-500 ppm acetaldehyde and moderately severe lesions seen in the highest exposure group, with increased lesion severity and extent as the exposure duration increased. Acetaldehyde exposure was also associated with inflammation, hyperplasia, and squamous metaplasia of the respiratory epithelium. These responses were seen in animals exposed to > or = 500 ppm acetaldehyde. Acetaldehyde exposure was not associated with increased DNA-protein cross-link formation in the respiratory or olfactory epithelium. A model of acetaldehyde pharmacokinetics in the nose was used to derive an inhalation reference concentration (RfC) of 0.4 ppm, based on the no-observed-adverse-effect level (NOAEL) of 50 ppm for the nasal pathology seen in this study.

A PBPK model for evaluating the impact of aldehyde dehydrogenase polymorphisms on comparative rat and human nasal tissue acetaldehyde dosimetry

Acetaldehyde is an important intermediate in the chemical synthesis and normal oxidative metabolism of several industrially important compounds, including ethanol, ethyl acetate, and vinyl acetate. Chronic inhalation of acetaldehyde leads to degeneration of the olfactory and respiratory epithelium in rats at concentrations > 50 ppm (90 day exposure) and respiratory and olfactory nasal tumors at concentrations > or = 750 ppm, the lowest concentration tested in the 2-yr chronic bioassay. Differences in the anatomy and biochemistry of the rodent and human nose, including polymorphisms in human high-affinity acetaldehyde dehydrogenase (ALDH2), are important considerations for interspecies extrapolations in the risk assessment of acetaldehyde. A physiologically based pharmacokinetic model of rat and human nasal tissues was constructed for acetaldehyde to support a dosimetry-based risk assessment for acetaldehyde (Dorman et al., 2008). The rodent model was developed using published metabolic constants and calibrated using upper-respiratory-tract acetaldehyde extraction data. The human nasal model incorporates previously published tissue volumes, blood flows, and acetaldehyde metabolic constants. ALDH2 polymorphisms were represented in the human model as reduced rates of acetaldehyde metabolism. Steady-state dorsal olfactory epithelial tissue acetaldehyde concentrations in the rat were predicted to be 409, 6287, and 12,634 microM at noncytotoxic (50 ppm), and cytotoxic/tumorigenic exposure concentrations (750 and 1500 ppm), respectively. The human equivalent concentration (HEC) of the rat no-observed-adverse-effect level (NOAEL) of 50 ppm, based on steady-state acetaldehyde concentrations from continual exposures, was 67 ppm. Respiratory and olfactory epithelial tissue acetaldehyde and H(+) (pH) concentrations were largely linear functions of exposure in both species. The impact of presumed ALDH2 polymorphisms on human olfactory tissue concentrations was negligible; the high-affinity, low-capacity ALDH2 does not contribute significantly to acetaldehyde metabolism in the nasal tissues. The human equivalent acetaldehyde concentration for homozygous low activity was 66 ppm, 1.5% lower than for the homozygous full activity phenotype. The rat and human acetaldehyde PBPK models developed here can also be used as a bridge between acetaldehyde dose-response and mode-of-action data as well as between similar databases for other acetaldehyde-producing nasal toxicants.

Modeling approaches for estimating the dosimetry of inhaled toxicants in children

Risk assessment of inhaled toxicants has typically focused upon adults, with modeling used to extrapolate dosimetry and risks from lab animals to humans. However, behavioral factors such as time spent playing outdoors may lead to more exposure to inhaled toxicants in children. Depending on the inhaled agent and the age and size of the child, children may receive a greater internal dose than adults because of greater ventilation rate per body weight or lung surface area, or metabolic differences may result in different tissue burdens. Thus, modeling techniques need to be adapted to children in order to estimate inhaled dose and risk in this potentially susceptible life stage. This paper summarizes a series of inhalation dosimetry presentations from the U.S. EPA's Workshop on Inhalation Risk Assessment in Children held on June 8-9, 2006 in Washington, DC. These presentations demonstrate how existing default models for particles and gases may be adapted for children, and how more advanced modeling of toxicant deposition and interaction in respiratory airways takes into account children's anatomy and physiology. These modeling efforts identify child-adult dosimetry differences in respiratory tract regions that may have implications for children's vulnerability to inhaled toxicants. A decision framework is discussed that considers these different approaches and modeling structures including assessment of parameter values, supporting data, reliability, and selection of dose metrics.

mRNA-induction and cytokine release during in vitro exposure of human nasal respiratory epithelia to methyl methacrylate

BACKGROUND

Methyl methacrylate (MMA) has been reported to cause histopathological changes in rodent nasal epithelium after inhalation challenges. Data in humans are lacking.

METHODS

In this in vitro design 22 primary cell cultures taken from inferior turbinate tissue of healthy individuals were exposed to MMA concentrations of 50 ppm (German MAK-value) and 200 ppm. mRNA expression and cytokine release of inflammatory mediators were quantified after 4h and after 24h. Controls were exposed to synthetic air. Q-PCR analysis was performed for TNF-alpha, IL-1beta, IL-6, IL-8, MCP-1, GMCSF, Cox-1 and Cox-2. ELISA assays were performed from culture supernatants for TNF-alpha, IL-1beta, IL-6, IL-8, MCP-1 and GMCSF.

RESULTS

Acute inductions of mRNA after 4h were observed for TNF-alpha, IL-1beta, IL-6, IL-8 and MCP-1 at 50 ppm. ELISA analysis of the described parameters did not reveal any significant upregulations at both concentrations after both 4h and 24h.

CONCLUSIONS

The obtained data suggest that exposure of human respiratory epithelia in vitro to 50 ppm and to 200 ppm of MMA does not induce lasting upregulation of the inflammatory mediators measured in this study. The exposure limit of 50 ppm appears safe following these results obtained from human respiratory epithelia.

A physiological toxicokinetic model for inhaled propylene oxide in rat and human with special emphasis on the nose

Chronic exposure to high concentrations of PO induced inflammation in the respiratory nasal mucosa (RNM) of rodents and, for concentrations >or= 300 ppm, caused nasal tumors. Considering the nose to be the most relevant target organ for PO-induced tumorigenicity, we developed a physiological toxicokinetic model for PO in rats and humans. It includes compartments for arterial, venous, and pulmonary blood, liver, muscle, fat, richly perfused tissues, lung, and nose. It simulates inhalation of PO, its distribution into tissues by blood flow, and its elimination by exhalation and metabolism. In nose, lung, and liver of rats, PO conjugation with glutathione (GSH), PO-induced GSH depletion, and formation of PO adducts to DNA are described. Also modeled are PO adducts to hemoglobin of rats and humans. Required partition coefficients and metabolic parameters were derived experimentally or from publications. In rats, simulated PO concentrations in blood and GSH levels in tissues agreed with measured data. If compared with reported values, levels of adducts with hemoglobin were underpredicted up to a factor of about 2. Adducts with DNA differed up to a factor of 3. Hemoglobin adducts predicted for PO-exposed workers were 1.5-1.9 times higher than the reported ones. Considering identical conditions of PO exposure, similar PO concentrations in RNM were modeled for rats and humans. Also, PO concentrations in blood, about 1/30th of those in RNM, were similar in both species. Since the model was evaluated on all available data in rats and humans, we consider it to be useful for estimating the risk from inhalation exposure to PO.

Incorporation of Tissue Reaction Kinetics in a Computational Fluid Dynamics Model for Nasal Extraction of Inhaled Hydrogen Sulfide in Rats

Rodents exposed to hydrogen sulfide (H2S) develop olfactory neuronal loss. This lesion has been used by the risk assessment community to develop occupational and environmental exposure standards. A correlation between lesion locations and areas of high H2S flux to airway walls has been previously demonstrated, but a quantitative dose assessment is needed to extrapolate dose at lesion sites to humans. In this study, nasal extraction (NE) of 10, 80, and 200 ppm H2S was measured in the isolated upper respiratory tract of anesthetized rats under constant unidirectional inspiratory flow rates of 75, 150, and 300 ml/min. NE was dependent on inspired H2S concentration and air flow rate: increased NE was observed when H2S exposure concentrations or inspiratory air flow rates were low. An anatomically accurate, three-dimensional computational fluid dynamics (CFD) model of rat nasal passages was used to predict NE of inhaled H2S. To account for the observed dependence of NE on H2S exposure concentration, the boundary condition used at airway walls incorporated first-order and saturable kinetics in nasal tissue to govern mass flux at the air:tissue interface. Since the kinetic parameters cannot be obtained using the CFD model, they were estimated independently by fitting a well-mixed, two-compartment pharmacokinetic (PK) model to the NE data. Predicted extraction values using this PK-motivated CFD approach were in good agreement with the experimental measurements. The CFD model provides estimates of localized H2S flux to airway walls and can be used to calibrate lesion sites by dose.

Dose-Incidence Modeling: Consequences of Linking Quantal Measures of Response to Depletion of Critical Tissue Targets

In developing mechanistic PK-PD models, incidence of toxic responses in a population has to be described in relation to measures of biologically effective dose (BED). We have developed a simple dose-incidence model that links incidence with BED for compounds that cause toxicity by depleting critical cellular target molecules. The BED in this model was the proportion of target molecule adducted by the dose of toxic compound. Our modeling approach first estimated the proportion depleted for each dose and then calculated the tolerance distribution for toxicity in relation to either administered dose or log of administered dose. We first examined cases where the mean of the tolerance distribution for toxicity occurred when a significant proportion of target had been adducted (i.e., more than half). When a normal distribution was assumed to exist for the relationship of incidence and BED, the tolerance distribution based on administered dose for these cases becomes asymmetrical and logarithmic transformations of the administered dose axis lead to a more symmetrical distribution. These linked PK-PD models for tissue reactivity, consistent with conclusions from other work for receptor binding models (Lutz et al., 2005), indicate that log normal distributions with administered dose may arise from normal distributions for BED and nonlinear kinetics between BED and administered dose. These conclusions are important for developing biologically based dose response (BBDR) models that link incidences of toxicity or other biological responses to measures of BED.

Validation of human physiologically based pharmacokinetic model for vinyl acetate against human nasal dosimetry data

Vinyl acetate has been shown to induce nasal lesions in rodents in inhalation bioassays. A physiologically based pharmacokinetic (PBPK) model for vinyl acetate has been used in human risk assessment, but previous in vivo validation was conducted only in rats. Controlled human exposures to vinyl acetate were conducted to provide validation data for the application of the model in humans. Five volunteers were exposed to 1, 5, and 10 ppm 13C1,13C2 vinyl acetate via inhalation. A probe inserted into the nasopharyngeal region sampled both 13C1,13C2 vinyl acetate and the major metabolite 13C1,13C2 acetaldehyde during rest and light exercise. Nasopharyngeal air concentrations were analyzed in real time by ion trap mass spectrometry (MS/MS). Experimental concentrations of both vinyl acetate and acetaldehyde were then compared to predicted concentrations calculated from the previously published human model. Model predictions of vinyl acetate nasal extraction compared favorably with measured values of vinyl acetate, as did predictions of nasopharyngeal acetaldehyde when compared to measured acetaldehyde. The results showed that the current PBPK model structure and parameterization are appropriate for vinyl acetate. These analyses were conducted from 1 to 10 ppm vinyl acetate, a range relevant to workplace exposure standards but which would not be expected to saturate vinyl acetate metabolism. Risk assessment based on this model further concluded that 24 h per day exposures up to 1 ppm do not present concern regarding cancer or non-cancer toxicity. Validation of the vinyl acetate human PBPK model provides support for these conclusions.

Physiologically-based kinetic modeling of vapours toxic to the respiratory tract

The respiratory tract is frequently identified as a site of toxicity for inhaled xenobiotic chemicals. Usually, these observations come from controlled animal studies. For these studies to be of quantitative value to human health risk assessment, species-specific factors governing dosimetry of inhaled substances must be taken into account. Toxicokinetics of vapours in the respiratory tract are defined by absorption, distribution, metabolism, and excretion, as they are in other tissues; however, these concepts take on new dimensions when considering respiratory tract toxicants, especially those that elicit portal of entry effects by directly interacting with the tissue lining the respiratory tract. Species-specific factors related to anatomy, physiology and biochemistry govern inter-species extrapolation of toxicokinetics. This article discusses critical factors of respiratory tract kinetics that should be considered when developing physiological-based toxicokinetic (PBTK) models for inhaled vapours. Important considerations such as impact of regional airflow-delivery, water solubility, reactivity, and rates of local biotransformation on respiratory tract tissue dosimetry are highlighted. These factors can be accounted for only to a limited extent when using default approaches to extrapolate dosimetry of inhaled substances across species. On the other hand, PBTK modeling has the flexibility to accommodate many of the critical determinants of respiratory tract toxicity. PBTK models can also help identify the most critical toxicokinetic data necessary to replace defaults. PBTK approaches have led to more informed estimates of human target tissue dose, and therefore human health risk, especially where these risk assessments have been based on extrapolation of animal dosimetry studies. Experience derived from the development of more intensive case studies have, in turn, enabled simplified approaches to the use of PBTK modeling for respiratory tract toxicants. Whether simplified or highly complex, PBTK modeling approaches are proven to be of great utility to risk assesors interested in applying quantitative information to informed risk assessment evaluations.

Toxicokinetic modeling and its applications in chemical risk assessment

In recent years physiologically based pharmacokinetic (PBPK) modeling has found frequent application in risk assessments where PBPK models serve as important adjuncts to studies on modes of action of xenobiotics. In this regard, studies on mode of action provide insight into both the sites/mechanisms of action and the form of the xenobiotic associated with toxic responses. Validated PBPK models permit calculation of tissue doses of xenobiotics and metabolites for a variety of conditions, i.e. at low-doses, in different animal species, and in different members of a human population. In this manner, these PBPK models support the low-dose and interspecies extrapolations that are important components of current risk assessment methodologies. PBPK models are sometimes referred to as physiological toxicokinetic (PT) models to emphasize their application with compounds causing toxic responses. Pharmacokinetic (PK) modeling in general has a rich history. Data-based PK compartmental models were developed in the 1930's when only primitive tools were available for solving sets of differential equations. These models were expanded in the 1960's and 1970's to accommodate new observations on dose-dependent elimination and flow-limited metabolism. The application of clearance concepts brought many new insights about the disposition of drugs in the body. In the 1970's PBPK/PT models were developed to evaluate metabolism of volatile compounds of occupational importance, and, for the first time, dose-dependent processes in toxicology were included in PBPK models in order to assess the conditions under which saturation of metabolic and elimination processes lead to non-linear dose response relationships. In the 1980's insights from chemical engineers and occupational toxicology were combined to develop PBPK/PT models to support risk assessment with methylene chloride and other solvents. The 1990's witnessed explosive growth in risk assessment applications of PBPK/PT models and in applying sensitivity and variability methods to evaluate model performance. Some of the compounds examined in detail include butadiene, styrene, glycol ethers, dioxins and organic esters/aids. This paper outlines the history of PBPK/PT modeling, emphasizes more recent applications of PBPK/TK models in health risk assessment, and discusses the risk assessment perspective provided by modern uses of these modeling approaches.

Toxicokinetics and physiologically based toxicokinetics in toxicology and risk assessment

Toxicokinetics is the study of kinetics of absorption, distribution, metabolism, and excretion of a xenobiotic under the conditions of toxicity evaluation. Conventional toxicokinetics uses the hypothetical compartments, and the model is composed of rate equations that describe the time course of drug and chemical disposition. The utility of toxicokinetics in toxicity evaluation and interpretation of animal toxicology data is emerging as an important tool in product discovery and development. With implementation of the International Conference on Harmonization (ICH) guidelines on systemic exposure and dose selection, toxicokinetics have been integrated in routine toxicity evaluations. Although traditional compartmental/noncompartmental models are generally adequate for assessing systemic exposure, they are unable to the predict time course of drug disposition in target tissues and often fail to relate systemic drug levels to a biological response. Physiologically based toxicokinetic (PB-TK) models address this deficiency of traditional compartmental models. PB-TK models are the kinetic models of the uptake and disposition of chemicals based on rates of biochemical reactions, physiological and anatomical characteristics. These models, when developed appropriately, can predict target organ drug distribution in different species under variety of conditions. This minireview discusses the basic principles, and applications of traditional compartmental toxicokinetic and physiologically based toxicokinetics (PB-TK) models in drug development and risk assessment. Special emphasis will be placed on discussion related to interpretation of the ICH guidelines related to toxicokinetics and the utility of toxicokinetics data in dose selection for toxicity and carcinogenicity studies. The utility of PB-TK models in risk assessment of methylene chloride, vinyl chloride, retinoic acid, dioxin, and inhaled organic esters is discussed.

Physiologically based pharmacokinetic (PBPK) models for nasal tissue dosimetry of organic esters: assessing the state-of-knowledge and risk assessment applications with methyl methacrylate and vinyl acetate

Mathematical models have been developed to describe nasal epithelial tissue dosimetry with two compounds, vinyl acetate (VA) and methyl methacrylate (MMA), that cause toxicity in these tissues These models couple computational fluid dynamics (CFD) calculations that map airflow patterns within the nose with physiologically based pharmacokinetic (PBPK) models that integrate diffusion, metabolism, and tissue interactions of these compounds. Dose metrics estimated in these models for MMA and VA, respectively, were rates of MMA metabolism per volume of tissue and alterations in pH in target tissues associated with VA hydrolysis and metabolism. In this article, four scientists who have contributed significantly to development of these models describe the many similarities and relatively few differences between the MMA and VA models. Some differences arise naturally because of differences in target tissues, in the calculated measures of tissue dose, and in the modes of action for highly extracted vapors (VA) compared with poorly extracted vapors (MMA). A difference in the approach used to estimate metabolic parameters from human tissues provides insights into interindividual extrapolation and identifies opportunities for studies with human nasal tissues to enhance current risk assessments. In general, the differences in model structure for these two esters were essential for describing the biology of the observed responses and in accounting for the different measures of target tissue dose. This article is intended to serve as a guide for understanding issues of optimum model structure and optimal data sources for these nasal tissue dosimetry models. We also hope that it leads to greater international acceptance of these hybrid CFD/PBPK modeling approaches for improving risk assessment for many nasal toxicants. In general, these models predict either equivalent (VA) or lower (MMA) nasal tissue doses in humans compared with tissue doses at equivalent exposure concentrations in rats.

Dosimetry modeling of inhaled formaldehyde: comparisons of local flux predictions in the rat, monkey, and human nasal passages

Formaldehyde-induced nasal squamous cell carcinomas in rats and squamous metaplasia in rats and rhesus monkeys occur in specific regions of the nose with species-specific distribution patterns. Experimental approaches addressing local differences in formaldehyde uptake patterns and dose are limited by the resolution of dissection techniques used to obtain tissue samples and the rapid metabolism of absorbed formaldehyde in the nasal mucosa. Anatomically accurate, 3-dimensional computational fluid dynamics models of F344 rat, rhesus monkey, and human nasal passages were used to estimate and compare regional inhaled formaldehyde uptake patterns predicted among these species. Maximum flux values, averaged over a breath, in nonsquamous epithelium were estimated to be 2620, 4492, and 2082 pmol/(mm(2)-h-ppm) in the rat, monkey, and human respectively. Flux values predicted in sites where cell proliferation rates were measured as similar in rats and monkeys were also similar, as were fluxes predicted in a region of high tumor incidence in the rat nose and the anterior portion of the human nose. Regional formaldehyde flux estimates are directly applicable to clonal growth modeling of formaldehyde carcinogenesis to help reduce uncertainty in human cancer risk estimates.

Dosimetry modeling of inhaled formaldehyde: binning nasal flux predictions for quantitative risk assessment

Interspecies extrapolations of tissue dose and tumor response have been a significant source of uncertainty in formaldehyde cancer risk assessment. The ability to account for species-specific variation of dose within the nasal passages would reduce this uncertainty. Three-dimensional, anatomically realistic, computational fluid dynamics (CFD) models of nasal airflow and formaldehyde gas transport in the F344 rat, rhesus monkey, and human were used to predict local patterns of wall mass flux (pmol/[mm(2)-h-ppm]). The nasal surface of each species was partitioned by flux into smaller regions (flux bins), each characterized by surface area and an average flux value. Rat and monkey flux bins were predicted for steady-state inspiratory airflow rates corresponding to the estimated minute volume for each species. Human flux bins were predicted for steady-state inspiratory airflow at 7.4, 15, 18, 25.8, 31.8, and 37 l/min and were extrapolated to 46 and 50 l/min. Flux values higher than half the maximum flux value (flux median) were predicted for nearly 20% of human nasal surfaces at 15 l/min, whereas only 5% of rat and less than 1% of monkey nasal surfaces were associated with fluxes higher than flux medians at 0.576 l/min and 4.8 l/min, respectively. Human nasal flux patterns shifted distally and uptake percentage decreased as inspiratory flow rate increased. Flux binning captures anatomical effects on flux and is thereby a basis for describing the effects of anatomy and airflow on local tissue disposition and distributions of tissue response. Formaldehyde risk models that incorporate flux binning derived from anatomically realistic CFD models will have significantly reduced uncertainty compared with risk estimates based on default methods.

Mode of action and tissue dosimetry in current and future risk assessments

Two fundamental concepts have emerged to organize contemporary approaches to chemical risk assessment - mode of action and tissue dosimetry. Mode of action specifies the nature of the interactions between the chemical and the body that lead to toxic responses and should, under optimal circumstances, also specify the form of the tissue dose that leads to these effects. This paper highlights recent development of biologically based dose response (BBDR) models for specific toxic endpoints that use knowledge on mode of action to specify measures of dose. These dose measures then are used to support low dose and interspecies extrapolations. We first focus on a series of dose response models developed for several compounds that produce nasal toxicity. These examples demonstrate a range of model structures from simple dosimetry models (methylmethacrylate) to linkage of dosimetry with specific biological processes involved in carcinogenesis (formaldehyde). Two BBDR models with dioxin illustrate the organization of biological and dosimetry information into specific testable hypotheses that could distinguish these different models and lead to a more uniform approach to risk assessment for this compound. A final section discusses the impact of molecular biology and the genomic revolution in relation to development of BBDR models for specific toxic endpoints.

Using Human Data to Protect the Public's Health

The value of using human data in the assessment and management of risk is evaluated. Although the use of such data has a long and successful history with environmental contaminants and the development of drugs and commercial chemicals, recent deliberations within the Environmental Protection Agency (EPA) have questioned this practice in part. Specifically, we evaluate the degree to which reference doses (RfDs) and reference concentrations (RfCs) derived from human data on EPA's Integrated Risk Information System (IRIS) differ with RfDs and RfCs that we estimate from experimental animal data. We also use several minimal risk levels of the Agency for Toxic Substances and Disease Registry (ATSDR) and tolerable intakes of Health Canada in this comparison. Human-based RfDs are more than threefold lower than the corresponding animal-based RfDs for 23% of the comparisons. Human- based RfDs or RfCs are lower than corresponding animal-based RfDs or RfCs for 36% of the comparisons. Furthermore, for 10 of 43 possible comparisons, insufficient experimental animal data are readily available or data are inappropriate to estimate either RfDs or RfCs. We also discuss human pharmacokinetic data from volunteer studies and mechanistic studies with human tissues in vitro and demonstrate through a series of case discussions that utilization of such data is important when making decisions to protect exposed individuals. Moreover, physiologically based pharmacokinetic (PBPK) modeling evaluates critical information in assessing interindividual variability and identifying at-risk populations. Within the limits of our analysis, we conclude that the direct use and interpretation of human data, in conjunction with data gathered from experimental animals, are public health protective policies that should be encouraged.