Simulated pharmacokinetics of inhaled caffeine and melatonin from existing products indicate the lack of dosimetric considerations

Numerous commercially available inhalable products claim to improve sleep-wake cycle-related target indications by delivering a wide variety of chemicals like caffeine and melatonin. The resulting exposure-responses from inhaling different doses are unknown and obtaining early understanding of resulting pharmacokinetics is beneficial. This study applied a physiologically based pharmacokinetic modeling approach to predict the inhalation pharmacokinetics of caffeine and melatonin for different target indications related to the sleep-wake cycle. The model predicted rapid systemic delivery of caffeine and melatonin based on airway regional deposition of inhaled aerosol. A low inhaled dose of 1 mg of caffeine resulted in a 72.3-times lower plasma maximal concentration and was predicted to not improve cognitive performance task outcomes compared to oral consumption of coffee containing 80 mg of caffeine. Conversely, 2-mg oral and inhaled doses of melatonin under recommended directions of use result in more than 25.1- and 645-times higher plasma concentrations compared to endogenous melatonin, respectively. The recommended doses for inhalation products for potential improvement in the target indications vary widely. Additional research is needed to evaluate the human pharmacokinetics, efficacy, and safety of inhaled products. Given the lack of assessments, inhaled caffeine and melatonin must be consumed with caution as the toxicological concerns are not known and could outweigh the potential beneficial effects.

Aerosol delivery and spatiotemporal tissue distribution of hydroxychloroquine in rat lung

Inhalation enables the delivery of drugs directly to the lung, increasing the retention for prolonged exposure and maximizing the therapeutic index. However, the differential regional lung exposure kinetics and systemic pharmacokinetics are not fully known, and their estimation is critical for pulmonary drug delivery. The study evaluates the pharmacokinetics of hydroxychloroquine in different regions of the respiratory tract for multiple routes of administration. We also evaluated the influence of different inhaled formulations on systemic and lung pharmacokinetics by identifying suitable nebulizers followed by early characterization of emitted aerosol physicochemical properties. The salt- and freebase-based formulations required different nebulizers and generated aerosol with different physicochemical properties. An administration of hydroxychloroquine by different routes resulted in varied systemic and lung pharmacokinetics, with oral administration resulting in low tissue concentrations in all regions of the respiratory tract. A nose-only inhalation exposure resulted in higher and sustained lung concentrations of hydroxychloroquine with a lung parenchyma-to-blood ratio of 386 after 1440 min post-exposure. The concentrations of hydroxychloroquine in different regions of the respiratory tract (i.e., nasal epithelium, larynx, trachea, bronchi, and lung parenchyma) varied over time, indicating different retention kinetics. The spatiotemporal distribution of hydroxychloroquine in the lung is different due to the heterogeneity of cell types, varying blood perfusion rate, clearance mechanisms, and deposition of inhaled aerosol along the respiratory tract. In addition to highlighting the varied lung physiology, these results demonstrate the ability of the lung to retain increased levels of inhaled lysosomotropic drugs. Such findings are critical for the development of future inhalation-based therapeutics, aiming to optimize target site exposure, enable precision medicine, and ultimately enhance clinical outcomes.

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.

Use of micro-CT to determine tracheobronchial airway geometries in three strains of mice used in inhalation toxicology as disease models

Aerosol dosimetry estimates for mouse strains used as models for human disease are not available, primarily because of the lack of tracheobronchial airway morphometry data. By using micro-CT scans of in-situ prepared lung casts, tracheobronchial airway morphometry for four strains of mice were obtained: Balb/c, AJ, C57BL/6, and Apoe^-/- . The automated tracheobronchial airway morphometry algorithms for airway length and diameter were successfully verified against previously published manual and automated tracheobronchial airway morphometry data derived from two identical in-situ Balb/c mouse lung casts. There was also excellent agreement in tracheobronchial airway length and diameter between the automated and manual airway data for the AJ, C57BL/6, and Apoe^-/- mice. Differences in branch angle measurements were partially due to the differences in definition between the automated algorithms and manual morphometry techniques. Unlike the manual airway morphometry techniques, the automated algorithms were able to provide a value for inclination to gravity for each airway. Inclusion of an inclination to gravity angle for each airway along with airway length, diameter, and branch angle make the current automated tracheobronchial airway data suitable for use in dosimetry programs that can provide dosimetry estimates for inhaled material. The significant differences in upper tracheobronchial airways between Balb/c mice and between C57BL/6 and Apoe^-/- mice highlight the need for mouse strain-specific aerosol dosimetry estimates.

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.

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.

State-of-the-art methods and devices for generation, exposure, and collection of aerosols from e-vapor products

E-vapor products (EVP) have become popular alternatives for cigarette smokers who would otherwise continue to smoke. EVP research is challenging and complex, mostly because of the numerous and rapidly evolving technologies and designs as well as the multiplicity of e-liquid flavors and solvents available on the market. There is an urgent need to standardize all stages of EVP assessment, from the production of a reference product to e-vapor generation methods and from physicochemical characterization methods to nonclinical and clinical exposure studies. The objective of this review is to provide a detailed description of selected experimental setups and methods for EVP aerosol generation and collection and exposure systems for their in vitro and in vivo assessment. The focus is on the specificities of the product that constitute challenges and require development of ad hoc assessment frameworks, equipment, and methods. In so doing, this review aims to support further studies, objective evaluation, comparison, and verification of existing evidence, and, ultimately, formulation of standardized methods for testing EVPs.

Physicochemical studies of direct interactions between lung surfactant and components of electronic cigarettes liquid mixtures

Direct physicochemical interactions between the major components of electronic cigarette liquids (e-liquids): glycerol (VG) and propylene glycol (PG), and lung surfactant (LS) were studied by determining the dynamic surface tension under a simulated breathing cycle using drop shape method. The studies were performed for a wide range of concentrations based on estimated doses of e-liquid aerosols (up to 2500 × the expected nominal concentrations) and for various VG/PG ratios. The results are discussed as relationships among mean surface tension, surface tension amplitude, and surface rheological properties (dilatational elasticity and viscosity) versus concentration and composition of e-liquid. The results showed that high local concentrations (>200 × higher than the estimated average dose after a single puffing session) may induce measurable changes in biophysical activity of LS; however, only ultra-high e-liquid concentrations inactivated the surfactant. Physiochemical characterization of e-liquids provide additional insights for the safety assessment of electronic nicotine delivery systems (ENDS).

Nonlocal modulation of the energy cascade in broadband-forced turbulence

Classically, large-scale forced turbulence is characterized by a transfer of energy from large to small scales via nonlinear interactions. We have investigated the changes in this energy transfer process in broadband forced turbulence where an additional perturbation of flow at smaller scales is introduced. The modulation of the energy dynamics via the introduction of forcing at smaller scales occurs not only in the forced region but also in a broad range of length scales outside the forced bands due to nonlocal triad interactions. Broadband forcing changes the energy distribution and energy transfer function in a characteristic manner leading to a significant modulation of the turbulence. We studied the changes in this transfer of energy when changing the strength and location of the small-scale forcing support. The energy content in the larger scales was observed to decrease, while the energy transport power for scales in between the large and small scale forcing regions was enhanced. This was investigated further in terms of the detailed transfer function between the triad contributions and observing the long-time statistics of the flow. The energy is transferred toward smaller scales not only by wave numbers of similar size as in the case of large-scale forced turbulence, but by a much wider extent of scales that can be externally controlled.