In situ casting and imaging of the rat airway tree for accurate 3D reconstruction

Dynamic airway constriction in rats: heterogeneity and response to deep inspiration

Airway narrowing due to hyperresponsiveness severely limits gas exchange in patients with asthma. Imaging studies in humans and animals have shown that bronchoconstriction causes patchy patterns of ventilation defects throughout the lungs, and several computational models have predicted that these regions are due to constriction of smaller airways. However, these imaging approaches are often limited in their ability to capture dynamic changes in small airways, and the patterns of constriction are heterogeneous. To directly investigate regional variations in airway narrowing and the response to deep inspirations (DIs), we utilized tantalum dust and microfocal X-ray imaging of rat lungs to obtain dynamic images of airways in an intact animal model. Airway resistance was simultaneously measured using the flexiVent system. Custom-developed software was used to track changes in airway diameters up to generation 19 (~0.3-3 mm). Changes in diameter during bronchoconstriction were then measured in response to methacholine (MCh) challenge. In contrast with the model predictions, we observed significantly greater percent constriction in larger airways in response to MCh challenge. Although there was a dose-dependent increase in total respiratory resistance with MCh, the percent change in airway diameters was similar for increasing doses. A single DI following MCh caused a significant reduction in resistance but did not cause a significant increase in airway diameters. Multiple DIs did, however, cause significant increases in airway diameters. These measurements allowed us to directly quantify dynamic changes in airways during bronchoconstriction and demonstrated greater constriction in larger airways.

Computational fluid dynamics modeling of Bacillus anthracis spore deposition in rabbit and human respiratory airways

Three-dimensional computational fluid dynamics and Lagrangian particle deposition models were developed to compare the deposition of aerosolized Bacillus anthracis spores in the respiratory airways of a human with that of the rabbit, a species commonly used in the study of anthrax disease. The respiratory airway geometries for each species were derived respectively from computed tomography (CT) and μCT images. Both models encompassed airways that extended from the external nose to the lung with a total of 272 outlets in the human model and 2878 outlets in the rabbit model. All simulations of spore deposition were conducted under transient, inhalation-exhalation breathing conditions using average species-specific minute volumes. Two different exposure scenarios were modeled in the rabbit based upon experimental inhalation studies. For comparison, human simulations were conducted at the highest exposure concentration used during the rabbit experimental exposures. Results demonstrated that regional spore deposition patterns were sensitive to airway geometry and ventilation profiles. Due to the complex airway geometries in the rabbit nose, higher spore deposition efficiency was predicted in the nasal sinus compared to the human at the same air concentration of anthrax spores. In contrast, higher spore deposition was predicted in the lower conducting airways of the human compared to the rabbit lung due to differences in airway branching pattern. This information can be used to refine published and ongoing biokinetic models of inhalation anthrax spore exposures, which currently estimate deposited spore concentrations based solely upon exposure concentrations and inhaled doses that do not factor in species-specific anatomy and physiology for deposition.

Non-uniform central airways ventilation model based on vascular segmentation

Improvements in the understanding of the physiology of the central airways require an appropriate representation of the non-uniform ventilation at its terminal branches. This paper proposes a new technique for estimating the non-uniform ventilation at the terminal branches by modelling the volume change of their distal peripheral airways, based on vascular segmentation. The vascular tree is used for sectioning the dynamic CT-based 3D volume of the lung at 11 time points over the breathing cycle of a research animal. Based on the mechanical coupling between the vascular tree and the remaining lung tissues, the volume change of each individual lung segment over the breathing cycle was used to estimate the non-uniform ventilation of its associated terminal branch. The 3D lung sectioning technique was validated on an airway cast model of the same animal pruned to represent the truncated dynamic CT based airway geometry. The results showed that the 3D lung sectioning technique was able to estimate the volume of the missing peripheral airways within a tolerance of 2%. In addition, the time-varying non-uniform ventilation distribution predicted by the proposed sectioning technique was validated against CT measurements of lobar ventilation and showed good agreement. This significant modelling advance can be used to estimate subject-specific non-uniform boundary conditions to obtain subject-specific numerical models of the central airway flow.

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.

Respiratory tract lung geometry and dosimetry model for male Sprague-Dawley rats

While inhalation toxicological studies of various compounds have been conducted using a number of different strains of rats, mechanistic dosimetry models have only had tracheobronchial (TB) structural data for Long-Evans rats, detailed morphometric data on the alveolar region of Sprague-Dawley rats and limited alveolar data on other strains. Based upon CT imaging data for two male Sprague-Dawley rats, a 15-generation, symmetric typical path model was developed for the TB region. Literature data for the alveolar region of Sprague-Dawley rats were analyzed to develop an eight-generation model, and the two regions were joined to provide a complete lower respiratory tract model for Sprague-Dawley rats. The resulting lung model was used to examine particle deposition in Sprague-Dawley rats and to compare these results with predicted deposition in Long-Evans rats. Relationships of various physiologic variables and lung volumes were either developed in this study or extracted from the literature to provide the necessary input data for examining particle deposition. While the lengths, diameters and branching angles of the TB airways differed between the two Sprague-Dawley rats, the predicted deposition patterns in the three major respiratory tract regions were very similar. Between Sprague-Dawley and Long-Evans rats, significant differences in TB and alveolar predicted deposition fractions were observed over a wide range of particle sizes, with TB deposition fractions being up to 3- to 4-fold greater in Sprague-Dawley rats and alveolar deposition being significantly greater in Long-Evans rats. Thus, strain-specific lung geometry models should be used for particle deposition calculations and interspecies dose comparisons.

Casting materials and their application in research and teaching

From a biological point of view, casting refers to filling of anatomical and/or pathological spaces with extraneous material that reproduces a three-dimensional replica of the space. Casting may be accompanied by additional procedures such as corrosion, in which the soft tissue is digested out, leaving a clean cast, or the material may be mixed with radiopaque substances to allow x-ray photography or micro computed topography (µCT) scanning. Alternatively, clearing of the surrounding soft tissue increases transparency and allows visualization of the casted cavities. Combination of casting with tissue fixation allows anatomical dissection and didactic surgical procedures on the tissue. Casting materials fall into three categories namely, aqueous substances (India ink, Prussian blue ink), pliable materials (gelatins, latex, and silicone rubber), or hard materials (methyl methacrylates, polyurethanes, polyesters, and epoxy resins). Casting has proved invaluable in both teaching and research and many phenomenal biological processes have been discovered through casting. The choice of a particular material depends inter alia on the targeted use and the intended subsequent investigative procedures, such as dissection, microscopy, or µCT. The casting material needs to be pliable where anatomical and surgical manipulations are intended, and capillary-passable for ultrastructural investigations.