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

Inhalation dosimetry of diacetyl and butyric acid, two components of butter flavoring vapors

Occupational exposure to butter flavoring vapors (BFV) is associated with significant pulmonary injury. The goal of the current study was to characterize inhalation dosimetric patterns of diacetyl and butyric acid, two components of BFV, and to develop a hybrid computational fluid dynamic-physiologically based pharmacokinetic model (CFD-PBPK) to describe these patterns. Uptake of diacetyl and butyric acid vapors, alone and in combination, was measured in the upper respiratory tract of anesthetized male Sprague-Dawley rats under constant velocity flow conditions and the uptake data were used to validate the CFD-PBPK model. Diacetyl vapor (100 or 300 ppm) was scrubbed from the airstream with 76-36% efficiency at flows of 100-400 ml/min. Butryic acid (30 ppm) was scrubbed with >90% efficiency. Concurrent exposure to butyric acid resulted in a small but significant reduction of diacetyl uptake (36 vs. 31%, p < 0.05). Diacetyl was metabolized in nasal tissues in vitro, likely by diacetyl reductase, an enzyme known to be inhibited by butyric acid. The CFD-PBPK model closely described diacetyl uptake; the reduction in diacetyl uptake by butyric acid could be explained by inhibition of diacetyl reductase. Extrapolation to the human via the model suggested that inspired diacetyl may penetrate to the intrapulmonary airways to a greater degree in the human than in the rat. Thus, based on dosimetric relationships, extrapulmonary airway injury in the rat may be predictive of intrapulmonary airway injury in humans. Butyric acid may modulate diacetyl toxicity by inhibiting its metabolism and/or altering its inhalation dosimetric patterns.

Approaches for applications of physiologically based pharmacokinetic models in risk assessment

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

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.

Framework for evaluation of physiologically-based pharmacokinetic models for use in safety or risk assessment

Proposed applications of increasingly sophisticated biologically-based computational models, such as physiologically-based pharmacokinetic models, raise the issue of how to evaluate whether the models are adequate for proposed uses, including safety or risk assessment. A six-step process for model evaluation is described. It relies on multidisciplinary expertise to address the biological, toxicological, mathematical, statistical, and risk assessment aspects of the modeling and its application. The first step is to have a clear definition of the purpose(s) of the model in the particular assessment; this provides critical perspectives on all subsequent steps. The second step is to evaluate the biological characterization described by the model structure based on the intended uses of the model and available information on the compound being modeled or related compounds. The next two steps review the mathematical equations used to describe the biology and their implementation in an appropriate computer program. At this point, the values selected for the model parameters (i.e., model calibration) must be evaluated. Thus, the fifth step is a combination of evaluating the model parameterization and calibration against data and evaluating the uncertainty in the model outputs. The final step is to evaluate specialized analyses that were done using the model, such as modeling of population distributions of parameters leading to population estimates for model outcomes or inclusion of early pharmacodynamic events. The process also helps to define the kinds of documentation that would be needed for a model to facilitate its evaluation and implementation.