Integrating In Vitro, Modeling, and In Vivo Approaches to Investigate Warfarin Bioequivalence

A physiologically-based pharmacokinetic/pharmacodynamic modeling approach for drug-drug-gene interaction evaluation of S-warfarin with fluconazole

Warfarin is a widely used anticoagulant, and its S-enantiomer has higher potency compared to the R-enantiomer. S-warfarin is mainly metabolized by cytochrome P450 (CYP) 2C9, and its pharmacological target is vitamin K epoxide reductase complex subunit 1 (VKORC1). Both CYP2C9 and VKORC1 have genetic polymorphisms, leading to large variations in the pharmacokinetics (PKs) and pharmacodynamics (PDs) of warfarin in the population. This makes dosage management of warfarin difficult, especially in the case of drug-drug interactions (DDIs). This study provides a whole-body physiologically-based pharmacokinetic/PD (PBPK/PD) model of S-warfarin for predicting the effects of drug-drug-gene interactions on S-warfarin PKs and PDs. The PBPK/PD model of S-warfarin was developed in PK-Sim and MoBi. Drug-dependent parameters were obtained from the literature or optimized. Of the 34 S-warfarin plasma concentration-time profiles used, 96% predicted plasma concentrations within twofold range compared to observed data. For S-warfarin plasma concentration-time profiles with CYP2C9 genotype, 364 of 386 predicted plasma concentration values (~94%) fell within the twofold of the observed values. This model was tested in DDI predictions with fluconazole as CYP2C9 perpetrators, with all predicted DDI area under the plasma concentration-time curve to the last measurable timepoint (AUClast) ratio within twofold of the observed values. The anticoagulant effect of S-warfarin was described using an indirect response model, with all predicted international normalized ratio (INR) within twofold of the observed values. This model also incorporates a dose-adjustment method that can be used for dose adjustment and predict INR when warfarin is used in combination with CYP2C9 perpetrators.

Predicting bioequivalence and developing dissolution bioequivalence safe space in vitro for warfarin using a Physiologically-Based pharmacokinetic absorption model

OBJECTIVE

Bioequivalence (BE) studies support the approval and clinical use of both new and generic drug products. Narrow therapeutic index (NTI) drugs have relatively high costs and low success rates of BE evaluation clinical trials as high-risk drugs. A physiologically-based pharmacokinetic (PBPK) model can be used to evaluate the BE of two preparations.

METHODS

This study inputs the basic physical and chemical property parameters of warfarin sodium available at the present stage into GastroPlus™ software, and combined it with the Advanced Compartmental Absorption and Transit (ACAT™) model built into the software. The PBPK model of Chinese individuals taking 2.5 mg of warfarin sodium orally while fasted condition was developed using the disposal parameters calculated from the clinically measured PK data of the reference preparations. The model was tested using the PK data of other reference preparations and tested preparations from different domestic manufacturers.

RESULTS

The results revealed that at least 30% of drugs are released in 30 min under a pH of 4.5 condition, and at least 80% are released in 30 min under a pH of 6.8 condition, which can be used as bioequivalent dissolution limits under fasted conditions. The risk of BE failure in the fed condition will be significantly reduced for the clinical study on the BE of warfarin sodium, which is a NTI drug if the fasted condition is bioequivalent.

CONCLUSION

The results revealed that the PBPK models were successfully developed for 2.5 mg of warfarin sodium tablets in Chinese individuals. Developing a PBPK model for NTI drugs based on in vitro dissolution data in software is a promising method for BE evaluation, which can provide great help for developing new drugs and the clinical trial research of BE of generic drugs.

Physiologically Based Pharmacokinetic Model Development and Verification for Bioequivalence Testing of Bempedoic Acid Oral Suspension and Reference Tablet Formulation

The bioequivalence of bempedoic acid oral suspension and commercial immediate release (IR) tablet formulations were assessed using a physiologically based pharmacokinetic (PBPK) model. The mechanistic model, developed from clinical mass balance results and in vitro intrinsic solubility, permeability, and dissolution data, was verified against observed clinical pharmacokinetics (PK) results. Model inputs included a fraction of a dose in solution (0.01%), viscosity (118.8 cps), and median particle diameter (50 µm) for the suspension and particle diameter (36.4 µm) for IR tablets. Dissolution was determined in the relevant media (pH 1.2-6.8) in vitro. Model simulations of bioequivalence predicted oral suspension (test) to IR tablet (reference) geometric mean ratio estimates of 96.9% (90% confidence interval [CI]: 92.6-101) for maximum concentration and 98.2% (90% CI: 87.3-111) for the area under the concentration-time curve. Sensitivity analyses showed gastric transit time had a minor impact on model predictions. Oral suspension biopharmaceutical safe space was defined by extremes of particle size and the percent of bempedoic acid in solution. PBPK model simulations predicted that the rate and extent of bempedoic acid absorption are unlikely to exhibit clinically meaningful differences when dosed as an oral suspension compared with an IR tablet without requiring a clinical bioequivalence study in adults.

Conjunction of semi-mechanistic in vitro-in vivo modeling and population pharmacokinetics as a tool for virtual bioequivalence analysis - a case study for a BCS class II drug

Virtual bioequivalence trial (VBE) simulations based on (semi)mechanistic in vitro-in vivo (IVIV) modeling have gained a huge interest in the pharmaceutical industry. Sophisticated commercially available software allows modeling variable drug fates in the gastrointestinal tract (GIT). Surprisingly, the between-subject and inter-occasion variability (IOV) of the distribution volumes and clearances are ignored or simplified, despite substantially contributing to varied plasma drug concentrations. The paper describes a novel approach for IVIV-based VBE by using population pharmacokinetics (popPK). The data from two bioequivalence trials with a poorly soluble BCS class II drug were analyzed retrospectively. In the first trial, the test drug product (biobatch 1) did not meet the bioequivalence criteria, but after a reformulation, the second trial succeeded (biobatch 2). The popPK model was developed in the Monolix software (Lixoft SAS, Simulation Plus) based on the originator's plasma concentrations. The modified Noyes-Whitney model was fitted to the results of discriminative biorelevant dissolution tests of the two biobatches and seven other reformulations. Then, the IVIV model was constructed by joining the popPK model with fixed drug disposition parameters, the drug dissolution model, and mechanistic approximation of the GIT transit. It was used to simulate the drug concentrations at different IOV levels of the primary pharmacokinetic parameters and perform the VBE. Estimated VBE success rates for both biobatches well reflected the outcomes of the bioequivalence trials. The predicted 90% confidence intervals for the area under the time-concentration curves were comparable with the observed values, and the 10% IOV allowed the closest approximation to the clinical results. Simulations confirmed that a significantly lower maximum drug concentration for biobatch 1 was responsible for the first clinical trial's failure. In conclusion, the proposed workflow might aid formulation screening in generic drug development.

Physiologically Based Absorption Modelling to Explore the Formulation and Gastric pH Changes on the Pharmacokinetics of Acalabrutinib

Acalabrutinib, a selective Bruton's tyrosine kinase inhibitor, is a biopharmaceutics classification system class II drug. The aim of this study was to develop a physiologically based pharmacokinetic (PBPK) model to mechanistically describe absorption of immediate release capsule formulation of acalabrutinib in humans. Integration of in vitro biorelevant measurements, dissolution studies and in silico modelling provided clinically relevant inputs for the mechanistic absorption PBPK model. The batch specific dissolution data were integrated in two ways, by fitting a diffusion layer model scalar to the drug product dissolution with integration of drug substance laser diffraction particle size data, or by fitting a product particle size distribution to the dissolution data. The latter method proved more robust and biopredictive. In both cases, the drug surface solubility was well predicted by the Simcyp simulator. The model using the product particle size distribution (P-PSD) for each clinical batch adequately captured the PK profiles of acalabrutinib and its active metabolite. Average fold errors were 0.89 for both Cmax and AUC, suggesting good agreement between predicted and observed PK values. The model also accurately predicted pH-dependent drug-drug interactions between omeprazole and acalabrutinib, which was similar across all clinical formulations. The model predicted acalabrutinib geometric mean AUC ratios (with omeprazole vs acalabrutinib alone) were 0.51 and 0.68 for 2 batches of formulations, which are close to observed values of 0.43 and 0.51~0.63, respectively. The mechanistic absorption PBPK model could be potentially used for future applications such as optimizing formulations or predicting the PK for different batches of the drug product.

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.

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.

Structural Changes of Sodium Warfarin in Tablets Affecting the Dissolution Profiles and Potential Safety of Generic Substitution

At present, the risk of generic substitutions in warfarin tablets is still being discussed. The aim of this study was to assess whether API interactions with commonly used excipients may affect the safety of generic replacement of warfarin sodium tablets. These interactions were observed during an accelerated stability study, and the effect of the warfarin solid phase (crystalline/amorphous form) as well as the API particle size distribution was studied. Commercial tablets and prepared tablets containing crystalline warfarin or amorphous warfarin were used. In addition, binary mixtures of warfarin with various excipients were prepared. The structural changes before and after the stability study were monitored by dissolution test in different media, solid-state NMR spectroscopy and Raman microscopy. During the stability study, the conversion of the sodium in warfarin to its acid form was demonstrated by some excipients (e.g., calcium phosphate). This change in the solid phase of warfarin leads to significant changes in dissolution, especially with the different particle sizes of the APIs in the tablet. Thus, the choice of suitable excipients and particle sizes are critical factors influencing the safety of generic warfarin sodium tablets.

In Silico Modeling and Simulation to Guide Bioequivalence Testing for Oral Drugs in a Virtual Population

Model-informed drug discovery and development (MID3) shows great advantages in facilitating drug development. A physiologically based pharmacokinetic model is one of the powerful computational approaches of MID3, and the emerging field of virtual bioequivalence is well recognized to be the future of the physiologically based pharmacokinetic model. Based on the translational link between in vitro, in silico, and in vivo, virtual bioequivalence study can evaluate the similarity and potential difference of pharmacokinetic and clinical performance between test and reference formulations. With the aid of virtual bioequivalence study, the pivotal information of clinical trials can be provided to streamline the development for both new and generic drugs. However, a regulatory framework of virtual bioequivalence study has not reached its full maturity. Therefore, this article aims to present an overview of the current status of bioequivalence study, identify the framework of virtual bioequivalence studies for oral drugs, and also discuss the future opportunities of virtual bioequivalence in supporting the waiver and optimization of in vivo clinical trials.

Molecular Insights into Warfarin Sodium 2-Propanol Solvate Solid Form Changes and Disproportionation Using a Low Volume Two-Stage Dissolution Approach

The current research work focuses on understanding the reported discrepancies and our observations in the dissolution profiles of warfarin sodium tablets and potential patient-based failure modes during oral warfarin therapy. It was hypothesized that freely soluble crystalline warfarin sodium (WAR_C) at first transforms into noncrystalline warfarin sodium (WAR_NC) under stress conditions. The WAR_C → WAR_NC conversion facilitates the rapid formation of the poorly soluble unionized form, which could lead to dissolution failures and potential poor in vivo performance. Depressed warfarin concentrations locally in the gastrointestinal tract (GIT) may in turn lead to inadequate absorption and thereby affect bioavailability. A low volume two-stage dissolution method was developed to mimic in vivo GIT conditions. Warfarin sodium tablets exposed to room temperature and 75% relative humidity for 1 week showed approximately 23% decrease in drug release. The decline in drug release supports the hypothesis that WAR_NC is converted to the unionized form faster than WAR_C does under the same conditions. Solid state characterization (powder X-ray diffractometry and differential scanning calorimetry) data demonstrated the disproportionation of warfarin sodium to unionized warfarin after solubility and dissolution studies. The findings support the hypothesis and a possible failure mode of warfarin sodium tablets. This work is a second case study from our laboratory on narrow therapeutic index drug products in which the instability of the solid state of the drug substance is potentially responsible for observed clinical failures.

Physiologically Based Pharmacokinetic/Pharmacodynamic Modeling to Support Waivers of In Vivo Clinical Studies: Current Status, Challenges, and Opportunities

Physiologically based pharmacokinetic/pharmacodynamic (PBPK/PD) modeling has been extensively applied to quantitatively translate in vitro data, predict the in vivo performance, and ultimately support waivers of in vivo clinical studies. In the area of biopharmaceutics and within the context of model-informed drug discovery and development (MID3), there is a rapidly growing interest in applying verified and validated mechanistic PBPK models to waive in vivo clinical studies. However, the regulatory acceptance of PBPK analyses for biopharmaceutics and oral drug absorption applications, which is also referred to variously as "PBPK absorption modeling" [Zhang et al. CPT: Pharmacometrics Syst. Pharmacol. 2017, 6, 492], "physiologically based absorption modeling", or "physiologically based biopharmaceutics modeling" (PBBM), remains rather low [Kesisoglou et al. J. Pharm. Sci. 2016, 105, 2723] [Heimbach et al. AAPS J. 2019, 21, 29]. Despite considerable progress in the understanding of gastrointestinal (GI) physiology, in vitro biopharmaceutic and in silico tools, PBPK models for oral absorption often suffer from an incomplete understanding of the physiology, overparameterization, and insufficient model validation and/or platform verification, all of which can represent limitations to their translatability and predictive performance. The complex interactions of drug substances and (bioenabling) formulations with the highly dynamic and heterogeneous environment of the GI tract in different age, ethnic, and genetic groups as well as disease states have not been yet fully elucidated, and they deserve further research. Along with advancements in the understanding of GI physiology and refinement of current or development of fully mechanistic in silico tools, we strongly believe that harmonization, interdisciplinary interaction, and enhancement of the translational link between in vitro, in silico, and in vivo will determine the future of PBBM. This Perspective provides an overview of the current status of PBBM, reflects on challenges and knowledge gaps, and discusses future opportunities around PBPK/PD models for oral absorption of small and large molecules to waive in vivo clinical studies.

A physiologically based pharmacokinetic/pharmacodynamic modeling approach for drug-drug interaction evaluation of warfarin enantiomers with sorafenib

Sorafenib was suggested to cause drug-drug interaction (DDI) with the common anticoagulant, warfarin based on published studies. The inhibition on CYP2C9 enzyme was thought to be the mechanism, but further studies are warranted. Thus, a mechanistic PBPK/PD model for warfarin enantiomers was developed to predict DDI potential with sorafenib, aiming at providing reference for the rational use of both drugs. PBPK models of warfarin enantiomers were constructed by Simcyp software. A mechanistic PK/PD model was built in NONMEM software. PBPK model of sorafenib was fitted via a top-down method. The final PBPK/PD model of warfarin enantiomers was verified and validated by different dosing regimens, ethnicities and genetic polymorphisms, and used to perform DDI simulations between warfarin racemate and sorafenib among general populations and sub-populations with various CYP2C9 and VKORC1 genotypes. Results suggested low DDI risk between warfarin and sorafenib for general populations. Potentially serious consequence was seen for those carrying both CYP2C9 ∗2 and ∗3 and VKORC1 A/A genotypes. This PBPK/PD modeling approach for warfarin enantiomers enabled DDI evaluation with sorafenib. Close monitoring and warfarin dosage adjustment were recommended for patients carrying mutant genotypes. The novel model could be applied to investigate other drugs that may interact with warfarin.

Evaluation of the Impacts of Formulation Parameters on the Pharmacokinetics and Bioequivalence of Risperidone Orodispersible Film: a Physiologically Based Pharmacokinetic Modeling Approach

The purpose of this study was to investigate the impacts of the formulation parameters on the pharmacokinetics and bioequivalence of risperidone orodispersible film (ODF) using physiologically based pharmacokinetic model. The pharmacokinetic profiles of two risperidone ODFs, which exhibit different in vitro dissolution, were examined in Beagle dogs after supralingual administration. Subsequently, a physiologically based pharmacokinetic (PBPK) model was constructed to evaluate the in vivo performance of risperidone ODF. The parameter sensitivity analysis (PSA) was used to access the impacts of formulation parameters on the pharmacokinetics of risperidone. Moreover, the validated PBPK model was applied to predict human pharmacokinetic profiles and examine the bioequivalence of these two ODFs. These two ODFs displayed similar risperidone pharmacokinetic profiles in dogs. The parameter sensitivity analysis indicated that the changes in the solubility, particle size, particle density, and diffusion coefficient did not have obvious influence on the in vivo properties of risperidone ODF. Alternation of the in vitro complete dissolution time in water from 15 to 30 min led to a 30% decrease in C_max and 20% of increase in T_max. AUC_0-∞ would be decreased if risperidone was not fully released within 1 h. As both ODFs completely released risperidone within 15 min, the difference in the extent of in vivo absorption, intestinal regional absorption location, and plasma concentration-time curves between these two ODFs was almost negligible. Consequently, a bioequivalence was foreseen in humans. The in vitro cumulative dissolution percentage in water at 15 min was found to be the major determinant on the in vivo properties of risperidone ODF. PBPK modeling appears to be an innovative strategy to guide the development of risperidone ODF.

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.

Development and validation of an ultra-high-performance liquid chromatography-tandem mass spectrometry method to determine the bioavailability of warfarin and its major metabolite 7-hydroxy warfarin in rats dosed with oral formulations containing different polymorphic forms

A simple, sensitive and rapid ultra-high-performance liquid chromatography-electrospray ionization-tandem mass spectrometry method was developed and validated for the quantification of warfarin and 7-hydroxy warfarin in Sprague Dawley (SD) rats. Animals were administered a single dose of warfarin sodium formulations (crystalline and amorphous) at 12 mg/kg via oral gavage and blood was drawn over a 96-h time course. Sample process recoveries, matrix effect and analyte stability were determined. The linearity for warfarin and 7-hydroxy warfarin was from 5 to 2000 ng/mL in blank SD rat plasma. Correlation coefficients (r² ) for standard calibration curves were >.98 and analytes quantified within ±15% of target at all calibrator concentrations. The average percent accuracy and precision for intra- and inter-day were 93.7%-113.8% and ≤12.1%, respectively, for warfarin and 7-hydroxy warfarin, across the quality control standards (5, 10, 500, 1800 and 2000 ng/mL). Acceptable analytical recovery (>55%) was achieved with process efficiencies >41.5% and matrix effects <139.9% over the analytical range. Both analytes were stable in stock solution, autosampler, benchtop and three cycles of freeze-thaw with percent accuracy ≥90.2% and precision (percent relative standard deviation) ≤14%. The validated method was successfully applied to a pre-clinical bioavailability study of crystalline and amorphous warfarin sodium formulations in SD rats.