Particle dynamics and deposition in true-scale pulmonary acinar models

Substance deposition assessment in obstructed pulmonary system through numerical characterization of airflow and inhaled particles attributes

BACKGROUND

Chronic obstructive pulmonary disease (COPD) and asthma are considered as the two most widespread obstructive lung diseases, whereas they affect more than 500 million people worldwide. Unfortunately, the requirement for detailed geometric models of the lungs in combination with the increased computational resources needed for the simulation of the breathing did not allow great progress to be made in the past for the better understanding of inflammatory diseases of the airways through detailed modelling approaches. In this context, computational fluid dynamics (CFD) simulations accompanied by fluid particle tracing (FPT) analysis of the inhaled ambient particles are deemed critical for lung function assessment. Also they enable the understanding of particle depositions on the airways of patients, since these accumulations may affect or lead to inflammations. In this direction, the current study conducts an initial investigation for the better comprehension of particle deposition within the lungs. More specifically, accurate models of the airways obstructions that relate to pulmonary disease are developed and a thorough assessment of the airflow behavior together with identification of the effects of inhaled particle properties, such as size and density, is conducted. Our approach presents a first step towards an effective personalization of pulmonary treatment in regards to the geometric characteristics of the lungs and the in depth understanding of airflows within the airways.

METHODS

A geometry processing technique involving contraction algorithms is established and used to employ the different respiratory arrangements associated with lung related diseases that exhibit airways obstructions. Apart from the normal lung case, two categories of obstructed cases are examined, i.e. models with obstructions in both lungs and models with narrowings in the right lung only. Precise assumptions regarding airflow and deposition fraction (DF) over various sections of the lungs are drawn by simulating these distinct incidents through the finite volume method (FVM) and particularly the CFD and FPT algorithms. Moreover, a detailed parametric analysis clarifies the effects of the particles size and density in terms of regional deposition upon several parts of the pulmonary system. In this manner, the deposition pattern of various substances can be assessed.

RESULTS

For the specific case of the unobstructed lung model most particles are detected on the right lung (48.56% of total, when the air flowrate is 12.6 L/min), a fact that is also true when obstructions arise symmetrically in both lungs (51.45% of total, when the air flowrate is 6.06 L/min and obstructions occur after the second generation). In contrast, when narrowings are developed on the right lung only, most particles are pushed on the left section (68.22% of total, when the air flowrate is 11.2 L/min) indicating that inhaled medication is generally deposited away from the areas of inflammation. This observation is useful when designing medical treatment of lung diseases. Furthermore, particles with diameters from 1 μm to 10 μm are shown to be mainly deposited on the lower airways, whereas particles with diameters of 20 μm and 30 μm are mostly accumulated in the upper airways. As a result, the current analysis indicates increased DF levels in the upper airways when the particle diameter is enlarged. Additionally, when the particles density increases from 1000 Kg/m³ to 2000 Kg/m³, the DF is enhanced on every generation and for all cases investigated herein. The results obtained by our simulations provide an accurate and quantitative estimation of all important parameters involved in lung modeling.

CONCLUSIONS

The treatment of respiratory diseases with inhaled medical substances can be advanced by the clinical use of accurate CFD and FPT simulations and specifically by evaluating the deposition of inhaled particles in a regional oriented perspective in regards to different particle sizes and particle densities. Since a drug with specific characteristics (i.e. particle size and density) exhibits maximum deposition on particular lung areas, the current study provides initial indications to a qualified physician for proper selection of medication.

Capturing complexity in pulmonary system modelling

Respiratory disease is a significant problem worldwide, and it is a problem with increasing prevalence. Pathology in the upper airways and lung is very difficult to diagnose and treat, as response to disease is often heterogeneous across patients. Computational models have long been used to help understand respiratory function, and these models have evolved alongside increases in the resolution of medical imaging and increased capability of functional imaging, advances in biological knowledge, mathematical techniques and computational power. The benefits of increasingly complex and realistic geometric and biophysical models of the respiratory system are that they are able to capture heterogeneity in patient response to disease and predict emergent function across spatial scales from the delicate alveolar structures to the whole organ level. However, with increasing complexity, models become harder to solve and in some cases harder to validate, which can reduce their impact clinically. Here, we review the evolution of complexity in computational models of the respiratory system, including successes in translation of models into the clinical arena. We also highlight major challenges in modelling the respiratory system, while making use of the evolving functional data that are available for model parameterisation and testing.

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).

Transport of ellipsoid fibers in oscillatory shear flows: Implications for aerosol deposition in deep airways

It is widely acknowledged that inhaled fibers, e.g. air pollutants and anthropogenic particulate matter, hold the ability to deposit deep into the lungs reaching the distal pulmonary acinar airways as a result of their aerodynamic properties; these particles tend to align with the flow and thus stay longer airborne relative to their spherical counterpart, due to higher drag forces that resist sedimentation. Together with a high surface-to-volume ratio, such characteristics may render non-spherical particles, and fibers in particular, potentially attractive airborne carriers for drug delivery. Until present, however, our understanding of the dynamics of inhaled aerosols in the distal regions of the lungs has been mostly limited to spherical particles. In an effort to unravel the fate of non-spherical aerosols in the pulmonary depths, we explore through numerical simulations the kinematics of ellipsoid-shaped fibers in a toy model of a straight pipe as a first step towards understanding particle dynamics in more intricate acinar geometries. Transient translational and rotational motions of micron-sized ellipsoid particles are simulated as a function of aspect ratio (AR) for laminar oscillatory shear flows mimicking various inhalation maneuvers under the influence of aerodynamic (i.e. drag and lift) and gravitational forces. We quantify transport and deposition metrics for such fibers, including residence time and penetration depth, compared with spherical particles of equivalent mass. Our findings underscore how deposition depth is largely independent of AR under oscillatory conditions, in contrast with previous works where AR was found to influence deposition depth under steady inspiratory flow. Overall, our efforts underline the importance of modeling oscillatory breathing when predicting fiber deposition in the distal lungs, as they are inhaled and exhaled during a full inspiratory cycle. Such physical insight helps further explore the potential of fiber particles as attractive carriers for deep airway targeting.

Development of the lung

To fulfill the task of gas exchange, the lung possesses a huge inner surface and a tree-like system of conducting airways ventilating the gas exchange area. During lung development, the conducting airways are formed first, followed by the formation and enlargement of the gas exchange area. The latter (alveolarization) continues until young adulthood. During organogenesis, the left and right lungs have their own anlage, an outpouching of the foregut. Each lung bud starts a repetitive process of outgrowth and branching (branching morphogenesis) that forms all of the future airways mainly during the pseudoglandular stage. During the canalicular stage, the differentiation of the epithelia becomes visible and the bronchioalveolar duct junction is formed. The location of this junction stays constant throughout life. Towards the end of the canalicular stage, the first gas exchange may take place and survival of prematurely born babies becomes possible. Ninety percent of the gas exchange surface area will be formed by alveolarization, a process where existing airspaces are subdivided by the formation of new walls (septa). This process requires a double-layered capillary network at the basis of the newly forming septum. However, in parallel to alveolarization, the double-layered capillary network of the immature septa fuses to a single-layered network resulting in an optimized setup for gas exchange. Alveolarization still continues, because, at sites where new septa are lifting off preexisting mature septa, the required second capillary layer will be formed instantly by angiogenesis. The latter confirms a lifelong ability of alveolarization, which is important for any kind of lung regeneration.

Do nanoparticles provide a new opportunity for diagnosis of distal airspace disease?

There is a need for efficient techniques to assess abnormalities in the peripheral regions of the lungs, for example, for diagnosis of pulmonary emphysema. Considerable scientific efforts have been directed toward measuring lung morphology by studying recovery of inhaled micron-sized aerosol particles (0.4-1.5 µm). In contrast, it is suggested that the recovery of inhaled airborne nanoparticles may be more useful for diagnosis. The objective of this work is to provide a theoretical background for the use of nanoparticles in measuring lung morphology and to assess their applicability based on a review of the literature. Using nanoparticles for studying distal airspace dimensions is shown to have several advantages over other aerosol-based methods. 1) Nanoparticles deposit almost exclusively by diffusion, which allows a simpler breathing maneuver with minor artifacts from particle losses in the oropharyngeal and upper airways. 2) A higher breathing flow rate can be utilized, making it possible to rapidly inhale from residual volume to total lung capacity (TLC), thereby eliminating the need to determine the TLC before measurement. 3) Recent studies indicate better penetration of nanoparticles than micron-sized particles into poorly ventilated and diseased regions of the lungs; thus, a stronger signal from the abnormal parts is expected. 4) Changes in airspace dimensions have a larger impact on the recovery of nanoparticles. Compared to current diagnostic techniques with high specificity for morphometric changes of the lungs, computed tomography and magnetic resonance imaging with hyperpolarized gases, an aerosol-based method is likely to be less time consuming, considerably cheaper, simpler to use, and easier to interpret (providing a single value rather than an image that has to be analyzed). Compared to diagnosis by carbon monoxide (DL,CO), the uptake of nanoparticles in the lung is not affected by blood flow, hemoglobin concentration or alterations of the alveolar membranes, but relies only on lung morphology.

Streamline crossing: An essential mechanism for aerosol dispersion in the pulmonary acinus

The dispersion of inhaled microparticles in the pulmonary acinus of the lungs is often attributed to the complex interplay between convective mixing, due to irreversible flows, and intrinsic particle motion (i.e. gravity and diffusion). However, the role of each mechanism, the exact nature of such interplay between them and their relative importance still remain unclear. To gain insight into these dispersive mechanisms, we track liquid-suspended microparticles and extract their effective diffusivities inside an anatomically-inspired microfluidic acinar model. Such results are then compared to experiments and numerical simulations in a straight channel. While alveoli of the proximal acinar generations exhibit convective mixing characteristics that lead to irreversible particle trajectories, this local effect is overshadowed by a more dominant dispersion mechanism across the ductal branching network that arises from small but significant streamline crossing due to intrinsic diffusional motion in the presence of high velocity gradients. We anticipate that for true airborne particles, which exhibit much higher intrinsic motion, streamline crossing would be even more significant.

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.

Computational Models of Inhalation Therapy in Early Childhood: Therapeutic Aerosols in the Developing Acinus

BACKGROUND

Inhalation therapy targeted to the deep alveolated regions holds great promise, specifically in pediatric populations. Yet, inhalation devices and medical protocols are overwhelmingly derived from adult guidelines, with very low therapeutic efficiency in young children. During the first years of life, airway remodeling and changing ventilation patterns are anticipated to alter aerosol deposition with underachieving outcomes in infants. As past research is still overwhelmingly focused on adults or limited to models of upper airways, a fundamental understanding of inhaled therapeutic transport and deposition in the acinar regions is needed to shed light on delivering medication to the developing alveoli.

METHODS

Using computational fluid dynamics (CFD), we simulated inhalation maneuvers in anatomically-inspired models of developing acinar airways, covering the distinct phases of lung development, from underdeveloped, saccular pulmonary architectures in infants, to structural changes in toddlers, ultimately mimicking space-filling morphologies of a young child, representing scaled-down adult lungs. We model aerosols whose diameters span the range of sizes acknowledged to reach the alveolar regions and examine the coupling between morphological changes, varying ventilation patterns and particle characteristics on deposition outcomes.

RESULTS

Spatial distributions of deposited particles point to noticeable changes in the patterns of aerosol deposition with age, in particular in the youngest age group examined (3 month). Total deposition efficiency, as well as deposition dispersion, vary not only with the phases of lung development but also and critically with aerosol diameter.

CONCLUSIONS

Given the various challenges when prescribing inhalation therapy to a young infant, our findings underline some mechanistic aspects to consider when targeting medication to the developing alveoli. Not only does the intricate coupling between acinar morphology and ventilation patterns need to be considered, but the physical properties (i.e., aerodynamic size) of therapeutic aerosols also closely affect the anticipated success rates of the inhaled medication.

Aerosols in healthy and emphysematous in silico pulmonary acinar rat models

There has been relatively little attention given on predicting particle deposition in the respiratory zone of the diseased lungs despite the high prevalence of chronic obstructive pulmonary disease (COPD). Increased alveolar volume and deterioration of alveolar septum, characteristic of emphysema, may alter the amount and location of particle deposition compared to healthy lungs, which is particularly important for toxic or therapeutic aerosols. In an attempt to shed new light on aerosol transport and deposition in emphysematous lungs, we performed numerical simulations in models of healthy and emphysematous acini motivated by recent experimental lobar-level data in rats (Oakes et al., 2014a). Compared to healthy acinar structures, models of emphysematous subacini were created by removing inter-septal alveolar walls and enhancing the alveolar volume in either a homogeneous or heterogeneous fashion. Flow waveforms and particle properties were implemented to match the experimental data. The occurrence of flow separation and recirculation within alveolar cavities was found in proximal generations of the healthy zones, in contrast to the radial-like airflows observed in the diseased regions. In agreement with experimental data, simulations point to particle deposition concentrations that are more heterogeneously distributed in the diseased models compared with the healthy one. Yet, simulations predicted less deposition in the emphysematous models in contrast to some experimental studies, a likely consequence due to the shallower penetration depths and modified flow topologies in disease compared to health. These spatial-temporal particle transport simulations provide new insight on deposition in the emphysematous acini and shed light on experimental observations.