Computing normalised prediction distribution errors to evaluate nonlinear mixed-effect models: the npde add-on package for R

Population Pharmacokinetics of Palbociclib in aReal-World Situation

Palbociclib is an oral cyclin-dependent kinase inhibitor that is used in combination with aromatase inhibitors in the treatment of postmenopausal women with metastatic breast cancer. Its metabolism profile is associated with an important interpatient variability. We performed a population pharmacokinetics study of palbociclib in women routinely followed in a cancer center. One hundred and fifty-one samples were analyzed. The sampling times after administration ranged from 0.9 to 75 h and the samples were taken between 1 and 21 days after the beginning of the palbociclib cycle. Palbociclib was determined using a validated mass spectrometry method. The best model that described the concentrations was a one-compartment model with first-order absorption and an absorption lag time. Interindividual variability could only be estimated on the clearance and the first-order absorption. Creatinine clearance was found to be a significant covariate for the apparent clearance. No significant covariates could be observed with the first-order absorption. First-order absorption and absorption lag times were difficult to assess because of the constraints linked to the real-world setting due to the small number of samples used during the absorption process. However, palbociclib apparent clearance was satisfactorily estimated. Population pharmacokinetics (POP PK) with palbociclib could help to optimize dosing.

Pharmacokinetics and safety of pegylated recombinant human granulocyte colony-stimulating factor in children with acute leukaemia

AIMS

This open-label, phase I study evaluated the pharmacokinetics and safety of pegylated recombinant human granulocyte colony-stimulating factor (PEG-rhG-CSF) for the treatment of chemotherapy-induced neutropenia in children with acute leukaemia.

METHODS

PEG-rhG-CSF was administered as a single 100 mcg/kg (3 mg maximum dose) subcutaneous injection at the end of each chemotherapy period when neutropenia occurred. Blood samples were obtained from patients treated with PEG-rhG-CSF. PEG-rhG-CSF serum concentrations were determined by an enzyme-linked immunosorbent assay. Population pharmacokinetic (PPK) analysis was implemented using the nonlinear mixed-effects model. Short-term safety was evaluated through adverse events collection (registered at clinicaltrials.gov identifier: 03844360).

RESULTS

A total of 16 acute leukaemia patients (1.8-13.6 years) were included, of whom two (12.5%) had grade 3 neutropenia, six (37.5%) had grade 4 neutropenia, and eight (50.0%) had severe neutropenia. For PPK modelling, 64 PEG-rhG-CSF serum concentrations were obtainable. A one-compartment model with first-order elimination was used for pharmacokinetic data modelling. The current weight was a significant covariate. The median (range) of clearance (CL) and area under the serum concentration-time curve (AUC) were 5.65 (1.49-14.45) mL/h/kg and 16514.75 (6632.45-54423.30) ng·h/mL, respectively. Bone pain, pyrexia, anaphylaxis and nephrotoxicity were not observed. One patient died 13 days after administration, and the objective assessment of causality was that an association with PEG-rhG-CSF was "possible".

CONCLUSIONS

The AUC of PEG-rhG-CSF (100 mcg/kg, 3 mg maximum dose) in paediatric patients with acute leukaemia were similar to those of PEG-rhG-CSF (100 mcg/kg) in children with sarcoma. PEG-rhG-CSF is safe, representing an important therapeutic option for chemotherapy-induced neutropenia in paediatric patients with acute leukaemia.

Determination of vancomycin exposure target and individualised dosing recommendations for neonates: model-informed precision dosing

INTRODUCTION

Few studies incorporating population pharmacokinetic/pharmacodynamic (Pop-PK/PD) modelling have been conducted to quantify the exposure target of vancomycin in neonates. A retrospective observational cohort study was undertaken in neonates to determine this target and dosing recommendations (chictr.org.cn, ChiCTR1900027919).

METHODS

A Pop-PK model was developed to estimate PK parameters. Causalities between acute kidney injury (AKI) occurrence and vancomycin use were verified using Naranjo criteria. Thresholds of vancomycin exposure in predicting AKI or efficacy were identified via classification and regression tree analysis. Associations between exposure thresholds and clinical outcomes, including AKI and efficacy, were analysed by logistic regression. Dosing recommendations were designed using Monte Carlo simulations based on the optimised exposure target.

RESULTS

Pop-PK modelling included 182 neonates with 411 observations. On covariate analysis, neonatal physiological maturation, renal function and concomitant use of vasoactive agents (VAS) significantly affected vancomycin PK. Seven cases of vancomycin-induced AKI were detected. Area under the concentration-time curve from 0-24 hours (AUC_0-24) ≥ 485 mg•h/L was an independent risk factor for AKI after adjusting for VAS co-administration. The clinical efficacy of vancomycin was analysed in 42 patients with blood culture-proven staphylococcal sepsis. AUC_0-24 to minimum inhibitory concentration (AUC_0-24/MIC) ≥ 234 was the only significant predictor of clinical effectiveness. Monte Carlo simulations indicated that regimens in Neonatal Formulary 7 and Red Book (2018) were unsuitable for all neonates.

CONCLUSION

An AUC_0-24 of 240-480 (assuming MIC = 1 mg/L) is a recommended exposure target of vancomycin in neonates. Model-informed dosing regimens are valuable in clinical practice.

Comparison of Precision and Accuracy of Five Methods to Analyse Total Score Data

Total score (TS) data is generated from composite scales consisting of several questions/items, such as the Movement Disorder Society-Unified Parkinson's Disease Rating Scale (MDS-UPDRS). The analysis method that most fully uses the information gathered is item response theory (IRT) models, but these are complex and require item-level data which may not be available. Therefore, the TS is commonly analysed with standard continuous variable (CV) models, which do not respect the bounded nature of data. Bounded integer (BI) models do respect the data nature but are not as extensively researched. Mixed models for repeated measures (MMRM) are an alternative that requires few assumptions and handles dropout without bias. If an IRT model exists, the expected mean and standard deviation of TS can be computed through IRT-informed functions-which allows CV and BI models to estimate parameters on the IRT scale. The fit, performance on external data and parameter precision (when applicable) of CV, BI and MMRM to analyse simulated TS data from the MDS-UPDRS motor subscale are investigated in this work. All models provided accurate predictions and residuals without trends, but the fit of CV and BI models was improved by IRT-informed functions. The IRT-informed BI model had more precise parameter estimates than the IRT-informed CV model. The IRT-informed models also had the best performance on external data, while the MMRM model was worst. In conclusion, (1) IRT-informed functions improve TS analyses and (2) IRT-informed BI models had more precise IRT parameter estimates than IRT-informed CV models.

Theophylline dosing and pharmacokinetics for renal protection in neonates with hypoxic-ischemic encephalopathy undergoing therapeutic hypothermia

BACKGROUND

Theophylline, a non-selective adenosine receptor antagonist, improves renal perfusion in the setting of hypoxia-ischemia and may offer therapeutic benefit in neonates with hypoxic-ischemic encephalopathy (HIE) undergoing hypothermia. We evaluated the pharmacokinetics and dose-exposure relationships of theophylline in this population to guide dosing strategies.

METHODS

A population pharmacokinetic analysis was performed in 22 neonates with HIE undergoing hypothermia who were part of a prospective study or retrospective chart review. Aminophylline (intravenous salt form of theophylline) was given per institutional standard of care for low urine output and/or rising serum creatinine (5 mg/kg intravenous (i.v.) load then 1.8 mg/kg i.v. q6h). The ability of different dosing regimens to achieve target concentrations (4-10 mg/L) associated with clinical response was examined.

RESULTS

Birth weight was a significant predictor of theophylline clearance and volume of distribution (p < 0.05). The median half-life was 39.5 h (range 27.2-50.4). An aminophylline loading dose of 7 mg/kg followed by 1.6 mg/kg q12h was predicted to achieve target concentrations in 84% of simulated neonates.

CONCLUSIONS

In neonates with HIE undergoing hypothermia, theophylline clearance was low with a 50% longer half-life compared to full-term normothermic neonates without HIE. Dosing strategies need to consider the unique pharmacokinetic needs of this population.

IMPACT

Theophylline is a potential renal-protective therapy in neonates with HIE undergoing therapeutic hypothermia; however, the pharmacokinetics and dose needs in this population are not known. Theophylline clearance was low in neonates with HIE undergoing therapeutic hypothermia with a 50% longer half-life compared to full-term normothermic neonates without HIE. As theophylline is advanced in clinical development, dosing strategies will need to consider the unique pharmacokinetic needs of neonates with HIE undergoing therapeutic hypothermia.

Population pharmacokinetics of azathioprine active metabolite in patients with inflammatory bowel disease and dosage regimens optimisation

Azathioprine is a first-line drug used to maintain the remission of inflammatory bowel disease (IBD). As a prodrug, azathioprine is metabolised to produce active 6-thioguanine nucleotides (6-TGN). There are large individual variations in the pharmacokinetics/pharmacodynamics of 6-TGN in patients with IBD. Here, we aimed to develop a model to quantitatively investigate factors that affect 6-TGN pharmacokinetics to formulate a dosage guideline for azathioprine. Data were collected prospectively from 100 adult patients with IBD who were receiving azathioprine. Patients were genotyped for two single-nucleotide polymorphisms (TPMT*3C c.719A > G and NUDT15 c.415C > T). Using high-performance liquid chromatography, we measured 156 steady-state trough concentrations of 6-TGN within the range 0.09 to 1.16 mg/L (ie 133-1733 pmol per 8 × 10⁸ RBC). The covariates analysed included sex, age, body-weight, laboratory tests and concomitant medications. A population pharmacokinetic model was established using "non-linear mixed-effects modelling" software and the "first-order conditional estimation method with interaction." Body-weight, TPMT*3C polymorphisms and co-therapy with mesalazine were found to be important factors influencing the clearance of 6-TGN. A dosage guideline for azathioprine was developed based on the PPK model that enables individualised azathioprine dosing in adult patients with different body-weights, TPMT*3C genotypes and co-administration with mesalazine.

Population pharmacokinetics and dosing optimization of metformin in Chinese patients with type 2 diabetes mellitus

Approximately 35% of patients fail to attain ideal initial blood glucose control under metformin monotherapy. The objective of this observational study is to simulate the optimal protocol of metformin according to the different renal function.The population pharmacokinetics of metformin was performed in 125 subjects with type 2 diabetes mellitus. Plasma concentrations of metformin were quantified by high-performance liquid chromatography. A population pharmacokinetic model of metformin was developed using NONMEN (version 7.2, Icon Development Solutions, USA). Monte Carlo simulation was used to simulate the concentration-time profiles for doses of metformin for 1000 times at different stages of renal function.The mean population pharmacokinetic parameters were apparent clearance 53.0 L/h, apparent volume of distribution 438 L, absorption rate constant 1.4 hour and lag-time 0.91 hour. Covariate analyses revealed that estimated glomerular filtration rate (eGFR) and bodyweight as individual factors influencing the apparent oral clearance: CL/F = 53.0 × ( bodyweight/75) × (eGFR/102.5)EXP(0.1797). The results of the simulation showed that patients should be prescribed metformin 2550 mg/d (t.i.d.) vs 3000 mg/d (b.i.d.) as the minimum doses for patients with augmented renal clearance.eGFR had a significant impact on metformin pharmacokinetics. Patients administered metformin twice a day require higher total daily doses than those with a regimen of 3 times a day at each stage of kidney function.

Non-Linear Pharmacokinetics of Oral Roscovitine (Seliciclib) in Cystic Fibrosis Patients Chronically Infected with Pseudomonas aeruginosa: A Study on Population Pharmacokinetics with Monte Carlo Simulations

Roscovitine (Seliciclib), a new protein kinase inhibitor, was administered orally to adult patients with cystic fibrosis for the first time in the ROSCO-CF trial, a dose-escalation, phase IIa, randomized, controlled trial. Extensive pharmacokinetic sampling was performed up to 12 h after the first oral dose. Roscovitine and its main metabolite M3 were quantified by liquid chromatography coupled with tandem mass spectrometry. The pharmacokinetics analyses were performed by non-linear mixed effects modelling. Monte Carlo simulations were performed to assess the impact of dose on the pharmacokinetics of oral roscovitine. Twenty-three patients received oral doses ranging from 200 to 800 mg of roscovitine and 138 data points were available for both roscovitine and M3 concentrations. The pharmacokinetics was best described by a two-compartment parent-metabolite model, with a complex saturable absorption process modelled as the sum of Gaussian inverse density functions. The Monte Carlo simulations showed a dose-dependent and saturable first-pass effect leading to pre-systemic formation of M3. The treatment with proton-pump inhibitors reduced the rate of absorption of oral roscovitine. The pharmacokinetics of oral roscovitine in adult patients with cystic fibrosis was non-linear and showed significant inter-individual variability. A repeat-dose study will be required to assess the inter-occasional variability of its pharmacokinetics.

Severity parameters for asphyxia or hypoxic-ischemic encephalopathy do not explain inter-individual variability in the pharmacokinetics of phenobarbital in newborns treated with therapeutic hypothermia

BACKGROUND

The current study uses a population modeling approach to evaluate and quantify the impact of severity of asphyxia and hypoxic-ischemic encephalopathy (HIE) on the pharmacokinetics of phenobarbital in asphyxiated newborns treated with therapeutic hypothermia.

METHODS

Included newborns received phenobarbital (the TOBY trial protocol). 120 plasma samples were available from 50 newborns, median (IQR) weight 3.3 (2.8-3.5) kg and gestational age 39 (39-40) weeks. NONMEM^® version 7.2 was used for the data analysis. Age, body weight, sex, concomitant medications, kidney and liver function markers, as well as severity parameters of asphyxia and HIE were tested as potential covariates of pharmacokinetics of phenobarbital. Severe asphyxia was defined as pH of arterial umbilical cord blood ≤7.1 and Apgar 5 ≤5, and severe HIE was defined as time to normalization of amplitude-integrated electroencephalography (aEEG) >24 h.

RESULTS

Weight was found to be the only statistically significant covariate for the volume of distribution. At weight of 1 kg volume of distribution was 0.91 L and for every additional kg it increased in 0.91 L. Clearance was 0.00563 L/h. No covariates were statistically significant for the clearance of phenobarbital.

CONCLUSIONS

Phenobarbital dose adjustments are not indicated in the studied population, irrespective of the severity of asphyxia or HIE.

Impact of enantiomer-specific changes in pharmacokinetics between infants and adults on the target concentration of racemic ketorolac: A pooled analysis

Aims

Ketorolac is a nonsteroidal anti‐inflammatory racemic drug with analgesic effects only attributed to its S‐enantiomer. The aim of this study is to quantify enantiomer‐specific maturational pharmacokinetics (PK) of ketorolac and investigate if the contribution of both enantiomers to the total ketorolac concentration remains equal between infants and adults or if a change in target racemic concentration should be considered when applied to infants.

Methods

Data were pooled from 5 different studies in adults, children and infants, with 1020 plasma concentrations following single intravenous ketorolac administration. An allometry‐based enantiomer‐specific population PK model was developed with NONMEM 7.3. Simulations were performed in typical adults and infants to investigate differences in S‐ and R‐ketorolac exposure.

Results

S‐ and R‐ketorolac PK were best described with a 3‐ and a 2‐compartment model, respectively. The allometry‐based PK parameters accounted for changes between populations. No maturation function of ketorolac clearance could be identified. All model parameters were estimated with adequate precision (relative standard error <50%). Single dose simulations showed that a previously established analgesic concentration at half maximal effect in adults of 0.37 mg/L, had a mean S‐ketorolac concentration of 0.057 mg/L, but a mean S‐ketorolac concentration of 0.046 mg/L in infants. To match the effective adult S‐ketorolac‐concentration (0.057 mg/L) in typical infants, the EC_50‐racemic should be increased to 0.41 mg/L.

Conclusion

Enantiomer‐specific changes in ketorolac PK yield different concentrations and S‐ and R‐ketorolac ratios between infants and adults at identical racemic concentrations. These PK findings should be considered when studies on maturational pharmacodynamics are considered.

Delay differential equations for the description of Irbesartan pharmacokinetics: A population approach to model absorption complexities leading to dual peaks

Irbesartan is a poorly soluble BCS class II compound with weak acidic properties. After oral administration, dual peaks are noted in its concentration (C) - time (t) profile, a phenomenon that may be attributed to enterohepatic recirculation, gastric emptying and/or other absorption complexities related to its pH- and buffer capacity-dependent dissolution behavior. A population pharmacokinetic model, encompassing delay differential equations, was found the most appropriate approach to describe dual peaks in irbesartan's C-t profiles. Parameters estimated were: the absorption rate constant in the central compartment (ka = 0.304 h^-1), the constant time delay between the administration and the absorption (T=1.68 h), the apparent volume of distribution of the central (V1/F = 13.8 L) and peripheral (V2/F = 85.8 L) compartment, the apparent clearance from the central compartment (CL/F = 13.5 L/h), and the inter-compartmental clearance (Q/F = 17.7 L/h). Using simulations, it was made evident that changing the time delay results in significant changes of peak plasma concentrations but not of its blood pressure-lowering effect. In conclusion, delay differential equations may be useful to model dual peaks arising from absorption complexities, while changes of the time delay that reflect physiological processes that take place before absorption may have significant implications in proving bioequivalence.

Exposure-response analysis of Raltitrexed assessing liver toxicity

AIM

Raltitrexed (RTX) is a thymidylate synthase inhibitor with large pharmacokinetics (PK) variability that can be administered in case of 5-fluorouracil (5FU) intolerance or dihydropyrimidine dehydrogenase deficiency. While it is a more potent thymidylate synthase inhibitor than 5FU, RTX failed to replace this drug for colorectal cancer patients, mainly due to its toxicity at the recommended dose of 3 mg/m² every 3 weeks. However, every 2 weeks administration at 2 mg/m² demonstrated a favourable toxicity profile.

METHOD

We performed a randomized crossover comparative population PK study between every 2 weeks TOMOX (RTX 2 mg/m² ) and every 3 weeks TOMOX (RTX 3 mg/m² ).

RESULTS

A three-compartment model and a proportional error model best describe the data. Creatinine clearance and sex, but not body surface area (BSA), were covariates of RTX clearance leading to decrease of its interindividual variability of 28%. Weight and body surface area were covariates of central and peripheral volumes of distribution, respectively, leading to decreases of interindividual variability of 34.6% and 100%, respectively. In contrast to the dose, AUC was a good predictor of liver toxicity (P = 0.006, OR = 3.91, 95%CI = [1.48-10.34]). Using covariates to compute individual clearance and a threshold AUC (1.639, determined in this study), a covariates-based dose was calculated, leading to less variability in AUC than observed with the actual BSA-based or fixed doses.

CONCLUSION

These results advocate for the use of creatinine clearance and sex to determine the RTX dose instead of BSA.

Rapid Increase in Clearance of Phenobarbital in Neonates on Extracorporeal Membrane Oxygenation: A Pilot Retrospective Population Pharmacokinetic Analysis

OBJECTIVES

This study characterizes the changes in the pharmacokinetics of phenobarbital associated with extracorporeal membrane oxygenation treatment in neonates, to illustrate our findings and provide guidance on dosing.

DESIGN

Retrospective pilot population pharmacokinetic analysis.

SETTING

Neonatal ICU.

PATIENTS

Thirteen critically ill neonates (birth body weight, 3.21 kg [2.65-3.72 kg]; postnatal age at start of treatment: 2 d [0-7 d]; gestational age: 38 wk [38-41 wk]) receiving venovenous or venoarterial extracorporeal membrane oxygenation.

INTERVENTIONS

Phenobarbital administered in a loading dose of 7.5 mg/kg (8.5-16 mg/kg) and maintenance dose of 6.9 mg/kg/d (4.5-8.5 mg/kg/d).

MEASUREMENTS AND MAIN RESULTS

Therapeutic drug monitoring data were available, yielding 5, 31, and 19 phenobarbital concentrations before, during, and after extracorporeal membrane oxygenation, respectively. Population pharmacokinetic analysis was performed using NONMEM 7.3.0 (ICON Development Solutions, Ellicott City, MD). Maturation functions for clearance and volume of distribution were obtained from literature. In a one-compartment model, clearance and volume of distribution for a typical neonate off extracorporeal membrane oxygenation and with a median birth body weight (3.21 kg) at median postnatal age (2 d) were 0.0096 L/hr (relative SE = 11%)) and 2.72 L (16%), respectively. During extracorporeal membrane oxygenation, clearance was found to linearly increase with time. Upon decannulation, phenobarbital clearance initially decreased and subsequently increased slowly driven by maturation. Extracorporeal membrane oxygenation-related changes in volume of distribution could not be identified, possibly due to sparse data collection shortly after extracorporeal membrane oxygenation start. According to the model, target attainment is achieved in the first 12 days of extracorporeal membrane oxygenation with a regimen of a loading dose of 20 mg/kg and a maintenance dose of 4 mg/kg/d divided in two doses with an increase of 0.25 mg/kg every 12 hours during extracorporeal membrane oxygenation treatment.

CONCLUSIONS

We found a time-dependent increase in phenobarbital clearance during the first 12 days of extracorporeal membrane oxygenation treatment in neonates, which results in continuously decreasing phenobarbital exposure and increases the risk of therapeutic failure over time. Due to high unexplained variability, frequent and repeated therapeutic drug monitoring should be considered even with the model-derived regimen.

A Randomized, Controlled Pharmacokinetic and Pharmacodynamics Trial of Ambrisentan After Fontan Surgery

OBJECTIVES

To determine the pharmacokinetics, pharmacodynamics, and safety of the hepatically metabolized endothelin receptor antagonist, ambrisentan in children after Fontan surgery.

DESIGN

Prospective, randomized, double-blind, placebo-controlled pharmacokinetic/pharmacodynamics and safety trial.

SETTING

Single-center, postoperative cardiac ICU.

PATIENTS

Children undergoing elective Fontan surgery.

INTERVENTIONS

Subjects randomized on postoperative day number 1 to short-term (3 d) treatment with oral ambrisentan (2.5 mg in suspension, daily) versus placebo (4:1 randomization).

MEASUREMENTS AND MAIN RESULTS

Plasma drug concentrations were measured at 0.5, 1, 2, 4, and 18-36 hours after the first dose. We developed a population pharmacokinetic model in NONMEM 7.2 (Icon Solutions, Ellicott City, MD) and applied the model to dose-exposure simulations. Pharmacodynamics endpoints were assessed at baseline and 3 hours after study drug administration, using postoperative hemodynamic monitoring lines. The analysis included 16 patients, 13 on ambrisentan (77 plasma samples); median age 36 months (range, 26-72 mo), weight 13.3 kg (11.1-17.6 kg), and nine males. There were no differences in baseline characteristics between ambrisentan and controls. A one-compartment model with first-order absorption and lag-time characterized the data well. Allometrically scaled weight was the only covariate retained in the final model. Typical values for clearance and volume of distribution were lower than previously reported in adults, 1 L/hr/70 kg and 13.7 L/70 kg, respectively. Simulated exposures with doses of 0.1-0.2 mg/kg approximated therapeutic exposures in adults with pulmonary arterial hypertension receiving 5 mg or 10 mg doses. Ambrisentan lowered plasma brain natriuretic peptide concentrations (452 ± 479 to 413 ± 462; p = 0.046), Fontan pressures (16.8 ± 2.9 to 15.6 ± 2.9; p = 0.01), and indexed pulmonary vascular resistance (2.3 ± 0.9 to 1.8 ± 0.6; p = 0.01) with no drug-related adverse events.

CONCLUSIONS

Ambrisentan clearance is reduced following Fontan surgery, perhaps reflecting abnormal hepatic metabolism in this population. The observed safety profile appears favorable and hemodynamic effects of ambrisentan may be beneficial for Fontan patients.

Personalized Tuberculosis Treatment Through Model-Informed Dosing of Rifampicin

BACKGROUND AND OBJECTIVE

This study proposes a model-informed approach for therapeutic drug monitoring (TDM) of rifampicin to improve tuberculosis (TB) treatment.

METHODS

Two datasets from pulmonary TB patients were used: a pharmacokinetic study (34 patients, 373 samples), and TDM data (96 patients, 391 samples) collected at Radboud University Medical Center, The Netherlands. Nine suitable population pharmacokinetic models of rifampicin were identified in the literature and evaluated on the datasets. A model developed by Svensson et al. was found to be the most suitable based on graphical goodness of fit, residual diagnostics, and predictive performance. Prediction of individual area under the concentration-time curve from time zero to 24 h (AUC₂₄) and maximum concentration (C_max) employing various sampling strategies was compared with a previously established linear regression TDM strategy, using sampling at 2, 4, and 6 h, in terms of bias and precision (mean error [ME] and root mean square error [RMSE]).

RESULTS

A sampling strategy using 2- and 4-h blood collection was selected to be the most suitable. The bias and precision of the two strategies were comparable, except that the linear regression strategy was more biased in prediction of the AUC₂₄ than the model-informed approach (ME of 9.9% and 1.5%, respectively). A comparison of resulting dose advice, using predictions on a simulated dataset, showed no significant difference in sensitivity or specificity between the two methods. The model was successfully implemented in the InsightRX precision dosing platform.

CONCLUSION

Blood sampling at 2 and 4 h, combined with model-based prediction, can be used instead of the currently used linear regression strategy, shortening the sampling by 2 h and one sampling point without performance loss while simultaneously offering flexibility in sampling times.

Modeling the Effect of DAV132, a Novel Colon-Targeted Adsorbent, on Fecal Concentrations of Moxifloxacin and Gut Microbiota Diversity in Healthy Volunteers

To prevent antibiotic-induced perturbations on gut microbiota, DAV132, a novel colon-targeted adsorbent, which sequesters antibiotic residues in the lower gastrointestinal tract, was developed. We built an integrated pharmacological model of how DAV132 reduces fecal free moxifloxacin and preserves gut microbiota. We used plasma and fecal free moxifloxacin concentrations, and Shannon diversity index from 16S ribosomal RNA gene metagenomics analysis of fecal microbiota, of 143 healthy volunteers assigned randomly to receive moxifloxacin only, or with 10 DAV132 dose regimens, or to a control group. We modeled reduced fecal moxifloxacin concentrations using a transit model for DAV132 kinetics and a Michaelis-Menten model with an effect of the amount of activated charcoal on adsorption efficacy. Changes in moxifloxacin-induced perturbations on gut microbiota diversity were then quantified through a turnover model with the Emax model. With the developed model, the efficiency of pharmacokinetic antagonism and its consequences on gut microbiota diversity were quantified.

Population pharmacokinetics of ceftazidime in critically ill children: impact of cystic fibrosis

BACKGROUND

Pharmacokinetics data on ceftazidime are sparse for the paediatric population, particularly for children with cystic fibrosis (CF) or severe infections.

OBJECTIVES

To characterize the population pharmacokinetics of ceftazidime in critically ill children, identify covariates that affect drug disposition and evaluate the current dosing regimens.

METHODS

The study was registered with Clinicaltrials.gov (NCT01344512). Children receiving ceftazidime were selected in 13 French hospitals. Plasma concentrations were determined by UPLC-MS/MS. Population pharmacokinetic analyses were performed using NONMEN software.

RESULTS

One hundred and eight patients, aged 28 days to 12 years, with CF (n = 32), haematology and/or oncology disorders (n = 47) or severe infection (n = 29) were included. Ceftazidime was administered by continuous or intermittent infusions; 271 samples were available for analysis. A two-compartment model with first-order elimination and allometric scaling was developed and covariate analysis showed that ceftazidime pharmacokinetics were also significantly affected by CLCR and CF. Ceftazidime clearance was 82% higher in CF than in non-CF patients. Monte Carlo simulations showed that the percentage of target attainment (PTA) for the target of T>MIC = 65% was (i) lower in CF than in non-CF children with intermittent infusions and (ii) higher with continuous than intermittent infusion in all children.

CONCLUSIONS

The population pharmacokinetics model for ceftazidime in children was influenced by body weight, CLCR and CF. A higher PTA was obtained with continuous versus intermittent infusions. Further studies should explore the benefits of continuous versus intermittent infusion of ceftazidime, including current versus increased doses in CF children.

Aminoglycosides in critically ill patients: which dosing regimens for which pathogens?

Modifications of antibiotic pharmacokinetic parameters have been reported in critically ill patients, resulting in a risk of treatment failure. We aimed to determine optimised amikacin (AMK), gentamicin (GEN) and tobramycin (TOB) intravenous dosing regimens in this patient population. Patients admitted to the medical ICU and treated with AMK, GEN or TOB were included. Analyses were performed using a parametric population approach. Monte Carlo simulations were performed and the probability of target attainment (PTA) was calculated using C_max/MIC ≥ 8 and trough concentrations as targets. A total of 117 critically ill hospitalised patients were studied. Median values (interindividual variability, ɷ²) of clearance were 3.51 (0.539), 3.53 (0.297), 2.70 (0.339) and 5.07 (0.339) L/h for AMK, GEN, TOB, and TOB in cystic fibrosis (CF), respectively. Median values (ɷ²) of central volume of distribution were 30.2 (0.215), 20.0 (0.109) and 25.6 (0.177) L for AMK, GEN and TOB, respectively. Simulations showed that doses should be adjusted to actual body weight and creatinine clearance (CLCR) for AMK and GEN, and according to CLCR and presence of CF for TOB. In conclusion, our recommendations for treating Pseudomonas aeruginosa infections in this population include using initial doses of 35 mg/kg for AMK or 10 mg/kg for TOB (CF and non-CF patients). GEN demonstrated the best rates of target attainment against Staphylococcus aureus infections with a dose of 5 mg/kg. As high aminoglycoside doses are required in this population, efficacy and safety targets are conflicting and therapeutic drug monitoring remains an important tool to manage this issue.

Developmental population pharmacokinetics of caffeine in Chinese premature infants with apnoea of prematurity: A post-marketing study to support paediatric labelling in China

AIMS

The aim of the study was to evaluate the suitability of the current caffeine dosing regimen for the Chinese population using modelling and simulation approach.

METHODS

Pharmacokinetic samples were collected from 99 Chinese newborns with premature apnoea. The median (range) of gestational age and postmenstrual age were 28.3 (25.0-33.4) weeks and 31.1 (26.4-38.0) weeks, respectively. Newborns were receiving caffeine citrate at a loading dose of 20 mg/kg/d and a maintenance dose of 5-10 mg/kg/d. Caffeine concentrations and CYP1A2 polymorphisms were investigated. Population pharmacokinetic modelling of caffeine in Chinese preterm newborn on a population-wide scale was conducted using NONMEM.

RESULTS

A 1-compartment model with first-order elimination was used to describe population pharmacokinetic. With current weight implemented using 0.75 allometric scaling, clearance (CL) was positively related to current weight and postmenstrual age, but a negative relationship was observed with serum creatinine concentration. Eight genotypes of CYP1A2 were tested and none of them had a significant impact on caffeine pharmacokinetic parameters. Interindividual variability of CL and volume of distribution was 7.70 and 65.9%. The median (range) of 95% confidence intervals of CL were 0.0128 (0.0128-0.0131) L/h/kg. Monte Carlo simulation demonstrated that 80% (loading dose) and 98% (maintenance dose) of premature infants treated with a labelled dosing regimen attained the concentration target range of 5-20 mg/L.

CONCLUSION

A population PK model of caffeine was developed in Chinese newborns. Body weight-implemented allometric scaling, postmenstrual age and serum creatinine concentration markedly affected caffeine clearance. The labelled dosing regimen is suitable for Chinese premature infants.

Risperidone plasma concentrations are associated with side effects and effectiveness in children and adolescents with autism spectrum disorder

AIM

Risperidone is the most commonly prescribed antipsychotic drug to children and adolescents worldwide, but it is associated with serious side effects, including weight gain. This study assessed the relationship of risperidone and 9-hydroxyrisperidone trough concentrations, maximum concentrations and 24-hour area under the curves (AUCs) with body mass index (BMI) z-scores in children and adolescents with autism spectrum disorder (ASD) and behavioural problems. Secondary outcomes were metabolic, endocrine, extrapyramidal and cardiac side effects and effectiveness.

METHODS

Forty-two children and adolescents (32 males) aged 6-18 years were included in a 24-week prospective observational trial. Drug plasma concentrations, side effects and effectiveness were measured at several time points during follow-up. Relevant pharmacokinetic covariates, including medication adherence and CYP2D6, CYP3A4, CYP3A5 and P-glycoprotein (ABCB1) genotypes, were measured. Nonlinear mixed-effects modelling (NONMEM®) was used for a population pharmacokinetic analysis with 205 risperidone and 205 9-hydroxyrisperidone concentrations. Subsequently, model-based trough concentrations, maximum concentrations and 24-hour AUCs were analysed to predict outcomes using generalized and linear mixed-effects models.

RESULTS

A risperidone two-compartment model combined with a 9-hydroxyrisperidone one-compartment model best described the measured concentrations. Of all the pharmacokinetic parameters, higher risperidone sum trough concentrations best predicted higher BMI z-scores during follow-up (P < .001). Higher sum trough concentrations also predicted more sedation (P < .05), higher prolactin levels (P < .001) and more effectiveness measured with Aberrant Behavior Checklist irritability score (P < .01).

CONCLUSION

Our results indicate a therapeutic window exists, which suggests that therapeutic drug monitoring of risperidone might increase safety and effectiveness in children and adolescents with ASD and behavioural problems.

Exploration of Reduced Doses and Short-Cycle Therapy for Darunavir/Cobicistat in Patients with HIV Using Population Pharmacokinetic Modeling and Simulations

BACKGROUND AND OBJECTIVES

Protease inhibitors such as darunavir are an important therapeutic option in the anti-human immunodeficiency virus arsenal. Current dosage guidelines recommend using cobicistat- or ritonavir-boosted darunavir 800 mg every 24 h (q24h) in protease inhibitor-naïve patients, or ritonavir-boosted darunavir 600 mg q12h in experienced patients. However, darunavir displays a large, poorly characterized, inter-individual pharmacokinetic variability. The objectives of this study were to investigate the pharmacokinetics of darunavir and to elucidate the sources of its inter-individual variability using population pharmacokinetic modeling. Then, to determine the appropriateness of current treatment guidelines and the feasibility of alternative dosing regimens in a representative cohort of adult patients using simulations.

METHODS

Sparse pharmacokinetic samples were collected in 127 patients with human immunodeficiency virus type 1 infection, then supplemented with rich sampling data from a subset of 12 individuals. Data were analyzed using the nonlinear mixed-effects modeling software NONMEM. The effect of reduced doses (600 mg q24h and 400 mg q24h) or reduced frequency of administration (800 mg q24h for 5 days followed by 2 days of treatment interruption) was simulated.

RESULTS

Our model adequately described the pharmacokinetics of darunavir. Predictors of individual exposure were CYP3A5*3 and SLCO3A1 rs8027174 genotypes, sex, and alpha-1 acid glycoprotein level. No relationship was apparent between darunavir area under the curve and treatment efficacy or safety. For reduced dose regimens, darunavir concentrations remained above the protein binding-corrected EC₅₀ in the majority of subjects. More stringent pharmacokinetic targets were not reached in a significant proportion of patients.

CONCLUSIONS

These results add to the growing body of evidence that darunavir-based therapy could be simplified to reduce costs and toxicity, as well as to improve patient compliance. However, the heterogeneity in pharmacokinetic response should be considered when assessing whether individual patients could benefit from a particular regimen, for instance through the use of population pharmacokinetic models.

CLINICAL TRIAL REGISTRATION

ClinicalTrials.gov identifier: NCT03101644, date of registration: 5 April, 2017.

Reappraisal of the Optimal Dose of Meropenem in Critically Ill Infants and Children: a Developmental Pharmacokinetic-Pharmacodynamic Analysis

Data of developmental pharmacokinetics (PK) of meropenem in critically ill infants and children with severe infections are limited. We assessed the population PK and defined the appropriate regimen to optimize treatment in this population based on developmental PK-pharmacodynamic (PD) analysis. Blood samples were collected from pediatric intensive care unit patients with severe infection treated with standard dosage regimens for meropenem. Population PK data were analyzed using NONMEM software. Fifty-seven patients (mean age, 2.96 years [range, 0.101 to 14.4]; mean body weight, 15.8 kg [range, 5.0 to 65.0]) were included. A total of 135 meropenem concentrations were obtainable for population PK modeling. The median number of samples per patients was 2 (range, 1 to 4). A two-compartment model with first-order elimination was optimal for PK modeling. Weight and creatinine clearance (estimated by the Schwartz formula) were significantly correlated with the PK parameters of meropenem. The probabilities of target attainment for pathogens with low MICs of 1 and 2 μg/ml were 87.5% and 68.6% following administration of 40 mg/kg/dose (every 8 h [q8h]) as a 4-h infusion and 98.0% and 73.3% with high MICs of 4 and 8 μg/ml following administration of 110 mg/kg/day as a continuous infusion in critically ill infants and children under 70% fT _>MIC (the free time during which the plasma concentration of meropenem exceeds the MIC), respectively. The standard dosage regimens for meropenem did not meet an appropriate PD target, and an optimal dosing regimen was established in critically ill infants and children. (This study has been registered at ClinicalTrials.gov under identifier NCT03643497.).

Population Pharmacokinetics of Fludarabine in Children and Adults during Conditioning Prior to Allogeneic Hematopoietic Cell Transplantation

Background

Fludarabine is often used as an important drug in reduced toxicity conditioning regimens prior to hematopoietic cell transplantation (HCT). As no definitive pharmacokinetic (PK) basis for HCT dosing for the wide age and weight range in HCT is available, linear body surface area (BSA)-based dosing is still used.

Objective

We sought to describe the population PK of fludarabine in HCT recipients of all ages.

Methods

From 258 HCT recipients aged 0.3–74 years, 2605 samples were acquired on days 1 (42%), 2 (17%), 3 (4%) and 4 (37%) of conditioning. Herein, the circulating metabolite of fludarabine was quantified, and derived concentration-time data were used to build a population PK model using non-linear mixed-effects modelling.

Results

Variability was extensive where area under the curve ranged from 10 to 66 mg h/L. A three-compartment model with first-order kinetics best described the data. Actual body weight (BW) with standard allometric scaling was found to be the best body-size descriptor for all PK parameters. Estimated glomerular filtration rate (eGFR) was included as a descriptor of renal function. Thus, clearance was differentiated into a non-renal (3.24 ± 20% L/h/70 kg) and renal (eGFR × 0.782 ± 11% L/h/70 kg) component. The typical volumes of distribution of the central (V1), peripheral (V2), and second peripheral (V3) compartments were 39 ± 8%, 20 ± 11%, and 50 ± 9% L/70 kg respectively. Intercompartmental clearances between V1 and V2, and V1 and V3, were 8.6 ± 8% and 3.8 ± 13% L/h/70 kg, respectively.

Conclusion

BW and eGFR are important predictors of fludarabine PK. Therefore, current linear BSA-based dosing leads to highly variable exposure, which may lead to variable treatment outcome.

Population pharmacokinetics and dose optimization of ceftriaxone for children with community-acquired pneumonia

PURPOSE

To assess ceftriaxone population pharmacokinetics in a large pediatric population and describe the proper dose for establishing an optimized antibiotic regimen.

METHODS

From pediatric patients using ceftriaxone, blood samples were obtained and the concentration was measured using high-performance liquid chromatography ultraviolet detection. The NONMEM software program was used for population pharmacokinetic analysis, for which data from 99 pediatric patients (2 to 12 years old) was collected and 175 blood concentrations were obtained.

RESULTS

The best fit with the data was shown by the one-compartment model with first-order elimination. According to covariate analysis, weight had a significant impact on the clearance of ceftriaxone. Using Monte Carlo simulation, in a pediatric population with community-acquired pneumonia, a dose regimen of 100 mg/kg every 24 h produced satisfactory target attainment rates while remaining within the required minimum inhibitory concentration (2 mg/L).

CONCLUSION

Population pharmacokinetics of ceftriaxone was evaluated in children and an optimum dosing regimen was constructed on the basis of the pharmacokinetics-pharmacodynamics model-based approach.

Use of normalized prediction distribution errors for assessing population physiologically-based pharmacokinetic model adequacy

Currently employed methods for qualifying population physiologically-based pharmacokinetic (Pop-PBPK) model predictions of continuous outcomes (e.g., concentration-time data) fail to account for within-subject correlations and the presence of residual error. In this study, we propose a new method for evaluating Pop-PBPK model predictions that account for such features. The approach focuses on deriving Pop-PBPK-specific normalized prediction distribution errors (NPDE), a metric that is commonly used for population pharmacokinetic model validation. We describe specific methodological steps for computing NPDE for Pop-PBPK models and define three measures for evaluating model performance: mean of NPDE, goodness-of-fit plots, and the magnitude of residual error. Utility of the proposed evaluation approach was demonstrated using two simulation-based study designs (positive and negative control studies) as well as pharmacokinetic data from a real-world clinical trial. For the positive-control simulation study, where observations and model simulations were generated under the same Pop-PBPK model, the NPDE-based approach denoted a congruency between model predictions and observed data (mean of NPDE =  - 0.01). In contrast, for the negative-control simulation study, where model simulations and observed data were generated under different Pop-PBPK models, the NPDE-based method asserted that model simulations and observed data were incongruent (mean of NPDE =  - 0.29). When employed to evaluate a previously developed clindamycin PBPK model against prospectively collected plasma concentration data from 29 children, the NPDE-based method qualified the model predictions as successful (mean of NPDE = 0). However, when pediatric subpopulations (e.g., infants) were evaluated, the approach revealed potential biases that should be explored.

Population Pharmacokinetics of Vemurafenib in Children With Recurrent/Refractory BRAF Gene V600E-Mutant Astrocytomas

Vemurafenib (Zelboraf) is an orally available BRAF^V600E inhibitor approved for the treatment of unresectable or metastatic BRAF^V600E -mutant melanoma. The primary objective of this work was to characterize the pharmacokinetics (PK) of vemurafenib in pediatric patients with recurrent/refractory astrocytomas harboring the BRAF^V600E mutation. The study was also designed to evaluate the feasibility of replacing whole vemurafenib tablets with crushed tablets in young children unable to swallow tablets. Twenty-five pediatric patients (median age, 8.8 years; range, 3.3-19.2) with recurrent/refractory BRAF^V600E -mutant astrocytomas received whole (n = 19) or crushed (n = 6) vemurafenib tablets twice daily. Plasma samples were collected on days 1, 15, and 22 in cycle 1 of vemurafenib treatment. Descriptive PK analyses demonstrated significant variability (approximately 6-fold) in drug exposure. A 1-compartment model with first-order absorption and elimination was developed by adjusting the vemurafenib PK model previously validated in adults with mutant BRAF^V600E melanoma. After inclusion of allometric scaling on total body weight, the model adequately described the PK of vemurafenib in children between a wide age range of 3 to 19 years old. In the crushed-tablet cohort, relative bioavailability was approximately 96% (95% confidence interval, 49%-142%) compared to that seen in pediatric patients receiving whole tablets based on the preliminary comparison analysis results. Moderate intrapatient variability (48%) of vemurafenib clearance was observed. There was significant correlation (R² = 0.83) between area under the plasma concentration-time curve and trough concentration at steady state. These results will help increase the number of pediatric patients for whom vemurafenib is accessible and facilitate improved dosing in pediatric patients with recurrent/refractory BRAF^V600E astrocytomas.

First dose in neonates: pharmacokinetic bridging study from juvenile mice to neonates for drugs metabolized by CYP3A

First dose prediction is challenging in neonates. Our objective in this proof-of-concept study was to perform a pharmacokinetic (PK) bridging study from juvenile mice to neonates for drugs metabolized by CYP3A. We selected midazolam and clindamycin as model drugs. We developed juvenile mice population PK models using NONMEM. The PK parameters of these two drugs in juvenile mice were used to bridge PK parameters in neonates using different correction methods. The bridging results were evaluated by the fold-error of 0.5- to 1.5-fold. Simple allometry with and without a correction factor for maximum lifespan potential could be used for a bridging of clearance (CL) and volume of distribution (V_d), respectively, from juvenile mice to neonates. Simulation results demonstrated that for midazolam, 100% of clinical studies for which both the predictive CL and V_d were within 0.5- to 1.5-fold of the observed. For clindamycin, 75% and 100% of clinical studies for which the predictive CL and V_d were within 0.5- to 1.5-fold of the observed. A PK bridging of drugs metabolized by CYP3A is feasible from juvenile mice to neonates. It could be a complement to the ADE and PBPK models to support the first dose in neonates.

Enantiomer specific pharmacokinetics of ibuprofen in preterm neonates with patent ductus arteriosus

Aims

Racemic ibuprofen is widely used for the treatment of preterm neonates with patent ductus arteriosus. Currently used bodyweight‐based dosing guidelines are based on total ibuprofen, while only the S‐enantiomer of ibuprofen is pharmacologically active. We aimed to optimize ibuprofen dosing for preterm neonates of different ages based on an enantiomer‐specific population pharmacokinetic model.

Methods

We prospectively collected 210 plasma samples of 67 preterm neonates treated with ibuprofen for patent ductus arteriosus (median gestational age [GA] 26 [range 24–30] weeks, median body weight 0.83 [0.45–1.59] kg, median postnatal age [PNA] 3 [1–12] days), and developed a population pharmacokinetic model for S‐ and R‐ibuprofen.

Results

We found that S‐ibuprofen clearance (CL_S, 3.98 mL/h [relative standard error {RSE} 8%]) increases with PNA and GA, with exponents of 2.25 (RSE 6%) and 5.81 (RSE 15%), respectively. Additionally, a 3.11‐fold higher CL_S was estimated for preterm neonates born small for GA (RSE 34%). Clearance of R‐ibuprofen was found to be high compared to CL_S (18 mL/h [RSE 24%]), resulting in a low contribution of R‐ibuprofen to total ibuprofen exposure. Current body weight was identified as covariate on both volume of distribution of S‐ibuprofen and R‐ibuprofen.

Conclusion

S‐ibuprofen clearance shows important maturation, especially with PNA, resulting in an up to 3‐fold increase in CL_S during a 3‐day treatment regimen. This rapid increase in clearance needs to be incorporated in dosing guidelines by adjusting the dose for every day after birth to achieve equal ibuprofen exposure.

Population pharmacokinetics and target attainment of ciprofloxacin in critically ill patients

Purpose

To develop and validate a population pharmacokinetic model of ciprofloxacin intravenously in critically ill patients, and determine target attainment to provide guidance for more effective regimens.

Methods

Non-linear mixed-effects modelling was used for the model development and covariate analysis. Target attainment of an ƒAUC_0–24/MIC ≥ 100 for different MICs was calculated for standard dosing regimens. Monte Carlo simulations were performed to define the probability of target attainment (PTA) of several dosing regimens.

Results

A total of 204 blood samples were collected from 42 ICU patients treated with ciprofloxacin 400–1200 mg/day, with median values for age of 66 years, APACHE II score of 22, BMI of 26 kg/m², and eGFR of 58.5 mL/min/1.73 m². The median ƒAUC_0–24 and ƒC_max were 29.9 mg•h/L and 3.1 mg/L, respectively. Ciprofloxacin pharmacokinetics were best described by a two-compartment model. We did not find any significant covariate to add to the structural model. The proportion of patients achieving the target ƒAUC_0–24/MIC ≥ 100 were 61.9% and 16.7% with MICs of 0.25 and 0.5 mg/L, respectively. Results of the PTA simulations suggest that a dose of ≥ 1200 mg/day is needed to achieve sufficient ƒAUC_0–24/MIC ratios.

Conclusions

The model described the pharmacokinetics of ciprofloxacin in ICU patients adequately. No significant covariates were found and high inter-individual variability of ciprofloxacin pharmacokinetics in ICU patients was observed. The poor target attainment supports the use of higher doses such as 1200 mg/day in critically ill patients, while the variability of inter-individual pharmacokinetics parameters emphasizes the need for therapeutic drug monitoring to ensure optimal exposure.

Electronic supplementary material

The online version of this article (10.1007/s00228-020-02873-5) contains supplementary material, which is available to authorized users.

Population Pharmacokinetic Analysis of Tiropramide in Healthy Korean Subjects

The purpose of this study was to perform population pharmacokinetic (PPK) analysis of tiropramide in healthy Korean subjects, as well as to investigate the possible effects of various covariates on pharmacokinetic (PK) parameters of tiropramide. Although tiropramide is commonly used in digestive system-related diseases as an antispasmodic, PPK reporting and factors affecting PKs are not clearly reported. Thus, this study for healthy subjects is very significant because it could find new covariates in patients that had not been reported before or predict PPK for patients in the clinic by establishing PPK in healthy adults. By using Phoenix NLME, PK, demographic, and genetic data (collected to explain PK diversity of tiropramide in population) analyses were performed. As a basic model, a one-compartment with first-order absorption and lag-time was established and extended to include covariates that influenced the inter-subject variability. The total protein significantly influenced the distribution volume and systemic clearance of tiropramide, but genetic factors such as ABCB1 (1236C>T, 2677G>T/A, and 3435C>T), CYP2D6 (*1 and *10), OCT2 (808G>T), and PEPT1 (1287G>C) genes did not show any significant association with PK parameters of tiropramide. The final PPK model of tiropramide was validated, and suggested that some of the PK diversity in the population could be explained. Herein, we first describe the establishment of the PPK model of tiropramide for healthy Korean subjects, which may be useful as a dosing algorithm for the diseased population.

Population pharmacokinetics of riociguat and its metabolite in patients with chronic thromboembolic pulmonary hypertension from routine clinical practice

Pharmacokinetic data for riociguat in patients with chronic thromboembolic pulmonary hypertension (CTEPH) have previously been reported from randomized clinical trials, which may not fully reflect the population encountered in routine practice. The aim of the current study was to characterize the pharmacokinetic of riociguat and its metabolite M1 in the patients from routine clinical practice. A population pharmacokinetic model was developed in NONMEM 7.3, based on riociguat and its metabolite plasma concentrations from 49 patients with CTEPH. One sample with riociguat and M1 concentrations was available from each patient obtained at different time points after last dose. Age, bodyweight, sex, smoking status, concomitant medications, kidney and liver function markers were tested as potential covariates of pharmacokinetic of riociguat and its metabolite. Riociguat and M1 disposition was best described with one-compartment models. Apparent volume of distribution (Vd/F) for riociguat and M1 were assumed to be the same. Total bilirubin and creatinine clearance were the most predictive covariates for apparent riociguat metabolic clearance to M1 (CL_f,M1/F) and for apparent riociguat clearance through remaining pathways (CL_e,r/F), respectively. CL_f,M1/F, CL_e,r/F, Vd/F of riociguat and M1, and clearance of M1 (CL_e,M1/F) for a typical individual with 70 mL/min creatinine clearance and 0.69 mg/dL total bilirubin were 0.665 L/h (relative standard error = 17%)), 0.66 (18%) L/h, 3.63 (15%) L and 1.47 (19%) L/h, respectively. Upon visual identification of six outlying individuals, an absorption lag-time of 2.95 (6%) h was estimated for these patients. In conclusion, the only clinical characteristics related to riociguat exposure in patients with CTEPH from routine clinical practice are total bilirubin and creatinine clearance. This confirms the findings of the previous population pharmacokinetic studies based on data from randomized clinical trials.

Intrapulmonary concentrations of meropenem administered by continuous infusion in critically ill patients with nosocomial pneumonia: a randomized pharmacokinetic trial

Background

Optimal antimicrobial drug exposure in the lung is required for successful treatment outcomes for nosocomial pneumonia. Little is known about the intrapulmonary pharmacokinetics (PK) of meropenem when administered by continuous infusion (CI). The aim of this study was to evaluate the PK of two dosages of meropenem (3 g vs 6 g/day by CI) in the plasma and epithelial lining fluid (ELF) in critically ill patients with nosocomial pneumonia.

Methods

Thirty-one patients (81% male, median (IQR) age 72 (22) years) were enrolled in a prospective, randomized, clinical trial. Sixteen patients received 1 g/8 h and 15 2 g/8 h by CI (8 h infusion). Plasma and ELF meropenem concentrations were modeled using a population methodology, and Monte Carlo simulations were performed to estimate the probability of attaining (PTA) a free ELF concentration of 50% of time above MIC (50% fT>MIC), which results in logarithmic killing and the suppression of resistance in experimental models of pneumonia.

Results

The median (IQR) of meropenem AUC_0–24 h in the plasma and ELF was 287.6 (190.2) and 84.1 (78.8) mg h/L in the 1 g/8 h group vs 448.1 (231.8) and 163.0 (201.8) mg h/L in the 2 g/8 h group, respectively. The penetration ratio was approximately 30% and was comparable between the dosage groups. In the Monte Carlo simulations, only the highest approved dose of meropenem of 2 g/8 h by CI allowed to achieve an optimal PTA for all isolates with a MIC < 4 mg/L.

Conclusions

An increase in the dose of meropenem administered by CI achieved a higher exposure in the plasma and ELF. The use of the highest licensed dose of 6 g/day may be necessary to achieve an optimal coverage in ELF for all susceptible isolates (MIC ≤ 2 mg/L) in patients with conserved renal function. An alternative therapy should be considered when the presence of microorganisms with a MIC greater than 2 mg/L is suspected.

Trial registration

The trial was registered in the European Union Drug Regulating Authorities Clinical Trials Database (EudraCT-no. 2016-002796-10). Registered on 27 December 2016.

Population Pharmacokinetics of Intravenous Isavuconazole in Solid-Organ Transplant Recipients

Isavuconazole (ISA) is a triazole antifungal with activity against yeasts and molds. We established a population pharmacokinetic (pop PK) model of intravenous (i.v.) ISA to identify covariates that affect pharmacokinetics, using plasma samples from solid-organ transplant (SOT) recipients receiving peritransplant prophylaxis. Samples (n = 471) from 79 SOT recipients were utilized for pop PK analysis using nonlinear mixed-effect modeling NONMEM software. ISA (i.v.) PK parameters were best described by a two-compartment model. Significant covariates were sex on clearance (∼53% higher in women than men) and body mass index on peripheral volume of distribution. Incorporating drug exposure into Monte Carlo simulations, we demonstrated that standard ISA dosing is likely to attain the PK-pharmacodynamic (PD) target (area under the concentration curve/MIC ratio [AUC/MIC]) for treatment effectiveness against almost all infections caused by Aspergillus fumigatus isolates exhibiting MICs of ≤0.5 μg/ml (modal MIC). In contrast, standard dosing is predicted to attain PK-PD targets against <20% of infections caused by Candida albicans and Candida glabrata isolates exhibiting MICs of ≥0.016 and ≥0.5 μg/ml, respectively (modal MICs). These data are consistent with our SOT program's experience with ISA breakthrough infections, where 3 of 4 were caused by C. glabrata for which probabilities of PK-PD target attainment (PTA) were only 70% and 0% for isolates with MICs of 0.25 μg/ml and 1 μg/ml. Our findings provide important new insights into how ISA use might be optimized following SOT.

Population Pharmacokinetics of Alemtuzumab (Campath) in Pediatric Hematopoietic Cell Transplantation: Towards Individualized Dosing to Improve Outcome

BACKGROUND AND OBJECTIVE

Alemtuzumab (Campath^®) is used to prevent graft-versus-host disease and graft failure following pediatric allogeneic hematopoietic cell transplantation. The main toxicity includes delayed immune reconstitution, subsequent viral reactivations, and leukemia relapse. Exposure to alemtuzumab is highly variable upon empirical milligram/kilogram dosing.

METHODS

A population pharmacokinetic (PK) model for alemtuzumab was developed based on a total of 1146 concentration samples from 206 patients, aged 0.2-19 years, receiving a cumulative intravenous dose of 0.2-1.5 mg/kg, and treated between 2003 and 2015 in two centers.

RESULTS

Alemtuzumab PK were best described using a two-compartment model with a parallel saturable and linear elimination pathway. The linear clearance pathway, central volume of distribution, and intercompartmental distribution increased with body weight. Blood lymphocyte counts, a potential substrate for alemtuzumab, did not impact clearance.

CONCLUSION

The current practice with uniform milligram/kilogram doses leads to highly variable exposures in children due to the non-linear relationship between body weight and alemtuzumab PK. This model may be used for individualized dosing of alemtuzumab.

Population pharmacokinetics of sildenafil in extremely premature infants

AIMS

To characterize the population pharmacokinetics (PK) of sildenafil and its active metabolite, N-desmethyl sildenafil (DMS), in premature infants.

METHODS

We performed a multicentre, open-label trial to characterize the PK of sildenafil in infants ≤28 weeks gestation and < 365 postnatal days (cohort 1) or < 32 weeks gestation and 3-42 postnatal days (cohort 2). In cohort 1, we obtained PK samples from infants receiving sildenafil as ordered per the local standard of care (intravenous [IV] or enteral). In cohort 2, we administered a single IV dose of sildenafil and performed PK sampling. We performed a population PK analysis and dose-exposure simulations using the software NONMEM®.

RESULTS

We enrolled 34 infants (cohort 1 n = 25; cohort 2 n = 9) and collected 109 plasma PK samples. Sildenafil was given enterally (0.42-2.09 mg/kg) in 24 infants in cohort 1 and via IV (0.125 or 0.25 mg/kg) in all infants in cohort 2. A 2-compartment PK model for sildenafil and 1-compartment model for DMS, with presystemic conversion of sildenafil to DMS, characterized the data well. Coadministration of fluconazole (n = 4), a CYP3A inhibitor, resulted in an estimated 59% decrease in sildenafil clearance. IV doses of 0.125, 0.5 and 1 mg/kg every 8 hours (in the absence of fluconazole) resulted in steady-state maximum sildenafil concentrations that were generally within the range of those reported to inhibit phosphodiesterase type 5 activity in vitro.

CONCLUSIONS

We successfully characterized the PK of sildenafil and DMS in premature infants and applied the model to inform dosing for a follow-up, phase II study.

Pharmacokinetic modeling of intravenous sildenafil in newborns with congenital diaphragmatic hernia

PURPOSE

We developed a pharmacokinetic model of intravenous sildenafil in newborns with congenital diaphragmatic hernia (CDH) to achieve a target plasma concentration of over 50 μg/l.

METHODS

Twenty-three CDH newborns with pulmonary hypertension (64 blood samples) received intravenous sildenafil. Patients received a loading dose of 0.35 mg/kg (IQR 0.16 mg/kg) for 3 h, followed by a continuous infusion of 1.5 mg/kg/day (IQR 0.1 mg/kg/day). For model development, non-linear mixed modeling was used. Inter-individual variability (IIV) and inter-occasion variability were tested. Demographic and laboratory parameters were evaluated as covariates. Normalized prediction distribution errors (NPDE) and visual predictive check (VPC) were used for model validation.

RESULTS

A two-compartment disposition model of sildenafil and a one-compartment disposition model of desmethyl sildenafil (DMS) was observed with IIV in sildenafil and DMS clearance and volume of distribution of sildenafil. NPDE and VPC revealed adequate predictability. Only postnatal age increased sildenafil clearance. This was partly compensated by a higher DMS concentration, which also has a therapeutic effect. In this small group of patients, sildenafil was tolerated well.

CONCLUSIONS

This model for sildenafil in CDH patients shows that concentration-targeted sildenafil dosing of 0.4 mg/kg in 3 h, followed by 1.6 mg/kg/day continuous infusion achieves appropriate sildenafil plasma levels.

Rapidly maturing fentanyl clearance in preterm neonates

BACKGROUND

Fentanyl is frequently used off-label in preterm newborns. Due to very limited pharmacokinetic and pharmacodynamic data, fentanyl dosing is mostly based on bodyweight. This study describes the maturation of the pharmacokinetics in preterm neonates born before 32 weeks of gestation.

METHODS

442 plasma samples from 98 preterm neonates (median gestational age: 26.9 (range 23.9-31.9) weeks, postnatal age: 3 (range 0-68) days, bodyweight 1.00 (range 0.39-2.37) kg) were collected in an opportunistic trial and fentanyl plasma levels were determined. NONMEM V.7.3 was used to develop a population pharmacokinetic model and to perform simulations.

RESULTS

Fentanyl pharmacokinetics was best described by a two-compartment model. A pronounced non-linear influence of postnatal and gestational age on clearance was identified. Clearance (L/hour/kg) increased threefold, 1.3-fold and 1.01-fold in the first, second and third weeks of life, respectively. In addition, clearance (L/hour/kg) was 1.4-fold and 1.7-fold higher in case of a gestational age of 28 and 31 weeks, respectively, compared with 25 weeks. Volume of distribution changed linearly with bodyweight and was 8.7 L/kg. To achieve similar exposure across the entire population, a continuous infusion (µg/kg/hour) dose should be reduced by 50% and 25% in preterm neonates with a postnatal age of 0-4 days and 5-9 days in comparison to 10 days and older.

CONCLUSION

Because of low clearance, bodyweight-based dosages may result in fentanyl accumulation in neonates with the lowest postnatal and gestational ages which may require dose reduction. Together with additional information on the pharmacodynamics, the results of this study can be used to guide dosing.

Temocillin dosing in haemodialysis patients based on population pharmacokinetics of total and unbound concentrations and Monte Carlo simulations

Objectives

To develop a population model describing temocillin pharmacokinetics (PK) in patients undergoing haemodialysis and investigate how pharmacokinetic/pharmacodynamic (PD) targets can be met with different dosage regimens.

Patients and methods

Sixteen patients received the currently licenced dosing of 1, 2 or 3 g of temocillin (total of 61 doses) corresponding to an inter-dialytic period of 20, 44 or 68 h, respectively, and a dialysis period of 4 h. A non-linear mixed-effects model was developed jointly for total and unbound temocillin serum concentrations. The performance of clinically feasible dosing regimens was evaluated using a 5000-subject Monte Carlo (MC) simulation for determining the highest MIC for which the PK/PD target of 40%ƒT>MIC would be reached in 90% of patients [probability of target attainment (PTA)]. This PK study was registered at ClinicalTrials.gov (NCT02285075).

Results

Temocillin unbound and total serum concentrations (429 samples) were used to fit an open two-compartment model with non-linear albumin binding and first-order elimination. In addition to total body clearance, dialysis clearance was modelled using the Michaels function. The currently licenced dosing achieved a 90% PTA for an MIC up to 8 mg/L. A new temocillin dosage regimen was designed that would achieve a 90% PTA for an MIC of 16 mg/L (MIC90 of target organisms) adjusted to patient weight and inter-dialytic period.

Conclusions

Currently licensed dosage regimen is suboptimal for MICs >8 mg/L (frequently found in clinical isolates). Model-based simulations allowed suggestion of a new dosage regimen with improved probability of microbiological success, applicability in routine clinical practice and more appropriate for empirical therapy.

Nonlinear Mixed-Effects Model Development and Simulation Using nlmixr and Related R Open-Source Packages

nlmixr is a free and open-source R package for fitting nonlinear pharmacokinetic (PK), pharmacodynamic (PD), joint PK-PD, and quantitative systems pharmacology mixed-effects models. Currently, nlmixr is capable of fitting both traditional compartmental PK models as well as more complex models implemented using ordinary differential equations. We believe that, over time, it will become a capable, credible alternative to commercial software tools, such as NONMEM, Monolix, and Phoenix NLME.

Development and comparison of population pharmacokinetic models of vancomycin in neurosurgical patients based on two different renal function markers

WHAT IS KNOWN AND OBJECTIVES

Some previous studies have indicated that serum cystatin C (Cys C) is a better marker than serum creatinine (SCR) for assessing the glomerular filtering rate (GFR). However, in almost all population pharmacokinetic models of vancomycin, the GFR is usually estimated from SCR. Therefore, the aim of this study was to compare the GFR estimated from SCR (sGFR) with the GFR estimated from Cys C (cGFR) and investigate which one can describe the characteristics of vancomycin population pharmacokinetics better in Chinese neurosurgical adult patients.

METHODS

Patients from the Neurosurgery Department aged ≥18 years were enrolled retrospectively. Among these patients, the data from 222 patients were used to establish two population pharmacokinetic models based on sGFR and cGFR, separately. The data from another 95 patients were used for the external validation of these two models. Non-linear mixed-effect modelling (NONMEM) 7.4.3 was used for the population pharmacokinetic analysis.

RESULTS

We developed two one-compartment models with first-order absorption based on Cys C and SCR, separately. In the Cys C model, age, body weight and cGFR were significant covariates on the clearance rate (CL) of vancomycin (typical value, 6.4 L/hour). In the SCR model, age and sGFR were significant covariates on the CL (typical value, 6.46 L/hour). The external validation results showed that the predictive performance of the two models was similar.

WHAT IS NEW AND CONCLUSION

In this study, the predictive performance of two models was similar in neurosurgical patients. We did not find a significant improvement in the predictive performance of the model when GFR was estimated from Cys C.

Impact of antiretroviral treatment on height evolution of HIV infected children

BACKGROUND

Antiretroviral treatment (ART) has been shown to have a beneficial effect on the weight evolution but its effect on height remains unclear. We described patterns of height evolution and identified predictors of catch-up growth in HIV-infected children on ART.

METHODS

To describe the height evolution from birth to adulthood, we developed a nonlinear mixed effect model using data from perinatally HIV-infected children who initiated ART from 1999 to 2013 in a prospective cohort study in Thailand. The main covariates of interest were: sex, ART regimen (dual nucleoside reverse-transcriptase inhibitor, non-nucleoside reverse transcriptase inhibitor (NNRTI)-, or protease inhibitor (PI)-based), baseline CD4 percentage, HIV-RNA load and CDC HIV Classification stage and occurrence of AIDS-defining events.

RESULTS

A total 477 children (43% boys) contributed 18,596 height measurements over a median duration of 6.3 years on ART (interquartile range, 3.0 to 8.3). At ART initiation, median age was 6.2 years (1.8 to 9.6), 16% of children were underweight (weight-for-age z-score < - 2), 49% presented stunting (height-for-age z-score < - 2), and 7% wasting (weight-for-height z-score < - 2). The most frequent regimen at ART initiation was NNRTI-based (79%). A model with 4 components, birth length and 3 exponential functions of age accounting for the 3 growth phases was developed and show that the height-growth velocity was inversely associated with the age at ART initiation, the adult height was significantly lower in those who had experienced at least one AIDS-defining event while, as expected, the model found that adult height in females was lower than in males. Age at ART initiation, type of ART regimen, CDC stage, CD4 percentages, and HIV-RNA load were not associated with the final height.

CONCLUSIONS

The younger the children at ART initiation, the greater the effect on height-growth velocity, supporting the World Health Organization's recommendation to start ART as early as possible. However, final adult height was not linked to the age at ART initiation.

NONMEM Tutorial Part II: Estimation Methods and Advanced Examples

In this second tutorial on NONMEM, the examples of typical pharmacokinetic/pharmacodynamic modeling problems that occur in the pharmaceutical field will be presented, which the reader can use as a template for his or her own modeling endeavors. Each of the problems presented is challenging in some way, and the logic behind setting up each problem is discussed. Logical concepts of the problem itself as well as the technical aspect of how to set it up in NONMEM are described and demonstrated. The concepts behind the various estimation algorithms will first be described to allow the user a better understanding of how to use them.

A Population Pharmacokinetic Model to Predict the Individual Starting Dose of Tacrolimus Following Pediatric Renal Transplantation

Background

Multiple clinical, demographic, and genetic factors affect the pharmacokinetics of tacrolimus in children, yet in daily practice, a uniform body-weight based starting dose is used. It can take weeks to reach the target tacrolimus pre-dose concentration.

Objectives

The objectives of this study were to determine the pharmacokinetics of tacrolimus immediately after kidney transplantation and to find relevant parameters for dose individualization using a population pharmacokinetic analysis.

Methods

A total of 722 blood samples were collected from 46 children treated with tacrolimus over the first 6 weeks after renal transplantation. Non-linear mixed-effects modeling (NONMEM^®) was used to develop a population pharmacokinetic model and perform a covariate analysis. Simulations were performed to determine the optimal starting dose and to develop dosing guidelines.

Results

The data were accurately described by a two-compartment model with allometric scaling for bodyweight. Mean tacrolimus apparent clearance was 50.5 L/h, with an inter-patient variability of 25%. Higher bodyweight, lower estimated glomerular filtration rate, and higher hematocrit levels resulted in lower total tacrolimus clearance. Cytochrome P450 3A5 expressers and recipients who received a kidney from a deceased donor had a significantly higher tacrolimus clearance. The model was successfully externally validated. In total, these covariates explained 41% of the variability in clearance. From the significant covariates, the cytochrome P450 3A5 genotype, bodyweight, and donor type were useful to adjust the starting dose to reach the target pre-dose concentration. Dosing guidelines range from 0.27 to 1.33 mg/kg/day.

Conclusion

During the first 6 weeks after transplantation, the tacrolimus weight-normalized starting dose should be higher in pediatric kidney transplant recipients with a lower bodyweight, those who express the cytochrome P450 3A5 genotype, and those who receive a kidney from a deceased donor.

A Pediatric Covariate Function for CYP3A-Mediated Midazolam Clearance Can Scale Clearance of Selected CYP3A Substrates in Children

Recently a framework was presented to assess whether pediatric covariate models for clearance can be extrapolated between drugs sharing elimination pathways, based on extraction ratio, protein binding, and other drug properties. Here we evaluate when a pediatric covariate function for midazolam clearance can be used to scale clearance of other CYP3A substrates. A population PK model including a covariate function for clearance was developed for midazolam in children aged 1–17 years. Commonly used CYP3A substrates were selected and using the framework, it was assessed whether the midazolam covariate function accurately scales their clearance. For eight substrates, reported pediatric clearance values were compared numerically and graphically with clearance values scaled using the midazolam covariate function. For sildenafil, clearance values obtained with population PK modeling based on pediatric concentration-time data were compared with those scaled with the midazolam covariate function. According to the framework, a midazolam covariate function will lead to systemically accurate clearance scaling (absolute prediction error (PE) < 30%) for CYP3A substrates binding to albumin with an extraction ratio between 0.35 and 0.65 when binding < 10% in adults, between 0.05 and 0.55 when binding > 90%, and with an extraction ratio ranging between these values when binding between 10 and 90%. Scaled clearance values for eight commonly used CYP3A substrates were reasonably accurate (PE < 50%). Scaling of sildenafil clearance was accurate (PE < 30%). We defined for which CYP3A substrates a pediatric covariate function for midazolam clearance can accurately scale plasma clearance in children. This scaling approach may be useful for CYP3A substrates with scarce or no available pediatric PK information.

Body Surface Area-Based Dosing Regimen of Caspofungin in Children: a Population Pharmacokinetics Confirmatory Study

We evaluated the population pharmacokinetics of caspofungin in children (2 to 12 years of age). The real-world data from 48 children were best fit by a two-compartment model with first-order elimination. Subsequent covariate analysis demonstrated that body surface area had a significant correlation with caspofungin pharmacokinetics, compared to body weight. The population pharmacokinetics of caspofungin confirmed that adjustment of caspofungin dosage based on body surface area is most appropriate for pediatric use.

Population Pharmacokinetics and Dosing Optimization of Imipenem in Children with Hematological Malignancies

Imipenem is widely used for the treatment of children with serious infections. Currently, studies on the pharmacokinetics of imipenem in children with hematological malignancies are lacking. Given the significant impact of disease on pharmacokinetics and increased resistance, we aimed to conduct a population pharmacokinetic study of imipenem and optimize the dosage regimens for this vulnerable population. After children were treated with imipenem-cilastatin (IMP-CS), blood samples were collected from the children and the concentrations of imipenem were quantified using high-performance liquid chromatography with UV detection. Then, a population-level pharmacokinetic analysis was conducted using NONMEM software. Data were collected from 56 children (age range, 2.03 to 11.82 years) with hematological malignancies to conduct a population pharmacokinetic analysis. In this study, a two-compartment model that followed first-order elimination was found to be the most suitable. The parameters of current weight, age, and creatinine elimination rate were significant covariates that influenced imipenem pharmacokinetics. As a result, 41.4%, 56.1%, and 67.1% of the children reached the pharmacodynamic target (the percentage of the time during the total dosing interval that the free drug concentration remains above the MIC of 70%) against sensitive pathogens with an MIC of 0.5 mg/liter with imipenem at 15, 20, and 25 mg/kg of body weight every 6 h (q6h), respectively. However, only 11.1% of the children achieved the pharmacodynamic target against Pseudomonas aeruginosa isolates with an MIC of 2 mg/liter at a dose of 25 mg/kg q6h. The population pharmacokinetics of imipenem were assessed in children. The current dosage regimens of imipenem result in underdosing against resistant pathogens, including Pseudomonas aeruginosa and Acinetobacter baumannii However, for sensitive pathogens, imipenem has an acceptable pharmacodynamic target rate at a dosage of 25 mg/kg q6h. (The study discussed in this paper has been registered at ClinicalTrials.gov under identifier NCT03113344.).