Development of characteristic upper tracheobronchial airway models for testing pharmaceutical aerosol delivery

Leveraging Numerical Simulation Technology to Advance Drug Preparation: A Comprehensive Review of Application Scenarios and Cases

BACKGROUND/OBJECTIVES

Numerical simulation plays an important role in pharmaceutical preparation recently. Mechanistic models, as a type of numerical model, are widely used in the study of pharmaceutical preparations. Mechanistic models are based on a priori knowledge, i.e., laws of physics, chemistry, and biology. However, due to interdisciplinary reasons, pharmacy researchers have greater difficulties in using computer models.

METHODS

In this paper, we highlight the application scenarios and examples of mechanistic modelling in pharmacy research and provide a reference for drug researchers to get started.

RESULTS

By establishing a suitable model and inputting preparation parameters, researchers can analyze the drug preparation process. Therefore, mechanistic models are effective tools to optimize the preparation parameters and predict potential quality problems of the product. With product quality parameters as the ultimate goal, the experiment design is optimized by mechanistic models. This process emphasizes the concept of quality by design.

CONCLUSIONS

The use of numerical simulation saves experimental cost and time, and speeds up the experimental process. In pharmacy experiments, part of the physical information and the change processes are difficult to obtain, such as the mechanical phenomena during tablet compression and the airflow details in the nasal cavity. Therefore, it is necessary to predict the information and guide the formulation with the help of mechanistic models.

Development of an effective two-equation turbulence modeling approach for simulating aerosol deposition across a range of turbulence levels

Pharmaceutical aerosol systems present a significant challenge to computational fluid dynamics (CFD) modeling based on the need to capture multiple levels of turbulence, frequent transition between laminar and turbulent flows, anisotropic turbulent particle dispersion, and near-wall particle transport phenomena often within geometrically complex systems over multiple time scales. Two-equation turbulence models, such as the k - ω family of approximations, offer a computationally efficient solution approach, but are known to require the use of near-wall (NW) corrections and eddy interaction model (EIM) modifications for accurate predictions of aerosol deposition. The objective of this study was to develop an efficient and effective two-equation turbulence modeling approach that enables accurate predictions of pharmaceutical aerosol deposition across a range of turbulence levels. Key systems considered were the traditional aerosol deposition benchmark cases of a 90-degree bend (R e = 6,000 ) and a vertical straight section of pipe (R e = 10,000 ), as well as a highly complex case of direct-to-infant (D2I) nose-to-lung pharmaceutical aerosol delivery from an air-jet dry powder inhaler (DPI) including a patient interface and infant nasal geometry through mid-trachea (500 < R e < 7,000 ). Of the k - ω family of models, the low Reynolds number (LRN) shear stress transport (SST) approach was determined to provide the best agreement with experimental aerosol deposition data in the D2I system, based on an improved simulation of turbulent jet flow that frequently occurs in DPIs. Considering NW corrections, a new correlation was developed to quantitatively predict best regional values of the y + l i m i t , within which anisotropic NW turbulence is approximated. Considering EIM modifications, a previously described drift correction approach was implemented in pharmaceutical aerosol simulations for the first time. Considering all model corrections and modifications applied to the D2I system, regional relative errors in deposition fractions between CFD predictions and new experimental data were improved from 19-207% (no modifications) to 2-15% (all modifications) with a notable decrease in computational time (up to ∼15%). In conclusion, the highly efficient two-equation k - ω models with physically realistic corrections and modifications provided a viable, efficient and accurate approach to simulate the transport and deposition of pharmaceutical aerosols in complex airway systems that include laminar, turbulent and transitional flows.

Inhalation Profiles Through a Dry Powder Inhaler: Relation Between Inhalation Technique and Spirometric Measures

Background: The understanding of the real flow profiles through a dry powder inhaler (DPI), generated by asthma patients, is a prerequisite for satisfactory drug delivery to the lungs. The aims of the study were to assess the relationship between spirometric measures and inhalation profiles through a low-resistance DPI, and to compare parameters of those profiles between optimal and suboptimal inhalation technique type. Methods: Both healthy adult volunteers and patients with asthma were included in the study. Spirometry was conducted along with modified flow-volume test to detect expiratory levels (maximum "100%" exhalation to residual volume [RV] and halfway "50%" to RV). These were the reference levels of the depth of exhalation for each patient to simulate the effect of incomplete exhalation. Individual inhalation profiles were recorded using spirometry in-house software as the volumetric airflow through the inhaler versus time. Inspiratory flow parameters were extracted: time to peak inspiratory flow through inhaler (PIF_inh), time at which peak inspiratory flow occurs (tPIF_inh), total inhalation time (T), and inhaled volume during maneuver (V). Results and Conclusions: There are significant relationships between spirometric indices and parameters of inhalation through a low-resistance, cyclohaler-type DPI (assessed by single-factor analysis of Spearman's rank correlation coefficient). Multiple regression models were constructed, predicting inspiratory flow parameters (including spirometric indices, demographic parameters, and inhaler's usage history as determinants). The exhalation halfway to RV before inhalation did not affect significantly PIF_inh and tPIF_inh (and, thus, initial flow dynamics) in asthma patients. T and V parameters were then significantly decreased, but seemed sufficient for successful DPI performance. Both exhalation to RV and incomplete exhalation halfway to RV preceding inhalation allow for effective usage of low-resistance DPI.

Influences of puff protocols and upper airway anatomy on cannabis pharmacokinetics: A CFPD-PK study

Predicting the optimal administration doses of the inhaled Δ9-tetrahydrocannabinol (THC), i.e., one of the major natural compounds in cannabis, is critical for maximizing the therapeutic outcomes and minimizing the toxic side effects. Thus, it is essential to developing an aerosol dosimetry model to simulate the transport, deposition, and translocation of inhaled THC aerosols from the human respiratory system to the systemic region. In this study, a computational fluid-particle dynamics (CFPD) plus pharmacokinetics (PK) model was developed and validated to quantify the localized vapor and particle uptake rates of THC and the resultant THC-plasma concentrations using two human upper airway geometries. In addition, two different puff protocols (4.0/10.0 s and 1.6/11.4 s as the inhalation/holding time ratios) were employed, associated with two different inhaled THC doses (2.0 mg and 8.82 mg, respectively). The computational results demonstrated that multiple parameters had noticeable influences on THC particle deposition and vapor absorption in the upper airways, as well as the resultant pharmacokinetic behaviors. These factors include anatomical features of the upper airway, puff flow rate, duration, and holding time. The results indicated that puff protocol with 4.0/10.0 s inhalation/holding time ratio would be recommended if the treatment needs THC delivery to the deeper lung. Furthermore, the inhaled THC dose had a dominant effect on the THC-plasma PK profiles, which could override the influences of anatomical variability and puff protocols. The developed CFPD-PK modeling framework has the potential to provide localized lung absorption data and PK profiles for in vitro-in vivo correlation, as well as supporting the development and assessment of drug products containing cannabis or cannabis-derived compounds.

Role of CFD based in silico modelling in establishing an in vitro-in vivo correlation of aerosol deposition in the respiratory tract

Effective evaluation and prediction of aerosol transport deposition in the human respiratory tracts are critical to aerosol drug delivery and evaluation of inhalation products. Establishment of an in vitro-in vivo correlation (IVIVC) requires the understanding of flow and aerosol behaviour and underlying mechanisms at the microscopic scale. The achievement of the aim can be facilitated via computational fluid dynamics (CFD) based in silico modelling which treats the aerosol delivery as a two-phase flow. CFD modelling research, in particular coupling with discrete phase model (DPM) and discrete element method (DEM) approaches, has been rapidly developed in the past two decades. This paper reviews the recent development in this area. The paper covers the following aspects: geometric models of the respiratory tract, CFD turbulence models for gas phase and its coupling with DPM/DEM for aerosols, and CFD investigation of the effects of key factors associated with geometric variations, flow and powder characteristics. The review showed that in silico study based on CFD models can effectively evaluate and predict aerosol deposition pattern in human respiratory tracts. The review concludes with recommendations on future research to improve in silico prediction to achieve better IVIVC.

Focused targeting of inhaled magnetic aerosols in reconstructed in vitro airway models

The pulmonary tract is an attractive route for topical treatments of lung diseases. Yet, our ability to confine the deposition of inhalation aerosols to specific lung regions, or local airways, remains still widely beyond reach. It has been hypothesized that by coupling magnetic particles to inhaled therapeutics the ability to locally target airway sites can be substantially improved. Although the underlying principle has shown promise in seminal in vivo animal experiments as well as in vitro and in silico studies, its practical implementation has come short of delivering efficient localized airway targeting. Here, we demonstrate in an in vitro proof-of-concept an inhalation framework to leverage magnetically-loaded aerosols for airway targeting in the presence of an external magnetic field. By coupling the delivery of a short pulsed bolus of sub-micron (~500 nm diameter) droplet aerosols with a custom ventilation machine that tracks the volume of air inhaled past the bolus, focused targeting can be maximized during a breath hold maneuver. Specifically, we visualize the motion of the pulsed SPION-laden (superparamagnetic iron oxide nanoparticles) aerosol bolus and quantify under microscopy ensuing deposition patterns in reconstructed 3D airway models. Our aerosol inhalation platform allows for the first time to deposit inhaled particles to specific airway sites while minimizing undesired deposition across the remaining airspace, in an effort to significantly augment the targeting efficiency (i.e. deposition ratio between targeted and untargeted regions). Such inhalation strategy may pave the way for improved treatment outcomes, including reducing side effects in chemotherapy.

Computational fluid dynamic models as tools to predict aerosol distribution in tracheobronchial airways

Aerosol and pollutants, in form of particulates 5-8 μm in main size face every day our respiratory system as natural suspension in air or forced to be inhaled as a coadjutant in a medical therapy for respiratory diseases. This inhalation happens in children to elderly, women and men, healthy or sick and disable people. In this paper we analyzed the inhalation of aerosol in conditions assimilable to the thermal therapy. We use a computational fluid dynamic 3D model to compute and visualize the trajectories of aerosol (3-7-10-25 µm) down to the sixth generation of bronchi, in a steady and dynamic condition (7 µm) set as breath cycle at rest. Results, compared to a set of milestone experimental studies published in literature, allow the comprehension of particles behavior during the inhalation from mouth to bronchi sixth generation, the visualization of jet at larynx constriction and vortices, in an averaged characteristic rigorous geometrical model including tracheal rings. Results on trajectories and deposition show the importance of the including transient physiological breath cycle on aerosol deposition analyses. Numerical and graphical results, may enable the design of medical devices and protocols to make the inhalations more effective in all the users' population.

In silico optimization of targeted aerosol delivery in upper airways via Inhaled Volume Tracking

BACKGROUND

Despite the widespread use of aerosol inhalation as a drug delivery method, targeted delivery to the upper airways remains an ongoing challenge in the quest for improved clinical response in respiratory disease.

METHODS

Here, we examine in silico flow and particle dynamics when using an oral Inhaled Volume Tracking manoeuvre. A short pulsed aerosol bolus is injected during slow inhalation flow rates followed by clean air, and a breath-hold is initiated once it reaches the desired depth. We explore the fate of a broad particle size range (1-40 μm) for both upright and supine positions.

FINDINGS

Our findings illustrate that despite attempts to mitigate dispersion using slower flow rates, the laryngeal jet disperses the aerosol bolus and thus remains a hurdle for efficient targeted delivery. Nevertheless, we show a decrease in extra-thoracic deposition; large aerosols in the range of 10-30 μm potentially outperform existing inhalation methods, showing deposition fractions of up to 80% in an upright orientation.

INTERPRETATION

The improved deposition during Inhaled Volume Tracking shows promise for clinical applications and could be leveraged to deliver larger payloads to the upper airways.

CFD Guided Optimization of Nose-to-Lung Aerosol Delivery in Adults: Effects of Inhalation Waveforms and Synchronized Aerosol Delivery

PURPOSE

The objective of this study was to optimize nose-to-lung aerosol delivery in an adult upper airway model using computational fluid dynamics (CFD) simulations in order to guide subsequent human subject aerosol delivery experiments.

METHODS

A CFD model was developed that included a new high-flow nasal cannula (HFNC) and pharmaceutical aerosol delivery unit, nasal cannula interface, and adult upper airway geometry. Aerosol deposition predictions in the system were validated with existing and new experimental results. The validated CFD model was then used to explore aerosol delivery parameters related to synchronizing aerosol generation with inhalation and inhalation flow rate.

RESULTS

The low volume of the new HFNC unit minimized aerosol transit time (0.2 s) and aerosol bolus spread (0.1 s) enabling effective synchronization of aerosol generation with inhalation. For aerosol delivery correctly synchronized with inhalation, a small particle excipient-enhanced growth delivery strategy reduced nasal cannula and nasal depositional losses each by an order of magnitude and enabled ~80% of the nebulized dose to reach the lungs. Surprisingly, nasal deposition was not sensitive to inhalation flow rate due to use of a nasal cannula interface with co-flow inhaled air and the small initial particle size.

CONCLUSIONS

The combination of correct aerosol synchronization and small particle size enabled high efficiency nose-to-lung aerosol delivery in adults, which was not sensitive to inhalation flow rate.

In Silico Optimization of Fiber-Shaped Aerosols in Inhalation Therapy for Augmented Targeting and Deposition across the Respiratory Tract

Motivated by a desire to uncover new opportunities for designing the size and shape of fiber-shaped aerosols towards improved pulmonary drug delivery deposition outcomes, we explore the transport and deposition characteristics of fibers under physiologically inspired inhalation conditions in silico, mimicking a dry powder inhaler (DPI) maneuver in adult lung models. Here, using computational fluid dynamics (CFD) simulations, we resolve the transient translational and rotational motion of inhaled micron-sized ellipsoid particles under the influence of aerodynamic (i.e., drag, lift) and gravitational forces in a respiratory tract model spanning the first seven bifurcating generations (i.e., from the mouth to upper airways), coupled to a more distal airway model representing nine generations of the mid-bronchial tree. Aerosol deposition efficiencies are quantified as a function of the equivalent diameter (dp) and geometrical aspect ratio (AR), and these are compared to outcomes with traditional spherical particles of equivalent mass. Our results help elucidate how deposition patterns are intimately coupled to dp and AR, whereby high AR fibers in the narrow range of dp = 6-7 µm yield the highest deposition efficiency for targeting the upper- and mid-bronchi, whereas fibers in the range of dp= 4-6 µm are anticipated to cross through the conducting regions and reach the deeper lung regions. Our efforts underscore previously uncovered opportunities to design the shape and size of fiber-like aerosols towards targeted pulmonary drug delivery with increased deposition efficiencies, in particular by leveraging their large payloads for deep lung deposition.

Ultrafine particle deposition in a realistic human airway at multiple inhalation scenarios

The scarcity of regional deposition data in distal respiratory airways represents an important challenge for current toxicology and pharmacology research. To bridge this gap, a realistic airway model extending from nasal and oral openings to distal bronchial airways with varying pathway length was built in this study. Transport and deposition characteristics of naturally inhaled ultrafine particles (UFPs) ranging from 1 to 100 nm were numerically investigated, and effects of different inhalation scenarios were considered. To enable intercase particle deposition comparison, an adjusted parameter, unified deposition enhancement factor (UDEF), was proposed for quantifying the localised deposition concentration. Results show that compartment particle deposition peaked around the ultrafine end of the considered size range, and it dropped rapidly with the increase of particle size. Different inhalation modes caused notable deposition changes in the extrathoracic region, while its effects in the TB airway are much less. For UFPs larger than 10 nm, predicted deposition efficiencies in all compartments are all at lowest levels among considered particle size range, implying UFPs ranging from 10 to 100 nm can travel through the whole respiratory airway model and escape to the alveolar region. Furthermore, high enhancement factors were observed at the vicinity of most bifurcation apexes, and more even UDEF distribution was observed from 1-nm particle cases. While for 100-nm cases, the deposited particles tend to concentrate at few "hot spots" (areas of high deposition concentration in relation to surrounding surfaces) with greater UDEF in the tracheobronchial airway.

In Silico Methods for Development of Generic Drug-Device Combination Orally Inhaled Drug Products

The development of generic, single‐entity, drug–device combination products for orally inhaled drug products is challenging in part because of the complex nature of device design characteristics and the difficulties associated with establishing bioequivalence for a locally acting drug product delivered to the site of action in the lung. This review examines in silico models that may be used to support the development of generic orally inhaled drug products and how model credibility may be assessed.

Innovation for Generic Drugs: Science and Research Under the Generic Drug User Fee Amendments of 2012

Regulatory science is science and research intended to improve decision making in a regulatory framework. Improvements in decision making can be in both accuracy (making better decisions) and in efficiency (making faster decisions). Science and research supported by the Generic Drug User Fee Amendments of 2012 (GDUFA) have focused on two innovative methodologies that work together to enable new approaches to development and review of generic drugs: quantitative models and advanced in vitro product characterization. Quantitative models faithfully represent current scientific understanding. They are tools pharmaceutical scientists and clinical pharmacologists use for making better and faster product development decisions. Advances in the in vitro product comparisons provide the measurements of product differences that are the critical input into the models. This paper outlines four areas where science and research funded by GDUFA support synergistic use of models and characterization at critical decision points during generic drug product development and review.

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.

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.

Development of an infant complete-airway in vitro model for evaluating aerosol deposition

A complete-airway in vitro model would be very useful for toxicological dosimetry testing and for developing targeted inhaled medications in cases where conducting in vivo experiments are exceedingly difficult, as with infants. The objective of this study was to determine whether packed bed in vitro models, which contain spheres as the primary repeating unit, provide a realistic representation of aerosol deposition in the tracheobronchial region of infant lungs based on computational fluid dynamics (CFD) predictions. The packed bed (PB) CFD model contained an inlet consistent with airway bifurcation B3 (∼lobar bronchi) leading to a spherical array with voids between the spheres forming a divided flow pathway. The hydrodynamic diameter of the voids was approximately matched to the diameter of bifurcations in various lung regions. For comparison, a CFD stochastic individual pathway (SIP) geometry with realistic bifurcations extending from B4-B15 (terminal bronchioles) was selected as an anatomically accurate model. The CFD-SIP model predictions were benchmarked with existing algebraic correlations for aerosol deposition in the lungs and found to be reasonable. Unfortunately, the CFD-PB model did not provide a good representation of aerosol deposition in the tracheobronchial region of human lungs. Through careful selection of the PB sphere size and inlet conditions, total deposition in the CFD-PB model matched CFD-SIP deposition within 10% absolute error across a range of relevant aerosol sizes. However, regional deposition within the CFD-PB model was very different from the CFD-SIP case. Therefore, the PB approach cannot be recommended for determining spatial or temporal distribution of aerosol transport and impaction deposition through the lungs.

Recommendations for Simulating Microparticle Deposition at Conditions Similar to the Upper Airways with Two-Equation Turbulence Models

The development of a CFD model, from initial geometry to experimentally validated result with engineering insight, can be a time-consuming process that often requires several iterations of meshing and solver set-up. Applying a set of guidelines in the early stages can help to streamline the process and improve consistency between different models. The objective of this study was to determine both mesh and CFD solution parameters that enable the accurate simulation of microparticle deposition under flow conditions consistent with the upper respiratory airways including turbulent flow. A 90° bend geometry was used as a characteristic model that occurs throughout the airways and for which high-quality experimental aerosol deposition data is available in the transitional and turbulent flow regimes. Four meshes with varying degrees of near-wall resolution were compared, and key solver settings were applied to determine the parameters that minimize sensitivity to the near-wall (NW) mesh. The Low Reynolds number (LRN) k-ω model was used to resolve the turbulence field, which is a numerically efficient two-equation turbulence model, but has recently been considered overly simplistic. Some recent studies have used more complex turbulence models, such as Large Eddy Simulation (LES), to overcome the perceived weaknesses of two-equation models. Therefore, the secondary objective was to determine whether the more computationally efficient LRN k-ω model was capable of providing deposition results that were comparable to LES. Results show how NW mesh sensitivity is reduced through application of the Green-Gauss Node-based gradient discretization scheme and physically realistic near-wall corrections. Using the newly recommended meshing parameters and solution guidelines gives an excellent match to experimental data. Furthermore, deposition data from the LRN k-ω model compares favorably with LES results for the same characteristic geometry. In summary, this study provides a set of meshing and solution guidelines for simulating aerosol deposition in transitional and turbulent flows found in the upper respiratory airways using the numerically efficient LRN k-ω approach.

Aerosol Drug Delivery During Noninvasive Positive Pressure Ventilation: Effects of Intersubject Variability and Excipient Enhanced Growth

BACKGROUND

Nebulized aerosol drug delivery during the administration of noninvasive positive pressure ventilation (NPPV) is commonly implemented. While studies have shown improved patient outcomes for this therapeutic approach, aerosol delivery efficiency is reported to be low with high variability in lung-deposited dose. Excipient enhanced growth (EEG) aerosol delivery is a newly proposed technique that may improve drug delivery efficiency and reduce intersubject aerosol delivery variability when coupled with NPPV.

MATERIALS AND METHODS

A combined approach using in vitro experiments and computational fluid dynamics (CFD) was used to characterize aerosol delivery efficiency during NPPV in two new nasal cavity models that include face mask interfaces. Mesh nebulizer and in-line dry powder inhaler (DPI) sources of conventional and EEG aerosols were both considered.

RESULTS

Based on validated steady-state CFD predictions, EEG aerosol delivery improved lung penetration fraction (PF) values by factors ranging from 1.3 to 6.4 compared with conventional-sized aerosols. Furthermore, intersubject variability in lung PF was very high for conventional aerosol sizes (relative differences between subjects in the range of 54.5%-134.3%) and was reduced by an order of magnitude with the EEG approach (relative differences between subjects in the range of 5.5%-17.4%). Realistic in vitro experiments of cyclic NPPV demonstrated similar trends in lung delivery to those observed with the steady-state simulations, but with lower lung delivery efficiencies. Reaching the lung delivery efficiencies reported with the steady-state simulations of 80%-90% will require synchronization of aerosol administration during inspiration and reducing the size of the EEG aerosol delivery unit.

CONCLUSIONS

The EEG approach enabled high-efficiency lung delivery of aerosols administered during NPPV and reduced intersubject aerosol delivery variability by an order of magnitude. Use of an in-line DPI device that connects to the NPPV mask appears to be a convenient method to rapidly administer an EEG aerosol and synchronize the delivery with inspiration.

Validating Whole-Airway CFD Predictions of DPI Aerosol Deposition at Multiple Flow Rates

BACKGROUND

The objective of this study was to compare aerosol deposition predictions of a new whole-airway CFD model with available in vivo data for a dry powder inhaler (DPI) considered across multiple inhalation waveforms, which affect both the particle size distribution (PSD) and particle deposition.

METHODS

The Novolizer DPI with a budesonide formulation was selected based on the availability of 2D gamma scintigraphy data in humans for three different well-defined inhalation waveforms. Initial in vitro cascade impaction experiments were conducted at multiple constant (square-wave) particle sizing flow rates to characterize PSDs. The whole-airway CFD modeling approach implemented the experimentally determined PSDs at the point of aerosol formation in the inhaler. Complete characteristic airway geometries for an adult were evaluated through the lobar bronchi, followed by stochastic individual pathway (SIP) approximations through the tracheobronchial region and new acinar moving wall models of the alveolar region.

RESULTS

It was determined that the PSD used for each inhalation waveform should be based on a constant particle sizing flow rate equal to the average of the inhalation waveform's peak inspiratory flow rate (PIFR) and mean flow rate [i.e., AVG(PIFR, Mean)]. Using this technique, agreement with the in vivo data was acceptable with <15% relative differences averaged across the three regions considered for all inhalation waveforms. Defining a peripheral to central deposition ratio (P/C) based on alveolar and tracheobronchial compartments, respectively, large flow-rate-dependent differences were observed, which were not evident in the original 2D in vivo data.

CONCLUSIONS

The agreement between the CFD predictions and in vivo data was dependent on accurate initial estimates of the PSD, emphasizing the need for a combination in vitro-in silico approach. Furthermore, use of the AVG(PIFR, Mean) value was identified as a potentially useful method for characterizing a DPI aerosol at a constant flow rate.

Particle transport in the human respiratory tract: formulation of a nodal inverse distance weighted Eulerian-Lagrangian transport and implementation of the Wind-Kessel algorithm for an oral delivery

This paper is the first in a series wherein efficient computational methods are developed and implemented to accurately quantify the transport, deposition, and clearance of the microsized particles (range of interest: 2 to 10 µm) in the human respiratory tract. In particular, this paper (part I) deals with (i) development of a detailed 3D computational finite volume mesh comprising of the NOPL (nasal, oral, pharyngeal and larynx), trachea and several airway generations; (ii) use of CFD Research Corporation's finite volume Computational Biology (CoBi) flow solver to obtain the flow physics for an oral inhalation simulation; (iii) implement a novel and accurate nodal inverse distance weighted Eulerian-Lagrangian formulation to accurately obtain the deposition, and (iv) development of Wind-Kessel boundary condition algorithm. This new Wind-Kessel boundary condition algorithm allows the 'escaped' particles to reenter the airway through the outlets, thereby to an extent accounting for the drawbacks of having a finite number of lung generations in the computational mesh. The deposition rates in the NOPL, trachea, the first and second bifurcation were computed, and they were in reasonable accord with the Typical Path Length model. The quantitatively validated results indicate that these developments will be useful for (i) obtaining depositions in diseased lungs (because of asthma and COPD), for which there are no empirical models, and (ii) obtaining the secondary clearance (mucociliary clearance) of the deposited particles. Copyright © 2015 John Wiley & Sons, Ltd.

Current Inhalers Deliver Very Small Doses to the Lower Tracheobronchial Airways: Assessment of Healthy and Constricted Lungs

To evaluate the regional delivery of conventional aerosol medications, a new whole-lung computational fluid dynamics modeling approach was applied for metered dose inhaler (MDI) and dry powder inhaler (DPI) aerosols delivered to healthy and constricted airways. The computational fluid dynamics approach included complete airways through the third respiratory bifurcation (B3) and applied the new stochastic individual pathway modeling technique beyond B3 through the remainder of the conducting airways together with a new model of deposition in the alveolar region. Bronchiolar (B8-B15) deposition fraction values were low (∼1%) for both MDI and DPI aerosols with the healthy geometry, whereas delivery to the constricted model was even lower, with deposition fraction values of 0.89% and 0.81% for the MDI and DPI, respectively. Calculating dose per unit surface area for the commercial MDI and DPI products resulted in approximately 10(-3) μg/cm(2) in the lower tracheobronchial region of B8-B15 and 10(-4) μg/cm(2) in the alveolar region. Across the lung, dose per unit surface area varied by 2 orders of magnitude, which increased to 4 orders of magnitude when the mouth-throat region was included. The MDI and DPI both provided very low drug dose per unit surface area to the small tracheobronchial and alveolar airways.

Validating CFD Predictions of Pharmaceutical Aerosol Deposition with In Vivo Data

PURPOSE

CFD provides a powerful approach to evaluate the deposition of pharmaceutical aerosols; however, previous studies have not compared CFD results of deposition throughout the lungs with in vivo data.

METHODS

The in vivo datasets selected for comparison with CFD predictions included fast and slow clearance of monodisperse aerosols as well as 2D gamma scintigraphy measurements for a dry powder inhaler (DPI) and softmist inhaler (SMI). The CFD model included the inhaler, a characteristic model of the mouth-throat (MT) and upper tracheobronchial (TB) airways, stochastic individual pathways (SIPs) representing the remaining TB region, and recent CFD-based correlations to predict pharmaceutical aerosol deposition in the alveolar airways.

RESULTS

For the monodisperse aerosol, CFD predictions of total lung deposition agreed with in vivo data providing a percent relative error of 6% averaged across aerosol sizes of 1-7 μm. With the DPI and SMI, deposition was evaluated in the MT, central airways (bifurcations B1-B7), and intermediate plus peripheral airways (B8 through alveoli). Across these regions, CFD predictions produced an average relative error <10% for each inhaler.

CONCLUSIONS

CFD simulations with the SIP modeling approach were shown to accurately predict regional deposition throughout the lungs for multiple aerosol types and different in vivo assessment methods.

Efficient Nose-to-Lung (N2L) Aerosol Delivery with a Dry Powder Inhaler

PURPOSE

Delivering aerosols to the lungs through the nasal route has a number of advantages, but its use has been limited by high depositional loss in the extrathoracic airways. The objective of this study was to evaluate the nose-to-lung (N2L) delivery of excipient enhanced growth (EEG) formulation aerosols generated with a new inline dry powder inhaler (DPI). The device was also adapted to enable aerosol delivery to a patient simultaneously receiving respiratory support from high flow nasal cannula (HFNC) therapy.

METHODS

The inhaler delivered the antibiotic ciprofloxacin, which was formulated as submicrometer combination particles containing a hygroscopic excipient prepared by spray-drying. Nose-to-lung delivery was assessed using in vitro and computational fluid dynamics (CFD) methods in an airway model that continued through the upper tracheobronchial region.

RESULTS

The best performing device contained a 2.3 mm flow control orifice and a 3D rod array with a 3-4-3 rod pattern. Based on in vitro experiments, the emitted dose from the streamlined nasal cannula had a fine particle fraction <5 μm of 95.9% and mass median aerodynamic diameter of 1.4 μm, which was considered ideal for nose-to-lung EEG delivery. With the 2.3-343 device, condensational growth in the airways increased the aerosol size to 2.5-2.7 μm and extrathoracic deposition was <10%. CFD results closely matched the in vitro experiments and predicted that nasal deposition was <2%.

CONCLUSIONS

The developed DPI produced high efficiency aerosolization with significant size increase of the aerosol within the airways that can be used to enable nose-to-lung delivery and aerosol administration during HFNC therapy.

Development of a new technique for the efficient delivery of aerosolized medications to infants on mechanical ventilation

PURPOSE

To evaluate the efficiency of a new technique for delivering aerosols to intubated infants that employs a new Y-connector, access port administration of a dry powder, and excipient enhanced growth (EEG) formulation particles that change size in the airways.

METHODS

A previously developed CFD model combined with algebraic correlations were used to predict delivery system and lung deposition of typical nebulized droplets (MMAD = 4.9 μm) and EEG dry powder aerosols. The delivery system consisted of a Y-connector [commercial (CM); streamlined (SL); or streamlined with access port (SL-port)] attached to a 4-mm diameter endotracheal tube leading to the airways of a 6-month-old infant.

RESULTS

Compared to the CM device and nebulized aerosol, the EEG approach with an initial 0.9 μm aerosol combined with the SL and SL-port geometries reduced device depositional losses by factors of 3-fold and >10-fold, respectively. With EEG powder aerosols, the SL geometry provided the maximum tracheobronchial deposition fraction (55.7%), whereas the SL-port geometry provided the maximum alveolar (67.6%) and total lung (95.7%) deposition fractions, respectively.

CONCLUSIONS

Provided the aerosol can be administered in the first portion of the inspiration cycle, the proposed new method can significantly improve the deposition of pharmaceutical aerosols in the lungs of intubated infants.

Creation of an in vitro biomechanical model of the trachea using rapid prototyping

Previous in vitro models of the airways are either rigid or, if flexible, have not matched in vivo compliance characteristics. Rapid prototyping provides a quickly evolving approach that can be used to directly produce in vitro airway models using either rigid or flexible polymers. The objective of this study was to use rapid prototyping to directly produce a flexible hollow model that matches the biomechanical compliance of the trachea. The airway model consisted of a previously developed characteristic mouth-throat region, the trachea, and a portion of the main bronchi. Compliance of the tracheal region was known from a previous in vivo imaging study that reported cross-sectional areas over a range of internal pressures. The compliance of the tracheal region was matched to the in vivo data for a specific flexible resin by iteratively selecting the thicknesses and other dimensions of tracheal wall components. Seven iterative models were produced and illustrated highly non-linear expansion consisting of initial rapid size increase, a transition region, and continued slower size increase as pressure was increased. Thickness of the esophageal interface membrane and initial trachea indention were identified as key parameters with the final model correctly predicting all phases of expansion within a value of 5% of the in vivo data. Applications of the current biomechanical model are related to endotracheal intubation and include determination of effective mucus suctioning and evaluation of cuff sealing with respect to gases and secretions.

Targeted Lung Delivery of Nasally Administered Aerosols

Using the nasal route to deliver pharmaceutical aerosols to the lungs has a number of advantages including co-administration during non-invasive ventilation. The objective of this study was to evaluate the growth and deposition characteristics of nasally administered aerosol throughout the conducting airways based on delivery with streamlined interfaces implementing two forms of controlled condensational growth technology. Characteristic conducting airways were considered including a nose-mouth-throat (NMT) geometry, complete upper tracheobronchial (TB) model through the third bifurcation (B3), and stochastic individual path (SIP) model to the terminal bronchioles (B15). Previously developed streamlined nasal cannula interfaces were used for the delivery of submicrometer particles using either enhanced condensational growth (ECG) or excipient enhanced growth (EEG) techniques. Computational fluid dynamics (CFD) simulations predicted aerosol transport, growth and deposition for a control (4.7 μm) and three submicrometer condensational aerosols with budesonide as a model insoluble drug. Depositional losses with condensational aerosols in the cannula and NMT were less than 5% of the initial dose, which represents an order-of-magnitude reduction compared to the control. The condensational growth techniques increased the TB dose by a factor of 1.1-2.6x, delivered at least 70% of the dose to the alveolar region, and produced final aerosol sizes ≥2.5 μm. Compared to multiple commercial orally inhaled products, the nose-to-lung delivery approach increased dose to the biologically important lower TB region by factors as large as 35x. In conclusion, nose-to-lung delivery with streamlined nasal cannulas and condensational aerosols was highly efficient and targeted deposition to the lower TB and alveolar regions.