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

Targeted delivery of inhalable drug particles in the tracheobronchial tree model of a pediatric patient with bronchopneumonia: A numerical study

Pneumonia is a common cause of hospitalization and death in children worldwide. Inhalation therapy is one of the methods treating pneumonia However, there are limited studies that distinguish between the physiology of children and adults, especially with respect to targeted drug delivery. A tracheobronchial (TB) tree model of an 11-year-old child with bronchopneumonia is selected as a testbed for in silico trials of targeted drug delivery. The airflow and particle transport are solved by the computational fluid dynamics method at an airflow rate of 15 LPM. The results indicate that the distribution of deposited particles shows aggregation on the particle release map. Point-source aerosol release (PSAR) method can significantly reduce the deposition efficiency (DE) of particles in the TB tree model. Specifically, the PSAR method can reduce the DE of large particles (i.e., 7.5 µm and 10 µm) by 7.57% and 9.61%, respectively. This enables rapid design of patient-specific treatment for different population age groups and different airway diseases.

Spatial aerosol deposition correlated to anatomic feature development in 6-year-old upper airway computational models

The upper airways of children undergo developmental changes around age 6, yielding differences between adult and pediatric anatomies. These differences include the cricoid ring area shape, the location of narrowest constriction, and the angle of the epiglottis, all of which are expected to alter local fluid dynamic profiles and subsequent upper airway deposition and downstream aerosol delivery of inhaled therapeutics. In this work, we quantify "pediatric"-like and "adult"-like geometric and fluid dynamic features of two computed tomography (CT)-scan derived models of 6-year-old upper airways in healthy subjects and compare to an idealized model. The two CT-scan models had a mixture of "adult"- and "pediatric"-like anatomic features, with Subject B exhibiting more "pediatric"-like features than Subject A, while the idealized model exhibited entirely "adult"-like features. By computational fluid-particle dynamics, these differences in anatomical features yielded distinct local fluid profiles with altered aerosol deposition between models. Notably, the idealized model better predicted deposition characteristics of Subject A, the more "adult"-like model, including the relationship between the impaction parameter, dp2Q and the fraction of deposition across a range of flow rates and particle diameters, as well as deposition of an approximate pharmaceutical particle size distribution model. Our results with even this limited dataset suggest that there are key personalized metrics that are influenced by anatomical development, which should be considered when developing pediatric inhalable therapeutics. Quantifying anatomical development and correlating to aerosol deposition has the potential for high-throughput developmental characterization and informing desired aerosol characteristics for pediatric applications.

Heat Wave and Bushfire Meteorology in New South Wales, Australia: Air Quality and Health Impacts

The depletion of air quality is a major problem that is faced around the globe. In Australia, the pollutants emitted by bushfires play an important role in making the air polluted. These pollutants in the air result in many adverse impacts on the environment. This paper analysed the air pollution from the bushfires from November 2019 to July 2020 and identified how it affects the human respiratory system. The bush fires burnt over 13 million hectares, destroying over 2400 buildings. While these immediate effects were devastating, the long-term effects were just as devastating, with air pollution causing thousands of people to be admitted to hospitals and emergency departments because of respiratory complications. The pollutant that caused most of the health effects throughout Australia was Particulate Matter (PM) PM_2.5 and PM₁₀. Data collection and analysis were covered in this paper to illustrate where and when PM_2.5 and PM_10, and other pollutants were at their most concerning levels. Susceptible areas were identified by analysing environmental factors such as temperature and wind speed. The study identified how these pollutants in the air vary from region to region in the same time interval. This study also focused on how these pollutant distributions vary according to the temperature, which helps to determine the relationship between the heatwave and air quality. A computational model for PM_2.5 aerosol transport to the realistic airways was also developed to understand the bushfire exhaust aerosol transport and deposition in airways. This study would improve the knowledge of the heat wave and bushfire meteorology and corresponding respiratory health impacts.

Polydisperse Aerosol Transport and Deposition in Upper Airways of Age-Specific Lung

A comprehensive understanding of airflow characteristics and particle transport in the human lung can be useful in modelling to inform clinical diagnosis, treatment, and management, including prescription medication and risk assessment for rehabilitation. One of the difficulties in clinical treatment of lung disorders lies in the patients’ variable physical lung characteristics caused by age, amongst other factors, such as different lung sizes. A precise understanding of the comparison between different age groups with various flow rates is missing in the literature, and this study aims to analyse the airflow and aerosol transport within the age-specific lung. ANSYS Fluent solver and the large-eddy simulation (LES) model were employed for the numerical simulation. The numerical model was validated with the available literature and the computational results showed airway size-reduction significantly affected airflow and particle transport in the upper airways. This study reports higher deposition at the mouth-throat region for larger diameter particles. The overall deposition efficiency (DE) increased with airway size reduction and flow rate. Lung aging effected the pressure distribution and a higher pressure drop was reported for the aged lung as compared to the younger lung. These findings could inform medical management through individualised simulation of drug-aerosol delivery processes for the patient-specific lung.

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.

Phlorotannins ameliorate extracellular matrix production in human vocal fold fibroblasts and prevent vocal fold fibrosis via aerosol inhalation in a laser-induced fibrosis model

Vocal fold fibrosis is an abnormal condition characterized by unfavorable changes in the organization of the extracellular matrix in vocal fold lamina propria. To prevent and treat vocal fold fibrosis, a number of synthetic drugs, such as mitomycin C and the glucocorticoid family, are used after surgery, but these are known to have some side effects. Therefore, using both in vitro and in vivo studies, this study investigated whether phlorotannins extracted from Ecklonia cava have the potential to prevent vocal fold fibrosis with minimal side effects. The results show that phlorotannins suppressed both the expression of the fibrotic phenotypic marker and cell migration by inhibiting the activation of the mitogen-activated protein kinase (MAPK) and Smad2/3 signaling pathways in human vocal fold fibroblasts stimulated by transforming growth factor-β. Additionally, phlorotannins exhibited antifibrotic efficacy without an excessive inflammatory response in a laser-induced fibrosis rabbit model when delivered as an aerosol via inhalation. Based on these results, phlorotannins should be considered a promising candidate for use in the prevention of vocal fold fibrosis.

Targeting inhaled aerosol delivery to upper airways in children: Insight from computational fluid dynamics (CFD)

Despite the prevalence of inhalation therapy in the treatment of pediatric respiratory disorders, most prominently asthma, the fraction of inhaled drugs reaching the lungs for maximal efficacy remains adversely low. By and large drug delivery devices and their inhalation guidelines are typically derived from adult studies with child dosages adapted according to body weight. While it has long been recognized that physiological (e.g. airway sizes, breathing maneuvers) and physical transport (e.g. aerosol dynamics) characteristics are critical in governing deposition outcomes, such knowledge has yet to be extensively adapted to younger populations. Motivated by such shortcomings, the present work leverages in a first step in silico computational fluid dynamics (CFD) to explore opportunities for augmenting aerosol deposition in children based on respiratory physiological and physical transport determinants. Using an idealized, anatomically-faithful upper airway geometry, airflow and aerosol motion are simulated as a function of age, spanning a five year old to an adult. Breathing conditions mimic realistic age-specific inhalation maneuvers representative of Dry Powder Inhalers (DPI) and nebulizer inhalation. Our findings point to the existence of a single dimensionless curve governing deposition in the conductive airways via the dimensionless Stokes number (Stk). Most significantly, we uncover the existence of a distinct deposition peak irrespective of age. For the DPI simulations, this peak (∼ 80%) occurs at Stk ≈ 0.06 whereas for nebulizer simulations, the corresponding peak (∼ 45%) occurs in the range of Stk between 0.03-0.04. Such dimensionless findings hence translate to an optimal window of micron-sized aerosols that evolves with age and varies with inhalation device. The existence of such deposition optima advocates revisiting design guidelines for optimizing deposition outcomes in pediatric inhalation therapy.

Deposition of bolus and continuously inhaled aerosols in rhythmically moving terminal alveoli

The particle dynamics in an oscillating alveolus under tidal breathing can be dramatically different from those in a static alveolus. Despite its close relevance to pulmonary drug delivery and health risk from airborne exposure, quantifications of alveolar deposition are scarce due to its inaccessibility to in vivo measurement instruments, tiny size to replicate in vitro, and dynamic wall motions to model. The objective of this study is to introduce a numerical method to quantify alveolar deposition with continuous particle release in a rhythmically oscillating alveolus by integrating the deposition curves for bolus aerosols and use this method to develop correlations applicable in assessing alveolar drug delivery efficiency or dosimetry of inhaled toxicants. An idealized blind-end terminal alveolus model was developed with rhythmically moving alveolar boundary conditions in phase with tidal breathing. The dynamic wall expansion mode and magnitude were based on experimentally measured chest wall motions and tidal volumes. A well-validated Lagrangian tracking model was used to simulate the transport and deposition of inhaled micrometer particles. Large differences were observed between dynamic and static alveoli in particle motion, deposition onset, and final alveolar deposition fraction. Alveolar deposition of bolus aerosols is highly sensitive to breath-holding duration, particle release time, and alveolar dimension. For 1 µm particles, there exists a cut-off release time (zero bolus deposition), which decreases with alveolar size (i.e., 1.0 s in a 0.2-mm-diameter alveolus and 0.56 s in a 0.8-mm-diameter alveolus). The cumulative alveolar deposition was predicted to be 39% for a 0.2-mm-diameter alveolus, 22% for a 0.4-mm-diameter alveolus, and 10% for a 0.8-mm-diameter alveolus. A cumulative alveolar deposition correlation was developed for inhalation delivery with a prescribed period of drug release and the second correlation for the time variation of alveolar deposition of ambient aerosols, both of which captured the relative dependence of the particle release time and alveolar dimension.

Biomimetics of the pulmonary environment in vitro: A microfluidics perspective

The entire luminal surface of the lungs is populated with a complex yet confluent, uninterrupted airway epithelium in conjunction with an extracellular liquid lining layer that creates the air-liquid interface (ALI), a critical feature of healthy lungs. Motivated by lung disease modelling, cytotoxicity studies, and drug delivery assessments amongst other, in vitro setups have been traditionally conducted using macroscopic cultures of isolated airway cells under submerged conditions or instead using transwell inserts with permeable membranes to model the ALI architecture. Yet, such strategies continue to fall short of delivering a sufficiently realistic physiological in vitro airway environment that cohesively integrates at true-scale three essential pillars: morphological constraints (i.e., airway anatomy), physiological conditions (e.g., respiratory airflows), and biological functionality (e.g., cellular makeup). With the advent of microfluidic lung-on-chips, there have been tremendous efforts towards designing biomimetic airway models of the epithelial barrier, including the ALI, and leveraging such in vitro scaffolds as a gateway for pulmonary disease modelling and drug screening assays. Here, we review in vitro platforms mimicking the pulmonary environment and identify ongoing challenges in reconstituting accurate biological airway barriers that still widely prevent microfluidic systems from delivering mainstream assays for the end-user, as compared to macroscale in vitro cell cultures. We further discuss existing hurdles in scaling up current lung-on-chip designs, from single airway models to more physiologically realistic airway environments that are anticipated to deliver increasingly meaningful whole-organ functions, with an outlook on translational and precision medicine.

Particle deposition in tracheobronchial airways of an infant, child and adult

BACKGROUND

Particle deposition in human airways is important for assessing both health effects of inhaled particles and therapeutic efficacy of inhaled drug aerosols, but is not well understood for infants and children.

OBJECTIVE

We investigate particle deposition in infants and children by using computational fluid dynamics (CFD), and compare this with particle deposition in adults.

METHODS

We chose three population age groups: 7-month infant, 4-year old child, and 20-year old adult. Both airway structures and breathing conditions are considered to vary as a human grows from infancy to adulthood. We investigated deposition of micron-size particles (1-10μm) in both the upper (G3-G6) and lower (G9-G12) tracheobronchial (TB) airways under sedentary conditions.

RESULTS

We found that particle deposition in both upper and lower airways is the highest in an infant, next in a child, and lowest in an adult. As age increases, particle deposition decreases in the upper airways but increases in the lower. For infants, inertial impaction is the dominant deposition mechanism, thus particles are deposited more in the upper airways than in the lower. However, particles are deposited more in the lower airways than in the upper in adults, as gravitational sedimentation is the dominant deposition mechanism.

CONCLUSION

Given the differences in the airway structure and particle deposition mechanisms, particle deposition in infants and children differs from that in adults, not only in the efficiency of deposition but also in the site. Our findings provide evidence that "children are not small adults".

Airflow Simulations in Infant, Child, and Adult Pulmonary Conducting Airways

The airway structure continuously evolves from birth to adulthood, influencing airflow dynamics and respiratory mechanics. We currently know very little about how airflow patterns change throughout early life and its impact on airway resistance, namely because of experimental limitations. To uncover differences in respiratory dynamics between age groups, we performed subject-specific airflow simulations in an infant, child, and adult conducting airways. Airflow throughout the respiration cycle was calculated by coupling image-based models of the conducting airways to the global respiratory mechanics, where flow was driven by a pressure differential. Trachea diameter was 19, 9, and 4.5 mm for the adult (36 years, female), child (6 years, male), and infant (0.25 years, female), respectively. Mean Reynolds number within the trachea was nearly the same for each subject (1100) and Womersley number was above unity for all three subjects and largest for the adult, highlighting the significance of transient effects. In general, air speeds and airway resistances within the conducting airways were inversely correlated with age; the 3D pressure drop was highest in the infant model. These simulations provide new insight into age-dependent flow dynamics throughout the respiration cycle within subject-specific airways.

Recent advances in our understanding of the mechanisms of late lung development and bronchopulmonary dysplasia

The objective of lung development is to generate an organ of gas exchange that provides both a thin gas diffusion barrier and a large gas diffusion surface area, which concomitantly generates a steep gas diffusion concentration gradient. As such, the lung is perfectly structured to undertake the function of gas exchange: a large number of small alveoli provide extensive surface area within the limited volume of the lung, and a delicate alveolo-capillary barrier brings circulating blood into close proximity to the inspired air. Efficient movement of inspired air and circulating blood through the conducting airways and conducting vessels, respectively, generates steep oxygen and carbon dioxide concentration gradients across the alveolo-capillary barrier, providing ideal conditions for effective diffusion of both gases during breathing. The development of the gas exchange apparatus of the lung occurs during the second phase of lung development-namely, late lung development-which includes the canalicular, saccular, and alveolar stages of lung development. It is during these stages of lung development that preterm-born infants are delivered, when the lung is not yet competent for effective gas exchange. These infants may develop bronchopulmonary dysplasia (BPD), a syndrome complicated by disturbances to the development of the alveoli and the pulmonary vasculature. It is the objective of this review to update the reader about recent developments that further our understanding of the mechanisms of lung alveolarization and vascularization and the pathogenesis of BPD and other neonatal lung diseases that feature lung hypoplasia.

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.