Numerical optimization of targeted delivery of charged nanoparticles to the ostiomeatal complex for treatment of rhinosinusitis

Characterizing regional drug delivery within the nasal airways

INTRODUCTION

The nose has been receiving increased attention as a route for drug delivery. As the site of deposition constitutes the first point of contact of the body with the drug, characterization of the regional deposition of intranasally delivered droplets or particles is paramount to formulation and device design of new products.

AREAS COVERED

This review article summarizes the recent literature on intranasal regional drug deposition evaluated in vivo, in vitro and in silico, with the aim of correlating parameters measured in vitro with formulation and device performance. We also highlight the relevance of regional deposition to two emerging applications: nose-to-brain drug delivery and intranasal vaccines.

EXPERT OPINION

As in vivo studies of deposition can be costly and time-consuming, researchers have often turned to predictive in vitro and in silico models. Variability in deposition is high due in part to individual differences in nasal geometry, and a complete predictive model of deposition based on spray characteristics remains elusive. Carefully selected or idealized geometries capturing population average deposition can be useful surrogates to in vivo measurements. Continued development of in vitro and in silico models may pave the way for development of less variable and more effective intranasal drug products.

Optimized gravity-driven intranasal drop administration delivers significant doses to the ostiomeatal complex and maxillary sinus

Chronic and allergic rhinosinusitis impacts approximately 12% of the global population. Challenges in rhinosinusitis treatment include paranasal sinus inaccessibility and variability in delivery efficiency among individuals. This study addresses these challenges of drug delivery by developing a high-efficiency, low-variability protocol for nasal drop delivery to the ostiomeatal complex (OMC) and maxillary sinus. Patient-specific nasal casts were dissected to reveal the configurations of conchae and meatus, providing insights into anatomical features amendable for sinus delivery. Fluorescent dye-enhanced videos visualized the dynamic liquid translocation in transparent nasal casts, allowing real-time assessment and quick adjustment to delivery parameters. Dosimetry to the OMC and maxillary sinus were quantified as drop count and mass using a precision scale. Key delivery factors, including the device type, formulation, and head-chin orientation, were systematically investigated in a cohort of ten nasal casts. Results show that both the squeeze bottle and soft-mist nasal pump yielded notably low doses to the OMC with high variability, and no dose from these two devices was detected within the maxillary sinuses. In contrast, the proposed approach, which included a curved nozzle surpassing the nasal valve and leveraged gravity-driven liquid translocation along the lateral nasal wall, delivered significant doses to the OMC and maxillary sinus. Iterative experimentations identified the optimal head tilt to be 40° and chin tilt to be° from the lateral recumbent position. Statistical analyses established the drop count required for effective OMC/sinus delivery. The proposed delivery protocol holds the potential to enhance chronic rhinosinusitis treatment outcomes with low variability. The dual role of nasal anatomy in posing challenges and offering opportunities highlights the need for future investigations using diverse formulations in a larger cohort of nasal models. Optimized gravity-driven intranasal drop administration delivers significant doses to the ostiomeatal complex and maxillary sinus.

Respiratory Fluid Mechanics of the Effect of Mouth Breathing on High-Arched Palate: Computational Fluid Dynamics Analyses

Computational fluid dynamics (CFD) was introduced into the study of palate growth and development to explain the mechanisms by which mouth breathing affects palate descent from an aerodynamic perspective. Cone beam computed tomography (CBCT) data were used to reconstruct a 3-dimensional model during natural mouth breathing of a volunteer. The model was imported into CFX 19.0 for numerical simulation of nasal breathing, mouth-nasal breathing, and mouth breathing. The pressure in the oronasal cavity was analyzed, and the pressure difference between the oral and nasal surfaces of hard palate under different breathing patterns was calculated. CFD can be used to simulate the stress on the oral and nasal surfaces of the palate under different breathing patterns. The pressure differences and resultant force between the oral and nasal surfaces of the hard palate during nasal inspiration, nasal expiration, mouth-nasal inspiration, mouth-nasal expiration, mouth inspiration, and mouth expiration were 0 Pa, 4 Pa (upward), 9 Pa (upward), 3 Pa (downward), 474 Pa (upward), 263 Pa (downward), respectively, and 87.99 N (upward), 88.03 N (upward), 88.01 N (upward), 88.01 N (upward), 88.05 N (upward), 87.94 N (upward), respectively. Therefore, CFD can be used to investigate the growth and development of the palate. When the volunteer opened his mouth, the pressure difference between the oral and nasal surfaces of the hard palate was about 88 N upward regardless of whether there was airflow in the mouth. The reversal of the direction of the force on the hard palate may be one of the factors affecting its descent of it.

Advancements in acoustic drug delivery for paranasal sinuses: A comprehensive review

Chronic rhinosinusitis (CRS) impacts patients' quality of life and healthcare costs. Traditional methods of drug delivery, such as nasal sprays and irrigation, have limited effectiveness. Acoustic Drug Delivery (ADD) using a nebulizer offers targeted delivery of drug to the sinuses, which may improve the treatment of CRS. This review examines the influence of aerosol particle characteristics, aero-acoustic parameters, inlet flow conditions, and acoustic waves on sinus drug delivery. Key findings reveal that smaller particles improve the ADD efficiency, whereas larger sizes or increased density impair it. The oscillation amplitude of the air plug in the ostium is crucial for the ADD efficiency. Introducing acoustic waves at the NC-sinus system's resonance frequency improves aerosol deposition within sinuses. Future research should address advanced models, optimizing particle characteristics, investigating novel acoustic waveforms, incorporating patient-specific anatomy, and evaluating long-term safety and efficacy. Tackling these challenges, ADD could offer more effective and targeted treatments for sinus-related conditions such as CRS.

Characterization of fungus ball CT-hyperdensities within maxillary and sphenoid sinuses

OBJECTIVES

CT-scan hyperdensities (HD) are described in more than 60% of all paranasal sinus fungus ball (FB) cases. Two types can be distinguished according to their density: calcium and metal types. We aimed to establish the prevalence and density of the HD observed in sphenoid and maxillary sinus FB and their relation to dental factors.

METHODS

This retrospective study included 64 patients operated in a tertiary referral center for unilateral maxillary or sphenoid FB diagnosed by histology or mycology. Pre-operative CT scans were analyzed by three independent observers (two ENT and one radiologist).

RESULTS

There were 45 maxillary FB and 19 sphenoid FB. 63 FB showed HD. Metal-type HD were observed in 28 maxillary FB but not in sphenoid sinuses. Among maxillary FB, the prevalence of endodontic treatment was significantly more significant on the FB side than on the healthy side (p = 0.02). The prevalence of endodontic treatment on the pathological side was more significant in the metal-type group than in the group without metal-type HD (p = 0.01). Isolated calcium-type HD were evidenced in 17 maxillary FB and 18 sphenoid FB (p = 0.019).

CONCLUSION

This study highlights the existence of two different types of HD in FBs of the paranasal sinuses with an association between metal-type HD and endodontic treatments.

Two-way coupling and Kolmogorov scales on inhaler spray plume evolutions from Ventolin, ProAir, and Qvar

Previous numerical studies of pulmonary drug delivery using metered-dose inhalers (MDIs) often neglected the momentum transfer from droplets to fluid. However, Kolmogorov length scales in MDI flows can be comparable to the droplet sizes in the orifice vicinity, and their interactions can modify the spray behaviors. This study aimed to evaluate the two-way coupling effects on spray plume evolutions compared to one-way coupling. The influences from the mass loading, droplet size, and inhaler type were also examined. Large-eddy simulation and Lagrangian approach were used to simulate the flow and droplet motions. Two-way coupled predictions appeared to provide significantly improved predictions of the aerosol behaviors close to the Ventolin orifice than one-way coupling. Increasing the applied MDI dose mass altered both the fluid and aerosol dynamics, notably bending the spray plume downward when applying a dose ten times larger. The droplet size played a key role in spray dynamics, with the plume being suppressed for 2-µm aerosols and enhanced for 20-µm aerosols. The Kolmogorov length scale ratio dp/η correlated well with the observed difference in spray plumes, with suppressed plumes when dp/η < 0.1 and enhanced plumes when dp/η > 0.1. For the three inhalers considered (Ventolin, ProAir, and Qvar), significant differences were predicted using two-way and one-way coupling despite the level and manifestation of these differences varied. Two-way coupling effects were significant for MDI sprays and should be considered in future numerical studies.

Count- and mass-based dosimetry of MDI spray droplets with polydisperse and monodisperse size distributions

Most previous numerical studies of inhalation drug delivery used monodisperse aerosols or quantified deposition as the ratio of deposited particle number over the total number of released particles (i.e., count-based). These practices are reasonable when the aerosols have a sufficiently narrow size range. However, spray droplets from metered-dose inhalers (MDIs) are often polydisperse with a wide size range, so using monodisperse aerosols and/or count-based deposition quantification may lead to significant errors. The objective of this study was to develop a mass-based dosimetry method and evaluate its performance in lung delivery in a mouth-lung (G9) geometry with an albuterol-CFC inhaler. The conventional practices (monodisperse and polydisperse-count-based) were also simulated for comparison purposes. The MDI actuation in the open space was studied using both high-speed imaging and LES-Lagrangian simulations. Experimentally measured spray velocities and size distribution were implemented in the computational model as boundary conditions. Good agreement was achieved between recorded and simulated spray plume evolution spatially and temporally. The polydisperse-mass-based predictions of MDI doses compared favorably with the measurements in all three regions considered (device, mouth-throat, and lung). Significant errors in MDI regional deposition were predicted using the monodisperse and count-based methods. The new polydisperse-mass-based method also predicted local deposition hot spots that were one order of magnitude higher in intensity than the two conventional methods. The results of this study highlighted that a presentative polydisperse size distribution and appropriate deposition quantification method should be applied to reliably predict the MDI drug delivery in the human respiratory tract.

Lower Inspiratory Breathing Depth Enhances Pulmonary Delivery Efficiency of ProAir Sprays

Effective pulmonary drug delivery using a metered-dose inhaler (MDI) requires a match between the MDI sprays, the patient’s breathing, and respiratory physiology. Different inhalers generate aerosols with distinct aerosol sizes and speeds, which require specific breathing coordination to achieve optimized delivery efficiency. Inability to perform the instructed breathing maneuver is one of the frequently reported issues during MDI applications; however, their effects on MDI dosimetry are unclear. The objective of this study is to systemically evaluate the effects of breathing depths on regional deposition in the respiratory tract using a ProAir-HFA inhaler. An integrated inhaler mouth-throat-lung geometry model was developed that extends to the ninth bifurcation (G9). Large-eddy simulation (LES) was used to compute the airflow dynamics due to concurrent inhalation and orifice flows. The discrete-phase Lagrangian model was used to track droplet motions. Experimental measurements of ProAir spray droplet sizes and speeds were used as initial and boundary conditions to develop the computational model for ProAir-pulmonary drug delivery. The time-varying spray plume from a ProAir-HFA inhaler into the open air was visualized using a high-speed imaging system and was further used to validate the computational model. The inhalation dosimetry of ProAir spray droplets in the respiratory tract was compared among five breathing depths on a regional, sub-regional, and local basis. The results show remarkable differences in airflow dynamics within the MDI mouthpiece and the droplet deposition distribution in the oral cavity. The inhalation depth had a positive relationship with the deposition in the mouth and a negative relationship with the deposition in the five lobes beyond G9 (small airways). The highest delivery efficiency to small airways was highest at 15 L/min and declined with an increasing inhalation depth. The drug loss inside the MDI was maximal at 45–60 L/min. Comparisons to previous experimental and numerical studies revealed a high dosimetry sensitivity to the inhaler type and patient breathing condition. Considering the appropriate inhalation waveform, spray actuation time, and spray properties (size and velocity) is essential to accurately predict inhalation dosimetry from MDIs. The results highlight the importance of personalized inhalation therapy to match the patient’s breathing patterns for optimal delivery efficiencies. Further complimentary in vitro or in vivo experiments are needed to validate the enhanced pulmonary delivery at 15 L/min.

Effect of MDI Actuation Timing on Inhalation Dosimetry in a Human Respiratory Tract Model

Accurate knowledge of the delivery of locally acting drug products, such as metered-dose inhaler (MDI) formulations, to large and small airways is essential to develop reliable in vitro/in vivo correlations (IVIVCs). However, challenges exist in modeling MDI delivery, due to the highly transient multiscale spray formation, the large variability in actuation–inhalation coordination, and the complex lung networks. The objective of this study was to develop/validate a computational MDI-releasing-delivery model and to evaluate the device actuation effects on the dose distribution with the newly developed model. An integrated MDI–mouth–lung (G9) geometry was developed. An albuterol MDI with the chlorofluorocarbon propellant was simulated with polydisperse aerosol size distribution measured by laser light scatter and aerosol discharge velocity derived from measurements taken while using a phase Doppler anemometry. The highly transient, multiscale airflow and droplet dynamics were simulated by using large eddy simulation (LES) and Lagrangian tracking with sufficiently fine computation mesh. A high-speed camera imaging of the MDI plume formation was conducted and compared with LES predictions. The aerosol discharge velocity at the MDI orifice was reversely determined to be 40 m/s based on the phase Doppler anemometry (PDA) measurements at two different locations from the mouthpiece. The LES-predicted instantaneous vortex structures and corresponding spray clouds resembled each other. There are three phases of the MDI plume evolution (discharging, dispersion, and dispensing), each with distinct features regardless of the actuation time. Good agreement was achieved between the predicted and measured doses in both the device, mouth–throat, and lung. Concerning the device–patient coordination, delayed MDI actuation increased drug deposition in the mouth and reduced drug delivery to the lung. Firing MDI before inhalation was found to increase drug loss in the device; however, it also reduced mouth–throat loss and increased lung doses in both the central and peripheral regions.

Recent In Vitro and In Silico Advances in the Understanding of Intranasal Drug Delivery

BACKGROUND

Many drugs are delivered intranasally for local or systemic effect, typically in the form of droplets or aerosols. Due to the high cost of in vivo studies, drug developers and researchers often turn to in vitro or in silico testing when first evaluating the behavior and properties of intranasal drug delivery devices and formulations. Recent advances in manufacturing and computer technologies have allowed for increasingly realistic and sophisticated in vitro and in silico reconstructions of the human nasal airways.

OBJECTIVE

The study aims to perform a summary of advances in the understanding of intranasal drug delivery based on recent in vitro and in silico studies.

CONCLUSION

The turbinates are a common target for local drug delivery applications, and while nasal sprays are able to reach this region, there is currently no broad consensus across the in vitro and in silico literature concerning optimal parameters for device design, formulation properties and patient technique which would maximize turbinate deposition. Nebulizers can more easily target the turbinates, but come with the disadvantage of significant lung deposition. Targeting of the olfactory region of the nasal cavity has been explored for the potential treatment of central nervous system conditions. Conventional intranasal devices, such as nasal sprays and nebulizers, deliver very little dose to the olfactory region. Recent progress in our understanding of intranasal delivery will be useful in the development of the next generation of intranasal drug delivery devices.

The impact of geometrical parameters on acoustically driven drug delivery to maxillary sinuses

Acoustically driven nebulized drug delivery (acoustic aerosol delivery) is the most efficient noninvasive technique for drug delivery to maxillary sinuses (MS). This method is based on the oscillation of the air plug inside the ostium to transport drug particles from the nasal cavity (NC) to the MS. The larger the wavelength of the air plug oscillation in the ostium, the greater the penetration of drug particles to the MS. However, using this technique, the maximum drug delivery efficiency achieved to date is 5%, which means 95% of the aerosolized drugs do not enter the MS and are wasted. Since the largest amplitude of the air plug oscillation occurs at its resonance frequency, to achieve an improved MS drug delivery efficiency, it is important to determine the resonance frequency of the nose-sinus combination accurately. This paper aims to investigate the impact of geometrical parameters on the resonance frequency of the nose-sinus model. Both experimental and computational acoustic models, along with the theoretical analysis, were conducted to determine the resonance frequency of an idealized nose-sinus model. The computational modeling was carried out using computational fluid dynamics (CFD) and finite element analysis (FEA), whereas in the analytical solution, the mathematical relationships developed for a conventional Helmholtz resonator were employed. A series of experiments were also conducted to measure the resonance frequency of a realistic NC-MS combination. The results demonstrated a good agreement between the experimental and CFD modeling, while the FEA and theoretical analysis showed a significant deviation from the experimental data. Also, it was shown that the resonance frequency of the idealized nose-sinus model increases by up to twofold with increasing the ostium diameter from 3 to 9 mm; however, it has an inverse relationship with the ostium length and sinus volume. It was also reported that the resonance frequency of the nose-sinus model is independent of the NC width and MS shape.

Numerical simulation of thermal water delivery in the human nasal cavity

This work describes an extensive numerical investigation of thermal water delivery for the treatment of inflammatory disorders in the human nasal cavity. The numerical simulation of the multiphase air-droplets flow is based upon the Large Eddy Simulation (LES) technique, with droplets of thermal water described via a Lagrangian approach. Droplet deposition is studied for different sizes of water droplets, corresponding to two different thermal treatments, i.e. aerosol and inhalation. Numerical simulations are conducted on a patient-specific anatomy, employing two different grid sizes, under steady inspiration at two breathing intensities. The results are compared with published in vivo and in vitro data. The effectiveness of the various thermal treatments is then assessed qualitatively and quantitatively, by a detailed analysis of the deposition patterns of the droplets. Discretization effects on the deposition dynamics are addressed. The level of detail of the present work, together with the accuracy afforded by the LES approach, leads to an improved understanding of how the mixture of air-water droplets is distributed within the nose and the paranasal sinuses.

Understanding the mechanisms underlying pulsating aerosol delivery to the maxillary sinus: In vitro tests and computational simulations

BACKGROUND

Pulsating aerosol delivery has been demonstrated in depositing medications into paranasal sinuses. However, its mechanisms are not fully understood. Influences of the nasal anatomy and sound frequency on intrasinus delivery are not yet clear.

OBJECTIVES

This study aimed to gain a better understanding of the mechanisms for enhanced intrasinus delivery with pulsating sound. Specifically, effects of the pulsation frequency, ostium size, and sinus shape on the intrasinus dosage and resonance frequency would be examined.

METHODS AND MATERIALS

Both experiments and computational modeling were conducted to understand the pulsating aerosol delivery in both idealized (two-bottle) and realistic nose-sinus models. A computational model of intrasinus pulsation delivery was developed using COMSOL and was cross-validated with both experimental and theoretical results.

RESULTS

In contrast to previous studies, seemingly erratic relations between the intrasinus dosage and ostium diameter were observed in experiments, which suggested a more complicated particle transport mechanism. Improved agreement was achieved when grouping the ostium size and sinus volume into the resonance frequency, and therefore, validated the hypothesis that intrasinus deposition strongly depends on the resonance frequency. Extensive computational simulations revealed that the deposition was highest at the resonance frequency and decreased gradually at off-resonance frequencies. The resonance frequency depended on the ostium and sinus morphology, but was independent of the nasal cavity.

CONCLUSION

Results of this study verified the hypothesis of resonance being the mechanism for enhanced particle deposition in the maxillary sinus. A better knowledge of the relationship between sinus dosages, pulsating frequency, and nasal morphometry is essential for improving the design of intrasinus delivery devices.

Modeling and Simulations of Olfactory Drug Delivery with Passive and Active Controls of Nasally Inhaled Pharmaceutical Aerosols

There are many advantages of direct nose-to-brain drug delivery in the treatment of neurological disorders. However, its application is limited by the extremely low delivery efficiency (< 1%) to the olfactory mucosa that directly connects the brain. It is crucial to develop novel techniques to deliver neurological medications more effectively to the olfactory region. The objective of this study is to develop a numerical platform to simulate and improve intranasal olfactory drug delivery. A coupled image-CFD method was presented that synthetized the image-based model development, quality meshing, fluid simulation, and magnetic particle tracking. With this method, performances of three intranasal delivery protocols were numerically assessed and compared. Influences of breathing maneuvers, magnet layout, magnetic field strength, drug release position, and particle size on the olfactory dosage were also numerically studied. From the simulations, we found that clinically significant olfactory dosage (up to 45%) were feasible using the combination of magnet layout and selective drug release. A 64 -fold higher delivery of dosage was predicted in the case with magnetophoretic guidance compared to the case without it. However, precise guidance of nasally inhaled aerosols to the olfactory region remains challenging due to the unstable nature of magnetophoresis, as well as the high sensitivity of olfactory dosage to patient-, device-, and particle-related factors.

Anatomical Details of the Rabbit Nasal Passages and Their Implications in Breathing, Air Conditioning, and Olfaction

The rabbit is commonly used as a laboratory animal for inhalation toxicology tests and detail knowledge of the rabbit airway morphometry is needed for outcome analysis or theoretical modeling. The objective of this study is to quantify the morphometric dimension of the nasal airway of a New Zealand white rabbit and to relate the morphology and functions through analytical and computational methods. Images of high-resolution MRI scans of the rabbit were processed to measure the axial distribution of the cross-sectional areas, perimeter, and complexity level. The lateral recess, which has functions other than respiration or olfaction, was isolated from the nasal airway and its dimension was quantified separately. A low Reynolds number turbulence model was implemented to simulate the airflow, heat transfer, vapor transport, and wall shear stress. Results of this study provide detailed morphological information of the rabbit that can be used in the studies of olfaction, inhalation toxicology, drug delivery, and physiology-based pharmacokinetics modeling. For the first time, we reported a spiral nasal vestibule that splits into three paths leading to the dorsal meatus, maxilloturbinate, and ventral meatus, respectively. Both non-dimensional functional analysis and CFD simulations suggested that the airflow in the rabbit nose is laminar and the unsteady effect is only significantly during sniffing. Due to the large surface-to-volume ratio, the maxilloturbinate is highly effective in warming and moistening the inhaled air to body conditions. The unique anatomical structure and respiratory airflow pattern may have important implications for designing new odorant detectors or electronic noses. Anat Rec, 299:853-868, 2016. © 2016 Wiley Periodicals, Inc.

Simulation study of electric-guided delivery of 0.4µm monodisperse and polydisperse aerosols to the ostiomeatal complex

Despite the high prevalence of rhinosinusitis, current inhalation therapy shows limited efficacy due to extremely low drug delivery efficiency to the paranasal sinuses. Novel intranasal delivery systems are needed to enhance targeted delivery to the sinus with therapeutic dosages. An optimization framework for intranasal drug delivery was developed to target polydisperse charged aerosols to the ostiomeatal complex (OMC) with electric guidance. The delivery efficiency of a group of charged aerosols recently reported in the literature was numerically assessed and optimized in an anatomically accurate nose-sinus model. Key design variables included particle charge number, particle size and distribution, electrode strength, and inhalation velocity. Both monodisperse and polydisperse aerosol profiles were considered. Results showed that the OMC delivery efficiency was highly sensitive to the applied electric field and electrostatic charges carried by the particles. Through the synthesis of electric-guidance and point drug release, focused deposition with significantly enhanced dosage in the OMC can be achieved. For 0.4 µm charged aerosols, an OMC delivery efficiency of 51.6% was predicted for monodisperse aerosols and 34.4% for polydisperse aerosols. This difference suggested that the aerosol profile exerted a notable effect on intranasal deliveries. Sensitivity analysis indicated that the OMC deposition fraction was highly sensitive to the charge and size of particles and was less sensitive to the inhalation velocity considered in this study. Experimental studies are needed to validate the numerically optimized designs. Further studies are warranted to investigate the targeted OMC delivery with both electric and acoustics controls, the latter of which has the potential to further deliver the drug particles into the sinus cavity.