Pharmacokinetic-Pharmacodynamic Modeling of the Anticonvulsant and Electroencephalogram Effects of Phenytoin in Rats

• Published in: The Journal of pharmacology and experimental therapeutics 1998-02, Vol.284 (2), p.460-466
• Main Authors: Della Paschoa, O E, Mandema, J W, Voskuyl, R A, Danhof, M
• Format: Article
• Language: English
• Subjects: Animals; Anticonvulsants - pharmacokinetics; Anticonvulsants - pharmacology; Blood Proteins - metabolism; Dose-Response Relationship, Drug; Electroencephalography; Male; Models, Biological; Phenytoin - blood; Phenytoin - pharmacokinetics; Phenytoin - pharmacology; Protein Binding; Rats; Rats, Wistar
• ISSN: 0022-3565; 1521-0103

In this study a pharmacokinetic-pharmacodynamic model is proposed for drugs with nonlinear elimination kinetics. We applied such an integrated approach to characterize the pharmacokinetic-pharmacodynamic relationship of phenytoin. In parallel, the anticonvulsant effect and the electroencephalogram (EEG) effect were used to determine the pharmacodynamics. Male Wistar-derived rats received a single intravenous dose of 40 mg · kg −1 phenytoin. The increase in the threshold for generalized seizure activity (TGS) was used as the anticonvulsant effect and the increase in the total number of waves in the 11.5 to 30 Hz frequency band was taken as the EEG effect measure. Phenytoin pharmacokinetics was described by a saturation kinetics model with Michaelis-Menten elimination. V max and K m values were, respectively, 386 ± 31 μg · min −1 and 15.4 ± 2.2 μg · ml −1 for the anticonvulsant effect in the cortical stimulation model and 272 ± 31 μg · min −1 and 5.9 ± 0.7 μg · ml −1 for the EEG effect. In both groups, a delay to the onset of the effect was observed relative to plasma concentrations. The relationship between phenytoin plasma concentrations and effect site was estimated by an equilibration kinetics routine, yielding mean k e0 values of 0.108 and 0.077 min −1 for the anticonvulsant and EEG effects, respectively. The EEG changes in the total number of waves could be fitted by the sigmoid E max model, but E max values could not be estimated for the nonlinear relationship between concentration and the increase in TGS. An exponential equation ( E = E 0 + B n  ·  C n ) derived from the sigmoid E max model was applied to describe the concentration-anticonvulsant effect relationship, under the assumption that E max values cannot be reached within acceptable electric stimulation levels. This approach yielded a coefficient ( B ) of 2.0 ± 0.4 μA · ml · μg −1 and an exponent ( n ) of 2.7 ± 0.9. The derived EC 50 value of 12.5 ± 1.3 μg · ml −1 for the EEG effect coincides with the “therapeutic range” in humans.