Development of a hybrid CFD-PBPK model to predict the transport of xenon gas around a human respiratory system to systemic regions

Reconstruction of exposure to methylene diphenyl-4,4'-diisocyanate (MDI) aerosol using computational fluid dynamics, physiologically based toxicokinetics and statistical modeling

OBJECTIVES

This study employed computational fluid dynamics (CFD), physiologically based toxicokinetics (PBTK), and statistical modeling to reconstruct exposure to methylene diphenyl-4,4'-diisocyanate (MDI) aerosol. By utilizing a validated CFD model, human respiratory deposition of MDI aerosol in different workload conditions was investigated, while a PBTK model was calibrated using experimental rat data. Biomonitoring data and Markov Chain Monte Carlo (MCMC) simulation were utilized for exposure assessment.

RESULTS

Deposition fraction of MDI in the respiratory tract at the light, moderate, and heavy activity were 0.038, 0.079, and 0.153, respectively. Converged MCMC results as the posterior means and prior values were obtained for several PBTK model parameters. In our study, we calibrated a rat model to investigate the transport, absorption, and elimination of 4,4'-MDI via inhalation exposure. The calibration process successfully captured experimental data in the lungs, liver, blood, and kidneys, allowing for a reasonable representation of MDI distribution within the rat model. Our calibrated model also represents MDI dynamics in the bloodstream, facilitating the assessment of bioavailability. For human exposure, we validated the model for recent and long-term MDI exposure using data from relevant studies.

CONCLUSION

Our computational models provide reasonable insights into MDI exposure, contributing to informed risk assessment and the development of effective exposure reduction strategies.

Deconvolution of Systemic Pharmacokinetics Predicts Inhaled Aerosol Dosimetry of Nicotine

Absorption of inhaled compounds can occur from multiple sites based on upper and lower respiratory tract deposition, and clearance mechanisms leading to differential local and systemic pharmacokinetics. Deriving inhaled aerosol dosimetry and local tissue concentrations for nose-only exposure in rodents and inhaled products in humans is challenging. In this study we use inhaled nicotine as an example to identify regional respiratory tract deposition, absorption fractions, and their contribution toward systemic pharmacokinetics in rodents and humans. A physiologically based pharmacokinetic (PBPK) model was constructed to describe the disposition of nicotine and its major metabolite, cotinine. The model description for the lungs was simplified to include an upper respiratory tract region with active mucociliary clearance and a lower respiratory tract region. The PBPK model parameters such as rate of oral absorption, metabolism and clearance were fitted to the published nicotine and cotinine plasma concentrations post systemic administration and oral dosing. The fractional deposition of inhaled aerosol in the upper and lower respiratory tract regions was estimated by fitting the plasma concentrations. The model predicted upper respiratory tract deposition was 63.9% for nose-only exposure to nicotine containing nebulized aqueous aerosol in rats and 60.2% for orally inhaled electronic vapor product in humans. A marked absorption of nicotine from the upper respiratory tract and the gastrointestinal tract for inhaled aqueous aerosol contributed to the differential systemic pharmacokinetics in rats and humans. The PBPK model derived dosimetry shows that the current aerosol dosimetry models with their posteriori application using independent aerosol physicochemical characterization to predict aerosol deposition are insufficient and will need to consider complex interplay of inhaled aerosol evolutionary process. While the study highlights the needs for future research, it provides a preliminary framework for interpreting pharmacokinetics of inhaled aerosols to facilitate the analysis of in vivo exposure-responses for pharmacological and toxicological assessments.

Airborne transmission of biological agents within the indoor built environment: a multidisciplinary review

The nature and airborne dispersion of the underestimated biological agents, monitoring, analysis and transmission among the human occupants into building environment is a major challenge of today. Those agents play a crucial role in ensuring comfortable, healthy and risk-free conditions into indoor working and leaving spaces. It is known that ventilation systems influence strongly the transmission of indoor air pollutants, with scarce information although to have been reported for biological agents until 2019. The biological agents' source release and the trajectory of airborne transmission are both important in terms of optimising the design of the heating, ventilation and air conditioning systems of the future. In addition, modelling via computational fluid dynamics (CFD) will become a more valuable tool in foreseeing risks and tackle hazards when pollutants and biological agents released into closed spaces. Promising results on the prediction of their dispersion routes and concentration levels, as well as the selection of the appropriate ventilation strategy, provide crucial information on risk minimisation of the airborne transmission among humans. Under this context, the present multidisciplinary review considers four interrelated aspects of the dispersion of biological agents in closed spaces, (a) the nature and airborne transmission route of the examined agents, (b) the biological origin and health effects of the major microbial pathogens on the human respiratory system, (c) the role of heating, ventilation and air-conditioning systems in the airborne transmission and (d) the associated computer modelling approaches. This adopted methodology allows the discussion of the existing findings, on-going research, identification of the main research gaps and future directions from a multidisciplinary point of view which will be helpful for substantial innovations in the field.

Predicting Regional Respiratory Tissue and Systemic Concentrations of Orally Inhaled Drugs through a Novel PBPK Model

Oral inhalation (OI) of drugs is the route of choice to treat respiratory diseases or for recreational drug use (e.g., cannabis). After OI, the drug is deposited in and systemically absorbed from various regions of the respiratory tract. Measuring regional respiratory tissue drug concentrations at the site of action is important for evaluating the efficacy and safety of orally inhaled drugs (OIDs). Because such a measurement is routinely not possible in humans, the only alternative is to predict these concentrations, for example by physiologically based pharmacokinetic (PBPK) modeling. Therefore, we developed an OI-PBPK model to integrate the interplay between regional respiratory drug deposition and systemic absorption to predict regional respiratory tissue and systemic drug concentrations. We validated our OI-PBPK model by comparing the simulated and observed plasma concentration-time profiles of two OIDs, morphine and nicotine. Furthermore, we performed sensitivity analyses to quantitatively demonstrate the impact of key parameters on the extent and pattern of regional respiratory drug deposition, absorption, and the resulting regional respiratory tissue and systemic plasma concentrations. Our OI-PBPK model can be applied to predict regional respiratory tissue and systemic drug concentrations to optimize OID formulations, delivery systems, and dosing regimens. Furthermore, our model could be used to establish the bioequivalence of generic OIDs for which systemic plasma concentrations are not measurable or are not a good surrogate of the respiratory tissue drug concentrations. SIGNIFICANCE STATEMENT: Our OI-PBPK model is the first comprehensive model to predict regional respiratory deposition, as well as systemic and regional tissue concentrations of OIDs, especially at the drug's site of action, which is difficult to measure in humans. This model will help optimize OID formulations, delivery systems, dosing regimens, and bioequivalence assessment of generic OID. Furthermore, this model can be linked with organs-on-chips, pharmacodynamic and quantitative systems pharmacology models to predict and evaluate the safety and efficacy of OID.

TRANSITION AND LAMINAR FLOWS IN A REALISTIC GEOMETRY OF HUMAN UPPER AIRWAY

In this study, a realistic respiratory airway model extending from oral to the end of the trachea including all the key details of the passage was produced. A series of CT scan images were used to generate the topological data of airway cross-sections that were used to generate the computational model, as well as the three-dimensional (3D) printed model of the passage for experimental study. The airflow velocity field and pressure drop in the airway for different breathing rates of 5, 7.5, 10, and 12.5[Formula: see text]L/min were investigated numerically (by laminar and transition models) and experimentally. The velocity distributions, pressure variation, and streamlines along the oral–trachea airway model were studied. The maximum pressure drop was shown to occur in the narrowest part of the larynx region. It was also concluded that the laryngeal jet could significantly influence the airway flow patterns in the trachea. A comparison between the numerical results and experimental data showed that the transition [Formula: see text]–kl–[Formula: see text] model can give better predictions of pressure losses, especially for flow rates higher than 10[Formula: see text]L/min. The simulation results for the velocity profiles in the trachea were also compared with the available particle image velocimetry (PIV) data and earlier simulations. Despite inter-personal variability and difference in the flow regime, the qualitative agreement was found.

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

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

Investigation of airflow at different activity conditions in a realistic model of human upper respiratory tract

In the present study, the turbulent flows inside a realistic model of the upper respiratory tract were investigated numerically and experimentally. The airway model included the geometrical details of the oral cavity to the end of the trachea that was based on a series of CT-scan images. The topological data of the respiratory tract were used for generating the computational model as well as the 3D-printed model that was used in the experimental pressure drop measurement. Different airflow rates of 30, 45, and 60 L/min, which correspond to the light, semi-light, and heavy activity breathing conditions, were investigated numerically using turbulence and transition models, as well as experimentally. Simulation results for airflow properties, including velocity vectors, pressure drops, streamlines, eddy viscosity, and turbulent kinetic energy contours in the oral-trachea airway model, were presented. The simulated pressure drop was compared with the experimental data, and reasonable agreement was found. The obtained results showed that the maximum pressure drop occurs in the narrowest part of the larynx region. A comparison between the numerical results and experimental data showed that the transition (γ-Reθ) SST model predicts higher pressure losses, especially at higher breathing rates. Formations of the secondary flows in the oropharynx and trachea regions were also observed. In addition, the simulation results showed that in the trachea region, the secondary flow structures dissipated faster for the flow rate of 60 L/min compared to the lower breathing rates of 30 and 45 L/min.

Bridging inhaled aerosol dosimetry to physiologically based pharmacokinetic modeling for toxicological assessment: nicotine delivery systems and beyond

One of the challenges for toxicological assessment of inhaled aerosols is to accurately predict their deposited and absorbed dose. Transport, evolution, and deposition of liquid aerosols are driven by complex processes dominated by convection-diffusion that depend on various factors related to physics and chemistry. These factors include the physicochemical properties of the pure substance of interest and associated mixtures, the physical and chemical properties of the aerosols generated, the interplay between different factors during transportation and deposition, and the subject-specific inhalation topography. Several inhalation-based physiologically based pharmacokinetic (PBPK) models have been developed, but the applicability of these models for aerosols has yet to be verified. Nicotine is among several substances that are often delivered via the pulmonary route, with varied kinetics depending upon the route of exposure. This was used as an opportunity to review and discuss the current knowledge and state-of-the-art tools combining aerosol dosimetry predictions with PBPK modeling efforts. A validated tool could then be used to perform for toxicological assessment of other inhaled therapeutic substances. The Science Panel from the Alliance of Risk Assessment have convened at the "Beyond Science and Decisions: From Problem Formulation to Dose-Response Assessment" workshop to evaluate modeling approaches and address derivation of exposure-internal dose estimations for inhaled aerosols containing nicotine or other substances. The discussion involved PBPK model evaluation criteria, challenges, and choices that arise in such a model design, development, and application as a computational tool for use in human toxicological assessments.