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

Fixed parameters in the population pharmacokinetic modeling of valproic acid might not be suitable: external validation in Chinese adults with epilepsy or after neurosurgery

PURPOSE

This study aims to assess the predictive performance of published valproic acid (VPA) population pharmacokinetic (PPK) models using an external data set in Chinese adults with epilepsy or after neurosurgery.

METHODS

A total of 384 concentrations from 290 Chinese adults with epilepsy or after neurosurgery were used for external validation. Data on published VPA PPK models were extracted from the literature. Prediction-based diagnostics (such as F20 and F30), simulation-based diagnostics, and Bayesian forecasting were used to evaluate the predictability of models.

RESULTS

The results of prediction-based diagnostics of all models were unsatisfactory. Models B, F, and H showed the best prediction performance in simulation-based diagnostics and Bayesian forecasting, demonstrating superior precision and accuracy. Bayesian forecasting demonstrated significant improvements in the model predictability.

CONCLUSION

The published PPK models showed extensive variation in predictive performance for extrapolation among Chinese adults with epilepsy or after neurosurgery patients. Fixed parameters of Vd and Ka in the PPK modeling of VPA might be the reason for the unsatisfied predictive performance. Bayesian forecasting significantly improved model predictability and may help to individualize VPA dosing.

Optimal dosing regimen of caspofungin in adolescents with allogeneic haematopoietic stem cell transplantation

OBJECTIVES

The optimal dosing regimen of caspofungin in adolescents undergoing allogeneic haematopoietic stem cell transplantation against Candida spp. is unknown. The study aimed to compare body surface area (BSA)-based and fixed dosing regimens through population pharmacokinetic (PPK) analysis and to optimize dosing regimens likely to achieve therapeutic exposures.

METHODS

Opportunistic sampling was used to collect plasma concentrations through a prospective observational pharmacokinetic study. PPK analysis and Monte Carlo simulations (n = 1000) were performed using NONMEM.

RESULTS

A total of 86 samples of 30 adolescents (12-17 years old) were best described by a two-compartment pharmacokinetic model. BSA is the only covariate on clearance and central volume of distribution. For Candida glabrata and Candida albicans, a standard dosing regimen could achieve at least a 90% probability of target attainment for the indicator of AUC0-24/MIC90. Dosing regimen simulations identified a BSA cut-off value of 1.3 m2, where a fixed loading dose (LD) is preferred when BSA ≥ 1.3 m2 and a BSA-based LD is preferred when BSA < 1.3 m2. For maintenance dose (MD), however, the BSA-based dose was proposed, regardless of BSA. The current maximum dosing regimen of LD 70 mg/day and MD 70 mg/day could not result in sufficient antifungal exposure for Candida parapsilosis with MIC90 of 1 mg/L. Furthermore, an LD of 70 mg/day and MD of 60 mg/m2/day rendered 90.4% steady-state trough concentration (Ctrough) over 1 mg/L in the virtual population.

CONCLUSIONS

Our study proposed optimized dosing regimens of caspofungin based on AUC0-24/MIC90 or Ctrough, which may support further individualized treatment.

Using Real-World Data to Externally Evaluate Population Pharmacokinetic Models of Dexmedetomidine in Children and Infants

Dexmedetomidine is a sedative used in both adults and off-label in children with considerable reported pharmacokinetic (PK) interindividual variability affecting drug exposure across populations. Several published models describe the population PKs of dexmedetomidine in neonates, infants, children, and adolescents, though very few have been externally evaluated. A prospective PK dataset of dexmedetomidine plasma concentrations in children and young adults aged 0.01-19.9 years was collected as part of a multicenter opportunistic PK study. A PubMed search of studies reporting dexmedetomidine PK identified five population PK models developed with data from demographically similar children that were selected for external validation. A total of 168 plasma concentrations from 102 children were compared with both population (PRED) and individualized (IPRED) predicted values from each of the five published models by quantitative and visual analyses using NONMEM (v7.3) and R (v4.1.3). Mean percent prediction errors from observed values ranged from -1% to 120% for PRED, and -24% to 60% for IPRED. The model by James et al, which was developed using similar "real-world" data, nearly met the generalizability criteria from IPRED predictions. Other models developed using clinical trial data may have been limited by inclusion/exclusion criteria and a less racially diverse population than this study's opportunistic dataset. The James model may represent a useful, but limited tool for model-informed dosing of hospitalized children.

Population pharmacokinetics of total and unbound valemetostat and platelet dynamics in healthy volunteers and patients with non-Hodgkin lymphoma

Valemetostat is an EZH2/1 inhibitor that has been approved in Japan for the treatment of patients with relapsed/refractory adult T-cell leukemia/lymphoma, based mainly on results from a single-arm phase II trial. It is currently under investigation worldwide for the treatment of other non-Hodgkin lymphomas (NHLs), including peripheral T-cell lymphoma, and for solid tumors. Semi-mechanistic population pharmacokinetic modeling of total and unbound valemetostat and an analysis of the platelet time course during treatment with valemetostat were conducted using data from five clinical trials (two in patients with NHL and three in healthy volunteers). Pharmacokinetic data, including 3162 total/1871 unbound valemetostat observations from 102 patients and 72 healthy volunteers, were described by a three-compartment model with sequential zero-/first-order absorption and saturable binding in the central compartment. Alpha-1-acid glycoprotein (AAG) was the most influential covariate for total valemetostat exposure, yet had little impact on unbound exposure, meaning no dose adjustment was warranted based on AAG levels. The longitudinal platelet data from 101 patients (2313 observations) were adequately described by a modified Friberg model with two proliferation compartments, which characterized unique spontaneous recovery of platelet counts without dose modifications. A model-based simulation quantitatively assessed the proposed dose-adjustment guidance in case of platelet count decreased by comparing the probability of treatment discontinuation due to platelet count decreased with or without the dose adjustment. In summary, the models described observed total and unbound valemetostat concentrations and a unique time course of platelets during treatment, which can justify the clinical dose and provide dose-adjustment guidance.

Population pharmacokinetics of olanzapine in pediatric patients with psychiatric disorders

OBJECTIVE

To develop and validate a population pharmacokinetic (PPK) model of oral olanzapine in pediatric Chinese patients in order to individualize therapy in this population.

METHODS

A total of 897 serum concentrations from 269 pediatric patients taking oral olanzapine (ages 8-17 years) were collected. Demographic parameters, biological characteristics and concomitant medications were investigated as covariates. The data were analyzed using a nonlinear mixed-effects modeling approach. Bootstrapping (1000 runs), normalized prediction distribution error (NPDE), and external validation of 62 patients were employed. Simulations were performed to explore the individualized dosing regimens in various situations.

RESULTS

The one-compartment model with first-order absorption and elimination had an apparent clearance (CL/F) of 10.38 L/h, a distribution volume (V/F) of 9.41 L/kg and an absorption rate constant (Ka) fixed at 0.3 h-1. The equation was CL∕F (L∕h) = 10.38 × (body weight∕60)0.25 ×1.33 (if male) × 0.71 (if co-occurrence of infection) × 0.51 (if co-therapy with fluvoxamine) × 1.27 (if co-therapy with sertraline) × 1.43 (if co-therapy with valproate). The final model had satisfactory stability, robustness, and predictive ability. The results from a simulation suggested the oral olanzapine doses required for male and female pediatric patients weighing between 40 and 60 kg without co-medication were 10-15 mg/day and 7.5-10 mg/day, respectively, and dosage adjustments should be based on sex and body weight; and co-administrated with valproate, sertraline, or fluvoxamine.

CONCLUSION

This model may help individualize optimum dosing of oral olanzapine for pediatric patients.

Effect of Interindividual Variability in Metabolic Clearance and Relative Bioavailability on Rifampicin Exposure in Tuberculosis Patients with and without HIV Co-Infection: Does Formulation Quality Matter?

The present study aims to characterise the pharmacokinetics of rifampicin (RIF) in tuberculosis (TB) patients with and without HIV co-infection, considering the formation of 25-O-desacetyl-rifampicin (desRIF). It is hypothesised that the metabolite formation, HIV co-infection and drug formulation may further explain the interindividual variation in the exposure to RIF. Pharmacokinetic, clinical, and demographic data from TB patients with (TB-HIV+ group; n = 18) or without HIV (TB-HIV- group; n = 15) who were receiving RIF as part of a four-drug fixed-dose combination (FDC) regimen (RIF, isoniazid, pyrazinamide, and ethambutol) were analysed, along with the published literature data on the relative bioavailability of different formulations. A population pharmacokinetic model, including the formation of desRIF, was developed and compared to a model based solely on the parent drug. HIV co-infection does not alter the plasma exposure to RIF and the desRIF formation does not contribute to the observed variability in the RIF disposition. The relative bioavailability and RIF plasma exposure were significantly lower than previously reported for the standard regimen with FDC tablets. Furthermore, participants weighting less than 50 kg do not reach the same RIF plasma exposure as compared to those weighting >50 kg. In conclusion, as no covariate was identified other than body weight on CL/F and Vd/F, low systemic exposure to RIF is likely to be caused by the low bioavailability of the formulation.

Population pharmacokinetic modeling of sufentanil in adult Korean patients undergoing cardiopulmonary bypass surgery

Sufentanil is frequently used as an anesthetic agent in cardiac surgery owing to its cardiovascular safety and favorable pharmacokinetics. However, the pharmacokinetics profiles of sufentanil in patients undergoing cardiopulmonary bypass (CPB) surgery remain less understood, which is crucial for achieving the desired level of anesthesia and mitigating surgical complications. Therefore, this study aimed to develop a population pharmacokinetic model of sufentanil in patients undergoing CPB surgery and elucidate the clinical factors affecting its pharmacokinetic profile. Adult patients who underwent cardiac surgery with CPB and were administered sufentanil for anesthesia were enrolled. Arterial blood samples were collected to quantify plasma concentrations of sufentanil and clinical laboratory parameters, including inflammatory cytokines. A population pharmacokinetic model was established using nonlinear mixed-effects modeling. Simulations were performed using the pharmacokinetic parameters of the final model. Overall, 20 patients were included in the final analysis. Sufentanil pharmacokinetics were modeled using a two-compartment model, accounting for CPB effects. Sufentanil clearance increased 2.80-fold during CPB and warming phases, while the central compartment volume increased 2.74-fold during CPB. CPB was a significant covariate affecting drug clearance and distribution volume. No other significant covariates were identified despite increased levels of the inflammatory cytokines, including IL-6, IL-8, and TNF-α during CPB. The simulation indicated a 30 μg loading dose and 40 μg/h maintenance infusion for target-controlled infusion. Additionally, a bolus dose of 60 μg was added at CPB initiation to adjust for exposure changes during this phase. Considering the target sufentanil concentrations, a uniform dosing regimen was acceptable for effective analgesia.

Population pharmacokinetics of daptomycin in critically ill patients receiving extracorporeal membrane oxygenation

BACKGROUND

Daptomycin is widely used in critically ill patients for Gram-positive bacterial infections. Extracorporeal membrane oxygenation (ECMO) is increasingly used in this population and can potentially alter the pharmacokinetic (PK) behaviour of antibiotics. However, the effect of ECMO has not been evaluated in daptomycin. Our study aims to explore the effect of ECMO on daptomycin in critically ill patients through population pharmacokinetic (PopPK) analysis and to determine optimal dosage regimens based on both efficacy and safety considerations.

METHODS

A prospective, open-label PK study was carried out in critically ill patients with or without ECMO. The total concentration of daptomycin was determined by UPLC-MS/MS. NONMEM was used for PopPK analysis and Monte Carlo simulations.

RESULTS

Two hundred and ninety-three plasma samples were collected from 36 critically ill patients, 24 of whom received ECMO support. A two-compartment model with first-order elimination can best describe the PK of daptomycin. Creatinine clearance (CLCR) significantly affects the clearance of daptomycin while ECMO has no significant effect on the PK parameters. Monte Carlo simulations showed that, when the MICs for bacteria are  ≥1 mg/L, the currently recommended dosage regimen is insufficient for critically ill patients with CLCR > 30 mL/min. Our simulations suggest 10 mg/kg for patients with CLCR between 30 and 90 mL/min, and 12 mg/kg for patients with CLCR higher than 90 mL/min.

CONCLUSIONS

This is the first PopPK model of daptomycin in ECMO patients. Optimal dosage regimens considering efficacy, safety, and pathogens were provided for critical patients based on pharmacokinetic-pharmacodynamic analysis.

Does Sample Size, Sampling Strategy, or Handling of Concentrations Below the Lower Limit of Quantification Matter When Externally Evaluating Population Pharmacokinetic Models?

BACKGROUND AND OBJECTIVES

Precision dosing requires selecting the appropriate population pharmacokinetic model, which can be assessed through external evaluations (EEs). The lack of understanding of how different study design factors influence EE study outcomes makes it challenging to select the most suitable model for clinical use. This study aimed to evaluate the impact of sample size, sampling strategy, and handling of concentrations below the lower limit of quantification (BLQ) on the outcomes of EE for four population pharmacokinetic models using vancomycin and tobramycin as examples.

METHODS

Three virtual patient populations undergoing vancomycin or tobramycin therapy were simulated with varying sample size and sampling scenarios. The three approaches used to handle BLQ data were to (1) discard them, (2) impute them as LLOQ/2, or (3) use a likelihood-based approach. EEs were performed with NONMEM and R.

RESULTS

Sample size did not have an important impact on the EE results for a given scenario. Increasing the number of samples per patient did not improve predictive performance for two out of the three evaluated models. Evaluating a model developed with rich sampling did not result in better performance than those developed with regular therapeutic drug monitoring. A likelihood-based method to handle BLQ samples impacted the outcomes of the EE with lower bias for predicted troughs.

CONCLUSIONS

This study suggests that a large sample size may not be necessary for an EE study, and models selected based on TDM may be more generalizable. The study highlights the need for guidelines for EE of population pharmacokinetic models for clinical use.

Optimization of Vancomycin Initial Dosing Regimen in Neonates Using an Externally Evaluated Population Pharmacokinetic Model

BACKGROUND

Vancomycin therapeutic monitoring guidelines were revised in March 2020, and a population pharmacokinetics-guided Bayesian approach to estimate the 24-hour area under the concentration-time curve to the minimum inhibitory concentration ratio has since been recommended instead of trough concentrations. To comply with these latest guidelines, we evaluated published population pharmacokinetic models of vancomycin using an external dataset of neonatal patients and selected the most predictive model to develop a new initial dosing regimen.

METHODS

The models were identified from the literature and tested using a retrospective dataset of Canadian neonates. Their predictive performance was assessed using prediction- and simulation-based diagnostics. Monte Carlo simulations were performed to develop the initial dosing regimen with the highest probability of therapeutic target attainment.

RESULTS

A total of 144 vancomycin concentrations were derived from 63 neonates in the external population. Five of the 28 models retained for evaluation were found predictive with a bias of 15% and an imprecision of 30%. Overall, the Grimsley and Thomson model performed best, with a bias of -0.8% and an imprecision of 20.9%; therefore, it was applied in the simulations. A novel initial dosing regimen of 15 mg/kg, followed by 11 mg/kg every 8 hours should favor therapeutic target attainment.

CONCLUSIONS

A predictive population pharmacokinetic model of vancomycin was identified after an external evaluation and used to recommend a novel initial dosing regimen. The implementation of these model-based tools may guide physicians in selecting the most appropriate initial vancomycin dose, leading to improved clinical outcomes.

Population Pharmacokinetics of Cabozantinib in Metastatic Renal Cell Carcinoma Patients: Towards Drug Expenses Saving Regimens

INTRODUCTION

Cabozantinib is one of the preferred treatment options in the latest metastatic renal cell carcinoma (mRCC) guidelines. Cabozantinib is also associated with high drug expenses irrespective of the used dose, because a flat-prizing model has been implemented. In addition, concomitant intake with a high-fat meal increases its bioavailability on average by 57%. Combined with the long terminal half-life of cabozantinib (99 h), this creates possibilities to extend the dosing interval to reduce drug expenses whilst maintaining equivalent exposure.

OBJECTIVES

The primary objective was to evaluate the population pharmacokinetic (POPPK) model of cabozantinib developed for its registration using real-world patients' therapeutic drug monitoring (TDM) data. The secondary objective was to design, simulate, and evaluate alternative dose regimens with the aim to reduce drug expenses whilst maintaining comparable exposure.

METHODS

Retrospective TDM data from mRCC patients treated with cabozantinib were obtained. The data were evaluated using the published Food and Drug Administration (FDA) cabozantinib POPPK model, a two-compartment disposition model with a dual (fast and slow) lagged first-order absorption process derived from FDA registration documents, as a basis. Subsequently, simulations of alternative drug expenses saving regimens were evaluated.

RESULTS

Twenty-seven mRCC patients with 75 pharmacokinetic observations were included. Patients were treated for a median of 75 days with a median dose of 40 mg. Model evaluation results showed that the cabozantinib TDM concentrations were adequately predicted by the published FDA cabozantinib POPPK model, except for a slightly higher clearance (CL) of 3.11 L/h compared to the reported value (2.23 L/h). The simulation study indicated that an alternative dose regimen that consists of taking 60 mg of cabozantinib for 2 days and then skipping 1 day results in comparable average exposure when compared with a 40 mg daily dose, both without food interaction, while saving 33.3% of the total drug expenses per month. The food effect of a high-fat meal was also taken into account when simulating other alternative dose regimens; 40 mg every 72 h combined with a high-fat meal resulted in comparable exposure when compared with a 20 mg daily dose fasted, while saving 66.7% in drug expenses.

CONCLUSIONS

In this study, the optimized cabozantinib POPPK model resulted in adequate prediction of real-world cabozantinib pharmacokinetic data. Alternative dosing regimens with and without using known food interactions were proposed that resulted in potential strategies to significantly reduce cabozantinib drug expenses.

Meropenem Disposition in Neonatal and Pediatric Extracorporeal Membrane Oxygenation and Continuous Renal Replacement Therapy

This study aimed to characterize the impact of extracorporeal membrane oxygenation (ECMO) on the pharmacokinetics (PK) of meropenem in neonates and children and to provide recommendations for meropenem dosing in this specific population of patients. Therapeutic drug monitoring (152 meropenem plasma concentrations) data from 45 patients (38 received ECMO) with a body weight (BW) of 7.88 (3.62-11.97) kg (median (interquartile range)) and postnatal age of 3 (0-465) days were collected. The population PK analysis was performed using NONMEM V7.3.0. Monte Carlo simulations were performed to assess the probability of target achievement (PTA) for 40% of time the free drug remained above the minimum inhibitory concentration (fT > MIC) and 100% fT > MIC. BW was found to be a significant covariate for the volume of distribution (Vd) and clearance (CL). Additionally, continuous renal replacement therapy (CRRT) was associated with a two-fold increase in Vd. In the final model, the CL and Vd for a typical patient with a median BW of 7.88 kg that was off CRRT were 1.09 L/h (RSE = 8%) and 3.98 L (14%), respectively. ECMO did not affect meropenem PK, while superimposed CRRT significantly increased Vd. We concluded that current dosing regimens provide acceptably high PTA for MIC ≤ 4 mg/L for 40% fT > MIC, but individual dose adjustments are needed for 100% fT > MIC.

A systematic evaluation of population pharmacokinetic models for polymyxin B in patients with liver and/or kidney dysfunction

Polymyxin B (PMB) is considered a last-line treatment for multidrug-resistant (MDR) gram-negative bacterial infections. Model-informed precision dosing with population pharmacokinetics (PopPK) models could help to individualize PMB dosing regimens and improve therapy. However, the external prediction ability of the established PopPK models has not been fully elaborated. This study aimed to systemically evaluate eleven PMB PopPK models from ten published literature based on a new independent population, which was divided into four different populations, patients with liver dysfunction, kidney dysfunction, liver and kidney dysfunction, and normal liver and kidney function. The whole data set consisted of 146 patients with 391 PMB concentrations. The prediction- and simulation-based diagnostics and Bayesian forecasting were conducted to evaluate model predictability. In the overall evaluation process, none of the models exhibited satisfactory predictive ability in both prediction- and simulation-based diagnostic simultaneously. However, the evaluation of the models in the subgroup of patients with normal liver and kidney function revealed improved predictive performance compared to those with liver and/or kidney dysfunction. Bayesian forecasting demonstrated enhanced predictability with the incorporation of two to three prior observations. The external evaluation highlighted a lack of consistency between the prediction results of published models and the external validation dataset. Nonetheless, Bayesian forecasting holds promise in improving the predictive performance of the models, and feedback from therapeutic drug monitoring is crucial in optimizing individual dosing regimens.

External Evaluation of Population Pharmacokinetic Models for High-Dose Methotrexate in Adult Patients with Hematological Tumors

Currently, numerous population pharmacokinetic (popPK) models for methotrexate (MTX) have been published for estimating PK parameters and variability. However, it is unclear whether the accuracy of these models is sufficient for clinical application. The aim of this study is to evaluate published models and assess their predictive performance according to the standards of scientific research. A total of 237 samples from 74 adult patients who underwent high-dose MTX (HDMTX) treatment at Shanghai Changzheng Hospital were collected. The software package NONMEM was used to perform an external evaluation for each model, including prediction-based diagnosis, simulation-based diagnosis, and Bayesian forecasting. The simulation-based diagnosis includes normalized prediction distribution error (NPDE) and visual predictive check (VPC). Following screening, 7 candidate models suitable for external validation were identified for comparison. However, none of these models exhibited excellent predictive performance. Bayesian simulation results indicated that the prediction precision and accuracy of all models significantly improved when incorporating prior concentration information. The published popPK models for MTX exhibit significant differences in their predictive performance, and none of the models were able to accurately predict MTX concentrations in our data set. Therefore, before adopting any model in clinical practice, extensive evaluation should be conducted.

Extending the code in the open-source saemix package to fit joint models of longitudinal and time-to-event data

BACKGROUND AND OBJECTIVE

Joint modeling of longitudinal and time-to-event data has gained attention over recent years with extensive developments including nonlinear models for longitudinal outcomes and flexible time-to-event models for survival outcomes, possibly involving competing risks. However, in popular software such as R, the function used to describe the biomarker dynamic is mainly linear in the parameters, and the survival submodel relies on pre-implemented functions (exponential, Weibull, ...). The objective of this work is to extend the code from the saemix package (version 3.1 on CRAN) to fit parametric joint models where longitudinal submodels are not necessary linear in their parameters, with full user control over the model function.

METHODS

We used the saemix package, designed to fit nonlinear mixed-effects models (NLMEM) through the Stochastic Approximation Expectation Maximization (SAEM) algorithm, and extended the main functions to joint model estimation. To compute standard errors (SE) of parameter estimates, we implemented a recently developed stochastic algorithm. A simulation study was proposed to assess (i) the performances of parameter estimation, (ii) the SE computation and (iii) the type I error when testing independence between the two submodels. Four joint models were considered in the simulation study, combining a linear or nonlinear mixed-effects model for the longitudinal submodel, with a single terminal event or a competing risk model.

RESULTS

For all simulation scenarios, parameters were precisely and accurately estimated with low bias and uncertainty. For complex joint models (with NLMEM), increasing the number of chains of the algorithm was necessary to reduce bias, but earlier censoring in the competing risk scenario still challenged the estimation. The empirical SE of parameters obtained over all simulations were very close to those computed with the stochastic algorithm. For more complex joint models (involving NLMEM), some estimates of random effects variances had higher uncertainty and their SE were moderately under-estimated. Finally, type I error was controlled for each joint model.

CONCLUSIONS

saemix is a flexible open-source package and we adapted it to fit complex parametric joint models that may not be estimated using standard tools. Code and examples to help users get started are freely available on Github.

When will the Glomerular Filtration Rate in Former Preterm Neonates Catch up with Their Term Peers?

AIMS

Whether and when glomerular filtration rate (GFR) in preterms catches up with term peers is unknown. This study aims to develop a GFR maturation model for (pre)term-born individuals from birth to 18 years of age. Secondarily, the function is applied to data of different renally excreted drugs.

METHODS

We combined published inulin clearance values and serum creatinine (Scr) concentrations in (pre)term born individuals throughout childhood. Inulin clearance was assumed to be equal to GFR, and Scr to reflect creatinine synthesis rate/GFR. We developed a GFR function consisting of GFRbirth (GFR at birth), and an Emax model dependent on PNA (with GFRmax, PNA50 (PNA at which half ofGFR max is reached) and Hill coefficient). The final GFR model was applied to predict gentamicin, tobramycin and vancomycin concentrations.

RESULT

In the GFR model, GFRbirth varied with birthweight linearly while in the PNA-based Emax equation, GA was the best covariate for PNA50, and current weight for GFRmax. The final model showed that for a child born at 26 weeks GA, absolute GFR is 18%, 63%, 80%, 92% and 96% of the GFR of a child born at 40 weeks GA at 1 month, 6 months, 1 year, 3 years and 12 years, respectively. PopPK models with the GFR maturation equations predicted concentrations of renally cleared antibiotics across (pre)term-born neonates until 18 years well.

CONCLUSIONS

GFR of preterm individuals catches up with term peers at around three years of age, implying reduced dosages of renally cleared drugs should be considered below this age.

Towards precision medicine of long-acting aripiprazole through population pharmacokinetic modelling

Population pharmacokinetic (popPK) models constitute a valuable tool for characterizing the pharmacokinetic properties of once-monthly long-acting injectable aripiprazole (LAI aripiprazole) and quantifying the sources of variability in drug exposure. Our aim is to develop a popPK model of both aripiprazole and its metabolite dehydro-aripiprazole in patients treated with LAI aripiprazole, and to personalize the dosing regimen of aripiprazole across different sub-groups of patients. This is a prospective study investigating the pharmacokinetics of LAI aripiprazole. A total of 93 patients were included, 21 for model development and 71 for external model evaluation. A one-compartment model with linear absorption and elimination adequately described both aripiprazole and dehydro-aripiprazole concentrations. The weight of the patients has been shown to be the factor that most influences the absorption. However, the metabolizing phenotype for CYP2D6 and the concomitant treatment with strong inhibitors of this cytochrome have been shown to be the covariates that most influence total drug exposure. This is the first popPK model developed for LAI aripiprazole that includes aripiprazole and its main active metabolite, dehydroaripiprazole. It provides a personalized dosage recommendation that maximizes the probability of achieving optimal therapeutic concentrations and minimizes the difficulties associated with trial-and-error therapeutic strategies carried out in clinical practice.

Toward Model-Based Informed Precision Dosing of Vancomycin in Hematologic Cancer Patients: A First Step

BACKGROUND AND OBJECTIVE

There is no consensus on the optimal vancomycin dose to achieve pharmacokinetic/pharmacodynamic (PK/PD) target in patients with hematologic cancer or in hematopoietic stem cell transplant (HSCT) recipients. A 24-h area under the concentration-time curve (AUC) >400 mg*h/L must be achieved early for successful treatment of severe methicillin-resistant Staphylococcus aureus (MRSA) infections. Current nomograms derived from general population data are not sufficiently accurate to allow AUC-based model-informed precision dosing. The objective of this study was to characterize vancomycin PK in patients with hematologic cancer or in HSCT recipients and to develop a model-informed dosing tool based on PK/PD target requirements.

METHODS

Pooled retrospective and prospective vancomycin serum concentrations were analyzed using NONMEM® to evaluate the performance of previously published population PK (popPK) models built from hematologic cancer datasets and to develop a novel Bayesian PK model. Patients' characteristics and clinical data were tested as potential covariates. The popPK model was validated internally and externally. Predictions of vancomycin concentrations for different dosing regimens were made using Monte-Carlo simulations, and a nomogram strategy was proposed according to selected probability of target attainment (PTA).

RESULTS

The predictive performance of the published popPK models was found to be suboptimal for our population. A novel popPK model was developed using 240 vancomycin concentrations (60 patients). A two-compartment structural model with an additive error model best described the data. Ideal body weight and estimated glomerular filtration rate (eGFR) [Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI)] were selected as covariates for volume of distribution (V) and clearance (CL). Bootstrapping confirmed the stability and precision of the popPK parameters. The volume of distribution was V1 = 46.8 L and V2 = 56.1 L, while CL = 5.63 L/h. External validation using 107 vancomycin concentrations (24 patients) demonstrated the predictivity of the model. A nomogram was developed to reach minimally PTA >50% for 400 < AUC < 600 mg*h/L.

CONCLUSION

To our knowledge, this study provides the first model-informed AUC-based strategy in North American hematologic cancer patients with or without HSCT. The resulting nomogram generated provides a simplified approach to improving the accuracy of initial vancomycin dosing in this population.

Clinical pharmacokinetics and pharmacodynamics of ivosidenib in Chinese patients with relapsed or refractory IDH1-mutated acute myeloid leukemia

PURPOSE

This study aimed to assess the pharmacokinetics (PK) and pharmacodynamics (PD) of ivosidenib in Chinese patients with relapsed or refractory acute myeloid leukemia (R/R AML) carrying the mIDH1 mutation.

METHODS

A bridging study (NCT04176393) was conducted involving 29 Chinese patients who received a daily dose of ivosidenib 500 mg in 28-day cycles. Plasma concentrations of ivosidenib and D-2-hydroxyglutarate (2-HG) were measured before and after treatment. Non-compartmental analysis (NCA) was employed to evaluate the PK, and an established population pharmacokinetic (popPK) model developed from non-Chinese patients was externally validated.

RESULTS

The findings revealed comparable PD effects of ivosidenib in Chinese patients with mIDH1 R/R AML. After adjusting for concomitant drug effects, PK characteristics were similar between Chinese and non-Chinese patients. Furthermore, the popPK model offered additional insights into the possible causes of the apparent ethnic difference in PK exposure.

CONCLUSION

The study indicates that ivosidenib can be used effectively in Chinese patients, and the observed ethnic differences in PK exposure can be explained by concomitant drug effects. The popPK model contributes to a better understanding and optimization of personalized dosing in Chinese patients with mIDH1 R/R AML.

Ampicillin Pharmacokinetics in Peripartum and Laboring Women

OBJECTIVE

Ampicillin is used for multiple peripartum indications including prevention of neonatal group B streptococcus (GBS) and treatment of chorioamnionitis. Despite its widespread use in obstetrics, existing pharmacokinetic data for ampicillin do not address contemporary indications or dosing paradigms for this population. We sought to characterize the pharmacokinetic profile of ampicillin administered to laboring women.

STUDY DESIGN

Using whole blood dried blood spot sampling technique, maternal blood samples were collected at specified times from 31 women receiving intravenous (IV) ampicillin for peripartum indications. Women received either a 2-g loading dose with 1-g administered every 4 hours (GBS) or 2-g every 6 hours (chorioamnionitis). Pharmacokinetics were analyzed via a population approach with nonlinear mixed-effect modeling.

RESULTS

The data were best described by a two-compartment model with first-order elimination, with the following whole blood parameters: central volume of distribution (V1), 75.2 L (95% confidence interval [CI]: 56.3-93.6); clearance (CL), 82.4 L/h (95% CI: 59.7-95.7); intercompartmental clearance (Q), 20.9 L/h (95% CI: 16.2-38.2); and peripheral volume of distribution (V2), 61.1 L (95% CI: 26.1-310.5). Interpatient variation in CL and V1 was large (42.0 and 56.7%, respectively). Simulations of standard dosing strategies demonstrated over 98% of women are predicted to achieve an estimated free plasma concentration above mean inhibitory concentration (MIC) of 0.5 μg/mL for more than 50% of the dosing interval.

CONCLUSION

Although large variation in the pharmacokinetics of ampicillin in pregnant women exists, as predicted by our model, current standard dosing strategies achieve adequate exposure for GBS in nearly all patients.

KEY POINTS

· IV ampicillin is widely used in obstetrics.. · Pharmacokinetic studies are lacking.. · Ampicillin pharmacokinetics were established.. · Ampicillin clearance and volume of distribution are high.. · Current ampicillin dosing strategies are sufficient..

How to Dose Vancomycin in Overweight and Obese Patients with Varying Renal (Dys)function in the Novel Era of AUC 400-600 mg·h/L-Targeted Dosing

BACKGROUND AND OBJECTIVE

The latest vancomycin guideline recommends area under the curve (AUC)-targeted dosing and monitoring for efficacy and safety. However, guidelines for AUC-targeted starting dosing in patients with obesity and/or renal insufficiency are currently lacking. This study quantifies the pharmacokinetics (PK) of vancomycin in this population and provides AUC-targeted dosing recommendations.

METHODS

Vancomycin concentrations (n = 1188) from therapeutic drug monitoring of 210 overweight and obese patients with varying degrees of renal (dys)function from the ward (74.8%) and intensive care unit (ICU, 25.2%) were pooled with published rich concentration-time data (n = 207) from 20 (morbidly) obese subjects undergoing bariatric surgery. A population model was developed using NONMEM 7.4. Stochastic simulations were performed to design dosing guidelines targeting an AUC24 between 400-600 mg·h/L.

RESULTS

Vancomycin clearance (CL) was found to increase linearly with total bodyweight and with renal function (CKD-EPI) in a power relation. Additionally, CL proved 15.5% lower in ICU patients. Our model shows that, to reach the target AUC between 400 and 600 mg·h/L in the first 48 h, two loading doses are required for both continuous infusion and intermittent dosing regimens. Maintenance doses were found to require adjustment for total bodyweight, renal function, and ICU admission status. With this guideline, the median AUC24 is well within the target from the start of the treatment onwards.

CONCLUSIONS

To achieve safe and effective vancomycin exposure for maintenance doses in overweight and obese patients, renal function, total bodyweight, and ICU admission status should be taken into account.

TRIAL REGISTRATION

The AMIGO trial was registered in the Dutch Trial Registry [NTR6058].

New insights into the role of VKORC1 polymorphisms for optimal warfarin dose selection in Caribbean Hispanic patients through an external validation of a population PK/PD model

Warfarin, an oral anticoagulant, has been used for decades to prevent thromboembolic events. The complex interplay between CYP2C9 and VKORC1 genotypes on warfarin PK and PD properties is not fully understood in special sub-groups of patients. This study aimed to externally validate a population pharmacokinetic/pharmacodynamic (PK/PD) model for the effect of warfarin on international normalized ratio (INR) and to evaluate optimal dosing strategies based on the selected covariates in Caribbean Hispanic patients. INR, and CYP2C9 and VKORC1 genotypes from 138 patients were used to develop a population PK/PD model in NONMEM. The structural definition of a previously published PD model for INR was implemented. A numerical evaluation of the parameter-covariate relationship was performed. Simulations were conducted to determine optimal dosing strategies for each genotype combinations, focusing on achieving therapeutic INR levels. Findings revealed elevated IC50 for G/G, G/A, and A/A VKORC1 haplotypes (11.76, 10.49, and 9.22 mg/L, respectively), in this population compared to previous reports. The model-guided dosing analysis recommended daily warfarin doses of 3-5 mg for most genotypes to maintain desired INR levels, although subjects with combination of CYP2C9 and VKORC1 genotypes * 2/* 2-, * 2/* 3- and * 2/* 5-A/A would require only 1 mg daily. This research underscores the potential of population PK/PD modeling to inform personalized warfarin dosing in populations typically underrepresented in clinical studies, potentially leading to improved treatment outcomes and patient safety. By integrating genetic factors and clinical data, this approach could pave the way for more effective and tailored anticoagulation therapy in diverse patient groups.

CYP3A5*3 and CYP3A4*22 Cluster Polymorphism Effects on LCP-Tac Tacrolimus Exposure: Population Pharmacokinetic Approach

The aim of the study is to develop a population pharmacokinetic (PopPK) model and to investigate the influence of CYP3A5/CYP3A4 and ABCB1 single nucleotide polymorphisms (SNPs) on the Tacrolimus PK parameters after LCP-Tac formulation in stable adult renal transplant patients. The model was developed, using NONMEM v7.5, from full PK profiles from a clinical study (n = 30) and trough concentrations (C0) from patient follow-up (n = 68). The PK profile of the LCP-Tac formulation was best described by a two-compartment model with linear elimination, parameterized in elimination (CL/F) and distributional (CLD/F) clearances and central compartment (Vc/F) and peripheral compartment (Vp/F) distribution volumes. A time-lagged first-order absorption process was characterized using transit compartment models. According to the structural part of the base model, the LCP-Tac showed an absorption profile characterized by two transit compartments and a mean transit time of 3.02 h. Inter-individual variability was associated with CL/F, Vc/F, and Vp/F. Adding inter-occasion variability (IOV) on CL/F caused a statistically significant reduction in the model minimum objective function MOFV (p < 0.001). Genetic polymorphism of CYP3A5 and a cluster of CYP3A4/A5 SNPs statistically significantly influenced Tac CL/F. In conclusion, a PopPK model was successfully developed for LCP-Tac formulation in stable renal transplant patients. CYP3A4/A5 SNPs as a combined cluster including three different phenotypes (high, intermediate, and poor metabolizers) was the most powerful covariate to describe part of the inter-individual variability associated with apparent elimination clearance. Considering this covariate in the initial dose estimation and during the therapeutic drug monitoring (TDM) would probably optimize Tac exposure attainments.

Modeling the protein binding non-linearity in population pharmacokinetic model of valproic acid in children with epilepsy: a systematic evaluation study

Background: Several studies have investigated the population pharmacokinetics (popPK) of valproic acid (VPA) in children with epilepsy. However, the predictive performance of these models in the extrapolation to other clinical environments has not been studied. Hence, this study evaluated the predictive abilities of pediatric popPK models of VPA and identified the potential effects of protein binding modeling strategies. Methods: A dataset of 255 trough concentrations in 202 children with epilepsy was analyzed to assess the predictive performance of qualified models, following literature review. The evaluation of external predictive ability was conducted by prediction- and simulation-based diagnostics as well as Bayesian forecasting. Furthermore, five popPK models with different protein binding modeling strategies were developed to investigate the discrepancy among the one-binding site model, Langmuir equation, dose-dependent maximum effect model, linear non-saturable binding equation and the simple exponent model on model predictive ability. Results: Ten popPK models were identified in the literature. Co-medication, body weight, daily dose, and age were the four most commonly involved covariates influencing VPA clearance. The model proposed by Serrano et al. showed the best performance with a median prediction error (MDPE) of 1.40%, median absolute prediction error (MAPE) of 17.38%, and percentages of PE within 20% (F20, 55.69%) and 30% (F30, 76.47%). However, all models performed inadequately in terms of the simulation-based normalized prediction distribution error, indicating unsatisfactory normality. Bayesian forecasting enhanced predictive performance, as prior observations were available. More prior observations are needed for model predictability to reach a stable state. The linear non-saturable binding equation had a higher predictive value than other protein binding models. Conclusion: The predictive abilities of most popPK models of VPA in children with epilepsy were unsatisfactory. The linear non-saturable binding equation is more suitable for modeling non-linearity. Moreover, Bayesian forecasting with prior observations improved model fitness.

Towards precision dosing of aripiprazole in children and adolescents with autism spectrum disorder: Linking blood levels to weight gain and effectiveness

AIMS

Aripiprazole is one of the most commonly prescribed antipsychotic drugs to children and adolescents worldwide, but it is associated with serious side-effects, including weight gain. This study assessed the population pharmacokinetics of aripiprazole and its active metabolite and investigated the relationship between pharmacokinetic parameters and body mass index (BMI) in children and adolescents with autism spectrum disorder (ASD) and behavioural problems. Secondary outcomes were metabolic, endocrine, extrapyramidal and cardiac side-effects and drug effectiveness.

METHODS

Twenty-four children and adolescents (15 males, 9 females) aged 6-18 years were included in a 24-week prospective observational trial. Drug plasma concentrations, side-effects and drug effectiveness were measured at several time points during follow-up. Relevant pharmacokinetic covariates, including CYP2D6, CYP3A4, CYP3A5 and P-glycoprotein (ABCB1) genotypes, were determined. Nonlinear mixed-effects modelling (NONMEM®) was used for a population pharmacokinetic analysis with 92 aripiprazole and 91 dehydro-aripiprazole concentrations. Subsequently, model-based trough concentrations, maximum concentrations and 24-h area under the curves (AUCs) were analysed to predict outcomes using generalized and linear mixed-effects models.

RESULTS

For both aripiprazole and dehydro-aripiprazole, one-compartment models best described the measured concentrations, with albumin and BMI as significant covariates. Of all the pharmacokinetic parameters, higher sum (aripiprazole plus dehydro-aripiprazole) trough concentrations best predicted higher BMI z-scores (P < .001) and higher Hb1Ac levels (P = .03) during follow-up. No significant association was found between sum concentrations and effectiveness.

CONCLUSIONS

Our results indicate a threshold with regard to safety, which suggests that therapeutic drug monitoring of aripiprazole could potentially increase safety in children and adolescents with ASD and behavioural problems.

Towards Model-Informed Precision Dosing of Voriconazole: Challenging Published Voriconazole Nonlinear Mixed-Effects Models with Real-World Clinical Data

BACKGROUND AND OBJECTIVES

Model-informed precision dosing (MIPD) frequently uses nonlinear mixed-effects (NLME) models to predict and optimize therapy outcomes based on patient characteristics and therapeutic drug monitoring data. MIPD is indicated for compounds with narrow therapeutic range and complex pharmacokinetics (PK), such as voriconazole, a broad-spectrum antifungal drug for prevention and treatment of invasive fungal infections. To provide guidance and recommendations for evidence-based application of MIPD for voriconazole, this work aimed to (i) externally evaluate and compare the predictive performance of a published so-called 'hybrid' model for MIPD (an aggregate model comprising features and prior information from six previously published NLME models) versus two 'standard' NLME models of voriconazole, and (ii) investigate strategies and illustrate the clinical impact of Bayesian forecasting for voriconazole.

METHODS

A workflow for external evaluation and application of MIPD for voriconazole was implemented. Published voriconazole NLME models were externally evaluated using a comprehensive in-house clinical database comprising nine voriconazole studies and prediction-/simulation-based diagnostics. The NLME models were applied using different Bayesian forecasting strategies to assess the influence of prior observations on model predictivity.

RESULTS

The overall best predictive performance was obtained using the aggregate model. However, all NLME models showed only modest predictive performance, suggesting that (i) important PK processes were not sufficiently implemented in the structural submodels, (ii) sources of interindividual variability were not entirely captured, and (iii) interoccasion variability was not adequately accounted for. Predictive performance substantially improved by including the most recent voriconazole observations in MIPD.

CONCLUSION

Our results highlight the potential clinical impact of MIPD for voriconazole and indicate the need for a comprehensive (pre-)clinical database as basis for model development and careful external model evaluation for compounds with complex PK before their successful use in MIPD.

Population pharmacokinetics of oxcarbazepine active metabolite in Chinese children with epilepsy

Oxcarbazepine (OXC) is an antiepileptic drug whose efficacy is largely attributed to its monohydroxy derivative metabolite (MHD). Nevertheless, there exists significant inter-individual variability in both the pharmacokinetics and therapeutic response of this drug. The objective of this study is to explore the impact of patients' characteristics and genetic variants on MHD clearance in a population pharmacokinetic (PPK) model of Chinese pediatric patients with epilepsy. The PPK model was developed using a nonlinear mixed effects modeling method based on 231 MHD plasma concentrations obtained from 185 children with epilepsy. The one-compartment model and combined residual model were established to describe the pharmacokinetics of MHD. Forward addition and backward elimination were employed to evaluate the impact of covariates on the model parameters. The model was evaluated using goodness-of-fit, bootstrap, visual predictive checks, and normalized prediction distribution errors. In the two final PPK models, age, estimated glomerular filtration rate (eGFR), and a combined genotype of six variants (rs1045642, rs2032582, rs7668282, rs2396185, rs2304016, rs1128503) were found to significantly reduce inter-individual variability for MHD clearance. The inter-individual clearance equals to 1.38 × (Age/4.74)0.29 × (eGFR/128.66)0.25 × eθABCB-UGT-SCN-INSR for genetic variants included model and 1.30 × (Age/4.74)0.30 × (eGFR/128.66)0.23 for model without genetic variants. The precision of all parameters was deemed acceptable, and the model exhibited good predictability while remaining stable and effective.    Conclusion: Age, eGFR, and genotype may play a significant role in MHD clearance in children with epilepsy. The developed PPK models hold potential utility in facilitating oxcarbazepine dose adjustment in pediatric patients. What is Known: • The adjustment of the oxcarbazepine regimen remains difficult due to the considerable inter- and intra-individual variability of oxcarbazepine pharmacokinetics. • Body weight and co-administration with enzyme-inducing antiepileptic drugs emerge as the most influential factors contributing to the pharmacokinetics of MHD. What is New: • A positive correlation was observed between eGFR and the clearance of MHD in pediatric patients with epilepsy. • We explored the influence of genetic polymorphisms on MHD clearance and identified a combined genotype (ABCB-UGT-SCN-INSR) that exhibited a significant association with MHD concentration.

Population pharmacokinetics of topiramate in Chinese children with epilepsy

OBJECTIVE

Topiramate, a broad-spectrum antiepileptic drug, exhibits substantial inter-individual variability in both its pharmacokinetics and therapeutic response. The aim of this study was to investigate the influence of patient characteristics and genetic variants on topiramate clearance using population pharmacokinetic (PPK) models in a cohort of Chinese pediatric patients with epilepsy.

METHOD

The PPK model was constructed using a nonlinear mixed-effects modeling approach, utilizing a dataset comprising 236 plasma concentrations of topiramate obtained from 181 pediatric patients with epilepsy. A one-compartment model combined with a proportional residual model was employed to characterize the pharmacokinetics of topiramate. Covariate analysis was performed using forward addition and backward elimination to assess the influence of covariates on the model parameters. The model was thoroughly evaluated through goodness-of-fit analysis, bootstrap, visual predictive checks, and normalized prediction distribution errors. Monte Carlo simulations were utilized to devise topiramate dosing strategies.

RESULT

In the final PPK models of topiramate, body weight, co-administration with oxcarbazepine, and a combined genotype of GKIR1-UGT (GRIK1 rs2832407, UGT2B7 rs7439366, and UGT1A1 rs4148324) were identified as significant covariates affecting the clearance (CL). The clearance was estimated using the formulas CL (L/h) = 0.44 × (BW⁄11.7)0.82 × eOXC for the model without genetic variants and CL (L/h) = 0.49 × (BW⁄11.7)0.81 × eOXC × eGRIK1-UGT for the model incorporating genetic variants. The volume of distribution (Vd) was estimated using the formulas Vd (L) = 6.6 × (BW⁄11.7). The precision of all estimated parameters was acceptable. Furthermore, the model demonstrated good predictability, exhibiting stability and effectiveness in describing the pharmacokinetics of topiramate.

CONCLUSION

The clearance of topiramate in pediatric patients with epilepsy may be subject to the influence of factors such as body weight, co-administration with oxcarbazepine, and genetic polymorphism. In this study, PPK models were developed to better understand and account for these factors, thereby improving the precision and individualization of topiramate therapy in children with epilepsy.

Population pharmacokinetics of Amisulpride in Chinese patients with schizophrenia with external validation: the impact of renal function

Introduction: Amisulpride is primarily eliminated via the kidneys. Given the clear influence of renal clearance on plasma concentration, we aimed to explicitly examine the impact of renal function on amisulpride pharmacokinetics (PK) via population PK modelling and Monte Carlo simulations. Method: Plasma concentrations from 921 patients (776 in development and 145 in validation) were utilized. Results: Amisulpride PK could be described by a one-compartment model with linear elimination where estimated glomerular filtration rate, eGFR, had a significant influence on clearance. All PK parameters (estimate, RSE%) were precisely estimated: apparent volume of distribution (645 L, 18%), apparent clearance (60.5 L/h, 2%), absorption rate constant (0.106 h-1, 12%) and coefficient of renal function on clearance (0.817, 10%). No other significant covariate was found. The predictive performance of the model was externally validated. Covariate analysis showed an inverse relationship between eGFR and exposure, where subjects with eGFR= 30 mL/min/1.73 m2 had more than 2-fold increase in AUC, trough and peak concentration. Simulation results further illustrated that, given a dose of 800 mg, plasma concentrations of all patients with renal impairment would exceed 640 ng/mL. Discussion: Our work demonstrated the importance of renal function in amisulpride dose adjustment and provided a quantitative framework to guide individualized dosing for Chinese patients with schizophrenia.

Maturation of Paracetamol Elimination Routes in Preterm Neonates Born Below 32 Weeks of Gestation

PURPOSE

Despite being off-label, intravenous paracetamol (PCM) is increasingly used to control mild-to-moderate pain in preterm neonates. Here we aim to quantify the maturation of paracetamol elimination pathways in preterm neonates born below 32 weeks of gestation.

METHODS

Datasets after single dose (rich data) or multiple doses (sparse data) of intravenous PCM dose (median (range)) 9 (3-25) mg/kg were pooled, containing 534 plasma and 44 urine samples of PCM and metabolites (PCM-glucuronide, PCM-sulfate, PCM-cysteine, and PCM-mercapturate) from 143 preterm neonates (gestational age 27.7 (24.0-31.9) weeks, birthweight 985 (462-1,925) g, postnatal age (PNA) 5 (0-30) days, current weight 1,012 (462-1,959) g. Population pharmacokinetic analysis was performed using NONMEM® 7.4.

RESULTS

For a typical preterm neonate (birthweight 985 g; PNA 5 days), PCM clearance was 0.137 L/h, with glucuronidation, sulfation, oxidation and unchanged renal clearance accounting for 5.3%, 73.7%, 16.3% and 4.6%, respectively. Maturational changes in total PCM clearance and its elimination pathways were best described by birthweight and PNA. Between 500-1,500 g birthweight, total PCM clearance increases by 169%, with glucuronidation, sulfation and oxidation clearance increasing by 347%, 164% and 164%. From 1-30 days PNA for 985 g birthweight neonate, total PCM clearance increases by 167%, with clearance via glucuronidation and oxidation increasing by 551%, and sulfation by 69%.

CONCLUSION

Birthweight and PNA are the most important predictors for maturational changes in paracetamol clearance and its glucuronidation, sulfation and oxidation. As a result, dosing based on bodyweight alone will not lead to consistent paracetamol concentrations among preterm neonates.

External validation of population pharmacokinetic models of vancomycin in postoperative neurosurgical patients

OBJECTIVE

Vancomycin is commonly used in the prevention and treatment of intracranial infections in postoperative neurosurgical patients with narrow therapeutic window and large pharmacokinetic variations. Several population pharmacokinetic (PPK) models of vancomycin have been established for neurosurgical patients. But comprehensive external evaluation has not been performed for almost all models. The objective of this study was to evaluate the predictive ability of published vancomycin PPK models in adult postoperative neurosurgical patients using an independent dataset.

METHOD

PubMed, Embase and China National Knowledge Internet databases were searched to identify published vancomycin PPK models in adult postoperative neurosurgical patients. Prediction-based and simulation-based diagnostics were used to evaluate model predictability. Bayesian forecasting was used to assess the influence of prior concentration on model prediction performance.

RESULT

A total of 763 vancomycin plasma concentrations from 493 postoperative neurosurgical patients were included in the external dataset. Eight population pharmacokinetic models of vancomycin in postoperative neurosurgical patients were included and evaluated. The model by Zhang et al. exhibited the best predictive performance in prediction-based diagnostics and prediction-corrected visual predictive checks, followed by the model by Shen et al. The predictive performance of other models was not satisfactory. The normalized predictive distribution error test shows that none of the models is suitable to describe our data. The predictive performance of vancomycin models was obviously improved by maximum a posteriori Bayesian forecasting.

CONCLUSION

The published PPK models for adult postoperative neurosurgical patients show extensive variation in predictive performance in our patients. Although it is challenging to recommend initial doses of vancomycin from these predictive models, the combination of model-based prediction and therapeutic drug monitoring can be used for dose optimization.

Optimal Teicoplanin Dosage Regimens in Critically Ill Patients: Population Pharmacokinetics and Dosing Simulations Based on Renal Function and Infection Type

Purpose

To develop a population pharmacokinetic model describing teicoplanin concentrations in patients hospitalized in intensive care unit (ICU) and to perform Monte Carlo simulations to provide detailed dosing regimens of teicoplanin.

Methods

This single-center, prospective, observational study was conducted on 151 patients in ICU with 347 plasma samples. The population pharmacokinetics model was established and various covariates were evaluated. The probability of target attainment (PTA) of various proposal dosing regimens was calculated by Monte Carlo simulations.

Results

The two-compartment model adequately described teicoplanin concentration-time data. The estimated glomerular filtration rate (eGFR) associated with systemic clearance (CL) was the only covariate included in the final model. The estimate of CL was 0.838 L/h, with the eGFR adjustment factor of 0.00823. The volume of the central compartment (Vc), inter-compartmental clearance (Q) and volumes of the peripheral compartments (Vp) were 14.4 L, 3.08 L/h and 51.6 L, respectively. The simulations revealed that the standard dosage regimen was only sufficient for the patients with severe renal dysfunction (eGFR ≤ 30 mL/min/1.73 m2) to attain target trough concentration (Cmin, PTA 52.8%). When eGFR > 30 mL/min/1.73 m2, increasing dose and the administration times of loading doses were the preferred options to achieve target Cmin based on the renal function and types of infection.

Conclusion

The most commonly used standard dosage regimen was insufficient for all ICU patients. Our study provided detailed dosing regimens of teicoplanin stratified by eGFR and types of infection for ICU patients.

Population Pharmacokinetics of Posaconazole in Immune-Compromised Children and Assessment of Target Attainment in Invasive Fungal Disease

BACKGROUND AND OBJECTIVE

Posaconazole (PSZ) is a triazole antifungal for the management of invasive fungal disease (IFD) in adults and children. Although PSZ is available as an intravenous (IV) solution, oral suspension (OS) and delayed-release tablets (DRTs), OS is the preferred formulation for pediatric use because of potential safety concerns associated with an excipient in the IV formulation and difficulty in swallowing intact tablets by children. However, poor biopharmaceutical characteristics of the OS formulation leads to an unpredictable dose-exposure profile of PSZ in children, potentially risking therapeutic failure. The goal of this study was to characterize the population pharmacokinetics (PK) of PSZ in immunocompromised children and assess therapeutic target attainment.

METHODS

Serum concentrations of PSZ were collected retrospectively from records of hospitalized patients. A population PK analysis was performed in a nonlinear mixed-effects modeling framework with NONMEM (v7.4). The PK parameters were scaled to body weight, then potential covariate effects were assessed. The final PK model was used to evaluate recommended dosing schemes through simulation of target attainment (as a percentage of the population having steady-state trough concentrations above the recommended target) using Simulx (v2021R1).

RESULTS

Repeated measurement data of 202 serum concentrations of total PSZ were acquired from 47 immunocompromised patients between 1 and 21 years of age receiving PSZ either intravenously or orally, or both. A one-compartment PK model with first-order absorption and linear elimination best fit the data. The estimated absolute bioavailability (95% confidence interval) for suspension (Fs) was 16% (8-27%), which was significantly lower than the reported tablet bioavailability (Ft) [67%]. Fs was reduced by 62% and 75% upon concomitant administration with pantoprazole (PAN) and omeprazole (OME), respectively. Famotidine resulted in a reduction of Fs by only 22%. Both fixed dosing and weight-based adaptive dosing provided adequate target attainment when PAN or OME were not coadministered with the suspension.

CONCLUSIONS

The results of this study revealed that both fixed and weight-based adaptive dosing schemes can be appropriate for target attainment across all PSZ formulations, including suspension. Additionally, covariate analysis suggests that concomitant proton pump inhibitors should be contraindicated during PSZ suspension dosing.

Vancomycin population pharmacokinetics analysis in Chinese paediatric patients with varying degrees of renal function and ages: development of new practical dosing recommendations

OBJECTIVES

To describe the pharmacokinetics of vancomycin in a large Chinese paediatric cohort with varying degrees of renal function and ages and to develop practical dosing guidelines.

PATIENTS AND METHODS

We conducted a retrospective population pharmacokinetic study using data from paediatric patients who received vancomycin between June 2013 and June 2022. A non-linear mixed-effect modelling approach with a one-compartment model structure was applied. Monte Carlo simulations were used to stimulate an optimal dosage regimen to achieve the target of AUC24/MIC between 400 and 650.

RESULTS

We analysed a total of 673 paediatric patients and 1547 vancomycin serum concentrations. Covariate analysis revealed that physiological maturation, renal function, albumin and cardiothoracic surgery (CTS) significantly affected vancomycin pharmacokinetics. The typical clearance and volume of distribution, standardized to 70 kg, were 7.75 L/h (2.3% relative standard error, RSE) and 36.2 L (1.7% RSE), respectively. Based on the model, we proposed an optimal dosing regimen that considers the patient's age and estimate glomerular filtration rate (eGFR) to achieve a target AUC24/MIC for CTS and non-CTS patients. We also found that a loading dose of 20 mg/kg can help patients with an eGFR of <60 mL/min/1.73 m2 achieve the target AUC on the first day of treatment.

CONCLUSIONS

We established vancomycin pharmacokinetic parameters in Chinese paediatric patients and proposed a dosing guideline integrating eGFR, age and CTS status, potentially improving clinical outcomes and reducing nephrotoxicity risk.

Comparison of Three Renal Function Formulas for Ganciclovir/Valganciclovir Dose Individualization in CMV-Infected Solid Organ Transplantation Patients Using a Population Approach

BACKGROUND AND OBJECTIVE

The gold standard treatment of established cytomegalovirus infection or prevention in solid organ transplantation is the intravenous administration of ganciclovir (GCV) or oral administration of valganciclovir (VGCV), both adjusted to the renal function. In both instances, there is a high interindividual pharmacokinetic variability, mainly owing to the wide range of variation of both the renal function and body weight. Therefore, accurate estimation of the renal function is crucial for GCV/VGCV dose optimization. This study aimed to compare three different formulas for estimating the renal function in solid organ transplantation patients with cytomegalovirus infection, for individualizing antiviral therapy with GCV/VGCV, using a population approach.

METHODS

A population pharmacokinetic analysis was performed using NONMEM 7.4. A total of 650 plasma concentrations obtained after intravenous GCV and oral VGCV administrations were analyzed, from intensive and sparse sampling designs. Three different population pharmacokinetic models were built with the renal function given by Cockcroft-Gault, Modification of Diet in Renal Disease, or Chronic Kidney Disease EPIdemiology Collaboration (CKD-EPI) formulas. Pharmacokinetic parameters were allometrically scaled to body weight.

RESULTS

The CKD-EPI formula was identified as the best predictor of between-patient variability in GCV clearance. Internal and external validation techniques showed that the CKD-EPI model had better stability and performed better compared with the others.

CONCLUSIONS

The model based on the more accurate estimation of the renal function with the CKD-EPI formula and body weight as a size metric most used in the clinical practice can refine initial dose recommendations and contribute to GCV and VGCV dose individualization when required in the prevention or treatment of cytomegalovirus infection in solid organ transplantation patients.

Intraperitoneally Administered Vancomycin in Patients with Peritoneal Dialysis-Associated Peritonitis: Population Pharmacokinetics and Dosing Implications

Peritonitis is a limiting complication of peritoneal dialysis, which is treated by intraperitoneal administration of antibiotics. Various dosing strategies are recommended for intraperitoneally administered vancomycin, which leads to large differences in intraperitoneal vancomycin exposure. Based on data from therapeutic drug monitoring, we developed the first-ever population pharmacokinetic model for intraperitoneally administered vancomycin to evaluate intraperitoneal and plasma exposure after dosing schedules recommended by the International Society for Peritoneal Dialysis. According to our model, currently recommended dosing schedules lead to possible underdosing of a large proportion of patients. To prevent this, we suggest avoiding intermittent intraperitoneal vancomycin administration, and for the continuous dosing regimen, we suggest a loading dose of 20 mg/kg followed by maintenance doses of 50 mg/L in each dwell to improve the intraperitoneal exposure. Vancomycin plasma level measurement on the fifth day of treatment with subsequent dose adjustment would prevent it from reaching toxic levels in the few patients who are susceptible to overdose.

Vancomycin Clearance in Obese Adults is not Predictive of Clearance in Obese Adolescents

BACKGROUND AND OBJECTIVE

Contradictory pharmacokinetic (PK) results have been observed between obese adults and obese adolescents, with absolute clearance (CL) reported to be either unaltered, lower, or higher in obese adolescents compared to obese adults. This study investigates the PK of vancomycin in adolescents and adults who are overweight or obese.

METHODS

Data from 125 overweight and obese adolescents (aged 10-18 years, weight 28.3-188 kg) and 81 overweight and obese adults (aged 29-88 years, weight 66.7-143 kg) were analysed using population PK modelling. In addition to age, sex, renal function estimates, and regular weight descriptors, we evaluated standard weight (WT_standard, defined as weight for length, age, and sex in adolescents and weight for length in adults) and excess weight (WT_excess, defined as total body weight (TBW) minus WT_standard) as covariates in order to distinguish between weight resulting from length versus weight resulting from obesity.

RESULTS

Analyzing adolescents and adults together, vancomycin CL was found to increase with TBW and decrease with increasing age (p < 0.001). A covariate analysis investigating adolescents and adults separately found that vancomycin CL increased with WT_standard in adolescents and adults, albeit with different functions, with adolescents having a higher CL per WT_standard than adults. Moreover, in this separate model, adolescent males had 21% higher CL than adolescent females of the same WT_standard, while in adults, CL decreased with increasing age (p < 0.001).

CONCLUSION

There are apparent differences in vancomycin CL in overweight and obese adults versus overweight and obese adolescents, implying that dosing of vancomycin cannot be directly extrapolated between these populations.

New Recommendations of a Height-Based Dosing Regimen of Tobramycin for Cystic Fibrosis in Adults: A Population Pharmacokinetic Analysis

BACKGROUND

Acute pulmonary exacerbations (APEs) in patients with adult cystic fibrosis (CF) are treated with a beta-lactam and an aminoglycoside for activity against Pseudomonas aeruginosa (PA). Emerging drug resistance and changing pharmacokinetic profile in an aging population involve a reevaluation of tobramycin dosing recommendations. The objective of this study was to develop a population pharmacokinetic model and establish optimal dosing recommendations for tobramycin using Monte Carlo simulations.

METHODS

This retrospective clinical study and data collection were performed at the CF center of the McGill University Health Center (MUHC), Canada. Model development and simulations were performed using a nonlinear mixed-effect modeling approach (NONMEM, version 7.4.2). The ratios of maximal concentration (C max ) to the minimal inhibitory concentration (MIC) (C max /MIC ≥8 and ≥10) and area under the curve (AUC) to the MIC (AUC/MIC ≥70 and ≥100) were evaluated.

RESULTS

Adult patients with CF (n = 51) treated with tobramycin were included in the study. Plasma concentrations of tobramycin were obtained for 699 samples from the MUHC database. The two-compartmental model best described the pharmacokinetics of tobramycin. The association of patient height with the central volume of distribution significantly improved this model. Height, rather than weight, induced the best reduction in objective function. According to simulations, doses between 3.4 mg/cm and 4.4 mg/cm were necessary to achieve C max /MIC values of ≥8 and ≥10, respectively. However, higher doses were required to achieve the AUC/MIC targets.

CONCLUSIONS

This study demonstrated that height of the patients seems to be more suitable than their weight for dosing adjustments in adult patients with CF. According to this model, initial doses of tobramycin between 3.4 and 4.4 mg/cm should be recommended for patients with a median height of 164 cm and weight of 55 kg to achieve the target plasma concentrations.

Transferability of Published Population Pharmacokinetic Models for Apixaban and Rivaroxaban to Subjects with Obesity Treated for Venous Thromboembolism: A Systematic Review and External Evaluations

Apixaban and rivaroxaban have first-line use for many patients needing anticoagulation for venous thromboembolism (VTE). The pharmacokinetics of these drugs in non-obese subjects have been extensively studied, and, while changes in pharmacokinetics have been documented in obese patients, data remain scarce for these anticoagulants. The aim of this study was to perform an external validation of published population pharmacokinetic (PPK) models of apixaban and rivaroxaban in a cohort of obese patients with VTE. A literature search was conducted in the PubMed/MEDLINE, Scopus, and Embase databases following the PRISMA statement. External validation was performed using MonolixSuite software, using prediction-based and simulation-based diagnostics. An external validation dataset from the university hospitals of Brest and Rennes, France, included 116 apixaban pharmacokinetic samples from 69 patients and 121 rivaroxaban samples from 81 patients. Five PPK models of apixaban and 16 models of rivaroxaban were included, according to the inclusion criteria of the study. Two of the apixaban PPK models presented acceptable performances, whereas no rivaroxaban PPK model did. This study identified two published models of apixaban applicable to apixaban in obese patients with VTE. However, none of the rivaroxaban models evaluated were applicable. Dedicated studies appear necessary to elucidate rivaroxaban pharmacokinetics in this population.

Population Pharmacokinetics of Isavuconazole in Critical Care Patients with COVID-19-Associated Pulmonary Aspergillosis and Monte Carlo Simulations of High Off-Label Doses

Isavuconazole is a triazole antifungal agent recently recommended as first-line therapy for invasive pulmonary aspergillosis. With the COVID-19 pandemic, cases of COVID-19-associated pulmonary aspergillosis (CAPA) have been described with a prevalence ranging from 5 to 30%. We developed and validated a population pharmacokinetic (PKpop) model of isavuconazole plasma concentrations in intensive care unit patients with CAPA. Nonlinear mixed-effect modeling Monolix software were used for PK analysis of 65 plasma trough concentrations from 18 patients. PK parameters were best estimated with a one-compartment model. The mean of ISA plasma concentrations was 1.87 [1.29-2.25] mg/L despite prolonged loading dose (72 h for one-third) and a mean maintenance dose of 300 mg per day. Pharmacokinetics (PK) modeling showed that renal replacement therapy (RRT) was significantly associated with under exposure, explaining a part of clearance variability. The Monte Carlo simulations suggested that the recommended dosing regimen did not achieve the trough target of 2 mg/L in a timely manner (72 h). This is the first isavuconazole PKpop model developed for CAPA critical care patients underlying the need of therapeutic drug monitoring, especially for patients under RRT.

Alternative Pharmacokinetic Metrics in Single-Dose Studies to Ensure Bioequivalence of Prolonged-Release Products at Steady State-A Case Study

(1) Background: this article investigates which PK metrics in a single-dose study (concentration at the end of posology interval, C_τ, partial areas under the curve, pAUCs, or half-value duration, HVD) are more sensitive and less variable for predicting the failure of a prolonged-release product at steady-state that was the bioequivalent for C_max, AUC_0-t and AUC_0-inf, in the single-dose study; (2) Methods: a cross-over study was performed in 36 subjects receiving desvenlafaxine 100 mg prolonged-release tablets. Conventional (C_max, AUC_0-t and AUC_0-inf) and additional (C_τ, pAUCs and HVD) PK metrics were considered after single-dose conditions. Predicted PK metrics at steady state (AUC_0-τ, C_max,ss, and C_τ,ss) were derived using a population PK model approach; (3) Results: the existing differences in the shape of the concentration-time curves precluded to show equivalence for C_τ,ss in the simulated study at steady state. This failure to show equivalence at steady state was predicted by C_τ, pAUCs and HVD in the single-dose study. C_τ was the most sensitive metric for detecting the different shape, with a lower intra-subject variability than HVD; (4) Conclusions: conventional PK metrics for single-dose studies (C_max, AUC_0-t and AUC_0-inf) are not enough to guarantee bioequivalence at steady state for prolonged-release products.

Fractal Kinetic Implementation in Population Pharmacokinetic Modeling

Compartment modeling is a widely accepted technique in the field of pharmacokinetic analysis. However, conventional compartment modeling is performed under a homogeneity assumption that is not a naturally occurring condition. Since the assumption lacks physiological considerations, the respective modeling approach has been questioned, as novel drugs are increasingly characterized by physiological or physical features. Alternative approaches have focused on fractal kinetics, but evaluations of their application are lacking. Thus, in this study, a simulation was performed to identify desirable fractal-kinetics applications in conventional modeling. Visible changes in the profiles were then investigated. Five cases of finalized population models were collected for implementation. For model diagnosis, the objective function value (OFV), Akaike's information criterion (AIC), and corrected Akaike's information criterion (AICc) were used as performance metrics, and the goodness of fit (GOF), visual predictive check (VPC), and normalized prediction distribution error (NPDE) were used as visual diagnostics. In most cases, model performance was enhanced by the fractal rate, as shown in a simulation study. The necessary parameters of the fractal rate in the model varied and were successfully estimated between 0 and 1. GOF, VPC, and NPDE diagnostics show that models with the fractal rate described the data well and were robust. In the simulation study, the fractal absorption process was, therefore, chosen for testing. In the estimation study, the rate application yielded improved performance and good prediction-observation agreement in early sampling points, and did not cause a large shift in the original estimation results. Thus, the fractal rate yielded explainable parameters by setting only the heterogeneity exponent, which reflects true physiological behavior well. This approach can be expected to provide useful insights in pharmacological decision making.

Influence of Renal Function on Phosphoramide Mustard Exposure: A Nonlinear Mixed-Effects Analysis

Phosphoramide mustard (PM) is the final cytotoxic metabolite formed from the parent compound cyclophosphamide through a complex metabolic pathway, primarily through hepatic metabolism. Little is known about the effect of renal elimination on the disposition of PM. We evaluated the effect of renal function on PM exposure after single doses of cyclophosphamide in 85 patients undergoing allogeneic hematopoietic cell transplantation using nonlinear mixed-effects modeling. Mixed linear and nonlinear elimination pathways were required to adequately describe the disposition of PM. Creatinine clearance (CrCL) was incorporated as a covariate associated with first-order elimination, representing renal clearance (Cl_R ) of PM. For a 70-kg patient, Cl_R was 14.9 L/h, Volume of distribution was 525 L, maximum rate was 81.2 mg/h, and the concentration to achieve 50% of maximum rate was 0.51 mg/L. We conducted simulations to explore the impact of CrCL as a measure of renal function and observed that when CrCL decreases from 120 to 40 mL/min, PM area under the plasma concentration-time curve (AUC) from time 0 to 8 hours and AUC increases by 9.2% and 80.9% on average after a single dose, respectively. Our data suggest that renal function has limited influence on PM exposure during the first 8 hours after dosing but has a large impact on the total exposure. Dose adjustment of cyclophosphamide may not be necessary in hematopoietic cell transplant recipients with moderate to severe kidney dysfunction to attain targeted exposures based on AUC from time 0 to 8 hours. However, dose reduction may be necessary if demonstrated at some future time that total AUC is a better surrogate for safety or toxicity.

Population Pharmacokinetic Analysis for Model-Based Therapeutic Drug Monitoring of Tacrolimus in Chinese Han Heart Transplant Patients

BACKGROUND AND OBJECTIVE

Tacrolimus has become the first-line immunosuppressant for preventing rejection after heart transplantation. The present study aimed to investigate genetic variants and clinical factors affecting the variability of tacrolimus in Chinese Han heart transplant patients using a population pharmacokinetic approach.

METHODS

The retrospective study included 53 hospitalized patients with 547 tacrolimus concentrations for analysis. Nonlinear mixed-effects modeling was used to develop the population pharmacokinetics model for tacrolimus in patients with heart transplants, followed by Monte Carlo simulations to design initial dosing regimens.

RESULTS

In our study, the mutation rate of CYP3A4*18B (C>T) was 27.36%. An oral one-compartment model with first-order absorption and elimination was used to describe the pharmacokinetics of tacrolimus in heart transplant patients. In the final model, the estimated apparent clearance (CL/F) and volume of distribution (V/F) were 532.5 L/h [12.20% interindividual variability, IIV] and 16.87 L (23.16% IIV), respectively. Albumin, postoperative time, and rs2242480 (CYP3A4*18B) gene polymorphisms were the significant covariates affecting CL/F, and creatinine clearance had significant effects on the V/F.

CONCLUSION

The population pharmacokinetic model of tacrolimus in heart transplant patients can better estimate the population and individual pharmacokinetic parameters of patients and can provide a reference for the design of individualized dosing regimens.

Isavuconazole Pharmacokinetics and Pharmacodynamics in Children

Isavuconazole is a broad-spectrum azole anti-fungal not yet approved in children. We conducted a retrospective, single-center review of isavuconazole use and routine therapeutic drug monitoring in pediatric patients, extracting demographic, dosing, concentration, mortality and hepatoxicity data. We constructed a nonparametric population model using Pmetrics. Of 26 patients, 19 (73%) were male. The mean (SD) age and weight were 12.7 (5.5) years and 50.9 (26.8) kg. Eighty percent received between 9.7 and 10.6 mg/kg per dose. Ten (38%) subjects had proven fungal disease and eight (31%) had probable disease, mostly with Candida and Aspergillus spp. The predicted steady-state isavuconazole concentrations in our patients were similar to previous reports in children and adults, and simulations with the proposed dosing of 10 mg/kg/dose every 8 h for 2 days followed by once daily maintenance matched effective adult exposures. Attributable mortality (5 of 11 deaths) was associated with steady-state daily AUC < 60 mg∗h/L and higher AST/ALT with trough concentrations > 5 mg/L. Neither dose nor trough alone correlated well with AUC, but AUC can be estimated with one sample 10 h after the first maintenance dose or a trough concentration, if combined with a Bayesian approach or a peak and trough without a Bayesian approach.

Population Pharmacokinetic/Pharmacodynamic Modelling of Daptomycin for Schedule Optimization in Patients with Renal Impairment

The aims of this study are (i) to develop a population pharmacokinetic/pharmacodynamic model of daptomycin in patients with normal and impaired renal function, and (ii) to establish the optimal dose recommendation of daptomycin in clinical practice. Several structural PK models including linear and non-linear binding kinetics were evaluated. Monte Carlo simulations were conducted with a fixed combination of creatinine clearance (30-90 mL/min/1.73 m2) and body weight (50-100 kg). The final dataset included 46 patients and 157 daptomycin observations. A two-compartment model with first-order peripheral distribution and elimination kinetics assuming non-linear protein-binding kinetics was selected. The bactericidal effect for Gram+ strains with MIC ≤ 0.5 mg/L could be achieved with 5-12 mg/kg daily daptomycin based on body weight and renal function. The administration of 10-17 mg/kg q48 h daptomycin allows to achieve bactericidal effect for Gram+ strains with MIC ≤ 1 mg/L. Four PK samples were selected as the optimal sampling strategy for an accurate AUC estimation. A quantitative framework has served to characterize the non-linear binding kinetics of daptomycin in patients with normal and impaired renal function. The impact of different dosing regimens on the efficacy and safety outcomes of daptomycin treatment based on the unbound exposure of daptomycin and individual patient characteristics has been evaluated.

Impact of model misspecification on model-based tests in PK studies with parallel design: real case and simulation studies

This article evaluates the performance of pharmacokinetic (PK) equivalence testing between two formulations of a drug through the Two-One Sided Tests (TOST) by a model-based approach (MB-TOST), as an alternative to the classical non-compartmental approach (NCA-TOST), for a sparse design with a few time points per subject. We focused on the impact of model misspecification and the relevance of model selection for the reference data. We first analysed PK data from phase I studies of gantenerumab, a monoclonal antibody for the treatment of Alzheimer’s disease. Using the original rich sample data, we compared MB-TOST to NCA-TOST for validation. Then, the analysis was repeated on a sparse subset of the original data with MB-TOST. This analysis inspired a simulation study with rich and sparse designs. With rich designs, we compared NCA-TOST and MB-TOST in terms of type I error and study power. With both designs, we explored the impact of misspecifying the model on the performance of MB-TOST and adding a model selection step. Using the observed data, the results of both approaches were in general concordance. MB-TOST results were robust with sparse designs when the underlying PK structural model was correctly specified. Using the simulated data with a rich design, the type I error of NCA-TOST was close to the nominal level. When using the simulated model, the type I error of MB-TOST was controlled on rich and sparse designs, but using a misspecified model led to inflated type I errors. Adding a model selection step on the reference data reduced the inflation. MB-TOST appears as a robust alternative to NCA-TOST, provided that the PK model is correctly specified and the test drug has the same PK structural model as the reference drug.

Population pharmacokinetic modeling of caffeine in preterm infants with apnea of prematurity: New findings from concomitant erythromycin and AHR genetic polymorphisms

Current standard-dose caffeine therapy results in significant intersubject variability. The aims of this study were to develop and evaluate population pharmacokinetic (PPK) models of caffeine in preterm infants through comprehensive screening of covariates and then to propose model-informed precision dosing of caffeine for this population. A total of 129 caffeine concentrations from 96 premature neonates were incorporated into this study. Comprehensive medical record and genotype data of these neonates were collected for analysis. PPK modeling was performed by a nonlinear mixed effects modeling program (NONMEM). Final models based on the current weight (CW) or body surface area (BSA) were evaluated via multiple graphic and statistical methods. The model-informed dosing regimen was performed through Monte Carlo simulations. In addition to CW or BSA, postnatal age, coadministration with erythromycin (ERY), and aryl hydrocarbon receptor coding gene (AHR) variant (rs2158041) were incorporated into the final PPK models. Multiple evaluation results showed satisfactory prediction performance and stability of the CW- and BSA-based models. Monte Carlo simulations demonstrated that trough concentrations of caffeine in preterm infants would be affected by concomitant ERY therapy and rs2158041 under varying dose regimens. For the first time, ERY and rs2158041 were found to be associated with the clearance of caffeine in premature infants. Similar predictive performance and stability were obtained for both CW- and BSA-based PPK models. These findings provide novel insights into caffeine precision therapy for preterm infants.

Pharmacokinetic/Pharmacodynamic Evaluation of the Dipeptidyl Peptidase 1 Inhibitor Brensocatib for Non-cystic Fibrosis Bronchiectasis

Background and Objective

Brensocatib is an investigational, first-in-class, selective, and reversible dipeptidyl peptidase 1 inhibitor that blocks activation of neutrophil serine proteases (NSPs). The NSPs neutrophil elastase, cathepsin G, and proteinase 3 are believed to be central to the pathogenesis of several chronic inflammatory diseases, including bronchiectasis. In a phase II study, oral brensocatib 10 mg and 25 mg reduced sputum neutrophil elastase activity and prolonged the time to pulmonary exacerbation in patients with non-cystic fibrosis bronchiectasis (NCFBE). A population pharmacokinetic (PPK) model was developed to characterize brensocatib exposure, determine potential relationships between brensocatib exposure and efficacy and safety measures, and inform dose selection in clinical studies.

Methods

Pharmacokinetic (PK) data pooled from a phase I study of once-daily brensocatib (10, 25, and 40 mg) in healthy adults and a phase II study of once-daily brensocatib (10 mg and 25 mg) in adults with NCFBE were used to develop a PPK model and to evaluate potential covariate effects on brensocatib pharmacokinetics. PK–efficacy relationships for sputum neutrophil elastase below the level of quantification (BLQ) and reduction in pulmonary exacerbation and PK–safety relationships for adverse events of special interest (AESIs; periodontal disease, hyperkeratosis, and infections other than pulmonary infections) were evaluated based on model-predicted brensocatib exposure. A total of 1284 steady-state brensocatib concentrations from 225 individuals were included in the PPK data set; 241 patients with NCFBE from the phase II study were included in the pharmacodynamic (PD) population for the PK/PD analyses.

Results

The PPK model that best described the observed data consisted of two distributional compartments and linear clearance. Two significant covariates were found: age on volume of distribution and renal function on apparent oral clearance. PK–efficacy analysis revealed a threshold brensocatib exposure (area under the concentration–time curve) effect for attaining sputum neutrophil elastase BLQ and a strong relationship between sputum neutrophil elastase BLQ and reduction in pulmonary exacerbations. A PK–safety evaluation showed no noticeable trends between brensocatib exposure and the incidence of AESIs. Based on the predicted likelihood of clinical outcomes for sputum neutrophil elastase BLQ and pulmonary exacerbations, brensocatib doses of 10 mg and 25 mg once daily were selected for a phase III clinical trial in patients with NCFBE (ClinicalTrials.gov identifier: NCT04594369).

Conclusions

PPK results revealed that age and renal function have a moderate effect on brensocatib exposure. However, this finding does not warrant dose adjustments based on age or in those with mild or moderate renal impairment. The PK/PD evaluation demonstrated the clinically meaningful relationship between suppression of neutrophil elastase activity and reduction in exacerbations in brensocatib-treated patients with NCFBE, supporting further development of brensocatib for bronchiectasis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40262-022-01147-w.

Amiodarone/N-desethylamiodarone population pharmacokinetics in pediatric patients

The population pharmacokinetics (PK) of amiodarone and its active metabolite, N-desethylamiodarone (DEA) were investigated in pediatric patients with arrhythmias, mainly supraventricular tachycardias. A total of 55 patients from the department of pediatric intensive care and pediatric cardiology at Necker Enfants-malades Hospital (Paris, France) provided 72 concentrations for both amiodarone and DEA following repeated oral or intravenous administration. Blood samples drawn for biological analyses were used for drug concentrations. Plasma amiodarone concentrations were measured by a liquid chromatography method coupled with mass spectrometry (LC-MS/MS) and the data were modelled using the software Monolix 2019R2. Parent pharmacokinetics was described with a two-compartment open model and the metabolite formation was connected to the central parent compartment. Parameter estimates scaled allometrically on body weight (normalized to 70 kg) were respectively (% relative standard errors or RSE) 6.32 (31%) and 7.14 L/h (26%) for elimination (CL) and intercompartmental clearances and 167 (31%) and 3930 (32%) L for V₁ and V₂ . Oral bioavailability was 0.362 (21.5%). The clearance between subject variability (ω, square root of the variance) was 0.462 (RSE 21%). The proportional residual variabilities were respectively 0.453 (RSE 13%) and 0.423 (RSE 12%) for amiodarone and DEA respectively. The terminal half-lives were 34 and 14.5 days for amiodarone and DEA respectively. A dosage schedule was established for 3 weight bands along 2 time periods. The high pharmacokinetic variability suggests that therapeutic drug monitoring (TDM) might be useful to improve individual efficacy and safety.