A Bayesian approach for estimating detection times in horses : exploring the pharmacokinetics of a urinary acepromazine metabolite

We describe the population pharmacokinetics of an acepromazine (ACP) metabolite (2-(1-hydroxyethyl)promazine) (HEPS) in horses for the estimation of likely detection times in plasma and urine. Acepromazine (30 mg) was administered to 12 horses, and blood and urine samples were taken at frequent intervals for chemical analysis. A Bayesian hierarchical model was fitted to describe concentration-time data and cumulative urine amounts for HEPS. The metabolite HEPS was modelled separately from the parent ACP as the half-life of the parent was considerably less than that of the metabolite. The clearance ($Cl/F_{PM}$) and volume of distribution ($V/F_{PM}$), scaled by the fraction of parent converted to metabolite, were estimated as 769 L/h and 6874 L, respectively. For a typical horse in the study, after receiving 30 mg of ACP, the upper limit of the detection time was 35 hours in plasma and 100 hours in urine, assuming an arbitrary limit of detection of 1 $\mu$g/L, and a small ($\approx 0.01$) probability of detection. The model derived allowed the probability of detection to be estimated at the population level. This analysis was conducted on data collected from only 12 horses, but we assume that this is representative of the wider population.

Population pharmacokinetic and pharmacogenetic analysis of 6-mercaptopurine in paediatric patients with acute lymphoblastic leukaemia

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT center dot The cytotoxic effects of 6-mercaptopurine (6-MP) were found to be due to drug-derived intracellular metabolites (mainly 6-thioguanine nucleotides and to some extent 6-methylmercaptopurine nucleotides) rather than the drug itself. center dot Current empirical dosing methods for oral 6-MP result in highly variable drug and metabolite concentrations and hence variability in treatment outcome. WHAT THIS STUDY ADDS center dot The first population pharmacokinetic model has been developed for 6-MP active metabolites in paediatric patients with acute lymphoblastic leukaemia and the potential demographic and genetically controlled factors that could lead to interpatient pharmacokinetic variability among this population have been assessed. center dot The model shows a large reduction in interindividual variability of pharmacokinetic parameters when body surface area and thiopurine methyltransferase polymorphism are incorporated into the model as covariates. center dot The developed model offers a more rational dosing approach for 6-MP than the traditional empirical method (based on body surface area) through combining it with pharmacogenetically guided dosing based on thiopurine methyltransferase genotype. To investigate the population pharmacokinetics of 6-mercaptopurine (6-MP) active metabolites in paediatric patients with acute lymphoblastic leukaemia (ALL) and examine the effects of various genetic polymorphisms on the disposition of these metabolites. Data were collected prospectively from 19 paediatric patients with ALL (n = 75 samples, 150 concentrations) who received 6-MP maintenance chemotherapy (titrated to a target dose of 75 mg m(-2) day(-1)). All patients were genotyped for polymorphisms in three enzymes involved in 6-MP metabolism. Population pharmacokinetic analysis was performed with the nonlinear mixed effects modelling program (nonmem) to determine the population mean parameter estimate of clearance for the active metabolites. The developed model revealed considerable interindividual variability (IIV) in the clearance of 6-MP active metabolites [6-thioguanine nucleotides (6-TGNs) and 6-methylmercaptopurine nucleotides (6-mMPNs)]. Body surface area explained a significant part of 6-TGNs clearance IIV when incorporated in the model (IIV reduced from 69.9 to 29.3%). The most influential covariate examined, however, was thiopurine methyltransferase (TPMT) genotype, which resulted in the greatest reduction in the model's objective function (P

Population pharmacokinetic model of canrenone after intravenous administration of potassium canrenoate to paediatric patients

Objectives: To characterize the population pharmacokinetics of canrenone following administration of potassium canrenoate in paediatric patients. Patients and Methods: Data were collected prospectively from 23 paediatric patients (2 days to 10 years of age; median weight 4 kg, range 2.16-28.0 kg) who received intravenous potassium canrenoate (K-canrenoate) as part of their intensive care therapy for removal of retained fluids e.g. in pulmonary oedema due to chronic lung disease and for the management of congestive heart failure. Plasma samples were analysed by HPLC for determination of canrenone (the major metabolite and pharmacologically active moiety) and the data subjected to pharmacokinetic analysis using NONMEM. Results: A one-compartment model best described the data. The only significant covariate was weight (WT). The final population models for canrenone clearance (CL/F) and volume of distribution (V/F) were CL/F (L/hr) = 11.4 × (WT /70.0)(0.75) and V/F (L) = 374.2 × (WT/70) where WT is in kg. The values of CL/F and V/F in a 4 kg child would be 1.33 L/hr and 21.4 L, respectively, resulting in an elimination half-life of 11.2 hr. Conclusions: The range of estimated CL/F in the study population was 0.67-7.38 L/hr. The data suggest that adjustment of K-canrenoate dosage according to body weight is appropriate in paediatric patients

Prophylactic ranitidine treatment in critically ill children - a population pharmacokinetic study

Aims: To characterize the population pharmacokinetics of ranitidine in critically ill children and to determine the influence of various clinical and demographic factors on its disposition. Methods: Data were collected prospectively from 78 paediatric patients (n = 248 plasma samples) who received oral or intravenous ranitidine for prophylaxis against stress ulcers, gastrointestinal bleeding or the treatment of gastro-oesophageal reflux. Plasma samples were analysed using high-performance liquid chromatography, and the data were subjected to population pharmacokinetic analysis using nonlinear mixed-effects modelling. Results: A one-compartment model best described the plasma concentration profile, with an exponential structure for interindividual errors and a proportional structure for intra-individual error. After backward stepwise elimination, the final model showed a significant decrease in objective function value (-12.618; P <0.001) compared with the weight-corrected base model. Final parameter estimates for the population were 32.1lh for total clearance and 285l for volume of distribution, both allometrically modelled for a 70kg adult. Final estimates for absorption rate constant and bioavailability were 1.31h and 27.5%, respectively. No significant relationship was found between age and weight-corrected ranitidine pharmacokinetic parameters in the final model, with the covariate for cardiac failure or surgery being shown to reduce clearance significantly by a factor of 0.46. Conclusions: Currently, ranitidine dose recommendations are based on children's weights. However, our findings suggest that a dosing scheme that takes into consideration both weight and cardiac failure/surgery would be more appropriate in order to avoid administration of higher or more frequent doses than necessary.

Potassium canrenoate treatment in paediatric patients: a population pharmacokinetic study using novel dried blood spot sampling

Objective: To characterize the population pharmacokinetics of canrenone following administration of potassium canrenoate (K-canrenoate) in paediatric patients.

Methods: Data were collected prospectively from 37 paediatric patients (median weight 2.9?kg, age range 2 days–0.85 years) who received intravenous K-canrenoate for management of retained fluids, for example in heart failure and chronic lung disease. Dried blood spot (DBS) samples (n?=?213) from these were analysed for canrenone content and the data subjected to pharmacokinetic analysis using nonlinear mixed-effects modelling. Another group of patients (n?=?16) who had 71 matching plasma and DBS samples was analysed separately to compare canrenone pharmacokinetic parameters obtained using the two different matrices.

Results: A one-compartment model best described the DBS data. Significant covariates were weight, postmenstrual age (PMA) and gestational age. The final population models for canrenone clearance (CL/F) and volume of distribution (V/F) in DBS were CL/F (l/h)?=?12.86?×? (WT/70.0)0.75?×?e [0.066?×? (PMA?-?40]) and V/F (l)?=?603.30?×? (WT/70)?×?(GA/40)1.89 where weight is in kilograms. The corresponding values of CL/F and V/F in a patient with a median weight of 2.9?kg are 1.11?l/h and 20.48?l, respectively. Estimated half-life of canrenone based on DBS concentrations was similar to that based on matched plasma concentrations (19.99 and 19.37?h, respectively, in 70?kg patient).

Conclusion: The range of estimated CL/F in DBS for the study population was 0.12–9.62?l/h; hence, bodyweight-based dosage adjustment of K-canrenoate appears necessary. However, a dosing scheme that takes into consideration both weight and age (PMA/gestational age) of paediatric patients seems more appropriate.

Pharmacokinetics of melatonin in preterm infants

Aims: Preterm infants are deprived of the normal intra-uterine exposure to maternal melatonin and may benefit from replacement therapy. We conducted a pharmacokinetic study to guide potential therapeutic trials. Methods: Melatonin was administered to 18 preterm infants in doses ranging from 0.04-0.6μgkg^-1 over 0.5-6h. Pharmacokinetic profiles were analyzed individually and by population methods. Results: Baseline melatonin was largely undetectable. Infants receiving melatonin at 0.1μgkg^-1h^-1 for 2h showed a median half-life of 15.82h and median maximum plasma concentration of 203.3pgml^-1. On population pharmacokinetics, clearance was 0.045lh^-1, volume of distribution 1.098l and elimination half-life 16.91h with gender (P = 0.047) and race (P < 0.0001) as significant covariates. Conclusions: A 2h infusion of 0.1μgkg^-1h^-1 increased blood melatonin from undetectable to approximately peak adult concentrations. Slow clearance makes replacement of a typical maternal circadian rhythm problematic. The pharmacokinetic profile of melatonin in preterm infants differs from that of adults so dosage of melatonin for preterm infants cannot be extrapolated from adult studies. Data from this study can be used to guide therapeutic clinical trials of melatonin in preterm infants. © 2013 The British Pharmacological Society.

Pantoprazole intraveineux aux soins intensifs pédiatriques: un modèle de pharmacocinétique de population

Objectifs : Définir les paramètres pharmacocinétiques du pantoprazole intraveineux en soins intensifs pédiatriques et déterminer l’influence qu’exercent sur ceux-ci les facteurs démographiques, le syndrome de réponse inflammatoire systémique (SRIS), la dysfonction hépatique et l’administration d’un inhibiteur du cytochrome (CYP) 2C19. Méthode : Cent cinquante-six concentrations plasmatiques de pantoprazole provenant d’une population de 20 patients (âgés de 10 jours à 16.4 ans) à risque ou atteints d’une hémorragie gastroduodénale de stress, ayant reçu des doses quotidiennes de pantoprazole de 19.9 à 140.6 mg/1.73m2, ont été analysées selon les méthodes non compartimentale et de modélisation non linéaire à effets mixtes. Résultats : Une clairance médiane (CL) de 0.14 L/h/kg, un volume apparent de distribution de 0.20 L/kg et une demi-vie d’élimination de 1.7 h ont été déterminés via l’approche non compartimentale. Le modèle populationnel à deux compartiments avec une infusion d’ordre zéro et une élimination d’ordre un représentait fidèlement la cinétique du pantoprazole. Le poids, le SRIS, la dysfonction hépatique et l’administration d’un inhibiteur du CYP2C19 constituaient les covariables significatives rendant compte de 75 % de la variabilité interindividuelle observée pour la CL. Seul le poids influençait significativement le volume central de distribution (Vc). Selon les estimations du modèle final, un enfant de cinq ans pesant 20 kg avait une CL de 5.28 L/h et un Vc de 2.22 L. La CL du pantoprazole augmentait selon l’âge et le poids tandis qu’elle diminuait respectivement de 62.3%, 65.8% et 50.5% en présence d’un SRIS, d’un inhibiteur du CYP2C19 ou d’une dysfonction hépatique. Conclusion : Ces résultats permettront de guider les cliniciens dans le choix d’une dose de charge et dans l’ajustement des posologies du pantoprazole en soins intensifs pédiatriques dépendamment de facteurs fréquemment rencontrés dans cette population.

Optimisation de l’administration des médicaments chez les enfants transplantés grâce à la pharmacocinétique de population

Ce travail de thèse porte sur l’application de la pharmacocinétique de population dans le but d’optimiser l’utilisation de certains médicaments chez les enfants immunosupprimés et subissant une greffe. Parmi les différents médicaments utilisés chez les enfants immunosupprimés, l’utilisation du busulfan, du tacrolimus et du voriconazole reste problématique, notamment à cause d’une très grande variabilité interindividuelle de leur pharmacocinétique rendant nécessaire l’individualisation des doses par le suivi thérapeutique pharmacologique. De plus, ces médicaments n’ont pas fait l’objet d’études chez les enfants et les doses sont adaptées à partir des adultes. Cette dernière pratique ne prend pas en compte les particularités pharmacologiques qui caractérisent l’enfant tout au long de son développement et rend illusoire l’extrapolation aux enfants des données acquises chez les adultes. Les travaux effectués dans le cadre de cette thèse ont étudié successivement la pharmacocinétique du busulfan, du voriconazole et du tacrolimus par une approche de population en une étape (modèles non-linéaires à effets mixtes). Ces modèles ont permis d’identifier les principales sources de variabilités interindividuelles sur les paramètres pharmacocinétiques. Les covariables identifiées sont la surface corporelle et le poids. Ces résultats confirment l’importance de tenir en compte l’effet de la croissance en pédiatrie. Ces paramètres ont été inclus de façon allométrique dans les modèles. Cette approche permet de séparer l’effet de la mesure anthropométrique d’autres covariables et permet la comparaison des paramètres pharmacocinétiques en pédiatrie avec ceux des adultes. La prise en compte de ces covariables explicatives devrait permettre d’améliorer la prise en charge a priori des patients. Ces modèles développés ont été évalués pour confirmer leur stabilité, leur performance de simulation et leur capacité à répondre aux objectifs initiaux de la modélisation. Dans le cas du busulfan, le modèle validé a été utilisé pour proposer par simulation une posologie qui améliorerait l’atteinte de l’exposition cible, diminuerait l’échec thérapeutique et les risques de toxicité. Le modèle développé pour le voriconazole, a permis de confirmer la grande variabilité interindividuelle dans sa pharmacocinétique chez les enfants immunosupprimés. Le nombre limité de patients n’a pas permis d’identifier des covariables expliquant cette variabilité. Sur la base du modèle de pharmacocinétique de population du tacrolimus, un estimateur Bayesien a été mis au point, qui est le premier dans cette population de transplantés hépatiques pédiatriques. Cet estimateur permet de prédire les paramètres pharmacocinétiques et l’exposition individuelle au tacrolimus sur la base d’un nombre limité de prélèvements. En conclusion, les travaux de cette thèse ont permis d’appliquer la pharmacocinétique de population en pédiatrie pour explorer les caractéristiques propres à cette population, de décrire la variabilité pharmacocinétique des médicaments utilisés chez les enfants immunosupprimés, en vue de l’individualisation du traitement. Les outils pharmacocinétiques développés s’inscrivent dans une démarche visant à diminuer le taux d'échec thérapeutique et l’incidence des effets indésirables ou toxiques chez les enfants immunosupprimés suite à une transplantation.

Quantitative justification for target concentration intervention - parameter variability and predictive performance using population pharmacokinetic models for aminoglycosides

Aims [1] To quantify the random and predictable components of variability for aminoglycoside clearance and volume of distribution [2] To investigate models for predicting aminoglycoside clearance in patients with low serum creatinine concentrations [3] To evaluate the predictive performance of initial dosing strategies for achieving an aminoglycoside target concentration. Methods Aminoglycoside demographic, dosing and concentration data were collected from 697 adult patients (>=20 years old) as part of standard clinical care using a target concentration intervention approach for dose individualization. It was assumed that aminoglycoside clearance had a renal and a nonrenal component, with the renal component being linearly related to predicted creatinine clearance. Results A two compartment pharmacokinetic model best described the aminoglycoside data. The addition of weight, age, sex and serum creatinine as covariates reduced the random component of between subject variability (BSVR) in clearance (CL) from 94% to 36% of population parameter variability (PPV). The final pharmacokinetic parameter estimates for the model with the best predictive performance were: CL, 4.7 l h(-1) 70 kg(-1); intercompartmental clearance (CLic), 1 l h(-1) 70 kg(-1); volume of central compartment (V-1), 19.5 l 70 kg(-1); volume of peripheral compartment (V-2) 11.2 l 70 kg(-1). Conclusions Using a fixed dose of aminoglycoside will achieve 35% of typical patients within 80-125% of a required dose. Covariate guided predictions increase this up to 61%. However, because we have shown that random within subject variability (WSVR) in clearance is less than safe and effective variability (SEV), target concentration intervention can potentially achieve safe and effective doses in 90% of patients.

Some considerations on the design of population pharmacokinetic studies

The goal of this manuscript is to introduce a framework for consideration of designs for population pharmacokinetic orpharmacokinetic-pharmacodynamic studies. A standard one compartment pharmacokinetic model with first-order input and elimination is considered. A series of theoretical designs are considered that explore the influence of optimizing the allocation of sampling times, allocating patients to elementary designs, consideration of sparse sampling and unbalanced designs and also the influence of single vs. multiple dose designs. It was found that what appears to be relatively sparse sampling (less blood samples per patient than the number of fixed effects parameters to estimate) can also be highly informative. Overall, it is evident that exploring the population design space can yield many parsimonious designs that are efficient for parameter estimation and that may not otherwise have been considered without the aid of optimal design theory.

Pharmacokinetics of melatonin in preterm infants

Aims - Preterm infants are deprived of the normal intra-uterine exposure to maternal melatonin and may benefit from replacement therapy. We conducted a pharmacokinetic study to guide potential therapeutic trials. Methods - Melatonin was administered to 18 preterm infants in doses ranging from 0.04–0.6 μg kg−1 over 0.5–6 h. Pharmacokinetic profiles were analyzed individually and by population methods. Results - Baseline melatonin was largely undetectable. Infants receiving melatonin at 0.1 μg kg−1 h−1 for 2 h showed a median half-life of 15.82 h and median maximum plasma concentration of 203.3 pg ml−1. On population pharmacokinetics, clearance was 0.045 l h−1, volume of distribution 1.098 l and elimination half-life 16.91 h with gender (P = 0.047) and race (P < 0.0001) as significant covariates. Conclusions - A 2 h infusion of 0.1 μg kg−1 h−1 increased blood melatonin from undetectable to approximately peak adult concentrations. Slow clearance makes replacement of a typical maternal circadian rhythm problematic. The pharmacokinetic profile of melatonin in preterm infants differs from that of adults so dosage of melatonin for preterm infants cannot be extrapolated from adult studies. Data from this study can be used to guide therapeutic clinical trials of melatonin in preterm infants.

Prophylactic ranitidine treatment in critically ill children - a population pharmacokinetic study

Aims - To characterize the population pharmacokinetics of ranitidine in critically ill children and to determine the influence of various clinical and demographic factors on its disposition. Methods - Data were collected prospectively from 78 paediatric patients (n = 248 plasma samples) who received oral or intravenous ranitidine for prophylaxis against stress ulcers, gastrointestinal bleeding or the treatment of gastro-oesophageal reflux. Plasma samples were analysed using high-performance liquid chromatography, and the data were subjected to population pharmacokinetic analysis using nonlinear mixed-effects modelling. Results - A one-compartment model best described the plasma concentration profile, with an exponential structure for interindividual errors and a proportional structure for intra-individual error. After backward stepwise elimination, the final model showed a significant decrease in objective function value (−12.618; P < 0.001) compared with the weight-corrected base model. Final parameter estimates for the population were 32.1 l h−1 for total clearance and 285 l for volume of distribution, both allometrically modelled for a 70 kg adult. Final estimates for absorption rate constant and bioavailability were 1.31 h−1 and 27.5%, respectively. No significant relationship was found between age and weight-corrected ranitidine pharmacokinetic parameters in the final model, with the covariate for cardiac failure or surgery being shown to reduce clearance significantly by a factor of 0.46. Conclusions - Currently, ranitidine dose recommendations are based on children's weights. However, our findings suggest that a dosing scheme that takes into consideration both weight and cardiac failure/surgery would be more appropriate in order to avoid administration of higher or more frequent doses than necessary.

Population pharmacokinetic and pharmacogenetic analysis of 6-mercaptopurine in paediatric patients with acute lymphoblastic leukaemia

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT • The cytotoxic effects of 6-mercaptopurine (6-MP) were found to be due to drug-derived intracellular metabolites (mainly 6-thioguanine nucleotides and to some extent 6-methylmercaptopurine nucleotides) rather than the drug itself. • Current empirical dosing methods for oral 6-MP result in highly variable drug and metabolite concentrations and hence variability in treatment outcome. WHAT THIS STUDY ADDS • The first population pharmacokinetic model has been developed for 6-MP active metabolites in paediatric patients with acute lymphoblastic leukaemia and the potential demographic and genetically controlled factors that could lead to interpatient pharmacokinetic variability among this population have been assessed. • The model shows a large reduction in interindividual variability of pharmacokinetic parameters when body surface area and thiopurine methyltransferase polymorphism are incorporated into the model as covariates. • The developed model offers a more rational dosing approach for 6-MP than the traditional empirical method (based on body surface area) through combining it with pharmacogenetically guided dosing based on thiopurine methyltransferase genotype. AIMS - To investigate the population pharmacokinetics of 6-mercaptopurine (6-MP) active metabolites in paediatric patients with acute lymphoblastic leukaemia (ALL) and examine the effects of various genetic polymorphisms on the disposition of these metabolites. METHODS - Data were collected prospectively from 19 paediatric patients with ALL (n = 75 samples, 150 concentrations) who received 6-MP maintenance chemotherapy (titrated to a target dose of 75 mg m−2 day−1). All patients were genotyped for polymorphisms in three enzymes involved in 6-MP metabolism. Population pharmacokinetic analysis was performed with the nonlinear mixed effects modelling program (nonmem) to determine the population mean parameter estimate of clearance for the active metabolites. RESULTS - The developed model revealed considerable interindividual variability (IIV) in the clearance of 6-MP active metabolites [6-thioguanine nucleotides (6-TGNs) and 6-methylmercaptopurine nucleotides (6-mMPNs)]. Body surface area explained a significant part of 6-TGNs clearance IIV when incorporated in the model (IIV reduced from 69.9 to 29.3%). The most influential covariate examined, however, was thiopurine methyltransferase (TPMT) genotype, which resulted in the greatest reduction in the model's objective function (P < 0.005) when incorporated as a covariate affecting the fractional metabolic transformation of 6-MP into 6-TGNs. The other genetic covariates tested were not statistically significant and therefore were not included in the final model. CONCLUSIONS - The developed pharmacokinetic model (if successful at external validation) would offer a more rational dosing approach for 6-MP than the traditional empirical method since it combines the current practice of using body surface area in 6-MP dosing with a pharmacogenetically guided dosing based on TPMT genotype.

Population pharmacokinetic model of canrenone after intravenous administration of potassium canrenoate to paediatric patients

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT • Little is known about the pharmacokinetics of potassium canrenoate/canrenone in paediatric patients WHAT THIS STUDY ADDS • A population pharmacokinetic model has been developed to evaluate the pharmacokinetics of canrenone in paediatric patients who received potassium canrenoate as part of their therapy in the intensive care unit. AIMS To characterize the population pharmacokinetics of canrenone following administration of potassium canrenoate to paediatric patients. METHODS Data were collected prospectively from 23 paediatric patients (2 days to 10 years of age; median weight 4 kg, range 2.16–28.0 kg) who received intravenous potassium canrenoate (K-canrenoate) as part of their intensive care therapy for removal of retained fluids, e.g. in pulmonary oedema due to chronic lung disease and for the management of congestive heart failure. Plasma samples were analyzed by HPLC for determination of canrenone (the major metabolite and pharmacologically active moiety) and the data subjected to pharmacokinetic analysis using NONMEM. RESULTS A one compartment model best described the data. The only significant covariate was weight (WT). The final population models for canrenone clearance (CL/F) and volume of distribution (V/F) were CL/F (l h−1) = 11.4 × (WT/70.0)0.75 and V/F (l) = 374.2 × (WT/70) where WT is in kg. The values of CL/F and V/F in a 4 kg child would be 1.33 l h−1 and 21.4 l, respectively, resulting in an elimination half-life of 11.2 h. CONCLUSIONS The range of estimated CL/F in the study population was 0.67–7.38 l h−1. The data suggest that adjustment of K-canrenoate dosage according to body weight is appropriate in paediatric patients.