Comparison of different approaches to predict metabolic drug–drug interactions

Prediction of herb-drug interactions involving consumption of furanocoumarin-mixtures and cytochrome P450 1A2-mediated caffeine metabolism inhibition in humans

Herb-drug interactions (HDI) has become important due to the increasing popularity of natural health product consumption worldwide. HDI is difficult to predict as botanical drugs usually contain complex phytochemical-mixtures, which interact with drug metabolism. Currently, there is no specific pharmacological tool to predict HDI since almost all in vitro-in vivo-extrapolation (IVIVE) Drug-Drug Interaction (DDI) models deal with one inhibitor-drug and one victim-drug. The objectives were to modify-two IVIVE models for the prediction of in vivo interaction between caffeine and furanocoumarin-containing herbs, and to confirm model predictions by comparing the DDI predictive results with actual human data. The models were modified to predict in vivo herb-caffeine interaction using the same set of inhibition constants but different integrated dose/concentration of furanocoumarin mixtures in the liver. Different hepatic inlet inhibitor concentration ([I]_H) surrogates were used for each furanocoumarin. In the first (hybrid) model, the [I]_H was predicted using the concentration-addition model for chemical-mixtures. In the second model, the [I]_H was calculated by adding individual furanocoumarins together. Once [I]_H values were determined, the models predicted an area-under-curve-ratio (AUCR) value of each interaction. The results indicate that both models were able to predict the experimental AUCR of herbal products reasonably well. The DDI model approaches described in this study may be applicable to health supplements and functional foods also.

Beyond the Michaelis-Menten: Accurate Prediction of Drug Interactions through Cytochrome P450 3A4 Induction

The US Food and Drug Administration (FDA) guidance has recommended several model-based predictions to determine potential drug-drug interactions (DDIs) mediated by cytochrome P450 (CYP) induction. In particular, the ratio of substrate area under the plasma concentration-time curve (AUCR) under and not under the effect of inducers is predicted by the Michaelis-Menten (MM) model, where the MM constant ( K m $$ {K}_{\mathrm{m}} $$ ) of a drug is implicitly assumed to be sufficiently higher than the concentration of CYP enzymes that metabolize the drug ( E T $$ {E}_{\mathrm{T}} $$ ) in both the liver and small intestine. Furthermore, the fraction absorbed from gut lumen ( F a $$ {F}_{\mathrm{a}} $$ ) is also assumed to be one because F a $$ {F}_{\mathrm{a}} $$ is usually unknown. Here, we found that such assumptions lead to serious errors in predictions of AUCR. To resolve this, we propose a new framework to predict AUCR. Specifically, F a $$ {F}_{\mathrm{a}} $$ was re-estimated from experimental permeability values rather than assuming it to be one. Importantly, we used the total quasi-steady-state approximation to derive a new equation, which is valid regardless of the relationship between K m $$ {K}_{\mathrm{m}} $$ and E T $$ {E}_{\mathrm{T}} $$ , unlike the MM model. Thus, our framework becomes much more accurate than the original FDA equation, especially for drugs with high affinities, such as midazolam or strong inducers, such as rifampicin, so that the ratio between K m $$ {K}_{\mathrm{m}} $$ and E T $$ {E}_{\mathrm{T}} $$ becomes low (i.e., the MM model is invalid). Our work greatly improves the prediction of clinical DDIs, which is critical to preventing drug toxicity and failure.

Application of physiologically-based pharmacokinetic model approach to predict pharmacokinetics and drug-drug interaction of rivaroxaban: A case study of rivaroxaban and carbamazepine

Rivaroxaban (RIV; Xarelto; Janssen Pharmaceuticals, Beerse, Belgium) is one of the direct oral anticoagulants. The drug is a strong substrate of cytochrome P450 (CYP) enzymes and efflux transporters. This study aimed to develop a physiologically-based pharmacokinetic (PBPK) model for RIV. It contained three hepatic metabolizing enzyme reactions (CYP3A4, CYP2J2, and CYP-independent) and two active transporter-mediated transfers (P-gp and BCRP transporters). To illustrate the performance of the developed RIV PBPK model on the prediction of drug-drug interactions (DDIs), carbamazepine (CBZ) was selected as a case study due to the high DDI potential. Our study results showed that CBZ significantly reduces the exposure of RIV. The area under the concentration-time curve from zero to infinity (AUC_inf ) of RIV was reduced by 35.2% (from 2221.3 to 1438.7 ng*h/ml) and by 25.5% (from 2467.3 to 1838.4 ng*h/ml) after the first dose and at the steady-state, respectively, whereas the maximum plasma concentration (C_max ) of RIV was reduced by 37.7% (from 266.3 to 166.1 ng/ml) and 36.4% (from 282.3 to 179.5 ng/ml), respectively. The developed PBPK model of RIV could be paired with PBPK models of other interested perpetrators to predict DDI profiles. Further studies investigating the extent of DDI between CBZ and RIV should be conducted in humans to gain a full understanding of their safety and effects.

Physiologically-based pharmacokinetic modeling to predict drug interactions of lemborexant with CYP3A inhibitors

Lemborexant, a recently approved dual orexin receptor antagonist for treatment of adults with insomnia, is eliminated primarily by cytochrome P450 (CYP)3A metabolism. The recommended dose of lemborexant is 5 mg once per night, with a maximum recommended dose of 10 mg once daily. A physiologically-based pharmacokinetic (PBPK) model for lemborexant was developed and applied to integrate data obtained from in vivo drug-drug interaction (DDI) assessments, and to further explore lemborexant interaction with CYP3A inhibitors and inducers. The model predictions were in good agreement with observed pharmacokinetic data and with DDI results from clinical studies with CYP3A inhibitors, itraconazole and fluconazole. The model further predicted that DDI effects of weak CYP3A inhibitors (fluoxetine and ranitidine) are weak, and effects of moderate inhibitors (erythromycin and verapamil) are moderate. Based on the PBPK simulations and clinical efficacy and safety data, the maximum daily recommended lemborexant dose when administered with weak CYP3A inhibitors is 5 mg; co-administration of moderate and strong inhibitors should be avoided except in countries where 2.5 mg has been approved.

Inter-phenotypic differences in CYP2C9 and CYP2C19 metabolism: Bayesian meta-regression of human population variability in kinetics and application in chemical risk assessment

Quantifying variability in pharmacokinetics (PK) and toxicokinetics (TK) provides a science-based approach to refine uncertainty factors (UFs) for chemical risk assessment. In this context, genetic polymorphisms in cytochromes P450 (CYPs) drive inter-phenotypic differences and may result in reduction or increase in metabolism of drugs or other xenobiotics. Here, an extensive literature search was performed to identify PK data for probe substrates of the human polymorphic isoforms CYP2C9 and CYP2C19. Relevant data from 158 publications were extracted for markers of chronic exposure (clearance and area under the plasma concentration-time curve) and analysed using a Bayesian meta-regression model. Enzyme function (EF), driven by inter-phenotypic differences across a range of allozymes present in extensive and poor metabolisers (EMs and PMs), and fraction metabolised (Fm), were identified as exhibiting the highest impact on the metabolism. The Bayesian meta-regression model provided good predictions for such inter-phenotypic differences. Integration of population distributions for inter-phenotypic differences and estimates for EF and Fm allowed the derivation of CYP2C9- and CYP2C19-related UFs which ranged from 2.7 to 12.7, and were above the default factor for human variability in TK (3.16) for PMs and major substrates (Fm >60%). These results provide population distributions and pathway-related UFs as conservative in silico options to integrate variability in CYP2C9 and CYP2C19 metabolism using in vitro kinetic evidence and in the absence of human data. The future development of quantitative extrapolation models is discussed with particular attention to integrating human in vitro and in vivo PK or TK data with pathway-related variability for chemical risk assessment.

Application of unbound liver-to-plasma concentration ratio to quantitative projection of cytochrome P450-mediated drug-drug interactions using physiologically based pharmacokinetic modelling approach

1. This study evaluated the prediction accuracy of cytochrome P450 (CYP)-mediated drug-drug interaction (DDI) using minimal physiologically-based pharmacokinetic (PBPK) modelling incorporating the hepatic accumulation factor of an inhibitor (i.e. unbound liver/unbound plasma concentration ratio [K_p,uu,liver]) based on 22 clinical DDI studies. 2. K_p,uu,liver values were estimated using three methods: (1) ratio of cell-to-medium ratio in human cryopreserved hepatocytes (C/M_u) at 37 °C to that on ice (K_p,uu,C/M), (2) multiplication of total liver/unbound plasma concentration ratio (K_p,u,liver) estimated from C/M_u at 37 °C with unbound fraction in human liver homogenate (K_p,uu,cell) and (3) observed K_p,uu,liver in rats after intravenous infusion (K_p,uu,rat). 3. PBPK model using each K_p,uu,liver projected the area under the curve (AUC) increase of substrates more accurately than the model assuming a K_p,uu,liver of 1 for the average fold error and root mean square error did. Particularly, the model with a K_p,uu,liver of 1 underestimated the AUC increase of triazolam following co-administration with CYP3A4 inhibitor itraconazole by five-fold, whereas the AUC increase projected using the model incorporating the K_p,uu,C/M, K_p,uu,cell, or K_p,uu,rat of itraconazole and hydroxyitraconazole was within approximately two-fold of the actual value. 4. The results indicated that incorporating K_p,uu,liver into the PBPK model improved the accuracy of DDI projection.

Quantitative prediction of the extent of drug-drug interaction using a physiologically based pharmacokinetic model that includes inhibition of drug metabolism determined in cryopreserved hepatocytes

1. A physiologically based pharmacokinetic (PBPK) model that includes inhibition constant evaluated in cryopreserved hepatocytes was used to predict drug-drug interactions (DDIs) between orally administered nifedipine, a CYP substrate, and fluconazole or ketoconazole, CYP inhibitors, in rats. 2. The Kp,uu, ratio of unbound inhibitor concentration in liver ([I]liver,u) to that in plasma ([I]sys,u), of fluconazole and ketoconazole was 1.0 and 13.0, indicating that ketoconazole accumulates in liver. The ratios of inhibition constants in rat liver microsomes (Ki,mic,u) to that in rat cryopreserved hepatocytes (Ki,hep,u) for fluconazole and ketoconazole were 1.5 and 25.5, which were similar to the Kp,uu and suggested that cryopreserved hepatocytes could mimic the hepatic accumulation of inhibitors. 3. The increases in AUC of nifedipine predicted by the minimal PBPK model using [I]liver,u/Ki,mic,u and [I]sys,u/Ki,hep,u were within 1.5-fold of the observed values for both inhibitors, whereas the model using [I]sys,u/Ki,mic,u underestimated the AUC increase caused by ketoconazole 21-fold. 4. These results indicated that hepatic accumulation factor of an inhibitor is required for a precise DDI projection and that cryopreserved hepatocytes would be useful to obtain the Ki including hepatic accumulation factor. It was demonstrated that PBPK model using Ki,hep,u could be a valuable approach for quantitative DDI projection.

Comparison of the static in vivo approach to a physiologically based pharmacokinetic approach for metabolic drug–drug interactions prediction

Background: The in vivo mechanistic static model (IMSM) and the physiologically based pharmacokinetic (PBPK) model are two approaches used to predict the magnitude of drug–drug interactions (DDIs). The aim of this study was to evaluate the performance of IMSM and to compare IMSM with the PBPK approach implemented in Simcyp. Methods: The predictive performances of IMSM were evaluated on a panel of 628 DDIs. Subsequently, the IMSM and PBPK approaches were compared on a set of 104 DDIs. Results: The IMSM yielded 85% of predictions within 1.5-fold of the observed value on the 628 DDIs panel. The predictive performances of IMSM were better than those of the PBPK approach (median fold error 1 vs 0.86 on 104 studies; p = 0.02). Conclusion: The IMSM approach is an alternative tool for metabolic DDIs prediction.

Quantitative Systems Pharmacology: A Framework for Context

Quantitative Systems Pharmacology (QSP) is receiving increased attention. As the momentum builds and the expectations grow it is important to (re)assess and formalize the basic concepts and approaches. In this short review, I argue that QSP, in addition to enabling the rational integration of data and development of complex models, maybe more importantly, provides the foundations for developing an integrated framework for the assessment of drugs and their impact on disease within a broader context expanding the envelope to account in great detail for physiology, environment and prior history. I articulate some of the critical enablers, major obstacles and exciting opportunities manifesting themselves along the way. Charting such overarching themes will enable practitioners to identify major and defining factors as the field progressively moves towards personalized and precision health care delivery.

Identification and Mechanistic Investigation of Drug-Drug Interactions Associated With Myopathy: A Translational Approach

Myopathy is a group of muscle diseases that can be induced or exacerbated by drug–drug interactions (DDIs). We sought to identify clinically important myopathic DDIs and elucidate their underlying mechanisms. Five DDIs were found to increase the risk of myopathy based on analysis of observational data from the Indiana Network of Patient Care. Loratadine interacted with simvastatin (relative risk 95% confidence interval [CI] = [1.39, 2.06]), alprazolam (1.50, 2.31), ropinirole (2.06, 5.00), and omeprazole (1.15, 1.38). Promethazine interacted with tegaserod (1.94, 4.64). In vitro investigation showed that these DDIs were unlikely to result from inhibition of drug metabolism by CYP450 enzymes or from inhibition of hepatic uptake via the membrane transporter OATP1B1/1B3. However, we did observe in vitro synergistic myotoxicity of simvastatin and desloratadine, suggesting a role in loratadine–simvastatin interaction. This interaction was epidemiologically confirmed (odds ratio 95% CI = [2.02, 3.65]) using the data from the US Food and Drug Administration Adverse Event Reporting System.

Use of physiologically based pharmacokinetic modeling for assessment of drug-drug interactions

Interactions between co-administered medicines can reduce efficacy or lead to adverse effects. Understanding and managing such interactions is essential in bringing safe and effective medicines to the market. Ideally, interaction potential should be recognized early and minimized in compounds that reach late stages of drug development. Physiologically based pharmacokinetic models combine knowledge of physiological factors with compound-specific properties to simulate how a drug behaves in the human body. These software tools are increasingly used during drug discovery and development and, when integrating relevant in vitro data, can simulate drug interaction potential. This article provides some background and presents illustrative examples. Physiologically based models are an integral tool in the discovery and development of drugs, and can significantly aid our understanding and prediction of drug interactions.

Application of a systems approach to the bottom-up assessment of pharmacokinetics in obese patients: expected variations in clearance

BACKGROUND AND OBJECTIVES

The maintenance dose of a drug is dependent on drug clearance, and thus any biochemical and physiological changes in obesity that affect parameters such as cardiac output, renal function, expression of drug-metabolizing enzymes and protein binding may result in altered clearance compared with that observed in normal-weight subjects (corrected or uncorrected for body weight). Because of the increasing worldwide incidence of obesity, there is a need for more information regarding the optimal dosing of drug therapy to be made available to prescribers. This is usually provided via clinical studies in obese people; however, such studies are not available for all drugs that might be used in obese subjects. Incorporation of the relevant physiological and biochemical changes into predictive bottom-up pharmacokinetic models in order to optimize dosage regimens may offer a logical way forward for the cases where no clinical data exist. The aims of the current report are to apply such a 'systems approach' to identify the likelihood of observing variations in the clearance of drugs in obesity and morbid obesity for a set of compounds for which clinical data, as well as the necessary in vitro information, are available, and to provide a framework for assessing other drugs in the future.

METHODS

The population-specific changes in demographic, physiological and biochemical parameters that are known to be relevant to obese and morbidly obese subjects were collated and incorporated into two separate population libraries. These libraries, together with mechanistic in vitro-in vivo extrapolations (IVIVE) within the Simcyp Population-based Simulator™, were used to predict the clearance of oral alprazolam, oral caffeine, oral chlorzoxazone, oral ciclosporin, intravenous and oral midazolam, intravenous phenytoin, oral theophylline and oral triazolam. The design of the simulated studies was matched as closely as possible with that of the clinical studies. Outcome was measured by the predicted ratio of the clearance of the drug in obese and lean subjects ± its 90% confidence interval, compared with observed values. The overall statistical measures of the performance of the model to detect differences in compound clearance between obese and lean populations were investigated by measuring sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV). A power calculation was carried out to investigate the impact of the sample size on the overall outcome of clinical studies.

RESULTS

The model was successful in predicting clearance in obese subjects, with the degree to which simulations could mimic the outcome of in vivo studies being greater than 60% for six of the eight drugs. A clear difference in the clearance of chlorzoxazone was correctly picked up via simulation. The overall statistical measures of the performance of the Simcyp Simulator were 100% sensitivity, 66% specificity, 60% PPV and 100% NPV. Studies designed on the basis of the ratio of the absolute values required substantial numbers of participants in order to detect a significant difference, except for phenytoin and chlorzoxazone, where the ratios of the weight-normalized clearances generally showed statistically significant differences with a smaller number of subjects.

CONCLUSION

Extension of a mechanistic predictive pharmacokinetic model to accommodate physiological and biochemical changes associated with obesity and morbid obesity allowed prediction of changes in drug clearance on the basis of in vitro data, with reasonable accuracy across a range of compounds that are metabolized by different enzymes. Prediction of the effects of obesity on drug clearance, normalized by various body size scalars, is of potential value in the design of clinical studies during drug development and in the introduction of dosage adjustments that are likely to be needed in clinical practice.

In vitro-to-in vivo predictions of drug-drug interactions involving multiple reversible inhibitors

INTRODUCTION

Predictions of drug-drug interactions (DDIs) are commonly performed for single inhibitors, but interactions involving multiple inhibitors also frequently occur. Predictions of such interactions involving stereoisomer pairs, parent/metabolite combinations and simultaneously administered multiple inhibitors are increasing in importance. This review provides the framework for predicting inhibitory DDIs of multiple inhibitors with any combination of reversible inhibition mechanism.

AREAS COVERED

The review provides an overview of the reliability of the in vitro determined reversible inhibition mechanism. Furthermore, the article provides a method to predict DDIs for multiple reversible inhibitors that allows substituting the inhibition constant (K(i)) with an inhibitor affinity (IC(50)) value determined at S << K(M).

EXPERT OPINION

A better understanding and the prediction methods of DDIs, resulting from multiple inhibitors, are important. The inhibition mechanism of a reversible inhibitor is often equivocal across studies and unreliable. Determination of the K(i) requires the assignment of reversible inhibition mechanism but in vitro-to-in vivo prediction of DDI risk can be achieved for multiple inhibitors from estimates of the inhibitor affinity (IC(50)) only, regardless of the inhibition mechanism.