A Physiologically-Based Pharmacokinetic (PBPK) Model Network for the Prediction of CYP1A2 and CYP2C19 Drug-Drug-Gene Interactions with Fluvoxamine, Omeprazole, S-mephenytoin, Moclobemide, Tizanidine, Mexiletine, Ethinylestradiol, and Caffeine

In-Depth Analysis of the Selection of PBPK Modeling Tools: Bibliometric and Social Network Analysis of the Open Systems Pharmacology Community

Since the Open Source Initiative laid the foundation for the open source software environment in 1998, the popularity of free and open source software has been steadily increasing. Model-informed drug discovery and development (MID3), a key component of pharmaceutical research and development, heavily makes use of computational models which can be developed using various software including the Open Systems Pharmacology (OSP) software (PK-Sim/MoBi), a free and open source software tool for physiologically based pharmacokinetic (PBPK) modeling. In this study, we aimed to investigate the impact, application areas, and reach of the OSP software as well as the relationships and collaboration patterns between organizations having published OSP-related articles between 2017 and 2023. Therefore, we conducted a bibliometric analysis of OSP-related publications and a social network analysis of the organizations with which authors of OSP-related publications were affiliated. On several levels, we found evidence for a significant growth in the size of the OSP community as well as its visibility in the MID3 community since OSP's establishment in 2017. Specifically, the annual publication rate of PubMed-indexed PBPK-related articles using the OSP software outpaced that of PBPK-related articles using any software. Our bibliometric analysis and network analysis demonstrated that the expansion of the OSP community was predominantly driven by new authors and organizations without prior connections to the community involving the generation of research clusters de novo and an overall diversification of the network. These findings suggest an ongoing evolution of the OSP community toward a more segmented, diverse, and inclusive network.

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.

Physiologically Based Pharmacokinetic Modeling to Simulate CYP3A4-Mediated Drug-Drug Interactions for Pyrotinib

INTRODUCTION

Pyrotinib is a newly developed tyrosine kinase inhibitor whose in vivo clearance relies heavily on cytochrome P450 3A4 (CYP3A4) activity. Clinical trials are ongoing to explore the effects of coadministration with CYP3A4 perpetrators on pyrotinib exposure. The present study aims to utilize physiologically based pharmacokinetic (PBPK) modeling to predict CYP3A4-based drug interactions of pyrotinib.

METHODS

Pyrotinib PBPK model was developed in the PK-Sim® multicompartmental physiology structure. Physiochemical parameters were obtained from the literature, and clearance-related parameters were optimized by fitting clinical single-dose pharmacokinetic data. Pharmacokinetic parameters from the model output were compared with the observed data to validate the model predictive performance. Using validated CYP3A4 perpetrator models, we conducted PBPK simulations for drug interactions in a virtual population to explore the impacts of comedication with these perpetrators.

RESULTS

The PBPK model accurately describes pyrotinib single- and multi-dose pharmacokinetics. The model also predicts dramatic exposure change of pyrotinib in the presence of itraconazole and rifampicin, though the impact of rifampicin is somewhat underestimated. According to model predictions, coadministration with typical potent or moderate CYP3A4 perpetrators increases pyrotinib concentration by over sixfold, extinguishing the possibility of dose adjustment for pyrotinib. A weak CYP3A4 inhibitor has minimal influence on pyrotinib pharmacokinetics.

CONCLUSION

PBPK modeling provides valuable information to avoid irrational medication when receiving pyrotinib chemotherapy.

Physiologically-based pharmacokinetic modeling of quinidine to establish a CYP3A4, P-gp, and CYP2D6 drug-drug-gene interaction network

The antiarrhythmic agent quinidine is a potent inhibitor of cytochrome P450 (CYP) 2D6 and P-glycoprotein (P-gp) and is therefore recommended for use in clinical drug-drug interaction (DDI) studies. However, as quinidine is also a substrate of CYP3A4 and P-gp, it is susceptible to DDIs involving these proteins. Physiologically-based pharmacokinetic (PBPK) modeling can help to mechanistically assess the absorption, distribution, metabolism, and excretion processes of a drug and has proven its usefulness in predicting even complex interaction scenarios. The objectives of the presented work were to develop a PBPK model of quinidine and to integrate the model into a comprehensive drug-drug(-gene) interaction (DD(G)I) network with a diverse set of CYP3A4 and P-gp perpetrators as well as CYP2D6 and P-gp victims. The quinidine parent-metabolite model including 3-hydroxyquinidine was developed using pharmacokinetic profiles from clinical studies after intravenous and oral administration covering a broad dosing range (0.1-600 mg). The model covers efflux transport via P-gp and metabolic transformation to either 3-hydroxyquinidine or unspecified metabolites via CYP3A4. The 3-hydroxyquinidine model includes further metabolism by CYP3A4 as well as an unspecific hepatic clearance. Model performance was assessed graphically and quantitatively with greater than 90% of predicted pharmacokinetic parameters within two-fold of corresponding observed values. The model was successfully used to simulate various DD(G)I scenarios with greater than 90% of predicted DD(G)I pharmacokinetic parameter ratios within two-fold prediction success limits. The presented network will be provided to the research community and can be extended to include further perpetrators, victims, and targets, to support investigations of DD(G)Is.

Development and Evaluation of a Physiologically Based Pharmacokinetic Model for Predicting Haloperidol Exposure in Healthy and Disease Populations

The physiologically based pharmacokinetic (PBPK) approach can be used to develop mathematical models for predicting the absorption, distribution, metabolism, and elimination (ADME) of administered drugs in virtual human populations. Haloperidol is a typical antipsychotic drug with a narrow therapeutic index and is commonly used in the management of several medical conditions, including psychotic disorders. Due to the large interindividual variability among patients taking haloperidol, it is very likely for them to experience either toxic or subtherapeutic effects. We intend to develop a haloperidol PBPK model for identifying the potential sources of pharmacokinetic (PK) variability after intravenous and oral administration by using the population-based simulator, PK-Sim. The model was initially developed and evaluated to predict the PK of haloperidol and its reduced metabolite in adult healthy population after intravenous and oral administration. After evaluating the developed PBPK model in healthy adults, it was used to predict haloperidol–rifampicin drug–drug interaction and was extended to tuberculosis patients. The model evaluation was performed using visual assessments, prediction error, and mean fold error of the ratio of the observed-to-predicted values of the PK parameters. The predicted PK values were in good agreement with the corresponding reported values. The effects of the pathophysiological changes and enzyme induction associated with tuberculosis and its treatment, respectively, on haloperidol PK, have been predicted precisely. For all clinical scenarios that were evaluated, the predicted values were within the acceptable two-fold error range.

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

Physiologically‐based pharmacokinetic (PBPK) modeling is being increasingly used in drug development to avoid unnecessary clinical drug–drug interaction (DDI) studies and inform drug labels. Thus, regulatory agencies are recommending, or indeed requesting, more rigorous demonstration of the prediction accuracy of PBPK platforms in the area of their intended use. We describe a framework for qualification of the Simcyp Simulator with respect to competitive and mechanism‐based inhibition (MBI) of CYP1A2, CYP2D6, CYP2C8, CYP2C9, CYP2C19, and CYP3A4/5. Initially, a DDI matrix, consisting of a range of weak, moderate, and strong inhibitors and substrates with varying fraction metabolized by specific CYP enzymes that were susceptible to different degrees of inhibition, were identified. Simulations were run with 123 clinical DDI studies involving competitive inhibition and 78 clinical DDI studies involving MBI. For competitive inhibition, the overall prediction accuracy was good with an average fold error (AFE) of 0.91 and 0.92 for changes in the maximum plasma concentration (C_max) and area under the plasma concentration (AUC) time profile, respectively, as a consequence of the DDI. For MBI, an AFE of 1.03 was determined for both C_max and AUC. The prediction accuracy was generally comparable across all CYP enzymes, irrespective of the isozyme and mechanism of inhibition. These findings provide confidence in application of the Simcyp Simulator (V19 R1) for assessment of the DDI potential of drugs in development either as inhibitors or victim drugs of CYP‐mediated interactions. The approach described herein and the identified DDI matrix can be used to qualify subsequent versions of the platform.

Applications, Challenges, and Outlook for PBPK Modeling and Simulation: A Regulatory, Industrial and Academic Perspective

Several regulatory guidances on the use of physiologically based pharmacokinetic (PBPK) analyses and physiologically based biopharmaceutics model(s) (PBBM(s)) have been issued. Workshops are routinely held, demonstrating substantial interest in applying these modeling approaches to address scientific questions in drug development. PBPK models and PBBMs have remarkably contributed to model-informed drug development (MIDD) such as anticipating clinical PK outcomes affected by extrinsic and intrinsic factors in general and specific populations. In this review, we proposed practical considerations for a "base" PBPK model construction and development, summarized current status, challenges including model validation and gaps in system models, and future perspectives in PBPK evaluation to assess a) drug metabolizing enzyme(s)- or drug transporter(s)- mediated drug-drug interactions b) dosing regimen prediction, sampling timepoint selection and dose validation in pediatric patients from newborns to adolescents, c) drug exposure in patients with renal and/or and hepatic organ impairment, d) maternal-fetal drug disposition during pregnancy, and e) pH-mediated drug-drug interactions in patients treated with proton pump inhibitors/acid-reducing agents (PPIs/ARAs) intended for gastric protection. Since PBPK can simulate outcomes in clinical studies with enrollment challenges or ethical issues, the impact of PBPK models on waivers and how to strengthen study waiver is discussed.

Physiologically Based Pharmacokinetic (PBPK) Modeling of Clopidogrel and Its Four Relevant Metabolites for CYP2B6, CYP2C8, CYP2C19, and CYP3A4 Drug–Drug–Gene Interaction Predictions

The antiplatelet agent clopidogrel is listed by the FDA as a strong clinical index inhibitor of cytochrome P450 (CYP) 2C8 and weak clinical inhibitor of CYP2B6. Moreover, clopidogrel is a substrate of—among others—CYP2C19 and CYP3A4. This work presents the development of a whole-body physiologically based pharmacokinetic (PBPK) model of clopidogrel including the relevant metabolites, clopidogrel carboxylic acid, clopidogrel acyl glucuronide, 2-oxo-clopidogrel, and the active thiol metabolite, with subsequent application for drug–gene interaction (DGI) and drug–drug interaction (DDI) predictions. Model building was performed in PK-Sim^® using 66 plasma concentration-time profiles of clopidogrel and its metabolites. The comprehensive parent-metabolite model covers biotransformation via carboxylesterase (CES) 1, CES2, CYP2C19, CYP3A4, and uridine 5′-diphospho-glucuronosyltransferase 2B7. Moreover, CYP2C19 was incorporated for normal, intermediate, and poor metabolizer phenotypes. Good predictive performance of the model was demonstrated for the DGI involving CYP2C19, with 17/19 predicted DGI AUC_last and 19/19 predicted DGI C_max ratios within 2-fold of their observed values. Furthermore, DDIs involving bupropion, omeprazole, montelukast, pioglitazone, repaglinide, and rifampicin showed 13/13 predicted DDI AUC_last and 13/13 predicted DDI C_max ratios within 2-fold of their observed ratios. After publication, the model will be made publicly accessible in the Open Systems Pharmacology repository.

Current Practices, Gap Analysis, and Proposed Workflows for PBPK Modeling of Cytochrome P450 Induction: An Industry Perspective

The International Consortium for Innovation and Quality (IQ) Physiologically Based Pharmacokinetic (PBPK) Modeling Induction Working Group (IWG) conducted a survey across participating companies around general strategies for PBPK modeling of induction, including experience with its utility to address various questions, regulatory interactions, and regulatory acceptance. The results highlight areas where PBPK modeling is used with high confidence and identifies opportunities where confidence is lower and further evaluation is needed. To enhance the survey results, the PBPK-IWG also collected case studies and analyzed recent literature examples where PBPK models were applied to predict CYP3A induction-mediated drug-drug interactions. PBPK modeling of induction has evolved and progressed significantly, proving to have great potential to accelerate drug discovery and development. With the aim of enabling optimal use for new molecular entities that are either substrates and/or inducers of CYP3A, the PBPK-IWG proposes initial workflows for PBPK application, discusses future trends, and identifies gaps that need to be addressed.

Predictive In Vitro-In Vivo Extrapolation for Time Dependent Inhibition of CYP1A2, CYP2C8, CYP2C9, CYP2C19 and CYP2D6 Using Pooled Human Hepatocytes, Human Liver Microsomes, and a Simple Mechanistic Static Model

Inactivation of Cytochrome P450 (CYP450) enzymes can lead to significant increases in exposure of co-medicants. The majority of reported in vitro to in vivo extrapolation (IVIVE) data have historically focused on CYP3A4 leaving the assessment of other CYP isoforms insubstantial. To this end, the utility of human hepatocytes (HHEP) and microsome (HLM) to predict clinically relevant DDIs was investigated with a focus on CYP1A2, CYP2C8, CYP2C9, CYP2C19 and CYP2D6. Evaluation of IVIVE for CYP2B6 was limited to only weak inhibition. A search of the University of Washington Drug-Drug Interaction Database was conducted to identify a clinically relevant weak, moderate and strong inhibitor for selective substrates of CYP1A2, CYP2C8, CYP2C9, CYP2C19 and CYP2D6, resulting in 18 inhibitors for in vitro characterization against 119 clinical interaction studies. Pooled human hepatocytes and HLM were pre-incubated with increasing concentrations of inhibitors for designated timepoints. Time dependent inhibition (TDI) was detected in HLM for four moderate/strong inhibitors suggesting that some optimization of incubation conditions (i.e. lower protein concentrations) is needed to capture weak inhibition. Clinical risk assessment was conducted by incorporating the in vitro derived kinetic parameters k_inact and K_I into static equations recommended by regulatory authorities. Significant overprediction was observed when applying the basic models recommended by regulatory agencies. Mechanistic static models (MSM), which consider the fraction of metabolism through the impacted enzyme, using the unbound hepatic inlet concentration lead to the best overall prediction accuracy with 92% and 85% of data from HHEPs and HLM, respectively, within 2-fold of the observed value. Significance Statement Collectively, the data demonstrate that coupling time-dependent inactivation parameters derived from pooled human hepatocytes and HLM with a mechanistic static model provides an easy and quantitatively accurate means to determine clinical DDI risk from in vitro data. Weak and moderate inhibitors did not show TDI under standard incubation conditions using HLM and optimization of incubation conditions is warranted. Recommendations are made with respect to input parameters for IVIVE of TDI with non-CYP3A enzymes using available data from HLM and HHEPs.

Evaluation of the drug-drug interaction potential of treosulfan using a physiologically-based pharmacokinetic modelling approach

AIMS

The aim of this work is the development of a mechanistic physiologically-based pharmacokinetic (PBPK) model using in vitro to in vivo extrapolation to conduct a drug-drug interaction (DDI) assessment of treosulfan against two cytochrome p450 (CYP) isoenzymes and P-glycoprotein (P-gp) substrates.

METHODS

A PBPK model for treosulfan was developed de novo based on literature and unpublished clinical data. The PBPK DDI analysis was conducted using the U.S. Food and Drug Administration (FDA) DDI index drugs (probe substrates) midazolam, omeprazole and digoxin for CYP3A4, CYP2C19 and P-gp, respectively. Qualified and documented PBPK models of the probe substrates have been adopted from an open-source online model database.

RESULTS

The PBPK model for treosulfan, based on both in vitro and in vivo data, was able to predict the plasma concentration-time profiles and exposure levels of treosulfan applied for a standard conditioning treatment. Medium and low potentials for DDI on CYP3A4 (maximum area under the concentration-time curve ratio (AUCR_max  = 2.23) and CYP2C19 (AUCR_max  = 1.6) were predicted, respectively, using probe substrates midazolam and omeprazole. Treosulfan was not predicted to cause a DDI on P-gp.

CONCLUSION

Medicinal products with a narrow therapeutic index (eg, digoxin) that are substrates for CYP3A4, CYP2C19 or P-gp should not be given during treatment with treosulfan. However, considering the comprehensive treosulfan-based conditioning treatment schedule and the respective pharmacokinetic properties of the concomitantly used drugs (eg, half-life), the potential for interaction on all evaluated mechanisms would be low (AUCR < 1.25), if concomitantly administered drugs are dosed either 2 hours before or 8 hours after the 2-hour intravenous infusion of treosulfan.

Drug Interactions Involving 17α-Ethinylestradiol: Considerations Beyond Cytochrome P450 3A Induction and Inhibition

It is widely acknowledged that drug-drug interactions (DDIs) involving estrogen (17α-ethinylestradiol (EE))-containing oral contraceptives (OCs) are important. Consequently, sponsors of new molecular entities (NMEs) often conduct clinical studies with priority given to OCs as victims of cytochrome P450 (CYP) 3A (CYP3A) induction and inhibition. Such scenarios are reflected in the US Food and Drug Administration-issued guidance documentation related to OC DDI studies. Although CYP3A is important, OCs such as EE are metabolized by sulfotransferase 1E1 and UDP-glucuronosyltransferase (UGT) 1A1, expressed in the gut and liver, and so both can also serve as loci of victim OC DDI. Therefore, for any NME, one should carefully consider its induction and inhibition profile involving CYP3A4/5, UGT1A1, and SULT1E1. As DDI perpetrators, available clinical DDI data indicate that EE-containing OCs can induce (e.g., UGT1A4 and CYP2A6) and inhibit (CYP1A2 ≥ CYP2C19 > CYP3A4/5 > CYP2C8, CYP2B6, CYP2D6, and CYP2C9) various CYP forms. Although available in vitro CYP inhibition data do not explain such a graded inhibitory effect in vivo, it is hypothesized that EE differentially modulates CYP expression via potent agonism of the estrogen receptor expressed in the gut and liver. From the standpoint of the NME as potential OC DDI victim, therefore, it is important to assess its projected (pre-phase I) or known therapeutic index and pharmacokinetic profile (fraction absorbed, absolute oral bioavailability, clearance/extraction class, fraction metabolized by CYP1A2, CYP2C19, CYP2A6, and UGT1A4). Such information can enable the prioritization, design, and interpretation of NME-OC DDI studies.

A generic framework for the physiologically-based pharmacokinetic platform qualification of PK-Sim and its application to predicting cytochrome P450 3A4-mediated drug-drug interactions

The success of applications of physiologically‐based pharmacokinetic (PBPK) modeling in drug development and drug labeling has triggered regulatory agencies to demand rigorous demonstration of the predictive capability of the specific PBPK platform for a particular intended application purpose. The effort needed to comply with such qualification requirements exceeds the costs for any individual PBPK application. Because changes or updates of a PBPK platform would require (re‐)qualification, a reliable and efficient generic qualification framework is needed. We describe the development and implementation of an agile and sustainable technical framework for automatic PBPK platform (re‐)qualification of PK‐Sim^® embedded in the open source and open science GitHub landscape of Open Systems Pharmacology. The qualification approach enables the efficient assessment of all aspects relevant to the qualification of a particular purpose and provides transparency and traceability for all stakeholders. As a showcase example for the power and versatility of the qualification framework, we present the qualification of PK‐Sim^® for the intended purpose of predicting cytochrome P450 3A4 (CYP3A4)–mediated drug–drug interactions (DDIs). Several perpetrator PBPK models featuring various degrees of CYP3A4 modulation and different types of mechanisms (competitive inhibition, mechanism‐based inactivation, and induction) were coupled with a set of PBPK models of sensitive CYP3A4 victim drugs. Simulations were compared to a comprehensive data set of 135 observations from published clinical DDI studies. The platform's overall predictive performance showed reasonable accuracy and precision (geometric mean fold error of 1.4 for both area under the plasma concentration‐time curve ratios and peak plasma concentration ratios with/without perpetrator) and suggests that PK‐Sim^® can be applied to quantitatively assess CYP3A4‐mediated DDI in clinically untested scenarios.

Disease-drug and drug-drug interaction in COVID-19: Risk and assessment

COVID-19 is announced as a global pandemic in 2020. Its mortality and morbidity rate are rapidly increasing, with limited medications. The emergent outbreak of COVID-19 prompted by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) keeps spreading. In this infection, a patient's immune response plays pivotal role in the pathogenesis. This inflammatory factor was shown by its mediators that, in severe cases, reach the cytokine at peaks. Hyperinflammatory state may sparks significant imbalances in transporters and drug metabolic machinery, and subsequent alteration of drug pharmacokinetics may result in unexpected therapeutic response. The present scenario has accounted for the requirement for therapeutic opportunities to relive and overcome this pandemic. Despite the diminishing developments of COVID-19, there is no drug still approved to have significant effects with no side effect on the treatment for COVID-19 patients. Based on the evidence, many antiviral and anti-inflammatory drugs have been authorized by the Food and Drug Administration (FDA) to treat the COVID-19 patients even though not knowing the possible drug-drug interactions (DDI). Remdesivir, favipiravir, and molnupiravir are deemed the most hopeful antiviral agents by improving infected patient’s health. Dexamethasone is the first known steroid medicine that saved the lives of seriously ill patients. Some oligopeptides and proteins have also been using. The current review summarizes medication updates to treat COVID-19 patients in an inflammatory state and their interaction with drug transporters and drug-metabolizing enzymes. It gives an opinion on the potential DDI that may permit the individualization of these drugs, thereby enhancing the safety and efficacy.