Validating CFD Predictions of Pharmaceutical Aerosol Deposition with In Vivo Data

A quasi-realistic computational model development and flow field study of the human upper and central airways

The impact of drug delivery and particulate matter exposure on the human respiratory tract is influenced by various anatomical and physiological factors, particularly the structure of the respiratory tract and its fluid dynamics. This study employs computational fluid dynamics (CFD) to investigate airflow in two 3D models of the human air conducting zone. The first model uses a combination of CT-scan images and geometrical data from human cadaver to extract the upper and central airways down to the ninth generation, while the second model develops the lung airways from the first Carina to the end of the ninth generation using Kitaoka's deterministic algorithm. The study examines the differences in geometrical characteristics, airflow rates, velocity, Reynolds number, and pressure drops of both models in the inhalation and exhalation phases for different lobes and generations of the airways. From trachea to the ninth generation, the average air flowrates and Reynolds numbers exponentially decay in both models during inhalation and exhalation. The steady drop is the case for the average air velocity in Kitaoka's model, while that experiences a maximum in the 3rd or 4th generation in the quasi-realistic model. Besides, it is shown that the flow field remains laminar in the upper and central airways up to the total flow rate of 15 l/min. The results of this work can contribute to the understanding of flow behavior in upper respiratory tract.

Designing inhaled small molecule drugs for severe respiratory diseases: an overview of the challenges and opportunities

INTRODUCTION

Inhaled drugs offer advantages for the treatment of respiratory diseases over oral drugs by delivering the drug directly to the lung, thus improving the therapeutic index. There is an unmet medical need for novel therapies for lung diseases, exacerbated by a multitude of challenges for the design of inhaled small molecule drugs.

AREAS COVERED

The authors review the challenges and opportunities for the design of inhaled drugs for respiratory diseases with a focus on new target discovery, medicinal chemistry, and pharmacokinetic, pharmacodynamic, and toxicological evaluation of drug candidates.

EXPERT OPINION

Inhaled drug discovery is facing multiple unique challenges. Novel biological targets are scarce, as is the guidance for medicinal chemistry teams to design compounds with inhalation-compatible features. It is exceedingly difficult to establish a PK/PD relationship given the complexity of pulmonary PK and the impact of physical properties of the drug substance on PK. PK, PD and toxicology studies are technically challenging and require large amounts of drug substance. Despite the current challenges, the authors foresee that the design of inhaled drugs will be facilitated in the future by our increasing understanding of pathobiology, emerging medicinal chemistry guidelines, advances in drug formulation, PBPK models, and in vitro toxicology assays.

Enhancing Inhalation Drug Delivery: A Comparative Study and Design Optimization of a Novel Valved Holding Chamber

This paper presents an innovative approach to the design optimization of valved holding chambers (VHCs), crucial devices for aerosol drug delivery. We present the design of an optimal cylindrical VHC body and introduce a novel valve based on particle impaction theory. The research combines computational simulations and physical experiments to assess the performance of various VHCs, with a special focus on the deposition patterns of medication particles within these devices. The methodology incorporates both experimental and simulation approaches to validate the reliability of the simulation. Emphasis is placed on the deposition patterns observed on the VHC walls and the classification of fine and large particles for salbutamol sulfate particles. The study reveals the superior efficacy of our valve design in separating particles compared to commercially available VHCs. In standard conditions, our valve design allows over 95% of particles under 7 μm to pass through while effectively filtering those larger than 8 μm. The optimized body design accomplishes a 60% particle mass flow fraction at the outlet and an average particle size reduction of 58.5%. When compared numerically in terms of size reduction, the optimal design outperforms the two commercially available VHCs selected. This study provides valuable insights into the optimization of VHC design, offering significant potential for improved aerosol drug delivery. Our findings demonstrate a new path forward for future studies, aiming to further optimize the design and performance of VHCs for enhanced pulmonary drug delivery.

Validation of computational fluid dynamics models for airway deposition with SPECT data of the same population

This study compared computational fluid dynamic (CFD) model predictions on aerosol deposition in six asthmatic patients to the in-vivo results of the same group. Patient-specific ventilation and internal air distribution were prescribed using inspiratory and expiratory CT scans of each patient, accounting for individual lobar air flow distribution. Moreover, the significant influence of realistic mouth and throat geometries on regional deposition is demonstrated. The in-vivo data were obtained from single photon emission computed tomography (SPECT) in 6 subjects with mild asthma selected from a database of historical clinical trials. The governing flow and particle tracking equations were solved numerically using a commercial CFD tool, and the modeled deposition results were compared to the SPECT data. Good agreement was found between the CFD model, applying k-ω SST turbulence model, and SPECT in terms of aerosol deposition. The average difference for the lobar deposition obtained from CFD model and SPECT/CT data was 2.1%. The high level of agreement is due to applying patient specific airway geometries and inspiratory/expiratory CT images, anatomical upper airways, and realistic airway trees. The results of this study show that CFD is a powerful tool to simulate patient-specific deposition if correct boundary conditions are applied and can generate similar information obtained with functional imaging tools, such as SPECT.

Validated respiratory drug deposition predictions from 2D and 3D medical images with statistical shape models and convolutional neural networks

For the one billion sufferers of respiratory disease, managing their disease with inhalers crucially influences their quality of life. Generic treatment plans could be improved with the aid of computational models that account for patient-specific features such as breathing pattern, lung pathology and morphology. Therefore, we aim to develop and validate an automated computational framework for patient-specific deposition modelling. To that end, an image processing approach is proposed that could produce 3D patient respiratory geometries from 2D chest X-rays and 3D CT images. We evaluated the airway and lung morphology produced by our image processing framework, and assessed deposition compared to in vivo data. The 2D-to-3D image processing reproduces airway diameter to 9% median error compared to ground truth segmentations, but is sensitive to outliers of up to 33% due to lung outline noise. Predicted regional deposition gave 5% median error compared to in vivo measurements. The proposed framework is capable of providing patient-specific deposition measurements for varying treatments, to determine which treatment would best satisfy the needs imposed by each patient (such as disease and lung/airway morphology). Integration of patient-specific modelling into clinical practice as an additional decision-making tool could optimise treatment plans and lower the burden of respiratory diseases.

Use of the Same Model or Modeling Strategy Across Multiple Submissions: Focus on Complex Drug Products

Evidence shows that there is an increasing use of modeling and simulation to support product development and approval for complex generic drug products in the USA, which includes the use of mechanistic modeling and model-integrated evidence (MIE). The potential for model reuse was the subject of a workshop session summarized in this review, where the session included presentations and a panel discussion from members of the U.S. Food and Drug Administration (FDA), academia, and the generic drug product industry. Concepts such as platform performance assessment and MIE standardization were introduced to provide potential frameworks for model reuse related to mechanistic models and MIE, respectively. The capability of models to capture formulation and product differences was explored, and challenges with model validation were addressed for drug product classes including topical, orally inhaled, ophthalmic, and long-acting injectable drug products. An emphasis was placed on the need for communication between FDA and the generic drug industry to continue to foster maturation of modeling and simulation that may support complex generic drug product development and approval, via meetings and published guidance from FDA. The workshop session provided a snapshot of the current state of modeling and simulation for complex generic drug products and offered opportunities to explore the use of such models across multiple drug products.

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.

Effects of vocal fold adduction on the particle deposition in the glottis: A numerical analysis and in vitro assessment

BACKGROUND

The efficacy of inhalation therapy depends on the drug deposition in the human respiratory tract. This study investigates the effects of vocal fold adduction on the particle deposition in the glottis.

METHODS

A realistic mouth-throat (MT) geometry was built based on CT images of a healthy adult (MT-A). Mild (MT-B) and great (MT-C) vocal fold (VF) adduction were incorporated in the original model. Monodisperse particles range in size from 3 to 12 μm were simulated at inspiration flow rates of 15, 30 and 45 L per minute (LPM). The regional deposition of drug aerosols was performed in 3D-printed models and quantified using high-performance liquid chromatography.

RESULTS

Both the numerical analysis and in vitro experiments show that most particles are deposited in the mouth, pharynx and supraglottis, while few are deposited in the glottis and subglottis. For most cases in MT-A, the particle quantity in glottis is lower than 0.02 N/mm2 at 15 and 30 LPM while they increase dramatically at 45 LPM. It peaked at 0.347 N/mm2 for 5-μm particles at 45 LPM in MT-B and 2.324 N/mm2 for 6-μm particles at 30 LPM in MT-C. The lowest drug mass faction in the glottis in vitro were found at 15 LPM for MT-A and MT-C, and at 30 LPM for MT-B, whereas it peaked at 45 LPM for all MT models, 0.71% in MT-A, 1.16% in MT-B, and 2.53% in MT-C, respectively.

CONCLUSION

Based on the results of this study, larger particles are more likely to be deposited in the oral cavity, oropharynx, and supraglottis than in the glottis. However, particle deposition in the glottis generally increases with VF adduction and greater inspiratory flow rates.

The Pharmacokinetics of Inhaled Drugs

The pharmacokinetic (PK) profile of a drug after inhalation may differ quite markedly from that seen after dosing by other routes of administration. Drugs may be administered to the lung to elicit a local action or as a portal for systemic delivery of the drug to its site of action elsewhere in the body. Some knowledge of PK is important for both locally- and systemically-acting drugs. For a systemically-acting drug, the plasma concentration-time profile shares some similarities with drug given by the oral or intravenous routes, since the plasma concentrations (after the distribution phase) will be in equilibrium with concentrations at the site of action. For a locally-acting drug, however, the plasma concentrations reflect its fate after it has been absorbed and removed from the airways, and not what is available to its site of action in the lung. Consequently, those typical PK parameters which are determined from plasma concentration measurements, e.g., area under the curve (AUC), Cmax, tmax and post-peak t1/2 may provide information on the deposition and absorption of drugs from the lung; however, the information from these parameters becomes more complicated to decipher for those drugs which are locally-acting in the lung. Additionally, the plasma concentration profile for both locally- and systemically-acting drugs will not only reflect drug absorbed from the lung but also that absorbed from the gastrointestinal (GI) tract from the portion of the dose which is swallowed. This absorption from the GI tract adds a further complication to the interpretation of plasma concentrations, particularly for locally-acting drugs. The influence of physiological and pathological factors needs to be considered in the absorption of some inhaled drugs. The absorption of some hydrophilic drugs is influenced by the inspiratory maneuver used during initial inhalation of the drug, and at later times after deposition. Similarly, the effects of smoking have been shown to increase lung permeability and increase the absorption of certain hydrophilic drugs. The effects of different disease states of the lung have less defined influences on absorption into the systemic circulation.

Dry Powder Formulations for Inhalation Require a Smaller Aerodynamic Diameter for Usage at High Altitude

BACKGROUND

High Altitude Pulmonary Edema (HAPE) seriously threatens the health of people at high altitudes. There are drug treatments for HAPE, and dry powder formulations (DPFs) represent a rapid and accessible delivery vehicle for these drugs. However, there are presently no reports on the inhalability of DPFs in low-pressure environments. Given the reduced atmospheric pressure typical at high altitudes, conventional DPFs might not be suitable for inhalation. Therefore, it is necessary to elucidate the deposition behaviors of dry powder in the respiratory tract at low pressure, as well as to improve their pulmonary deposition efficiency via adjustments to their formulation and design.

METHODS

The effect of air pressure, inspiratory velocity, and particle properties (such as size, density, and aerodynamic diameter) on pulmonary deposition of DPFs was calculated by a computational fluid dynamics (CFD)-coupled discrete phase model. DPFs of various aerodynamic diameters were prepared by spray drying, and the inhalability of these DPFs in a low-pressure environment was evaluated in mice. Finally, a mouse model of HAPE was established, and the treatment of HAPE by nifedipine-loaded DPFs with small aerodynamic diameter was validated.

RESULTS

CFD results showed that low pressure decreased the deposition of DPFs in the lungs. At 0.5 standard atmosphere, DPFs with aerodynamic diameter of ∼2.0 μm could not enter the lower respiratory tract; however, a decrease in the physical diameter, density, and, consequently, the aerodynamic diameter of the DPFs was able to enhance pulmonary deposition of these powders. To validate the CFD results, three kinds of dry powder with aerodynamic diameters of 0.66, 0.98, and 2.00 μm were prepared by spray drying. Powders with smaller aerodynamic diameter could be inhaled into the lungs of mice more effectively, and, consequently could ameliorate the progression of HAPE more effectively than conventional powders. These results were consistent with the CFD results.

CONCLUSIONS

Low atmospheric pressure can prevent the pulmonary deposition of DPFs at high altitudes. Compared with conventional DPFs, powders with smaller aerodynamic diameter can be effectively inhaled at these pressures and thus might be more suitable for the treatment the HAPE.

Predicting Inhaled Drug Dose Generated by Mesh Nebulizers

Background: The lung dose of nebulized drugs for spontaneous breathing is influenced by breathing patterns and nebulizer performance. This study aimed to develop a system for measuring breath patterns and a formula for estimating inhaled drugs, and then to validate the hypothesized prediction formula. Methods: An in vitro model was first used to determine correlations among the delivered dose, breath patterns, and doses deposited on the accessories and reservoirs testing with a breathing simulator to generate 12 adult breathing patterns (n = 5). A pressure sensor was developed to measure breathing parameters and used along with a prediction formula that accounted for the initial charge dose, respiratory pattern, and dose on the accessory and reservoir of a nebulizer. Three brands of nebulizers were tested by placing salbutamol (5.0 mg/2.5 mL) in the drug holding chamber. Ten healthy individuals participated in the ex vivo study to validate the prediction formula. The agreement between the predicted and inhaled doses was analyzed using the Bland-Altman plot. Results: The in vitro model showed that the inspiratory time to total respiratory cycle time (Ti/Ttotal; %) was significantly directly correlated with the delivered dose among the respiratory factors, followed by inspiratory flow, respiratory rate, and tidal volume. The ex vivo model showed that Ti/Ttotal was significantly directly correlated with the delivered dose among the respiratory factors, in addition to the nebulization time and accessory dose. The Bland-Altman plots for the ex vivo model showed similar results between the two methods. Large differences in inhaled dose measured at the mouth were observed among the subjects, ranging from 12.68% to 21.68%; however, the difference between the predicted dose and inhaled dose was lower, at 3.98%-5.02%. Conclusions: The inhaled drug dose could be predicted with the hypothesized estimation formula, which was validated by the agreement between the inhaled and predicted doses of breathing patterns of healthy individuals.

A path forward in the development of new aerosol drug delivery devices for pediatrics

Inhaled medications are widely accepted as being the optimal route for treating pediatric respiratory diseases, a leading cause of hospitalization and death. Despite jet nebulizers being the preferred inhalation device for neonates and infants, current devices face performance issues with most of the drug never reaching the target lung location. Previous work has aimed to improve pulmonary drug deposition, yet nebulizer efficiency remains low. The development of an inhalant therapy that is efficacious and safe for pediatrics depends on a well-designed delivery system and formulation. To accomplish this, the field needs to rethink the current practice of basing pediatric treatments on adult studies. The rapidly evolving pediatric patient (i.e. neonates to eighteen) needs to be considered because they are different from adults with respect to airway anatomy, breathing patterns, and adherence. Previous research approaches to improve deposition efficiency have been limited due to the complexity of combining physics, which drives aerosol transport and deposition, and biology, especially within the area of pediatrics. To address these critical knowledge gaps, we need a better understanding of how patient age and disease state affect deposition of aerosolized drugs. The complexity of the multiscale respiratory system makes scientific investigation very challenging. The authors have simplified the complex problem into five components with these three areas as ones to address first: how the aerosol is (i) generated in a medical device, (ii) delivered to the patient, and (iii) deposited inside the lung. In this review, we discuss the technological advances and innovations made from experiments, simulations, and predictive models in each of these areas. In addition, we discuss the impact on patient treatment efficacy and recommend a clinical direction, with a focus on pediatrics. In each area, a series of research questions are posed and steps for future research to improve efficacy in aerosol drug delivery are outlined.

Advancements in the Design and Development of Dry Powder Inhalers and Potential Implications for Generic Development

Dry powder inhalers (DPIs) are drug-device combination products where the complexity of the formulation, its interaction with the device, and input from users play important roles in the drug delivery. As the landscape of DPI products advances with new powder formulations and novel device designs, understanding how these advancements impact performance can aid in developing generics that are therapeutically equivalent to the reference listed drug (RLD) products. This review details the current understanding of the formulation and device related principles driving DPI performance, past and present research efforts to characterize these performance factors, and the implications that advances in formulation and device design may present for evaluating bioequivalence (BE) for generic development.

Scientific and regulatory activities initiated by the U.S. Food and drug administration to foster approvals of generic dry powder inhalers: Bioequivalence perspective

Regulatory science for generic dry powder inhalers (DPIs) in the United States (U.S.) has evolved over the last decade. In 2013, the U.S. Food and Drug Administration (FDA) published the draft product-specific guidance (PSG) for fluticasone propionate and salmeterol xinafoate inhalation powder. This was the first PSG for a DPI available in the U.S., which provided details on a weight-of-evidence approach for establishing bioequivalence (BE). A variety of research activities including in vivo and in vitro studies were used to support these recommendations, which have led to the first approval of a generic DPI in the U.S. for fluticasone propionate and salmeterol xinafoate inhalation powder in January of 2019. This review describes the scientific and regulatory activities that have been initiated by FDA to support the current BE recommendations for DPIs that led to the first generic DPI approvals, as well as research with novel in vitro and in silico methods that may potentially facilitate generic DPI development and approval.

A computed tomography imaging-based subject-specific whole-lung deposition model

The respiratory tract is an important route for beneficial drug aerosol or harmful particulate matter to enter the body. To assess the therapeutic response or disease risk, whole-lung deposition models have been developed, but were limited by compartment, symmetry or stochastic approaches. In this work, we proposed an imaging-based subject-specific whole-lung deposition model. The geometries of airways and lobes were segmented from computed tomography (CT) lung images at total lung capacity (TLC), and the regional air-volume changes were calculated by registering CT images at TLC and functional residual capacity (FRC). The geometries were used to create the structure of entire subject-specific conducting airways and acinar units. The air-volume changes were used to estimate the function of subject-specific ventilation distributions among acinar units and regulate flow rates in respiratory airway models. With the airway dimensions rescaled to a desired lung volume and the airflow field simulated by a computational fluid dynamics model, particle deposition fractions were calculated using deposition probability formulae adjusted with an enhancement factor to account for the effects of secondary flow and airway geometry in proximal airways. The proposed model was validated in silico against existing whole-lung deposition models, three-dimensional (3D) computational fluid and particle dynamics (CFPD) for an acinar unit, and 3D CFPD deep lung model comprising conducting and respiratory regions. The model was further validated in vivo against the lobar particle distribution and the coefficient of variation of particle distribution obtained from CT and single-photon emission computed tomography (SPECT) images, showing good agreement. Subject-specific airway structure increased the deposition fraction of 10.0-μm particles and 0.01-μm particles by approximately 10%. An enhancement factor increased the overall deposition fractions, especially for particle sizes between 0.1 and 1.0 μm.

Characterizing the Effects of Nasal Prong Interfaces on Aerosol Deposition in a Preterm Infant Nasal Model

The objective of this study was to characterize the effects of multiple nasal prong interface configurations on nasal depositional loss of pharmaceutical aerosols in a preterm infant nose-throat (NT) airway model. Benchmark in vitro experiments were performed in which a spray-dried powder formulation was delivered to a new preterm NT model with a positive-pressure infant air-jet dry powder inhaler using single- and dual-prong interfaces. These results were used to develop and validate a computational fluid dynamics (CFD) model of aerosol transport and deposition in the NT geometry. The validated CFD model was then used to explore the NT depositional characteristic of multiple prong types and configurations. The CFD model highlighted a turbulent jet effect emanating from the prong(s). Analysis of NT aerosol deposition efficiency curves for a characteristic particle size and delivery flowrate (3 µm and 1.4 L/min (LPM)) revealed little difference in NT aerosol deposition fraction (DF) across the prong insertion depths of 2-5 mm (DF = 16-24%) with the exception of a single prong with 5-mm insertion (DF = 36%). Dual prongs provided a modest reduction in deposition vs. a single aerosol delivery prong at the same flow for insertion depths < 5 mm. The presence of the prongs increased nasal depositional loss by absolute differences in the range of 20-70% compared with existing correlations for ambient aerosols. In conclusion, the use of nasal prongs was shown to have a significant impact on infant NT aerosol depositional loss prompting the need for prong design alterations to improve lung delivery efficiency.

Analytical method using SEM-EDS for metal elements present in particulate matter generated from stainless steel flux-cored arc welding process

Welding fumes (WFs) can cause occupational pneumonoconiosis and other diseases in workers. WFs have complex chemical composition and morphology depending on the welding conditions. The WF surface is a key factor affecting those diseases. The objective of this study was to establish an analytical method focused on characterizing individual WFs and welding slags (WSs) formed during CO₂ arc welding processes for knowledge acquisition of risk assessment. Especially, the characterization was focused on the elemental distributions near the surfaces obtained using fluxing agents and size of the WFs. WFs were collected using personal samplers. After welding, WS was also collected. The fluxing elemental distribution (e.g., Bi) near the surfaces WS and WFs were analyzed through scanning electron microscopy and energy-dispersive X-ray spectroscopy. As a result, some of the micron-sized spherical particles (SPs) grew by incorporating nanosized primary particles composed of other metal species. The fluxing agents formed elemental distribution patterns on the SP surface. Bi were dotted in an agglomerate. Mn amount in WS depends on Mn amount in the WFs. These results obtained through the analysis of both the WS and WF surface as well as the particle sizes will facilitate the establishment of exposure assessment models.

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.

Evaluation of the Polyhedral Mesh Style for Predicting Aerosol Deposition in Representative Models of the Conducting Airways

A critical factor affecting the accuracy of Computational Fluid Dynamic (CFD) simulations and the time required to conduct them is construction of the computational mesh. This study aimed to evaluate the relatively new polyhedral mesh style for simulating aerosol deposition in the upper conducting airways compared with established meshing techniques and experimental data. Hexahedral and polyhedral mesh solutions were compared in two benchmark geometries: 1) a 90°-bend with flow characteristics similar to the extrathoracic airways of an adolescent child, and 2) a double bifurcation representing bifurcations B3-B5 in an adult. Both 4-block and 5-block hexahedral meshes were used in the 90°-bend to capture the potential of fully-structured hexahedral meshes. In the 90°-bend, polyhedral elements matched polydisperse in vitro deposition data with 20% relative error (RE; averaged across the particle sizes considered), which is an improvement on the accuracy of the 4-block hexahedral mesh (35% RE) and is similar to the accuracy of the 5-block hexahedral mesh (19% RE). In the double bifurcation, deposition fraction relative differences evaluated between polyhedral and hexahedral meshes ranged from 0.3% to 28.6% for the different particle sizes assessed, which is an order of magnitude improvement compared with previous studies that considered hexahedral vs. hybrid tetrahedral-prism meshes for the same flow field. Solution convergence time with polyhedral elements was found to be 50% to 140% higher than with hexahedral meshes of comparable size. While application dependent, the increase in simulation time observed with polyhedral meshes will likely be outweighed by the ease and convenience of polyhedral mesh construction. It was concluded that the polyhedral mesh style, with sufficient resolution especially near the walls, is an excellent alternative to the highly regarded hexahedral mesh style for predicting upper airway aerosol transport and deposition and provides a powerful new tool in the assessment of respiratory aerosol dosimetry.

Reconciling Oxygen and Aerosol Delivery with a Hood on In Vitro Infant and Paediatric Models

This study aimed to evaluate optimal aerosol and oxygen delivery with a hood on an infant model and a paediatric model. A facemask and a hood with three inlets, with or without a front cover, were used. A small-volume nebuliser with a unit-dose of salbutamol was used for drug delivery and an air entrainment nebuliser was used to deliver oxygen at 35%. Infant and paediatric breathing patterns were mimicked; a bacterial filter was connected to the end of a manikin trachea for aerosol drug collection, and an oxygen analyser was used to measure the oxygen concentration. For the infant model, inhaled drug dose was significantly higher when the nebuliser was placed in the back of the hood and with a front cover. This was verified by complementary computational simulations in a comparable infant-hood model. For the paediatric model, the inhaled dose was greater with a facemask than with a hood. Oxygen delivery with a facemask and a hood with a front cover achieved a set concentration in both models, yet a hood without a front cover delivered oxygen at far lower concentrations than the set concentration.

Development of a rat capnoperitoneum phantom to study drug aerosol deposition in the context of anticancer research on peritoneal carcinomatosis

Pressurized Intraperitoneal Aerosol Chemotherapy (PIPAC) is a promising approach with a high optimization potential for the treatment of peritoneal carcinomatosis. To study the efficacy of PIPAC and drugs, first rodent cancer models were developed. But inefficient drug aerosol supply and knowledge gaps concerning spatial drug distribution can limit the results based on such models. To study drug aerosol supply/deposition, computed tomography scans of a rat capnoperitoneum were used to deduce a virtual and a physical phantom of the rat capnoperitoneum (RCP). RCP qualification was performed for a specific PIPAC method, where the capnoperitoneum is continuously purged by the drug aerosol. In this context, also in-silico analyses by computational fluid dynamic modelling were conducted on the virtual RCP. The physical RCP was used for ex-vivo granulometric analyses concerning drug deposition. Results of RCP qualification show that aerosol deposition in a continuous purged rat capnoperitoneum depends strongly on the position of the inlet and outlet port. Moreover, it could be shown that the droplet size and charge condition of the drug aerosol define the deposition efficiency. In summary, the developed virtual and physical RCP enables detailed in-silico and ex-vivo analyses on drug supply/deposition in rodents.

Computational pharmaceutics - A new paradigm of drug delivery

In recent decades pharmaceutics and drug delivery have become increasingly critical in the pharmaceutical industry due to longer time, higher cost, and less productivity of new molecular entities (NMEs). However, current formulation development still relies on traditional trial-and-error experiments, which are time-consuming, costly, and unpredictable. With the exponential growth of computing capability and algorithms, in recent ten years, a new discipline named "computational pharmaceutics" integrates with big data, artificial intelligence, and multi-scale modeling techniques into pharmaceutics, which offered great potential to shift the paradigm of drug delivery. Computational pharmaceutics can provide multi-scale lenses to pharmaceutical scientists, revealing physical, chemical, mathematical, and data-driven details ranging across pre-formulation studies, formulation screening, in vivo prediction in the human body, and precision medicine in the clinic. The present paper provides a comprehensive and detailed review in all areas of computational pharmaceutics and "Pharma 4.0", including artificial intelligence and machine learning algorithms, molecular modeling, mathematical modeling, process simulation, and physiologically based pharmacokinetic (PBPK) modeling. We not only summarized the theories and progress of these technologies but also discussed the regulatory requirements, current challenges, and future perspectives in the area, such as talent training and a culture change in the future pharmaceutical industry.

Salbutamol Transport and Deposition in the Upper and Lower Airway with Different Devices in Cats: A Computational Fluid Dynamics Approach

Simple Summary

Administration of inhaled salbutamol via metered-dose inhalers can effectively treat bronchoconstriction. Different devices are used for the delivery of this drug in cats, either in the hospital or at home, for long-term treatment. Effective drug administration may depend on the drug delivery device as well as patient cooperation. By using non-invasive computational fluid dynamics techniques, the impact of these devices on the deposition and transport of salbutamol particles in the cat airways was simulated and assessed. The results confirm a variable drug distribution depending on the device used. The percentage of particles reaching the lung was reduced when using spacers and increased when applied directly into an endotracheal tube.

Abstract

Pressurized metered-dose inhalers (pMDI) with or without spacers are commonly used for the treatment of feline inflammatory airway disease. During traditional airways treatments, a substantial amount of drugs are wasted upstream of their target. To study the efficiency of commonly used devices in the transport of inhaled salbutamol, different computational models based on two healthy adult client-owned cats were developed. Computed tomographic images from one cat were used to generate a three-dimensional geometry, and two masks (spherical and conical shapes) and two spacers (10 and 20 cm) completed the models. A second cat was used to generate a second model having an endotracheal tube (ETT) with and without the same spacers. Airflow, droplet spray transport, and deposition were simulated and studied using computational fluid dynamics techniques. Four regions were evaluated: device, upper airways, primary bronchi, and downstream lower airways/parenchyma (“lung”). Regardless of the model, most salbutamol is deposited in devices and/or upper airways. In general, particles reaching the lung varied between 5.8 and 25.8%. Compared with the first model, pMDI application through the ETT with or without a spacer had significantly higher percentages of particles reaching the lung (p = 0.006).

Aerosol deposition and airflow dynamics in healthy and asthmatic human airways during inhalation

Inhalation of aerosols such as pharmaceutical aerosols or virus aerosol uptake is of great concern to the human population. To elucidate the underlying aerosol dynamics, the deposition fractions (DFs) of aerosols in healthy and asthmatic human airways of generations 13-15 are predicted. The Navier-stokes equations governing the gaseous phase and the discrete phase model for particles' motion are solved using numerical methods. The main forces responsible for deposition are inertial impaction forces and complex secondary flow velocities. The curvatures and sinusoidal folds in the asthmatic geometry lead to the formation of complex secondary flows and hence higher DFs. The intensities of complex secondary flows are strongest at the generations affected by asthma. The DF in the healthy airways is 0%, and it ranges from 1.69% to 52.93% in the asthmatic ones. From this study, the effects of the pharmaceutical aerosol particle diameters in the treatment of asthma patients can be established, which is conducive to inhibiting the inflammation of asthma airways. Furthermore, with the recent development of COVID-19 which causes pneumonia, the predicted physics and effective simulation methods of bioaerosols delivery to asthma patients are vital to prevent the exacerbation of the chronic ailment and the epidemic.

Deposition of Smoke Particles in Human Airways with Realistic Waveform

Exposure to toxic particles from smoke generated either from bush fire, stable burning, or direct smoking is very harmful to our health. The tiny particles easily penetrate deep into the lungs after exposure and damage the airways. Tobacco smoking causes the direct emission of 2.6 million tons of CO2 and 5.2 million tons of methane annually into the atmosphere. Nevertheless, it is one of the significant contributors to various respiratory diseases leading to lung cancer. These particles’ deposition in the human airway is computed in the present article for refining our understanding of the adverse health effects due to smoke particle inhalation, especially cigarette smoke. Until recently, little work has been reported to account for the transient flow pattern of cigarette smoking. Consideration of transient flow may change the deposition pattern of the particle. A high-resolution CT scan image of the respiratory tract model consisting of the oral cavity, throat, trachea, and first to sixth generations of the lungs helps predict cigarette smoke particle (CSP) deposition. With the same scan, a realistic geometric model of the human airways of an adult subject is used to simulate the transport of air and particle. The CSP deposition is determined at different locations from the oral cavity to the sixth generation of the bronchi. In addition, an unsteady breathing curve indicative of realistic smoking behavior is utilized to represent the breathing conditions accurately. The discrete phase model (DPM) technique is used to determine smoke particle deposition in the human airways. It is found that the deposition increases with the size of the smoke particle. Particles tend to deposit in the oral cavity around the bifurcation junction of the airways. The deposition fraction of CSP with the realistic waveform of smoking is found to be smaller compared to that during the stable flow condition. It is also observed that the fine particles (0.1–1.0 micron) escape to lower generations, leading to higher deposition of fine particles in the deeper airways. The outcome of the study is helpful for understanding smoke-related pulmonary complications.

A quasi-3D model of the whole lung: airway extension to the tracheobronchial limit using the constrained constructive optimization and alveolar modeling, using a sac-trumpet model

Existing computational models used for simulating the flow and species transport in the human airways are zero-dimensional (0D) compartmental, three-dimensional (3D) computational fluid dynamics (CFD), or the recently developed quasi-3D (Q3D) models. Unlike compartmental models, the full CFD and Q3D models are physiologically and anatomically consistent in the mouth and the upper airways, since the starting point of these models is the mouth-lung surface geometry, typically created from computed tomography (CT) scans. However, the current resolution of CT scans limits the airway detection between the 3rd-4th and 7th-9th generations. Consequently, CFD and the Q3D models developed using these scans are generally limited to these generations. In this study, we developed a method to extend the conducting airways from the end of the truncated Q3D lung to the tracheobronchial (TB) limit. We grew the lung generations within the closed lung lobes using the modified constrained constructive optimization, creating an aerodynamically optimized network aiming to produce equal pressure at the distal ends of the terminal segments. This resulted in a TB volume and lateral area of ∼165 cc and ∼2000 cm², respectively. We created a "sac-trumpet" model at each of the TB outlets to represent the alveoli. The volumes of the airways and the individual alveolar generations match the anatomical values by design: with the functional residual capacity at 2611 cc. Lateral surface areas were scaled to match the physiological values. These generated Q3D whole lung models can be efficiently used for conducting multiple breathing cycles of drug transport and deposition simulations.

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.

Septal destruction enhances chaotic mixing and increases cellular doses of nanoparticles in emphysematous acinus

One hallmark of emphysema is the breakdown of inter-alveolar septal walls in pulmonary acini. How the acinar dosimetry of environmental aerosols varies at different stages of emphysema remains unclear; this is specifically pertinent to users of tobacco products, which is the leading cause of emphysema. The objective of this study is to systematically assess the impacts of septal destruction on the behavior and fate of nanoparticles (1–800 nm) in a pyramid-shaped sub-acinar model consisting of 496 alveoli. Four diseased geometry variants were created by gradually removing the septal walls from the base model. Particle motions within the acinar region were tracked for particles raging 1–800 nm at four emphysema stages using a well-tested Lagrangian tracking model. Both spatial profile and temporal variation of particle deposition were predicted in healthy and diseased sub-acinar geometries on both a total and regional basis. Results show large differences in airflow and particle dynamics among different emphysema stages. Large differences in particle dynamics are also observed among different particle sizes, with one order of magnitude’s variation in the speeds of particles of 1, 10, and 200 nm. The destruction of septal walls also changed the deposition mechanisms, shifting from connective diffusion to chaotic mixing with emphysema progression. The sub-acinar dosimetry became less sensitive to particle size variation with more septal destructions. The lowest retention rate was found at 200–500 nm in the healthy sub-acinar geometry, but at 800 nm in all emphysematous models considered. The acinus-averaged dose for nanoparticles (1–800 nm) increases with aggravating septal destructions, indicating an even higher risk to the acinus at later emphysema stages.

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.

A multiscale absorption and transit model for oral delivery of hydroxychloroquine: Pharmacokinetic modeling and intestinal concentration prediction to assess toxicity and drug-induced damage in healthy subjects

Hydroxychloroquine (HCQ) is commonly used in the treatment of malaria and rheumatic diseases. Recently it has also been identified as possible therapeutic option in combating COVID-19. However, the use of HCQ is known to induce cytotoxicity. In 2020, we developed a multiscale absorption and transit (MAT) toolkit to simulate the dissolution, transport, absorption, distribution, metabolism, and elimination of orally administered drugs in the human GIT at multiple levels. MAT was constructed by integrating the spatially accurate first-principles driven high-fidelity drug transport, dissolution, and absorption model in the human stomach and GIT using the recently published Quasi-3D framework. The computational results showed that MAT was able to match the experimental concentration results better than the traditional compartmental models. In this study, we adapted MAT, to predict the pharmacokinetics of orally delivered HCQ in healthy subjects. The computational results matched the experimental concentration results. The simulated stomach and intestinal fluid and enterocyte concentrations were compared with the in vitro CC50 values. While the peak enterocyte concentrations were several orders lower than the in vitro CC50 values, the peak stomach and the intestinal fluid concentrations were only one order smaller than the in vitro CC50 values. In particular, the peak stomach and the duodenum fluid concentrations were just 3× smaller than the in vitro CC50 values. This implies that the lumen walls are much more susceptible to cytotoxicity-based damage than the enterocyte layers. We envision that MAT can be used to optimize the dosing regimen of HCQ by maximizing its bioavailability, while simultaneously minimizing the cytotoxic damage.

Computational Fluid Dynamics (CFD)-Based Optimization of Injection Process during Endoscopic Mucosal Therapy

Creation of a submucosal plane to separate the lesion from the deeper muscle layer in gastrointestinal tract is an integral and essential part of endoscopic resection therapies such as endoscopic mucosal resection (EMR) and endoscopic submucosal dissection (ESD). Thereby, an optimized submucosal injection technique is required to ensure a successful process. In this study, the computational fluid dynamics (CFD) technique is employed as a foundational step towards the development of a framework that can provide useful directions to optimize the injection process. Three different lifting agents, including Glycerol, Eleview^®, and ORISE^® gel have been used for this study. The role of different injection angles, injection dynamics, and effect of temperature are studied to understand the lifting characteristic of each agent. The study shows that Eleview^® provides the highest lifting effect, including the initial injection period. To evaluate the impact of the injection process, two cases are simulated, termed static injection and dynamic injection. Under static injection, the injection angle is investigated from lower to higher angles of injection. In the dynamic injection, two cases are modulated, where a continuous change of injection angle from lower to higher degrees (denoted as clockwise) and vice-versa in the anti-clockwise direction are investigated. Increased lifting characteristics are observed at decreasing/lower angle of injection. Further, the correlation between temperature of the lifting agents and their lifting characteristics is investigated.

Contemporary Formulation Development for Inhaled Pharmaceuticals

Pulmonary delivery has gained increased interests over the past few decades. For respiratory conditions, targeted drug delivery directly to the site of action can achieve a high local concentration for efficacy with reduced systemic exposure and adverse effects. For systemic conditions, the unique physiology of the lung evolutionarily designed for rapid gaseous exchange presents an entry route for systemic drug delivery. Although the development of inhaled formulations has come a long way over the last few decades, many aspects of it remain to be elucidated. In particular, a reliable and well-understood method for in vitro-in vivo correlations remains to be established. With the rapid and ongoing advancement of technology, there is much potential to better utilise computational methods including different types of modelling and simulation approaches to support inhaled formulation development. This review intends to provide an introduction on some fundamental concepts in pulmonary drug delivery and inhaled formulation development followed by discussions on some challenges and opportunities in the translation of inhaled pharmaceuticals from preclinical studies to clinical development. The review concludes with some recent advancements in modelling and simulation approaches that could play an increasingly important role in modern formulation development of inhaled pharmaceuticals.

Development of Dry Powder Inhaler Patient Interfaces for Improved Aerosol Delivery to Children

The objective of this study was to explore different internal flow passages in the patient interface region of a new air-jet-based dry powder inhaler (DPI) in order to minimize device and extrathoracic aerosol depositional losses using computational fluid dynamics (CFD) simulations. The best-performing flow passages were used for oral and nose-to-lung (N2L) aerosol delivery in pediatric extrathoracic airway geometries consistent with a 5-year-old child. Aerosol delivery conditions were based on a previously developed and tested air-jet DPI device and included a base flow rate of 13.3 LPM (delivered from a small ventilation bag) and an inhaled air volume of 750 mL. Initial CFD models of the system clearly established that deposition on either the back of the throat or nasal cannula bifurcation was strongly correlated with the maximum velocity exiting the flow passage. Of all designs tested, the combination of a 3D rod array and rapid expansion of the flow passage side walls was found to dramatically reduce interface and device deposition and improve lung delivery of the aerosol. For oral aerosol administration, the optimal flow passage compared with a base case reduced device, mouthpiece, and mouth-throat deposition efficiencies by factors of 8-, 3-, and 2-fold, respectively. For N2L aerosol administration, the optimal flow pathway compared with a base case reduced device, nasal cannula, and nose-throat deposition by 16-, 6-, and 1.3-fold, respectively. In conclusion, a new patient interface design including a 3D rod array and rapid expansion dramatically improved transmission efficiency of a dry powder aerosol.

Anthropometry-based generation of personalized and population-specific human airway models

Understanding aerosol deposition in the human lung is of great significance in pulmonary toxicology and inhalation pharmacology. Adverse effects of inhaled environmental aerosols and pharmacological efficacy of inhaled therapeutics are dependent on aerosol properties as well as person-specific respiratory tract anatomy and physiology. Anatomical geometry and physiological function of human airways depend on age, gender, weight, fitness, health, and disease status. Tools for the generation of the population- and subject-specific virtual airway anatomical geometry based on anthropometric data and physiological vitals are invaluable in respiratory diagnostics, personalized pulmonary pharmacology, and model-based management of chronic respiratory diseases. Here we present a novel protocol and software framework for the generation of subject-specific airways based on anthropometric measurements of the subject's body, using the anatomical input, and the conventional spirometry, providing the functional (physiological) data. This model can be used for subject-specific simulations of respiration physiology, gas exchange, and aerosol inhalation and deposition.

AVATREE: An open-source computational modelling framework modelling Anatomically Valid Airway TREE conformations

This paper presents AVATREE, a computational modelling framework that generates Anatomically Valid Airway tree conformations and provides capabilities for simulation of broncho-constriction apparent in obstructive pulmonary conditions. Such conformations are obtained from the personalized 3D geometry generated from computed tomography (CT) data through image segmentation. The patient-specific representation of the bronchial tree structure is extended beyond the visible airway generation depth using a knowledge-based technique built from morphometric studies. Additional functionalities of AVATREE include visualization of spatial probability maps for the airway generations projected on the CT imaging data, and visualization of the airway tree based on local structure properties. Furthermore, the proposed toolbox supports the simulation of broncho-constriction apparent in pulmonary diseases, such as chronic obstructive pulmonary disease (COPD) and asthma. AVATREE is provided as an open-source toolbox in C++ and is supported by a graphical user interface integrating the modelling functionalities. It can be exploited in studies of gas flow, gas mixing, ventilation patterns and particle deposition in the pulmonary system, with the aim to improve clinical decision making.

A multiscale absorption and transit model for oral drug delivery: Formulation and applications during fasting conditions

Most Food and Drug Administration (FDA)-approved drugs are administered orally, despite the complex process of oral drug absorption that is difficult to analyze experimentally. Oral bioavailability is dependent on the drug compound as well as the physiological and anatomical states of the user. Thus, computational models have emerged to mechanistically capture and predict the oral absorption process. The current models are generally 0D compartmental models and are limited by (a) simplified physiological characteristics of the gastrointestinal tract (GIT), (b) semiempirical/analytical dissolution profiles of the tested drugs, (c) incorrect absorption for some drug BCS classes (class IIa, for example), (d) GITs size variability among population, (e) incorrectly predicting the absorption of drugs that are GIT target specific, and (f) erroneous mixing in the domain. In this study, we have developed a multiscale absorption and transit (MAT) toolkit to simulate the dissolution, transport, absorption, distribution, metabolism, and elimination of orally administered drugs in the human GIT at multiple levels. MAT was constructed by integrating the spatially accurate first-principles driven high-fidelity drug transport, dissolution, and absorption model in the human stomach and GIT using our recently published quasi-3D (Q3D) framework. The process integrated the multilayer intestine physiologically based pharmacokinetics models with the whole-body compartmental models to predict the systemic pharmacokinetics of oral drugs. The computational results showed that this multiscale tool was able to match the experimental concentration results (individual and population) better than the traditional compartmental models. Ultimately, MAT will be developed into a commercial product to meet urgent demands from pharmaceutical and biomedical industries.

Prediction of Aerosol Deposition in the Human Respiratory Tract via Computational Models: A Review with Recent Updates

The measurement of deposited aerosol particles in the respiratory tract via in vivo and in vitro approaches is difficult due to those approaches’ many limitations. In order to overcome these obstacles, different computational models have been developed to predict the deposition of aerosol particles inside the lung. Recently, some remarkable models have been developed based on conventional semi-empirical models, one-dimensional whole-lung models, three-dimensional computational fluid dynamics models, and artificial neural networks for the prediction of aerosol-particle deposition with a high accuracy relative to experimental data. However, these models still have some disadvantages that should be overcome shortly. In this paper, we take a closer look at the current research trends as well as the future directions of this research area.

A Review of Respiratory Anatomical Development, Air Flow Characterization and Particle Deposition

The understanding of complex inhalation and transport processes of pollutant particles through the human respiratory system is important for investigations into dosimetry and respiratory health effects in various settings, such as environmental or occupational health. The studies over the last few decades for micro- and nanoparticle transport and deposition have advanced the understanding of drug-aerosol impacts in the mouth-throat and the upper airways. However, most of the Lagrangian and Eulerian studies have utilized the non-realistic symmetric anatomical model for airflow and particle deposition predictions. Recent improvements to visualization techniques using high-resolution computed tomography (CT) data and the resultant development of three dimensional (3-D) anatomical models support the realistic representation of lung geometry. Yet, the selection of different modelling approaches to analyze the transitional flow behavior and the use of different inlet and outlet conditions provide a dissimilar prediction of particle deposition in the human lung. Moreover, incorporation of relevant physical and appropriate boundary conditions are important factors to consider for the more accurate prediction of transitional flow and particle transport in human lung. This review critically appraises currently available literature on airflow and particle transport mechanism in the lungs, as well as numerical simulations with the aim to explore processes involved. Numerical studies found that both the Euler–Lagrange (E-L) and Euler–Euler methods do not influence nanoparticle (particle diameter ≤50 nm) deposition patterns at a flow rate ≤25 L/min. Furthermore, numerical studies demonstrated that turbulence dispersion does not significantly affect nanoparticle deposition patterns. This critical review aims to develop the field and increase the state-of-the-art in human lung modelling.

Devices for Improved Delivery of Nebulized Pharmaceutical Aerosols to the Lungs

Nebulizers have a number of advantages for the delivery of inhaled pharmaceutical aerosols, including the use of aqueous formulations and the ability to deliver process-sensitive proteins, peptides, and biological medications. A frequent disadvantage of nebulized aerosols is poor lung delivery efficiency, which wastes valuable medications, increases delivery times, and may increase side effects of the medication. A focus of previous development efforts and previous nebulizer reviews, has been an improvement of the underlying nebulization technology controlling the breakup of a liquid into droplets. However, for a given nebulization technology, a wide range of secondary devices and strategies can be implemented to significantly improve lung delivery efficiency of the aerosol. This review focuses on secondary devices and technologies that can be implemented to improve the lung delivery efficiency of nebulized aerosols and potentially target the region of drug delivery within the lungs. These secondary devices may (1) modify the aerosol size distribution, (2) synchronize aerosol delivery with inhalation, (3) reduce system depositional losses at connection points, (4) improve the patient interface, or (5) guide patient inhalation. The development of these devices and technologies is also discussed, which often includes the use of computational fluid dynamic simulations, three-dimensional printing and rapid prototype device and airway model construction, realistic in vitro experiments, and in vivo analysis. Of the devices reviewed, the implementation of streamlined components may be the most direct and lowest cost approach to enhance aerosol delivery efficiency within nonambulatory nebulizer systems. For applications involving high-dose medications or precise dose administration, the inclusion of active devices to control aerosol size, guide inhalation, and synchronize delivery with inhalation hold considerable promise.

Biological Obstacles for Identifying In Vitro-In Vivo Correlations of Orally Inhaled Formulations

Oral inhalation of drugs is the classic therapy of obstructive lung diseases. In contrast to the oral route, the link between in vitro and in vivo findings is less well defined and predictive models and parameters for in vitro-in vivo correlations are missing. Frequently used in vitro models and problems in obtaining in vivo values to establish such models and to identify the action of formulations in vivo are discussed. It may be concluded that major obstacles to link in vitro parameters on in vivo action include lack of treatment adherence and incorrect use of inhalers by patients, variation in inhaler performance, changes by humidity, uncertainties about lung deposition, and difficulties to measure drug levels in epithelial lining fluid and tissue. Physiologically more relevant in vitro models, improvement in inhaler performance, and better techniques for in vivo measurements may help to better understand importance and interactions between individual in vitro parameters in pulmonary delivery.

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.

High-Efficiency Nose-to-Lung Aerosol Delivery in an Infant: Development of a Validated Computational Fluid Dynamics Method

Background: Computational fluid dynamics (CFD) provides a powerful tool for developing new high-efficiency aerosol delivery strategies, such as nose-to-lung (N2L) aerosol administration to infants and children using correctly sized aerosols. The objective of this study was to establish numerically efficient CFD solution methods and guidelines for simulating N2L aerosol administration to an infant based on comparisons with concurrent in vitro experiments. Materials and Methods: N2L administration of a micrometer-sized aerosol (mass median aerodynamic diameter [MMAD] = 1.4 μm) was evaluated using concurrent CFD simulations and in vitro experiments. Aerosol transport and deposition was assessed in a new nasal airway geometry of a 6-month-old infant with a streamlined nasal cannula interface, which was constructed as a CFD mesh and three-dimensionally printed to form an identical physical prototype. CFD meshes explored were a conventional tetrahedral approach with near-wall (NW) prism elements and a new polyhedral mesh style with an equally refined NW layer. The presence of turbulence in the model was evaluated using a highly efficient low-Reynolds number (LRN) k-ω turbulence model, with previously established NW corrections that accounted for anisotropic wall-normal turbulence as well as improved NW velocity interpolations and hydrodynamic particle damping. Results: Use of the new polyhedral mesh was found to improve numerical efficiency by providing more rapid convergence and requiring fewer control volumes. Turbulent flow was found in the nasal geometry, generated by the inlet jets from the nasal cannula interface. However, due to the small particle size, turbulent dispersion was shown to have little effect on deposition. Good agreement was established between the CFD predictions using the numerically efficient LRN k-ω model with appropriate NW corrections and in vitro deposition data. Aerosol transmission efficiencies through the delivery tube, nasal cannula, and infant nasal model, based on experimental and CFD predictions, were 93.0% and 91.5%, respectively. Conclusions: A numerically efficient CFD approach was established to develop transnasal aerosol administration to infants and children. Small particle aerosols with aerodynamic diameters of ∼1.5 μm were confirmed to have low inertial depositional loss, and have low deposition from turbulent dispersion, making them ideal for high-efficiency lung delivery through an infant nasal cannula interface.