Development of Mechanistic In Vitro-In Vivo Extrapolation to Support Bioequivalence Assessment of Long-Acting Injectables.

Long-acting injectable (LAI) formulations provide sustained drug release over an extended period ranging from weeks to several months to improve efficacy, safety, and compliance. Nevertheless, many challenges arise in the development and regulatory assessment of LAI drug products due to a limited understanding of the tissue response to injected particles (e.g., inflammation) impacting in vivo performance. Mechanism-based in silico methods may support the understanding of LAI-physiology interactions. The objectives of this study were as follows: (1) to use a mechanistic modeling approach to delineate the in vivo performance of DepoSubQ Provera® and formulation variants in preclinical species; (2) to predict human exposure based on the knowledge gained from the animal model. The PBPK model evaluated different elements involved in LAI administration and showed that (1) the effective in vivo particle size is potentially larger than the measured in vitro particle size, which could be due to particle aggregation at the injection site, and (2) local inflammation is a key process at the injection site that results in a transient increase in depot volume. This work highlights how a mechanistic modeling approach can identify critical physiological events and product attributes that may affect the in vivo performance of LAIs.

Predicting Human Dermal Drug Concentrations Using PBPK Modeling and Simulation: Clobetasol Propionate Case Study.

Quantitative in silico tools may be leveraged to mechanistically predict the dermato-pharmacokinetics of compounds delivered from topical and transdermal formulations by integrating systems of rate equations that describe permeation through the formulation and layers of skin and pilo-sebaceous unit, and exchange with systemic circulation via local blood flow. Delivery of clobetasol-17 propionate (CP) from DermovateTM cream was simulated using the Transdermal Compartmental Absorption & Transit (TCATTM) Model in GastroPlus®. The cream was treated as an oil-in-water emulsion, with model input parameters estimated from publicly available information and quantitative structure-permeation relationships. From the ranges of values available for model input parameters, a set of parameters was selected by comparing model outputs to CP dermis concentration-time profiles measured by dermal open-flow microperfusion (Bodenlenz et al. Pharm Res. 33(9):2229-38, 2016). Predictions of unbound dermis CP concentrations were reasonably accurate with respect to time and skin depth. Parameter sensitivity analyses revealed considerable dependence of dermis CP concentration profiles on drug solubility in the emulsion, relatively less dependence on dispersed phase volume fraction and CP effective diffusivity in the continuous phase of the emulsion, and negligible dependence on dispersed phase droplet size. Effects of evaporative water loss from the cream and corticosteroid-induced vasoconstriction were also assessed. This work illustrates the applicability of computational modeling to predict sensitivity of dermato-pharmacokinetics to changes in thermodynamic and transport properties of a compound in a topical formulation, particularly in relation to rate-limiting steps in skin permeation. Where these properties can be related to formulation composition and processing, such a computational approach may support the design of topically applied formulations.

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.

PBPK Modeling Approach to Predict the Behavior of Drugs Cleared by Metabolism in Pregnant Subjects and Fetuses.

This study aimed to develop a physiologically based pharmacokinetic (PBPK) model that simulates metabolically cleared compounds' pharmacokinetics (PK) in pregnant subjects and fetuses. This model accounts for the differences in tissue sizes, blood flow rates, enzyme expression levels, plasma protein binding, and other physiological factors affecting the drugs' PK in both the pregnant woman and the fetus. The PBPKPlus™ module in GastroPlus® was used to model the PK of metoprolol, midazolam, and metronidazole for both non-pregnant and pregnant groups. For each of the three compounds, the model was first developed and validated against PK data in healthy non-pregnant volunteers and then applied to predict the PK in the pregnant groups. The model accurately described the PK in both the non-pregnant and pregnant groups and explained well the differences in the plasma concentration due to pregnancy. When available, the fetal plasma concentration, placenta, and fetal tissue concentrations were also predicted reasonably well at different stages of pregnancy. The work described the use of a PBPK approach for drug development and demonstrates the ability to predict differences in PK in pregnant subjects and fetal exposure for metabolically cleared compounds.

Mechanistic modeling of ophthalmic, nasal, injectable, and implant generic drug products: A workshop summary report.

For approval, a proposed generic drug product must demonstrate it is bioequivalent (BE) to the reference listed drug product. For locally acting drug products, conventional BE approaches may not be feasible because measurements in local tissues at the sites of action are often impractical, unethical, or cost-prohibitive. Mechanistic modeling approaches, such as physiologically-based pharmacokinetic (PBPK) modeling, may integrate information from drug product properties and human physiology to predict drug concentrations in these local tissues. This may allow clinical relevance determination of critical drug product attributes for BE assessment during the development of generic drug products. In this regard, the Office of Generic Drugs of the US Food and Drug Administration has recently established scientific research programs to accelerate the development and assessment of generic products by utilizing model-integrated alternative BE approaches. This report summarizes the presentations and panel discussion from a public workshop that provided research updates and information on the current state of the use of PBPK modeling approaches to support generic product development for ophthalmic, injectable, nasal, and implant drug products.

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.

Developing Clinically Relevant Dissolution Specifications (CRDSs) for Oral Drug Products: Virtual Webinar Series.

A webinar series that was organised by the Academy of Pharmaceutical Sciences Biopharmaceutics focus group in 2021 focused on the challenges of developing clinically relevant dissolution specifications (CRDSs) for oral drug products. Industrial scientists, together with regulatory and academic scientists, came together through a series of six webinars, to discuss progress in the field, emerging trends, and areas for continued collaboration and harmonisation. Each webinar also hosted a Q&A session where participants could discuss the shared topic and information. Although it was clear from the presentations and Q&A sessions that we continue to make progress in the field of CRDSs and the utility/success of PBBM, there is also a need to continue the momentum and dialogue between the industry and regulators. Five key areas were identified which require further discussion and harmonisation.

In silico predictions of absorption of MDI substances after dermal or inhalation exposures to support a category based read-across assessment.

Methylenediphenyl diisocyanate (MDI) substances used polyurethane production can range from their simplest monomeric forms (e.g., 4,4'-MDI) to mixtures of the monomers with various homologues, homopolymer, and prepolymer derivatives. The relative dermal or inhalation absorption of 39 constituents of these substances in human were predicted using the GastroPlus® program. Predicted dermal uptake and absorption of the three MDI monomers from an acetone vehicle was 84-86% and 1.4-1.5%, respectively, with lower uptake and absorption predicted for the higher MW analogs. Lower absorption was predicted from exposures in a more lipophilic vehicle (1-octanol). Modeled inhalation exposures afforded the highest pulmonary absorption for the MDI monomers (38-54%), with 3-27% for the MW range of 381-751, and <0.1% for the remaining, higher MW derivatives. Predicted oral absorption, representing mucociliary transport, ranged from 5 to 10% for the MDI monomers, 10-25% for constituents of MW 381-751, and ≤3% for constituents with MW > 900. These in silico evaluations should be useful in category-based, worst-case, Read-Across assessments for MDI monomers and modified MDI substances for potential systemic effects. Predictions of appreciable mucociliary transport may also be useful to address data gaps in oral toxicity testing for this category of compounds.

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.

PBPK Modeling Approach to Predict the Behavior of Drugs Cleared by Kidney in Pregnant Subjects and Fetus.

The purpose of this study was to develop a physiologically based pharmacokinetic (PBPK) model predicting the pharmacokinetics (PK) of different compounds in pregnant subjects. This model considers the differences in tissue sizes, blood flow rates, enzyme expression levels, glomerular filtration rates, plasma protein binding, and other factors affected during pregnancy in both the maternal and fetal models. The PBPKPlus™ module in GastroPlus® was used to model the PK of cefuroxime and cefazolin. For both compounds, the model was first validated against PK data in healthy non-pregnant volunteers and then applied to predict pregnant groups PK. The model accurately described the PK in both non-pregnant and pregnant groups and explained well differences in the plasma concentration due to pregnancy. The fetal plasma and amniotic fluid concentrations were also predicted reasonably well at different stages of pregnancy. This work describes the use of a PBPK approach for drug development and demonstrates the ability to predict differences in PK in pregnant subjects and fetal exposure for compounds excreted renally. The prediction for pregnant groups is also improved when the model is calibrated with postpartum or non-pregnant female group if such data are available.

Physiologically Based Pharmacokinetic Modeling Approach to Identify the Drug-Drug Interaction Mechanism of Nifedipine and a Proton Pump Inhibitor, Omeprazole.

BACKGROUND AND OBJECTIVES: Proton pump inhibitors (PPIs) can affect the intragastric release of other drugs from their dosage forms by elevating the gastric pH. They may also influence drug absorption and metabolism by interacting with P-glycoprotein or with the cytochrome P450 (CYP) enzyme system. Nifedipine is a Biopharmaceutics Classification System (BCS) class II drug with low solubility across physiologic pH and high permeability. Previous studies have demonstrated that drug-drug interaction (DDI) existed between omeprazole and nifedipine with significantly increased systemic exposure of nifedipine in subjects after pre-treatment for 7 days with omeprazole compared to the subjects without omeprazole treatment. It was shown that omeprazole not only induced an increase in intragastric pH, but also inhibited the CYP3A4 activity, while CYP3A4-mediated oxidation is the main metabolic pathway of nifedipine. The purpose of this study is to apply a physiologically based pharmacokinetic (PBPK) modeling approach to investigate the DDI mechanism for an immediate release formulation of nifedipine with omeprazole. METHODS: A previously published model for omeprazole was modified to integrate metabolites and to update CYP inhibition based on the most updated published in vitro data. We simulated the nifedipine pharmacokinetics in healthy subjects with or without the multiple-dose pretreatment of omeprazole (20 mg) following oral administrations of immediate-release (IR) (10 mg) nifedipine. Nifedipine solubility at different pHs was used to simulate the nifedipine pharmacokinetics for both clinical arms. Multiple sensitivity analyses were performed to understand the impact of gastric pH and the CYP3A4-mediated gut and liver first pass metabolism on the overall nifedipine pharmacokinetics. RESULTS: The developed PBPK model properly described the pharmacokinetics of nifedipine and predicted the inhibitory effect of multiple-dose omeprazole on CYP3A4 activity. With the incorporation of the physiologic effect of omeprazole on both gastric pH and CYP3A4 to the PBPK model, the verified PBPK model allows evaluating the impact of the increase in gastric pH and/or CYP3A4 inhibition. The simulated results show that the nifedipine metabolic inhibition by omeprazole may play an important role in the DDI between nifedipine and omeprazole for IR nifedipine formulation. CONCLUSION: The developed full PBPK model with the capability to simulate DDI by considering gastric pH change and metabolic inhibition provides a mechanistic understanding of the observed DDI of nifedipine with a PPI, omeprazole.

Ocular Physiologically Based Pharmacokinetic Modeling for Ointment Formulations.

PURPOSE: The purpose of this study is to show how the Ocular Compartmental Absorption & Transit (OCAT™) model in GastroPlus® can be used to characterize ocular drug pharmacokinetic performance in rabbits for ointment formulations. METHODS: A newly OCAT™ model developed for fluorometholone, as well as a previously verified model for dexamethasone, were used to characterize the aqueous humor (AH) concentration following the administration of multiple ointment formulations to rabbit. The model uses the following parameters: application surface area (SA), a fitted application time, and the fitted Higuchi release constant to characterize the rate of passage of the active pharmaceutical ingredient from the ointment formulations into the tears in vivo. RESULTS: Parameter sensitivity analysis was performed to understand the impact of ointment formulation changes on ocular exposure. While application time was found to have a significant impact on the time of maximal concentration in AH, both the application SA and the Higuchi release constant significantly influenced both the maximum concentration and the ocular exposure. CONCLUSIONS: This initial model for ointment ophthalmic formulations is a first step to better understand the interplay between physiological factors and ophthalmic formulation physicochemical properties and their impact on in vivo ocular drug pharmacokinetic performance in rabbits.

Physiologically Based Pharmacokinetic Model to Support Ophthalmic Suspension Product Development.

FDA's Orange Book lists 17 currently marketed active pharmaceutical ingredients (API) formulated within ophthalmic suspensions in which a majority has 90% or more of the API undissolved. We used an ocular physiologically based pharmacokinetic (O-PBPK) model to compare a suspension with a solution for ophthalmic products with dexamethasone (Dex) as the model drug. Simulations with a Dex suspension O-PBPK model previously verified in rabbit were used to characterize the consequences of drug clearance mechanism in the precorneal compartment on pharmacokinetic (PK) exposure and to assess the ocular and systemic PK characteristics of ophthalmic suspensions with different strengths or magnitudes of viscosity. O-PBPK-based simulations show that (1) Dex suspension 0.05% has a 2.5- and 5-fold higher AUC in aqueous humor and plasma, respectively, than the Dex saturated solution; (2) strength increase by 5- and 10-fold induces a respective 2.2- and 3.3-fold increase in aqueous humor and 4.4- and 8.6-fold increase in plasma Cmax and AUC; and (3) increasing formulation viscosity (from 1.6 to 75 cP) causes an overall increase in API available for absorption in the cornea resulting in a higher ocular Cmax and AUC with no significant impact on systemic exposure. This research demonstrates that solid particles present in a suspension can not only help to achieve a higher ocular exposure but also unfavorably raise systemic exposure. A model able to correlate formulation changes to both ocular and plasma exposure is a necessary tool to support ocular product development taking into consideration both local efficacy and systemic safety aspects.

Application of Mechanistic Ocular Absorption Modeling and Simulation to Understand the Impact of Formulation Properties on Ophthalmic Bioavailability in Rabbits: a Case Study Using Dexamethasone Suspension.

Developing mathematical models to predict changes in ocular bioavailability and pharmacokinetics due to differences in the physicochemical properties of complex topical ophthalmic suspension formulations is important in drug product development and regulatory assessment. Herein, we used published FDA clinical pharmacology review data, in-house, and literature rabbit pharmacokinetic data generated for dexamethasone ophthalmic suspensions to demonstrate how the mechanistic Ocular Compartmental Absorption and Transit model by GastroPlus™ can be used to characterize ocular drug pharmacokinetic performance in rabbits for suspension formulations. This model was used to describe the dose-dependent (0.01 to 0.1%) non-linear pharmacokinetic in ocular tissues and characterize the impact of viscosity (1.67 to 72.9 cP) and particle size (5.5 to 22 µm) on in vivo ocular drug absorption and disposition. Parameter sensitivity analysis (hypothetical suspension particle size: 1 to 10 µm, viscosity: 1 to 100 cP) demonstrated that the interplay between formulation properties and physiological clearance through drainage and tear turnover rates in the pre-corneal compartment drives the ocular drug bioavailability. The quick removal of drug suspended particles from the pre-corneal compartment renders the impact of particle size inconsequential relative to viscosity modification. The in vivo ocular absorption is (1) viscosity non-sensitive when the viscosity is high and the impact of viscosity on the pre-corneal residence time reaches the maximum physiological system capacity or (2) viscosity sensitive when the viscosity is below a certain limit. This study reinforces our understanding of the interplay between physiological factors and ophthalmic formulation physicochemical properties and their impact on in vivo ocular drug PK performance in rabbits.