The Use of Physiologically Based Pharmacokinetic Analyses-in Biopharmaceutics Applications -Regulatory and Industry Perspectives

More Than a Gut Feeling─A Combination of Physiologically Driven Dissolution and Pharmacokinetic Modeling as a Tool for Understanding Human Gastric Motility

In vivo studies of formulation performance with in vitro and/or in silico simulations are often limited by significant gaps in our knowledge of the interaction between administered dosage forms and the human gastrointestinal tract. This work presents a novel approach for the investigation of gastric motility influence on dosage form performance, by combining biopredictive dissolution tests in an innovative PhysioCell apparatus with mechanistic physiology-based pharmacokinetic modeling. The methodology was based on the pharmacokinetic data from a large (n = 118) cohort of healthy volunteers who ingested a capsule containing a highly soluble and rapidly absorbed drug under fasted conditions. The developed dissolution tests included biorelevant media, varied fluid flows, and mechanical stress events of physiological timing and intensity. The dissolution results were used as inputs for pharmacokinetic modeling that led to the deduction of five patterns of gastric motility and their prevalence in the studied population. As these patterns significantly influenced the observed pharmacokinetic profiles, the proposed methodology is potentially useful to other in vitro-in vivo predictions involving immediate-release oral dosage forms.

Parameterization of Physiologically Based Biopharmaceutics Models: Workshop Summary Report

This Article shares the proceedings from the August 29th, 2023 (day 1) workshop "Physiologically Based Biopharmaceutics Modeling (PBBM) Best Practices for Drug Product Quality: Regulatory and Industry Perspectives". The focus of the day was on model parametrization; regulatory authorities from Canada, the USA, Sweden, Belgium, and Norway presented their views on PBBM case studies submitted by industry members of the IQ consortium. The presentations shared key questions raised by regulators during the mock exercise, regarding the PBBM input parameters and their justification. These presentations also shed light on the regulatory assessment processes, content, and format requirements for future PBBM regulatory submissions. In addition, the day 1 breakout presentations and discussions gave the opportunity to share best practices around key questions faced by scientists when parametrizing PBBMs. Key questions included measurement and integration of drug substance solubility for crystalline vs amorphous drugs; impact of excipients on apparent drug solubility/supersaturation; modeling of acid-base reactions at the surface of the dissolving drug; choice of dissolution methods according to the formulation and drug properties with a view to predict the in vivo performance; mechanistic modeling of in vitro product dissolution data to predict in vivo dissolution for various patient populations/species; best practices for characterization of drug precipitation from simple or complex formulations and integration of the data in PBBM; incorporation of drug permeability into PBBM for various routes of uptake and prediction of permeability along the GI tract.

Conventional vs Mechanistic IVIVC: A Comparative Study in Establishing Dissolution Safe Space for Extended Release Formulations

The use of in vitro-in vivo correlation (IVIVC) for extended release oral dosage forms is an important technique that can avoid potential clinical studies. IVIVC has been a topic of discussion over the past two decades since the inception of USFDA guidance. It has been routinely used for biowaivers, establishment of dissolution safe space and clinically relevant dissolution specifications, for supporting site transfers, scale-up and post approval changes. Although conventional or mathematical IVIVC is routinely used, other approach such as mechanistic IVIVC can be of attractive choice as it integrates all the physiological aspects. In the present study, we have performed comparative evaluation of mechanistic and conventional IVIVC for establishment of dissolution safe space using divalproex sodium and tofacitinib extended release formulations as case examples. Conventional IVIVC was established using Phoenix and mechanistic IVIVC was set up using Gastroplus physiologically based biopharmaceutics model (PBBM). Virtual dissolution profiles with varying release rates were constructed around target dissolution profile using Weibull function. After internal and external validation, the virtual dissolution profiles were integrated into mechanistic and conventional IVIVC and safe space was established by absolute error and T/R ratio's methods. The results suggest that mechanistic IVIVC yielded wider safe space as compared to conventional IVIVC. The results suggest that a mechanistic approach of establishing IVIVC may be a flexible approach as it integrates physiological aspects. These findings suggest that mechanistic IVIVC has wider potential as compared to conventional IVIVC to gain wider dissolution safe space and thus can avoid potential clinical studies.

Physiologically Based Biopharmaceutics Modeling (PBBM): Best Practices for Drug Product Quality, Regulatory and Industry Perspectives: 2023 Workshop Summary Report

Physiologically based biopharmaceutics modeling (PBBM) is used to elevate drug product quality by providing a more accurate and holistic understanding of how drugs interact with the human body. These models are based on the integration of physiological, pharmacological, and pharmaceutical data to simulate and predict drug behavior in vivo. Effective utilization of PBBM requires a consistent approach to model development, verification, validation, and application. Currently, only one country has a draft guidance document for PBBM, whereas other major regulatory authorities have had limited experience with the review of PBBM. To address this gap, industry submitted confidential PBBM case studies to be reviewed by the regulatory agencies; software companies committed to training. PBBM cases were independently and collaboratively discussed by regulators, and academic colleagues participated in some of the discussions. Successful bioequivalence "safe space" industry case examples are also presented. Overall, six regulatory agencies were involved in the case study exercises, including ANVISA, FDA, Health Canada, MHRA, PMDA, and EMA (experts from Belgium, Germany, Norway, Portugal, Spain, and Sweden), and we believe this is the first time such a collaboration has taken place. The outcomes were presented at this workshop, together with a participant survey on the utility and experience with PBBM submissions, to discuss the best scientific practices for developing, validating, and applying PBBMs. The PBBM case studies enabled industry to receive constructive feedback from global regulators and highlighted clear direction for future PBBM submissions for regulatory consideration.

Research Landscape of Physiologically Based Pharmacokinetic Model Utilization in Different Fields: A Bibliometric Analysis (1999-2023)

PURPOSE

The physiologically based pharmacokinetic (PBPK) modeling has received increasing attention owing to its excellent predictive abilities. However, there has been no bibliometric analysis about PBPK modeling. This research aimed to summarize the research development and hot points in PBPK model utilization overall through bibliometric analysis.

METHODS

We searched for publications related to the PBPK modeling from 1999 to 2023 in the Web of Science Core Collection (WoSCC) database. The Microsoft Office Excel, CiteSpace and VOSviewers were used to perform the analyses.

RESULTS

A total of 4,649 records from 1999 to 2023 were identified, and the largest number of publications focused in the period 2018-2023. The United States was the leading country, and the Environmental Protection Agency (EPA) was the leading institution. The journal Drug Metabolism and Disposition published and co-cited the most articles. Drug-drug interactions, special populations, and new drug development are the main topics in this research field.

CONCLUSION

We first visualize the research landscape and hotspots of the PBPK modeling through bibliometric methods. Our study provides a better understanding for researchers, especially beginners about the dynamization of PBPK modeling and presents the relevant trend in the future.

Predicting Surface pH in Unbuffered Conditions for Acids, Bases, and Their Salts - A Review of Modeling Approaches and Their Performance

Dissolution of ionizable drugs and their salts is a function of drug surface solubility driven by the surface pH, i.e., the microenvironmental pH at the solid/liquid interface, which will deviate from bulk pH when there is an acid-base reaction occurring at the solid/liquid interface. In this work, we first present a brief overview of the modeling approaches available in the literature, classified according to the rate-determining step assumed in the dissolution process. In the second part, we present and evaluate the prediction performance of two different modeling approaches for surface pH. The first method relies only on thermodynamic equilibria, while the second method accounts for transport phenomena of charged compounds through the diffusional boundary layer using the Nernst - Planck equation. Model outcomes are compared with experimental data taken from the literature and obtained during this work. In terms of surface pH predictions, the models provide identical values for weak acids or weak bases. The models' outcomes for bases are in good agreement with experimental data in acidic conditions (bulk pH 1-4), while overpredictions are observed in the 5-7 bulk pH range in a system-dependent manner. Deviations can be related to the effect of surface dissolution (also referred to as surface reaction), which may become a controlling mechanism and slow the replenishment of the unionized drug at the surface of the crystal. Surface pH predictions for acids are generally in good agreement with experiments, with a slight underestimation for some drug examples, which could be related to errors in intrinsic solubility determination or to the assumption of thermodynamic equilibrium at the surface of the drug. A good agreement is also observed for salts with the thermodynamic model except for mesylate salts, suggesting that other phenomena, not currently included in the thermodynamic equilibrium model, may determine the surface pH.

Role of Physiologically Based Biopharmaceutics Modeling (PBBM) in Fed Bioequivalence Study Waivers: Regulatory Outlook, Case Studies and Future Perspectives

Over the past few decades, physiologically based biopharmaceutics modeling (PBBM) has demonstrated its utility in both new drug and generic product development. Applications of PBBM for fed bioequivalence study waivers is an upcoming area. Recently Innovation & Quality (IQ) consortium demonstrated utility of PBBM to avoid repeat food effect studies for new drugs. In the similar lines, the current manuscript aims to discuss role of PBBM in generic fed bioequivalence study waivers. Generic industry practices related to PBBM model development to predict fed bioequivalence was portrayed with special emphasis on fed bio-predictive media. Media that can simulate fed bioequivalence study outcome were discussed from practical perspective. In-depth analysis, collating the data from 36 products was performed to understand predictability of PBBM for fed bioequivalence. Cases where PBBM was successful to predict fed bioequivalence was correlated with BCS class, formulation category and type of food effect. Further, two case studies were presented wherein fed bioequivalence study waiver obtained with PBBM approach. Lastly, future direction in terms of fed bioequivalence study waivers, regulatory perspectives and best practices for PBBM were portrayed. Overall, this article paves a way to utilize PBBM for generic fed bioequivalence study waivers.

Physiologically Based Pharmacokinetic Modeling of the Drug-Drug Interaction Between CYP3A4 Substrate Glasdegib and Moderate CYP3A4 Inducers in Lieu of a Clinical Study

Glasdegib (DAURISMO) is a hedgehog pathway inhibitor approved for the treatment of acute myeloid leukemia (AML). Cytochrome P450 3A4 (CYP3A4) has been identified as a major metabolism and clearance pathway for glasdegib. The role of CYP3A4 in the clearance of glasdegib has been confirmed with clinical drug-drug interaction (DDI) studies following the coadministration of glasdegib with the strong CYP3A4 inhibitor ketoconazole and the strong inducer rifampin. To evaluate potential drug interactions with CYP3A4 modulators, the coadministration of glasdegib with a moderate CYP3A4 inducer, efavirenz, was evaluated using physiologically based pharmacokinetic (PBPK) modeling using the Simcyp simulator. The glasdegib compound file was developed using measured physicochemical properties, data from human intravenous and oral pharmacokinetics, absorption, distribution, metabolism, and excretion studies, and in vitro reaction phenotyping results. The modeling assumptions, model parameters, and assignments of fractional CYP3A4 metabolism were verified using results from clinical pharmacokinetics (PK) and DDI studies with ketoconazole and rifampin. The verified glasdegib and efavirenz compound files, the latter of which was available in the Simcyp simulator, were used to estimate the potential impact of efavirenz on the PK of glasdegib. PBPK modeling predicted a glasdegib area under the concentration-time curve ratio of 0.45 and maximum plasma concentration ratio of 0.75 following coadministration with efavirenz. The PBPK results, in lieu of a formal clinical study, informed the drug label, with the recommendation to double the clinical dose of glasdegib when administered in conjunction with a moderate CYP3A4 inducer, followed by a resumption of the original dose 7 days post-discontinuation.

Approaches to Account for Colon Absorption in Physiologically Based Biopharmaceutics Modeling of Extended-Release Drug Products

The rate and extent of colon absorption are important determinants of the in vivo performance of extended-release (ER) drug products. The ability to appropriately predict this at different stages of development using mechanistic physiologically based biopharmaceutic modeling (PBBM) is highly desirable. This investigation aimed to evaluate the prediction performance of three different approaches to account for colon absorption in predictions of the in vivo performance of ER drug product variants with different in vitro release profiles. This was done by mechanistic predictions of the absorption and plasma exposure of the ER drug products using GastroPlus and GI-Sim for five drugs with different degrees of colon absorption limitations in humans. Colon absorption was accounted for in the predictions using three different approaches: (1) by an a priori approach using the default colon models, (2) by fitting the colon absorption scaling factors to the observed plasma concentration-time profiles after direct administration to the colon in humans, or (3) from the ER drug product variant with the slowest in vitro release profile. The prediction performance was evaluated based on the percentage prediction error and the average absolute prediction error (AAPE). Two levels of acceptance criteria corresponding to highly accurate (AAPE ≤ 20%) and accurate (AAPE 20-50%) predictions were defined prior to the evaluation. For the a priori approach, the relative bioavailability (Frel), AUC0-t, and Cmax of the ER drug product variants for the low to medium colon absorption limitation risk drugs was accurately predicted with an AAPE range of 11-53 and 8-59% for GastroPlus and GI-Sim, respectively. However, the prediction performance was poor for the high colon absorption limitation risk drugs. Moreover, accounting for the human regional colon absorption data in the models did not improve the prediction performance. In contrast, using the colon absorption scaling factors derived from the slowest ER variant significantly improved the prediction performance regardless of colon absorption limitation, with a majority of the predictions meeting the high accuracy criteria. For the slowest ER approach, the AAPE ranges were 5-24 and 5-32% for GastroPlus and GI-Sim, respectively, excluding the low permeability drug. In conclusion, the a priori PBBM can be used during candidate selection and early product design to predict the in vivo performance of ER drug products for low to medium colon absorption limitation risk drugs with sufficient accuracy. The results also indicate a limited value in performing human regional absorption studies in which the drug is administered to the colon as a bolus to support PBBM development for ER drug products. Instead, by performing an early streamlined relative bioavailability study with the slowest relevant ER in vitro release profile, a highly accurate PBBM suitable for ER predictions for commercial and regulatory applications can be developed, except for permeability-limited drugs.

Adapting physiologically-based pharmacokinetic models for machine learning applications

Both machine learning and physiologically-based pharmacokinetic models are becoming essential components of the drug development process. Integrating the predictive capabilities of physiologically-based pharmacokinetic (PBPK) models within machine learning (ML) pipelines could offer significant benefits in improving the accuracy and scope of drug screening and evaluation procedures. Here, we describe the development and testing of a self-contained machine learning module capable of faithfully recapitulating summary pharmacokinetic (PK) parameters produced by a full PBPK model, given a set of input drug-specific and regimen-specific information. Because of its widespread use in characterizing the disposition of orally administered drugs, the PBPK model chosen to demonstrate the methodology was an open-source implementation of a state-of-the-art compartmental and transit model called OpenCAT. The model was tested for drug formulations spanning a large range of solubility and absorption characteristics, and was evaluated for concordance against predictions of OpenCAT and relevant experimental data. In general, the values predicted by the ML models were within 20% of those of the PBPK model across the range of drug and formulation properties. However, summary PK parameter predictions from both the ML model and full PBPK model were occasionally poor with respect to those derived from experiments, suggesting deficiencies in the underlying PBPK model.

Biopharmaceutics Risk Assessment-Connecting Critical Bioavailability Attributes with In Vitro, In Vivo Properties and Physiologically Based Biopharmaceutics Modeling to Enable Generic Regulatory Submissions

Quality risk assessment following ICH Q9 principles is an important activity to ensure optimal clinical efficacy and safety of a drug product. Typically, risk assessment is focused on product performance wherein critical material attributes, formulation variables, and process parameters are evaluated from a manufacturing perspective. Extending ICH Q9 principles to biopharmaceutics risk assessment to identify factors that can impact in vivo performance is an upcoming area. This is evident by recent regulatory trends wherein a new term critical bioavailability attributes (CBA) has been coined to identify such factors. Although significant work has been performed for biopharmaceutics risk assessment for new molecules, there is a need for harmonized biopharmaceutics risk assessment workflow for generic submissions. In this manuscript, we attempted to provide a framework for performing biopharmaceutics risk assessment for generic regulatory submissions. A detailed workflow for performing biopharmaceutics risk assessment includes identification of initial CBA (iCBA), their confirmatory evaluation followed by definition of the control strategy. Tools for biopharmaceutics risk assessment, i.e., bio-discriminatory dissolution method and physiologically based biopharmaceutics modeling (PBBM) were discussed from a practical perspective. Furthermore, a case study for CBA evaluation using PBBM modeling for an extended-release product for regulatory submission has been described using the proposed workflow. Finally, future directions of integrating CBA evaluation, biopharmaceutics risk assessment to the FDA Knowledge Aided Structured Assessment (KASA) initiative, the necessity of risk assessment templates, and knowledge sharing between industry and academia are discussed. Overall, the work described in this manuscript can facilitate and provide guidance for biopharmaceutics risk assessment for generic submissions.

Modelling Based Approaches to Support Generic Drug Regulatory Submissions-Practical Considerations and Case Studies

Model informed drug development (MiDD) is useful to predict in vivo exposure of drugs during various stages of the drug development process. This approach employs a variety of quantitative tools to assess the risks during the drug development process. One important tool in the MiDD tool kit is the Physiologically Based Pharmacokinetic Modelling (PBPK). This tool is extensively used to reduce the development cost and to accelerate the access of medicines to the patients. In this work, we provide an overview of PBPK modelling approaches in the generic drug development process, with a special emphasis on the bio-waiver applications. We describe herein approaches and common pitfalls while submitting model based justifications as a response to the regulatory deficiencies during the generic drug development process. With some in-house case studies, we have attempted to provide a clear path for PBPK model based justifications for bio-waivers. With this review, the gap between theoretical knowledge and practical application of modelling and simulation tools for generic drug product development could be potentially reduced.

Application of a novel PhysioCell apparatus for biopredictive dissolution tests of oral immediate release formulations - a case study workflow for in vitro-in vivo predictions

Biorelevant dissolution tests of oral solid dosage forms open the gate to valid in vitro-in vivo predictions (IVIVP). A recently developed apparatus, PhysioCell, allows mimicking the fluid flow and pressure waves occurring in the human fasted stomach. In this work, we used the PhysioCell to perform IVIVP for vortioxetine immediate-release (IR) tablets: the originator (Brintellix) and generic product candidates (VORTIO). The dissolved drug was monitored in the gastric (StressCell) and intestinal (Collection Vessel) compartments that contained biorelevant media. Simulated intermittent gastric stress at 15 min and 'housekeeping wave' at 30 min increased the dissolution of Brintellix formulations only. A mechanistic model that best described the observations involved the first-order tablet disintegration with a stress-induced enhancement for Brintellix, dissolution of solid particles in the StressCell, and drug transfer to the Collection Vessel. Then, a semi-mechanistic pharmacokinetic model with dissolution parameters as inputs simulated vortioxetine plasma concentrations in healthy volunteers after single and multiple dosing of Brintellix. Despite different dissolution characteristics, VORTIO provided similar concentration profiles to the originator. In conclusion, PhysioCell dissolution tests, combined with semi-mechanistic IVIVP, can be successfully used to develop IR dosage forms exhibiting gastric stress-related effects.

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.

Acalabrutinib Maleate Tablets: The Physiologically Based Biopharmaceutics Model behind the Drug Product Dissolution Specification

Acalabrutinib maleate tablets correspond to an improved formulation compared to acalabrutinib capsules as they can be dosed with and without acid reducing agents and therefore benefit more cancer patients. The dissolution specification for the drug product was determined using all the information available on the drug safety, efficacy, and in vitro performance. In addition, a physiologically based biopharmaceutics model was developed for acalabrutinib maleate tablets on the back of a previously published model for acalabrutinib capsules to establish that the proposed drug product dissolution specification would ensure safe and effective products for all patients including those under acid reducing agent treatment. The model was built, validated, and used to predict the exposure of virtual batches where the dissolution was slower than that of the clinical target. A combination of exposure prediction and the use of a PK-PD model allowed it to be demonstrated that the proposed drug product dissolution specification was acceptable. This combination of models enabled a larger safe space than would have been granted by consideration of bioequivalence only.

Best Practices for Integration of Dissolution Data into Physiologically Based Biopharmaceutics Models (PBBM): A Biopharmaceutics Modeling Scientist Perspective

Dissolution is considered as a critical input into physiologically based biopharmaceutics models (PBBM) as it governs in vivo exposure. Despite many workshops, initiatives by academia, industry, and regulatory, wider practices are followed for dissolution data input into PBBM models. Due to variety of options available for dissolution data input into PBBM models, it is important to understand pros, cons, and best practices while using specific dissolution model. This present article attempts to summarize current understanding of various dissolution models and data inputs in PBBM software's and aims to discuss practical challenges and ways to overcome such scenarios. Different approaches to incorporate dissolution data for immediate, modified, and delayed release formulations are discussed in detail. Common challenges faced during fitting of z-factor are discussed along with novel approach of dissolution data incorporation using P-PSD model. Ways to incorporate dissolution data for MR formulations using Weibull and IVIVR approaches were portrayed with examples. Strategies to incorporate dissolution data for DR formulations was depicted along with practical aspects. Approaches to generate virtual dissolution profiles, using Weibull function, DDDPlus, and time scaling for defining dissolution safe space, and strategies to generate virtual dissolution profiles for justifying single and multiple dissolution specifications were discussed. Finally, novel ways to integrate dissolution data for complex products such as liposomes, data from complex dissolution systems, importance of precipitation, and bio-predictive ability of QC media for evaluation of CBA's impact were discussed. Overall, this article aims to provide an easy guide for biopharmaceutics modeling scientist to integrate dissolution data effectively into PBBM models.

A mechanistic physiologically based model to assess the effect of study design and modified physiology on formulation safe space for virtual bioequivalence of dermatological drug products

Physiologically based pharmacokinetic (PBPK) models are widely accepted tools utilised to describe and predict drug pharmacokinetics (PK). This includes the use of dermal PBPK models at the regulatory level including virtual bioequivalence (VBE) studies. The current work considers the Topicort® Spray formulation, which contains 0.25% desoximetasone (DSM), as an example formulation. Quantitative formulation composition and in vitro permeation testing (IVPT) data were obtained from the public literature to develop a mechanistic model using the multi-phase, multi-layer (MPML) MechDermA IVPT module in the Simcyp Simulator. In vitro-in vivo extrapolation functionality was used to simulate in vivo PK for various scenarios and predict a 'safe space' for formulation bioequivalence using the VBE module. The potential effect of vasoconstriction, impaired barrier function, and various dosing scenarios on the formulation safe space was also assessed. The model predicted 'safe space' for formulation solubility suggesting that a 50% change in solubility may cause bio-in-equivalence, whereas viscosity could deviate by orders of magnitude and the formulation may still remain bioequivalent. Evaporation rate and fraction of volatile components showed some sensitivity, suggesting that large changes in the volume or composition of the volatile fraction could cause bio-in-equivalence. The tested dosing scenarios showed decreased sensitivity for all formulation parameters with a decreased dose. The relative formulation bioequivalence was insensitive to vasoconstriction, but the safe space became wider with decreased barrier function for all parameters, except viscosity that was unaffected.