Current Scientific Considerations to Verify Physiologically-Based Pharmacokinetic Models and Their Implications for Locally Acting Products

Clinical Ocular Exposure Extrapolation for a Complex Ophthalmic Suspension Using Physiologically Based Pharmacokinetic Modeling and Simulation

The development of generic ophthalmic drug products with complex formulations is challenging due to the complexity of the ocular system and a lack of sensitive testing to evaluate the interplay of its physiology with ophthalmic drugs. New methods are needed to facilitate the development of ophthalmic generic drug products. Ocular physiologically based pharmacokinetic (O-PBPK) models can provide insight into drug partitioning in eye tissues that are usually not accessible and/or are challenging to sample in humans. This study aims to demonstrate the utility of an ocular PBPK model to predict human exposure following the administration of ophthalmic suspension. Besifloxacin (Bes) suspension is presented as a case study. The O-PBPK model for Bes ophthalmic suspension (Besivance® 0.6%) accounts for nasolacrimal drainage, suspended particle dissolution in the tears, ocular absorption, and distribution in the rabbit eye. A topical controlled release formulation was used to integrate the effect of Durasite® on Bes ocular retention. The model was subsequently used to predict Bes exposure after its topical administration in humans. Drug-specific parameters were used as validated for rabbits. The physiological parameters were adjusted to match human ocular physiology. Simulated human ocular pharmacokinetic profiles were compared with the observed ocular tissue concentration data to assess the OCAT models' ability to predict human ocular exposure. The O-PBPK model simulations adequately described the observed concentrations in the eye tissues following the topical administration of Bes suspension in rabbits. After adjustment of physiological parameters to represent the human eye, the extrapolation of clinical ocular exposure following a single ocular administration of Bes suspension was successful.

Update on the advances and challenges in bioequivalence testing methods for complex topical generic products

Most of the government regulatory agencies, including the United States Food and Drug Administration and the European Medicine Agency, demand that the generic complex topical products prove pharmaceutical and bioequivalence. The evaluation of bioequivalence for complex topical dermatological formulations is a challenging task that requires careful consideration of several factors. Although comparative clinical studies are still considered the gold standard approach for establishing bioequivalence in most formulations, these studies can be costly and insensitive to detect formulation differences. Therefore, significant efforts have been made to develop and validate alternative approaches that demonstrate bioequivalence and expedite the availability of high-quality generic topical dermatological products. This article reviews the current methods for determining the bioequivalence of topical formulations in humans, with particular emphasis on recent advances in these methodologies. Most of the alternative methods are sensitive and reproducible, with the capability to ease the financial burden of comparative clinical studies within a short delivery time. The limitations associated with each technique are reviewed in detail.

Physiologically Based Biopharmaceutics Model (PBBM) of Minimally Absorbed Locally Acting Drugs in the Gastrointestinal Tract-Case Study: Tenapanor

A physiologically based biopharmaceutics model (PBBM) was developed to predict stool and urine sodium content in response to tenapanor administration in healthy subjects. Tenapanor is a minimally absorbed small molecule that inhibits the sodium/hydrogen isoform 3 exchanger (NHE3). It is used to treat irritable bowel syndrome with constipation (IBS-C). Its mode of action in the gastrointestinal tract reduces the uptake of sodium, resulting in an increase in water secretion in the intestinal lumen and accelerating intestinal transit time. The strategy employed was to perform drug-drug interaction (DDI) modelling between sodium and tenapanor, with sodium as the "victim" administered as part of daily food intake and tenapanor as the "perpetrator" altering sodium absorption. Food effect was modelled, including meal-induced NHE3 activity using sodium as an inducer by normalising the induction kinetics of butyrate to sodium equivalents. The presented model successfully predicted both urine and stool sodium content in response to tenapanor dosed in healthy subjects (within 1.25-fold error) and provided insight into the clinical observations of tenapanor dosing time relative to meal ingestion. The PBBM model was applied retrospectively to assess the impact of different forms of tenapanor (free base vs. HCl salt) on its pharmacodynamic (PD) effect. The developed modelling strategy can be effectively adopted to increase confidence in using PBBM models for the prediction of the in vivo behaviour of minimally absorbed, locally acting drugs in the gastrointestinal tract, when other approaches (e.g., biomarkers or PD data) are not available.

Applications of Modeling and Simulation Approaches in Support of Drug Product Development of Oral Dosage Forms and Locally Acting Drug Products: a Symposium Summary

The number of modeling and simulation applications, including physiologically based pharmacokinetic (PBPK) models, physiologically based biopharmaceutics modeling (PBBM), and empirical models, has been constantly increasing along with the regulatory acceptance of these methodologies. While aiming at minimizing unnecessary human testing, these methodologies are used today to support the development and approval of novel drug products and generics. Modeling approaches are leveraged today for assessing drug-drug interaction, informing dose adjustments in renally or hepatically impaired patients, perform dose selection in pediatrics and pregnant women and diseased populations, and conduct biopharmaceutics-related assessments such as establish clinically relevant specifications for drug products and achieve quality assurance throughout the product life cycle. In the generics space, PBPK analyses are utilized toward virtual bioequivalence assessments within the scope of alternative bioequivalence approaches, product-specific guidance development, and food effect assessments among others. Case studies highlighting the evolving and expanding role of modeling and simulation approaches within the biopharmaceutics space were presented at the symposium titled "Model Informed Drug Development (MIDD): Role in Dose Selection, Vulnerable Populations, and Biowaivers - Chemical Entities" and Prologue "PBPK/PBBM to inform the Bioequivalence Safe Space, Food Effects, and pH-mediated DDIs" at the American Association of Pharmaceutical Scientists (AAPS) PharmSci 360 Annual Meeting in Boston, MA, on October 16-19, 2022, and are summarized here.

Clinical Ocular Exposure Extrapolation for Ophthalmic Solutions Using PBPK Modeling and Simulation

BACKGROUND

The development of generic ophthalmic drug products is challenging due to the complexity of the ocular system, and a lack of sensitive testing to evaluate the interplay of physiology with ophthalmic formulations. While measurements of drug concentration at the site of action in humans are typically sparse, these measurements are more easily obtained in rabbits. The purpose of this study is to demonstrate the utility of an ocular physiologically based pharmacokinetic (PBPK) model for translation of ocular exposure from rabbit to human.

METHOD

The Ocular Compartmental Absorption and Transit (OCAT™) model within GastroPlus® v9.8.2 was used to build PBPK models for levofloxacin (Lev), moxifloxacin (Mox), and gatifloxacin (Gat) ophthalmic solutions. in the rabbit eye. The models were subsequently used to predict Lev, Mox, and Gat exposure after ocular solution administrations in humans. Drug-specific parameters were used as fitted and validated in the rabbit OCAT model. The physiological parameters were scaled to match human ocular physiology.

RESULTS

OCAT model simulations for rabbit well described the observed concentrations in the eye compartments following Lev, Mox, and Gat solution administrations of different doses and various administration schedules. The clinical ocular exposure following ocular administration of Lev, Mox, and Gat solutions at different doses and various administration schedules was well predicted.

CONCLUSION

Even though additional case studies for different types of active pharmaceutical ingredients (APIs) and formulations will be needed, the current study represents an important step in the validation of the extrapolation method to predict human ocular exposure for ophthalmic drug products using PBPK models.

Modelling and Simulation Approaches to Support Formulation Optimization, Clinical Development and Regulatory Assessment of the Topically Applied Formulations- Nimesulide Solution Gel Case Study

The objective of the study was to show how mechanistic modelling can be used to characterize the skin absorption of Nimesulide (NIM) in both in vitro systems and in vivo subjects. A basic PBPK model for oral absorption to characterize the systemic disposition of NIM and MPML MechDermA^TM models for in vitro permeation and in vivo, topical absorption was developed and verified using published data. The developed models utilize drug physicochemical properties, formulation attributes and physiology information either collected from literature and/or from Simcyp databases (systems' data). Following the verification of the PBPK models virtual bioequivalence (VBE) trials were performed both at systemic and local exposure levels (dermis concentrations) to compare these formulations. A parameter sensitivity analysis was conducted to understand the impact of vehicle-related attributes on IVPT (in vitro permeation test) data. The vehicle-stratum corneum lipids partition coefficient in the formulation layer (Kp_{sc_lip:vehicle}) was identified to be an appropriate parameter to take into account the differences in dermal absorption of marketed preparations based on the qualitative composition. Thus, this parameter was optimized for each marketed product based on the published in vitro data. After verification of the IVPT model, IVIVE was performed to assess the predictability of the model for studying the in vivo pharmacokinetics of NIM. The VBE analysis concluded that these formulations are bioequivalent at the level of systemic and local dermis exposure. To summarize, the study shows the use of modelling and simulation (M&S) tools to better understand the behaviour of formulations and their interaction with human physiology.

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.

Optimization of topical formulations using a combination of in vitro methods to quantify the transdermal passive diffusion of drugs

This paper describes a new approach to the early-stage optimization of topical products and selection of lead formulation candidates. It demonstrates the application of open flow microperfusion in vitro in conjunction with the Franz diffusion cell to compare time-resolved, 24-hour profiles of diclofenac passive diffusion through all skin layers (including the skin barrier, dermis, and subcutis) resulting from nine topical formulations of different composition. The technique was successfully validated for in vitro sampling of diclofenac in interstitial fluid. A multi-compartmental model integrating the two datasets was analyzed and revealed that the passive diffusion of diclofenac through the dermis and subcutis does not correlate with its diffusion through the skin barrier and cannot be predicted using Franz diffusion cell data alone. The combined application of the two techniques provides a new, convenient tool for product development and selection enabling the comparison of topical formulation candidates and their impact on drug delivery through all skin layers. This approach can also generate the experimental data required to improve the robustness of mechanistic PBPK models, and when combined with clinical sampling via open flow microperfusion - for the development of better in vivo-in vitro correlative models.

Development of PBPK model for intra-articular injection in human: methotrexate solution and rheumatoid arthritis case study

A physiologically based model describing the dissolution, diffusion, and transfer of drug from the intra-articular (IA) space to the plasma, was developed for GastroPlus® v9.8. The model is subdivided into compartments representing the synovial fluid, synovium, and cartilage. The synovium is broken up into two sublayers. The intimal layer acts as a diffusion barrier between the synovial fluid and the subintimal layer. The subintimal layer of the synovium has fenestrated capillaries that allow the free drug to be transported into systemic circulation. The articular cartilage is broken up into 10 diffusion sublayers as it is much thicker than the synovium. The cartilage acts as a depot tissue for the drug to diffuse into from synovial fluid. At later times, the drug will diffuse from the cartilage back into synovial fluid once a portion of the dose enters systemic circulation. In this study, a listing of all relevant details and equations for the model is presented. Methotrexate was chosen as a case study to show the application and utility of the model, based on the availability of intravenous (IV), oral (PO) and IA administration data in patients presenting rheumatoid arthritis (RA) symptoms. Systemic disposition of methotrexate in RA patients was described by compartmental pharmacokinetic (PK) model with PK parameters extracted using the PKPlus™ module in GastroPlus®. The systemic PK parameters were validated by simulating PO administration of methotrexate before being used for simulation of IA administration. For methotrexate, the concentrations of drug in the synovial fluid and plasma were well described after adjustments of physiological parameters to account for RA disease state, and with certain assumptions about binding and diffusion. The results indicate that the model can correctly describe PK profiles resulting from administration in the IA space, however, additional cases studies will be required to evaluate ability of the model to scale between species and/or doses.

Physiologically-Based Pharmacokinetic Modeling to Support Determination of Bioequivalence for Dermatological Drug Products: Scientific and Regulatory Considerations

Physiologically-based pharmacokinetic (PBPK) modeling and simulation provides mechanism-based predictions of the pharmacokinetics of an active ingredient following its administration in humans. Dermal PBPK models describe the skin permeation and disposition of the active ingredient following the application of a dermatological product on the skin of virtual healthy and diseased human subjects. These models take into account information on product quality attributes, physicochemical properties of the active ingredient and skin (patho)physiology, and their interplay with each other. Regulatory and product development decision makers can leverage these quantitative tools to identify factors impacting local and systemic exposure. In the realm of generic drug products, the number of US Food and Drug Administratioin (FDA) interactions that use dermal PBPK modeling to support alternative bioequivalence (BE) approaches is increasing. In this report, we share scientific considerations on the development, verification and validation (V&V), and application of PBPK models within the context of a virtual BE assessment for dermatological drug products. We discuss the challenges associated with model V&V for these drug products stemming from the fact that target-site active ingredient concentrations are typically not measurable. Additionally, there are no established relationships between local and systemic PK profiles, when the latter are quantifiable. To that end, we detail a multilevel model V&V approach involving validation for the model of the drug product of interest coupled with the overall assessment of the modeling platform in use while leveraging in vitro and in vivo data related to local and systemic bioavailability.

Physiologically-based pharmacokinetic model for alectinib, ruxolitinib, and panobinostat in the presence of cancer, renal impairment, and hepatic impairment

Renal (RIP) and hepatic (HIP) impairments are prevalent conditions in cancer patients. They can cause changes in gastric emptying time, albumin levels, hematocrit, glomerular filtration rate, hepatic functional volume, blood flow rates, and metabolic activity that can modify drug pharmacokinetics. Performing clinical studies in such populations has ethical and practical issues. Using predictive physiologically-based pharmacokinetic (PBPK) models in the evaluation of the PK of alectinib, ruxolitinib, and panobinostat exposures in the presence of cancer, RIP, and HIP can help in using optimal doses with lower toxicity in these populations. Verified PBPK models were customized under scrutiny to account for the pathophysiological changes induced in these diseases. The PBPK model-predicted plasma exposures in patients with different health conditions within average 2-fold error. The PBPK model predicted an area under the curve ratio (AUCR) of 1, and 1.8, for ruxolitinib and panobinostat, respectively, in the presence of severe RIP. On the other hand, the severe HIP was associated with AUCR of 1.4, 2.9, and 1.8 for alectinib, ruxolitinib, and panobinostat, respectively, in agreement with the observed AUCR. Moreover, the PBPK model predicted that alectinib therapeutic cerebrospinal fluid levels are achieved in patients with non-small cell lung cancer, moderate HIP, and severe HIP at 1-, 1.5-, and 1.8-fold that of healthy subjects. The customized PBPK models showed promising ethical alternatives for simulating clinical studies in patients with cancer, RIP, and HIP. More work is needed to quantify other pathophysiological changes induced by simultaneous affliction by cancer and RIP or HIP.

Physiologically-based pharmacokinetic modeling to support bioequivalence and approval of generic products: A case for diclofenac sodium topical gel, 1

Establishing bioequivalence (BE) for dermatological drug products by conducting comparative clinical end point studies can be costly and the studies may not be sufficiently sensitive to detect certain formulation differences. Quantitative methods and modeling, such as physiologically‐based pharmacokinetic (PBPK) modeling, can support alternative BE approaches with reduced or no human testing. To enable PBPK modeling for regulatory decision making, models should be sufficiently verified and validated (V&V) for the intended purpose. This report illustrates the US Food and Drug Administration (FDA) approval of a generic diclofenac sodium topical gel that was based on a totality of evidence, including qualitative and quantitative sameness and physical and structural similarity to the reference product, an in vivo BE study with PK end points, and, more importantly, for the purposes of this report, a virtual BE assessment leveraging dermal PBPK modeling and simulation instead of a comparative clinical end point study in patients. The modeling approach characterized the relationship between systemic (plasma) and local (skin and synovial fluid) diclofenac exposure and demonstrated BE between the generic and reference products at the presumed site of action. Based on the fit‐for‐purpose modeling principle, the V&V process involved assessing observed data of diclofenac concentrations in skin tissues and plasma, and the overall performance of the modeling platform for relevant products. Using this case as an example, this report provides current scientific considerations on good practices for model V&V and the establishment of BE for dermatological drug products when leveraging PBPK modeling and simulation for regulatory decision making.

Development and Qualification of a Physiologically Based Pharmacokinetic Model of Finasteride and Minoxidil Following Scalp Application

In this study, we aimed to develop and qualify a PBPK model for scalp application using two drugs with marked differences in physicochemical properties and PK profiles. The parameters related to scalp physiology, drug PK, and formulations were incorporated into a Multi-Phase and Multi-Layer (MPML) Mechanistic Dermal Absorption (MechDermA) model within the Simcyp® Simulator V17. The finasteride PBPK model was linked to its effect on dihydrotestosterone (DHT) levels in plasma and scalp using an indirect response model. Predicted PK (and PD for finasteride) profiles and parameters were compared against the clinically reported data and justified by visual predictive checks and two-fold error criteria for model verification. The PBPK/PD model for finasteride reasonably demonstrated an ability to predict its respective PK and PD profiles, and parameters following scalp application under various clinical scenarios. Using the same scalp physiological input parameters, the minoxidil PBPK model was then developed and satisfactorily qualified with independent clinical datasets. Collectively, these results suggested that the established PBPK model may have broader utility for other topical formulations intended for scalp application.

In Vitro Dissolution and in Silico Modeling Shortcuts in Bioequivalence Testing

Purpose: To review in vitro testing and simulation platforms that are in current use to predict in vivo performances of generic products as well as other situations to provide evidence for biowaiver and support drug formulations development. Methods: Pubmed and Google Scholar databases were used to review published literature over the past 10 years. The terms used were “simulation AND bioequivalence” and “modeling AND bioequivalence” in the title field of databases, followed by screening, and then reviewing. Results: A total of 22 research papers were reviewed. Computer simulation using software such as GastroPlus™, PK-Sim^® and SimCyp^® find applications in drug modeling. Considering the wide use of optimization for in silico predictions to fit observed data, a careful review of publications is required to validate the reliability of these platforms. For immediate release (IR) drug products belonging to the Biopharmaceutics Classification System (BCS) classes I and III, difference factor (ƒ₁) and similarity factor (ƒ₂) are calculated from the in vitro dissolution data of drug formulations to support biowaiver; however, this method can be more discriminatory and may not be useful for all dissolution profiles. Conclusions: Computer simulation platforms need to improve their mechanistic physiologically based pharmacokinetic (PBPK) modeling, and if prospectively validated within a small percentage of error from the observed clinical data, they can be valuable tools in bioequivalence (BE) testing and formulation development.