Proof of Concept in Assignment of Within-Subject Variability During Virtual Bioequivalence Studies: Propagation of Intra-Subject Variation in Gastrointestinal Physiology Using Physiologically Based Pharmacokinetic Modeling

Integrating Clinical Variability into PBPK Models for Virtual Bioequivalence of Single and Multiple Doses of Tofacitinib Modified-Release Dosage Form

Tofacitinib is a potent, selective inhibitor of the Janus kinase (JAK) family of kinases with a high degree of selectivity within the human genome's set of protein kinases. Currently approved formulations for tofacitinib citrate are immediate-release (IR) tablets, modified-release (MR) tablets, and IR solution. A once daily MR microsphere formulation was developed for use in pediatric patients. Demonstration of bioequivalence (BE) between the 10 mg once daily (q.d.) MR microsphere formulation and 5 mg twice daily (b.i.d.) IR solution is needed to enable the exposure-response analyses-based bridging to support regulatory approval. To assess BE between MR microsphere and IR solution, an innovative approach was utilized with physiologically-based pharmacokinetic (PBPK) virtual BE trials (VBE) in lieu of a clinical BE trial. A PBPK model was developed to characterize the absorption of different formulations of tofacitinib using Simcyp ADAM module. VBE trials were conducted by simulating PK profiles using the verified PBPK model and integrating the clinically observed intrasubject coefficient of variation (ICV) where BE was assessed with a predetermined sample size and prespecified criteria. The VBE trials demonstrated BE between IR solution 5 mg b.i.d. and MR microsphere 10 mg q.d. after a single dose on day 1 and after multiple doses on day 5. This research presents an innovative approach that incorporates clinically observed ICV in PBPK model-based VBE trials, which could reduce unnecessary drug exposure to healthy volunteers and streamline new formulation development strategies.

Virtual Bioequivalence Assessment of Ritlecitinib Capsules with Incorporation of Observed Clinical Variability Using a Physiologically Based Pharmacokinetic Model

Ritlecitinib, an orally available Janus kinase 3 and tyrosine kinase inhibitor being developed for the treatment of alopecia areata (AA), is highly soluble across the physiological pH range at the therapeutic dose. As such, it is expected to dissolve rapidly in any in vitro dissolution conditions. However, in vitro dissolution data showed slower dissolution for 100-mg capsules, used for the clinical bioequivalence (BE) study, compared with proposed commercial 50-mg capsules. Hence, a biowaiver for the lower 50-mg strength using comparable multimedia dissolution based on the f2 similarity factor was not possible. The in vivo relevance of this observed in vitro dissolution profile was evaluated with a physiologically based pharmacokinetic (PBPK) model. This report describes the development, verification, and application of the ritlecitinib PBPK model to translate observed in vitro dissolution data to an in vivo PK profile for ritlecitinib capsule formulations. Virtual BE (VBE) trials were conducted using the Simcyp VBE module, including the model-predicted within-subject variability or intra-subject coefficient of variation (ICV). The results showed the predicted ICV was predicted to be smaller than observed clinical ICV, resulting in a more optimistic BE risk assessment. Additional VBE assessment was conducted by incorporating clinically observed ICV. The VBE trial results including clinically observed ICV demonstrated that proposed commercial 50-mg capsules vs clinical 100-mg capsules were bioequivalent, with > 90% probability of success. This study demonstrates a PBPK model-based biowaiver for a clinical BE study while introducing a novel method to integrate clinically observed ICV into VBE trials with PBPK models. Trial registration: NCT02309827, NCT02684760, NCT04004663, NCT04390776, NCT05040295, NCT05128058.

A Rational Approach to Predicting Immediate Release Formulation Behavior in Multiple Gastric Motility Patterns: A Combination of a Biorelevant Apparatus, Design of Experiments, and Machine Learning

Gastric mechanical stress often impacts drug dissolution from solid oral dosage forms, but in vitro experiments cannot recreate the substantial variability of gastric motility in a reasonable time. This study, for the first time, combines a novel dissolution apparatus with the design of experiments (DoE) and machine learning (ML) to overcome this obstacle. The workflow involves the testing of soft gelatin capsules in a set of fasted-state biorelevant dissolution experiments created with DoE. The dissolution results are used by an ML algorithm to build the classification model of the capsule's opening in response to intragastric stress (IS) within the physiological space of timing and magnitude. Next, a random forest algorithm is used to model the further drug dissolution. The predictive power of the two ML models is verified with independent dissolution tests, and they outperform a polynomial-based DoE model. Moreover, the developed tool reasonably simulates over 50 dissolution profiles under varying IS conditions. Hence, we prove that our method can be utilized for the simulation of dissolution profiles related to the multiplicity of individual gastric motility patterns. In perspective, the developed workflow can improve virtual bioequivalence trials and the patient-centric development of immediate-release oral dosage forms.

Integration of Biorelevant Pediatric Dissolution Methodology into PBPK Modeling to Predict In Vivo Performance and Bioequivalence of Generic Drugs in Pediatric Populations: a Carbamazepine Case Study

This study investigated the impact of gastro-intestinal fluid volume and bile salt (BS) concentration on the dissolution of carbamazepine (CBZ) immediate release (IR) 100 mg tablets and to integrate these in vitro biorelevant dissolution profiles into physiologically based pharmacokinetic modelling (PBPK) in pediatric and adult populations to determine the biopredictive dissolution profile. Dissolution profiles of CBZ IR tablets (100 mg) were generated in 50-900 mL biorelevant adult fasted state simulated gastric and intestinal fluid (Ad-FaSSGF and Ad-FaSSIF), also in three alternative compositions of biorelevant pediatric FaSSGF and FaSSIF medias at 200 mL. This study found that CBZ dissolution was poorly sensitive to changes in the composition of the biorelevant media, where dissimilar dissolution (F2 = 46.2) was only observed when the BS concentration was changed from 3000 to 89 μM (Ad-FaSSIF vs Ped-FaSSIF 50% 14 BS). PBPK modeling demonstrated the most predictive dissolution volume and media composition to forecast the PK was 500 mL of Ad-FaSSGF/Ad-FaSSIF media for adults and 200 mL Ped-FaSSGF/FaSSIF media for pediatrics. A virtual bioequivalence simulation was conducted by using Ad-FaSSGF and/or Ad-FaSSIF 500 mL or Ped-FaSSGF and/or Ped-FaSSIF 200 mL dissolution data for CBZ 100 mg (reference and generic test) IR product. The CBZ PBPK models showed bioequivalence of the product. This study demonstrates that the integration of biorelevant dissolution data can predict the PK profile of a poorly soluble drug in both populations. Further work using more pediatric drug products is needed to verify biorelevant dissolution data to predict the in vivo performance in pediatrics.

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.

Physiologically Based Pharmacokinetics Modeling in Biopharmaceutics: Case Studies for Establishing the Bioequivalence Safe Space for Innovator and Generic Drugs

For successful oral drug development, defining a bioequivalence (BE) safe space is critical for the identification of newer bioequivalent formulations or for setting of clinically relevant in vitro specifications to ensure drug product quality. By definition, the safe space delineates the dissolution profile boundaries or other drug product quality attributes, within which the drug product variants are anticipated to be bioequivalent. Defining a BE safe space with physiologically based biopharmaceutics model (PBBM) allows the establishment of mechanistic in vitro and in vivo relationships (IVIVR) to better understand absorption mechanism and critical bioavailability attributes (CBA). Detailed case studies on how to use PBBM to establish a BE safe space for both innovator and generic drugs are described. New case studies and literature examples demonstrate BE safe space applications such as how to set in vitro dissolution/particle size distribution (PSD) specifications, widen dissolution specification to supersede f2 tests, or application toward a scale-up and post-approval changes (SUPAC) biowaiver. A workflow for detailed PBBM set-up and common clinical study data requirements to establish the safe space and knowledge space are discussed. Approaches to model in vitro dissolution profiles i.e. the diffusion layer model (DLM), Takano and Johnson models or the fitted PSD and Weibull function are described with a decision tree. The conduct of parameter sensitivity analyses on kinetic dissolution parameters for safe space and virtual bioequivalence (VBE) modeling for innovator and generic drugs are shared. The necessity for biopredictive dissolution method development and challenges with PBBM development and acceptance criteria are described.

Concept and application of the probability of pharmacological success (PoPS) as a decision tool in drug development: a position paper

BACKGROUND

In drug development, few molecules from a large pool of early candidates become successful medicines after demonstrating a favourable benefit-risk ratio. Many decisions are made along the way to continue or stop the development of a molecule. The probability of pharmacological success, or PoPS, is a tool for informing early-stage decisions based on benefit and risk data available at the time.

RESULTS

The PoPS is the probability that most patients can achieve adequate pharmacology for the intended indication while minimising the number of subjects exposed to safety risk. This probability is usually a function of dose; hence its computation typically requires exposure-response models for pharmacology and safety. The levels of adequate pharmacology and acceptable risk must be specified. The uncertainties in these levels, in the exposure-response relationships, and in relevant translation all need to be identified. Several examples of different indications are used to illustrate how this approach can facilitate molecule progression decisions for preclinical and early clinical development. The examples show that PoPS assessment is an effective mechanism for integrating multi-source data, identifying knowledge gaps, and forcing transparency of assumptions. With its application, translational modelling becomes more meaningful and dose prediction more rigorous. Its successful implementation calls for early planning, sound understanding of the disease-drug system, and cross-discipline collaboration. Furthermore, the PoPS evolves as relevant knowledge grows.

CONCLUSION

The PoPS is a powerful evidence-based framework to formally capture multiple uncertainties into a single probability term for assessing benefit-risk ratio. In GSK, it is now expected for governance review at all early-phase decision gates.

Delineating gene-environment effects using virtual twins of patients treated with clozapine

Studies that focus on individual covariates, while ignoring their interactions, may not be adequate for model-informed precision dosing (MIPD) in any given patient. Genetic variations that influence protein synthesis should be studied in conjunction with environmental covariates, such as cigarette smoking. The aim of this study was to build virtual twins (VTs) of real patients receiving clozapine with interacting covariates related to genetics and environment and to delineate the impact of interacting covariates on predicted clozapine plasma concentrations. Clozapine-treated patients with schizophrenia (N = 42) with observed clozapine plasma concentrations, demographic, environmental, and genotype data were used to construct VTs in Simcyp. The effect of increased covariate virtualization was assessed by performing simulations under three conditions: "low" (demographic), "medium" (demographic and environmental interaction), and "high" (demographic and environmental/genotype interaction) covariate virtualization. Increasing covariate virtualization with interaction improved the coefficient of variation (R² ) from 0.07 in the low model to 0.391 and 0.368 in the medium and high models, respectively. Whereas R² was similar between the medium and high models, the high covariate virtualization model had improved accuracy, with systematic bias of predicted clozapine plasma concentration improving from -138.48 ng/ml to -74.65 ng/ml. A high level of covariate virtualization (demographic, environmental, and genotype) may be required for MIPD using VTs.

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.