Study of the variability in upper and lower airway morphology in Sprague-Dawley rats using modern micro-CT scan-based segmentation techniques

Computational, In Vitro, and In Vivo Models for Nose-to-Brain Drug Delivery Studies

Direct nose-to-brain drug delivery offers the opportunity to treat central nervous system disorders more effectively due to the possibility of drug molecules reaching the brain without passing through the blood-brain barrier. Such a delivery route allows the desired anatomic site to be reached while ensuring drug effectiveness, minimizing side effects, and limiting drug losses and degradation. However, the absorption of intranasally administered entities is a complex process that considerably depends on the interplay between the characteristics of the drug delivery systems and the nasal mucosa. Various preclinical models (in silico, in vitro, ex vivo, and in vivo) are used to study the transport of drugs after intranasal administration. The present review article attempts to summarize the different computational and experimental models used so far to investigate the direct delivery of therapeutic agents or colloidal carriers from the nasal cavity to the brain tissue. Moreover, it provides a critical evaluation of the data available from different studies and identifies the advantages and disadvantages of each model.

Effect of MDI Actuation Timing on Inhalation Dosimetry in a Human Respiratory Tract Model

Accurate knowledge of the delivery of locally acting drug products, such as metered-dose inhaler (MDI) formulations, to large and small airways is essential to develop reliable in vitro/in vivo correlations (IVIVCs). However, challenges exist in modeling MDI delivery, due to the highly transient multiscale spray formation, the large variability in actuation–inhalation coordination, and the complex lung networks. The objective of this study was to develop/validate a computational MDI-releasing-delivery model and to evaluate the device actuation effects on the dose distribution with the newly developed model. An integrated MDI–mouth–lung (G9) geometry was developed. An albuterol MDI with the chlorofluorocarbon propellant was simulated with polydisperse aerosol size distribution measured by laser light scatter and aerosol discharge velocity derived from measurements taken while using a phase Doppler anemometry. The highly transient, multiscale airflow and droplet dynamics were simulated by using large eddy simulation (LES) and Lagrangian tracking with sufficiently fine computation mesh. A high-speed camera imaging of the MDI plume formation was conducted and compared with LES predictions. The aerosol discharge velocity at the MDI orifice was reversely determined to be 40 m/s based on the phase Doppler anemometry (PDA) measurements at two different locations from the mouthpiece. The LES-predicted instantaneous vortex structures and corresponding spray clouds resembled each other. There are three phases of the MDI plume evolution (discharging, dispersion, and dispensing), each with distinct features regardless of the actuation time. Good agreement was achieved between the predicted and measured doses in both the device, mouth–throat, and lung. Concerning the device–patient coordination, delayed MDI actuation increased drug deposition in the mouth and reduced drug delivery to the lung. Firing MDI before inhalation was found to increase drug loss in the device; however, it also reduced mouth–throat loss and increased lung doses in both the central and peripheral regions.

Leveraging statistical shape modeling in computational respiratory dynamics: Nanomedicine delivery in remodeled airways

BACKGROUND AND OBJECTIVE

Accurate knowledge of the delivered doses to the diseased site in the respiratory tract is crucial to elicit desired therapeutic outcomes. However, such information is still difficult to obtain due to inaccessibility for measurement or visualization, complex network structure, and challenges in reconstructing lung geometries with disease-invoked airway remodeling. This study presents a novel method to simulate the airway remodeling in a mouth-lung geometry extending to G9.

METHODS

Statistical shape modeling was used to extract morphological features from a lung geometry database and four new models (i.e., M1-M4) were generated with parameter-controlled dilated/constricted bronchioles in the left-lower (LL) lung. The variations in airflow and particle deposition due to the airway remodeling were simulated using a well-tested k-ω turbulence model and a Lagrangian tracking approach.

RESULTS

Significant variations in flow partitions between the lower and upper lobes of the left lung, as well as between the left and right lungs. The flow partition into the LL lobe varied by 10-fold between the most dilated and constricted models in this study. Significantly lower doses were also predicted on the surface of the constricted LL bronchioles G4-G9, as well as into the peripheral airways beyond G9. However, the total dosimetry in the mouth-lung geometry (up to G9) exhibited low sensitivity to the LL lobar remodeling. Results in this study suggest that the optimal nanomedicine should be 2-10 nm in diameter if targeted at the constricted bronchioles G4-G9 as in topical inhalation therapy but should be larger than 20 nm if targeted at the alveolar region as in systemic therapy.

CONCLUSIONS

This study highlights the large dose variability from local airway remodeling and the need to consider these variations in the treatment planning for pneumonia and other obstructive respiratory diseases.

Inhalation dosimetry of nasally inhaled respiratory aerosols in the human respiratory tract with locally remodeled conducting lungs

Objective: Respiratory diseases are often accompanied by alterations to airway morphology. However, inhalation dosimetry data in remodeled airways are scarce due to the challenges in reconstructing diseased respiratory morphologies. This study aims to study the airway remodeling effects on the inhalation dosimetry of nasally inhaled nanoparticles in a nose-lung geometry that extends to G9 (ninth generation).Materials and methods: Statistical shape modeling was used to develop four diseased lung models with varying levels of bronchiolar dilation/constriction in the left-lower (LL) lobe (i.e. M1-M4). Respiratory airflow and particle deposition were simulated using a low Reynolds number k-ω turbulence model and a Lagrangian tracking approach.Results: Significant discrepancies were observed in the flow partitions between the left and right lungs, as well as between the lower and upper lobes of the left lung, which changed by 10-fold between the most dilated and constricted models.Much lower doses were predicted on the surface of the constricted LL bronchioles G4-G9, as well as into the peripheral airways beyond G9 of the LL lung. However, the LL lobar remodeling had little effect on the dosimetry in the nasopharynx, as well as on the total dosimetry in the nose-lung geometry (up to G9).Conclusion: It is suggested that airway remodeling may pose a higher viral infection risk to the host by redistributing the inhaled viruses to healthy lung lobes. Airway remodeling effects should also be considered in the treatment planning of inhalation therapies, not only because of the dosimetry variation from altered lung morphology but also its evolution as the disease progresses.

Simple low dose radiography allows precise lung volume assessment in mice

X-ray based lung function (XLF) as a planar method uses dramatically less X-ray dose than computed tomography (CT) but so far lacked the ability to relate its parameters to pulmonary air volume. The purpose of this study was to calibrate the functional constituents of XLF that are biomedically decipherable and directly comparable to that of micro-CT and whole-body plethysmography (WBP). Here, we developed a unique set-up for simultaneous assessment of lung function and volume using XLF, micro-CT and WBP on healthy mice. Our results reveal a strong correlation of lung volumes obtained from radiographic XLF and micro-CT and demonstrate that XLF is superior to WBP in sensitivity and precision to assess lung volumes. Importantly, XLF measurement uses only a fraction of the radiation dose and acquisition time required for CT. Therefore, the redefined XLF approach is a promising tool for preclinical longitudinal studies with a substantial potential of clinical translation.

[Functional respiratory imaging after neostigmine- or sugammadex-enhanced recovery from neuromuscular blockade in the anesthetised rat: a randomised controlled pilot study]

OBJECTIVES

Reductions in diaphragm activity are associated with the postoperative development of atelectasis. Neostigmine reversal is also associated with increased atelectasis. We assessed the effects of neostigmine, sugammadex, and spontaneous reversal on regional lung ventilation and airway flow.

METHODS

Six Sprague-Dawley rats were paralysed with rocuronium and mechanically ventilated until recovery of the train-of-four ratio to 0.5. We administered neostigmine (0.06mg.kg^-1), sugammadex (15mg.kg^-1), or saline (n=2 per group). Computed tomography scans were obtained during the breathing cycle. Three-dimensional models of lung lobes were generated using functional respiratory imaging technology, and lobar volumes were calculated during the breathing cycle. The diaphragmatic surface was segmented for the end-expiratory and end-inspiratory scans. The total change in volume was reported by the lung volume change from the end-expiratory scan to the end-inspiratory scan. Chest wall movement was defined as the lung volume change minus the volume change that resulted from diaphragm excursion.

RESULTS

The two rats that received neostigmine exhibited a smaller relative contribution of diaphragm movement to the total change in lung volume compared with the two rats that received sugammadex or saline (chest wall contribution (%): 26.69 and 25.55 for neostigmine; -2.77 and 15.98 for sugammadex; 18.82 and 10.30 for saline).

CONCLUSION

This pilot study in rats demonstrated an increased relative contribution of chest wall expansion after neostigmine compared with sugammadex or saline. This smaller relative contribution of diaphragm movement may be explained by a neostigmine-induced decrease in phrenic nerve activity or by remaining occupied acetylcholine receptors after neostigmine.

Longitudinal assessment of bleomycin-induced lung fibrosis by Micro-CT correlates with histological evaluation in mice

Background

The intratracheal instillation of bleomycin in mice induces early damage to alveolar epithelial cells and development of inflammation followed by fibrotic tissue changes and represents the most widely used model of pulmonary fibrosis to investigate human IPF.

Histopathology is the gold standard for assessing lung fibrosis in rodents, however it precludes repeated and longitudinal measurements of disease progression and does not provide information on spatial and temporal distribution of tissue damage.

Here we investigated the use of the Micro-CT technique to allow the evaluation of disease onset and progression at different time-points in the mouse bleomycin model of lung fibrosis. Micro-CT was throughout coupled with histological analysis for the validation of the imaging results.

Methods

In bleomycin-instilled and control mice, airways and lung morphology changes were assessed and reconstructed at baseline, 7, 14 and 21 days post-treatment based on Micro-CT images. Ashcroft score, percentage of collagen content and percentage of alveolar air area were detected on lung slides processed by histology and subsequently compared with Micro-CT parameters.

Results

Extent (%) of fibrosis measured by Micro-CT correlated with Ashcroft score, the percentage of collagen content and the percentage of alveolar air area (r² = 0.91; 0.77; 0.94, respectively). Distal airway radius also correlated with the Ashcroft score, the collagen content and alveolar air area percentage (r² = 0.89; 0.78; 0.98, respectively).

Conclusions

Micro-CT data were in good agreement with histological read-outs as micro-CT was able to quantify effectively and non-invasively disease progression longitudinally and to reduce the variability and number of animals used to assess the damage. This suggests that this technique is a powerful tool for understanding experimental pulmonary fibrosis and that its use could translate into a more efficient drug discovery process, also helping to fill the gap between preclinical setting and clinical practice.

Airway tree segmentation in serial block-face cryomicrotome images of rat lungs

A highly automated method for the segmentation of airways in the serial block-face cryomicrotome images of rat lungs is presented. First, a point inside of the trachea is manually specified. Then, a set of candidate airway centerline points is automatically identified. By utilizing a novel path extraction method, a centerline path between the root of the airway tree and each point in the set of candidate centerline points is obtained. Local disturbances are robustly handled by a novel path extraction approach, which avoids the shortcut problem of standard minimum cost path algorithms. The union of all centerline paths is utilized to generate an initial airway tree structure, and a pruning algorithm is applied to automatically remove erroneous subtrees or branches. Finally, a surface segmentation method is used to obtain the airway lumen. The method was validated on five image volumes of Sprague-Dawley rats. Based on an expert-generated independent standard, an assessment of airway identification and lumen segmentation performance was conducted. The average of airway detection sensitivity was 87.4% with a 95% confidence interval (CI) of (84.9, 88.6)%. A plot of sensitivity as a function of airway radius is provided. The combined estimate of airway detection specificity was 100% with a 95% CI of (99.4, 100)%. The average number and diameter of terminal airway branches was 1179 and 159 μm, respectively. Segmentation results include airways up to 31 generations. The regression intercept and slope of airway radius measurements derived from final segmentations were estimated to be 7.22 μm and 1.005, respectively. The developed approach enables the quantitative studies of physiology and lung diseases in rats, requiring detailed geometric airway models.

Rat airway morphometry measured from in situ MRI-based geometric models

Rodents have been widely used to study the environmental or therapeutic impact of inhaled particles. Knowledge of airway morphometry is essential in assessing geometric influence on aerosol deposition and in developing accurate lung models of aerosol transport. Previous morphometric studies of the rat lung performed ex situ provided high-resolution measurements (50-125 μm). However, it is unclear how the overall geometry of these casts might have differed from the natural in situ appearance. In this study, four male Wistar rat (268 ± 14 g) lungs were filled sequentially with perfluorocarbon and phosphate-buffered saline before being imaged in situ in a 7-T magnetic resonance (MR) scanner at a resolution of 0.2 × 0.2 × 0.27 mm. Airway length, diameter, gravitational, bifurcation, and rotational angles were measured for the first four airway generations from 3D geometric models built from the MR images. Minor interanimal variability [expressed by the relative standard deviation RSD (=SD/mean)] was found for length (0.18 ± 0.07), diameter (0.15 ± 0.15), and gravitational angle (0.12 ± 0.06). One rat model was extended to 16 airway generations. Organization of the airways using a diameter-defined Strahler ordering method resulted in lower interorder variability than conventional generation-based grouping for both diameter (RSD = 0.12 vs. 0.42) and length (0.16 vs. 0.67). Gravitational and rotational angles averaged 82.9 ± 37.9° and 53.6 ± 24.1°, respectively. Finally, the major daughter branch bifurcated at a smaller angle (19.3 ± 14.6°) than the minor branch (60.5 ± 19.4°). These data represent the most comprehensive set of rodent in situ measurements to date and can be used readily in computational studies of lung function and aerosol exposure.

Development of an in vitro model to assess deposition of aerosol particles in a representative replica of the rat's respiratory tract

BACKGROUND

Small rodents continue to be the mainstay for the assessment of pharmacological and toxicological data of inhaled therapeutics. For meaningful interpretation of the results information about deposition of aerosol particles in the respiratory tract is warranted, but not trivial to obtain for animals with nose-only aerosol exposure. The purpose of this study was to develop and evaluate a general method to characterize the deposition of inhaled test particles in an in vitro model (IVR) of the rat's respiratory tract.

METHODS

A highly detailed, realistic and representative image using micro-CT scanning technology was obtained and the generated morphological data was used to construct a plastic replica of the average rat respiratory tract. The model was connected to a rodent ventilator, which allowed the breathing frequency (f, min(-1)) and tidal volume (V(T), mL) to be varied as required. Polydisperse fluorescent microsphere particles with an average mass median aerodynamic diameter (MMAD) of 3.1 μm and geometric standard deviation (GSD) of 2.2 μm were used as model compound.

RESULTS

Comparison of the experimental data for total and regional deposition levels with predicted outputs using the in silico MPPD model showed reasonably good relative agreement between the two models. The predictions were closest to the experimental values when default respiratory conditions of f=102 breaths/min and V(T) of 2.0 mL were used. Moreover, the IVR model revealed good correlation with published in vivo data.

CONCLUSION

Using the IVR model allows an easy, fast and reasonably precise estimation of the inhaled dose in rodent inhalation studies. The IVR has the potential to be used along with live rats in inhalation exposure studies, and thus provides the unique possibility to run an internal standard for dose deposition in the respiratory tract in each inhalation experiment. This should contribute to enable a greater understanding of drug pharmacokinetics and dynamics in rats and may improve dose extrapolation from animals to humans.

A case series on lung deposition analysis of inhaled medication using functional imaging based computational fluid dynamics in asthmatic patients: effect of upper airway morphology and comparison with in vivo data

CONTEXT

Asthma affects 20 million Americans resulting in an economic burden of approximately $18 billion in the US alone (Allergies and Asthma Foundation 2000; National Center for Environmental Health (NCEH) 1999). Research studies based on differences in patient-specific airway morphology for asthma and the associated effect on deposition of inhaled aerosols are currently not available in the literature. Therefore, the role of morphological variations such as upper airway (extrathoracic) occlusion is not well documented.

OBJECTIVE

Functional imaging based computational fluid dynamics (CFD) of the respiratory airways for five asthmatic subjects is performed in this study using computed tomography (CT) based patient-specific airway models and boundary conditions.

METHODS

CT scans for 5 asthma patients were used to reconstruct 3D lung models using segmentation software. An averaged inhalation profile and patient-specific lobar flow distribution were used to perform the simulation. The simulations were used to obtain deposition for BDP/Formoterol® HFA pMDI in the patient-specific airway models.

RESULTS

The lung deposition obtained using CFD was in excellent agreement with available in vivo data using the same product. Specifically, CFD resulted in 30% lung deposition, whereas in vivo lung deposition was reported to be approximately 31%.

CONCLUSION

It was concluded that a combination of patient-specific airway models and lobar boundary conditions can be used to obtain accurate lung deposition estimates. Lower lung deposition can be expected for patients with higher extrathoracic resistance. Novel respiratory drug delivery devices need to accommodate population sub-groups based on these morphological and anatomical differences in addition to subject age.

Functional imaging using computer methods to compare the effect of salbutamol and ipratropium bromide in patient-specific airway models of COPD

BACKGROUND

Salbutamol and ipratropium bromide improve lung function in patients with chronic obstructive pulmonary disease (COPD). However, their bronchodilating effect has not yet been compared in the central and distal airways. Functional imaging using computational fluid dynamics offers the possibility of making such a comparison. The objective of this study was to assess the effects of salbutamol and ipratropium bromide on the geometry and computational fluid dynamics-based resistance of the central and distal airways.

METHODS

Five patients with Global Initiative for Chronic Obstructive Lung Disease Stage III COPD were randomized to a single dose of salbutamol or ipratropium bromide in a crossover manner with a 1-week interval between treatments. Patients underwent lung function testing and a multislice computed tomography scan of the thorax that was used for functional imaging. Two hours after dosing, the patients again underwent lung function tests and repeat computed tomography.

RESULTS

Lung function parameters, including forced expiratory volume in 1 second, vital capacity, overall airway resistance, and specific airway resistance, changed significantly after administration of each product. On functional imaging, the bronchodilating effect was greater in the distal airways, with a corresponding drop in airway resistance, compared with the central airways. Salbutamol and ipratropium bromide were equally effective at first glance when looking at lung function tests, but when viewed in more detail with functional imaging, hyporesponsiveness could be shown for salbutamol in one patient. Salbutamol was more effective in the other patients.

CONCLUSION

This pilot study gives an innovative insight into the modes of action of salbutamol and ipratropium bromide in patients with COPD, using the new techniques of functional imaging and computational fluid dynamics.

The effects of long-term noninvasive ventilation in hypercapnic COPD patients: a randomized controlled pilot study

Introduction

Noninvasive ventilation (NIV) is a well-established treatment for acute-on- chronic respiratory failure in hypercapnic COPD patients. Less is known about the effects of a long-term treatment with NIV in hypercapnic COPD patients and about the factors that may predict response in terms of improved oxygenation and lowered CO₂ retention.

Methods

In this study, we randomized 15 patients to a routine pharmacological treatment (n = 5, age 66 [standard deviation ± 6] years, FEV₁ 30.5 [±5.1] %pred, PaO₂ 65 [±6] mmHg, PaCO₂ 52.4 [±6.0] mmHg) or to a routine treatment and NIV (using the Synchrony BiPAP device [Respironics, Inc, Murrsville, PA]) (n = 10, age 65 [±7] years, FEV₁ 29.5 [±9.0] %pred, PaO₂ 59 [±13] mmHg, PaCO₂ 55.4 [±7.7] mmHg) for 6 months. We looked at arterial blood gasses, lung function parameters and performed a low-dose computed tomography of the thorax, which was later used for segmentation (providing lobe and airway volumes, iVlobe and iVaw) and post-processing with computer methods (providing airway resistance, iRaw) giving overall a functional image of the separate airways and lobes.

Results

In both groups there was a nonsignificant change in FEV₁ (NIV group 29.5 [9.0] to 38.5 [14.6] %pred, control group 30.5 [5.1] to 36.8 [8.7] mmHg). PaCO₂ dropped significantly only in the NIV group (NIV: 55.4 [7.7] → 44.5 [4.70], P = 0.0076; control: 52.4 [6.0] → 47.6 [8.2], NS). Patients actively treated with NIV developed a more inhomogeneous redistribution of mass flow than control patients. Subsequent analysis indicated that in NIV-treated patients that improve their blood gases, mass flow was also redistributed towards areas with higher vessel density and less emphysema, indicating that flow was redistributed towards areas with better perfusion. There was a highly significant correlation between the % increase in mass flow towards lobes with a blood vessel density of >9% and the increase in PaO₂. Improved ventilation–perfusion match and recruitment of previously occluded small airways can explain the improvement in blood gases.

Conclusion

We can conclude that in hypercapnic COPD patients treated with long-term NIV over 6 months, a mass flow redistribution occurs, providing a better ventilation–perfusion match and hence better blood gases and lung function. Control patients improve homogeneously in iVaw and iRaw, without improvement in gas exchange since there is no improved ventilation/perfusion ratio or increased alveolar ventilation. These differences in response can be detected through functional imaging, which gives a more detailed report on regional lung volumes and resistances than classical lung function tests do. Possibly only patients with localized small airway disease are good candidates for long-term NIV treatment. To confirm this and to see if better arterial blood gases also lead to better health related quality of life and longer survival, we have to study a larger population.

Comparison of hyperpolarized (3)He MRI with Xe-enhanced computed tomography imaging for ventilation mapping of rat lung

Lung ventilation was mapped in five healthy Brown Norway rats (210-377 g) using both hyperpolarized (3)He MRI and Xe-enhanced computed tomography (Xe-CT) under similar ventilator conditions. Whole-lung measurements of ventilation r obtained with (3)He MRI were not significantly different from those obtained from Xe-CT (p = 0.1875 by Wilcoxon matched pairs test). The ventilation parameter r is defined as the fraction of refreshed gas per unit volume per breath. Regional ventilation was also measured in four regions of the lung using both methods. A two-tailed paired t-test was performed for each region, yielding p > 0.05 for all but the upper portion of the right lung. The distribution of regional ventilation was evaluated by calculating ventilation gradients in the superior/inferior (S/I) direction. The average S/I gradient obtained using the (3)He MRI method was found to be 0.17 ± 0.04 cm(-1) , whereas the average S/I gradient obtained using the Xe-CT method was found to be 0.016 ± 0.005 cm(-1) . In general, S/I ventilation gradients obtained from both methods were significantly different from each other (p = 0.0019 by two-tailed paired t-test). These regional differences in ventilation measurements may be caused by the manner in which the gas contrast agents distribute physiologically and/or by the imaging modality.