Comparison of nonlinear mixed-effect and non-parametric expectation maximisation modelling for Bayesian estimation of carboplatin clearance in children

Parametric and Nonparametric Population Pharmacokinetic Models to Assess Probability of Target Attainment of Imipenem Concentrations in Critically Ill Patients

Population pharmacokinetic modeling and simulation (M&S) are used to improve antibiotic dosing. Little is known about the differences in parametric and nonparametric M&S. Our objectives were to compare (1) the external validation of parametric and nonparametric models of imipenem in critically ill patients and (2) the probability of target attainment (PTA) calculations using simulations of both models. The M&S software used was NONMEM 7.2 (parametric) and Pmetrics 1.5.2 (nonparametric). The external predictive performance of both models was adequate for eGFRs ≥ 78 mL/min but insufficient for lower eGFRs, indicating that the models (developed using a population with eGFR ≥ 60 mL/min) could not be extrapolated to lower eGFRs. Simulations were performed for three dosing regimens and three eGFRs (90, 120, 150 mL/min). Fifty percent of the PTA results were similar for both models, while for the other 50% the nonparametric model resulted in lower MICs. This was explained by a higher estimated between-subject variability of the nonparametric model. Simulations indicated that 1000 mg q6h is suitable to reach MICs of 2 mg/L for eGFRs of 90-120 mL/min. For MICs of 4 mg/L and for higher eGFRs, dosing recommendations are missing due to largely different PTA values per model. The consequences of the different modeling approaches in clinical practice should be further investigated.

Diversity of dose-individualization and therapeutic drug monitoring practices of platinum compounds: a review

Introduction: Platinum-derived drugs are commonly used for the treatment of solid tumors. The differences in chemical structures of these molecules lead to different pharmacological properties, in terms of indication, efficacy, and toxicity. Their pharmacokinetics (PK) differ according to their respective renal elimination and have led to many studies investigating their dose optimization. Area covered: This review attempts to summarize and compare PK and pharmacodynamics of cisplatin, carboplatin, and oxaliplatin, with an emphasis on differences of dose calculations and opportunities for therapeutic drug monitoring (TDM) in various patient populations. Expert opinion: Although cisplatin and carboplatin can be considered as analogs since they share the same DNA interacting properties, the slower hydrolysis of the latter results in a better safety profile. Carboplatin is the only drug in oncology to be administrated according to a target area under the curve of concentration versus time, considering that its PK variability is almost fully explained by renal function, not by body size. This enables individual dosing based on predicted carboplatin clearance (along with patients renal characteristics) or on actual clearance with TDM, especially in a high-dose protocol.

Population Pharmacokinetics of Imipenem in Critically Ill Patients: A Parametric and Nonparametric Model Converge on CKD-EPI Estimated Glomerular Filtration Rate as an Impactful Covariate

BACKGROUND

Population pharmacokinetic (popPK) models for antibiotics are used to improve dosing strategies and individualize dosing by therapeutic drug monitoring. Little is known about the differences in results of parametric versus nonparametric popPK models and their potential consequences in clinical practice. We developed both parametric and nonparametric models of imipenem using data from critically ill patients and compared their results.

METHODS

Twenty-six critically ill patients treated with intravenous imipenem/cilastatin were included in this study. Median estimated glomerular filtration rate (eGFR) measured by the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation was 116 mL/min/1.73 m2 (interquartile range 104-124) at inclusion. The usual dosing regimen was 500 mg/500 mg four times daily. On average, five imipenem levels per patient (138 levels in total) were drawn as peak, intermediate, and trough levels. Imipenem concentration-time profiles were analyzed using parametric (NONMEM 7.2) and nonparametric (Pmetrics 1.5.2) popPK software.

RESULTS

For both methods, data were best described by a model with two distribution compartments and the CKD-EPI eGFR equation unadjusted for body surface area as a covariate on the elimination rate constant (Ke). The parametric population parameter estimates were Ke 0.637 h-1 (between-subject variability [BSV]: 19.0% coefficient of variation [CV]) and central distribution volume (Vc) 29.6 L (without BSV). The nonparametric values were Ke 0.681 h-1 (34.0% CV) and Vc 31.1 L (42.6% CV).

CONCLUSIONS

Both models described imipenem popPK well; the parameter estimates were comparable and the included covariate was identical. However, estimated BSV was higher in the nonparametric model. This may have consequences for estimated exposure during dosing simulations and should be further investigated in simulation studies.

A limited sampling strategy to estimate individual pharmacokinetic parameters of methotrexate in children with acute lymphoblastic leukemia

PURPOSE

The objectives of this study were to characterize the population pharmacokinetics of MTX in patients with acute lymphoblastic leukemia (ALL) with ages ranging from 2 to 16 years and to propose a limited sampling strategy to estimate individual pharmacokinetic parameters.

METHODS

Seventy-nine children were enrolled in this study; they received 1-4 courses of chemotherapy. MTX was administered at a dose of 5 g/m2. MTX population parameters were estimated from 61 patients (231 courses; age range: 2-16 years). The data were analyzed by nonlinear mixed-effect modeling with use of a two-compartment structural model. The interoccasion variability was taken into account in the model. Eighteen additional patients (70 courses) were used to evaluate the predictive performances of the Bayesian approach and to devise a limited sampling strategy.

RESULTS

The following population parameters were obtained: total clearance (CL) = 8.8 l/h (inter-individual variability: 43%), initial volume of distribution (V1) = 17.3 l (48%), k12 = 0.0225 h(-1) (41%), and k21 = 0.0629 h(-1) (24%). The inter-individual variability in the initial volume of distribution was partially explained by the fact that this parameter was weight-dependent. Intercourse variability was limited, with a mean variation of 13.2%. The protocol involving two sampling times, 24 and 48 h after the beginning of infusion, allows precise and accurate determination of individual pharmacokinetic parameters and consequently, it was possible to predict the time at which the MTX concentration reached the predicted threshold (0.2 microM) below which the administration of folinic acid could be stopped.

CONCLUSION

The results of this study combine the relationships between the pharmacokinetic parameters of MTX and patient covariates that may be useful for dose adjustment, with a convenient sampling procedure that may aid in optimizing pediatric patient care.

Modelling and clinical trials in paediatrics

Clinical trials are more difficult to conduct in children, but they are even more necessary than in adults their scarcity is an ethical scandal. Mathematical models can be built that can describe both the disease process and the mechanism of action of drugs. These models can then be used to simulate the outcome of clinical trials. Inspection of the simulated results then facilitates optimisation of the trial design and proposed methods of analysis. Validation is a crucial issue for the good practice of modelling and simulation. The participants of Round Table No. 6 recommend: (i) that modelling be systematically employed; (ii) that all the required professional personnel be involved, at all phases; (iii) that all data needed are made accessible; (iv) that clinicians be trained; (v) that specialists develop training tool kits; and (vi) that universities provide appropriate training.

Population pharmacokinetics of ciprofloxacin in pediatric and adolescent patients with acute infections

The aim of the present study was to characterize the population pharmacokinetics of ciprofloxacin in patients with and without cystic fibrosis ranging in age from 1 day to 24 years and to propose a limited sampling strategy to estimate individual pharmacokinetic parameters. Patients were divided into four groups according to the treatment schedule. They received ciprofloxacin by intravenous infusion (30 min) or by the oral route. The number of samples collected from each patient ranged from 1 to 12. The population parameters were computed for an initial group of 37 patients. The data were analyzed by nonlinear mixed-effect modeling by use of a two-compartment structural model. The interindividual variability in clearance (CL) was partially explained by a dependence on age and the patient's clinical status. In addition, a significant relationship was found between weight and the initial volume of distribution. Eighteen additional patients were used for model validation and evaluation of limited sampling strategies. When ciprofloxacin was administered intravenously, sampling at a single point (12 h after the start of infusion) allowed the precise and accurate estimation of CL and the elimination half-life, as well as the ciprofloxacin concentration at the end of the infusion. It should be noted that to take into account the presence of a lag time after oral administration, a schedule based on two sampling times of 1 and 12 h is needed. The results of this study combine relationships between ciprofloxacin pharmacokinetic parameters and patient covariates that may be useful for dose adjustment and a convenient sampling procedure that can be used for further studies.

Pharmacology of cytotoxic agents: a helpful tool for building dose adjustment guidelines in the elderly

Aging is associated with multidimensional changes, including alterations in physiological functions, co-morbidities and poly-medications. These changes may lead to modifications in the absorption, distribution, metabolism and excretion of drugs. The lack of a scientific basis for optimal drug dosing in the elderly is a major problem. The development and validation of guidelines are therefore essential to improve treatment administration and monitoring in elderly patients. Even though it has been widely demonstrated that standard therapies used in adults may be of great benefit in the elderly, there may be a higher incidence of toxicity. This could be avoided by using dosage individualization based on a sound knowledge of the physiological factors implicated in the pharmacokinetic (PK) characteristics of the drugs administered and in their observed pharmacodynamic (PD) effects in each patient. The so-called "population modeling" approach renders such studies feasible by allowing the analysis of PK-PD relationships from sparse observational data.

Individual dosing of carboplatin based on drug monitoring in children receiving high-dose chemotherapy

Individual dosing of carboplatin based on drug monitoring was performed within a multi-centric phase I study based on high AUC-levels in children. Twelve patients (aged 3-17 years old) have been included: 3, 5, and 4 patients at the overall target ultrafilterable carboplatin AUC of 20, 25, or 30 mg/ml x min, respectively. Carboplatin was administered as a daily 60-min infusion, repeated on five consecutive days. The initial daily dose corresponding to the three first days was calculated according to the carboplatin clearance (CL) predicted from patients' characteristics (body weight, serum creatinine and nephrectomy status). Three blood samples were taken per patient. The individual CL were estimated by MAP (maximum a posteriori approach) Bayesian method implemented in the MP-K program. The doses for day 4 and 5 was adjusted in order to obtain the overall target AUC. Drug monitoring led to a change in the carboplatin dose (overall administered dose versus overall dose planned) ranging from -41% to +45%. Pharmacokinetics were performed at day 5 for 7/12 children: mean relative change between day 1 and day 5 was -11% showing a statistically significant, but limited, decrease of CL from day 1 to day 5. The percentage of difference between the observed and target overall AUC ranged between -7% and +14%. Three patients (one at each AUC level) who were previously treated with cisplatin experienced dose-limiting hearing loss. In conclusion, drug monitoring and dose adjustment is needed for the control of carboplatin plasma exposure when administering high doses of carboplatin in children.

NONMEM and NPEM2 population modeling: a comparison using tobramycin data in neonates

Nonlinear mixed effects modeling (NONMEM) and nonparametric expectation maximization (NPEM2) have both been used in population modeling of tobramycin. We compared both methods for differences in population pharmacokinetic parameters in relation to error models used. Predictive performance was compared between models. A group of 470 neonates who had received tobramycin according to a gestational age (GA)-dependent dosing interval was analyzed according to a one-compartment model with NONMEM and NPEM2. Additional models were constructed where the assay error pattern in NPEM2 mimicked NONMEM residual error and vice versa. Individual pharmacokinetic parameter estimates were compared. Predictive performance was evaluated in a separate group of 61 patients. Population estimates and variation coefficients (CV) for optimal models were NONMEM K(el) 0.071 h(-1) (27%), V(d) 0.59 L/kg (9%); NPEM2 K(el) 0.079 h(-1) (42%), V(d) 0.65 L/kg (48%). Forcing NONMEM to use the NPEM2 error pattern as residual error or vice versa resulted in smaller differences in CVs of the estimates. NONMEM gave less bias (P < 0.05) than NPEM2 and comparable precision with this approach. In conclusion NONMEM and NPEM2 are dissimilar in population estimates. Differences in ranges of pharmacokinetic parameter estimates between NONMEM and NPEM2 are largely determined by the method of incorporating error patterns in both programs.