Richard E. Jacob

Comparative computational modeling of airflows and vapor dosimetry in the respiratory tracts of rat, monkey, and human

Computational fluid dynamics (CFD) models are useful for predicting site-specific dosimetry of airborne materials in the respiratory tract and elucidating the importance of species differences in anatomy, physiology, and breathing patterns. We improved the imaging and model development methods to the point where CFD models for the rat, monkey, and human now encompass airways from the nose or mouth to the lung. A total of 1272, 2172, and 135 pulmonary airways representing 17±7, 19±9, or 9±2 airway generations were included in the rat, monkey and human models, respectively. A CFD/physiologically based pharmacokinetic model previously developed for acrolein was adapted for these anatomically correct extended airway models. Model parameters were obtained from the literature or measured directly. Airflow and acrolein uptake patterns were determined under steady-state inhalation conditions to provide direct comparisons with prior data and nasal-only simulations. Results confirmed that regional uptake was sensitive to airway geometry, airflow rates, acrolein concentrations, air:tissue partition coefficients, tissue thickness, and the maximum rate of metabolism. Nasal extraction efficiencies were predicted to be greatest in the rat, followed by the monkey, and then the human. For both nasal and oral breathing modes in humans, higher uptake rates were predicted for lower tracheobronchial tissues than either the rat or monkey. These extended airway models provide a unique foundation for comparing material transport and site-specific tissue uptake across a significantly greater range of conducting airways in the rat, monkey, and human than prior CFD models.

19F MR imaging of ventilation and diffusion in excised lungs.

Perfluorinated gases, particularly C2F6, are potentially suitable alternatives to hyperpolarized noble gases for pulmonary airspace spin density and diffusion MRI. This work focuses mainly on 19F imaging of C2F6 gas in healthy and emphysematous explanted lungs, avoiding regulatory issues of human in vivo measurements. Three‐dimensional gradient echo and spin echo spin density images of human lungs can be made in 10 s with adequate signal‐to‐noise, demonstrating the feasibility for breathing dynamics to be captured during a succession of short breath holds. As expected, the spin echo images have much smaller susceptibility artifacts than the gradient echo images. 19F and 3He images of the same lungs are compared. The apparent diffusion coefficient (ADC) of C2F6 is sensitive to restrictions imposed by the lung microstructure: the average ADC is measured to be 0.018 cm2/s in healthy lungs versus 0.031 cm2/s in emphysematous lungs at a diffusion time Δ = 2.2 ms. The low free diffusivity of pure C2F6 (D0 = 0.033 cm2/s) places it in a regime where the ADC measurement allows the surface‐to‐volume ratio to be determined in each voxel, a potentially valuable quantitative characterization of regional lung tissue destruction in emphysema. Magn Reson Med, 2005.

Effects of diffusion time on short‐range hyperpolarized 3He diffusivity measurements in emphysema

To characterize the effect of diffusion time on short‐range hyperpolarized 3He magnetic resonance imaging (MRI) diffusion measurements across a wide range of emphysema severity.

Application of Magnetic Resonance (MR) Imaging for the Development and Validation of Computational Fluid Dynamic (CFD) Models of the Rat Respiratory System

Computational fluid dynamic (CFD) models of the respiratory system provide a quantitative basis for extrapolating the localized dose of inhaled materials and improving human health risk assessments based upon inhalation studies conducted in animals. Nevertheless, model development and validation have historically been tedious and time-consuming tasks. In recognition of this, we previously reported on the use of proton (1H) magnetic resonance (MR) imaging for visualizing nasal-sinus passages in the rat, and for speeding computational mesh generation. Here, the generation and refinement of meshes for rat nasal airways are described in more detail and simulated airflows are presented. To extend the CFD models to the complete respiratory tract, three-dimensional (3D) 1H MR imaging of rat pulmonary casts was also utilized to construct pulmonary airway meshes using procedures developed for the nasal airways. Furthermore, the feasibility of validating CFD predictions with MR was tested by imaging hyperpolarized 3He gas at physiological flow rates in a straight pipe with a diameter comparable to the rat trachea. Results from these diverse studies highlight the potential utility of MR imaging not only for speeding CFD development but also possibly for model validation.

In vivo MRI of altered proton signal intensity and T2 relaxation in a bleomycin model of pulmonary inflammation and fibrosis

To investigate the ability of proton (1H) magnetic resonance imaging (MRI) to distinguish between pulmonary inflammation and fibrosis.

Comparison of two quantitative methods of discerning airspace enlargement in smoke-exposed mice.

In this work, we compare two methods for evaluating and quantifying pulmonary airspace enlargement in a mouse model of chronic cigarette smoke exposure. Standard stereological sample preparation, sectioning, and imaging of mouse lung tissues were performed for semi-automated acquisition of mean linear intercept (Lm) data. After completion of the Lm measurements, D2, a metric of airspace enlargement, was measured in a blinded manner on the same lung images using a fully automated technique developed in-house. An analysis of variance (ANOVA) shows that although Lm was able to separate the smoke-exposed and control groups with statistical significance (p = 0.034), D2 was better able to differentiate the groups (p<0.001) and did so without any overlap between the control and smoke-exposed individual animal data. In addition, the fully automated implementation of D2 represented a time savings of at least 24x over semi-automated Lm measurements. Although D2 does not provide 3D stereological metrics of airspace dimensions as Lm does, results show that it has higher sensitivity and specificity for detecting the subtle airspace enlargement one would expect to find in mild or early stage emphysema. Therefore, D2 may serve as a more accurate screening measure for detecting early lung disease than Lm.

Phase-Contrast MRI and CFD Modeling of Apparent 3He Gas Flow in Rat Pulmonary Airways

Phase-contrast (PC) magnetic resonance imaging (MRI) with hyperpolarized ³He is potentially useful for developing and testing patient-specific models of pulmonary airflow. One challenge, however, is that PC-MRI provides apparent values of local ³He velocity that not only depend on actual airflow but also on gas diffusion. This not only blurs laminar flow patterns in narrow airways but also introduces anomalous airflow structure that reflects gas-wall interactions. Here, both effects are predicted in a live rat using computational fluid dynamics (CFD), and for the first time, simulated patterns of apparent ³He gas velocity are compared with in vivo PC-MRI. Results show (1) that correlations (R²) between measured and simulated airflow patterns increase from 0.23 to 0.79 simply by accounting for apparent ³He transport, and (2) that remaining differences are mainly due to uncertain airway segmentation and partial volume effects stemming from relatively coarse MRI resolution. Higher-fidelity testing of pulmonary airflow predictions should therefore be possible with future imaging improvements.

Development of a rhesus monkey lung geometry model and application to particle deposition in comparison to humans

The exposure-dose-response characterization of an inhalation hazard established in an animal species needs to be translated to an equivalent characterization in humans relative to comparable doses or exposure scenarios. Here, the first geometry model of the conducting airways for rhesus monkeys is developed based upon CT images of the conducting airways of a 6-month-old male, rhesus monkey. An algorithm was developed for adding the alveolar region airways using published rhesus morphometric data. The resultant lung geometry model can be used in mechanistic particle or gaseous dosimetry models. Such dosimetry models require estimates of the upper respiratory tract volume of the animal and the functional residual capacity, as well as of the tidal volume and breathing frequency of the animal. The relationship of these variables to rhesus monkeys of differing body weights was established by synthesizing and modeling published data as well as modeling pulmonary function measurements on 121 rhesus control animals. Deposition patterns of particles up to 10 µm in size were examined for endotracheal and and up to 5 µm for spontaneous breathing in infant and young adult monkeys and compared to those for humans. Deposition fraction of respirable size particles was found to be higher in the conducting airways of infant and young adult rhesus monkeys compared to humans. Due to the filtering effect of the conducting airways, pulmonary deposition in rhesus monkeys was lower than that in humans. Future research areas are identified that would either allow replacing assumptions or improving the newly developed lung model.

An automated self-similarity analysis of the pulmonary tree of the Sprague-Dawley rat.

We present the results of an automated analysis of the morphometry of the pulmonary airway trees of the Sprague–Dawley rat. Our work is motivated by a need to inform lower‐dimensional mathematical models to prescribe realistic boundary conditions for multiscale hybrid models of rat lung mechanics. Silicone casts were made from three age‐matched, male Sprague–Dawley rats, immersed in a gel containing a contrast agent and subsequently imaged with magnetic resonance (MR). From a segmentation of this data, we extracted a connected graph, representing the airway centerline. Segment statistics (lengths and diameters) were derived from this graph. To validate this MR imaging/digital analysis method, airway segment measurements were compared with nearly 1,000 measurements collected by hand using an optical microscope from one of the rat lung casts. To evaluate the reproducibility of the MR imaging/digital analysis method, two lung casts were each imaged three times with randomized orientations in the MR bore. Diameters and lengths of randomly selected airways were compared among each of the repeated imaging datasets to estimate the variability. Finally, we analyzed the morphometry of the airway tree by assembling individual airway segments into structures that span multiple generations, which we call branches. We show that branches not segments are the fundamental repeating unit in the rat lung and develop simple mathematical relationships describing these structures for the entire lung. Our analysis shows that airway diameters and lengths have both a deterministic and stochastic character. Anat Rec, 2008.

Comparative Risks of Aldehyde Constituents in Cigarette Smoke Using Transient Computational Fluid Dynamics/Physiologically Based Pharmacokinetic Models of the Rat and Human Respiratory Tracts

Computational fluid dynamics (CFD) modeling is well suited for addressing species-specific anatomy and physiology in calculating respiratory tissue exposures to inhaled materials. In this study, we overcame prior CFD model limitations to demonstrate the importance of realistic, transient breathing patterns for predicting site-specific tissue dose. Specifically, extended airway CFD models of the rat and human were coupled with airway region-specific physiologically based pharmacokinetic (PBPK) tissue models to describe the kinetics of 3 reactive constituents of cigarette smoke: acrolein, acetaldehyde and formaldehyde. Simulations of aldehyde no-observed-adverse-effect levels for nasal toxicity in the rat were conducted until breath-by-breath tissue concentration profiles reached steady state. Human oral breathing simulations were conducted using representative aldehyde yields from cigarette smoke, measured puff ventilation profiles and numbers of cigarettes smoked per day. As with prior steady-state CFD/PBPK simulations, the anterior respiratory nasal epithelial tissues received the greatest initial uptake rates for each aldehyde in the rat. However, integrated time- and tissue depth-dependent area under the curve (AUC) concentrations were typically greater in the anterior dorsal olfactory epithelium using the more realistic transient breathing profiles. For human simulations, oral and laryngeal tissues received the highest local tissue dose with greater penetration to pulmonary tissues than predicted in the rat. Based upon lifetime average daily dose comparisons of tissue hot-spot AUCs (top 2.5% of surface area-normalized AUCs in each region) and numbers of cigarettes smoked/day, the order of concern for human exposures was acrolein > formaldehyde > acetaldehyde even though acetaldehyde yields were 10-fold greater than formaldehyde and acrolein.