Department of Pharmaceutics, Faculty of Pharmaceutical
Sciences, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan.
 |
I. Introduction |
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 |
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:
|
(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 |
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
).
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:
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):
|
(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:
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):
|
(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:
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):
|
(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
):
|
(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:
|
(6)
|
|
(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:
|
(8)
|
Equation (7) can be rearranged to give the following equation:
|
(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:
|
(10)
|
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:
|
(11)
|
Combining equations (10) and (11) yields the following equation:
|
(12)
|
Because Cu,ss encountered clinically is usually much
less than Km, CLint,1 and CLint',1
can be expressed as follows:
where Iu is the unbound concentration of the
inhibitor. Therefore,
|
(13)
|
can be derived. Combining equations (12) and (13) yields the
following equation:
|
(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:
|
(15)
|
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/
/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):
|
(16)
|
Furthermore, as Cu,ss encountered clinically is
usually much less than Km, CLint,1 and
CLint',1 can be expressed as follows:
Therefore,
|
(17)
|
can be derived. Combining equations (16) and (17) yields the
following equation:
|
(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:
|
(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,
).
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:
|
(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:
|
(21)
|
As
Qh = Qa + Qpv,
the following equation can be derived:
|
(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:
|
(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 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 ).
: Tolbutamide (500 mg p.o.) alone; : 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):
|
(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
-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
- 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.
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
-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:
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.
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