Development of a population pharmacokinetic model for carbamazepine based on sparse therapeutic monitoring data from pediatric patients with epilepsy.
BACKGROUND
Population models can be important extensions of therapeutic drug monitoring (TDM), as they allow estimation of individual pharmacokinetic parameters based on a small number of measured drug concentrations.


OBJECTIVE
This study used a Bayesian approach to explore the utility of routinely collected and sparse TDM data (1 sample per patient) for carbamazepine (CBZ) monotherapy in developing a population pharmacokinetic (PPK) model for CBZ in pediatric patients that would allow prediction of CBZ concentrations for both immediate- and controlled-release formulations.


METHODS
Patient and TDM data were obtained from a pediatric neurology outpatient database. Data were analyzed using an iterative 2-stage Bayesian algorithm and a nonparametric adaptive grid algorithm. Models were compared by final log likelihood, mean error (ME) as a measure of bias, and root mean squared error (RMSE) as a measure of precision.


RESULTS
Fifty-seven entries with data on CBZ monotherapy were identified from the database and used in the analysis (36 from males, 21 from females; mean [SD] age, 9.1 [4.4] years [range, 2-21 years]). Preliminary models estimating clearance (Cl) or the elimination rate constant (K(el)) gave good prediction of serum concentrations compared with measured serum concentrations, but estimates of Cl and K(el) were highly correlated with estimates of volume of distribution (V(d)). Different covariate models were then tested. The selected model had zero-order input and had age and body weight as covariates. Cl (L/h) was calculated as K(el) . V(d), where K(el) = [K(i) - (K(s) . age)] and V(d) = [V(i) + (V(s) . body weight)]. Median parameter estimates were V(i) (intercept) = 11.5 L (fixed); V(s) (slope) = 0.3957 L/kg (range, 0.01200-1.5730); K(i) (intercept) = 0.173 h(-1) (fixed); and K(s) (slope) = 0.004487 h(-1) . y(-1) (range, 0.0001800-0.02969). The fit was good for estimates of steady-state serum concentrations based on prior values (population median estimates) (R = 0.468; R(2) = 0.219) but was even better for predictions based on individual Bayesian posterior values (R(2) = 0.991), with little bias (ME = -0.079) and good precision (RMSE = 0.055).


CONCLUSIONS
Based on the findings of this study, sparse TDM data can be used for PPK modeling of CBZ clearance in children with epilepsy, and these models can be used to predict Cl at steady state in pediatric patients. However, to estimate additional pharmacokinetic model parameters (eg, the absorption rate constant and V(d)), it would be necessary to combine sparse TDM data with additional well-timed samples. This would allow development of more informative PPK models that could be used as part of Bayesian dose-individualization strategies.
School of Pharmacy, Faculty of Mathematics and Natural Sciences, University of Oslo, 0316 Oslo, Norway. k.c.carlsson@farmasi.uio.no
Kristin Cecilie
KC
Carlsson
Nils Ove
NO
Hoem
Tracy
T
Glauser
Alexander A
AA
Vinks
D000293:Adolescent; D000328:Adult; D000927:Anticonvulsants / Q000097:blood / Q000493:pharmacokinetics* / Q000627:therapeutic use; D001499:Bayes Theorem; D002220:Carbamazepine / Q000097:blood / Q000493:pharmacokinetics* / Q000627:therapeutic use; D002648:Child; D002675:Child, Preschool; D003692:Delayed-Action Preparations; D016903:Drug Monitoring / Q000379:methods*; D004827:Epilepsy / Q000097:blood / Q000188:drug therapy*; D005260:Female; D006207:Half-Life; D006801:Humans; D008297:Male; D008657:Metabolic Clearance Rate; D008954:Models, Biological
D000927:Anticonvulsants; D003692:Delayed-Action Preparations; D002220:Carbamazepine