Aerosols in healthy and emphysematous in silico pulmonary acinar rat models

Design of a multilayer lung chip with multigenerational alveolar ducts to investigate the inhaled particle deposition

We present the development and application of a multilayer microfluidic lung chip designed to accurately replicate the human respiratory bronchi, providing an innovative platform for controlled particle deposition in the lung. By employing a quantitative control method of fluid velocity through the deformation of an elastic PDMS membrane, this platform mimics the passive breathing process in humans and allows for precise simulation of the respiration cycle. We utilized time-lapse photography of fluorescent particles in a water/glycerol solution to qualitatively observe fluid morphology in the channel, while a chip-aerosol exposure device combined with microscopy imaging was employed to visualise aerosol deposition. Both experimental and numerical simulation results showed that particle concentration decreased towards the distal generations of the lung, and that changes in breathing pattern significantly affected particle deposition trends. Furthermore, we found that increasing the residence time of particles in the channel facilitated deeper particle deposition, achievable by adjusting parameters such as breath-hold time, exhalation time, respiration cycle length, and tidal volume. The proposed microfluidic lung chip device has significant potential for future research in respiratory health and inhaled drug delivery, providing an efficient, cost-effective, and ethical alternative to traditional in vivo studies.

Inhalation toxicity of polyhexamethylene guanidine-phosphate in rats: A 4-week inhalation exposure and 24-week recovery period study

Humidifier disinfectant (HD) is a causative agent of atypical lung injury reported in 2011 in South Korea, and various diseases caused by HD after exposure cessation have been reported to date. However, there is limited research on most of the reported diseases in terms of their association with HD exposure, and information on the progression of diseases caused by HD exposure is also limited. Therefore, we investigated the effects of HD inhalation on the body in rats. Rats were exposed to 0.15, 0.50, and 1.60 mg/m³ polyhexamethylene guanidine-phosphate (PHMG-p), which is the major component of HDs and most closely related to HD-associated lung injury. We conducted necropsy four times during the recovery period (0, 4, 12, and 24 weeks) and evaluated general systemic toxicities. There were significant changes in the mortality rate, body weight, and food consumption in the PHMG-p-exposed groups. Hematology revealed changes in hemoglobin level, hematocrit, red blood cell, reticulocyte, and white blood cell counts until 12 weeks of the recovery period. PHMG-p induced a delay in prothrombin time until 12 weeks of the recovery period. The aspartate aminotransferase, alanine aminotransferase, total bilirubin, and triglyceride levels were higher in the PHMG-p-exposed groups than in the control group at week 4 of the recovery period, and these parameters normalized after 12 weeks of the recovery period. Histopathological examination in PHMG-p exposed groups revealed several changes in the lungs, including the presence of alveolar macrophages, chronic inflammation, squamous metaplasia, alveolar emphysema, and pulmonary fibrosis. The severity of these symptoms was maintained or exacerbated till 24 weeks. Overall, PHMG-p inhalation can induce irreversible histological changes in the lungs and cause various types of damage throughout the body, even after exposure ends.

Microflows in two-generation alveolar cells at an acinar bifurcation

The alveolus is a basic functional unit of the human respiratory system, and the airflow in the alveoli plays an important role in determining the transport and deposition of particulate matter, which is crucial for inhaled disease diagnosis and drug delivery. In the present study, taking advantage of the precise control ability of the microfluidic technique, a rhythmically expanding alveolar chip with multiple alveoli in two generations is designed and both the geometric and kinematic similarities are matched with the real human respiration system. With the help of a micro-PIV measurement system, the microflow patterns inside each alveolus can be studied. The observed vortex and radial flow patterns and the discovery of stagnant saddle points are similar to those captured in our previous platform with only one alveolus [Lv et al., Lab Chip 20, 2394-2402 (2020)]. However, the interactions between multiple alveoli also uncover new phenomena, such as the finding of stagnant saddle points in non-vortex flow patterns and significant differences in the flow pattern around the points between the time of T/4 and 3T/4. The obtained results could enrich the understanding of microflow in a whole alveolar tree with multiple generations.

Recent advances in the understanding of alveolar flow

Understanding the dynamics of airflow in alveoli and its effect on the behavior of particle transport and deposition is important for understanding lung functions and the cause of many lung diseases. The studies on these areas have drawn substantial attention over the last few decades. This Review discusses the recent progress in the investigation of behavior of airflow in alveoli. The information obtained from studies on the structure of the lung airway tree and alveolar topology is provided first. The current research progress on the modeling of alveoli is then reviewed. The alveolar cell parameters at different generation of branches, issues to model real alveolar flow, and the current numerical and experimental approaches are discussed. The findings on flow behavior, in particular, flow patterns and the mechanism of chaotic flow generation in the alveoli are reviewed next. The different flow patterns under different geometrical and flow conditions are discussed. Finally, developments on microfluidic devices such as lung-on-a-chip devices are reviewed. The issues of current devices are discussed.

Effects of Pulmonary Fibrosis and Surface Tension on Alveolar Sac Mechanics in Diffuse Alveolar Damage

Diffuse alveolar damage (DAD) is a characteristic histopathologic pattern in most cases of acute respiratory distress syndrome and severe viral pneumonia, such as COVID-19. DAD is characterized by an acute phase with edema, hyaline membranes, and inflammation followed by an organizing phase with pulmonary fibrosis and hyperplasia. The degree of pulmonary fibrosis and surface tension is different in the pathological stages of DAD. The effects of pulmonary fibrosis and surface tension on alveolar sac mechanics in DAD are investigated by using the fluid-structure interaction (FSI) method. The human pulmonary alveolus is idealized by a three-dimensional honeycomb-like geometry, with alveolar geometries approximated as closely packed 14-sided polygons. A dynamic compression-relaxation model for surface tension effects is adopted. Compared to a healthy model, DAD models are created by increasing the tissue thickness and decreasing the concentration of the surfactant. The FSI results show that pulmonary fibrosis is more influential than the surface tension on flow rate, volume, P-V loop, and resistance. The lungs of the disease models become stiffer than those of the healthy models. According to the P-V loop results, the surface tension plays a more important role in hysteresis than the material nonlinearity of the lung tissue. Our study demonstrates the differences in air flow and lung function on the alveolar sacs between the healthy and DAD models.

A whole lung in silico model to estimate age dependent particle dosimetry

Anatomical and physiological changes alter airflow characteristics and aerosol distribution in the developing lung. Correlation between age and aerosol dosimetry is needed, specifically because youth are more susceptible to medication side effects. In this study, we estimate aerosol dosages (particle diameters of 1, 3, and 5 [Formula: see text]m) in a 3 month-old infant, a 6 year-old child, and a 36 year-old adult by performing whole lung subject-specific particle simulations throughout respiration. For 3 [Formula: see text]m diameter particles we estimate total deposition as 88, 73, and [Formula: see text] and the conducting versus respiratory deposition ratios as 4.0, 0.5, and 0.4 for the infant, child, and adult, respectively. Due to their lower tidal volumes and functional residual capacities the deposited mass is smaller while the tissue concentrations are larger in the infant and child subjects, compared to the adult. Furthermore, we find that dose cannot be predicted by simply scaling by tidal volumes. These results highlight the need for additional clinical and computational studies that investigate the efficiency of treatment, while optimizing dosage levels in order to alleviate side effects, in youth.

Septal destruction enhances chaotic mixing and increases cellular doses of nanoparticles in emphysematous acinus

One hallmark of emphysema is the breakdown of inter-alveolar septal walls in pulmonary acini. How the acinar dosimetry of environmental aerosols varies at different stages of emphysema remains unclear; this is specifically pertinent to users of tobacco products, which is the leading cause of emphysema. The objective of this study is to systematically assess the impacts of septal destruction on the behavior and fate of nanoparticles (1–800 nm) in a pyramid-shaped sub-acinar model consisting of 496 alveoli. Four diseased geometry variants were created by gradually removing the septal walls from the base model. Particle motions within the acinar region were tracked for particles raging 1–800 nm at four emphysema stages using a well-tested Lagrangian tracking model. Both spatial profile and temporal variation of particle deposition were predicted in healthy and diseased sub-acinar geometries on both a total and regional basis. Results show large differences in airflow and particle dynamics among different emphysema stages. Large differences in particle dynamics are also observed among different particle sizes, with one order of magnitude’s variation in the speeds of particles of 1, 10, and 200 nm. The destruction of septal walls also changed the deposition mechanisms, shifting from connective diffusion to chaotic mixing with emphysema progression. The sub-acinar dosimetry became less sensitive to particle size variation with more septal destructions. The lowest retention rate was found at 200–500 nm in the healthy sub-acinar geometry, but at 800 nm in all emphysematous models considered. The acinus-averaged dose for nanoparticles (1–800 nm) increases with aggravating septal destructions, indicating an even higher risk to the acinus at later emphysema stages.

Aerosol delivery into small anatomical airway model through spontaneous engineered breathing

Pulmonary administration is a noninvasive drug delivery method that, in contrast to systemic administration, reduces drug dosage and possible side effects. Numerous testing models, such as impingers and impactors, have previously been developed to evaluate the fate of inhaled drugs. However, such models are limited by the lack of information regarding several factors, such as pulmonary morphology and breathing motion, which are required to fully interpret actual inhaled-drug deposition profiles within the human respiratory tract. In this study, a spontaneous breathing-lung model that integrates branched morphology and deformable alveolar features was constructed using a multilayered fabrication technology to mimic the complex environment of the human lower respiratory tract. The developed model could emulate cyclic and spontaneous breathing motions to inhale and exhale aerosols generated by a nebulizer under diseaselike conditions. Results of this research demonstrate that aerosols (4.2 μm) could reach up to the deeper lung regions (generation 19 of the branched lung structure) within the obstructivelike model, whereas lesser penetration (generation 17) was observed when using the restrictivelike model. The proposed breathing-lung model can serve as a testing platform to provide a comprehensive understanding of the pharmacokinetics of pulmonary drugs within the lower lungs.

Characterization of air flow and lung function in the pulmonary acinus by fluid-structure interaction in idiopathic interstitial pneumonias

Background and objective

The idiopathic interstitial pneumonias (IIPs) are diffuse parenchymal lung disorders that are associated with substantial morbidity and mortality. Early diagnosis and disease stratification of IIP patients are important because these are related with the treatment and prognosis. Idiopathic pulmonary fibrosis (IPF) and nonspecific interstitial pneumonia (NSIP) are two major distinctive pathologic patterns of pulmonary fibrosis. We researched the application of the fluid-structure interaction (FSI) to the respiratory system and compared the pulmonary acinus mechanics and functions in healthy and IIP models.

Methods

The human pulmonary alveolus is idealized by a three-dimensional honeycomb-like geometry, and a fluid-structure interaction analysis is performed to study the normal and diseased breathing mechanics. The computational domain consists of two generations of alveolar ducts within the pulmonary acinus, with alveolar geometries approximated as closely packed 14-sided polygons.

Findings

In a normal breathing cycle, the flow rate of the healthy model is significantly larger than that of the NSIP and IPF models. Similar trends are observed for the volume change and the maximum pressure drop. The flow rate and the volume change of the NSIP are almost the same as those of IPF. The maximum pressure drop of NSIP is 5.5% larger than that of IPF. There is a 47% decrease in the pulmonary acinus compliance for the NSIP and IPF compared with that of the healthy model. The acinus resistances of NSIP and IPF are higher than those of the healthy lung by 6.4~11.2%. In particular, the pulmonary acinus resistance of the NSIP lung is higher than that of the IPF lung by 4.5%.

Conclusions

Our study demonstrates the differences of air flow and lung function in the pulmonary acinus between the healthy and the IIP models. These changes in the lung are important considerations for early diagnosis and disease stratification in patients. Patient-based geometry can to be included in the computational models in future studies.

Modeling Airflow and Particle Deposition in a Human Acinar Region

The alveolar region, encompassing millions of alveoli, is the most vital part of the lung. However, airflow behavior and particle deposition in that region are not fully understood because of the complex geometrical structure and intricate wall movement. Although recent investigations using 3D computer simulations have provided some valuable information, a realistic analysis of the air-particle dynamics in the acinar region is still lacking. So, to gain better physical insight, a physiologically inspired whole acinar model has been developed. Specifically, air sacs (i.e., alveoli) were attached as partial spheroids to the bifurcating airway ducts, while breathing-related wall deformation was included to simulate actual alveolar expansion and contraction. Current model predictions confirm previous notions that the location of the alveoli greatly influences the alveolar flow pattern, with recirculating flow dominant in the proximal lung region. In the midalveolar lung generations, the intensity of the recirculating flow inside alveoli decreases while radial flow increases. In the distal alveolar region, the flow pattern is completely radial. The micron/submicron particle simulation results, employing the Euler–Lagrange modeling approach, indicate that deposition depends on the inhalation conditions and particle size. Specifically, the particle deposition rate in the alveolar region increases with higher inhalation tidal volume and particle diameter. Compared to previous acinar models, the present system takes into account the entire acinar region, including both partially alveolated respiratory bronchioles as well the fully alveolated distal airways and alveolar sacs. In addition, the alveolar expansion and contraction have been calculated based on physiological breathing conditions which make it easy to compare and validate model results with in vivo lung deposition measurements. Thus, the current work can be readily incorporated into human whole-lung airway models to simulate/predict the flow dynamics of toxic or therapeutic aerosols.

Use of computational fluid dynamics deposition modeling in respiratory drug delivery

INTRODUCTION

Respiratory drug delivery is a surprisingly complex process with a number of physical and biological challenges. Computational fluid dynamics (CFD) is a scientific simulation technique that is capable of providing spatially and temporally resolved predictions of many aspects related to respiratory drug delivery from initial aerosol formation through respiratory cellular drug absorption.

AREAS COVERED

This review article focuses on CFD-based deposition modeling applied to pharmaceutical aerosols. Areas covered include the development of new complete-airway CFD deposition models and the application of these models to develop a next-generation of respiratory drug delivery strategies.

EXPERT OPINION

Complete-airway deposition modeling is a valuable research tool that can improve our understanding of pharmaceutical aerosol delivery and is already supporting medical hypotheses, such as the expected under-treatment of the small airways in asthma. These complete-airway models are also being used to advance next-generation aerosol delivery strategies, like controlled condensational growth. We envision future applications of CFD deposition modeling to reduce the need for human subject testing in developing new devices and formulations, to help establish bioequivalence for the accelerated approval of generic inhalers, and to provide valuable new insights related to drug dissolution and clearance leading to microdosimetry maps of drug absorption.

Simulation of Airflow in an Idealized Emphysematous Human Acinus

Emphysema is the permanent enlargement of air spaces in the respiratory regions of the lung due to destruction of the inter-alveolar septa. The progressive coalescence of alveoli and alveolar ducts into larger airspaces leads to the disruption of normal airway wall motion and airflow rates within the pulmonary acinus. To contribute to the understanding of the individual effects of emphysema during its earliest stages, computational fluid dynamics (CFD) simulations of airflow in mathematically derived models of the pulmonary acinus were performed. The here generated computational domain consists of two generations of alveolar ducts within the pulmonary acinus, with alveolar geometries approximated as closely packed, 14-sided polygons. Physiologically realistic airflow rates and wall motions were used to study airflow patterns within subsequent generations of alveolar ducts during the inspiratory and expiratory phases of the breathing cycle. The effects of progressive emphysema on the airway wall motion and flow rates were simulated by sequentially removing all alveolar septa within each alveolar duct. Parametric studies were presented to independently assess the relative influence of progressive septal destruction of airway motion and flow rates. The results illustrate that septal destruction lowers the flow resistance through the alveolar ducts but has little influence on the mass transport of oxygen into the alveoli. Septal destruction has a net effect on the flow field by favoring the development of recirculatory flow patterns in individual alveoli.

One (sub-)acinus for all: Fate of inhaled aerosols in heterogeneous pulmonary acinar structures

Computational Fluid Dynamics (CFD) have offered an attractive gateway to investigate in silico respiratory flows and aerosol transport in the depths of the lungs. Yet, not only do existing models lack sufficient anatomical realism in capturing the heterogeneity and morphometry of the acinar environment, numerical simulations have been widely restricted to domains capturing a mere few percent of a single acinus. Here, we present to the best of our knowledge the most detailed and comprehensive in silico simulations to date on the fate of aerosols in the acinar depths. Our heterogeneous acinar domains represent complete sub-acinar models (i.e. 1/8th of a full acinus) based on the recent algorithm of Koshiyama & Wada (2015), capturing statistics of human acinar morphometry (Ochs et al. 2004). Our simulations deliver high-resolution, 3D spatial-temporal data on aerosol transport and deposition, emphasizing how variances in acinar heterogeneity only play a minor role in determining general deposition outcomes. With such tools at hand, we revisit whole-lung deposition predictions (i.e. ICRP) based on past 1D lung models. While our findings under quiet breathing substantiate general deposition trends obtained with past predictions in the alveolar regions, we underscore how deposition fractions are anticipated to increase, in particular during deep inhalation. For such inhalation maneuver, our simulations support the notion of significantly augmented deposition for all aerosol sizes (0.005-5.0μm). Overall, our efforts not only help consolidate our mechanistic understanding of inhaled aerosol transport in the acinar depths but also continue to bridge the gap between "bottom-up" in silico models and regional deposition predictions from whole-lung models. Such quantifications provide what is deemed more accurate deposition predictions in morphometrically-faithful models and are particularly useful in assessing inhalation strategies for deep airway deposition (e.g. systemic delivery).

Aerosol transport throughout inspiration and expiration in the pulmonary airways

Little is known about transport throughout the respiration cycle in the conducting airways. It is challenging to appropriately describe the time-dependent number of particles entering back into the model during exhalation. Modeling the entire lung is not feasible; therefore, multidomain methods must be used. Here, we present a new framework that is designed to simulate particles throughout the respiration cycle, incorporating realistic airway geometry and respiration. This framework is applied for a healthy rat lung exposed to  ∼ 1μm diameter particles, chosen to facilitate parameterization and validation. The flow field is calculated in the conducting airways (3D domain) by solving the incompressible Navier-Stokes equations with experimentally derived boundary conditions. Particles are tracked throughout inspiration by solving a modified Maxey-Riley equation. Next, we pass the time-dependent particle concentrations exiting the 3D model to the 1D volume conservation and advection-diffusion models (1D domain). Once the 1D models are solved, we prescribe the time-dependent number of particles entering back into the 3D airways to again solve for 3D transport. The coupled simulations highlight that about twice as many particles deposit during inhalation compared to exhalation for the entire lung. In contrast to inhalation, where most particles deposit at the bifurcation zones, particles deposit relatively uniformly on the gravitationally dependent side of the 3D airways during exhalation. Strong agreement to previously collected regional experimental data is shown, as the 1D models account for lobe-dependent morphology. This framework may be applied to investigate dosimetry in other species and pathological lungs.

Aging effects on airflow dynamics and lung function in human bronchioles

Background and objective

The mortality rate for patients requiring mechanical ventilation is about 35% and this rate increases to about 53% for the elderly. In general, with increasing age, the dynamic lung function and respiratory mechanics are compromised, and several experiments are being conducted to estimate these changes and understand the underlying mechanisms to better treat elderly patients.

Materials and methods

Human tracheobronchial (G1 ~ G9), bronchioles (G10 ~ G22) and alveolar sacs (G23) geometric models were developed based on reported anatomical dimensions for a 50 and an 80-year-old subject. The aged model was developed by altering the geometry and material properties of the model developed for the 50-year-old. Computational simulations using coupled fluid-solid analysis were performed for geometric models of bronchioles and alveolar sacs under mechanical ventilation to estimate the airflow and lung function characteristics.

Findings

The airway mechanical characteristics decreased with aging, specifically a 38% pressure drop was observed for the 80-year-old as compared to the 50-year-old. The shear stress on airway walls increased with aging and the highest shear stress was observed in the 80-year-old during inhalation. A 50% increase in peak strain was observed for the 80-year-old as compared to the 50-year-old during exhalation. The simulation results indicate that there is a 41% increase in lung compliance and a 35%-50% change in airway mechanical characteristics for the 80-year-old in comparison to the 50-year-old. Overall, the airway mechanical characteristics as well as lung function are compromised due to aging.

Conclusion

Our study demonstrates and quantifies the effects of aging on the airflow dynamics and lung capacity. These changes in the aging lung are important considerations for mechanical ventilation parameters in elderly patients. Realistic geometry and material properties need to be included in the computational models in future studies.

The role of anisotropic expansion for pulmonary acinar aerosol deposition

Lung deformations at the local pulmonary acinar scale are intrinsically anisotropic. Despite progress in imaging modalities, the true heterogeneous nature of acinar expansion during breathing remains controversial, where our understanding of inhaled aerosol deposition still widely emanates from studies under self-similar, isotropic wall motions. Building on recent 3D models of multi-generation acinar networks, we explore in numerical simulations how different hypothesized scenarios of anisotropic expansion influence deposition outcomes of inhaled aerosols in the acinar depths. While the broader range of particles acknowledged to reach the acinar region (dp=0.005-5.0μm) are largely unaffected by the details of anisotropic expansion under tidal breathing, our results suggest nevertheless that anisotropy modulates the deposition sites and fractions for a narrow band of sub-micron particles (dp~0.5-0.75μm), where the fate of aerosols is greatly intertwined with local convective flows. Our findings underscore how intrinsic aerosol motion (i.e. diffusion, sedimentation) undermines the role of anisotropic wall expansion that is often attributed in determining aerosol mixing and acinar deposition.