A computational fluid dynamics analysis of the effects of size and shape of anterior nasal septal perforations.

BACKGROUND: Nasal septal perforations (NSPs) often cause bleeding, crusting, obstruction, and/or whistling. The objective was to analyze the impact of anterior NSP size and shape on nasal physiology using computational fluid dynamics (CFD). METHODS: A 3-dimensional model of the nasal cavity was constructed from a radiologically normal CT scan using imaging software. Anterior NSPs (ovoid (ONSP): 0.5, 1, 2, and 3 cm long anterior-to-posteriorly and round (RNSP, 0.5 and 1 cm)) were virtually created in the model and divided into ventral, dorsal, anterior, and posterior regions. Steady-state inspiratory airflow, heat, and water vapor transport were simulated using Fluent CFD software. Air crossover through the perforation, wall shear, heat flux, water vapor flux, resistance, and humidification were analyzed. RESULTS: Air crossover and wall shear increased with perforation size. Regionally, wall shear and heat and water vapor flux were highest posteriorly and lowest anteriorly, generally increasing with size in those regions. RNSPs had greater heat and water vapor flux compared to corresponding size ONSPs. Resistance decreased by 10% or more from normal only in the 3 cm ONSP. Maximum water content was achieved more posteriorly in larger NSP nasal cavities. CONCLUSIONS: High wall shear and heat and water vapor flux in posterior perforation regions may explain the crusting most commonly noted on posterior NSP edges. This preliminary study suggests that larger NSPs have a greater effect on nasal resistance and water content. Decrease in resistance with larger NSP size may be implicated in reported symptomatic improvement following enlargement of NSPs for treatment.

A computational study of functional endoscopic sinus surgery and maxillary sinus drug delivery.

BACKGROUND: Topical medication is increasingly used following functional endoscopic sinus surgery (FESS). Information on particle sizes that maximise maxillary sinus (MS) delivery is conflicting, and the effect of antrostomy size on delivery is unclear. The purpose of this study was to estimate antrostomy and particle size effects on topical MS drug delivery. METHODOLOGY: Sinonasal reconstructions were created from a pre- and a post-FESS CT scan in each of four chronic rhinosinusitis patients. Additional models were created from each post-FESS reconstruction representing four alternative antrostomy sizes. Airflow and particle deposition were simulated in each reconstruction using computational fluid dynamics for nebulised and sprayed delivery. RESULTS: MS ventilation and drug delivery increased following FESS, the largest virtual antrostomy led to greatest delivery, and MS delivery was sensitive to particle size. Particles within a 5-18 µm and 5-20 µm size range led to peak MS deposition for nebulised and sprayed particles, respectively. Post-FESS increases in drug delivery varied across individuals and within individuals by the type of antrostomy created. CONCLUSION: Our findings suggest that FESS, particularly with larger antrostomies, improves topical drug delivery, and that certain particle sizes improve this delivery. Further research is needed to contextualise these findings with other post-surgical effects.

A lung dosimetry model of vapor uptake and tissue disposition.

Inhaled vapors may be absorbed at the alveolar-capillary membrane and enter arterial blood flow to be carried to other organs of the body. Thus, the biological effects of inhaled vapors depend on vapor uptake in the lung and distribution to the rest of the body. A mechanistic model of vapor uptake in the human lung and surrounding tissues was developed for soluble and reactive vapors during a single breath. Lung uptake and tissue disposition of inhaled formaldehyde, acrolein, and acetaldehyde were simulated for different solubilities and reactivities. Formaldehyde, a highly reactive and soluble vapor, was estimated to be taken up by the tissues in the upper tracheobronchial airways with shallow penetration into the lung. Vapors with moderate solubility such as acrolein and acetaldehyde were estimated to penetrate deeper into the lung, reaching the alveolar region where absorbed vapors had a much higher probability of passing through the thin alveolar-capillary membrane to reach the blood. For all vapors, tissue concentration reached its maximum at the end of inhalation at the air-tissue interface. The depth of peak concentration moved within the tissue layer due to vapor desorption during exhalation. The proposed vapor uptake model offers a mechanistic approach for calculations of lung vapor uptake, air:tissue flux, and tissue concentration profiles within the respiratory tract that can be correlated to local biological response in the lung. In addition, the uptake model provides the necessary input for pharmacokinetic models of inhaled chemicals in the body, thus reducing the need for estimating requisite parameters.

Deviated nasal septum hinders intranasal sprays: a computer simulation study.

BACKGROUND: This study investigates how deviated nasal septum affects the quantity and distribution of spray particles, and examines the effects of inspiratory airflow and head position on particle transport. METHODS: Deposition of spray particles was analysed using a three-dimensional computational fluid dynamics model created from a computed tomography scan of a human nose with leftward septal deviation and a right inferior turbinate hypertrophy. Five simulations were conducted using FluentTM software, with particle sizes ranging from 20-110 μm, a spray speed of 3 m/s, plume angle of 68(deg), and with steady state inspiratory airflow either present (15.7 L/min) or absent at varying head positions. RESULTS: With inspiratory airflow present, posterior deposition on the obstructed side was approximately four times less than the contralateral side, regardless of head position, and was statistically significant. When airflow was absent, predicted deposition beyond the nasal valve on the left and right sides were between 16% and 69% lower and positively influenced by a dependent head position. CONCLUSION: Simulations predicted that septal deviation significantly diminished drug delivery on the obstructed side. Furthermore, increased particle penetration was associated with presence of nasal airflow. Head position is an important factor in particle deposition patterns when inspiratory airflow is absent.

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.

Correlation of regional and nonlinear formaldehyde-induced nasal cancer with proliferating populations of cells.

Formaldehyde induces nonlinear, concentration-related increases in nasal epithelial cell proliferation and squamous cell carcinomas (SCC) in rats. A formaldehyde carcinogenicity study was conducted in which a major end point was correlation of cell proliferation indices with sites of formaldehyde-induced SCC. A poor correlation in certain sites led to incorporation of the number of cells in each site into the correlation. Rats were exposed (6h/day, 5 days/week) to formaldehyde (0, 0.7, 2, 6, 10 or 15 ppm) for up to 24 months with interim sacrifice time points at 3, 6, 12, and 18 mo. A unit length labeling index (ULLI; S-phase nuclei/mm basement membrane) was determined for specific nasal regions in addition to a population-weighted ULLI (PWULLI). The PWULLI was defined as the product of regional ULLI and total number of nasal epithelial cells in the respective site. Nasal SCC sites of origin were mapped. Formaldehyde induced SCC in a highly nonlinear fashion, with no observed effect at the level of 2 ppm, a minimal response at 6 ppm, and a sharp increase at 10 and 15 ppm. The tumor incidence was 1, 22, and 47% at 6, 10 and 15 ppm, respectively. ULLI was significantly (P<0.05) increased at 10 and 15 ppm but not at the lower concentrations. There was a good correlation between PWULLI and regional tumor incidence (R(2) = 0.88), while the correlation of regional SCC with ULLI was relatively poor (R(2) = 0.46). We conclude that target cell population size and sustained increases of cell proliferation in these populations, determined by differences in regional airflow-driven formaldehyde binding to DNA dose to these sites, coupled with the known nonlinear kinetics of formaldehyde binding to DNA, can together account for the nonlinearity and site specificity of formaldehyde-induced nasal SCC in rats.

Simulation modeling of the tissue disposition of formaldehyde to predict nasal DNA-protein cross-links in Fischer 344 rats, rhesus monkeys, and humans.

Formaldehyde inhalation causes formation of DNA-protein cross-links (DPX) in the nasal mucosa of Fischer 344 (F344) rats and rhesus monkeys. DPX are considered to be part of the mechanism by which cytotoxic and carcinogenic effects of formaldehyde in laboratory animals are exerted, and DPX data have been used as a measure of tissue dose in cancer risk assessments for formaldehyde. Accurate prediction of DPX concentrations in humans is therefore desirable. The goal of this work was to increase confidence in the prediction of human DPX by refining earlier models of formaldehyde disposition and DPX kinetics in the nasal mucosa. Anatomically accurate, computational fluid dynamics models of the nasal airways of F344 rats, rhesus monkeys, and humans were used to predict the regional flux of formaldehyde to the respiratory and olfactory mucosa. A previously developed model of the tissue disposition of formaldehyde and of DPX kinetics was implemented in the graphical simulation tool SIMULINK and linked to the regional flux predictions. Statistical optimization was used to identify parameter values, and good simulations of the data were obtained. The parameter estimates for rats and monkeys were used to guide allometric scale-up to the human case. The relative levels of nasal mucosal DPX in rats, rhesus monkeys, and humans for a given inhaled concentration of formaldehyde were predicted by the model to vary with concentration. This modeling approach reduces uncertainty in the prediction of human nasal mucosal DPX resulting from formaldehyde inhalation.

Mass-transport models to predict toxicity of inhaled gases in the upper respiratory tract.

Mass transport (the movement of a chemical species) plays an important role in determining toxic responses of the upper respiratory tract (URT) to inhaled chemicals. Mathematical dosimetry models incorporate physical characteristics of mass transport and are used to predict quantitative uptake (absorption rate) and distribution of inhaled gases and vapors in the respiratory tract. Because knowledge of dose is an essential component of quantitative risk assessment, dosimetry modeling plays an important role in extrapolation of animal study results to humans. A survey of existing mathematical dosimetry models for the URT is presented, limitations of current models are discussed, and adaptations of existing models to produce a generally applicable model are suggested. Reviewed URT dosimetry models are categorized as early, lumped-parameter, and distributed- parameter models. Specific examples of other relevant modeling work are also presented.

Dosimetry modeling of inhaled formaldehyde: the human respiratory tract.

Formaldehyde (HCHO), which has been shown to be a nasal carcinogen in rats and mice, is used widely and extensively in various manufacturing processes. Studies in rhesus monkeys suggest that the lower respiratory tract may be at risk and some epidemiologic studies have reported an increase in lung cancer associated with HCHO; other studies have not. Thus, an assessment of possible human risk to HCHO exposure based on dosimetry information throughout the respiratory tract (RT) is desirable. To obtain dosimetry estimates for a risk assessment, two types of models were used. The first model (which is the subject of another investigation) used computational fluid dynamics (CFD) to estimate local fluxes in a 3-dimensional model of the nasal region. The subject of the present investigation (the second model) applied a 1-dimensional equation of mass transport to each generation of an adult human symmetric, bifurcating Weibel-type RT anatomical model, augmented by an upper respiratory tract. The two types of modeling approaches were made consistent by requiring that the 1-dimensional version of the nasal passages have the same inspiratory air-flow rate and uptake during inspiration as the CFD simulations for 4 daily human activity levels. Results obtained include the following: (1) More than 95% of the inhaled HCHO is predicted to be retained by the RT. (2) The CFD predictions for inspiration, modified to account for the difference in inspiration and complete breath times, are a good approximation to uptake in the nasal airways during a single breath. (3) In the lower respiratory tract, flux is predicted to increase for several generations and then decrease rapidly. (4) Compared to first pulmonary region generation fluxes, the first few tracheobronchial generations fluxes are over 1000 times larger. Further, there is essentially no flux in the alveolar sacs. (5) Predicted fluxes based on the 1-dimensional model are presented that can be used in a biologically based dose-response model for human carcinogenesis. Use of these fluxes will reduce uncertainty in a risk assessment for formaldehyde carcinogenicity.

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.

Analysis of vinyl acetate metabolism in rat and human nasal tissues by an in vitro gas uptake technique.

Physiologically based pharmacokinetic (PBPK) models require estimates of catalytic rate constants controlling the metabolism of xenobiotics. Usually, these constants are derived from whole tissue homogenates wherein cellular architecture and enzyme compartmentation are destroyed. Since the nasal cavity epithelium is composed of a heterogeneous cell population measurement of xenobiotic metabolizing enzymes using homogenates could yield artifactual results. In this article a method for measuring rates of metabolism of vinyl acetate, a metabolism-dependent carcinogen, is presented that uses whole-tissue samples and PBPK modeling techniques to estimate metabolic kinetic parameters in tissue compartments. The kinetic parameter estimates were compared to those derived from homogenate experiments using two methods of tissue normalization. When the in vitro gas uptake constants were compared to homogenate-derived values, using a normalization procedure that does not account for tissue architecture, there was poor agreement. Homogenate-derived values from rat nasal tissue were 3- to 23-fold higher than those derived using the in vitro gas uptake method. When the normalization procedure for the rat homogenate-derived values took into account tissue architecture, a good agreement was observed. Carboxylesterase activity in homogenates of human nasal tissues was undetectable. Using the in vitro gas uptake technique, however, carboxylesterase activity was detected. Rat respiratory carboxylesterase and aldehyde dehydrogenase activities were about three and two times higher than those of humans, respectively. Activities of the rat olfactory enzymes were about equivalent to those of humans. K(m) values did not differ between species. The results suggest that the in vitro gas uptake technique is useful for deriving enzyme kinetic constants where effects of tissue architecture are preserved. Furthermore, the results suggest that caution should be exercised when scaling homogenate-derived values to whole-organ estimates, especially in organs of cellular heterogeneity.

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.

Correlation of regional formaldehyde flux predictions with the distribution of formaldehyde-induced squamous metaplasia in F344 rat nasal passages.

Squamous epithelium lines the nasal vestibule of the rat, rhesus monkey, and human. Respiratory, transitional, and olfactory epithelia line most areas posterior to the nasal vestibule. Inhaled formaldehyde gas induces squamous metaplasia posterior to the nasal vestibule and does not induce lesions in the nasal vestibule in rats and rhesus monkeys, indicating that squamous epithelium is resistant to irritant effects of formaldehyde and that squamous metaplasia may be an adaptive response. If squamous metaplasia is determined by formaldehyde dosimetry rather than by tissue-specific factors, squamous epithelium may be protective by absorbing less formaldehyde than other epithelial types. In a previous study, a three-dimensional, anatomically accurate computational fluid dynamics (CFD) model of the anterior F344 rat nasal passages was used to simulate inspiratory airflow and inhaled formaldehyde transport. The present study consisted of two related parts. First, the rat CFD model was used to test the hypothesis that the distribution of formaldehyde-induced squamous metaplasia is related to the location of high-flux regions posterior to squamous epithelium. Regional formaldehyde flux into nonsquamous epithelium predicted by the CFD model correlated with regional incidence of formaldehyde-induced squamous metaplasia on the airway perimeter of one cross-sectional level of the noses of F344 rats exposed to 10 and 15 ppm formaldehyde gas for 6 months. Formaldehyde flux into nonsquamous epithelium was estimated to vary by an order of magnitude depending on the degree of formaldehyde absorption by squamous epithelium. These results indicate that the degree to which squamous epithelium absorbs formaldehyde strongly affects the rate and extent of the progression of squamous metaplasia with continued exposure to formaldehyde. In the second part of this study, the CFD model was used to predict squamous metaplasia progression. Data needs for verification of this model prediction are considered. These results indicate that information on the permeability of squamous epithelium in rats, monkeys, and humans is important for accurate prediction of uptake in regions posterior to the nasal vestibule.

Reconstruction of complex passageways for simulations of transport phenomena: development of a graphical user interface for biological applications.

Flow of fluids, such as blood, lymph and air, plays a major role in the normal physiology of all living organisms. Within individual organ systems, flow fields may significantly influence the transport of solutes, including nutrients and chemical toxicants, to and from the confining vessel walls (epithelia and endothelia). Computational fluid dynamics (CFD) provides a potentially useful tool for biologists and toxicologists investigating solute disposition in these flow fields in both normal and disease states. Application of CFD is dependent upon generation of accurate representations of the geometry of the system of interest in the form of a computational reconstruction. The present investigations, which were based on studies of the toxicology of inhaled reactive gases in the respiratory tract of rodents, provide computer programs for the generation of finite element meshes from serial tissue cross-sections. These programs, which interface with a commercial finite element fluid dynamics simulation package (FIDAP 7.05, Fluid Dynamics International, Evanston, IL), permit simulation of fluid flow in the complex geometries and local solute mass flux to the vessel walls of biological systems. The use of these programs and their application to studies of respiratory tract toxicology are described.

Changes in nasal airflow and heat transfer correlate with symptom improvement after surgery for nasal obstruction.

Surgeries to correct nasal airway obstruction (NAO) often have less than desirable outcomes, partly due to the absence of an objective tool to select the most appropriate surgical approach for each patient. Computational fluid dynamics (CFD) models can be used to investigate nasal airflow, but variables need to be identified that can detect surgical changes and correlate with patient symptoms. CFD models were constructed from pre- and post-surgery computed tomography scans for 10 NAO patients showing no evidence of nasal cycling. Steady-state inspiratory airflow, nasal resistance, wall shear stress, and heat flux were computed for the main nasal cavity from nostrils to posterior nasal septum both bilaterally and unilaterally. Paired t-tests indicated that all CFD variables were significantly changed by surgery when calculated on the most obstructed side, and that airflow, nasal resistance, and heat flux were significantly changed bilaterally as well. Moderate linear correlations with patient-reported symptoms were found for airflow, heat flux, unilateral allocation of airflow, and unilateral nasal resistance as a fraction of bilateral nasal resistance when calculated on the most obstructed nasal side, suggesting that these variables may be useful for evaluating the efficacy of nasal surgery objectively. Similarity in the strengths of these correlations suggests that patient-reported symptoms may represent a constellation of effects and that these variables should be tracked concurrently during future virtual surgery planning.

Derivation of mass transfer coefficients for transient uptake and tissue disposition of soluble and reactive vapors in lung airways.

Evaluation of vapor uptake by lung airways and subsequent dose to lung tissues provides the bridge connecting exposure episode to biological response. Respiratory vapor absorption depends on chemical properties of the inhaled material, including solubility, diffusivity, and metabolism/reactivity in lung tissues. Inter-dependent losses in the air and tissue phases require simultaneous calculation of vapor concentration in both phases. Previous models of lung vapor uptake assumed steady state, one-way transport into tissues with first-order clearance. A new approach to calculating lung dosimetry is proposed in which an overall mass transfer coefficient for vapor transport across the air-tissue interface is derived using air-phase mass transfer coefficients and analytical expressions for tissue-phase mass transfer coefficients describing unsteady transport by diffusion, first-order, and saturable pathways. Feasibility of the use of mass transfer coefficients was shown by calculating transient concentration levels of inhaled formaldehyde in the human tracheal airway and surrounding tissue. Formaldehyde tracheal air concentration and wall-flux declined throughout the breathing cycle. After the inhalation period, peak tissue concentration moved from the air-tissue interface into the tissue due to desorption into the air and continued diffusional transport across the tissue layer. While model predictions were performed for formaldehyde, which serves as a model of physiologically relevant, highly reactive vapors, the model is equally applicable to other soluble and reactive compounds.

Nasal dosimetry of inhaled gases and particles: where do inhaled agents go in the nose?

The anatomical structure of the nasal passages differs significantly among species, affecting airflow and the transport of inhaled gases and particles throughout the respiratory tract. Since direct measurement of local nasal dose is often difficult, 3-dimensional, anatomically accurate, computational models of the rat, monkey, and human nasal passages were developed to estimate regional transport and dosimetry of inhaled material. The computational models predicted that during resting breathing, a larger portion of inspired air passed through olfactory-lined regions in the rat than in the monkey or human. The models also predicted that maximum wall mass flux (mass per surface area per time) of inhaled formaldehyde in the nonsquamous epithelium was highest in monkeys (anterior middle turbinate) and similar in rats and humans (dorsal medial meatus in the rat and mid-septum in the human, near the squamous/nonsquamous epithelial boundary in both species). For particles that are 5 microm in aerodynamic diameter, preliminary simulations at minute volume flow rates predicted nasal deposition efficiencies of 92%, 11% and 25% in the rat, monkey, and human, respectively, with more vestibular deposition in the rat than in the monkey or human. Estimates such as these can be used to test hypotheses about mechanisms of toxicity and supply species-specific information for risk assessment, thus reducing uncertainty in extrapolating animal data to humans.

Computer simulation of cell growth governed by stochastic processes: application to clonal growth cancer models.

Cancer is a multistage process in which cell proliferation determines the growth of cells within stages and is associated with the transition of cells from one stage to the next. The usual model for cancer risk assessment, the linearized multistage model, does not explicitly include cell proliferation. More realistic cancer models are needed to reduce uncertainty in cancer risk assessment and to provide basic insights into the quantitative roles of cell proliferation and mutation. This report describes a simulation model for the transition of cells from one stage to the next and for clonal growth within stages. The model is intended to facilitate the use of experimental data on cell replication and preneoplastic lesions in risk assessment. When a population of cells is small its growth may be governed by stochastic processes. Such a population may disappear by chance even when the probability of cell division on a given time interval exceeds the probability of cell death. Procedures for estimating cell proliferation and mutation parameters from data for use in risk assessment should account for this random aspect of growth. The present model describes cell growth governed by stochastic processes, is consistent with earlier analytical expressions for such growth (Dewanjii et al., Risk Anal. 9, 179, 1989), and is flexible with respect to time-dependent data. A data set for spontaneous basophilic clones in male F344 rats (Popp et al., Fundam. Appl. Toxicol. 5, 314, 1985) is analyzed and predictions are made for (a) the probability of mutation to the basophilic genotype per division of a normal hepatocyte (3.5 x 10(-8)), (b) number of basophilic clones too small to be detected, and (c) number of basophilic clones that disappear by chance. This work illustrates the potential of computer simulation for quantitative analysis of the roles of cell division, cell death, and mutation in cancer.

Use of computational fluid dynamics models for dosimetry of inhaled gases in the nasal passages.

Computational fluid dynamics (CFD) models of the nasal passages of a rat, monkey, and human are being used (1) to determine important factors affecting nasal uptake, (2) to make interspecies dosimetric comparisons, (3) to provide detailed anatomical information for the rat, monkey, and human nasal passages, and (4) to provide estimates of regional air-phase mass transport coefficients (a measure of the resistance to gas transport from inhaled air to airway walls) in the nasal passages of all three species. For many inhaled materials, lesion location in the nose follows patterns that are both site and species specific. For reactive, water-soluble (Category 1) gases, regional uptake can be a major factor in determining lesion location. Since direct measurement of airflow and uptake is experimentally difficult, CFD models are used here to predict uptake patterns quantitatively in three-dimensional reconstructions of the F344 rat, rhesus monkey, and human nasal passages. In formaldehyde uptake simulations, absorption processes were assumed to be as rapid as possible, and regional flux (transport rate) of inhaled formaldehyde to airway walls was calculated for rats, primates, and humans. For uptake of gases like vinyl acetate and acrylic acid vapors, physiologically based pharmacokinetic uptake models incorporating anatomical and physical information from the CFD models were developed to estimate nasal tissue dose in animals and humans. The use of biologically based models in risk assessment makes sources of uncertainty explicit and, in doing so, allows quantification of uncertainty through sensitivity analyses. Limited resources can then be focused on reduction of important sources of uncertainty to make risk estimates more accurate.

Deposition of fine and coarse aerosols in a rat nasal mold.

The objective of this study was to investigate the deposition characteristics of large, inhalable particles in rat nasal passages by determining the deposition efficiencies of these particles in a nasal mold of an F344 rat for steady-state and pulsating flow conditions. Particles with geometric diameters ranging from 0.5 to 4 microm and flow rates ranging from 100 to 900 ml/min were employed for simulated inspiratory and expiratory flow situations. The optically clear acrylic mold was fabricated from a life-size metal cast that comprised the nares, nasal cavity, pharynx, and larynx. Deposition efficiencies were calculated for each flow situation and plotted as functions of particle inertia. Inspiratory and expiratory deposition efficiencies were similar for a given flow condition. Deposition efficiencies for the cases of pulsating flows were markedly higher than those of steady flows. The results for pulsating flows indicated higher deposition efficiencies than were found in previous studies performed with live rats. These differences may be due to uncertainties in particle inhalability, clearance, and flow rate in the previous studies, as well as differences between the nasal geometries of live rats and the geometry of the nasal mold made from a postmortem cast. The results suggest that the pulsating nature of breathing is an important consideration when determining the deposition of fine and coarse particles.