![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
School of Pharmacy and Pharmaceutical Sciences, University of Manchester, Manchester, United Kingdom (K.I., D.H., J.B.H.); and Department of Pharmacokinetics, Pharmacodynamics, and Drug Metabolism, Pfizer Global Research and Development, Groton, Connecticut (R.S.O.)
(Received January 14, 2005; Accepted March 9, 2005)
| Abstract |
|---|
|
|
|---|
A statistically significant increase in the area under the plasma concentration-time curve (AUC) of the victim drug in the presence of inhibitor (administered by multiple oral dosing) is the commonly used metric for an in vivo drug-drug interaction. Assuming that hepatic metabolism of the victim drug by a single pathway is the only route of elimination, the AUC ratio for an orally administered substrate in the presence and absence of inhibitor reflects the ratio of the control and inhibited intrinsic clearances. It is also assumed that the conditions of the "well stirred" liver model apply and that there is no impact of inhibitor on either the fraction absorbed (Fa) or the fraction unbound in plasma (fu) for the substrate. Thus, in theory, the change in clearance, and hence AUC, may be predicted from the inhibitor concentration available to the hepatic enzymes in vivo [I] and the in vitro Ki (Ito et al., 1998
, 2004
; Yao and Levy 2002
).
![]() | (1) |
![]() | (2) |
Although the use of eq. 2 to describe the degree of in vivo interaction between two drugs is widespread (von Moltke et al., 1994
; Ito et al., 1998
; Davit et al., 1999
; Lin, 2000
; Rodrigues, 2001; Tucker et al., 2001
; Yao and Levy, 2002
) and Ki values can be readily obtained from human liver microsomal studies, there is no general agreement over the selection of surrogate inhibitor concentration to reflect [I] (Tucker et al., 2001
; Ito et al., 2002
). In a recent analysis (Ito et al., 2004
), the use of hepatic input concentration gave the best resolution between positive and negative interactions and, most importantly, eliminated false negatives. It was also noted that there was a high incidence of false positive predictions, and many true positives were substantially over-predicted. One possible reason for these limitations is the assumption that there is a single metabolic pathway responsible for the elimination of the victim drug. In practice, this is frequently not the case; renal clearance of unchanged drug, parallel pathways of metabolism, and involvement of multiple P450s for the same reaction may all contribute to elimination. Thus, incorporation of the fraction metabolized by the particular P450 enzyme subject to inhibition (fmCYPi) may be a critically important consideration. Although the kinetic consequences of parallel elimination pathways have been alluded to by various investigators (Rowland and Matin, 1973
; Shaw and Houston, 1987
; Ito et al., 1998
; Yao and Levy, 2000; Rodrigues et al., 2001
), there has been no comprehensive evaluation of its impact. In general, accurate values of fmCYPi in vivo are not known, but in the case of CYP2D6 substrates, substantial information is often available from genetic polymorphism studies. By comparison of the poor and extensive phenotypes, a direct measure of the importance of CYP2D6 (that is, fmCYP2D6) can be made. Previously, we (Ito et al., 2004
) collated a number of in vivo studies involving CYP2D6 inhibition to assess the prediction of DDIs from the [I]/Ki ratio (see Fig. 1). These data are now used to explore the impact of parallel pathways on the prediction of DDIs. We also investigate the theoretical basis of the significance of parallel pathway and multiple P450 involvement and provide simulations of DDI based on a simple pharmacokinetic model. We demonstrate that prediction of DDIs from inhibitor characteristics (namely, [I] and Ki) is substantially improved by taking into account the fmCYPi for the victim drug.
|
| Materials and Methods |
|---|
|
|
|---|
![]() | (3) |
Case 1: Only One Pathway Subject to Inhibition. If the inhibitor does not affect both metabolic pathways, the intrinsic clearance in the presence of the inhibitor (CLintI) can be expressed as follows, again assuming linear metabolism, since the substrate concentration does not approach the Km and hence the distinction between competitive and noncompetitive inhibition mechanisms is not required:
![]() | (4) |
![]() | (5) |
Combining the eqs. 1 and 5, the following equation describes the AUC increase by the inhibitor:
![]() | (6) |
Case 2: Both Pathways Subject to Inhibition. If the inhibitor also inhibits the second pathway of substrate, then eq. 4 can be arranged as follows:
![]() | (7) |
![]() | (8) |
Simulations of AUC Ratio. First, the effect of a parallel elimination route was simulated according to eq. 6 for a range of values for fmCYP2D6 from 0.5 to 1. The effect of a parallel pathway subject to inhibition was simulated according to eq. 8, maintaining the fmCYP2D6 value at 0.8 and changing the relative affinity of the inhibitor for the two P450s, expressed as Ki,2/Ki, ratios of 1 to 10,000.
The CYP2D6 Paradigm: Estimation of fmCYP2D6 from Phenotype Studies and by Regression. CYP2D6 is known to show genetic polymorphism, and the clearance for CYP2D6 substrates in the extensive metabolizers (EMs) and poor metabolizers (PMs) can be expressed by eqs. 9 and 10, respectively, assuming that the CYP2D6 pathway is absent in the PMs:
![]() | (9) |
![]() | (10) |
![]() | (11) |
|
To obtain an alternative value of fmCYP2D6, the interaction data for each CYP2D6 substrate (see Fig. 1) were fitted to eq. 6 by nonlinear least-squares regression and the values obtained were compared to fmCYP2D6 from phenotype studies. For the interaction data involving desipramine as substrate, a regression analysis was also performed based on eq. 8 to obtain the best-fit for the Ki,2/Ki,1 ratio with the fixed value of fmCYP2D6 at 0.8. The total hepatic input concentration ([I]in) was used as the inhibitor concentration for all the predictions as this was found to be the most successful in qualitative zoning of DDIs (Ito et al., 2004
).
Data Source. Our original database (Ito et al., 2004
) contained 58 studies involving 14 substrates and 19 inhibitors. Overall, five substrates were excluded from the analysis due to the lack of phenotype studies or if they involved paroxetine, a mechanism-based inhibitor of CYP2D6 (Bertelsen et al., 2003
). Thus, 44 studies were used in the current analysis involving nine substrates and 12 reversible inhibitors.
The in vitro Ki values were previously collected (Ito et al., 2004
) and were all obtained against selective CYP2D6 probe reactions; in 25 cases the same probe was used in vitro and in vivo.
| Results |
|---|
|
|
|---|
![]() | (12) |
|
The shifts in both the position and the shape of the AUC ratio relationship to [I]/Ki illustrated in Fig. 2, A and B, result in a movement of the borderline between negative and positive DDIs. For the case of fmCYP2D6 equals 1 (eq. 2), an [I]/Ki ratio <1 would represent a low risk (resulting in AUC ratio <2), whereas the [I]/Ki >1 case would be a high risk (AUC ratio >2). However, the critical [I]/KI value that separates negative from positive interactions increases from 1 as fmCYP2D6 is reduced, as can be seen when eq. 6 is solved for an AUC ratio of 2.
![]() | (13) |
For fmCYP2D6 >0.75, the critical [1]/Ki value is relatively insensitive (between 1 and 2); however, as fmCYP2D6 decreases further, there is a marked increase in the borderline value as the limit of 0.5 is approached (see Fig. 2A, inset).
The CYP2D6 Paradigm. The effect of parallel pathways was first investigated in the interaction cases involving desipramine as the victim drug, since this was the substrate with the largest number of in vivo studies in our database (n = 16). Figure 3A includes studies showing an AUC increase in desipramine after both desipramine and imipramine dosing. The heavy broken line in Fig. 3A is the simulation line using eq. 6 with the fmCYP2D6 estimated from PM data. The light broken line is the regression line based on the same equation. The fmCYP2D6 values obtained from the PM data and regression analysis were comparable, with values of 0.88 and 0.8, respectively. A similar analysis using eq. 6 was performed for metoprolol (Fig. 4), the second most studied victim drug in our database. Although there were 10 studies, no major interactions were observed (AUC ratio ranged from 2.02 to 2.38). The fmCYP2D6 values obtained were 0.83 and 0.7 from PM data and regression approaches, respectively.
|
|
In addition, the impact of the inhibition of the second pathway was investigated using the desipramine data. The Ki ratio for the inhibition of both pathways was estimated by regression using eq. 8 and fixing the fmCYP2D6 value at 0.8. The line agrees well with the observed data and corresponds to inhibition of the second pathway that is 360 times weaker than the major pathway (Fig. 3B).
The effect of incorporating fmCYP2D6 in the prediction of AUC ratios for 44 DDI studies is shown in Fig. 5 as a three-dimensional surface. The fmCYP2D6 values for each substrate estimated by regression analysis were in good agreement with the PM-derived data (r = 0.748, p < 0.02; Table 1). For the nine substrates, a range of fmCYP2D6 values of 0.25 to 0.89 was found. The simulated surface is steepest as the fmCYP2D6 decreases from 1 to 0.8 and continues to decrease over the 0.8 to 0.5 range; below 0.5, all AUC ratios are less than 2. This indicates that many DDI over-predictions can be explained by taking parallel pathways into consideration and that false positives will occur when fmCYP2D6 is less than 0.5. This is clearly illustrated by the propranolol studies (fmCYP2D6 0.37).
|
|
| Discussion |
|---|
|
|
|---|
There have been numerous reports in the literature on the use of the [I]/Ki ratio approach (usually assuming a fmCYPi of 1), and they have been of mixed success. Recently, we (Ito et al., 2004
) compiled a database of 193 DDIs for which in vitro data were available in addition to in vivo inhibitor concentrations. The study showed good agreement on a qualitative level between the prediction from in vitro data and clinical observations. This systematic approach highlighted the importance of the selection of inhibitor concentration, whether this is an average systemic circulating concentration, a maximum hepatic portal vein concentration during the absorption process, a total drug concentration, or an unbound concentration. In terms of qualitative zoning of positive and negative DDI, the most successful estimate of [I] was the total input concentration, i.e., the total concentration of an inhibitor in the hepatic portal vein (Kanamitsu et al., 2000
). It was demonstrated that once mechanism-based inhibitors were excluded from the analysis, no false negative results were obtained. In terms of predicting positive interactions, more varied success was achieved, with a number of false positives and the majority of true positives being over-predicted. It was clear from this analysis that in order to predict DDIs quantitatively, additional considerations were necessary; hence the need for the present study that systematically explores the fraction metabolized by a particular cytochrome P450 enzyme (fmCYPi) as a substrate-dependent parameter.
The potential significance of parallel pathways of drug elimination in DDI prediction was pointed out as early as 1973 (Rowland and Matin, 1973
). In more recent years, reviewers of the area have also cited parallel elimination pathways and multiple P450 involvement as an important consideration (Yao and Levy 2000; Rodrigues et al., 2001
; Rostami-Hodjegan and Tucker, 2004
). However, there has been little substantial analysis to date to support this claim. We have selected a number of CYP2D6 DDI studies (n = 44) as appropriate test cases in view of the validity of the fmCYPi parameter obtained for these substrates. By using genetic phenotyping studies and comparing PM and EM data, we were able to obtain unequivocal estimates of the importance of CYP2D6 to the elimination of the substrates investigated. For the nine substrates studied, the fmCYP2D6 from PM studies varied from 0.37 (for propranolol) to 0.94 (for tolterodine). A series of simulations illustrated the importance of this range of fmCYP2D6 on the predicted AUC ratio. This finding was confirmed using information from 44 clinical studies, and we were able to increase success from 54% to 84% for positive and negative predictions and to reduce the number of over-predictions (as summarized in Fig. 5).
The present study highlights the critical importance of fmCYPi values >0.5. If the particular enzyme under examination contributes less than 50% of the total clearance, then it may be argued that there is little reason for concern over DDIs or, indeed, genetic polymorphisms. The maximum consequence of either phenomenon will be a doubling of AUC. Equation 12 indicates that the maximum effect can be predicted simply from the fmCYPi term. Although this outcome will be true for the vast majority of individuals, it is of interest to speculate on the possible impact on extreme individuals in the population for DDIs involving drugs whose metabolic clearance results from both a polymorphic enzyme (e.g., CYP2D6) and a nonpolymorphic enzyme (e.g., CYP3A4). Although inhibition of the polymorphic enzyme may be of minimal consequence when fmCYPi is low (say 25%), inhibition of the nonpolymorphic enzyme may cause substantial effects and, more importantly, these may be very varied. The latter situation will result from the dependence of clearance on the activity of the polymorphic enzyme in the absence of the major route of metabolite clearance. Hence the importance of genetic polymorphism in DDIs may not be limited to direct inhibition effects.
Certain inhibitors in the CYP2D6 DDI database used (e.g., fluoxetine, sertraline; see Ito et al., 2004
) have inhibitory metabolites; thus, the consequence of multiple inhibitors needs to be considered. As indicated (Rostami-Hodjegan and Tucker, 2004
), it is important to distinguish whether the metabolite, or indeed any second inhibitor, will act on the same or different enzyme. If the former situation occurs, then the cumulative effect will be considerably less than in the latter situation. In our simulations, we were able to identify the consequences of a secondary inhibitory effect on the clearance of desipramine (see Fig. 3B) despite only minor potency on the second enzyme (indicated by a Ki ratio of 360).
Application beyond CYP2D6 to other polymorphic P450s (e.g., CYPC19, CYP2C9) should be readily achievable either by the method used here (comparison of PM and EM data) or by use of selective inhibitors in EM subjects to "phenocopy" the PM status. Extending this work to other P450s that do not show polymorphic metabolism is more problematic due to the need for an estimate of fmCYPi. The total fraction metabolized can be readily determined once the fraction of drug excreted unchanged via either renal clearance or other mechanisms is known. Considerable effort, including the use of radiolabeled drugs, is routinely expended to resolve issues of mass balance, oral absorption, and fecal/biliary excretion. More specific information from urinary recovery of drug metabolites may be obtained provided the complexity of the metabolic pathways allows this resolution. If sequential metabolism occurs from more than one primary metabolite to a common secondary metabolite, it may be impossible to quantitatively assess the importance of particular drug clearance routes. However, the key information needed is assignment of particular enzymes to those drug pathways; that is, fmCYPi rather than fraction metabolized by a particular pathway. The importance of this specific information is being increasingly realized (Bjornsson et al., 2003
; Williams et al., 2003
). The use of data from reaction phenotyping may be very useful in this area, but to date, this has tended to be more of a qualitative measure than a quantitative one (McGinnity et al., 2000
; Lu et al., 2003
). However, the success apparent with CYP2D6 substrates would indicate that the above are worthy of detailed consideration, since even approximate fmCYPi values may markedly improve a prediction. The use of fuzzy set theory (Nestorov et al., 2002
) may provide a methodology to incorporate semiquantitative or qualitative estimates of fmCYPi.
The total hepatic input concentration was found to be the most useful estimate of [I] to predict the AUC ratio from the in vitro Ki value. This apparent need to ignore the consequences of plasma protein binding is surprising but similar to our previous experience with a larger database (Ito et al., 2004
). Incorporation of fmCYP2D6 improves the identification of true positive and negative DDIs from 54% to 84%. Consideration of the false positives and negatives (Table 2) are particularly informative, the former reducing from 20 of 44 to 4 of 44. However, there are some false negatives (3 of 44) when fmCYP2D6 is included in the prediction and these particular studies deserve further comment. One study involves inhibition of propranolol by propafenone (observed AUC ratio of 2.13) and two studies inhibition of imipramine by fluoxetine and fluvoxamine (observed AUC ratio of 3.33 and 3.63, respectively). For both substrates, fmCYP2D6 is less than 0.5, and hence, AUC ratios would not have been predicted to be greater than 2. It is of interest that the eight other propranolol studies and the three other imipramine studies that are part of this 2D6 DDI compilation are true negative DDIs. Therefore, inaccurate fmCYP2D6 values appear not to be the explanation for these three anomalies. Inhibition of the other P450s responsible for the metabolism of both substrates (CYP1A2, CYP2C19, CYP2C9, and CYP3A4) is the likely explanation. For fluoxetine and fluvoxamine, this has been comprehensively documented (Olesen and Linnet, 2000
; Hemeryck and Belpaire, 2002
) but propafenone studies appear to have been limited to CYP1A2 (Kobayashi et al., 1998
).
In conclusion, consideration of parallel elimination pathways and multiple P450 involvement provides extra information that improves prediction of CYP2D6 DDI. Incorporation of fmCYP2D6, readily available for CYP2D6 from PM and EM data, reduces the number of false positives and also the extent of quantitative over-predictions of true positives. This provides a valuable step forward in making quantitative predictions of DDIs, which, to date, have placed emphasis on delineating the inhibitor characteristics (namely [I] and Ki) and relatively little attention on the victim drug properties.
| Footnotes |
|---|
Article, publication date, and citation information can be found at http://dmd.aspetjournals.org.
ABBREVIATIONS: DDI, drug-drug interaction; AUC, area under the plasma concentration-time curve; CLint, intrinsic clearance; EM, extensive metabolizer; fmCYPi, fraction metabolized by a particular P450 enzyme i; [I], concentration of inhibitor available to the enzyme; PM, poor metabolizer.
Address correspondence to: Professor J. Brian Houston, School of Pharmacy and Pharmaceutical Sciences, University of Manchester, Oxford Road, Manchester M13 9PL, UK. E-mail: Brian.Houston{at}manchester.ac.uk
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. A. Rock, R. S. Foti, and J. T. Pearson The Combination of Chemical and Antibody Inhibitors for Superior P450 3A Inhibition in Reaction Phenotyping Studies Drug Metab. Dispos., December 1, 2008; 36(12): 2410 - 2413. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Kotsuma, T. Tokui, S. Freudenthaler, and K. Nishimura Effects of Ketoconazole and Quinidine on Pharmacokinetics of Pactimibe and Its Plasma Metabolite, R-125528, in Humans Drug Metab. Dispos., August 1, 2008; 36(8): 1505 - 1511. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Lu, P. Hatsis, C. Berg, F. W. Lee, and S. K. Balani Prediction of Pharmacokinetic Drug-Drug Interactions Using Human Hepatocyte Suspension in Plasma and Cytochrome P450 Phenotypic Data. II. In Vitro-in Vivo Correlation with Ketoconazole Drug Metab. Dispos., July 1, 2008; 36(7): 1255 - 1260. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Lu, C. Berg, S. R. Prakash, F. W. Lee, and S. K. Balani Prediction of Pharmacokinetic Drug-Drug Interactions Using Human Hepatocyte Suspension in Plasma and Cytochrome P450 Phenotypic Data. III. In Vitro-in Vivo Correlation with Fluconazole Drug Metab. Dispos., July 1, 2008; 36(7): 1261 - 1266. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J. Reese, R. M. Wurm, K. T. Muir, G. T. Generaux, L. St. John-Williams, and D. J. Mcconn An in Vitro Mechanistic Study to Elucidate the Desipramine/Bupropion Clinical Drug-Drug Interaction Drug Metab. Dispos., July 1, 2008; 36(7): 1198 - 1201. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. F. McGinnity, N. J. Waters, J. Tucker, and R. J. Riley Integrated in Vitro Analysis for the in Vivo Prediction of Cytochrome P450-Mediated Drug-Drug Interactions Drug Metab. Dispos., June 1, 2008; 36(6): 1126 - 1134. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Kotsuma, T. Tokui, T. Ishizuka-Ozeki, T. Honda, H. Iwabuchi, T. Murai, T. Ikeda, and H. Saji CYP2D6-Mediated Metabolism of a Novel Acyl Coenzyme A:Cholesterol Acyltransferase Inhibitor, Pactimibe, and Its Unique Plasma Metabolite, R-125528 Drug Metab. Dispos., March 1, 2008; 36(3): 529 - 534. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Fairman, C. Collins, and S. Chapple Progress Curve Analysis of CYP1A2 Inhibition: A More Informative Approach to the Assessment of Mechanism-Based Inactivation? Drug Metab. Dispos., December 1, 2007; 35(12): 2159 - 2165. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Nakashima, H. Takama, Y. Ogasawara, T. Kawakami, T. Nishitoba, S. Hoshi, E. Uchida, and H. Tanaka Effect of Cinacalcet Hydrochloride, a New Calcimimetic Agent, on the Pharmacokinetics of Dextromethorphan: In Vitro and Clinical Studies J. Clin. Pharmacol., October 1, 2007; 47(10): 1311 - 1319. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Obach, R. L. Walsky, and K. Venkatakrishnan Mechanism-Based Inactivation of Human Cytochrome P450 Enzymes and the Prediction of Drug-Drug Interactions Drug Metab. Dispos., February 1, 2007; 35(2): 246 - 255. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Lu, G. T. Miwa, S. R. Prakash, L.-S. Gan, and S. K. Balani A Novel Model for the Prediction of Drug-Drug Interactions in Humans Based on in Vitro Cytochrome P450 Phenotypic Data Drug Metab. Dispos., January 1, 2007; 35(1): 79 - 85. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. L. Walsky, A. V. Astuccio, and R. S. Obach Evaluation of 227 Drugs for In Vitro Inhibition of Cytochrome P450 2B6. J. Clin. Pharmacol., December 1, 2006; 46(12): 1426 - 1438. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. P. Gibbs, R. Hyland, and K. Youdim Minimizing Polymorphic Metabolism in Drug Discovery: Evaluation of the Utility of in Vitro Methods for Predicting Pharmacokinetic Consequences Associated with CYP2D6 Metabolism Drug Metab. Dispos., September 1, 2006; 34(9): 1516 - 1522. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Y. Chien, A. Lucksiri, C. S. Ernest II, J. C. Gorski, S. A. Wrighton, and S. D. Hall STOCHASTIC PREDICTION OF CYP3A-MEDIATED INHIBITION OF MIDAZOLAM CLEARANCE BY KETOCONAZOLE Drug Metab. Dispos., July 1, 2006; 34(7): 1208 - 1219. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Galetin, H. Burt, L. Gibbons, and J. B. Houston PREDICTION OF TIME-DEPENDENT CYP3A4 DRUG-DRUG INTERACTIONS: IMPACT OF ENZYME DEGRADATION, PARALLEL ELIMINATION PATHWAYS, AND INTESTINAL INHIBITION Drug Metab. Dispos., January 1, 2006; 34(1): 166 - 175. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Obach, R. L. Walsky, K. Venkatakrishnan, E. A. Gaman, J. B. Houston, and L. M. Tremaine The Utility of in Vitro Cytochrome P450 Inhibition Data in the Prediction of Drug-Drug Interactions J. Pharmacol. Exp. Ther., January 1, 2006; 316(1): 336 - 348. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. F. McGinnity, J. Tucker, S. Trigg, and R. J. Riley PREDICTION OF CYP2C9-MEDIATED DRUG-DRUG INTERACTIONS: A COMPARISON USING DATA FROM RECOMBINANT ENZYMES AND HUMAN HEPATOCYTES Drug Metab. Dispos., November 1, 2005; 33(11): 1700 - 1707. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Zhang, T. J. Chando, D. W. Everett, C. J. Patten, S. S. Dehal, and W. G. Humphreys IN VITRO INHIBITION OF UDP GLUCURONOSYLTRANSFERASES BY ATAZANAVIR AND OTHER HIV PROTEASE INHIBITORS AND THE RELATIONSHIP OF THIS PROPERTY TO IN VIVO BILIRUBIN GLUCURONIDATION Drug Metab. Dispos., November 1, 2005; 33(11): 1729 - 1739. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||