A Partial Differential Equation Approach to Inhalation Physiologically Based Pharmacokinetic Modeling

Assessment of Epithelial Lining Fluid Partitioning of Systemically Administered Monoclonal Antibodies in Rats

For systemically administered monoclonal antibodies (mAbs) with pharmacological targets in the epithelial lining fluid (ELF), information on the partitioning of mAb between plasma and ELF is instrumental for dose predictions. Bronchoalveolar lavage (BAL) combined with measurements of urea as indicator of sample dilution is often used to estimate ELF concentrations of a drug. However, unbalanced extraction of mAb and urea could potentially lead to a systematic bias in the back-calculated ELF concentration. In the present study 0.5, 1, or 4 mL phosphate-buffered saline was instilled to lungs of rats to obtain lavage samples after systemic dosing of mAb and tool small molecule (n≥4/group). Furthermore, extraction of urea, mAb and the small molecule was assessed by repeatedly lavaging the lung (n = 4). There was no statistically significant difference in the calculated partitioning into ELF between the evaluated instillation volumes. Repeated BAL demonstrated that urea and the small molecule were extracted from other sources than the ELF. In contrast, there was limited to none in-flow of mAb into the lavage fluid. The unbalanced extraction of urea and mAb could theoretically result in underestimated ELF concentrations and the calculated partitioning of 0.17±0.062 might therefore constitute a lower boundary for the true partitioning.

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.

Opportunities and challenges of physiologically based pharmacokinetic modeling in drug delivery

Physiologically based pharmacokinetic (PBPK) modeling is an important in silico tool to bridge drug properties and in vivo PK behaviors during drug development. Over the recent decade, the PBPK method has been largely applied to drug delivery systems (DDS), including oral, inhaled, transdermal, ophthalmic, and complex injectable products. The related therapeutic agents have included small-molecule drugs, therapeutic proteins, nucleic acids, and even cells. Simulation results have provided important insights into PK behaviors of new dosage forms, which strongly support drug regulation. In this review, we comprehensively summarize recent progress in PBPK applications in drug delivery, which shows large opportunities for facilitating drug development. In addition, we discuss the challenges of applying this methodology from a practical viewpoint.

A mechanistic framework for a priori pharmacokinetic predictions of orally inhaled drugs

The fate of orally inhaled drugs is determined by pulmonary pharmacokinetic processes such as particle deposition, pulmonary drug dissolution, and mucociliary clearance. Even though each single process has been systematically investigated, a quantitative understanding on the interaction of processes remains limited and therefore identifying optimal drug and formulation characteristics for orally inhaled drugs is still challenging. To investigate this complex interplay, the pulmonary processes can be integrated into mathematical models. However, existing modeling attempts considerably simplify these processes or are not systematically evaluated against (clinical) data. In this work, we developed a mathematical framework based on physiologically-structured population equations to integrate all relevant pulmonary processes mechanistically. A tailored numerical resolution strategy was chosen and the mechanistic model was evaluated systematically against data from different clinical studies. Without adapting the mechanistic model or estimating kinetic parameters based on individual study data, the developed model was able to predict simultaneously (i) lung retention profiles of inhaled insoluble particles, (ii) particle size-dependent pharmacokinetics of inhaled monodisperse particles, (iii) pharmacokinetic differences between inhaled fluticasone propionate and budesonide, as well as (iv) pharmacokinetic differences between healthy volunteers and asthmatic patients. Finally, to identify the most impactful optimization criteria for orally inhaled drugs, the developed mechanistic model was applied to investigate the impact of input parameters on both the pulmonary and systemic exposure. Interestingly, the solubility of the inhaled drug did not have any relevant impact on the local and systemic pharmacokinetics. Instead, the pulmonary dissolution rate, the particle size, the tissue affinity, and the systemic clearance were the most impactful potential optimization parameters. In the future, the developed prediction framework should be considered a powerful tool for identifying optimal drug and formulation characteristics.

Drug Interactions for Low-Dose Inhaled Nemiralisib: A Case Study Integrating Modeling, In Vitro, and Clinical Investigations

In vitro data for low-dose inhaled phosphoinositide 3-kinase delta inhibitor nemiralisib revealed that it was a substrate and a potent metabolism-dependent inhibitor of cytochrome P450 (P450) 3A4 and a P-glycoprotein (P-gp) substrate. An integrated in silico, in vitro, and clinical approach including a clinical drug interaction study as well as a bespoke physiologically based pharmacokinetic (PBPK) model was used to assess the drug-drug interaction (DDI) risk. Inhaled nemiralisib (100 µg, single dose) was coadministered with itraconazole, a potent P4503A4/P-gp inhibitor, following 200 mg daily administrations for 10 days in 20 male healthy subjects. Systemic exposure to nemiralisib (AUC_0-inf) increased by 2.01-fold versus nemiralisib alone. To extrapolate the clinical data to other P4503A4 inhibitors, an inhaled PBPK model was developed using Simcyp software. Retrospective simulation of the victim risk showed good agreement between simulated and observed data (AUC_0-inf ratio 2.3 vs. 2.01, respectively). Prospective DDI simulations predicted a weak but manageable drug interaction when nemiralisib was coadministered with other P4503A4 inhibitors, such as the macrolides clarithromycin and erythromycin (simulated AUC_0-inf ratio of 1.7), both common comedications in the intended patient populations. PBPK and static mechanistic models were also used to predict a negligible perpetrator DDI effect for nemiralisib on other P4503A4 substrates, including midazolam (a sensitive probe substrate of P4503A4) and theophylline (a narrow therapeutic index drug and another common comedication). In summary, an integrated in silico, in vitro, and clinical approach including an inhalation PBPK model has successfully discharged any potential patient DDI risks in future nemiralisib clinical trials. SIGNIFICANCE STATEMENT: This paper describes the integration of in silico, in vitro, and clinical data to successfully discharge potential drug-drug interaction risks for a low-dose inhaled drug. This work featured assessment of victim and perpetrator risks of drug transporters and cytochrome P450 enzymes, utilizing empirical and mechanistic approaches combined with clinical data (drug interaction and human absorption, metabolism, and pharmacokinetics) and physiologically based pharmacokinetic modeling approaches to facilitate bespoke risk assessment in target patient populations.

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.