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

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.

Augmenting regional and targeted delivery in the pulmonary acinus using magnetic particles.

BACKGROUND: It has been hypothesized that by coupling magnetic particles to inhaled therapeutics, the ability to target specific lung regions (eg, only acinar deposition), or even more so specific points in the lung (eg, tumor targeting), can be substantially improved. Although this method has been proven feasible in seminal in vivo studies, there is still a wide gap in our basic understanding of the transport phenomena of magnetic particles in the pulmonary acinar regions of the lungs, including particle dynamics and deposition characteristics. METHODS: Here, we present computational fluid dynamics-discrete element method simulations of magnetically loaded microdroplet carriers in an anatomically inspired, space-filling, multi-generation acinar airway tree. Breathing motion is modeled by kinematic sinusoidal displacements of the acinar walls, during which droplets are inhaled and exhaled. Particle dynamics are governed by viscous drag, gravity, and Brownian motion as well as the external magnetic force. In particular, we examined the roles of droplet diameter and volume fraction of magnetic material within the droplets under two different breathing maneuvers. RESULTS AND DISCUSSION: Our results indicate that by using magnetic-loaded droplets, 100% of the particles that enter are deposited in the acinar region. This is consistent across all particle sizes investigated (ie, 0.5-3.0 µm). This is best achieved through a deep inhalation maneuver combined with a breath-hold. Particles are found to penetrate deep into the acinus and disperse well, while the required amount of magnetic material is maintained low (<2.5%). Although particles in the size range of ~90-500 nm typically show the lowest deposition fractions, our results suggest that this feature could be leveraged to augment targeted delivery.

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.

Unsteady diffusional screening in 3D pulmonary acinar structures: from infancy to adulthood.

Diffusional screening in the lungs is a physical phenomenon where the specific topological arrangement of alveolated airways of the respiratory region leads to a depletion, or 'screening', of oxygen molecules with increasing acinar generation. Here, we revisit diffusional screening phenomena in anatomically-inspired pulmonary acinar models under realistic breathing maneuvers. By modelling 3D bifurcating alveolated airways capturing both convection and diffusion, unsteady oxygen transport is investigated under cyclic breathing motion. To evaluate screening characteristics in the developing lungs during growth, four representative stages of lung development were chosen (i.e. 3 months, 1 year and 9 months, 3 years and adulthood) that capture distinct morphological acinar changes spanning alveolarization phases to isotropic alveolar growth. Numerical simulations unveil the dramatic changes in O2 transport occurring during lung development, where young infants exhibit highest acinar efficiencies that rapidly converge with age to predictions at adulthood. With increased ventilatory effort, transient dynamics of oxygen transport is fundamentally altered compared to tidal breathing and emphasizes the augmented role of convection. Resolving the complex convective acinar flow patterns in 3D acinar trees allows for the first time a spatially-localized and time-resolved characterization of oxygen transport in the pulmonary acinus, from infancy to adulthood.

Particle dynamics and deposition in true-scale pulmonary acinar models.

Particle transport phenomena in the deep alveolated airways of the lungs (i.e. pulmonary acinus) govern deposition outcomes following inhalation of hazardous or pharmaceutical aerosols. Yet, there is still a dearth of experimental tools for resolving acinar particle dynamics and validating numerical simulations. Here, we present a true-scale experimental model of acinar structures consisting of bifurcating alveolated ducts that capture breathing-like wall motion and ensuing respiratory acinar flows. We study experimentally captured trajectories of inhaled polydispersed smoke particles (0.2 to 1 µm in diameter), demonstrating how intrinsic particle motion, i.e. gravity and diffusion, is crucial in determining dispersion and deposition of aerosols through a streamline crossing mechanism, a phenomenon paramount during flow reversal and locally within alveolar cavities. A simple conceptual framework is constructed for predicting the fate of inhaled particles near an alveolus by identifying capture and escape zones and considering how streamline crossing may shift particles between them. In addition, we examine the effect of particle size on detailed deposition patterns of monodispersed microspheres between 0.1-2 µm. Our experiments underline local modifications in the deposition patterns due to gravity for particles ≥0.5 µm compared to smaller particles, and show good agreement with corresponding numerical simulations.

Revisiting pulmonary acinar particle transport: convection, sedimentation, diffusion, and their interplay.

It is largely acknowledged that inhaled particles ranging from 0.001 to 10 m are able to reach and deposit in the alveolated regions of the lungs. To date, however, the bulk of numerical studies have focused mainly on micrometer sized particles whose transport kinematics are governed by convection and sedimentation, thereby capturing only a small fraction of the wider range of aerosols leading to acinar deposition. Too little is still known about the local acinar transport dynamics of inhaled (ultra)fine particles affected by diffusion and convection. Our study aims to fill this gap by numerically simulating the transport characteristics of particle sizes spanning three orders of magnitude (0.01-5 m) covering diffusive, convective, and gravitational aerosol motion across a multigenerational acinar network. By characterizing the deposition patterns as a function of particle size, we find that submicrometer particles [formulae see text (0.1 m)] reach deep into the acinar structure and are prone to deposit near alveolar openings; meanwhile, other particle sizes are restricted to accessing alveolar cavities in proximal generations. Our findings underline that a precise understanding of acinar aerosol transport, and ultrafine particles in particular, is contingent upon resolving the complex convective-diffusive interplay in determining their irreversible kinematics and local deposition sites.

Role of alveolar topology on acinar flows and convective mixing.

Due to experimental challenges, computational simulations are often sought to quantify inhaled aerosol transport in the pulmonary acinus. Commonly, these are performed using generic alveolar topologies, including spheres, toroids, and polyhedra, to mimic the complex acinar morphology. Yet, local acinar flows and ensuing particle transport are anticipated to be influenced by the specific morphological structures. We have assessed a range of acinar models under self-similar breathing conditions with respect to alveolar flow patterns, convective flow mixing, and deposition of fine particles (1.3 µm diameter). By tracking passive tracers over cumulative breathing cycles, we find that irreversible flow mixing correlates with the location and strength of the recirculating vortex inside the cavity. Such effects are strongest in proximal acinar generations where the ratio of alveolar to ductal flow rates is low and interalveolar disparities are most apparent. Our results for multi-alveolated acinar ducts highlight that fine 1 µm inhaled particles subject to alveolar flows are sensitive to the alveolar topology, underlining interalveolar disparities in particle deposition patterns. Despite the simplicity of the acinar models investigated, our findings suggest that alveolar topologies influence more significantly local flow patterns and deposition sites of fine particles for upper generations emphasizing the importance of the selected acinar model. In distal acinar generations, however, the alveolar geometry primarily needs to mimic the space-filling alveolar arrangement dictated by lung morphology.