Identifying clinically relevant sources of variability: The clopidogrel challenge.

High interindividual variability in clinical outcomes following clopidogrel's standard dosing regimen continues to be a challenge even two decades after its approval. CYP2C19 polymorphisms, obesity, older age, diabetes, and drug-drug interactions have been identified as risk factors for adverse events and treatment failure. We conducted a mechanism-based pharmacokinetic/pharmacodynamic analysis, where we integrated knowledge on in vitro enzyme kinetic, physiological, genetic, and demographic information to characterize changes in platelet reactivity from baseline following clopidogrel antiplatelet therapy. When considering the combined impact of these covariates, our analysis results indicate that higher maintenance doses are required for CYP2C19 intermediate metabolizers and poor metabolizers compared to extensive metabolizers and that respective maintenance doses have to be further increased for obese subjects for each of these CYP2C19 phenotypes. In addition, interindividual differences in the fraction absorbed and the CES1 activity were identified as sources of interindividual differences in clopidogrel's active metabolite concentrations and, thus, platelet reactivity.

Understanding the Behavior of Systems Pharmacology Models Using Mathematical Analysis of Differential Equations: Prolactin Modeling as a Case Study.

In this tutorial, we introduce basic concepts in dynamical systems analysis, such as phase-planes, stability, and bifurcation theory, useful for dissecting the behavior of complex and nonlinear models. A precursor-pool model with positive feedback is used to demonstrate the power of mathematical analysis. This model is nonlinear and exhibits multiple steady states, the stability of which is analyzed. The analysis offers insight into model behavior and suggests useful parameter regions, which simulations alone could not.

Mechanism-based modeling of the pharmacodynamic interaction of alphaxalone and midazolam in rats.

The objective of the present investigation was to characterize the pharmacodynamic interaction between the synthetic neuroactive steroid alphaxalone and the benzodiazepine midazolam. The time course of the electroencephalographic (EEG) effect (11.5-30 Hz) was determined in rats in conjunction with plasma concentrations. Alphaxalone was administered as a continuous intravenous infusion of 0, 1.2, 2.2, or 5.2 mg over 360 min. Midazolam was administered as a 5-min intravenous bolus infusion of 4 mg.kg-1. The pharmacokinetic profiles of both drugs were described by a two-compartment model. No pharmacokinetic interaction was observed. The EEG effect versus time profiles of midazolam and alphaxalone, when administered separately and in combination, were modeled on the basis of the recently proposed mechanism-based pharmacokinetic/pharmacodynamic model for GABAA receptor modulators, which contains separate expressions to describe the drug-receptor interaction and the stimulus-response relationship. The pharmacodynamic interaction between alphaxalone and midazolam was best characterized using an independent drug-drug interaction model without an expression for allosteric modulation of the effect of midazolam by alphaxalone. The final model contained an exponential expression to account for acute functional adaptation to the EEG effect upon continuous infusion of alphaxalone. The mechanism-based analysis showed that this functional adaptation is best explained by a change in the system-specific stimulus-response relationship, rather than the drug-receptor activation process. It is concluded that the pharmacodynamic interaction between alphaxalone and midazolam in vivo is best described using an independent interaction model without allosteric modulation.

Dose-dependent EEG effects of zolpidem provide evidence for GABA(A) receptor subtype selectivity in vivo.

Zolpidem is a nonbenzodiazepine GABA(A) receptor modulator that binds in vitro with high affinity to GABA(A) receptors expressing alpha(1) subunits but with relatively low affinity to receptors expressing alpha(2), alpha(3), and alpha(5) subunits. In the present study, it was investigated whether this subtype selectivity could be detected and quantified in vivo. Three doses (1.25, 5, and 25 mg) of zolpidem were administered to rats in an intravenous infusion over 5 min. The time course of the plasma concentrations was determined in conjunction with the change in the beta-frequency range of the EEG as pharmacodynamic endpoint. The concentration-effect relationship of the three doses showed a dose-dependent maximum effect and a dose-dependent potency. The data were analyzed for one- or two-site binding using two pharmacodynamic models based on 1) the descriptive model and 2) a novel mechanism-based pharmacokinetic/pharmacodynamic (PK/PD) model for GABA(A) receptor modulators that aims to separates drug- and system-specific properties, thereby allowing the estimation of in vivo affinity and efficacy. The application of two-site models significantly improved the fits compared with one-site models. Furthermore, in contrast to the descriptive model, the mechanism-based PK/PD model yielded dose-independent estimates for affinity (97 +/- 40 and 33,100 +/- 14,800 ng x ml(-1)). In conclusion, the mechanism-based PK/PD model is able to describe and explain the observed dose-dependent EEG effects of zolpidem and suggests the subtype selectivity of zolpidem in vivo.

Mechanism-based pharmacokinetic/pharmacodynamic modeling of the electroencephalogram effects of GABAA receptor modulators: in vitro-in vivo correlations.

A mechanism-based pharmacokinetic-pharmacodynamic (PK/PD) model for neuroactive steroids, comprising a separate characterization of 1) the receptor activation process and 2) the stimulus-response relationship, was applied to various nonsteroidal GABAA receptor modulators. The EEG effects of nine prototypical GABAA receptor modulators (six benzodiazepines, one imidazopyridine, one cyclopyrrolone, and one beta-carboline) were determined in rats in conjunction with plasma concentrations. Population PK/PD modeling revealed monophasic concentration-EEG effect relationships with large differences in potency (EC50) and intrinsic activity between the compounds. The data were analyzed on the basis of the mechanism-based PK/PD model for (synthetic) neuroactive steroids on the assumption of a single and unique stimulus-response relationship. The model converged yielding estimates of both the apparent in vivo receptor affinity (KPD) and the in vivo intrinsic efficacy (ePD). The values of KPD ranged from 0.41 +/- 0 ng.ml(-1) for bretazenil to 436 +/- 72 ng.ml(-1) for clobazam and the values for e(PD) from -0.27 +/- 0 for methyl 6,7-dimethoxy-4-ethyl-beta-carboline-3-carboxylate to 0.54 +/- 0.02 for diazepam. Significant linear correlations were observed between KPD for unbound concentrations and the affinity in an in vitro receptor bioassay (r = 0.93) and between e(PD) and the GABA-shift in vitro (r = 0.95). The findings of this investigation show that the in vivo effects of nonsteroidal GABAA receptor modulators and (synthetic) neuroactive steroids can be described on the basis of a single unique transducer function. In this paradigm, the nonsteroidal GABAA receptor modulators behave as partial agonists relative to neuroactive steroids.

Neuroactive steroids differ in potency but not in intrinsic efficacy at the GABA(A) receptor in vivo.

The objective of the present investigation was to characterize the in vivo EEG effects of (synthetic) neuroactive steroids on the basis of a recently proposed mechanism-based pharmacokinetic/pharmacodynamic (PK/PD) model. After intravenous administration, the time course of the EEG effect of pregnanolone, 2beta-3alpha-5alpha-3-hydroxy-2-(2,2-dimethylmorpholin-4-yl)-pregnan-11,20-dione (ORG 21465), 2beta-3alpha-5alpha-21-chloro-3-hydroxy-2-(4-morpholinyl)-pregnan-20-one (ORG 20599), and alphaxalone was determined in conjunction with plasma concentrations in rats. For each neuroactive steroid the PK/PD correlation was described on the basis of a two-compartment pharmacokinetic model with an effect compartment to account for hysteresis. The observed concentration EEG effect relationships were biphasic and characterized with a mechanism-based pharmacodynamic model, which is based on a separation between the receptor activation process and the stimulus-response relationship. A single unique biphasic stimulus-response relationship could be identified for all neuroactive steroids, which was successfully described by a parabolic function. The receptor activation process was described by a hyperbolic function. Estimates for the maximum activation (e(PD)) were similar for the different neuroactive steroids but values of the potency estimate (K(PD)) ranged from 157 +/- 16 ng. ml(-1) for pregnanolone, 221 +/- 83 ng. ml(-1) for ORG 20599, and 483 +/- 42 ng. ml(-1) for alphaxalone to 1619 +/- 208 ng. ml(-1) for ORG 21465. A statistically significant correlation was observed between the in vivo potency and the IC(50) in an in vitro [(35)S]t-butylbicyclophosphorothionate binding assay (r = 0.91). It is concluded that the new PK/PD model constitutes a new mechanism-based approach to the quantification of the effects of (synthetic) neuroactive steroids in vivo effects. The results show that the neuroactive steroids differ in potency but not in intrinsic efficacy at the GABA(A) receptor in vivo.

Mechanism-based pharmacokinetic-pharmacodynamic modeling of concentration-dependent hysteresis and biphasic electroencephalogram effects of alphaxalone in rats.

The neuroactive steroid alphaxalone reveals a complex biphasic concentration-effect relationship using the 11.5 to 30 Hz frequency band of the electroencephalogram (EEG) as biomarker. The purpose of the present investigation was to develop a mechanism-based pharmacokinetic-pharmacodynamic model to describe this observation. The proposed model is based on receptor theory and aims to separate the drug-receptor interaction from the transduction of the initial stimulus into the observed biphasic response. Individual concentration-time courses of alphaxalone were obtained in combination with continuous recording of the EEG parameter. Alphaxalone was administered intravenously in various dosages. The pharmacokinetics were described by a two-compartment model, and parameter estimates for clearance, intercompartmental clearance, volume of distribution 1 and 2 were 158 +/- 29 ml. min(-1). kg(-1), 143 +/- 31 ml. min(-1). kg(-1), 122 +/- 20 ml. kg(-1) and 606 +/- 48 ml. kg(-1), respectively. Concentration-effect relationships exhibited a biphasic pattern and delay in onset of effect. The hysteresis was described on the basis of an effect-compartment model with C(max) as covariate. The pharmacodynamic model consisted of a receptor model, featuring a monophasic saturable receptor activation model in combination with a biphasic stimulus-response model. The in vivo affinity (K(PD)) was estimated at 432 +/- 26 ng. ml(-1). Unique parameter estimates were obtained that were independent of the dose and the duration of the infusion. In conclusion, we have shown that this mechanism-based approach, which separates drug- and system-related properties in vivo, was successfully applied for the characterization of the biphasic effect versus time patterns of alphaxalone. The model should be of use in the characterization of other biphasic responses.

A set-point model with oscillatory behavior predicts the time course of 8-OH-DPAT-induced hypothermia.

Agonists for the 5-hydroxytryptamine (HT)(1A) receptor induce a hypothermic response that is believed to occur by lowering of the body's set-point temperature. We have developed a physiological model that can be used to predict the complex time course of the hypothermic response after administration of 5-HT(1A) agonists to rats. In the model, 5-HT(1A) agonists exert their effect by changing heat loss through a control mechanism with a thermostat signal that is proportional to the difference between measured and set-point temperature. Agonists exert their effect in a direct concentration-dependent manner, with saturation occurring at higher concentrations. On the basis of simulations, it is shown that, depending on the concentration and the intrinsic efficacy of a 5-HT(1A) agonist, the model shows oscillatory behavior. The model was successfully applied to characterize the complex hypothermic response profiles after administration of the reference 5-HT(1A) agonists R-8-hydroxy-2-(di-n-propylamino)tetralin (R-8-OH-DPAT) and S-8-OH-DPAT. This analysis revealed that the observed difference in effect vs. time profile for these two reference agonists could be explained by a difference in in vivo intrinsic efficacy.