Pharmacological Reviews Get Tables of Contents delivered automatically
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ito, K.
Right arrow Articles by Sugiyama, Y.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ito, K.
Right arrow Articles by Sugiyama, Y.

Vol. 50, Issue 3, 387-412, September 1998

Prediction of Pharmacokinetic Alterations Caused by Drug-Drug Interactions: Metabolic Interaction in the Liver

K. Ito, T. Iwatsubo, S. Kanamitsu, K. Ueda, H. Suzuki and Y. Sugiyamaa

Department of Pharmaceutics, Faculty of Pharmaceutical Sciences, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan.

I. Introduction
II. Drug-Drug Interactions Other Than Involving Metabolism
    A. Drug-drug interactions involving plasma protein binding
    B. Drug-drug interactions at the transport carrier
III. Drug-Drug Interactions Involving Metabolism in the Liver
    A. Examples of In Vivo Drug-Drug Interactions Involving P450 Metabolism
    B. Inhibition Mechanism of Drug Metabolism by P450
    C. Inhibition Patterns of Drug Metabolism
        1. Competitive Inhibition.
        2. Noncompetitive Inhibition.
        3. Uncompetitive Inhibition.
    D. Prediction of In Vivo Drug-Drug Interactions Based on In Vitro Data
        1. General Equations.
        2. The evaluation of the unbound concentration of the inhibitor in vivo.
    E. Examples of the Prediction of Drug-Drug Interactions Based on Literature Data
        1. Successful Cases of In Vitro/In Vivo Prediction.
        2. Interactions Predictable for the Objective Metabolic Pathway but not Predictable for the Overall Data.
        3. Interactions Not Predictable by In Vitro/In Vivo Scaling.
    F. Procedure for Predicting Inhibitory Effects of Coadministered Drugs on the Hepatic Metabolism of Other Drugs
    G. Mechanism-Based Inhibition
        1. Characteristics of Mechanism-Based Inhibition.
        2. Kinetic Analysis of Mechanism-Based Inhibition: Analysis of In Vitro Data.
        3. Prediction of In Vivo Interactions from In Vitro Data in the Case of Mechanism-Based Inhibition.
    H. Problems To Be Solved for the More Precise Prediction of Drug-Drug Interactions
        1. Estimation of the Tissue Unbound Concentration of the Inhibitor That Is Actively Transported into Hepatocytes.
        2. Evaluation of Drug-Drug Interactions Involving Drug Metabolism in the Gut.
References

    I. Introduction
Top
Next
References

Serious side-effects caused by drug interactions have attracted a great deal of attention and have become a social problem since the coadministration of ketoconazole and terfenadine was reported to cause potentially life-threatening ventricular arrhythmias (Monahan et al., 1990), and an interaction between sorivudine and fluorouracil resulted in fatal toxicity in Japan (Watabe, 1996; Okuda et al., 1997). The possible sites of drug-drug interaction which can change pharmacokinetic profiles include: (1) gastrointestinal absorption, (2) plasma and/or tissue protein binding, (3) carrier-mediated transport across plasma membranes (including hepatic or renal uptake and biliary or urinary secretion), and (4) metabolism. Pharmacodynamic interactions such as antagonism at the receptor may also increase or decrease the effects of a drug.

In this review, after brief comments on (2) and (3), we intend to focus on (4) and to discuss the possibility of the quantitative prediction of drug-drug interactions in vivo based on the analyses of data from literature obtained by in vitro experiments using human liver samples. Furthermore, strategic proposals for avoiding toxic interactions will be given from a pharmacokinetic point of view.

    II. Drug-Drug Interactions Other Than Involving Metabolism
Top
Previous
Next
References

A. Drug-drug interactions involving plasma protein binding

Although interactions involving plasma protein binding are well known, they rarely cause clinically serious problems (Rowland and Tozer, 1995; Rolan, 1994). The reasons are summarized below.

The unbound fraction (fu)b of a drug in plasma is increased when it is displaced by other drugs at the plasma protein binding sites. Subsequent alterations in plasma concentration profiles can be caused by changes in both clearance (CL) and volume of distribution (Vd) of the drug. The effect on the steady-state concentration (Css) and the area under concentration-time curve (AUC) can be predicted from the change in CL. It should be noted that the effect of protein binding replacement depends on the magnitude of CL and the route of administration. As shown in table 1, an analysis based on the well-stirred model has revealed that the protein binding replacement has little effect on the Css and AUC for unbound drugs (Cu,ss and AUCu) after oral administration, which are parameters directly related to the pharmacological and adverse effects, irrespective of the magnitude of CL. In the case of low clearance drugs, Cu,ss and AUCu after intravenous administration also are affected little by protein binding replacement. The only situation for a possible interaction is after the intravenous administration of a high clearance drug and there are few examples of this in clinical practice (Rolan, 1994).

                              
View this table:
[in this window]
[in a new window]
 
TABLE 1
Relationship between the area under concentration-time curve (AUC) or AUC for unbound drugs (AUCu) and the hepatic blood flow (Qh), the hepatic intrinsic clearance (CLint,h), and the blood unbound fraction (fb) based on the well-stirred model (Wilkinson, 1983)

The alteration of Vd caused by protein binding replacement also has an effect on the blood drug concentration (Rowland and Tozer, 1995). In the case of drugs with a relatively large Vd, Vd increases in parallel with fu. Although this leads to a transient reduction in total blood concentration caused by the redistribution of the drug into tissues, the unbound concentration is not affected. However, in the case of drugs with a small Vd, which depend on fu to a lesser extent, the total blood concentration is not affected so much by the change in fu, but the unbound concentration is greatly altered.

Figure 1 shows the simulation of the effects of protein binding replacement on the blood concentration profile during a constant intravenous infusion, where the protein binding and the tissue distribution of the drug are assumed to reach equilibrium rapidly, i.e., the concentration changes rapidly in response to a change in fu. In this simulation, changes in both CL and Vd associated with the change in fu were considered. As just described above, the steady-state unbound concentration is altered with the change in fu only for a high clearance drug. It is also clear from figure 1 that, in the case of drugs with a small Vd, a transient increase in the unbound concentration is observed even for a low clearance drug, and caution for the possible occurrence of side effects is needed.


View larger version (31K):
[in this window]
[in a new window]
 
Fig. 1.   Effects of protein binding replacement on the blood concentration profile during a constant intravenous infusion of (a) a high clearance drug with a high Vd, (b) a high clearance drug with a low Vd, (c) a low clearance drug with a high Vd, or (d) a low clearance drug with a low Vd. The protein binding and the tissue distribution of the drug are assumed to reach the equilibrium rapidly. An ideal situation is assumed where the concentration of the interacting drug (the displacer of the protein binding) immediately reaches a constant value at 10 min.

B. Drug-drug interactions at the transport carrier

Very few studies have focused on drug-drug interactions involving carrier-mediated transport across membranes, including the interactions involving renal secretion and reabsorption and those where p-glycoprotein (p-gp) plays a role (Tsuruo et al., 1981; Slater et al., 1986; Kusuhara et al., in press).

Along with metabolism, renal excretion is one of the most important processes affecting the total body clearance of a drug. Alterations in this process caused by drug-drug interactions should, therefore, be carefully considered. Secretion of drugs at the renal tubule is an active transport process, where organic anion transporters, organic cation transporters, and p-gp are known transport carriers (Hori et al., 1982; Takano et al., 1984; Tanigawara et al., 1992). The renal clearance of a drug is reduced by inhibition of these transport processes. It is known that both organic anion transporters and organic cation transporters exist on both the basolateral membrane (BLM) and the brush border membrane (BBM) and that they are different from each other, whereas p-gp is only present on the BBM. The inhibitors of these transporters interact with other drugs; for example, inhibition of the renal excretion of penicillin and other related drugs by probenecid (Hunter, 1951), methotrexate excretion by nonsteroidal anti-inflammatory drugs (Statkevich et al., 1993), and digoxin excretion by quinidine (Tanigawara et al., 1992) all involve this kind of interaction.

Most studies of pharmacokinetic drug-drug interactions reported so far have been limited to the analysis of hepatic metabolism. However, the hepatic clearance of many drugs has been found to be determined mainly by hepatic uptake (Yamazaki et al., 1995, 1996). The overall intrinsic clearance (CLint,all) can be expressed using the intrinsic clearance for metabolism (CLint) and that for membrane permeation (PS) as follows:
<UP>CL<SUB>int,all</SUB></UP>=<UP>PS<SUB>inf</SUB></UP> · <UP>CL<SUB>int</SUB>/</UP>(<UP>PS<SUB>eff</SUB></UP>+<UP>CL</UP><SUB><UP>int</UP></SUB>) (1)
where PSinf is intrinsic clearance for influx, and PSeff is intrinsic clearance for efflux.

It is clear from equation (1) that CLint,all equals CLint in the case of drugs with large (PS >>  CLint) and symmetrical (PSinf = PSeff) membrane permeability. Otherwise, hepatic clearance is affected by the membrane permeability of the drug. In such cases, it is important to evaluate drug-drug interactions involving not only metabolism but also membrane permeation. In our laboratory, several cases of drug-drug interactions were found in rats at the level of transporters involved in hepatobiliary transport as shown below. In the future, similar interactions at the transporter level possibly may be found in the clinical situation. The interactions found in rats include: inhibition of biliary excretion of glucuronides and sulfates of liquiritigenin, a flavonoid, by organic anions such as dibromosulfophthalein (DBSP) and glycyrrhizin, which has a glucuronide moiety (Shimamura et al., 1994); inhibition of biliary excretion of glycyrrhizin by DBSP (Shimamura et al., 1996); inhibition of biliary excretion of leukotriene C4, which has a glutathione moiety, by DBSP (Sathirakul et al., 1994); and reduction of plasma clearance, based on hepatic uptake and biliary excretion, of octreotide, a small octapeptide, by DBSP and taurocholate (Yamada et al., 1997). In vivo drug-drug interactions involving membrane transport remain to be predicted based on in vitro studies of the membrane permeability of drugs.

    III. Drug-Drug Interactions Involving Metabolism in the Liver
Top
Previous
Next
References

As a pharmacokinetic parameter directly related to the pharmacological and/or adverse effects of drugs, it is very important to predict the hepatic clearance. Because the use of animal scale-up is limited in the case of hepatic metabolic clearance due to large inherent interspecies differences, we have developed an alternative methodology to predict in vivo metabolic clearance in the liver; it is based on in vitro studies using mainly rat liver microsomes and isolated rat hepatocytes (Sugiyama and Ooie, 1993; Iwatsubo et al., 1996). Recently, with the greater availability of human liver samples, the method of in vitro/in vivo scaling can now be applied to human studies. We have already demonstrated that the method can be applied to P450 metabolism in humans based on in vitro and in vivo data obtained from the literature (Iwatsubo et al., 1997). However, the prediction of intrinsic clearance was not successful for some drugs, possibly because of the contribution of active transport into the liver and/or first-pass metabolism in the gut.

In order to prevent toxic drug-drug interactions, it is important to quantitatively predict pharmacokinetic changes caused by coadministration of drugs that are known to inhibit the hepatic metabolism of the drug under study (Sugiyama and Iwatsubo, 1996; Sugiyama et al., 1996). In this review, we have focused on the drug-drug interactions via inhibitory mechanisms and have tried to predict in vivo interactions from in vitro data on drug metabolism obtained from the literature.

A. Examples of In Vivo Drug-Drug Interactions Involving P450 Metabolism

Drug-drug interactions involving metabolism are one of the principal problems in clinical practice to evaluate the pharmacological and adverse effects of drugs. Parkinson (1996) summarized examples of substrates, inhibitors, and inducers of the major human liver microsomal P450 enzymes involved in drug metabolism. In the case of drugs that undergo metabolism by CYP3A4 and 2D6, particular attention should be paid to the interactions resulting in alterations in blood concentrations possibly accompanied by a change in its effects, because a number of drugs are metabolized by these enzymes (Bertz and Granneman, 1997). For example, blood concentrations of imipramine and desipramine, substrates for CYP2D6, are elevated several-fold by coadministration of fluoxetine, another substrate for CYP2D6 (Bergstrom et al., 1992). Similarly, concentrations of terfenadine, which is metabolized by CYP3A4, are increased in patients taking erythromycin, which is also a substrate for CYP3A4 (Honig et al., 1992). Quinidine is metabolized mainly by CYP3A4 but inhibits the metabolism of substrates for CYP2D6, such as sparteine, rather than those for CYP3A4 (Schellens et al., 1991). Furthermore, in the case of drugs whose metabolism is mediated by multiple isozymes (e.g., diazepam), drug-drug interaction may be complicated because of possible dose-dependent changes in the contribution of each isozyme to the overall metabolism (Iwatsubo et al., 1997).

B. Inhibition Mechanism of Drug Metabolism by P450

Drug metabolism by P450 can be inhibited by any of the following three mechanisms.

The first is mutual competitive inhibition caused by coadministration of drugs metabolized by the same P450 isozyme, such as the above-mentioned (see Sec. A.) combinations of imipramine or desipramine and fluoxetine (CYP2D6). In this case, as reported for metoprolol and propafenone (CYP2D6) (Wagner et al., 1987), blood concentrations of both drugs may be increased.

The second is the inactivation of P450 by the drug metabolite forming a complex with P450. This type of inhibition is designated as "mechanism-based inhibition" (Silverman, 1988). Inhibition by macrolide antibiotics, such as erythromycin, is a typical example of this type of interaction. As shown in figure 2, P450 demethylates and oxidizes the macrolide antibiotic into a nitrosoalkane that forms a stable, inactive complex with P450 (Periti et al., 1992).


View larger version (13K):
[in this window]
[in a new window]
 
Fig. 2.   Inhibition mechanism of P450 by macrolides (Periti et al., 1992).

The third is inhibition by the binding of imidazole or a hydrazine group to the haem portion of P450. In the case of cimetidine, the nitrogen in the imidazole ring binds to the haem portion of P450 causing nonselective inhibition of many P450 isozymes (Somogyi and Muirhead, 1987).

C. Inhibition Patterns of Drug Metabolism

The effects of inhibition of drug metabolism on in vivo pharmacokinetics are highly variable and depend on the properties of the drug, the route of administration, etc. (Rowland and Martin, 1973; Tucker, 1992). Except for the case of mechanism-based inhibition, inhibition of drug metabolism can be classified into the following three categories, and the equations corresponding to each inhibition type have been derived (Todhunter, 1979).

1. Competitive Inhibition. Competitive inhibition is a pattern of the inhibition where the inhibitor competes with the drug for the same binding site within an enzyme protein:
    <UP>E</UP>+<UP>S</UP><UP> ⇄ ES</UP> → <UP>E</UP>+<UP>P</UP>
<UP>K<SUB>m</SUB> </UP>(<UP>Michaelis constant for S</UP>)
<UP>E</UP>+<UP>I</UP><UP> ⇄ EI</UP>
<UP>K<SUB>i</SUB> </UP>(<UP>Inhibition constant for I</UP>)
where E is the enzyme, S is the substrate, ES is the enzyme-substrate complex, P is the product, I is the inhibitor, and EI is the enzyme-inhibitor complex. In the case of competitive inhibition, the metabolic rate (v) can be expressed by the following equation (2):
<UP>v = </UP><FR><NU><UP>V<SUB>max</SUB></UP> · <UP>S</UP></NU><DE><UP>K</UP><SUB><UP>m</UP></SUB>(<UP>1</UP>+<UP>I/K</UP><SUB><UP>i</UP></SUB>)+<UP>S</UP></DE></FR> (2)
where Vmax is the maximum metabolic rate.

It is clear from equation (2) that the inhibition by a given concentration of I is marked when the substrate concentration is low and becomes less marked with an increase in the substrate concentration.

2. Noncompetitive Inhibition. Noncompetitive inhibition is a pattern of inhibition where the inhibitor binds to the same enzyme as the drug but the binding site is different, resulting in a conformation change, etc., of the protein:
<AR><R><C><UP>E</UP>+<UP>S</UP></C><C><UP>⇄</UP></C><C><UP>ES → E</UP>+<UP>P</UP></C></R><R><C></C><C><UP>K</UP><SUB><UP>m</UP></SUB></C></R><R><C><UP>E</UP>+<UP>I</UP></C><C><UP>⇄</UP></C><C><UP>EI</UP></C></R><R><C></C><C><UP>K</UP><SUB><UP>i</UP></SUB></C></R><R><C><UP>EI</UP>+<UP>S</UP></C><C><UP>⇄</UP></C><C><UP>EIS</UP></C></R><R><C></C><C><UP>K</UP><SUB><UP>m</UP></SUB></C></R><R><C><UP>ES</UP>+<UP>I</UP></C><C><UP>⇄</UP></C><C><UP>EIS</UP></C></R><R><C></C><C><UP>K</UP><SUB><UP>i</UP></SUB></C></R></AR>
where EIS is the enzyme-inhibitor-substrate complex. It is assumed that the inhibitor binds to the free enzyme and the ES complex with the same affinity. In the case of noncompetitive inhibition, the metabolic rate can be expressed by the following equation (3):
<UP>v = </UP><FR><NU>{<UP>V<SUB>max</SUB>/</UP>(<UP>1</UP>+<UP>I/K</UP><SUB><UP>i</UP></SUB>)} · <UP>S</UP></NU><DE><UP>K<SUB>m</SUB></UP>+<UP>S</UP></DE></FR> (3)
It is clear from equation (3) that the degree of inhibition does not depend on the substrate concentration.

3. Uncompetitive Inhibition. Uncompetitive inhibition is a pattern of inhibition where the inhibitor binds only to the enzyme forming a complex with the drug:
<AR><R><C><UP>E</UP>+<UP>S</UP></C><C><UP>⇄</UP></C><C><UP>ES → E</UP>+<UP>P</UP></C></R><R><C></C><C><UP>K</UP><SUB><UP>m</UP></SUB></C></R><R><C><UP>ES</UP>+<UP>I</UP></C><C><UP>⇄</UP></C><C><UP>EIS</UP></C></R><R><C></C><C><UP>K</UP><SUB><UP>i</UP></SUB></C></R></AR>
Unlike competitive and noncompetitive inhibition, the inhibitor cannot bind to the free enzyme. In the case of uncompetitive inhibition, the metabolic rate can be expressed by the following equation (4):
<UP>v</UP> = <FR><NU>{<UP>V<SUB>max</SUB></UP>/(1+<UP>I/K<SUB>i</SUB></UP>)} · <UP>S</UP></NU><DE><UP>K<SUB>m</SUB></UP>/(1+<UP>I</UP>/<UP>K<SUB>i</SUB></UP>)+<UP>S</UP></DE></FR> (4)
It is clear from equation (4) that the inhibition becomes more marked with increasing substrate concentration.

The degree of inhibition depends on the inhibition pattern when the substrate concentration is high. However, when the substrate concentration is much lower than Km (Km >>  S), the degree of inhibition (R) is expressed by the following equation (5), independent of the inhibition pattern, except in the case of the uncompetitive inhibition (Tucker, 1992):
<UP>R</UP>=<FR><NU><UP>v</UP>(<UP>+inhibitor</UP>)</NU><DE><UP>v</UP>(<UP>−inhibitor</UP>)</DE></FR>=<FR><NU>1</NU><DE>1+<UP>I</UP>/<UP>K</UP><SUB><UP>i</UP></SUB></DE></FR> (5)
In clinical situations, the substrate concentration is usually much lower than Km. In this review, we will discuss the most frequently observed case in which equation (5) can be applied.

D. Prediction of In Vivo Drug-Drug Interactions Based on In Vitro Data

1. General Equations. The following factors determine the degree of change in Css and AUC caused by the drug-drug interaction in vivo:
1) The route of administration (intravenous or oral, i.e., whether the drug first passes through the liver or not).
2) Fraction (fh) of hepatic clearance (CLh) in total clearance (CLtot).
3) Fraction (fm) of the metabolic process subject to inhibition in CLh.
4) Unbound concentration of the inhibitor (Iu) around the enzyme.
5) Inhibition constant (Ki).
6) Plasma unbound concentration (Cu,ss) of the drug subject to inhibition.
7) Michaelis constant (Km) for the drug subject to inhibition.

fh and fm are expressed as follows:
<UP>f<SUB>h</SUB></UP>=<FR><NU><UP>CL</UP><SUB><UP>h</UP></SUB></NU><DE><UP>CL<SUB>h</SUB></UP>+<UP>CL</UP><SUB><UP>r</UP></SUB></DE></FR> (6)
<UP>f<SUB>m</SUB></UP>=<FR><NU><UP>CL</UP><SUB><UP>int,1</UP></SUB></NU><DE><UP>CL<SUB>int,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,2</UP></SUB></DE></FR> (7)
where CLh is hepatic clearance, CLr is renal clearance, and CLint,1 and CLint,2 represent the intrinsic clearance for the metabolic pathway inhibited and not inhibited by the inhibitor, respectively (CLint = CLint,1 + CLint,2). In equation (6), it is assumed that only the liver and kidney are the clearance organs. Equation (6) can be rearranged to give the following equation:
<UP>CL<SUB>r</SUB></UP>=<UP>CL<SUB>h</SUB></UP>(1/<UP>f<SUB>h</SUB></UP>−1) (8)
Equation (7) can be rearranged to give the following equation:
<UP>CL<SUB>int,2</SUB></UP>=<UP>CL<SUB>int,1</SUB></UP>(1/<UP>f<SUB>m</SUB></UP>−1) (9)
The fractional clearance for a particular metabolic pathway (CLh,m) can be expressed as fh multiplied by fm.

Rc, defined as the degree of increase in Css and AUC caused by the drug-drug interaction in vivo, can be calculated as shown below, depending on the route of administration.

a. INTRAVENOUS ADMINISTRATION. The change in AUC after intravenous bolus administration (AUCiv) and Css during intravenous infusion can be expressed by the following equation, if the dose or infusion rate is constant:

<UP>R<SUB>c</SUB></UP>=<FR><NU><UP>AUC</UP><SUB><UP>iv</UP></SUB>(<UP>+inhibitor</UP>)</NU><DE><UP>AUC</UP><SUB><UP>iv</UP></SUB>(<UP>−inhibitor</UP>)</DE></FR>=<FR><NU><UP>C</UP><SUB><UP>ss</UP></SUB>(<UP>+inhibitor</UP>)</NU><DE><UP>C</UP><SUB><UP>ss</UP></SUB>(<UP>−inhibitor</UP>)</DE></FR> (10)
=<FR><NU><UP>CL</UP><SUB><UP>tot</UP></SUB></NU><DE><UP>CL<SUB>tot</SUB>′</UP></DE></FR>=<FR><NU><UP>CL<SUB>h</SUB></UP>+<UP>CL</UP><SUB><UP>r</UP></SUB></NU><DE><UP>CL<SUB>h</SUB>′</UP>+<UP>CL</UP><SUB><UP>r</UP></SUB></DE></FR>=<FR><NU><UP>CL<SUB>h</SUB></UP>+<UP>CL</UP><SUB><UP>h</UP></SUB>(<UP>1/f<SUB>h</SUB> − 1</UP>)</NU><DE><UP>CL<SUB>h</SUB>′</UP>+<UP>CL</UP><SUB><UP>h</UP></SUB>(<UP>1/f<SUB>h</SUB> − 1</UP>)</DE></FR>
=<FR><NU><UP>CL<SUB>h</SUB>/f</UP><SUB><UP>h</UP></SUB></NU><DE><UP>CL<SUB>h</SUB>′</UP>+<UP>CL<SUB>h</SUB>/f<SUB>h</SUB></UP>−<UP>CL</UP><SUB><UP>h</UP></SUB></DE></FR>=<FR><NU><UP>1</UP></NU><DE><UP>f<SUB>h</SUB></UP> · <UP>CL<SUB>h</SUB>′/CL<SUB>h</SUB></UP>+<UP>1 − f</UP><SUB><UP>h</UP></SUB></DE></FR>
where ' represents the value after alteration by the drug-drug interaction.

i. High clearance drug. Because fb · CLint is much larger than the hepatic blood flow rate (Qh) (Qh <<  fb · CLint), CLh is rate-limited by the flow rate (CLh = Qh). When the altered CLh is still rate-limited by the flow rate (CLh' = Qh), i.e., Qh <<  fb. CLint', then CLh' equals CLh. Thus, Rc can be calculated to be unity by equation (10), indicating no change in AUCiv or Css. This is not the case when the inhibition is so extensive that CLh is not limited by the flow rate any more.

ii. Low clearance drug. In the case of a low clearance drug, CLh = fb · CLint and CLh' = fb. CLint'. If the protein binding is not altered by the inhibitor, the ratio (y) of CLh and CLh' can be calculated as follows:
  <UP>y</UP>=<FR><NU><UP>CL<SUB>h</SUB>′</UP></NU><DE><UP>CL</UP><SUB><UP>h</UP></SUB></DE></FR><UP> = </UP><FR><NU><UP>f<SUB>b</SUB></UP> · <UP>CL<SUB>int</SUB>′</UP></NU><DE><UP>f<SUB>b</SUB></UP> · <UP>CL</UP><SUB><UP>int</UP></SUB></DE></FR>=<FR><NU><UP>f<SUB>b</SUB></UP> · (<UP>CL<SUB>int</SUB>′<SUB>,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,2</UP></SUB>)</NU><DE><UP>f<SUB>b</SUB></UP> · (<UP>CL<SUB>int,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,2</UP></SUB>)</DE></FR> (11)
=<FR><NU><UP>CL<SUB>int</SUB>′<SUB>,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,1</UP></SUB>(<UP>1/f<SUB>m</SUB></UP>−<UP>1</UP>)</NU><DE><UP>CL<SUB>int,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,1</UP></SUB>(<UP>1/f<SUB>m</SUB></UP>−<UP>1</UP>)</DE></FR>
=<FR><NU><UP>CL<SUB>int</SUB>′<SUB>,1</SUB></UP>+<UP>CL<SUB>int,1</SUB>/f<SUB>m</SUB></UP>−<UP>CL</UP><SUB><UP>int,1</UP></SUB></NU><DE><UP>CL<SUB>int,1</SUB>/f</UP><SUB><UP>m</UP></SUB></DE></FR>
=<UP>f<SUB>m</SUB></UP> · <UP>CL<SUB>int</SUB>′<SUB>,1</SUB>/CL<SUB>int,1</SUB></UP>+1−<UP>f<SUB>m</SUB></UP>
Combining equations (10) and (11) yields the following equation:
  <UP>R<SUB>c</SUB></UP>=<FR><NU><UP>1</UP></NU><DE><UP>f</UP><SUB><UP>h</UP></SUB>(<UP>f<SUB>m</SUB></UP> · <UP>CL<SUB>int</SUB>′<SUB>,1</SUB>/CL<SUB>int,1</SUB></UP>+1−<UP>f</UP><SUB><UP>m</UP></SUB>)+1−<UP>f</UP><SUB><UP>h</UP></SUB></DE></FR> (12)
=<FR><NU><UP>1</UP></NU><DE><UP>f<SUB>h</SUB></UP> · <UP>f<SUB>m</SUB></UP> · <UP>CL<SUB>int</SUB>′<SUB>,1</SUB>/CL<SUB>int,1</SUB></UP>+1−<UP>f<SUB>h</SUB></UP> · <UP>f</UP><SUB><UP>m</UP></SUB></DE></FR>
Because Cu,ss encountered clinically is usually much less than Km, CLint,1 and CLint',1 can be expressed as follows:
<UP>CL<SUB>int,1</SUB></UP>=<UP>V<SUB>max</SUB>/K<SUB>m</SUB> and CL<SUB>int</SUB>′<SUB>,1</SUB></UP>=<UP>V<SUB>max</SUB>/K<SUB>m</SUB></UP>(1+<UP>I<SUB>u</SUB>/K<SUB>i</SUB></UP>)
where Iu is the unbound concentration of the inhibitor. Therefore,
<UP>CL<SUB>int</SUB>′<SUB>,1</SUB>/CL<SUB>int,1</SUB></UP>=<FR><NU><UP>1</UP></NU><DE>1+<UP>I<SUB>u</SUB>/K</UP><SUB><UP>i</UP></SUB></DE></FR> (13)
can be derived. Combining equations (12) and (13) yields the following equation:
<UP>R<SUB>c</SUB></UP>=<FR><NU><UP>1</UP></NU><DE><UP>f<SUB>h</SUB></UP> · <UP>f<SUB>m</SUB></UP> · {1/(1+<UP>I<SUB>u</SUB>/K<SUB>i</SUB></UP>)}+1−<UP>f<SUB>h</SUB></UP> · <UP>f</UP><SUB><UP>m</UP></SUB></DE></FR> (14)
It is clear from equation (14) that, in the case of the intravenous administration of a low clearance drug, the degree of increase in AUCiv is determined not by Km or Cu,ss but by Ki, Iu, fh, and fm, if Km >>  Cu,ss.

b. ORAL ADMINISTRATION. The change in AUCpo after a single oral administration and that in Css,av after repeated oral administration can be expressed by the following equation (15), if the dose and administration interval is constant:

  <UP>R<SUB>c</SUB></UP>=<FR><NU><UP>AUC</UP><SUB><UP>po</UP></SUB>(<UP>+inhibitor</UP>)</NU><DE><UP>AUC</UP><SUB><UP>po</UP></SUB>(<UP>−inhibitor</UP>)</DE></FR>=<FR><NU><UP>C</UP><SUB><UP>ss,av</UP></SUB>(<UP>+inhibitor</UP>)</NU><DE><UP>C</UP><SUB><UP>ss,av</UP></SUB>(<UP>−inhibitor</UP>)</DE></FR> (15)
=<FR><NU><UP>CL</UP><SUB><UP>oral</UP></SUB></NU><DE><UP>CL<SUB>oral</SUB>′</UP></DE></FR>=<FR><NU>(<UP>CL<SUB>h</SUB></UP>+<UP>CL</UP><SUB><UP>r</UP></SUB>)<UP>/</UP>(<UP>F<SUB>h</SUB></UP> · <UP>F</UP><SUB><UP>a</UP></SUB>)</NU><DE>(<UP>CL<SUB>h</SUB>′</UP>+<UP>CL</UP><SUB><UP>r</UP></SUB>)<UP>/</UP>(<UP>F<SUB>h</SUB>′</UP> · <UP>F</UP><SUB><UP>a</UP></SUB>)</DE></FR>
=<FR><NU>{<UP>CL<SUB>h</SUB></UP>+<UP>CL</UP><SUB><UP>h</UP></SUB>(<UP>1/f<SUB>h</SUB></UP>−<UP>1</UP>)}<UP>/F</UP><SUB><UP>h</UP></SUB></NU><DE>{<UP>CL<SUB>h</SUB>′</UP>+<UP>CL</UP><SUB><UP>h</UP></SUB>(<UP>1/f<SUB>h</SUB></UP>−<UP>1</UP>)}<UP>/F<SUB>h</SUB>′</UP></DE></FR>
=<FR><NU><UP>CL<SUB>h</SUB>/f</UP><SUB><UP>h</UP></SUB></NU><DE>(<UP>CL<SUB>h</SUB>′</UP>+<UP>CL<SUB>h</SUB>/f<SUB>h</SUB></UP>−<UP>CL</UP><SUB><UP>h</UP></SUB>)</DE></FR> · <FR><NU><UP>F<SUB>h</SUB>′</UP></NU><DE><UP>F</UP><SUB><UP>h</UP></SUB></DE></FR>
=<FR><NU><UP>1</UP></NU><DE><UP>f<SUB>h</SUB></UP> · <UP>CL<SUB>h</SUB>′/CL<SUB>h</SUB></UP>+1−<UP>f</UP><SUB><UP>h</UP></SUB></DE></FR> · <FR><NU><UP>F<SUB>h</SUB>′</UP></NU><DE><UP>F</UP><SUB><UP>h</UP></SUB></DE></FR>
where Fh is hepatic availability and Fa is the fraction absorbed from the gastrointestinal tract into the portal vein.

Some kind of mathematical model has to be introduced for the calculation of the hepatic intrinsic clearance (CLint) in vivo. In order to avoid a false negative prediction of drug-drug interactions, we tried to evaluate the maximum inhibitory effect expected. The well-stirred model was selected as one which can detect the maximum effect of the inhibitor. In the case of oral administration where D is dose, D/AUCpo = D/tau /Css,av = CLh/Fh = fb · CLint can be derived based on the well-stirred model, irrespective of the value of CLh relative to Qh, where D is dose. In this model, therefore, either AUCpo or Css,av is affected directly by the reduction in CLint without a contribution from the hepatic blood flow rate. For this reason, the well-stirred model can detect the maximum effect of an inhibitor. Thus, the well-stirred model was used in the following discussion of the prediction of drug-drug interactions after oral administration.

i. High clearance drug. Because fb · CLint is much larger than the hepatic blood flow rate (Qh <<  fb · CLint), CLh is rate-limited by the flow rate (CLh = Qh). When the altered CLh is still rate-limited by the flow rate (CLh' = Qh), i.e., Qh <<  fb · CLint', then CLh' equals CLh. On the other hand, Fh = Qh/(fb · CLint) and Fh' = Qh/(fb · CLint'). Therefore, the following equation (16) can be derived from equation (15):
   <UP>R<SUB>c</SUB></UP>=<FR><NU><UP>F<SUB>h</SUB>′</UP></NU><DE><UP>F</UP><SUB><UP>h</UP></SUB></DE></FR>=<FR><NU><UP>CL</UP><SUB><UP>int</UP></SUB></NU><DE><UP>CL<SUB>int</SUB>′</UP></DE></FR>=<FR><NU><UP>CL<SUB>int,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,2</UP></SUB></NU><DE><UP>CL<SUB>int</SUB>′<SUB>,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,2</UP></SUB></DE></FR> (16)
=<FR><NU><UP>CL<SUB>int,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,1</UP></SUB>(<UP>1/f<SUB>m</SUB></UP>−<UP>1</UP>)</NU><DE><UP>CL<SUB>int</SUB>′<SUB>,1</SUB></UP>+<UP>CL</UP><SUB><UP>int,1</UP></SUB>(<UP>1/f<SUB>m</SUB></UP>−<UP>1</UP>)</DE></FR>
=<FR><NU><UP>CL<SUB>int,1</SUB>/f</UP><SUB><UP>m</UP></SUB></NU><DE><UP>CL<SUB>int</SUB>′<SUB>,1</SUB></UP>+<UP>CL<SUB>int,1</SUB>/f<SUB>m</SUB></UP>−<UP>CL</UP><SUB><UP>int,1</UP></SUB></DE></FR>
=<FR><NU><UP>1</UP></NU><DE><UP>f<SUB>m</SUB></UP> · <UP>CL<SUB>int</SUB>′<SUB>,1</SUB>/CL<SUB>int,1</SUB></UP>+1−<UP>f</UP><SUB><UP>m</UP></SUB></DE></FR>
Furthermore, as Cu,ss encountered clinically is usually much less than Km, CLint,1 and CLint',1 can be expressed as follows:
<UP>CL<SUB>int,1</SUB></UP>=<UP>V<SUB>max</SUB>/K<SUB>m</SUB> and CL<SUB>int</SUB>′<SUB>,1</SUB></UP>=<UP>V<SUB>max</SUB>/K</UP><SUB><UP>m</UP></SUB>(<UP>1</UP>+<UP>I<SUB>u</SUB>/K</UP><SUB><UP>i</UP></SUB>)
Therefore,
<UP>CL<SUB>int</SUB>′<SUB>,1</SUB>/CL<SUB>int,1</SUB></UP>=<FR><NU><UP>1</UP></NU><DE><UP>1</UP>+<UP>I<SUB>u</SUB>/K</UP><SUB><UP>i</UP></SUB></DE></FR> (17)
can be derived. Combining equations (16) and (17) yields the following equation:
<UP>R<SUB>c</SUB></UP>=<FR><NU><UP>1</UP></NU><DE><UP>f<SUB>m</SUB> ·  </UP>{<UP>1/</UP>(<UP>1</UP>+<UP>I<SUB>u</SUB>/K</UP><SUB><UP>i</UP></SUB>)}+1−<UP>f</UP><SUB><UP>m</UP></SUB></DE></FR> (18)

ii. Low clearance drug. Since the first-pass hepatic availability is close to unity for low clearance drugs, the final equation (14) should be the same for intravenous and oral administration.

The effect of the inhibitor on the Cmax after oral administration also depends on the clearance of the drug. Assuming that the drug absorption from the gastrointestinal tract is sufficiently rapid, Cmax is proportional to Fh. Based on the well-stirred model, Fh can be expressed as follows:
<UP>F<SUB>h</SUB></UP>=<UP>Q<SUB>h</SUB>/</UP>(<UP>Q<SUB>h</SUB></UP>+<UP>f<SUB>b</SUB></UP> · <UP>CL</UP><SUB><UP>int</UP></SUB>) (19)
It is clear from equation (19) that Fh is minimally affected by the change in CLint in the case of a low clearance drug (Qh >>  fb · CLint: Fh = 1), but is inversely proportional to CLint in the case of a high clearance drug (Qh <<  fb · CLint: Fh = Qh/fb · CLint), in which case Cmax also changes in inverse proportion to CLint.

In summary, it is important to know the values of Ki, Iu, fh, and fm in order to predict in vivo drug-drug interactions. Approximated fh and fm values can be estimated from the urinary recovery of the parent compound and each metabolite. Ki values can be evaluated by kinetic analysis of in vitro data using human liver microsomes and recombinant systems and this has already been done for many compounds. The key, therefore, is the evaluation of Iu.

2. The evaluation of the unbound concentration of the inhibitor in vivo. Although Iu is the unbound concentration of the inhibitor around the metabolic enzyme in the liver, it is impossible to directly measure this in vivo. However, many drugs are transported into the liver by passive diffusion, allowing for the assumption that the unbound concentration in the liver equals that in the liver capillary at steady-state. This means that estimating the unbound concentration of the inhibitor in the liver capillary may be enough for some drugs. This assumption is not valid, however, in the case of drugs which are actively transported into or out of the liver; the unbound concentration in the liver may be higher in the former case or lower in the latter than in the liver capillary (fig. 3). In these cases, another experiment using human hepatocytes, human liver slices, etc., is required to estimate the kinetic parameters for the active transport. Furthermore, the unbound concentration in the liver capillary is always changing and a concentration gradient is formed from the entrance (portal vein) to the exit (hepatic vein). Which of these concentrations should be considered as Iu? An underestimation of Iu may lead to a "false negative" prediction of actually occurring in vivo drug interaction from in vitro data. In order to avoid a false negative prediction caused by underestimation of Iu, the plasma unbound concentration at the entrance to the liver, where the blood flow from the hepatic artery and portal vein meet, was considered the maximum value of Iu and was used in the prediction (Iin,u; fig. 4).


View larger version (18K):
[in this window]
[in a new window]
 
Fig. 3.   Relationship between unbound drug concentration in the liver capillary (Cf,p) and that in the liver (Cf,T). Cin and Cout represent the drug concentration at the entrance (portal vein side) and the exit (hepatic vein side) of the liver, respectively.


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 4.   Model for estimating inflow concentration of the inhibitor into the liver after oral administration (Iin). Iout, I, and Iblood represent the inhibitor concentration at the exit of liver (hepatic vein side), the inhibitor concentration at the liver capillary, and inhibitor concentration in the systemic circulation, respectively. Qa, Qpv and Qh (=Qa + Qpv) represent blood flow at hepatic artery, portal vein, and hepatic vein, respectively. ka, D, and Fa represent the absorption rate constant, dose, and the fraction absorbed from the gastrointestinal tract into the portal vein, respectively, of the inhibitor.

Practically, the maximum plasma concentration in the circulation (Imax) has been estimated for many inhibitors. When the value of Imax is not reported, it can be predicted from both the plasma concentration at a single time point after administration and the pharmacokinetic parameters such as the elimination half life (t1/2,beta ).

According to the model in figure 4, influx into the liver consists of contributions from the hepatic artery and portal vein (after gastrointestinal absorption). When the drug is absorbed from the gastrointestinal tract with a first-order rate constant (ka), the maximum influx rate into the liver (vin,max) can be expressed as follows:
<UP>v<SUB>in,max</SUB> ≦ Q<SUB>a</SUB>I<SUB>max</SUB></UP>+<UP>Q<SUB>pv</SUB>I<SUB>max</SUB></UP>+<UP>k<SUB>a</SUB></UP> · <UP>D · F<SUB>a</SUB></UP> · <UP>e<SUP>−kat′</SUP></UP> (20)
where Qa and Qpv represent the blood flow rate in the hepatic artery and the portal vein, respectively, Fa is the fraction absorbed from the gastrointestinal tract into the portal vein, and t' is the time after oral administration (after subtraction of the lag-time).

When the absorption rate is maximum (i.e., t' = 0), the final term in equation (20) can be expressed as ka · D · Fa and thus:
<UP>v<SUB>in,max</SUB> ≦ </UP>(<UP>Q<SUB>a</SUB></UP>+<UP>Q</UP><SUB><UP>pv</UP></SUB>)<UP>I<SUB>max</SUB></UP>+<UP>k<SUB>a</SUB></UP> · <UP>D</UP> · <UP>F<SUB>a</SUB></UP> (21)
As Qh = Qa + Qpv, the following equation can be derived:
<AR><R><C><UP>I<SUB>in,max</SUB></UP> = <UP>v<SUB>in,max</SUB>/Q<SUB>h</SUB>=</UP></C></R><R><C> </C></R><R><C> </C></R><R><C> </C></R><R><C> </C></R><R><C> </C></R></AR><AR><R><C>  <UNL><UP>I</UP><SUB><UP>max</UP></SUB></UNL></C></R><R><C> </C></R><R><C><UP>Contribution </UP></C></R><R><C><UP>from the</UP></C></R><R><C><UP>systemic </UP></C></R><R><C><UP>circulation</UP></C></R></AR><AR><R><C>+ </C></R><R><C> </C></R><R><C> </C></R><R><C> </C></R><R><C> </C></R><R><C> </C></R></AR><AR><R><C><UNL>(<UP>k<SUB>a</SUB></UP> · <UP>D/Q</UP><SUB><UP>h</UP></SUB>) · <UP>F</UP><SUB><UP>a</UP></SUB></UNL></C></R><R><C> </C></R><R><C><UP>Contribution </UP></C></R><R><C><UP>from the</UP></C></R><R><C><UP>absorption</UP></C></R><R><C> </C></R></AR> (22)
Therefore, Iin,max can be predicted if the parameters such as ka and Fa are available for the inhibitor. Taking the plasma protein binding into consideration, the unbound Iin,max (Iin,max,u) can be calculated as Iin,max · fu. Finally, comparing the value of Iin,max,u as Iin,u and that of Ki obtained in vitro allows the prediction of the maximum degree of in vivo drug-drug interaction caused by metabolic inhibition.

In general, the apparent absorption rate of the orally administered drug is maximum when the gastrointestinal absorption of the drug is so rapid that the rate-limiting step is the gastric emptying rate. A first order rate constant (ka) of about 0.1 min-1 is reported for gastric emptying in rats and humans (Oberle et al., 1990). On the other hand, the absorption rate constant in humans can be calculated from the time to reach the maximum concentration (Tmax) and the elimination constant (kel) as follows:
<UP>T<SUB>max</SUB></UP>=ln(<UP>k<SUB>a</SUB>/k<SUB>el</SUB></UP>)/(<UP>k<SUB>a</SUB></UP>−<UP>k<SUB>el</SUB></UP>) (23)
In practice, however, because of the possible estimation error of Tmax caused by interindividual differences etc., the calculated value of ka sometimes exceeds 0.1 min-1, though it should never exceed that, theoretically, for gastric emptying. Therefore, the theoretically maximum value of 0.1 min-1 was used for ka when it was calculated to be larger than 0.1 min-1. Moreover, in order to avoid a false negative prediction, the maximum ka of 0.1 min-1 was used to obtain the largest inhibitor concentration if ka was unknown.

E. Examples of the Prediction of Drug-Drug Interactions Based on Literature Data

The methodology described above (see Section III.D.) has been applied to the prediction of in vivo drug-drug interactions from in vitro data gathered from the literature.

1. Successful Cases of In Vitro/In Vivo Prediction. a. TOLBUTAMIDE-SULFAPHENAZOLE. Interactions between tolbutamide and sulfa-agents cause serious side effects such as hypoglycemic shock in patients (Christensen et al., 1963) and exhibit the marked interspecies differences in animals. Veronese et al. (1990) reported about a five-fold increase in both AUCpo and t1/2 of tolbutamide in humans following coadministration of 500 mg sulfaphenazole (table 2, fig. 5).

                              
View this table:
[in this window]
[in a new window]
 
TABLE 2
Inhibition of tolbutamide metabolism (CYP2C9) by coadministration of sulfaphenazole


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 5.   Effect of sulfaphenazole coadministration on plasma concentration of tolbutamide (Veronese et al., 1990). triangle : Tolbutamide (500 mg p.o.) alone; black-triangle: Tolbutamide (500 mg p.o.)+Sulfaphenazole (500 mg p.o., q12h).

The t1/2 of intravenous tolbutamide is prolonged and the CLtot is reduced markedly in rats, too, by sulfaphenazole (Sugita et al., 1981). On the contrary, the CLtot of tolbutamide is increased 15 to 30% in rabbits with little change in the t1/2 (fig. 6) (Sugita et al., 1984). Because tolbutamide is a low clearance drug with a small urinary excretion, the CLtot after intravenous administration is expressed by the following equation (24):
<UP>CL<SUB>tot</SUB></UP>=<UP>D/AUC</UP>=<UP>f<SUB>b</SUB></UP> · <UP>CL<SUB>int</SUB></UP> (24)
Sugita et al. (1984) tried to reconstruct the CLtot in vivo based on the values of unbound fraction in blood (fb) and CLint estimated separately by in vitro binding and metabolic studies. Sulfaphenazole inhibits both plasma protein binding and hepatic metabolism of tolbutamide, causing the increase in fb and the reduction in CLint, in both species. Although the CLint falls to about one-fourth and the fb increases about two-fold in rats, resulting in about a half-fold reduction in the CLtot, the CLint falls to about one-half and the fb increases about two-fold in rabbits resulting in little change in the CLtot (fig. 7). The effects on the CLint and fb of tolbutamide are not constant among sulfa-agents; sulfadimethoxine also reduces the CLint by about one-half in rabbits but increases the fb more than two-fold, resulting in an increase in the CLtot and a reduction in the AUC (Sugita et al., 1984).


View larger version (15K):
[in this window]
[in a new window]
 
Fig. 6.   Effect of sulfaphenazole (SP) coadministration on plasma concentration of tolbutamide (TB) in rabbits (A) and rats (B) (Sugita et al., 1981, 1984). Open and closed symbols represent plasma concentrations of sulfaphenazole (or its metabolite, N4-Ac SP) and tolbutamide, respectively.


View larger version (31K):
[in this window]
[in a new window]
 
Fig. 7.   Prediction of interspecies difference in the interaction of tolbutamide and sulfaphenazole (SP) or sulfadimethoxine (SDM) from in vitro data (Sugita et al., 1984).

The interaction between tolbutamide and sulfaphenazole involves both plasma protein binding and hepatic metabolism in humans, too. The main metabolic pathway of tolbutamide in vitro is CYP2C9-mediated hydroxylation, and the metabolite undergoes sequential metabolism forming a carboxylate in vivo (Thomas and Ikeda, 1966; Nelson and O'Reilly, 1961). The Ki of sulfaphenazole, a specific inhibitor of CYP2C9, for tolbutamide methyl hydroxylation in human liver microsomes in vitro is 0.1-0.2 µM (Miners et al., 1988; Back et al., 1988). The Imax of sulfaphenazole after a 500 mg dose was 70 µM in humans, and the absorption term [the second term in equation (22)] was calculated to be 8.0 µM using ka = 0.0095 min-1, D = 500 mg, Qh = 1610 ml/min, and Fa = 0.85. Iin,max was, therefore, calculated to be 78 µM, indicating that the contribution of systemic circulation is greater than that of absorption. Taking the fu value (0.32) of sulfaphenazole into consideration, Iin,u/Ki was calculated to be 125-250 (table 2). The plasma protein binding of tolbutamide is also inhibited by sulfaphenazole in humans, resulting in about a three-fold increase in fb (Christensen et al., 1963). However, the inhibition is considered almost complete in terms of the product of fb and CLint because the extent of inhibition of metabolism is much greater than that of its plasma protein binding. The contribution of the CYP2C9-related metabolic pathway of tolbutamide is about 80% of the total (fh · fm = CLh,m/CLtot = 0.8) (table 2). Therefore, complete inhibition of this metabolic pathway leads to an 80% reduction in CLint, and the AUCpo is predicted to be five times larger than the control value, which is consistent with the observed increase (table 2).

b. TRIAZOLAM-KETOCONAZOLE. Von Moltke et al. (1996) reported that plasma triazolam concentration after oral administration of 0.125 mg was greatly elevated by oral ketoconazole (200 mg), producing a nearly nine-fold reduction in the apparent oral clearance. They predicted this interaction based on in vitro studies using human liver microsomes (table 3). Triazolam is eliminated in humans mainly by CYP3A-mediated metabolism to alpha -hydroxy (OH)- and 4-OH-triazolam. Ketoconazole is a powerful inhibitor of both these metabolic pathways, with a mean Ki value of 0.006 and 0.023 µM, respectively (Von Moltke et al., 1996). In order to estimate ketoconazole concentrations in the liver, they conducted an in vitro study using mouse liver homogenates in human plasma spiked with ketoconazole; a liver/plasma partition ratio of 1.12 was obtained. On the other hand, the partition ratio was calculated to be 2.03 in the in vivo mouse study where the ketoconazole concentrations in plasma and liver were measured. The concentration of ketoconazole in the liver was estimated by multiplying this partition ratio by the total ketoconazole concentration in plasma in the clinical study (0.04-9.32 µM). Using the in vitro Ki values, ketoconazole concentration in the liver, and the contribution of both metabolic pathways (52.5% and 47.5% for alpha - and 4-OH-pathway, respectively), the predicted degree of reduction (>95%) in triazolam clearance in vivo was consistent with the 88% reduction actually observed in vivo (table 3). However, it should be noted that in this report, the total concentration of the inhibitor was used instead of unbound concentration in the liver. The unbound concentration needs to be estimated because the Ki values obtained in the in vitro studies are based on the concentration in the medium.

                              
View this table:
[in this window]
[in a new window]
 
TABLE 3
Comparison of the prediction of triazolam-ketoconazole interaction by von Moltke's and our method

Using our method proposed above (see Section III.D.), the ketoconazole-triazolam interaction would be predicted as follows: The Imax,ss of ketoconazole during administration of 200 mg×2/day was 6.6 µM (Daneshmend et al., 1981), and the absorption term [the second term in equation (22)] was calculated to be 1.4-2.5 µM using ka = 0.0099-0.018 min-1, D = 200 mg, Qh=1610 ml/min, and Fa = 0.59. The ka was calculated from equation (23) using the values of Tmax and t1/2 (=0.693/kel) (Daneshmend et al., 1984). The Iin,max is, therefore, calculated to be 8.0-9.1 µM. Since the fu of ketoconazole is 0.01, the Iin,max,u is calculated to be 0.080-0.091 µM and the obtained Iin,max,u/Ki value is 13-15 and 3.5-4.0 for the alpha -OH and 4-OH pathways, respectively, using a Ki value of 0.006 and 0.023 µM, respectively. Therefore, the reduction in the clearance can be estimated as follows, considering the contribution of each pathway to the total metabolism:
<UP>R</UP>=<FR><NU><UP>1</UP></NU><DE><UP>1</UP>+(<UP>13–15</UP>)</DE></FR>×<UP>0.525</UP>+<FR><NU><UP>1</UP></NU><DE><UP>1</UP>+(<UP>3.5–4.0</UP>)</DE></FR>×<UP>0.475</UP>
=<UP>0.128–0.143</UP>
Thus, an 85.7-87.2% reduction is predicted by this method, which is very close to the observed reduction (88%) (table 3). The degree of the inhibition should be larger if ketoconazole is actively transported into the liver.

Two cases have been shown here in which the interaction that had actually occurred in vivo was successfully predicted based on in vitro metabolism data. On the other hand, we believe that the ability to predict by the above-mentioned methods based on in vitro data should be very high in the case of drug combinations that do not interact with each other in vivo. In other words, the absence of in vivo drug-drug interactions should be successfully predicted, which has been partly confirmed in our study, though the data are not shown here.

2. Interactions Predictable for the Objective Metabolic Pathway but not Predictable for the Overall Data. a. SPARTEINE-QUINIDINE. Schellens et al. (1991) reported that the CLoral of sparteine (dose: 50 mg) fell from 979 to 341 ml/min (35% of the control value) after coadministration of 200 mg quinidine (table 4). The main metabolic pathway of sparteine is CYP2D6-mediated dehydration. Because quinidine is a specific inhibitor of CYP2D6, it is reasonable that metabolic inhibition is involved in this quinidine-induced reduction in the CLoral of sparteine. The Ki of quinidine for the CYP2D6-mediated metabolism in human liver microsomes in vitro is reported to be 0.06 µM. The Imax of quinidine after a dose of 200 mg was 4.1 µM, and the absorption term [the second term in equation (22)] was calculated to be 0.86-22 µM using ka = 0.0027-0.069 min-1, D = 200 mg, Qh = 1610 ml/min, and Fa = 0.83. Iin,max is, therefore, calculated to be 5-26 µM. Because the fu of quinidine is 0.15, the Iin,u and Iin,u/Ki are calculated to be 0.75-3.9 µM and 13-65, respectively (table 4). Thus, it was predicted that the dehydration pathway of sparteine would be almost completely inhibited by quinidine. The contribution of the dehydration pathway of sparteine to the total elimination is about 25% (fh · fm = CLh,m/CLtot = 0.25) (table 4). Therefore, the complete inhibition of this dehydration pathway will reduce the CLoral to 75% of the control value, which is about two-fold larger than the observed reduction to 35%. The reasons for this discrepancy may include the possibility that metabolic pathways other than dehydration may also be inhibited by quinidine.

                              
View this table:
[in this window]
[in a new window]
 
TABLE 4
Inhibition of sparteine metabolism (CYP2D6) by coadministration of quinidine

b. TERFENADINE-KETOCONAZOLE. Honig et al. (1993) reported that blood concentrations of terfenadine (dose: 60 mg), detectable only in one of six subjects when administered alone (Imax = 7 ng/ml), became detectable in all subjects following coadministration of 200 mg ketoconazole, with the highest Imax in the same subject elevated to 81 ng/ml (table 5). The main elimination pathway of terfenadine is CYP3A4-mediated N-dealkylation and hydroxylation yielding a carboxylate, which may undergo sequential N-dealkylation (Garteiz et al., 1982). Ki values of ketoconazole for terfenadine metabolism in vitro have been reported by two groups. Jurima-Romet et al. (1994) reported the Ki values of 3 and 10 µM in human liver microsomes and <1 and 3 µM in human hepatocytes (Jurima-Romet et al., 1996; Li et al., 1997) for the N-dealkylation and hydroxylation pathways, respectively. Based on the Ki values in human liver microsomes reported by Von Moltke et al. (1994), Iin,u/Ki value was calculated to be 3.3-3.8 and 0.33-0.38 for the N-dealkylation and hydroxylation pathways, respectively (table 5) [see our previous review (Ito et al., 1998) for the details]. As the estimated contribution (fh · fm = CLh,m/CLtot) of these two pathways to the total metabolism was about 0.13 and 0.45, respectively, the increase in the availability and Cmax caused by the metabolic inhibition was predicted to be about 1.3-fold according to equation (14). However, the Cmax was actually increased more than 10-fold. The possible reasons to explain this great discrepancy include: inaccurate measurements of clinical concentrations of terfenadine being around the detection limit of the assay, and the contribution of the other 50% of the metabolism of terfenadine neglected in the analysis. As shown later (see Section III.H.2.), interactions involving metabolism and/or efflux process in the gut may have some contributions in the in vitro/in vivo discrepancy. If the value of Iin,u/Ki</