Abstract
The pharmacodynamic interaction between midazolam and its active metabolite α-OH-midazolam was investigated to evaluate whether estimates of relevant pharmacodynamic parameters are possible after administration of a mixture of the two. Rats were administered 10 mg/kg of midazolam, 15 mg/kg of α-OH-midazolam, or a combination of 3.6 mg/kg of midazolam and 35 mg/kg of α-OH-midazolam. Increase in the 11.5- to 30-Hz frequency band of the electroencephalogram was used as the pharmacodynamic endpoint. The pharmacodynamics of midazolam and α-OH-midazolam after combined administration were first analyzed according to an empirical and a competitive interaction model to evaluate each model’s capability in retrieving the pharmacodynamic estimates of both compounds. Both models failed to accurately estimate the true pharmacodynamic estimates of midazolam and α-OH-midazolam. The pharmacodynamic interaction was subsequently analyzed according to a new mechanism-based model. This approach is based on classical receptor theory and allows estimation of the in vivo estimated receptor affinity and intrinsic in vivo drug efficacy. The relationship between stimulus and effect is characterized by a monotonically increasing function f, which is assumed to be identical for midazolam and α-OH-midazolam. The pharmacodynamic interaction is characterized by the classical equation for the competition between two substrates for a common receptor site. This mechanism-based interaction model was able to estimate the pharmacodynamic parameters of both midazolam and α-OH-midazolam with high accuracy. It is concluded that pharmacodynamic parameters of single drugs can be estimated after a combined administration when a mechanistically valid interaction model is applied.
In many instances drug effects are the result of the combined action of two or more chemical entities, when two or more drugs are given in combination, when a drug is converted into an active metabolite, or in the situation where, for an enantiomeric drug, a racemate of the two enantiomers is administered or metabolic interconversion occurs. In these situations the intensity of the pharmacological response can generally not be predicted by simply adding the effects of the two compounds separately; rather complex interaction models may be required. In this respect a clear distinction must be made between competitive interactions (where the two compounds interact at the same receptor) and noncompetitive interactions (where the two compounds bind to a separate receptor; Danhof and Mandema, 1992).
Many psychotropic drugs are converted into active metabolites (Caccia and Garratini, 1990) and these metabolites may contribute significantly to the pharmacological response (Danhof and Mandema, 1995). A well known example is benzodiazepines, where active metabolites have been demonstrated to contribute significantly to the effect (Crevoisier et al., 1983; Garzone and Kroboth, 1989; Mandema et al., 1992b).
When a drug is converted into an active metabolite, generally the situation arises where two active species are present which, as result of their close structural similarity, compete for the same receptor. Therefore, typically a competitive interaction model is required to characterize the combined effect of a drug and its active metabolite. Recently the modeling of the pharmacodynamic interactions between benzodiazepines has been the subject of a number of investigations. Specifically, empirical interaction models have been proposed to characterize the interactions between a benzodiazepine full agonist on one hand and benzodiazepine partial, competitive, or inverse (ant)agonists on the other (Mandema et al., 1992a,c). So far, however, the modeling of the pharmacodynamic interactions between two full agonists (or between parent drug and its active metabolite) has not been reported.
In this paper a new mechanism-based model is proposed for the modeling of competitive drug-drug interactions between full agonists in vivo that can also be used to characterize the pharmacodynamic interactions between a parent drug and its active metabolite. This model is based on receptor theory and contains the pharmacodynamic parameters in vivo estimated receptor affinity (K PD) and intrinsic in vivo drug efficacy (e PD). The relationship between stimulus and effect is characterized by a monotonically increasing function f, which is assumed to be identical for both parent drug and metabolite. The pharmacodynamic interaction is characterized by the classical equation for the competition between two ligands for a common receptor site. This new model is applied in a study on the pharmacodynamic interaction between midazolam and its active metabolite α-OH-midazolam in rats, using quantitative electroencephalogram (EEG) parameters as a pharmacodynamic endpoint. The performance of the new model is compared to that of the previously proposed, more empirical interaction models.
Materials and Methods
Animals.
Three groups of 6 to 8 male Wistar rats (Sylvius Laboratory Breeding Facility, Leiden, the Netherlands) weighing (mean ± S.E.) 261 ± 7 g were used in the study. The animals were kept individually in plastic cages with a normal 12-h light/dark cycle and were fed on a commercially available diet (Standard Laboratory Rat, Mouse and Hamster Diets, RMG-TM, Hope Farms, Woerden, the Netherlands) and water ad libitum. From the night before the experiment onward, the animals were deprived of food but had free access to water. For the measurement of EEG signals, chronic cortical EEG electrodes were implanted into the skull of the animals 1 week before the kinetic-dynamic experiments as described previously (Mandema et al., 1991b). One day before the experiment, indwelling cannulae were implanted in the right jugular vein (for drug administration) and the right femoral artery (for blood sample collection).
Drug Dosage and Blood Sampling.
Rats received 10 mg/kg midazolam, 15 mg/kg α-OH-midazolam, or a combined administration of 3.6 mg/kg midazolam and 35 mg/kg α-OH-midazolam during a 15-min infusion. Drugs were dissolved in 0.9% saline with the aid of an equimolar quantity of hydrochloric acid. To determine the pharmacokinetics of midazolam and α-OH-midazolam, blood samples of 100 or 200 μl (near the end of the experiment) were collected at fixed time intervals after drug administration over a period of 280 min. After the experiment the animals were sacrificed and a final blood sample was obtained by aortic puncture to be used for protein binding measurements. Heparinized blood samples were centrifuged and plasma was separated and stored at −35°C until the time of analysis.
EEG Measurements.
The output from bipolar EEG leads was continuously recorded using a Nihon Kohden EEG system consisting of a bioelectric input box JB-682G (Mihon Kohden Corporation, Tokyo, Japan), bioelectric amplifier AB-621G, and bioelectric input panel PB-680G. The low pass filter was set at 100 Hz, the time constant at 0.3 s. During the course of the experiment the animals were forced to walk in a slowly rotating drum to prevent spontaneous fluctuations in the level of vigilance (Mandema et al., 1991b). EEG recordings were commenced 15 min before drug administration for baseline determination. Two EEG leads, the fronto-central and central occipital lead on the left hemisphere, were quantified on-line by aperiodic analysis (Gregory and Pettus, 1986) as described previously (Mandema et al., 1991b). The amplitudes (μV/s) in the 11.5- to 30-Hz frequency band of the fronto-central lead as change over baseline were used as a measure of drug effect intensity.
Drug Analysis and Plasma Protein Binding.
Plasma concentrations of midazolam, α-OH-midazolam, and 4-OH-midazolam were determined by a gas chromatographic assay using electron-capture detection as described previously (Mandema et al., 1992b). Detection limits were 20 ng/ml for all three compounds. The extent of plasma protein binding of midazolam and α-OH-midazolam was determined for each individual animal by ultrafiltration at 37°C using the Amicon Micropartition System (Amicon Division, Danvers, MA) as described previously (Mandema et al., 1991b).
Pharmacokinetics.
Data were analyzed with the NONMEM computer program as developed by Beal and Sheiner (NONMEM project group, University of California at San Francisco, San Francisco, CA). The plasma concentration-time profiles of the drugs after i.v. infusion were described by a polyexponential equation for i.v. infusion:
Pharmacodynamics.
Pharmacodynamic data were analyzed with the NONMEM computer program as developed by Beal and Sheiner (NONMEM project group). Individual pharmacokinetic profiles were used to drive the pharmacodynamic fitting. Concentrations of midazolam and α-hydroxy-midazolam were directly linked to the EEG effect and characterized according to the sigmoidal maximal drug effect (Emax) model (Holford and Sheiner, 1982):
Empirical and Competitive Interaction Modeling.
The results obtained upon the combined administration of midazolam and α-OH-midazolam were analyzed on the basis of two different interaction models: an empirical interaction model and a pure competitive model. First, the data obtained upon administration of midazolam separately and upon administration of the combination were used to retrieve the pharmacodynamic parameters of midazolam and α-OH-midazolam. Next, the data obtained upon separate administration of α-OH-midazolam and the combination were used to estimate the pharmacodynamic parameters of both benzodiazepines.
An empirical model for competitive interaction between two ligands acting at the same receptor was proposed by Holford and Sheiner (1982):
To overcome the limitations of the above empirical model, a pure competitive interaction model can be derived, assuming that effects after single administration of the interacting compounds can be characterized on the basis of the sigmoidal Emaxmodel (eq. 2). Competitive models are of an additive nature. The interaction between two drugs A and B is additive when the isobole, i.e., the curve of drug concentration pairs that result in the same specific intensity of pharmacological effect, is a straight line. Mathematically this is represented by the additivity isobole:
Only in the special case that the pharmacodynamic parameters Emax and N of the interacting compounds are equal (EmaxA = EmaxB andN
A = N
B) the following algebraic solution for the competitive model can be derived:
Mechanism-Based Interaction Modeling.
Subsequently, the data were analyzed by a new mechanism-based model where effect is viewed to be a function of the benzodiazepine receptor binding-induced stimulus (S). This mechanism-based approach allows the event at the receptor level to be separated from the multitude of processes as result of the formation of the drug-receptor complex. This has the important advantage that empirical factors such as the Hill factor are omitted from the model. The details of the mechanism-based modeling approach are presented as an to this paper. When using this mechanism-based approach, the relationship between drug concentration and effect is characterized by the following equation:
Results
Pharmacokinetics and Pharmacodynamics of Midazolam and α-OH-Midazolam.
The plasma concentration-time profile and the EEG-time profile after 10 mg/kg of midazolam, 15 mg/kg α-OH-midazolam, or the combined administration of 3.6 mg/kg of midazolam and 35 mg/kg α-OH-midazolam are shown in Figs.1 and 2, respectively. After midazolam administration no measurable concentrations of the metabolites α-OH-midazolam and 4-OH-midazolam were observed. After α-OH-midazolam administration there were no measurable concentrations of midazolam and 4-OH-midazolam. For both compounds the plasma concentration versus time profiles were most adequately described by a biexponential equation. Pharmacokinetic parameters of midazolam and α-OH-midazolam were not different when estimated after separate or combined administration indicating the absence of a pharmacokinetic interaction. The values of the pharmacokinetic parameters and the free fraction in plasma are summarized in Table 1. The concentration EEG-effect relationships after separate administration of midazolam and α-OH-midazolam were analyzed on basis of the sigmoidal Emax pharmacodynamic model (Fig.3). The EC50 was 62 ± 10 ng/ml for midazolam and 406 ± 31 ng/ml for α-OH-midazolam (Table 2). Values (mean ± S.E.) for Emax and Hill factor were 191 ± 6% and 190 ± 4% and 0.87 ± 0.11 and 1.29 ± 0.12 for midazolam and α-OH-midazolam, respectively.
Empirical and Competitive Interaction Modeling.
The performance of the empirical interaction model (eq. 3) and the competitive interaction model (eq. 8) was evaluated by comparing the by each model derived pharmacodynamic estimates with those as estimated upon separate administration. First, the data obtained upon administration of midazolam separately and upon administration of the combination were used to retrieve the pharmacodynamic parameters of midazolam and α-OH-midazolam by fitting the pharmacodynamic models to both data sets simultaneously. Next, the data obtained upon separate administration of α-OH-midazolam and the combination were used to estimate the pharmacodynamic parameters of both benzodiazepines. The results of this analysis are shown in Table 3. Interestingly, upon visual inspection of the fits no remarkable differences were observed. Comparison with the estimated values obtained upon separate administration of the compounds (Table 2) shows, however, that the empirical interaction model is unable to accurately estimate values of the pharmacodynamic parameters of midazolam and α-OH-midazolam. With the competitive interaction model, accurate estimates of the pharmacodynamic parameters of midazolam but not of α-OH-midazolam were obtained. Here again by visual inspection the fits appeared good, without accurately estimating pharmacodynamic parameters. The competitive interaction model was also tested in a slightly modified form by assuming the Emax values and Hill factors of midazolam and α-OH-midazolam to be equal (eq. 8). This simplified model also failed to accurately estimate the pharmacodynamic parameters of midazolam and α-OH-midazolam.
Mechanism-Based Modeling.
The values of the pharmacodynamic parameters of midazolam and α-OH-midazolam according to the mechanism-based model are summarized in Table 3. Similar values of efficacy for midazolam and α-OH-midazolam were observed. The value of the K PD was 43 ± 5 ng/ml for midazolam and 258 ± 25 ng/ml for α-OH-midazolam. On the basis of the mechanism-based interaction model, realistic estimates of the pharmacodynamic parameters of both midazolam and α-OH-midazolam were obtained where using the data from the combined administration (Table3).
Discussion
The fact that many (psychotropic) drugs are converted into active metabolites has major implications for pharmacokinetic/pharmacodynamic modeling. Typically in this situation the effect cannot be predicted by adding the effects of the parent drug and the metabolite separately, and rather complex interaction models are required. In this investigation it is shown that in the situation where a drug is converted into an active metabolite, accurate prediction of the pharmacological response is only feasible if a mechanism-based interaction model is used.
In the present study midazolam was used as a model drug. It is well established that in humans this drug is converted into an active metabolite, α-OH-midazolam, and that particularly upon oral administration this metabolite contributes significantly to the pharmacological response (Crevoisier et al., 1983; Mandema et al., 1992b). In the rat no α-OH-midazolam is formed (Mandema et al., 1991b). This allows the pharmacodynamics of midazolam to be characterized in the absence of α-OH-midazolam and also offers the opportunity to characterize the pharmacodynamic interaction between parent drug and metabolite by separate and combined administration of the two. Quantitative EEG analysis was used as a pharmacodynamic endpoint as this provides a continuous and relevant measure of the effect on the central nervous system (Mandema et al., 1991a, 1992d). The separate administrations of midazolam and α-OH-midazolam allowed the true pharmacodynamic parameter estimates to be determined (Tables2 and 3). When analyzing the concentration-EEG effect relationships of both compounds on the basis of the sigmoidal Emaxpharmacodynamic model, a 7-fold difference in EC50 was observed whereas the values of the Emax were identical, indicating that both compounds act as full agonists at the GABA-benzodiazepine receptor complex in vivo. To characterize the interaction between the two compounds a number of different interaction models were studied. An important issue is how to determine which of the different models most accurately characterizes the pharmacodynamic interaction between midazolam and α-OH-midazolam (i.e., two full agonists). The most direct strategy would be to fit each model to the single drug data and then to determine how well the model predicts the combined drug effect. A limitation of this approach is, however, that it requires the comparison of the concentration-effect profiles. In general the various models did not differ to a large extent in their ability to forecast the response to the simultaneous administration. Furthermore the statistical power of existing techniques to detect differences in the overall concentration-effect profiles is rather poor. Thus by using this approach the power to detect differences in the performance of the various models is limited. Therefore we have followed another approach and determined how well the pharmacodynamic parameters of a single compound (either parent drug or metabolite) can be determined from a combined administration. An additional attractive feature of this approach is, that it may be generally applicable to determine the pharmacological activity of drug metabolites. Particularly for drugs which are subject to an extensive “first-pass” effect, high concentrations of both the parent compound and the metabolite are generally observed after oral administration, whereas upon i.v. administration only the parent compound is present in significant concentrations. The approach from the present investigation can then be used to estimate the pharmacodynamic parameters of the metabolite.
The empirical pharmacodynamic interaction model (eq. 3) failed to accurately estimate the EC50 for both midazolam and α-OH-midazolam (Table 2). Interestingly in other interaction studies with benzodiazepines, this model performed quite well (Mandema et al., 1992a,c). This can be explained by the fact that in those investigations, the benzodiazepines differed significantly in their intrinsic activity, which allowed for a different study design. In these previous investigations the concentration of one of the benzodiazepines (typically the full agonist) was brought to steady state by administering a continuous infusion. Subsequently the second benzodiazepine (partial/inverse agonist or competitive antagonist) was administered in a short infusion and the effect intensity was observed as it went from maximal to baseline (in the case of a competitive antagonist) or to the maximum effect of the second compound (in the case of a partial or inverse agonist), and subsequently, because of the decay in concentration, back to the maximal effect of the first compound. The model was then used to estimate the pharmacodynamic parameters of only the second compound. In the current investigation the empirical interaction model is used to characterize the pharmacodynamics of both compounds, which in addition have a very similar Emax value. Another important factor is that, on theoretical grounds, the model is not valid for compounds with a Hill factor that is different from 1. In that situation the model does not predict an additive interaction (as it should in the case of competition of two ligands for the same receptor) but, depending on the value of the Hill factor, either a synergistic or an antagonistic interaction. To further examine the role differences in the Hill factor between the compounds, the data were also analyzed assuming the values of either the Emax or the Hill factor to be identical for both compounds. These simplified forms of the empirical model however also failed to retrieve the accurate pharmacodynamic parameter estimates. This means that the empirical interaction model is of little value to characterize the pharmacodynamic interaction between two full agonists.
The data from the combined administration of midazolam and α-OH-midazolam were also analyzed according to the competitive interaction model on the basis of eq. 8, where it is assumed that both the values of Emax and the Hill factor are identical. This model has the advantage over the empirical interaction model in that it is correctly derived from the sigmoidal Emax model and that it predicts an additive interaction also when the value of the Hill factor is different than 1. This competitive interaction model was found to accurately estimate the pharmacodynamic parameters of midazolam, but failed to retrieve the pharmacodynamic parameters of α-OH-midazolam correctly (Table 2). Presumably this can be explained by the fact that the true values of the Hill factor are different for both compounds (Table 2). A more general form of the competitive interaction model (eq. 7) was also tested. This model differs from the previous model (eq. 8) in that it does not have the assumption of equal values of Emax and the Hill factor. This model also failed to accurately estimate the pharmacodynamic parameters of both midazolam and α-OH-midazolam. This is most likely due to the fact that this model, where both the Emax and the Hill factor are to be estimated for the two compounds, is overparameterized. This means that without the addition of extra information it is impossible to obtain realistic pharmacodynamic parameter estimates.
In the present investigation the mechanism-based interaction model was the only model that was able to retrieve accurate estimates of the pharmacodynamic parameters of both midazolam and α-OH-midazolam (Table 3). A unique feature of this mechanism-based model is that it does not need a Hill factor to account for a sigmoidal concentration-effect relationship as it incorporates sigmoidicity through the nonlinear stimulus-response relationship. As a result the interaction between the compounds at the receptor can be described by a theoretically correct competitive interaction model that is based on receptor theory. An interesting feature of this approach is that it in principle also allows characterization of the interaction between more than two drugs (Van den Brink, 1977). A limitation of the mechanism-based modeling approach is that it requires detailed knowledge with regard to the general shape of the stimulus-response relationship. This requires the simultaneous analysis of the concentration-effect relationships of at least two other compounds as is shown in the . This also implies that a considerable amount of extra information is incorporated in the interaction model. Another important feature is that the values of only two pharmacodynamic parameters are estimated (K PD ande PD.) instead of three (EC50, Emax, and Hill factor) in most other models. This explains why overparameterization is not a problem with this particular model. It seems unlikely, however, that the incorporation of additional information on the stimulus-effect relationship explains entirely the improved performance of the mechanism-based interaction model relative to the empirical and the competitive interaction model. The fact that the mechanism-based model utilizes a theoretically correct model to characterize the competitive interaction between midazolam and α-OH-midazolam at the receptor is probably the most important factor.
In conclusion, the results from the present investigation show that the pharmacodynamic interactions between two full agonists (i.e., a parent drug + an active metabolite) in vivo can only be characterized adequately on the basis of the mechanism-based pharmacodynamic interaction model.
Appendix
Previously published data on the concentration-EEG effect relationships of flunitrazepam, midazolam, clobazam, and oxazepam (Mandema et al., 1991) were re-analyzed according to a new model that is based on the receptor theory proposed by Stephenson (1956) andFurchgott (1966). The model combines two independent parts to describe drug action: a drug-specific part, reflecting the drug’s affinity to and intrinsic efficacy at the receptor and a system-specific part describing the relationship between receptor occupancy and effect. The drug-specific part is characterized by the following equation:
To apply this model to in vivo drug concentration effect relationships, the following simplifications must be made. First, the total amount of functional receptors cannot be measured in vivo, which means that only the product of ε and Rt can be estimated. Second, the absolute value of the product of ε and Rt can also not be determined. Only the relative value can be estimated by setting the value of εRt product of the drug reaching the highest maximum effect to one. The relationship between drug concentration and effect can then be characterized by the following equation (Kenakin, 1993b):
Simultaneous analysis of data on the concentration-EEG effect relationships of flunitrazepam, midazolam, oxazepam, and clobazam (Mandema et al., 1991) allowed independent estimation ofK PD,e PD, and f. Thereby the relationship f between initial stimulus and observed pharmacological effect was assumed to be identical for the four different benzodiazepines and characterized by a natural cubic spline (DeBoor, 1978). The knots of the spline were placed at equidistant intervals on the stimulus axis (i.e., equidistant intervals between 0 and 1). The number of knots was determined on basis of the Log-Likelihood criterion (p < .01) and visual inspection of the fits. The data were analyzed using a pooled data approach within the NONMEM computer program (NONMEM project group).
The mean values of the K PD ande PD of flunitrazepam, midazolam, oxazepam, and clobazam are given in Table4. The values of the binding constant for the interaction with the benzodiazepine receptor in vitro are presented as well. Table 4 shows that the in vivo estimatedK PD values based on free drug concentrations (K PD, u) are close to the in vitro estimated receptor affinityK i and that the ratio betweenK PD, u andK i is reasonably constant. This indicates that on basis of the mechanism-based model realistic estimates of the potency of the benzodiazepines are obtained. The relationship between stimulus and EEG effect is shown in Fig.4. It is nonlinear with relatively small increase in EEG effect at stimulus intensities less than 0.2 and an almost linear but steeper increase in EEG effect at stimulus intensities between 0.2 and 1.
Footnotes
-
Send reprint requests to: Dr. M. Danhof, Leiden/Amsterdam Center for Drug Research, Division of Pharmacology, University of Leiden, P.O. Box 9503, 2300 RA Leiden, the Netherlands. E-mail:m.danhof{at}lacdr.leidenuniv.nl
- Abbreviations:
- EEG
- electroencephalogram
- Cl
- clearance
- Vss
- volume of distribution at steady-state
- Emax
- maximal drug effect
- N
- constant expressing the sigmoidicity of the concentration-effect relationship
- KPD
- in vivo estimated receptor affinity
- ePD
- in vivo drug efficacy
- Received June 5, 1998.
- Accepted December 18, 1998.
- The American Society for Pharmacology and Experimental Therapeutics