Abstract
Models of the time course of the effect of P450 induction on substrate clearance have previously only considered induction through enhanced synthesis of protein. Induction of CYP2E1 does not always conform to this model, in that many chemicals induce the enzyme through stabilization of the protein apparently by binding to the active site. While such binding protects the enzyme from degradation, it also results in competitive inhibition of substrate clearance. We present a model based on experimental studies of chemical induction of CYP2E1 by ligand stabilization through which this mechanism of induction can be translated into its pharmacokinetic consequence with regard to clearance of substrate. CYP2E1 is considered to be localized in two pools: Pool 1 at which two mechanisms of degradation, fast and slow, operate and pool 2, at which only the slower mechanism operates. Binding of substrate to enzyme in pool 1 stabilizes it from degradation by the fast process, leaving only the slow process. Ligand stabilization therefore results in induction of CYP2E1 as enzyme accumulates as a consequence of unchanged synthesis. Binding of ligand to the active site results in competitive inhibition of the clearance of substrate. Model-based computer simulations show that the time course of interaction between inhibitor/inducer and substrate can be predicted from knowledge of I/Ki and S/Km and the synthesis and degradation kinetics of CYP2E1. The simulations demonstrate further that as long as inhibitor/inducer administration is not interrupted, the clearance of substrate will always be less than the value observed at low concentration of substrate even if the substrate concentration is raised to displace inhibitor/inducer from the active site. On the other hand, the degree of inhibition of clearance is less than would be seen if induction had not taken place. Clearance of substrate is observed to rise above the value observed in the absence of the inhibitor/inducer only after the inhibitor/inducer concentration declines low enough for substrate to gain access to the active site of the enzyme. The model-based simulations agree with reports of the interaction between isoniazid and acetaminophen in humans.
Cytochrome P450 2E1 (CYP2E1)1 is a constitutive and inducible mammalian form of the super family of cytochrome P450 enzymes (1-2). Physiologic and pathologic conditions associated with the induction of CYP2E1 include birth (3-5), fasting (6), obesity (7-8), and diabetes (9-11). These conditions modulate the level of mRNA via transcriptional activation, post-transcriptional activation, message stabilization, or some combination of these mechanisms (12-16). The biologic consequence of each of these mechanisms of induction is an increase in thede novo synthesis of CYP2E1 protein.
The pharmacokinetic consequences of induction via increasedde novo synthesis are well understood, as classical inducers such as 3-methylcholanthrene and phenobarbital fit into this group (17-20). This general mechanism results in a first-order rise in functional enzyme and substrate clearance to a new stable plateau that is maintained as long as the administration of the inducer is not altered (21-23).
Induction by Stabilization.
Chemical inducers of CYP2E1 in animals include isoniazid, imidazoles, pyridine, and other small molecular weight alcohols, as well as substances produced endogenously such as acetone (1-2, 15, 24-25). Chemical induction of CYP2E1 occurs either by increasing de novo synthesis or by decreasing degradation of the protein. Stabilization of CYP2E1 protein is of particular interest because the consequences of this mechanism with regard to the metabolism of substrates in vivo have only recently been systematically probed (26-30). The phenomenon of CYP2E1 stabilization was first shown with acetone in rats (14). In vivo labeling of protein and subsequent immunopurification demonstrated that native CYP2E1 protein turnover followed a biphasic pattern with respective half-lives of 7 and 37 hr. Treatment with acetone abolished the rapid component such that only the longer half-life was observed. Separate labeling data suggested that there was an apparent decrease in the de novosynthesis rate of CYP2E1 with acetone administration. This result indicated that induction of CYP2E1 by acetone is primarily due to stabilization of CYP2E1.
CYP2E1 induction through stabilization by ethanol has been demonstrated in rats at peak plasma ethanol concentrations ranging from 1.0 to 2.5 mg/ml (5, 31-34). In this concentration range, CYP2E1 activity and immunodetectable CYP2E1 in liver microsomes from ethanol-induced rats was increased in the absence of an increase in specific mRNA. In agreement with the results of CYP2E1 induction by acetone, CYP2E1 protein turnover followed a biphasic pattern, and ethanol treatment abolished the rapid CYP2E1 turnover phase (34-35).
Stabilization of CYP2E1 has also been demonstrated in primary cultures of rat hepatocytes with several small molecular weight ligands. A correlation between the binding affinities of ligands to the active site of CYP2E1 and their ability to prevent the loss of the enzyme has been demonstrated (36). CYP2E1 protein was maintained by ligand even when mRNA rapidly disappeared, showing that the maintenance of CYP2E1 in primary hepatocyte culture was independent of enzyme synthesis.
Taken together, the results of these studies indicate that binding of a ligand to the active site of CYP2E1 protects the protein from degradation, allowing it to accumulate as a consequence of unaltered synthesis.
Pharmacokinetic Consequence in Animals.
The consequence of induction by ligand stabilization with regard to drug metabolism was first illustrated in a study of the effect of chronic ethanol administration on the defluorination of enflurane (37), a substrate of CYP2E1 (38-39). In vivo ethanol administration increased defluorination ex vivo in microsomes prepared from rats in which the plasma concentration of ethanol was ≥ 0.6 mg/ml, as expected after induction of CYP2E1. However, at the same ethanol concentration, in vivodefluorination was inhibited. The in vivo enflurane defluorination increased above the pre-ethanol rate only when plasma ethanol concentration decreased to ≤ 0.5 mg/ml. These results are consistent with the ligand-stabilized mechanism of induction in that inhibition was observed in vivo at high concentrations of ethanol. Enhanced activity was observed in microsomes prepared at the same time since ethanol would be removed during preparation of the microsomes. In vivo defluorination of enflurane was enhanced (revealing the consequence of induction) as ethanol was eliminated and CYP2E1 became available for binding of substrate.
Similar results have been demonstrated with another CYP2E1 substrate, ethyl carbamate, in rats and mice induced with pyridine (40). CYP2E1 catalyzes the bioactivation of ethyl carbamate to a reactive epoxide and administration of pyridine inhibited the metabolism of [14C]-ethyl carbamate to its [14C]-epoxide in vivo. However, formation of this metabolite was enhanced in microsomes prepared at the same time after pyridine administration.
Pharmacokinetic Consequence in Humans.
The pharmacokinetic consequences of CYP2E1 induction have been studied in humans using isoniazid as the inducer and formation of the hepatotoxin N-acetyl-p-benzoquinone imine (NAPQI) from acetaminophen (41, 42) or chlorzoxazone 6-hydroxylation as a probe of CYP2E1 activity (43). That isoniazid stabilizes CYP2E1 has been suggested by an increase in CYP2E1 protein in microsomes of isoniazid-treated animals in which there was no increase in specific mRNA (16, 33).
In the first two of a series of studies (26-30), acetaminophen pharmacokinetics were determined on three occasions in healthy subjects administered isoniazid daily for several days: before isoniazid administration was initiated, together with isoniazid, and 48 hr after the last dose of isoniazid. When acetaminophen was ingested together with isoniazid, a 66% inhibition was observed in the oxidative metabolism of acetaminophen to NAPQI. At 48 hr after the withdrawal of isoniazid, NAPQI formation had returned to the pre-isoniazid value (26). In the second study (27), NAPQI formation was inhibited by 66% when acetaminophen was taken together with isoniazid. At 24 hr after the last dose of isoniazid, a 56% increase in the formation of NAPQI was observed, but it had returned to the pre-induced level 72 hr after the discontinuation of isoniazid. In this study the effect of isoniazid on chlorzoxazone clearance was also examined. When chlorzoxazone was taken together with isoniazid, chlorzoxazone clearance was inhibited by 58% in comparison to the pre-isoniazid control. At 48 hr after the last dose of isoniazid, there was a 50% increase in chlorzoxazone clearance, but this had returned to the pre-isoniazid level 96 hr after the last isoniazid dose.
The influence of isoniazid acetylation status on the interaction between isoniazid and acetaminophen in patients receiving a morning dose of isoniazid for a period of at least 6 months for prophylaxis of tuberculosis was consistent with expectations based on induction of CYP2E1 by stabilization (28). Acetaminophen was administered on separate occasions, either with the morning dose of isoniazid or in the evening 12 hr after the morning dose of isoniazid. At the conclusion of isoniazid therapy, a time-matched control phase (morningversus evening) of acetaminophen administration was performed. When acetaminophen and isoniazid were co-administered, the mean NAPQI formation clearance was inhibited by 64% in slow acetylators and by 52% in fast acetylators. When acetaminophen was given 12 hr after isoniazid, the mean NAPQI formation clearance in fast acetylators was increased 37% from the value observed in the absence of isoniazid therapy, but remained inhibited by 32% in slow acetylators, as expected if isoniazid was still present in those individuals in whom it would be eliminated more slowly.
A similar study was conducted in which chlorzoxazone was used as a probe of CYP2E1 activity (30). Co-administration of chlorzoxazone with isoniazid resulted in 80% inhibition of chlorzoxazone clearance in slow acetylators and 60% in fast acetylators. At 48 hr after the last dose of isoniazid, 6-hydroxy chlorzoxazone formation in slow acetylators was still increased by 60% and had returned to the pre-isoniazid value in fast acetylators.
Isoniazid inhibition and enhancement of NAPQI formation and chlorzoxazone clearance in humans are consistent with induction of CYP2E1 via ligand-stabilization, a mechanism likely to result in persistent inhibition of activity while the inhibitor/inducer is present and enhancement of activity after the inducer is eliminated. The time of reversal from inhibition to enhanced activity depends on the ability of the substrate to compete for the active site as the inhibitor/inducer ligand is eliminated. This mechanism is illustrated in fig. 1.
We have developed a pharmacokinetic model to describe the extent and consequences of enzyme inhibition/induction by ligand-stabilizationin vivo as a function of the affinities and concentrations of the inhibiting/inducing ligand and the substrate for the active site. Simulations conducted based on this model are consistent with the time course of interactions observed in humans and animals.
Theoretical
Model Characteristics.
There are conflicting putative mechanisms by which CYP2E1 is degraded. An ATP/ubiquitin-mediated proteolytic degradation in animals has been demonstrated by two groups of investigators (34, 44). Others suggest that the rapid phase of degradation is accounted for by a Mg/ATP-dependent phosphorylation, while an autophagosomal/autolysosomal process accounts for the slower phase (45-47). Although the biochemical mechanism(s) by which CYP2E1 is degraded are controversial, there is consensus from pulse labeling and immunopurification studies on the biphasic nature of its degradation (14, 34-35).
Several possible models that might result in biphasic loss of CYP2E1 can be proposed. The simplest model consists of a single pool of enzyme in which an endogenous ligand is assumed to slow the otherwise rapid degradation of CYP2E1. However, unless the ligand has a very slow dissociation constant (such that the half-life of dissociation would be several hours), thus making dissociation of the ligand the rate limiting step in the degradation of protein, the kinetics of CYP2E1 in this model are monophasic with a half-life that is a function of the dissociation constant and the elimination rate constants. Accordingly, an appropriate model must consider two kinetically separate pools of enzyme.
Kinetically distinct pools are most easily envisioned as chemically or physically distinct. Chemically distinct forms of the enzyme might be produced by post-translational modifications that do not compromise activity. While such a model would result in appropriate CYP2E1 degradation kinetics, there is no evidence for a post-translational modification of CYP2E1 that maintains catalytic activity. Physically distinct pools of catalytically competent CYP2E1 may exist. Catalytically competent CYP2E1 is present in the endoplasmic reticulum and the Golgi apparatus although the Golgi compartment accounts for relatively little enzyme (48). It is most likely that CYP2E1 in the Golgi arises from the smooth endoplasmic reticulum, while new CYP2E1 appears first in the rough endoplasmic reticulum. It is conceivable that the fast process of degradation operates at only one of these sites, while both processes can operate at the other. Thus, it is reasonable to consider a model in which the enzyme in the rapidly turned-over pool can be stabilized by ligand and thereby be protected from degradation by the fast process while it is still degraded by the slow process.
Several classes of model with physically distinct pools of catalytically competent enzyme can be envisioned. The first allows insertion of new enzyme into both pools without any transfer of enzyme between the pools. One pool has a rapid rate of degradation, which is slowed when ligand binds, while the other has only a slow rate. A second class that produces biphasic kinetics of CYP2E1 loss with physically distinct pools is characterized by insertion of newly synthesized protein into either the rapidly turned over pool or both of the pools and translocation of the protein in either or both directions between the compartments. If there is no translocation from the rapid to the slow pool, stabilization of the protein in the rapidly degraded pool will result in induction of the protein only within that pool. If translocation occurs from the rapidly turned over pool to the more slowly turned over pool, induction by stabilization in the rapidly turned over pool will generally result in an increase of enzyme in the more slowly turned over pool as well. Otherwise, the kinetic characteristics of these models are similar.
It is not clear which model is biochemically correct. However, initial simulations showed that to obtain induction to the extent observed in experimental animals, translocation from the rapidly turned over pool to the slowly turned over pool was required. Also, synthesis into one pool (viz. the rough endoplasmic reticulum pool) seemed most appropriate. Since the kinetic characteristics of the two pool class of models are similar, we present equations describing a model in which new enzyme is added to pool 1 only (viz. rough endoposamic reticulum). Enzyme in pool 1 is subject to a rapid degradation process when unbound by ligand and a slow degradation process when bound. The model is illustrated in fig. 2. Either free or bound enzyme can be translocated to pool 2, viz. smooth endoplasmic reticulum, as a consequence of either loss of ribosomes or by contact between the rough and smooth endoplasmic reticulum (48-52). Enzyme in pool 2 is degraded only by the slow process.
Mathematical model.
The rate of change in tissue content of total CYP2E1 prior to induction by ligand stabilization is described as the sum of two rate equations describing the respective pools of enzyme:
The rates can be written as the product of the respective rate constant and the concentration of the appropriate pool of the enzyme.
Upon induction with acetone or ethanol in vivo, the ligand-stabilized CYP2E1 is turned over at the same rate as the slowly turned over pool (pool 2) of enzyme. In contrast, the binding of the ligand to CYP2E1 in the slow pool does not alter the kinetics of degradation in this pool, although it can bind to CYP2E1 in pool 2. The rate equations describing the model when inhibiting/inducing ligand is introduced are:
The bound and free enzyme fractions are functions of the inducer concentration, I, and the equilibrium dissociation constant, Ki, of the (competitive2) inhibitor/inducer.
Clearance.
The basal clearance of a specific substrate in the absence of the inducer, Cl0, is described by the single-enzyme Michaelis-Menten equation:
An alternative case (“case 2”) can be constructed in which the basal clearance (Cl0) is the clearance in the absence of the inhibitor/inducer and when S ≪Km. In this case, Cl(t) is compared with clearance observed in the absence of inhibitor/inducer andat vanishingly low substrate concentration, regardless of the concentration of substrate influencing the value of Cl(t). Thus, the basal clearance, Cl0,(S/Km→0), is defined as E0·kcat/Km. The ratio of Cl(t) described in eq. 4 to this basal clearance thus becomes:
Enzyme at Steady-State.
As the enzyme concentration reaches steady-state at a steady-state inhibitor/inducer ligand concentration, dEss,total/dt = 0. The total amount of enzyme is the sum of eqs. 1 and 2 under the condition of steady-state:
Clearance at Steady-State.
An expression for clearance at steady-state and amount of enzyme is obtained by substituting eq. 7 into eq. 5 where the pre-induced basal clearance is Cl0 = E0·kcat/(Km+S).
Equation 8On the other hand, substituting eq. 7 into eq. 6, where the basal clearance is Cl0,(S/Km)→0 = E0·kcat/Km, yields an equation that describes the change in clearance relative to clearance at low substrate concentration:
Computer Simulations
Model Implementation.
In the time course simulations that follow, the inhibitor/inducer ligand concentration was simulated as a constant-rate infusion, followed by a first-order rate of elimination after the infusion was discontinued:
The initial division of CYP2E1 between the pools was estimated by simulating the time course of total enzyme with varying fractions in the pools. Since approximately 50% of the total area under the curve in the rat is described by pool 1 (14, 34), the appropriate basal division of total enzyme between the pools was the fraction that satisfied this condition: 85% in pool 1 and 15% pool 2. The degradation rate constants correspond to half-lives of 7 and 37 hr in the rat. ktrans was calculated from the steady-state pre-inhibitor/inducer condition for pool 2, ktrans = kslow·E0,2/E0,1, and corresponded to a half-life of 196 hr.
Eqs. 1, 2, 5, and 6 were used to simulate time course of enzyme level and substrate clearance. Eq. 7 was used to simulate steady-state enzyme level as a function of I/Ki and S/Km. Eqs. 8 through 11 were used to simulate the steady-state conditions of inhibitor/inducer, enzyme level, substrate clearance and inhibited clearance (Cl/E) as functions of I/Ki and S/Km. Simulations were performed using Mathmatica version 2.1.0.2 (Wolfram Research, Inc. Champaign, IL) implemented on a Macintosh computer.
Model Validation
The model was assessed by evaluating key expected results in simulations of (1) enzyme turnover after administration of a radio-labeled precursor amino acid, (2) the effect of ligand stabilization on the time course of total enzyme and the amount of total enzyme at steady-state, and (3) the effect of competitive inhibition alone on clearance of substrate, examined by dividing the clearance by enzyme amount.
Enzyme Turnover in the Two Pools.
Fig. 3 shows a simulation of the elimination of a pulse dose of radiolabeled enzyme formed by the administration of a radiolabeled precursor. This is accomplished by setting the rate of synthesis of enzyme to zero after enzyme is at steady-state. The simulation shows the characteristics required of the decline of enzyme. Decline of total enzyme is biphasic with half-lives of 7 and 37 hr, corresponding to the respective half-lives for pools 1 and 2 in the rat. Since newly synthesized enzyme appears first in pool 1, the amount of labeled enzyme in pool 2 rises to a maximum as it is transferred from pool 1 to pool 2. Also shown in fig. 3 is the time course of enzyme accumulation in the respective pools and the return to the pre-inhibitor/inducer condition after the end of the infusion. True steady-state is reached in the system by 200 hr after the beginning of the infusion and the pre-inhibitor/inducer condition is re-attained by 200 hr after the end of the infusion, times consistent with the longer of the CYP2E1 half-lives. However, steady-state and return to baseline is approximated much earlier in pool 1 and total CYP2E1 since the short half-life contributes to enzyme elimination unless the enzyme is completely bound.
Effect of Inhibitor/Inducer on Enzyme Amount.
Fig. 4 shows that infusion of inhibitor/inducer to steady-state (achieved at 24 hr after the start of the infusion) causes a rise in enzyme that is a function of I/Ki and not S/Km, as the model was constructed to allow evaluation of the effect of the inhibitor/inducer alone on enzyme behavior rather than to allow the contributions of inhibitor/inducer and substrate to interact.
The time course of induction of total enzyme from the beginning of the infusion of inhibitor/inducer to the return of enzyme to the basal state after the end of the infusion also is shown in fig. 4. At low steady-state inducer/inhibitor concentration, I/Ki = 0.1, enzyme induction is minimal and reaches steady-state in 50 hr. As I/Ki increases to 1, enzyme is induced by 50% over the initial value and reaches steady-state by 100 hr (the half-life becomes longer as a greater fraction of the total enzyme is stabilized). When administration of inducer is stopped, enzyme declines as inhibitor/inducer is removed and enzyme is degraded. The pre-induction level of enzyme is reached by 200 hr after the infusion is stopped, consistent with the slow half-life as CYP2E1 has been elevated in pool 2 as well as pool 1. Although, again, 50 hr after inhibitor/inducer infusion is stopped only a modest induction remains due to the influence of the short half-life on removal of CYP2E1. At I/Ki = 3, greater than 2-fold induction of the enzyme is attained at 168 hr. As I/Ki increases, the time required to reach steady-state increases because the greater the fraction of enzyme in pool 1 that is stabilized, the longer is the turnover time of total enzyme. This simulation shows that the time course of the response of enzyme is as expected based on the general scheme of the model.
Clearance Inhibition.
The interaction considered has two components that effect clearance: induction and inhibition. To assure that both are operating in the model as expected, the inhibition component was analyzed before their combined effects on clearance were considered.
The evaluation of the change in substrate clearance caused by the inhibitor/inducer requires a selection of the basal state to which the altered states of the enzyme and clearance are to be compared. As described in the mathematical description of the model, two basal states are relevant. The first, which is referred to as “case 1,” divides the clearance at any time (and any inhibitor/inducer concentration) by the clearance that would be observed if the inhibitor/inducer was suddenly removed and the substrate concentration was unchanged. Thus, in case 1, Cl0 (the basal clearance) is E·kcat/(Km+S).
The second case, “case 2,” uses as basal clearance the value of clearance in the absence of the inhibitor/inducer and at vanishingly low concentration of substrate, i.e. S/Km→0. In this case, Cl0is E·kcat/Km. This basal clearance is identified in the figures as Cl0,(S/Km→0). Case 2 allows evaluation of the net effect of inhibitor/inducer concentration and substrate concentration on substrate clearance in relation to the clearance that usually operates in therapeutics, low substrate concentration relative to Km and no interacting drug.
Case 1. The inhibition component of the interaction, Cl/E, as described by eq. 10 is shown in fig.5. For this case, at I/Ki = 0, Cl0/E0 is 1 by definition at all values of S/Km, since Cl0 is defined as E·kcat/(Km+S). At low S/Km, Cl/E decreases as I/Ki increases and inhibition becomes greater. At I/Ki = 3, the enzyme is nearly saturated with inhibitor/inducer. At this value of I/Ki, Cl/E is at the minimum at low S/Km. As S/Km increases, substrate displaces inhibitor/inducer and Cl/E increases in comparison to the clearance in the absence of inhibitor/inducer at the same concentration of substrate.
Case 2.
In fig. 6, the inhibition component of clearance, Cl/E, is shown as described by eq. 11. At low S/Km, Cl/E decreases as I/Ki increases (due to inhibition) and as S/Km increases (due to saturation of the enzyme with substrate) at low I/Ki. When I/Ki = 3, the enzyme is bound with inhibitor/inducer and Cl/E declines slightly as S/Km increases because the enzyme becomes saturated with substrate as inhibitor/inducer is displaced. Cl/E is minimum as both S/Km and I/Ki increase to 3.
Combined Effects of Inhibition and Induction on Substrate Clearance at Steady-State
The simulations presented in this section consider the situation in which inhibitor/inducer has been infused to steady-state and is continued as the combined effects of induction and inhibition on substrate clearance are evaluated. The two cases of expression of the alteration of clearance are included as above.
Case 1.
Fig. 7 shows the effect of the administration of inhibitor/inducer on the clearance of substrate expressed relative to clearance at the same concentration of substrate. When S/Km is approximately 0 and I/Ki increases, clearance decreases. Thus, at low substrate concentration the net effect of the inhibitor/inducer is inhibition of clearance as the low concentration of substrate cannot displace inhibitor/inducer from the enzyme. As S/Km is increased at I/Ki = 3, substrate can displace inhibitor/inducer and clearance rises. As S/Km approaches I/Ki, clearance becomes greater than that observed at the same concentration of substrate in the absence of inhibitor/inducer.
Case 2.
Fig. 8 shows the effect of the administration of inhibitor/inducer on the clearance of substrate expressed relative to clearance at vanishingly low substrate concentration. By definition in this case, Cl(t)/Cl0 declines as either I/Ki or S/Km increases. The most important result of this simulation is the illustration that the condition in which maximum substrate clearance will be observed is when inhibitor/inducer is absent and substrate concentration is low. In other words, as long as the inhibitor/inducer is present, clearance will not rise above the value observed in its absence at low substrate concentration even though the enzyme is induced. Although increasing substrate concentration will displace inhibitor/inducer from the enzyme (this result is apparent in case 1), the net effect of increasing substrate concentration is to saturate the enzyme and decrease clearance in comparison to the value observed at vanishingly low substrate concentration in the absence of inhibitor/inducer.
Time Course of Inhibition/Induction of Substrate Clearance as Inhibitor/Inducer is Added and Removed
The time course of substrate clearance over the entire period of exposure to inhibitor/inducer and the decline of enzyme to the basal level after inhibitor/inducer is removed is examined in this section. Fig. 9 shows the results for both cases of clearance for different values of S/Km at a steady-state I/Ki of 1.3
Case 1.
As inhibitor/inducer is introduced by constant rate infusion, clearance of substrate declines. The decline in substrate clearance is greater at the low values of S/Km, as competitive inhibition is greater when substrate concentration is lower. Substrate clearance increases over the period of infusion of inhibitor/inducer as the enzyme accumulates as a consequence of stabilization. At the end of the 350 hour infusion, substrate clearance rapidly increases as inhibitor/inducer concentration declines. The increase is greatest when S/Km is low because competitive inhibition is greatest at low S/Km and induction is a function only of inhibitor/inducer concentration, which did not vary among the simulations. The enzyme returns to the pre-induction value by 200 hr after the infusion is stopped.
Case 2.
As in case 1, substrate clearance declines as inhibitor/inducer is introduced. The decline is greatest when S/Km is 0.1 because of saturation of the enzyme by high substrate concentration. As the infusion of inhibitor/inducer is continued, clearance of substrate increases somewhat as enzyme accumulates as a consequence of stabilization. At the end of the infusion of inhibitor/inducer, clearance of substrate increases quickly as inhibitor/inducer is eliminated. The increase is most pronounced in this case when S/Km is small because the enzyme is not saturated with substrate and the degree of competitive inhibition is greatest when S/Km is low.
Fig. 10 shows the results for both expressions of the change in clearance for different values of I/Ki at S/Km of 0.1. In these simulations, the behavior of case 1 and case 2 are similar. As administration of inhibitor/inducer is initiated, substrate clearance falls as a consequence of inhibition, with the greatest decline when I/Ki is greatest. As administration of inhibitor/inducer continues, clearance rises as enzyme accumulates. At the end of the infusion of inhibitor/inducer, clearance rises abruptly as inhibitor/inducer is eliminated. The magnitude of the rise is greatest when I/Ki is greatest, as this is the condition in which induction was maximal.
Comparison of Model Predictions with Experimental Results
Fig. 11 shows the results of a simulation of the time course of the interaction between isoniazid and acetaminophen as predicted by the model in comparison to data from a series of clinical studies of the effect of isoniazid administration on the formation of NAPQI from acetaminophen and 6-hydroxychlorzoxazone from chlorzoxazone, both CYP2E1-modulated events (26-28). In each study, 300 mg isoniazid was administered once a day for several days and 500 mg acetaminophen or 750 mg chlorzoxazone was administered on three occasions, either before isoniazid was begun (26, 27) or well after the induction of CYP2E1 by isoniazid had reversed (1–2 weeks) (28), during isoniazid administration, and some time soon after the administration of isoniazid had been stopped (26-27). When given during the course of isoniazid administration, acetaminophen was usually given simultaneously with isoniazid (26, 27). In one study however (28), acetaminophen was given either with isoniazid in the morning or 12 hr after the morning dose of isoniazid. In that study, the time course of reversal of the interaction from inhibited to enhanced formation of NAPQI was evaluated as a function of isoniazid acetylator phenotype.
The simulations in fig. 11 take a value of S/Kmof 0.1. In humans, the dose of acetaminophen administered produces a maximum concentration in plasma of approximately 56 μM (28) and the dose of chlorzoxazone a maximum unbound concentration of 3 μM (55). The Km of acetaminophen for CYP2E1 is approximately 1 mM (56), and that for chlorzoxazone is 45 μM in human liver microsomes (42, 57). The half-lives of CYP2E1 have been established in the rat, but not in humans. In addition to conducting simulations with CYP2E1 half-lives of 7 and 37 hr as observed in the rat, we also performed simulations with values of 15 and 79 hr, as the half-life in humans is expected to be longer than rat. Also, the time for alcoholics to return to an uninduced CYP2E1 activity has been estimated to be 60 hr (58). The longer half-life values used provide the same distribution of enzyme between the two pools we estimated in the rat. Peak isoniazid concentration in plasma is reported to vary between 36 and 73 μM at a dose of 300 mg (59, 60), and a Ki of 40 μM for human CYP2E1 has been reported (53).
Fig. 11 shows data from the human studies and simulations with the literature parameters for I/Ki of 1.5 and 2, and CYP2E1 half-lives of 7 and 37 hr and 15 and 79 hr. The simulation with the longer half-lives and larger I/Ki best fit the data, although both agree qualitatively well with the time course of the interaction, inhibited clearance giving way to enhanced clearance.
Discussion
We have described a model of the time course of CYP2E1 inhibition/induction drug interactions that are a consequence of stabilization of the enzyme by bound ligand. The model accurately describes the biphasic degradation of enzyme that is observed in the absence of an inhibitor/inducer and the monophasic decline when the inhibitor/inducer is present. The degree of inhibition and induction is a function of inhibitor/inducer concentration and its affinity for the enzyme, while substrate concentration and Kminfluence the degree of inhibition as described by the competitive mechanism. Although the interaction is complex, the model suggests that the time course of the interaction and its magnitude can be anticipated from measurements made in vitro and knowledge of CYP2E1 half-lives and distribution between pools.
The effects of the inhibitor/inducer on the clearance of CYP2E1 substrates were examined by considering two separate baseline values of clearance, either that at the respective concentration of substrate in the absence of the inhibitor/inducer, or that at vanishingly low substrate concentration (S/Km≪1). It is clear from the former case that the magnitude of decrease in substrate clearance when the inhibitor/inducer is present is less than would be the case if induction had not occurred. The second case illustrates that at steady-state, in the presence of uninterrupted administration of the inhibitor/inducer, it is not possible to raise substrate concentration high enough to increase clearance above the value observed in the absence of the inhibitor/inducer and at low substrate concentration. While increasing substrate concentration reverses competitive inhibition it also saturates enzyme, resulting in diminished clearance in comparison to low substrate concentration.
On the other hand, it is also apparent that enhanced clearance is observed when inhibitor/inducer administration is discontinued and inhibitor/inducer is eliminated more quickly than CYP2E1. Thus, the model predictions are also generally consistent with the reversal from inhibited to enhanced NAPQI formation within the 24 hr dosing interval that we have observed in rapid acetylators of isoniazid in comparison with slow acetylators (28). Also, the model simulations are consistent with observations of the isoniazid-chlorzoxazone interaction reported by O’Shea et al. (30), in which an enhanced formation of 6-hydroxychlorzoxazone was observed in slow acetylators 2 days after the last dose of isoniazid but not in rapid acetylators. In the latter case, 6-hydroxychlorzoxazone formation at the same time had returned to the value observed in the absence of isoniazid. In both cases, the extent of inhibition and induction is generally underpredicted. The underprediction may indicate a contribution to the modulation of CYP2E1 by the hydrazine metabolites of isoniazid (53, 61-62).
Another interaction that occurs at least in part through the mechanism of substrate stabilization and seems to be of clinical significance is that between acetaminophen and ethanol (29). It is clear from these studies and the effects of ethanol on enflurane metabolism in rabbits that the formation clearance of NAPQI can be both inhibited and induced by ethanol. Thus, it would seem that simultaneous ingestion of ethanol with acetaminophen would tend to be protective against NAPQI-induced hepatotoxicity (63), while ingestion of acetaminophen when ethanol concentration has declined in the presence of CYP2E1 induction would be expected to enhance exposure to NAPQI. Derangements of glutathione homeostasis and UDP-glucuronic acid availability in the presence of fasting and alcoholic liver damage may alter these expectations (64,65).
To obtain biphasic elimination kinetics of CYP2E1 in the absence of the inhibitor/inducer in the model, it was necessary to incorporate two physically distinct pools of the enzyme into the model. The two pools share one mechanism (the slow mechanism) of enzyme degradation, but only one pool, the one into which new enzyme is synthesized, contains enzyme that is degraded by the rapid mechanism. Transfer of the enzyme from one pool to the other or binding of ligand protects enzyme from the rapid phase of degradation. Consideration of the known intracellular distribution suggested to us that the two pools might be the rough and smooth endoplasmic reticulum. However, it has not been demonstrated that a rapid degradation process operates at only one of these sites, while the slow process operates at both. Results of these simulations regarding the different time course and extent of induction in the two pools may be useful in identifying the actual pools.
The class of drug interactions that accrue through ligand stabilization of CYP2E1 may be unique to this enzyme.4 They represent an unusual case in which the effect of one drug on the other is not inhibition or induction, but both. The net result of the effect of the inhibitor/inducer on the clearance of the substrate is predictable from the synthesis and degradation kinetics of the enzyme and the affinities and concentrations of the inhibitor/inducer and substrate. The unusual feature of this interaction is that although induction of CYP2E1 begins as soon as the inhibitor/inducer is administered, enhanced clearance of substrate in relation to the pre-inhibitor/inducer condition is observed only after the inducer is removed. The isoniazid-acetaminophen interaction illustrates that it is possible for the interaction to cycle between inhibited and enhanced clearance within a dosing interval of the inhibitor/inducer. A detailed understanding of the pharmacokinetic consequences of this mechanism of interaction is likely to aid the interpretation of drug interaction studies involving CYP2E1.
Footnotes
-
Send reprint requests to: John T. Slattery, Ph.D., Department of Pharmaceutics, Box 357610, University of Washington, Seattle, WA 98195-7610. E-mail: jts{at}u.washington.edu.
-
This work was supported in part by grants GM 32165 and GM 48349 from the National Institutes of Health.
-
↵2 Since the inducing ligand induces a spin-shift consistent with occupation of the active site (36), the mechanism of inhibition most likely is competitive, as has been demonstrated with isoniazid in human liver microsomes and in a reconstituted system with expressed CYP2E1 (53). The equations for fu and fb are derived from Scatchard’s equilibrium binding equations (54).
-
↵3 For ease of comparison to steady-state simulations, the time course simulations are identified with values of I/Ki achieved at steady-state. Infusion rate of inhibitor/inducer was varied to achieve the indicated steady-state I/Ki; I/Ki did change over the course of infusion.
-
↵4 Induction of CYP3A isoforms byN-substituted imidazoles and macrolides may represent variants of this mechanism (66-68).
- Abbreviations used are::
- CYP2E1
- cytochrome P450 2E1
- NAPQI
- N-acetyl-p-benzoquinone imine
- Received March 31, 1997.
- Accepted June 19, 1997.
- The American Society for Pharmacology and Experimental Therapeutics