Comprehensive PBPK Model of Rifampicin for Quantitative Prediction of Complex Drug-Drug Interactions: CYP3A/2C9 Induction and OATP Inhibition Effects

Integrated PBPK-EO modeling of osimertinib to predict plasma concentrations and intracranial EGFR engagement in patients with brain metastases

The purpose of this study was to develop and validate a physiologically based pharmacokinetic (PBPK) model combined with an EGFR occupancy (EO) model for osimertinib (OSI) to predict plasma trough concentration (Ctrough) and the intracranial time-course of EGFR (T790M and L858R mutants) engagement in patient populations. The PBPK model was also used to investigate the key factors affecting OSI pharmacokinetics (PK) and intracranial EGFR engagement, analyze resistance to the target mutation C797S, and determine optimal dosing regimens when used alone and in drug-drug interactions (DDIs). A population PBPK-EO model of OSI was developed using physicochemical, biochemical, binding kinetic, and physiological properties, and then validated using nine clinical PK studies, observed EO study, and two clinical DDI studies. The PBPK-EO model demonstrated good consistency with observed data, with most prediction-to-observation ratios falling within the range of 0.7 to 1.3 for plasma AUC, Cmax, Ctrough and intracranial free concentration. The simulated time-course of C797S occupancy by the PBPK model was much lower than T790M and L858R occupancy, providing an explanation for OSI on-target resistance to the C797S mutation. The PBPK model identified ABCB1 CLint,u, albumin level, and EGFR expression as key factors affecting plasma Ctrough and intracranial EO for OSI. Additionally, PBPK-EO simulations indicated that the optimal dosing regimen for OSI in patients with brain metastases is either 80 mg once daily (OD) or 160 mg OD, or 40 mg or 80 mg twice daily (BID). When used concomitantly with CYP enzyme perpetrators, the PBPK-EO model suggested appropriate dosing regimens of 80 mg OD with fluvoxamine (FLUV) itraconazole (ITR) or fluvoxamine (FLUC) for co-administration and an increase to 160 mg OD with rifampicin (RIF) or efavirenz (EFA). In conclusion, the PBPK-EO model has been shown to be capable of simulating the pharmacokinetic concentration-time profiles and the time-course of EGFR engagement for OSI, as well as determining the optimum dosing in various clinical situations.

Beyond the basics: A deep dive into parameter estimation for advanced PBPK and QSP models

Physiologically-based pharmacokinetic (PBPK) models and quantitative systems pharmacology (QSP) models have contributed to drug development strategies. The parameters of these models are commonly estimated by capturing observed values using the nonlinear least-squares method. Software packages for PBPK and QSP modeling provide a range of parameter estimation algorithms. To choose the most appropriate method, modelers need to understand the basic concept of each approach. This review provides a general introduction to the key points of parameter estimation with a focus on the PBPK and QSP models, and the respective parameter estimation algorithms. The latter part assesses the performance of five parameter estimation algorithms - the quasi-Newton method, Nelder-Mead method, genetic algorithm, particle swarm optimization, and Cluster Gauss-Newton method - using three examples of PBPK and QSP modeling. The assessment revealed that some parameter estimation results were significantly influenced by the initial values. Moreover, the choice of algorithms demonstrating good estimation results heavily depends on factors such as model structure and the parameters to be estimated. To obtain credible parameter estimation results, it is advisable to conduct multiple rounds of parameter estimation under different conditions, employing various estimation algorithms.

Assessing Pleiotropic Effects of a Mixed-Mode Perpetrator Drug, Rifampicin, by Multiple Endogenous Biomarkers in Dogs

Rifampicin (RIF) is a mixed-mode perpetrator that produces pleiotropic effects on liver cytochrome P450 enzymes and drug transporters. To assess the complex drug-drug interaction liabilities of RIF in vivo, a known probe substrate, midazolam (MDZ), along with multiple endogenous biomarkers were simultaneously monitored in beagle dogs before and after a 7-day treatment period by RIF at 20 mg/kg per day. Confirmed by the reduced MDZ plasma exposure and elevated 4β-hydroxycholesterol (4β-HC, biomarker of CYP3A activities) level, CYP3A was significantly induced after repeated RIF doses, and such induction persisted for 3 days after cessation of the RIF administration. On the other hand, increased plasma levels of coproporphyrin (CP)-I and III [biomarkers of organic anion transporting polypeptides 1b (Oatp1b) activities] were observed after the first dose of RIF. Plasma CPs started to decline as RIF exposure decreased, and they returned to baseline 3 days after cessation of the RIF administration. The data suggested the acute (inhibitory) and chronic (inductive) effects of RIF on Oatp1b and CYP3A enzymes, respectively, and a 3-day washout period is deemed adequate to remove superimposed Oatp1b inhibition from CYP3A induction. In addition, apparent self-induction of RIF was observed as its terminal half-life was significantly altered after multiple doses. Overall, our investigation illustrated the need for appropriate timing of modulator dosing to differentiate between transporter inhibition and enzyme induction. As further indicated by the CP data, induction of Oatp1b activities was not likely after repeated RIF administration. SIGNIFICANCE STATEMENT: This investigation demonstrated the utility of endogenous biomarkers towards complex drug-drug interactions by rifampicin (RIF) and successfully determined the optimal timing to differentiate between transporter inhibition and enzyme induction. Based on experimental evidence, Oatp1b induction following repeated RIF administration was unlikely, and apparent self-induction of RIF elimination was observed.

Drug-drug interaction prediction of ziritaxestat using a physiologically based enzyme and transporter pharmacokinetic network interaction model

Ziritaxestat, an autotaxin inhibitor, was under development for the treatment of idiopathic pulmonary fibrosis. It is a substrate of cytochrome P450 3A4 (CYP3A4) and P-glycoprotein and a weak inhibitor of the CYP3A4 and OATP1B1 pathways. We developed a physiologically based pharmacokinetic (PBPK) network interaction model for ziritaxestat that incorporated its metabolic and transporter pathways, enabling prediction of its potential as a victim or perpetrator of drug-drug interactions (DDIs). Concurrently, we evaluated CYP3A4 autoinhibition, including time-dependent inhibition. In vitro information and clinical data from healthy volunteer studies were used for model building and validation. DDIs with rifampin, itraconazole, voriconazole, pravastatin, and rosuvastatin were predicted, followed by validation against a test dataset. DDIs of ziritaxestat as a victim or perpetrator were simulated using the final model. Predicted-to-observed DDI ratios for the maximum plasma concentration (Cmax ) and the area under the plasma concentration-time curve (AUC) were within a two-fold ratio for both the metabolic and transporter-mediated simulated DDIs. The predicted impact of autoinhibition/autoinduction or time-dependent inhibition of CYP3A4 was a 12% decrease in exposure. Model-based predictions for ziritaxestat as a victim of DDIs with a moderate CYP3A4 inhibitor (fluconazole) suggested a 2.6-fold increase in the AUC of ziritaxestat, while multiple doses of a strong inhibitor (voriconazole) would increase the AUC by 15-fold. Efavirenz would yield a three-fold decrease in the AUC of ziritaxestat. As a perpetrator, ziritaxestat was predicted to increase the AUC of the CYP3A4 index substrate midazolam by 2.7-fold. An overarching PBPK model was developed that could predict DDI liability of ziritaxestat for both CYP3A4 and the transporter pathways.

Physiologically-based pharmacokinetic modeling for investigating the effect of simeprevir on concomitant drugs and an endogenous biomarker of OATP1B

The orally available anti-hepatitis C virus (HCV) drug simeprevir exhibits nonlinear pharmacokinetics at the clinical doses due to saturation of cytochrome P450 (CYP) 3A4 metabolism and organic anion transporting peptide (OATP) 1B mediated hepatic uptake. Additionally, simeprevir increases exposures of concomitant drugs by CYP3A4 and OATP1B inhibition. The objective of this study was to develop physiologically-based pharmacokinetic (PBPK) models that could describe drug-drug interactions (DDIs) of simeprevir with concomitant drugs via CYP3A4 and OATP1B inhibition, and also to capture the effects on coproporphyrin-I (CP-I), an endogenous biomarker of OATP1B. PBPK modeling estimated unbound simeprevir inhibitory constant (Ki ) of 2.89 μM against CYP3A4 in the DDI results between simeprevir and midazolam in healthy volunteers. Then, we analyzed the DDIs between simeprevir and atorvastatin, a dual substrate of CYP3A4 and OATP1B, in healthy volunteers, and unbound Ki against OATP1B was estimated to be 0.00347 μM. Finally, we analyzed the increase in the blood level of CP-I by simeprevir to verify the Ki,OATP1B . Because CP-I was measured in subjects with HCV with various hepatic fibrosis state, Monte Carlo simulation was performed to involve the decreases in expression levels of hepatic CYP3A4 and OATP1B and their interindividual variabilities. The PBPK modeling coupled with Monte Carlo simulation using the Ki,OATP1B value obtained from atorvastatin study reasonably recovered the observed relationship between CP-I and simeprevir blood levels. In conclusion, the simeprevir PBPK model developed in this study can quantitatively describe the increase in exposures of concomitant drugs and an endogenous biomarker via inhibition of CYP3A4 and OATP1B.

Physiologically-based pharmacokinetic modeling for optimal dosage prediction of olaparib when co-administered with CYP3A4 modulators and in patients with hepatic/renal impairment

This study aimed to develop a physiologically-based pharmacokinetic (PBPK) model to predict the maximum plasma concentration (Cmax) and trough concentration (Ctrough) at steady-state of olaparib (OLA) in Caucasian, Japanese and Chinese. Furthermore, the PBPK model was combined with mean and 95% confidence interval to predict optimal dosing regimens of OLA when co-administered with CYP3A4 modulators and administered to patients with hepatic/renal impairment. The dosing regimens were determined based on safety and efficacy PK threshold Cmax (< 12,500 ng/mL) and Ctrough (772-2500 ng/mL). The population PBPK model for OLA was successfully developed and validated, demonstrating good consistency with clinically observed data. The ratios of predicted to observed values for Cmax and Ctrough fell within the range of 0.5 to 2.0. When OLA was co-administered with a strong or moderate CYP3A4 inhibitor, the recommended dosing regimens should be reduced to 100 mg BID and 150 mg BID, respectively. Additionally, the PBPK model also suggested that OLA could be not recommended with a strong or moderate CYP3A4 inducer. For patients with moderate hepatic and renal impairment, the dosing regimens of OLA were recommended to be reduced to 200 mg BID and 150 mg BID, respectively. In cases of severe hepatic and renal impairment, the PBPK model suggested a dosing regimen of 100 mg BID for OLA. Overall, this present PBPK model can determine the optimal dosing regimens for various clinical scenarios involving OLA.

A 20-Year Research Overview: Quantitative Prediction of Hepatic Clearance Using the In Vitro-In Vivo Extrapolation Approach Based on Physiologically Based Pharmacokinetic Modeling and Extended Clearance Concept

Understanding the extended clearance concept and establishing a physiologically based pharmacokinetic (PBPK) model are crucial for investigating the impact of changes in transporter and metabolizing enzyme abundance/functions on drug pharmacokinetics in blood and tissues. This mini-review provides an overview of the extended clearance concept and a PBPK model that includes transporter-mediated uptake processes in the liver. In general, complete in vitro and in vivo extrapolation (IVIVE) poses challenges due to missing factors that bridge the gap between in vitro and in vivo systems. By considering key in vitro parameters, we can capture in vivo pharmacokinetics, a strategy known as the top-down or middle-out approach. We present the latest progress, theory, and practice of the Cluster Gauss-Newton method, which is used for middle-out analyses. As examples of poor IVIVE, we discuss "albumin-mediated hepatic uptake" and "time-dependent inhibition" of OATP1Bs. The hepatic uptake of highly plasma-bound drugs is more efficient than what can be accounted for by their unbound concentration alone. This phenomenon is referred to as "albumin-mediated" hepatic uptake. IVIVE was improved by measuring hepatic uptake clearance in vitro in the presence of physiologic albumin concentrations. Lastly, we demonstrate the application of Cluster Gauss-Newton method-based analysis to the target-mediated drug disposition of bosentan. Incorporating saturable target binding and OATP1B-mediated hepatic uptake into the PBPK model enables the consideration of nonlinear kinetics across a wide dose range and the prediction of receptor occupancy over time. SIGNIFICANCE STATEMENT: There have been multiple instances where researchers' endeavors to unravel the underlying mechanism of poor in vitro-in vivo extrapolation have led to the discovery of previously undisclosed truths. These include 1) albumin-mediated hepatic uptake, 2) the target-mediated drug disposition in small molecules, and 3) the existence of a trans-inhibition mechanism by inhibitors for OATP1B-mediated hepatic uptake of drugs. Consequently, poor in vitro-in vivo extrapolation and the subsequent inquisitiveness of scientists may serve as a pivotal gateway to uncover hidden mechanisms.

Pharmacokinetic and Pharmacodynamic Drug-Drug Interactions: Research Methods and Applications

Because of the high research and development cost of new drugs, the long development process of new drugs, and the high failure rate at later stages, combining past drugs has gradually become a more economical and attractive alternative. However, the ensuing problem of drug-drug interactions (DDIs) urgently need to be solved, and combination has attracted a lot of attention from pharmaceutical researchers. At present, DDI is often evaluated and investigated from two perspectives: pharmacodynamics and pharmacokinetics. However, in some special cases, DDI cannot be accurately evaluated from a single perspective. Therefore, this review describes and compares the current DDI evaluation methods based on two aspects: pharmacokinetic interaction and pharmacodynamic interaction. The methods summarized in this paper mainly include probe drug cocktail methods, liver microsome and hepatocyte models, static models, physiologically based pharmacokinetic models, machine learning models, in vivo comparative efficacy studies, and in vitro static and dynamic tests. This review aims to serve as a useful guide for interested researchers to promote more scientific accuracy and clinical practical use of DDI studies.

Can Mechanistic Static Models for Drug-Drug Interactions Support Regulatory Filing for Study Waivers and Label Recommendations?

BACKGROUND AND OBJECTIVE

Mechanistic static and dynamic physiologically based pharmacokinetic models are used in clinical drug development to assess the risk of drug-drug interactions (DDIs). Currently, the use of mechanistic static models is restricted to screening DDI risk for an investigational drug, while dynamic physiologically based pharmacokinetic models are used for quantitative predictions of DDIs to support regulatory filing. As physiologically based pharmacokinetic model development by sponsors as well as a review of models by regulators require considerable resources, we explored the possibility of using mechanistic static models to support regulatory filing, using representative cases of successful physiologically based pharmacokinetic submissions to the US Food and Drug Administration under different classes of applications.

METHODS

Drug-drug interaction predictions with mechanistic static models were done for representative cases in the different classes of applications using the same data and modelling workflow as described in the Food and Drug Administration clinical pharmacology reviews. We investigated the hypothesis that the use of unbound average steady-state concentrations of modulators as driver concentrations in the mechanistic static models should lead to the same conclusions as those from physiologically based pharmacokinetic modelling for non-dynamic measures of DDI risk assessment such as the area under the plasma concentration-time curve ratio, provided the same input data are employed for the interacting drugs.

RESULTS

Drug-drug interaction predictions of area under the plasma concentration-time curve ratios using mechanistic static models were mostly comparable to those reported in the Food and Drug Administration reviews using physiologically based pharmacokinetic models for all representative cases in the different classes of applications.

CONCLUSIONS

The results reported in this study should encourage the use of models that best fit an intended purpose, limiting the use of physiologically based pharmacokinetic models to those applications that leverage its unique strengths, such as what-if scenario testing to understand the effect of dose staggering, evaluating the role of uptake and efflux transporters, extrapolating DDI effects from studied to unstudied populations, or assessing the impact of DDIs on the exposure of a victim drug with concurrent mechanisms. With this first step, we hope to trigger a scientific discussion on the value of a routine comparison of the two methods for regulatory submissions to potentially create a best practice that could help identify examples where the use of dynamic changes in modulator concentrations could make a difference to DDI risk assessment.

Analysis of Complex Absorption After Multiple Dosing: Application to the Interaction Between the P-glycoprotein Substrate Talinolol and Rifampicin

PURPOSE

In order to clarify the effect of rifampicin on the bioavailability of the P-glycoprotein substrate talinolol, its absorption kinetics was modeled after multiple-dose oral administration of talinolol in healthy subjects.

METHODS

A sum of two inverse Gaussian functions was used to calculate the time course of the input rate into the systemic circulation.

RESULTS

The estimated rate of drug entry into the systemic circulation revealed two distinct peaks at 1 and 3.5 h after administration. Rifampicin did not affect bioavailability of talinolol, but did shift the second peak of the input function by 1.3 h to later times. Elimination clearance and one of the intercompartmental distribution clearances increased significantly under rifampicin treatment.

CONCLUSIONS

Rifampicin changes the time course of absorption rate but not the fraction absorbed of talinolol. The model suggests the existence of two intestinal absorption windows for talinolol.

Alternatives to rifampicin: A review and perspectives on the choice of strong CYP3A inducers for clinical drug-drug interaction studies

N-Nitrosamine (NA) impurities are considered genotoxic and have gained attention due to the recall of several marketed drug products associated with higher-than-permitted limits of these impurities. Rifampicin is an index inducer of multiple cytochrome P450s (CYPs) including CYP2B6, 2C8, 2C9, 2C19, and 3A4/5 and an inhibitor of OATP1B transporters (single dose). Hence, rifampicin is used extensively in clinical studies to assess drug-drug interactions (DDIs). Despite NA impurities being reported in rifampicin and rifapentine above the acceptable limits, these critical anti-infective drugs are available for therapeutic use considering their benefit-risk profile. Reports of NA impurities in rifampicin products have created uncertainty around using rifampicin in clinical DDI studies, especially in healthy volunteers. Hence, a systematic investigation through a literature search was performed to determine possible alternative index inducer(s) to rifampicin. The available strong CYP3A inducers were selected from the University of Washington DDI Database and their in vivo DDI potential assessed using the data from clinical DDI studies with sensitive CYP3A substrates. To propose potential alternative CYP3A inducers, factors including lack of genotoxic potential, adequate safety, feasibility of multiple dose administration to healthy volunteers, and robust in vivo evidence of induction of CYP3A were considered. Based on the qualifying criteria, carbamazepine, phenytoin, and lumacaftor were identified to be the most promising alternatives to rifampicin for conducting CYP3A induction DDI studies. Strengths and limitations of the proposed alternative CYP3A inducers, the magnitude of in vivo CYP3A induction, appropriate study designs for each alternative inducer, and future perspectives are presented in this paper.

The Effect of Rifampicin on Darunavir, Ritonavir, and Dolutegravir Exposure within Peripheral Blood Mononuclear Cells: a Dose Escalation Study

Ritonavir-boosted darunavir (DRV/r) and dolutegravir (DTG) are affected by induction of metabolizing enzymes and efflux transporters caused by rifampicin (RIF). This complicates the treatment of people living with HIV (PLWH) diagnosed with tuberculosis. Recent data showed that doubling DRV/r dose did not compensate for this effect, and hepatic safety was unsatisfactory. We aimed to evaluate the pharmacokinetics of DRV, ritonavir (RTV), and DTG in the presence and absence of RIF in peripheral blood mononuclear cells (PBMCs). PLWH were enrolled in a dose-escalation crossover study with 6 treatment periods of 7 days. Participants started with DRV/r 800/100 mg once daily (QD), RIF and DTG were added before the RTV dose was doubled, and then they received DRV/r 800/100 twice daily (BD) and then 1,600/200 QD or vice versa. Finally, RIF was withdrawn. Plasma and intra-PBMC drug concentrations were measured through validated liquid chromatography-tandem mass spectrometry (LC-MS/MS) methods. Seventeen participants were enrolled but only 4 completed all study phases due to high incidence of liver toxicity. Intra-PBMC DRV trough serum concentration (C_trough) after the addition of RIF dropped from a median (interquartile range [IQR]) starting value of 261 ng/mL (158 to 577) to 112 ng/mL (18 to 820) and 31 ng/mL (12 to 331) for 800/100 BD and 1,600/200 QD DRV/r doses, respectively. The DRV intra-PBMC/plasma ratio increased significantly (P = 0.003). DTG and RIF intra-PBMC concentrations were in accordance with previous reports in the absence of RIF or DRV/r. This study showed a differential impact of enzyme and/or transporter induction on DRV/r concentrations in plasma and PBMCs, highlighting the usefulness of studying intra-PBMC pharmacokinetics with drug-drug interactions. (This study has been registered at ClinicalTrials.gov under registration no. NCT03892161.)

Understanding Inter-individual Variability in the Drug Interaction of a Highly Extracted CYP1A2 Substrate Tizanidine: Application of a Permeability-limited Multi-compartment Liver Model in a Population Based PBPK Framework

Tizanidine, a centrally acting skeletal muscle relaxant, is predominantly metabolised by CYP1A2 and undergoes extensive hepatic first-pass metabolism following oral administration. As a highly extracted drug, the systemic exposure to tizanidine exhibits considerable inter-individual variability and is altered substantially when co-administered with CYP1A2 inhibitors or inducers. The aim of the current study was to compare the performance of a permeability-limited multi-compartment liver (PerMCL) model, which operates as an approximation of the dispersion model (DM), and the well-stirred model (WSM) for predicting tizanidine DDIs. Physiologically-based pharmacokinetic (PBPK) models were developed for tizanidine, incorporating the PerMCL model and the WSM, respectively, to simulate the interaction of tizanidine with a range of CYP1A2 inhibitors and inducers. While the WSM showed a tendency to under-predict the fold change of tizanidine AUC (AUC ratio) in the presence of perpetrators, the use of PerMCL model increased precision (absolute average-fold error: 1.32 - 1.42 versus 1.58) and decreased bias (average-fold error: 0.97 - 1.25 versus 0.63) for the predictions of mean AUC ratios as compared to the WSM. The PerMCL model captured the observed range of individual AUC ratios of tizanidine as well as the correlation between individual AUC ratios and CYP1A2 activities without interactions, whereas the WSM was not able to capture these. The results demonstrate the advantage of using the PerMCL model over the WSM in predicting the magnitude and inter-individual variability of DDIs for a highly extracted sensitive substrate tizanidine. Significance Statement This study demonstrates the advantages of the permeability-limited multi-compartment liver (PerMCL) model, which operates as an approximation of the dispersion model (DM), in mitigating the tendency of the well-stirred model (WSM) to under-predict the magnitude and variability of DDIs of a highly extracted CYP1A2 substrate tizanidine when it is administered with CYP1A2 inhibitors or inducers. The PBPK modelling approach described herein is valuable to the understanding of drug interactions of highly extracted substrates and the source of its inter-individual variability.

Effect of CYP3A4 Inhibitors and Inducers on Pharmacokinetics and Pharmacodynamics of Saxagliptin and Active Metabolite M2 in Humans Using Physiological-Based Pharmacokinetic Combined DPP-4 Occupancy

We aimed to develop a physiological-based pharmacokinetic and dipepidyl peptidase 4 (DPP-4) occupancy model (PBPK-DO) characterized by two simultaneous simulations to predict pharmacokinetic (PK) and pharmacodynamic changes of saxagliptin and metabolite M2 in humans when coadministered with CYP3A4 inhibitors or inducers. Ketoconazole, delavirdine, and rifampicin were selected as a CYP3A4 competitive inhibitor, a time-dependent inhibitor, and an inducer, respectively. Here, we have successfully simulated PK profiles and DPP-4 occupancy profiles of saxagliptin in humans using the PBPK-DO model. Additionally, under the circumstance of actually measured values, predicted results were good and in line with observations, and all fold errors were below 2. The prediction results demonstrated that the oral dose of saxagliptin should be reduced to 2.5 mg when coadministrated with ketoconazole. The predictions also showed that although PK profiles of saxagliptin showed significant changes with delavirdine (AUC 1.5-fold increase) or rifampicin (AUC: a decrease to 0.19-fold) compared to those without inhibitors or inducers, occupancies of DPP-4 by saxagliptin were nearly unchanged, that is, the administration dose of saxagliptin need not adjust when there is coadministration with delavirdine or rifampicin.

Physiologically-Based Pharmacokinetic Modelling of Entrectinib Parent and Active Metabolite to Support Regulatory Decision-Making

BACKGROUND AND OBJECTIVE

Entrectinib is a selective inhibitor of ROS1/TRK/ALK kinases, recently approved for oncology indications. Entrectinib is predominantly cleared by cytochrome P450 (CYP) 3A4, and modulation of CYP3A enzyme activity profoundly alters the pharmacokinetics of both entrectinib and its active metabolite M5. We describe development of a combined physiologically based pharmacokinetic (PBPK) model for entrectinib and M5 to support dosing recommendations when entrectinib is co-administered with CYP3A4 inhibitors or inducers.

METHODS

A PBPK model was established in Simcyp^® Simulator. The initial model based on in vitro-in vivo extrapolation was refined using sensitivity analysis and non-linear mixed effects modeling to optimize parameter estimates and to improve model fit to data from a clinical drug-drug interaction study with the strong CYP3A4 inhibitor, itraconazole. The model was subsequently qualified against clinical data, and the final qualified model used to simulate the effects of moderate to strong CYP3A4 inhibitors and inducers on entrectinib and M5 pharmacokinetics.

RESULTS

The final model showed good predictive performance for entrectinib and M5, meeting commonly used predictive performance acceptance criteria in each case. The model predicted that co-administration of various moderate CYP3A4 inhibitors (verapamil, erythromycin, clarithromycin, fluconazole, and diltiazem) would result in an average increase in entrectinib exposure between 2.2- and 3.1-fold, with corresponding average increases for M5 of approximately 2-fold. Co-administration of moderate CYP3A4 inducers (efavirenz, carbamazepine, phenytoin) was predicted to result in an average decrease in entrectinib exposure between 45 and 79%, with corresponding average decreases for M5 of approximately 50%.

CONCLUSIONS

The model simulations were used to derive dosing recommendations for co-administering entrectinib with CYP3A4 inhibitors or inducers. PBPK modeling has been used in lieu of clinical studies to enable regulatory decision-making.

Population pharmacokinetic modeling and simulation to support qualification of pyridoxic acid as endogenous biomarker of OAT1/3 renal transporters

Renal clearance of many drugs is mediated by renal organic anion transporters OAT1/3 and inhibition of these transporters may lead to drug‐drug interactions (DDIs). Pyridoxic acid (PDA) and homovanillic acid (HVA) were indicated as potential biomarkers of OAT1/3. The objective of this study was to develop a population pharmacokinetic model for PDA and HVA to support biomarker qualification. Simultaneous fitting of biomarker plasma and urine data in the presence and absence of potent OAT1/3 inhibitor (probenecid, 500 mg every 6 h) was performed. The impact of study design (multiple vs. single dose of OAT1/3 inhibitor) and ability to detect interactions in the presence of weak/moderate OAT1/3 inhibitors was investigated, together with corresponding power calculations. The population models developed successfully described biomarker baseline and PDA/HVA OAT1/3‐mediated interaction data. No prominent effect of circadian rhythm on PDA and HVA individual baseline levels was evident. Renal elimination contributed greater than 80% to total clearance of both endogenous biomarkers investigated. Estimated probenecid unbound in vivo OAT inhibitory constant was up to 6.4‐fold lower than in vitro values obtained with PDA as a probe. The PDA model was successfully verified against independent literature reported datasets. No significant difference in power of DDI detection was found between multiple and single dose study design when using the same total daily dose of 2000 mg probenecid. Model‐based simulations and power calculations confirmed sensitivity and robustness of plasma PDA data to identify weak, moderate, and strong OAT1/3 inhibitors in an adequately powered clinical study to support optimal design of prospective clinical OAT1/3 interaction studies.

Inhibition and induction of CYP enzymes in humans: an update

The cytochrome P450 (CYP) enzyme family is the most important enzyme system catalyzing the phase 1 metabolism of pharmaceuticals and other xenobiotics such as herbal remedies and toxic compounds in the environment. The inhibition and induction of CYPs are major mechanisms causing pharmacokinetic drug–drug interactions. This review presents a comprehensive update on the inhibitors and inducers of the specific CYP enzymes in humans. The focus is on the more recent human in vitro and in vivo findings since the publication of our previous review on this topic in 2008. In addition to the general presentation of inhibitory drugs and inducers of human CYP enzymes by drugs, herbal remedies, and toxic compounds, an in-depth view on tyrosine-kinase inhibitors and antiretroviral HIV medications as victims and perpetrators of drug–drug interactions is provided as examples of the current trends in the field. Also, a concise overview of the mechanisms of CYP induction is presented to aid the understanding of the induction phenomena.

Improved Prediction of the Drug-Drug Interactions of Pemafibrate Caused by Cyclosporine A and Rifampicin via PBPK Modeling: Consideration of the Albumin-Mediated Hepatic Uptake of Pemafibrate and Inhibition Constants With Preincubation Against OATP1B

Pemafibrate (PMF) is highly albumin-bound (>99.8%) and a substrate for hepatic uptake transporters (OATP1B) and CYP enzymes. Here, we developed a PBPK model of PMF to capture drug-drug interactions (DDI) incurred by cyclosporine (CsA) and rifampicin (RIF), the two OATP1B inhibitors. Initial PBPK modeling of PMF utilized in vitro hepatic uptake clearance (PS_inf) obtained in the absence of albumin, but failed in capturing the blood PMF pharmacokinetic (PK) profiles. Based on the results that in vitro PS_inf of unbound PMF was enhanced in the presence of albumin, we applied the albumin-facilitated dissociation model and the resulting PS_inf parameters improved the prediction of the blood PMF PK profiles. In refining our PBPK model toward improved prediction of the observed DDI data (PMF co-administered with single dosing of CsA or RIF; PMF following multiple RIF dosing), we adjusted the previously obtained in vivo OATP1B inhibition constants (K_i,OATP1B) of CsA or RIF for pitavastatin by correcting for substrate-dependency. We also incorporated the induction of OATP1B and CYP enzymes after multiple RIF dosing. Sensitivity analysis informed that the higher gastrointestinal absorption rate constant could further improve capturing the observed DDI data, suggesting the possible inhibition of intestinal ABC transporter(s) by CsA or RIF.

Application of a microdosed cocktail of 3 oral factor Xa inhibitors to study drug-drug interactions with different perpetrator drugs

AIMS

Using 3 different perpetrators the impact of voriconazole, cobicistat and rifampicin (single dose), we evaluated the suitability of a microdose cocktail of factor Xa inhibitors (FXaI; rivaroxaban, apixaban and edoxaban; 100 μg in total) to study drug-drug interactions.

METHODS

Three cohorts of 6 healthy volunteers received 2 treatments with microdoses of rivaroxaban, apixaban and edoxaban alone and with coadministration of 1 of the perpetrators. Plasma and urine concentrations of microdosed apixaban, edoxaban and rivaroxaban were quantified using a validated ultra-performance liquid chromatography-tandem mass spectrometry with a lower limit of quantification of 2.5 pg/mL.

RESULTS

Voriconazole caused only a minor interaction with apixaban and rivaroxaban, none with edoxaban. Cobicistat significantly increased exposure of all 3 FXaI with area under the plasma concentration-time curve ratios of 1.67 (apixaban), 1.74 (edoxaban) and 2.0 (rivaroxaban). A single dose of rifampicin decreased the volume of distribution and elimination half-life of all 3 FXaI.

CONCLUSIONS

The microdosed FXaI cocktail approach is able to generate drug interaction data and can help elucidating the mechanism involved in the clearance of the different victim drugs. This is a safe approach to concurrently study drug-drug interactions with a drug class. (EudraCT 2016-003024-23).

Quantitative Assessment of Elagolix Enzyme-Transporter Interplay and Drug-Drug Interactions Using Physiologically Based Pharmacokinetic Modeling

INTRODUCTION

Elagolix is approved for the management of moderate-to-severe pain associated with endometriosis. The aim of this analysis was to develop a physiologically based pharmacokinetic (PBPK) model that describes the enzyme-transporter interplay involved in the disposition of elagolix and to predict the magnitude of drug-drug interaction (DDI) potential of elagolix as an inhibitor of P-glycoprotein (P-gp) and inducer of cytochrome P450 (CYP) 3A4.

METHODS

A PBPK model (SimCYP® version 15.0.86.0) was developed using elagolix data from in vitro, clinical PK and DDI studies. Data from DDI studies were used to quantify contributions of the uptake transporter organic anion transporting polypeptide (OATP) 1B1 and CYP3A4 in the disposition of elagolix, and to quantitatively assess the perpetrator potential of elagolix as a CYP3A4 inducer and P-gp inhibitor.

RESULTS

After accounting for the interplay between elagolix metabolism by CYP3A4 and uptake by OATP1B1, the model-predicted PK parameters of elagolix along with the DDI AUC∞ and Cmax ratios, were within 1.5-fold of the observed data. Based on model simulations, elagolix 200 mg administered twice daily is a moderate inducer of CYP3A4 (approximately 56% reduction in midazolam AUC∞). Simulations of elagolix 150 mg administered once daily with digoxin predicted an increase in digoxin Cmax and AUC∞ by 68% and 19%, respectively.

CONCLUSIONS

A PBPK model of elagolix was developed, verified, and applied to characterize the disposition interplay between CYP3A4 and OATP1B1, and to predict the DDI potential of elagolix as a perpetrator under dosing conditions that were not tested clinically. PBPK model-based predictions were used to support labeling language for DDI recommendations of elagolix.

Implementation of a Physiologically Based Pharmacokinetic Modeling Approach to Guide Optimal Dosing Regimens for Imatinib and Potential Drug Interactions in Paediatrics

Long-term use of imatinib is effective and well-tolerated in children with chronic myeloid leukaemia (CML) yet defining an optimal dosing regimen for imatinib in younger patients is a challenge. The potential interactions between imatinib and coadministered drugs in this “special” population also remains largely unexplored. This study implements a physiologically based pharmacokinetic (PBPK) modeling approach to investigate optimal dosing regimens and potential drug interactions with imatinib in the paediatric population. A PBPK model for imatinib was developed in the Simcyp Simulator (version 17) utilizing in silico, in vitro drug metabolism, and in vivo pharmacokinetic data and verified using an independent set of published clinical pharmacokinetic data. The model was then extrapolated to children and adolescents (aged 2–18 years) by incorporating developmental changes in organ size and maturation of drug-metabolising enzymes and plasma protein responsible for imatinib disposition. The PBPK model described imatinib pharmacokinetics in adult and paediatric populations and predicted drug interaction with carbamazepine, a cytochrome P450 (CYP)3A4 and 2C8 inducer, with a good accuracy (evaluated by visual inspections of the simulation results and predicted pharmacokinetic parameters that were within 1.25-fold of the clinically observed values). The PBPK simulation suggests that the optimal dosing regimen range for imatinib is 230–340 mg/m²/d in paediatrics, which is supported by the recommended initial dose for treatment of childhood CML. The simulations also highlighted that children and adults being treated with imatinib have similar vulnerability to CYP modulations. A PBPK model for imatinib was successfully developed with an excellent performance in predicting imatinib pharmacokinetics across age groups. This PBPK model is beneficial to guide optimal dosing regimens for imatinib and predict drug interactions with CYP modulators in the paediatric population.

Recent progress in in vivo phenotyping technologies for better prediction of transporter-mediated drug-drug interactions

Clinical reports on transporter-mediated drug-drug interactions (TP-DDIs) have rapidly accumulated and regulatory guidance/guidelines recommend that sponsors consider performing quantitative prediction of TP-DDI risks in the process of drug development. In vitro experiments for characterizing the function of drug transporters have been established and various parameters such as the inhibition constant (K_i) of drugs and the intrinsic uptake/efflux clearance for a certain transporter can be obtained. However, many reports have indicated large discrepancies between the parameters estimated from in vitro experiments and those rationally explaining drug pharmacokinetics. Thus, it is essential to evaluate directly the function of each transporter isoform in vivo in humans. At present, several transporter substrate drugs and endogenous compounds have been recognized as probe substrates for a specific transporter and transporter function was evaluated by monitoring the plasma and urine concentration of those probes; however, few compounds specifically transported via a single transporter isoform have been found. For monitoring the intraorgan concentration of drugs, positron emission tomography can be a powerful tool and clinical examples for quantification of in vivo transporter function have been published. In this review, novel methodologies for in vivo phenotyping of transporter function are summarized.

Quantitative evaluation of hepatic and intestinal induction of CYP3A in clinical practice

This is the first report quantitatively evaluating the clinical induction of CYP3A in the liver and the intestine.To evaluate hepatic induction, we collected literature data on endogenous biomarkers of hepatic CYP3A induction which we then used to calculate the fold-induction (inducer-mediated change in biomarker level). Literature data on decreases in the area under the curve (AUC) of alfentanil, a CYP3A substrate, caused by CYP3A inducers were also collected. We used the hepatic intrinsic clearance of alfentanil to calculate the hepatic induction ratio (inducer-mediated change in intrinsic clearance). For intestinal induction, the intestinal bioavailability (F_g) of alfentanil was used to calculate the intestinal induction ratio. We determined in vivo maximum induction (E_max) and the average unbound plasma concentration (C_av,u) required for half the maximum induction (EC₅₀) for inducers using an E_max model analysis.In our results, fold-induction was comparable to the induction ratio at several inducer concentrations, and almost the maximum induction was achieved by a therapeutic dose. Induction ratios in the intestine were similar to the liver.Our findings suggest that, by knowing only hepatic induction ratios for common inducers, we can quantitatively predict the decreases in the AUC of substrates by CYP3A induction.

Physiologically-Based Pharmacokinetic Models for Evaluating Membrane Transporter Mediated Drug-Drug Interactions: Current Capabilities, Case Studies, Future Opportunities, and Recommendations

Physiologically-based pharmacokinetic (PBPK) modeling has been extensively used to quantitatively translate in vitro data and evaluate temporal effects from drug-drug interactions (DDIs), arising due to reversible enzyme and transporter inhibition, irreversible time-dependent inhibition, enzyme induction, and/or suppression. PBPK modeling has now gained reasonable acceptance with the regulatory authorities for the cytochrome-P450-mediated DDIs and is routinely used. However, the application of PBPK for transporter-mediated DDIs (tDDI) in drug development is relatively uncommon. Because the predictive performance of PBPK models for tDDI is not well established, here, we represent and discuss examples of PBPK analyses included in regulatory submission (the US Food and Drug Administration (FDA), the European Medicines Agency (EMA), and the Pharmaceuticals and Medical Devices Agency (PMDA)) across various tDDIs. The goal of this collaborative effort (involving scientists representing 17 pharmaceutical companies in the Consortium and from academia) is to reflect on the use of current databases and models to address tDDIs. This challenges the common perceptions on applications of PBPK for tDDIs and further delves into the requirements to improve such PBPK predictions. This review provides a reflection on the current trends in PBPK modeling for tDDIs and provides a framework to promote continuous use, verification, and improvement in industrialization of the transporter PBPK modeling.

Development and Evaluation of Physiologically Based Pharmacokinetic Drug-Disease Models for Predicting Rifampicin Exposure in Tuberculosis and Cirrhosis Populations

The physiologically based pharmacokinetic (PBPK) approach facilitates the construction of novel drug–disease models by allowing incorporation of relevant pathophysiological changes. The aim of the present work was to explore and identify the differences in rifampicin pharmacokinetics (PK) after the application of its single dose in healthy and diseased populations by using PBPK drug–disease models. The Simcyp^® simulator was used as a platform for modeling and simulation. The model development process was initiated by predicting rifampicin PK in healthy population after intravenous (i.v) and oral administration. Subsequent to successful evaluation in healthy population, the pathophysiological changes in tuberculosis and cirrhosis population were incorporated into the developed model for predicting rifampicin PK in these populations. The model evaluation was performed by using visual predictive checks and the comparison of mean observed/predicted ratios (ratio_(Obs/pred)) of the PK parameters. The predicted PK parameters in the healthy population were in adequate harmony with the reported clinical data. The incorporation of pathophysiological changes in albumin concentration in the tuberculosis population revealed improved prediction of clearance. The developed PBPK drug–disease models have efficiently described rifampicin PK in tuberculosis and cirrhosis populations after administering single drug dose, as the ratio_(Obs/pred) for all the PK parameters were within a two-fold error range. The mechanistic nature of the developed PBPK models may facilitate their extension to other diseases and drugs.

Expanded Physiologically-Based Pharmacokinetic Model of Rifampicin for Predicting Interactions With Drugs and an Endogenous Biomarker via Complex Mechanisms Including Organic Anion Transporting Polypeptide 1B Induction

As rifampicin can cause the induction and inhibition of multiple metabolizing enzymes and transporters, it has been challenging to accurately predict the complex drug–drug interactions (DDIs). We previously constructed a physiologically‐based pharmacokinetic (PBPK) model of rifampicin accounting for the components for the induction of cytochrome P450 (CYP) 3A/CYP2C9 and the inhibition of organic anion transporting polypeptide 1B (OATP1B). This study aimed to expand and verify the PBPK model for rifampicin by incorporating additional components for the induction of OATP1B and CYP2C8 and the inhibition of multidrug resistance protein 2. The established PBPK model was capable of accurately predicting complex rifampicin‐induced alterations in the profiles of glibenclamide, repaglinide, and coproporphyrin I (an endogenous biomarker of OATP1B activities) with various dosing regimens. Our comprehensive rifampicin PBPK model may enable quantitative prediction of DDIs across diverse potential victim drugs and endogenous biomarkers handled by multiple metabolizing enzymes and transporters.

Physiologically-Based Pharmacokinetic Modeling Approach to Predict Rifampin-Mediated Intestinal P-Glycoprotein Induction

Physiologically‐based pharmacokinetic (PBPK) modeling is a powerful tool to quantitatively describe drug disposition profiles in vivo, thereby providing an alternative to predict drug–drug interactions (DDIs) that have not been tested clinically. This study aimed to predict effects of rifampin‐mediated intestinal P‐glycoprotein (Pgp) induction on pharmacokinetics of Pgp substrates via PBPK modeling. First, we selected four Pgp substrates (digoxin, talinolol, quinidine, and dabigatran etexilate) to derive in vitro to in vivo scaling factors for intestinal Pgp kinetics. Assuming unbound Michaelis‐Menten constant (K_m) to be intrinsic, we focused on the scaling factors for maximal efflux rate (J_max) to adequately recover clinically observed results. Next, we predicted rifampin‐mediated fold increases in intestinal Pgp abundances to reasonably recover clinically observed DDI results. The modeling results suggested that threefold to fourfold increases in intestinal Pgp abundances could sufficiently reproduce the DDI results of these Pgp substrates with rifampin. Hence, the obtained fold increases can potentially be applicable to DDI prediction with other Pgp substrates.

Perspectives from the Innovation and Quality Consortium Induction Working Group on Factors Impacting Clinical Drug-Drug Interactions Resulting from Induction: Focus on Cytochrome 3A Substrates

A recent publication from the Innovation and Quality Consortium Induction Working Group collated a large clinical data set with the goal of evaluating the accuracy of drug-drug interaction (DDI) prediction from in vitro data. Somewhat surprisingly, comparison across studies of the mean- or median-reported area under the curve ratio showed appreciable variability in the magnitude of outcome. This commentary explores the possible drivers of this range of outcomes observed in clinical induction studies. While recommendations on clinical study design are not being proposed, some key observations were informative during the aggregate analysis of clinical data. Although DDI data are often presented using median data, individual data would enable evaluation of how differences in study design, baseline expression, and the number of subjects contribute. Since variability in perpetrator pharmacokinetics (PK) could impact the overall DDI interpretation, should this be routinely captured? Maximal induction was typically observed after 5-7 days of dosing. Thus, when the half-life of the inducer is less than 30 hours, are there benefits to a more standardized study design? A large proportion of CYP3A4 inducers were also CYP3A4 inhibitors and/or inactivators based on in vitro data. In these cases, using CYP3A selective substrates has limitations. More intensive monitoring of changes in area under the curve over time is warranted. With selective CYP3A substrates, the net effect was often inhibition, whereas less selective substrates could discern induction through mechanisms not susceptible to inhibition. The latter included oral contraceptives, which raise concerns of reduced efficacy following induction. Alternative approaches for modeling induction, such as applying biomarkers and physiologically based pharmacokinetic modeling (PBPK), are also considered. SIGNIFICANCE STATEMENT: The goal of this commentary is to stimulate discussion on whether there are opportunities to optimize clinical drug-drug interaction study design. The overall aim is to reduce, understand and contextualize the variability observed in the magnitude of induction across reported clinical studies. A large clinical CYP3A induction dataset was collected and further analyzed to identify trends and gaps. Reporting individual victim PK data, characterizing perpetrator PK and including additional PK assessments for mixed-mechanism perpetrators may provide insights into how these factors impact differences observed in clinical outcomes. The potential utility of biomarkers and PBPK modeling are discussed in considering future directions.

From Endogenous Compounds as Biomarkers to Plasma-Derived Nanovesicles as Liquid Biopsy; Has the Golden Age of Translational Pharmacokinetics-Absorption, Distribution, Metabolism, Excretion-Drug-Drug Interaction Science Finally Arrived?

It is now established that a drug's pharmacokinetics (PK) absorption, distribution, metabolism, excretion (ADME) and drug-drug interaction (DDI) profile can be modulated by age, disease, and genotype. In order to facilitate subject phenotyping and clinical DDI assessment, therefore, various endogenous compounds (in plasma and urine) have been pursued as drug-metabolizing enzyme and transporter biomarkers. Compared with biomarkers, however, the topic of circulating extracellular vesicles as "liquid biopsy" has received little attention within the ADME community; most organs secrete nanovesicles (e.g., exosomes) into the blood that contain luminal "cargo" derived from the originating organ (proteins, messenger RNA, and microRNA). As such, ADME profiling of plasma exosomes could be leveraged to better define genotype-phenotype relationships and the study of ontogeny, disease, and complex DDIs. If methods to support the isolation of tissue-derived plasma exosomes are successfully developed and validated, it is envisioned that they will be used jointly with genotyping, biomarkers, and modeling tools to greatly progress translational PK-ADME-DDI science.

PBPK Modeling of Coproporphyrin I as an Endogenous Biomarker for Drug Interactions Involving Inhibition of Hepatic OATP1B1 and OATP1B3

The aim of the present study was to establish a physiologically based pharmacokinetic (PBPK) model for coproporphyrin I (CP-I), a biomarker supporting the prediction of drug-drug interactions (DDIs) involving hepatic organic anion transporting polypeptide 1B (OATP1B), using clinical DDI data with an OATP1B inhibitor rifampicin (300 and 600 mg, orally). The in vivo inhibition constants of rifampicin used as initial input parameters for OATP1Bs (K_i,u,OATP1Bs ) and multidrug resistance-associated protein two-mediated biliary excretion were estimated as 0.23 and 0.87 μM, respectively, from previous reports. Sensitivity analysis demonstrated that the K_i,u,OATP1Bs and biosynthesis rate of CP-I affected the magnitude of the interaction. K_i,u,OATP1Bs values optimized by nonlinear least-squares fitting were ~0.5-fold of the initial value. It was determined that the blood concentration-time profiles of four statins were well-predicted using corrected individual K_i,u,OATP1B values (ratio of in vitro K_i,u(statin) /in vitro K_i,u(CP-I) ). In conclusion, PBPK modeling of CP-I supports dynamic prediction of OATP1B-mediated DDIs.

PBPK Models for CYP3A4 and P-gp DDI Prediction: A Modeling Network of Rifampicin, Itraconazole, Clarithromycin, Midazolam, Alfentanil, and Digoxin

According to current US Food and Drug Administration (FDA) and European Medicines Agency (EMA) guidance documents, physiologically based pharmacokinetic (PBPK) modeling is a powerful tool to explore and quantitatively predict drug‐drug interactions (DDIs) and may offer an alternative to dedicated clinical trials. This study provides whole‐body PBPK models of rifampicin, itraconazole, clarithromycin, midazolam, alfentanil, and digoxin within the Open Systems Pharmacology (OSP) Suite. All models were built independently, coupled using reported interaction parameters, and mutually evaluated to verify their predictive performance by simulating published clinical DDI studies. In total, 112 studies were used for model development and 57 studies for DDI prediction. 93% of the predicted area under the plasma concentration‐time curve (AUC) ratios and 94% of the peak plasma concentration (C_max) ratios are within twofold of the observed values. This study lays a cornerstone for the qualification of the OSP platform with regard to reliable PBPK predictions of enzyme‐mediated and transporter‐mediated DDIs during model‐informed drug development. All presented models are provided open‐source and transparently documented.

Assessment of Pharmacokinetic Drug-Drug Interactions in Humans: In Vivo Probe Substrates for Drug Metabolism and Drug Transport Revisited

Pharmacokinetic parameters of selective probe substrates are used to quantify the activity of an individual pharmacokinetic process (PKP) and the effect of perpetrator drugs thereon in clinical drug-drug interaction (DDI) studies. For instance, oral caffeine is used to quantify hepatic CYP1A2 activity, and oral dagibatran etexilate for intestinal P-glycoprotein (P-gp) activity. However, no probe substrate depends exclusively on the PKP it is meant to quantify. Lack of selectivity for a given enzyme/transporter and expression of the respective enzyme/transporter at several sites in the human body are the main challenges. Thus, a detailed understanding of the role of individual PKPs for the pharmacokinetics of any probe substrate is essential to allocate the effect of a perpetrator drug to a specific PKP; this is a prerequisite for reliably informed pharmacokinetic models that will allow for the quantitative prediction of perpetrator effects on therapeutic drugs, also in respective patient populations not included in DDI studies.