Prediction of CYP‐mediated DDIs involving inhibition: Approaches to address the requirements for system qualification of the Simcyp Simulator

Evaluating drug interaction potential from cytokine release syndrome using a physiologically based pharmacokinetic model: A case study of teclistamab

Cytokine release syndrome (CRS) was associated with teclistamab treatment in the phase I/II MajesTEC-1 study. Cytokines, especially interleukin (IL)-6, are known suppressors of cytochrome P450 (CYP) enzymes' activity. A physiologically based pharmacokinetic model evaluated the impact of IL-6 serum levels on exposure of substrates of various CYP enzymes (1A2, 2C9, 2C19, 3A4, 3A5). Two IL-6 kinetics profiles were assessed, the mean IL-6 profile with a maximum concentration (Cmax) of IL-6 (21 pg/mL) and the IL-6 profile of the patient presenting the highest IL-6 Cmax (288 pg/mL) among patients receiving the recommended phase II dose of teclistamab in MajesTEC-1. For the mean IL-6 kinetics profile, teclistamab was predicted to result in a limited change in exposure of CYP substrates (area under the curve [AUC] mean ratio 0.87-1.20). For the maximum IL-6 kinetics profile, the impact on omeprazole, simvastatin, midazolam, and cyclosporine exposure was weak to moderate (mean AUC ratios 1.90-2.23), and minimal for caffeine and s-warfarin (mean AUC ratios 0.82-1.25). Maximum change in exposure for these substrates occurred 3-4 days after step-up dosing in cycle 1. These results suggest that after cycle 1, drug interaction from IL-6 effect has no meaningful impact on CYP activities, with minimal or moderate impact on CYP substrates. The highest risk of drug interaction is expected to occur during step-up dosing up to 7 days after the first treatment dose (1.5 mg/kg subcutaneously) and during and after CRS.

Prediction of Pediatric Pharmacokinetics for CYP3A4 Metabolized Drugs: Comparison of the Performance of Two Hepatic Ontogeny Within a Physiologically Based Pharmacokinetic Model

The rapid growth in the use of pediatric physiologically based pharmacokinetic (PBPK) models, particularly for regulatory applications, has focused emphasis on model verification and ensuring system parameters are robust, including how these change with age. Uncertainty remains regarding the ontogeny of some enzymes and transporters, in this study 2 published ontogeny profiles for hepatic CYP3A4 were compared. Clinical pharmacokinetic data on 4 intravenously administered CYP3A4 substrates (alfentanil, fentanyl, midazolam, and sildenafil) used across the pediatric age range was collected from the literature. The PBPK models were verified in the adult population and then used to compare the Salem and a modified Upreti ontogeny profiles for CYP3A4 in terms of parent drug clearance and area under the curve from birth onward. Overall, the modified Upreti ontogeny profile resulted in 15 out of 17 age-related predictions within 2-fold and 12 out of 17 predictions within 1.5-fold ranges of observed values, for the Salem ontogeny these values were 12 out of 17 and 8 out of 17, respectively. The Upreti ontogeny profile performed better than Salem, average fold error and absolute average fold error were 1.14 and 1.35 compared to 1.56 and 1.90, respectively. Identifying the optimal CYP3A4 ontogeny is important for regulatory use of PBPK especially given the number of drugs cleared by this enzyme. This study broadens the evidence from previous studies that Upreti is more favorable than Salem, but further work is needed especially in the neonatal and early infant age range.

Dose Optimization Informed by PBPK Modeling: State-of-the Art and Future

Model-informed drug development (MIDD) is a powerful quantitative approach that plays an integral role in drug development and regulatory review. While applied throughout the life cycle of the development of new drugs, a key application of MIDD is to inform clinical trial design including dose selection and optimization. To date, physiologically-based pharmacokinetic (PBPK) modeling, an established component of the MIDD toolkit, has mainly been used for assessment of drug-drug interactions (DDIs) and consequential dose adjustments in regulatory submissions. As a result of recent scientific advances and growing confidence in the utility of the approach, PBPK models are being increasingly applied to provide dose recommendations for subjects with differing ages, genetics, and disease states. In this review, we present our perspective on the current landscape of regulatory acceptance of PBPK applications supported by relevant case studies. We also discuss the recent progress and future challenges associated with expanding the utility of PBPK models into emerging areas for regulatory decision making, especially dose optimization in highly vulnerable and understudied populations and facilitating diversity in clinical trials.

Applied physiologically-based pharmacokinetic modeling to assess uridine diphosphate-glucuronosyltransferase-mediated drug-drug interactions for Vericiguat

Vericiguat (Verquvo; US: Merck, other countries: Bayer) is a novel drug for the treatment of chronic heart failure. Preclinical studies have demonstrated that the primary route of metabolism for vericiguat is glucuronidation, mainly catalyzed by uridine diphosphate-glucuronosyltransferase (UGT)1A9 and to a lesser extent UGT1A1. Whereas a drug-drug interaction (DDI) study of the UGT1A9 inhibitor mefenamic acid showed a 20% exposure increase, the effect of UGT1A1 inhibitors has not been assessed clinically. This modeling study describes a physiologically-based pharmacokinetic (PBPK) approach to complement the clinical DDI liability assessment and support prescription labeling. A PBPK model of vericiguat was developed based on in vitro and clinical data, verified against data from the mefenamic acid DDI study, and applied to assess the UGT1A1 DDI liability by running an in silico DDI study with the UGT1A1 inhibitor atazanavir. A minor effect with an area under the plasma concentration-time curve (AUC) ratio of 1.12 and a peak plasma concentration ratio of 1.04 was predicted, which indicates that there is no clinically relevant DDI interaction anticipated. Additionally, the effect of potential genetic polymorphisms of UGT1A1 and UGT1A9 was evaluated, which showed that an average modest increase of up to 1.7-fold in AUC may be expected in the case of concomitantly reduced UGT1A1 and UGT1A9 activity for subpopulations expressing non-wild-type variants for both isoforms. This study is a first cornerstone to qualify the PK-Sim platform for use of UGT-mediated DDI predictions, including PBPK models of perpetrators, such as mefenamic acid and atazanavir, and sensitive UGT substrates, such as dapagliflozin and raltegravir.

Challenges and Opportunities for In Vitro-In Vivo Extrapolation of Aldehyde Oxidase-Mediated Clearance: Toward a Roadmap for Quantitative Translation

Underestimation of aldehyde oxidase (AO)-mediated clearance by current in vitro assays leads to uncertainty in human dose projections, thereby reducing the likelihood of success in drug development. In the present study we first evaluated the current drug development practices for AO substrates. Next, the overall predictive performance of in vitro-in vivo extrapolation of unbound hepatic intrinsic clearance (CLint,u) and unbound hepatic intrinsic clearance by AO (CLint,u,AO) was assessed using a comprehensive literature database of in vitro (human cytosol/S9/hepatocytes) and in vivo (intravenous/oral) data collated for 22 AO substrates (total of 100 datapoints from multiple studies). Correction for unbound fraction in the incubation was done by experimental data or in silico predictions. The fraction metabolized by AO (fmAO) determined via in vitro/in vivo approaches was found to be highly variable. The geometric mean fold errors (gmfe) for scaled CLint,u (mL/min/kg) were 10.4 for human hepatocytes, 5.6 for human liver cytosols, and 5.0 for human liver S9, respectively. Application of these gmfe's as empirical scaling factors improved predictions (45%-57% within twofold of observed) compared with no correction (11%-27% within twofold), with the scaling factors qualified by leave-one-out cross-validation. A road map for quantitative translation was then proposed following a critical evaluation on the in vitro and clinical methodology to estimate in vivo fmAO In conclusion, the study provides the most robust system-specific empirical scaling factors to date as a pragmatic approach for the prediction of in vivo CLint,u,AO in the early stages of drug development. SIGNIFICANCE STATEMENT: Confidence remains low when predicting in vivo clearance of AO substrates using in vitro systems, leading to de-prioritization of AO substrates from the drug development pipeline to mitigate risk of unexpected and costly in vivo impact. The current study establishes a set of empirical scaling factors as a pragmatic tool to improve predictability of in vivo AO clearance. Developing clinical pharmacology strategies for AO substrates by utilizing mass balance/clinical drug-drug interaction data will help build confidence in fmAO.

A comprehensive evaluation in clinic and physiologically-based pharmacokinetic modeling and simulation to confirm lack of cytochrome P450-mediated drug-drug interaction potential for pomotrelvir

Pomotrelvir is a new chemical entity and potent direct-acting antiviral inhibitor of the main protease of coronaviruses. Here the cytochrome P450 (CYP)-mediated drug-drug interaction (DDI) potential of pomotrelvir was evaluated for major CYP isoforms, starting with in vitro assays followed by the basic static model assessment. The identified CYP3A4-mediated potential DDIs were evaluated clinically at a supratherapeutic dose of 1050 mg twice daily (b.i.d.) of pomotrelvir, including pomotrelvir coadministration with ritonavir (strong inhibitor of CYP3A4) or midazolam (sensitive substrate of CYP3A4). Furthermore, a physiologically-based pharmacokinetic (PBPK) model was developed within the Simcyp Population-based Simulator using in vitro and in vivo information and validated with available human pharmacokinetic (PK) data. The PBPK model was simulated to assess the DDI potential for CYP isoforms that pomotrelvir has shown a weak to moderate DDI in vitro and for CYP3A4 at the therapeutic dose of 700 mg b.i.d. To support the use of pomotrelvir in women of childbearing potential, the impact of pomotrelvir on the exposure of the representative oral hormonal contraceptive drugs ethinyl estradiol and levonorgestrel was assessed using the PBPK model. The overall assessment suggested weak inhibition of pomotrelvir on CYP3A4 and minimal impact of a strong CYP3A4 inducer or inhibitor on pomotrelvir PK. Therefore, pomotrelvir is not anticipated to have clinically meaningful DDIs at the clinical dose. These comprehensive in vitro, in clinic, and in silico efforts indicate that the DDI potential of pomotrelvir is minimal, so excluding patients on concomitant medicines in clinical studies would not be required.

Development and application of a PBPK modeling strategy to support antimalarial drug development

As part of a collaboration between Medicines for Malaria Venture (MMV), Certara UK and Monash University, physiologically-based pharmacokinetic (PBPK) models were developed for 20 antimalarials, using data obtained from standardized in vitro assays and clinical studies within the literature. The models have been applied within antimalarial drug development at MMV for more than 5 years. During this time, a strategy for their impactful use has evolved. All models are described in the supplementary material and are available to researchers. Case studies are also presented, demonstrating real-world development and clinical applications, including the assessment of the drug-drug interaction liability between combination partners or with co-administered drugs. This work emphasizes the benefit of PBPK modeling for antimalarial drug development and decision making, and presents a strategy to integrate it into the research and development process. It also provides a repository of shared information to benefit the global health research community.

A novel physiologically based pharmacokinetic model of rectal absorption, evaluated and verified using clinical data on 10 rectally administered drugs

We present a physiologically based pharmacokinetic (PBPK) model simulating systemic drug concentrations following administration to the human rectum. Rectum physiology is parameterized based on literature data. The model utilizes in vitro release (IVRT) profiles from which drug mass transfer through the rectal fluid and tissue and into the systemic circulation are predicted. Due to a lack of data, rectal fluid and tissue absorption parameters are predicted either from colon absorption, with modifications relevant to rectal physiology, or optimized. The PBPK model is evaluated by simulating 29 clinical studies for 10 drugs. For 8 drugs (diazepam, diclofenac, indomethacin, naproxen, paracetamol, pentobarbital, phenobarbital and theophylline) the bias (average fold error, AFE) and precision (absolute average fold error, AAFE) of Cmax, AUC0-t and AUC0-inf simulations range from 0.87 to 2.22, indicating good agreement with observed values. For prochlorperazine and promethazine, the AFEs and AAFEs of Cmax predictions are 1.30 and 2.52, respectively. TheAUC0-t and AUC0-inf are overpredicted for both compounds(AFEs and AAFEs from 2.66 to 4.90). This results from a lack of reliable elimination data for prochlorperazine and the relevance of the IVRT profiles used in the promethazine model. The model paves the way for more mechanistic rectal drug absorption studies and virtual bioequivalence methods for rectal drug products.

Clinical Pharmacology of Brigatinib: A Next-Generation Anaplastic Lymphoma Kinase Inhibitor

Brigatinib, a next-generation anaplastic lymphoma kinase (ALK) inhibitor designed to overcome mechanisms of resistance associated with crizotinib, is approved for the treatment of ALK-positive advanced or metastatic non-small cell lung cancer. After oral administration of single doses of brigatinib 30-240 mg, the median time to reach maximum plasma concentration ranged from 1 to 4 h. In patients with advanced malignancies, brigatinib showed dose linearity over the dose range of 60-240 mg once daily. A high-fat meal had no clinically meaningful effect on systemic exposures of brigatinib (area under the plasma concentration-time curve); thus, brigatinib can be administered with or without food. In a population pharmacokinetic analysis, a three-compartment pharmacokinetic model with transit absorption compartments was found to adequately describe brigatinib pharmacokinetics. In addition, the population pharmacokinetic analyses showed that no dose adjustment is required based on body weight, age, race, sex, total bilirubin (< 1.5× upper limit of normal), and mild-to-moderate renal impairment. Data from dedicated phase I trials have indicated that no dose adjustment is required for patients with mild or moderate hepatic impairment, while a dose reduction of approximately 40% (e.g., from 180 to 120 mg) is recommended for patients with severe hepatic impairment, and a reduction of approximately 50% (e.g., from 180 to 90 mg) is recommended when administering brigatinib to patients with severe renal impairment. Brigatinib is primarily metabolized by cytochrome P450 (CYP) 3A, and results of clinical drug-drug interaction studies and physiologically based pharmacokinetic analyses have demonstrated that coadministration of strong or moderate CYP3A inhibitors or inducers with brigatinib should be avoided. If coadministration with a strong or moderate CYP3A inhibitor cannot be avoided, the dose of brigatinib should be reduced by approximately 50% (strong CYP3A inhibitor) or approximately 40% (moderate CYP3A inhibitor), respectively. Brigatinib is a weak inducer of CYP3A in vivo; data from a phase I drug-drug interaction study showed that coadministration of brigatinib 180 mg once daily reduced the oral midazolam area under the plasma concentration-time curve from time zero to infinity by approximately 26%. Brigatinib did not inhibit CYP1A2, CYP2B6, CYP2C8, CYP2C9, CYP2C19, or CYP2D6 at clinically relevant concentrations in vitro. Exposure-response analyses based on data from the ALTA (ALK in Lung Cancer Trial of AP26113) and ALTA-1L pivotal trials of brigatinib confirm the favorable benefit versus risk profile of the approved titration dosing regimen of 180 mg once daily (after a 7-day lead-in at 90 mg once daily).

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.

Comparing the applications of machine learning, PBPK, and population pharmacokinetic models in pharmacokinetic drug-drug interaction prediction

The gold-standard approach for modeling pharmacokinetic mediated drug-drug interactions is the use of physiologically-based pharmacokinetic modeling and population pharmacokinetics. However, these models require extensive amounts of drug-specific data generated from a wide variety of in vitro and in vivo models, which are later refined with clinical data and system-specific parameters. Machine learning has the potential to be utilized for the prediction of drug-drug interactions much earlier in the drug discovery cycle, using inputs derived from, among others, chemical structure. This could lead to refined chemical designs in early drug discovery. Machine-learning models have many advantages, such as the capacity to automate learning (increasing the speed and scalability of predictions), improved generalizability by learning from multicase historical data, and highlighting statistical and potentially clinically significant relationships between input variables. In contrast, the routinely used mechanistic models (physiologically-based pharmacokinetic models and population pharmacokinetics) are currently considered more interpretable, reliable, and require a smaller sample size of data, although insights differ on a case-by-case basis. Therefore, they may be appropriate for later stages of drug-drug interaction assessment when more in vivo and clinical data are available. A combined approach of using mechanistic models to highlight features that can be used for training machine-learning models may also be exploitable in the future to improve the performance of machine learning. In this review, we provide concepts, strategic considerations, and compare machine learning to mechanistic modeling for drug-drug interaction risk assessment across the stages of drug discovery and development.