Model-informed Precision Dosing Research Articles

Development and application of a population pharmacokinetic model repository for caffeine dose tailoring in preterm infants

ABSTRACT Background Considerable interindividual variability for the pharmacokinetics of caffeine in preterm infants has been demonstrated, emphasizing the importance of personalized dosing. This study aimed to develop and apply a repository of currently published population pharmacokinetic (PopPK) models of caffeine in preterm infants to facilitate model-informed precision dosing (MIPD). Research design and methods Literature search was conducted using PubMed, Embase, Scopus, and Web of Science databases. Relevant publications were screened, and their quality was assessed. PopPK models were reestablished to develop the model repository. Covariate effects were evaluated and the concentration–time profiles were simulated. An online simulation and calculation tool was developed as an instance. Results Twelve PopPK models were finally included in the repository. Preterm infants’ age and body size, especially the postnatal age and current weight, were identified as the most clinically critical covariates. Simulated blood concentration–time profiles across these models were comparable. Caffeine citrate-dose regimen should be adjusted according to the age and body size of preterm infants. The developed online tool can be used to facilitate clinical decision-making. Conclusions The first developed repository of PopPK models for caffeine in preterm infants has a wide range of potential applications in the MIPD of caffeine.

Predictive performance of multi-model approaches for model-informed precision dosing of piperacillin in critically ill patients

Piperacillin (PIP)/tazobactam is a frequently prescribed antibiotic; however, over- or underdosing may contribute to toxicity, therapeutic failure, and development of antimicrobial resistance. An external evaluation of 24 published PIP-models demonstrated that model-informed precision dosing (MIPD) can enhance target attainment. Employing various candidate models, this study aimed to assess the predictive performance of different MIPD-approaches comparing (i) a single-model approach, (ii) a model selection algorithm (MSA) and (iii) a model averaging algorithm (MAA). Precision, accuracy and expected target attainment, considering either initial (B1) or initial and secondary (B2) therapeutic drug monitoring (TDM)-samples per patient, were assessed in a multicenter dataset (561 patients, 11 German centers, 3654 TDM-samples). The results demonstrated a slight superiority in predictive performance using MAA in B1, regardless of the candidate models, compared to MSA and the best single models (MAA, MSA, single models: inaccuracy ±3%, ±10%, ±8%; imprecision: <25%, <31%, <28%; expected target attainment >77%, >71%, >73%). The inclusion of a second TDM-sample notably improved precision and target attainment for all MIPD-approaches, particularly within the context of MSA and most of the single models. The expected target attainment is maximized (up to >90%) when the TDM-sample is integrated within 24 hours. In conclusion, MAA streamlines MIPD by reducing the risk of selecting an inappropriate model for specific patients. Therefore, MIPD of PIP using MAA implicates further optimization of antibiotic exposure in critically ill patients, by improving predictive performance with only one sample available for Bayesian forecasting, safety, and usability in clinical practice.

Optimization of oral isavuconazole dose for population in special physiological or pathological state: a physiologically based pharmacokinetics model-informed precision dosing.

To recommend precision dosing and improve therapeutic efficacy against invasive fungal disease, a physiologically based pharmacokinetic model (PBPK) of oral isavuconazole (ISA) was established and used to explore its disposition across populations in different physiological and pathological states. Twenty-five pharmacokinetic (PK) studies of oral ISA were identified through a systematic search of PubMed. Concentration-time data were extracted using WebPlotDigitizer. Physiochemical parameters were obtained from published literature and DrugBank. Model development and simulation used the Simcyp population-based simulator, and visual predictive check and predictive error were used for the model evaluation. Probability of target attainment and the cumulative fraction of response analyses were performed for dose optimization. The developed PBPK model was successfully validated in different populations. Most predicted concentration-time points aligned with the observed data, with acceptable predictive errors for the critical parameters. We predicted the PK profiles and parameters of ISA in a population with severe hepatic impairment (HI), a population with obesity and paediatric patients aged 1 to less than 6 years old. The probability of target attainment and cumulative fraction of response analyses indicated that the population with severe HI should have half the maintenance dose. The population with obesity and population with severe HI should have a loading dose of 300 mg every 8 h for 2 days. For paediatric patients aged 1 to less than 6 years old, a weight-based dosing regimen (5.38 mg/kg) of ISA was suggested. The predicted value aligns with observations, suggesting ISA's potential predictability in PK profiles for other populations. The recommended dosing regimens increase our understanding of the use of ISA in special populations.

The potential for therapeutic drug monitoring of belatacept and other biologicals in solid organ transplantation.

In solid organ transplantation (SOT), biologicals such as recombinant therapeutic proteins, monoclonal antibodies, fusion proteins and conjugates are increasingly used for immunosuppression, desensitization, ABO (blood group) incompatibility, antibody-mediated rejections and atypical haemolytic uremic syndrome. In this paper, we review the medical evidence available for biologicals used in SOT and the potential for improvement by the application of therapeutic drug monitoring (TDM) and model-informed precision dosing. Biologicals are used for off-label indications within the field of SOT, building on the experience from their use on labelled indications. Dosing is currently mostly standard, and experience vs. effect and toxicity is limited. Pharmacokinetic characteristics of these large, partly also immunogenic molecules differ from those of traditional small molecules. Individualization by concentration measurements and modelling has mostly been proof-of-concept or feasibility studies that lack the power to provide evidence for improvement in clinical outcome. For some drugs such as alemtuzumab, eculizumab, rituximab, tocilizumab and belatacept, studies have demonstrated significant interindividual variability in pharmacokinetics. Variability in absorption from subcutaneous administration may increase interindividual variability. There is also an economic aspect of appropriate dosing that needs to be pursued. Available assays and models to refine interpretation are in place, but trials of adequate size to document the usefulness of TDM and MIPD are scarce. Collaboration within the TDM community seems mandatory to establish studies of sufficient strength to provide evidence for the use of biologicals that are currently used off-label in SOT and furthermore to identify the settings where TDM may be beneficial.

External Evaluation of Population Pharmacokinetic Models of Piperacillin in Preterm and Term Patients from Neonatal Intensive Care.

Piperacillin/tazobactam is extensively used off-label to treat late-onset neonatal sepsis, but safety and pharmacokinetic data in this population are limited. Additionally, the organic immaturity of the newborns contributes to a high piperacillin pharmacokinetic variability. This affects the clinical efficacy of the antibiotic treatment and increases the probability of developing drug resistance. This study aimed to evaluate the predictive performance of reported piperacillin population pharmacokinetic models for their application in a model-informed precision dosing strategy in preterm and term Mexican neonatal intensive care patients. Published population pharmacokinetic models for piperacillin which included neonates in their study population were identified. From the reference models, structured models, population pharmacokinetic parameters, and interindividual and residual variability data were extracted to be replicated in pharmacokinetic software (NONMEM® version 7.4). For the clinical study, a sampling schedule was designed, and 2-3 blood samples of 250 µL were taken from neonates who met the inclusion criteria. Piperacillin plasma concentrations were determined by liquid chromatography/tandem mass spectrometry. The clinical treatment data were collected, and piperacillin plasma concentrations were estimated using reference pharmacokinetic models for an a priori or Bayesian approach. Statistical methods were used in terms of bias and precision to evaluate the differences between observed and estimated neonatal piperacillin plasma concentrations with the different approaches and to identify the pharmacokinetic model that best fits the neonatal data. A total of 70 plasma samples were collected from 25 neonatal patients, of which 15 were preterm neonates. The overall median value (range) postnatal age, gestational age, body weight, and serum creatinine at the sampling collecting day were 12 (3-26) days, 34.2 (26-41.1) weeks, 1.78 (0.08-3.90) Kg, 0.47 (0.20-0.90) mg/dL, respectively. Three population pharmacokinetic models for piperacillin in infants up to 2months were identified, and their predictive performance in neonatal data was evaluated. No pharmacokinetic model was suitable for our population using an a priori approach. The model published by Cohen-Wolkowiez et al. in 2014 with a Bayesian approach showed the best performance of the pharmacokinetic models evaluated in our neonatal data. The procedure requires two blood samples (predose and postdose), and, when applied, it predicted 66.6% of the observations with a relative median absolute predicted error of less than 30%. The population pharmacokinetic model developed by Cohen-Wolkowiez et al. in 2014 demonstrated superior performance in predicting the plasma concentration of piperacillin in preterm and term Mexican neonatal intensive care patients. The Bayesian approach, including two different piperacillin plasma concentrations, was clinically acceptable regarding bias and precision. Its application for model-informed precision dosing can be an option to optimize the piperacillin dosage in our population.

Personalized Dosing of Medicines for Children: A Primer on Pediatric Pharmacometrics for Clinicians.

The widespread use of drugs for unapproved purposes remains common in children, primarily attributable to practical, ethical, and financial constraints associated with pediatric drug research. Pharmacometrics, the scientific discipline that involves the application of mathematical models to understand and quantify drug effects, holds promise in advancing pediatric pharmacotherapy by expediting drug development, extending applications, and personalizing dosing. In this review, we delineate the principles of pharmacometrics, and explore its clinical applications and prospects. The fundamental aspect of any pharmacometric analysis lies in the selection of appropriate methods for quantifying pharmacokinetics and pharmacodynamics. Population pharmacokinetic modeling is a data-driven method ('top-down' approach) to approximate population-level pharmacokinetic parameters, while identifying factors contributing to inter-individual variability. Model-informed precision dosing is increasingly used to leverage population pharmacokinetic models and patient data, to formulate individualized dosing recommendations. Physiologically based pharmacokinetic models integrate physicochemical drug properties with biological parameters ('bottom-up approach'), and is particularly valuable in situations with limited clinical data, such as early drug development, assessing drug-drug interactions, or adapting dosing for patients with specific comorbidities. The effective implementation of these complex models hinges on strong collaboration between clinicians and pharmacometricians, given the pivotal role of data availability. Promising advancements aimed at improving data availability encompass innovative techniques such as opportunistic sampling, minimally invasive sampling approaches, microdialysis, and in vitro investigations. Additionally, ongoing research efforts to enhance measurement instruments for evaluating pharmacodynamics responses, including biomarkers and clinical scoring systems, are expected to significantly bolster our capacity to understand drug effects in children.

Systematic Evaluation of Osimertinib Population Pharmacokinetic Models in a Cohort of Dutch Adults with Non-Small Cell Lung Cancer.

Several population pharmacokinetic (popPK) studies have been reported that can guide the prediction of osimertinib plasma concentrations in individual patients. It is currently unclear which popPK model offers the best predictive performance and which popPK models are most suitable for nonadherence management and model-informed precision dosing. Therefore, the objective of this study was to externally validate all osimertinib popPK models available in the current literature. Published popPK models for osimertinib were constructed using NONMEM version 7.4.4. The predictive quality of the identified models was assessed with goodness-of-fit (GoF) plots, conditional weighted residuals (CWRES) plots and a prediction-corrected visual predictive check (pcVPC) for osimertinib and its active metabolite AZ5104. A subset from the Dutch OSIBOOST trial, where 11 patients with low osimertinib exposure were included, was used as evaluation cohort. The population GoF plots for all four models poorly followed the line of identity. For the individual GoF plots, all models performed comparable and were closely distributed among the line of identity. CWRES of the four models were skewed. The pcVPCs of all four models showed a similar trend, where all observed concentrations fell in the simulated shaded areas, but in the lower region of the simulated areas. All four popPK models can be used to individually predict osimertinib concentrations in patients with low osimertinib exposure. For population predictions, all four popPK models performed poorly in patients with low osimertinib exposure. A novel popPK model with good predictive performance should be developed for patients with low osimertinib exposure. Ideally, the cause for the relatively low osimertinib exposure in our evaluation cohort should be known. NCT03858491.

Population pharmacokinetics of TLD-1, a novel liposomal doxorubicin, in a phase I trial.

TLD-1 is a novel pegylated liposomal doxorubicin (PLD) formulation aiming to optimise the PLD efficacy-toxicity ratio. We aimed to characterise TLD-1's population pharmacokinetics using non-compartmental analysis and nonlinear mixed-effects modelling. The PK of TLD-1 was analysed by performing a non-compartmental analysis of longitudinal doxorubicin plasma concentration measurements obtained from a clinical trial in 30 patients with advanced solid tumours across a 4.5-fold dose range. Furthermore, a joint parent-metabolite PK model of doxorubicinentrapped, doxorubicinfree, and metabolite doxorubicinol was developed. Interindividual and interoccasion variability around the typical PK parameters and potential covariates to explain parts of this variability were explored. Medians standard deviations of dose-normalised doxorubicinentrapped+free Cmax and AUC0-∞ were 0.342 0.134mg/L and 40.1 18.9mg·h/L, respectively. The median half-life (95h) was 23.5h longer than the half-life of currently marketed PLD. The novel joint parent-metabolite model comprised a one-compartment model with linear release (doxorubicinentrapped), a two-compartment model with linear elimination (doxorubicinfree), and a one-compartment model with linear elimination for doxorubicinol. Body surface area on the volumes of distribution for free doxorubicin was the only significant covariate. The population PK of TLD-1, including its release and main metabolite, were successfully characterised using non-compartmental and compartmental analyses. Based on its long half-life, TLD-1 presents a promising candidate for further clinical development. The PK characteristics form the basis to investigate TLD-1 exposure-response (i.e., clinical efficacy) and exposure-toxicity relationships in the future. Once such relationships have been established, the developed population PK model can be further used in model-informed precision dosing strategies. ClinicalTrials.gov-NCT03387917-January 2, 2018.

Developing Parametric and Nonparametric Models for Model-Informed Precision Dosing: A Quality Improvement Effort in Vancomycin for Patients With Obesity.

Both parametric and nonparametric methods have been proposed to support model-informed precision dosing (MIPD). However, which approach leads to better models remains uncertain. Using open-source software, these 2 statistical approaches for model development were compared using the pharmacokinetics of vancomycin in a challenging subpopulation of class 3 obesity. Patients on vancomycin at the University of Vermont Medical Center from November 1, 2021, to February 14, 2023, were entered into the MIPD software. The inclusion criteria were body mass index (BMI) of at least 40 kg/m2 and 1 or more vancomycin levels. A parametric model was created using nlmixr2/NONMEM, and a nonparametric model was created using metrics. Then, a priori and a posteriori predictions were evaluated using the normalized root mean squared error (nRMSE) for precision and the mean percentage error (MPE) for bias. The parametric model was evaluated in a simulated MIPD context using an external validation dataset. In total, 83 patients were included in the model development, with a median age of 56.6 years (range: 24-89 years), and a median BMI of 46.3 kg/m2 (range: 40-70.3 kg/m2). Both parametric and nonparametric models were 2-compartmental, with creatinine clearance and fat-free mass as covariates to c clearance and volume parameters, respectively. The a priori MPE and nRMSE for the parametric versus nonparametric models were -6.3% versus 2.69% and 27.2% versus 30.7%, respectively. The a posteriori MPE and RMSE were 0.16% and 0.84%, and 13.8% and 13.1%. The parametric model matched or outperformed previously published models on an external validation dataset (n = 576 patients). Minimal differences were found in the model structure and predictive error between the parametric and nonparametric approaches for modeling vancomycin class 3 obesity. However, the parametric model outperformed several other models, suggesting that institution-specific models may improve pharmacokinetics management.

Model-Informed Precision Dosing Using Machine Learning for Levothyroxine in General Practice: Development, Validation and Clinical Simulation Trial.

Levothyroxine is one of the most prescribed drugs in the western world. Dosing is challenging due to high-interindividual differences in effective dosage and the narrow therapeutic window. Model-informed precision dosing (MIPD) using machine learning could assist general practitioners (GPs), but no such models exist for primary care. Furthermore, introduction of decision-support algorithms in healthcare is limited due to the substantial gap between developers and clinicians' perspectives. We report the development, validation, and a clinical simulation trial of the first MIPD application for primary care. Stable maintenance dosage of levothyroxine was the model target. The multiclass model generates predictions for individual patients, for different dosing classes. Random forest was trained and tested on a national primary care database (n = 19,004) with a final weighted AUC across dosing options of 0.71, even in subclinical hypothyroidism. TSH, fT4, weight, and age were most predictive. To assess the safety, feasibility, and clinical impact of MIPD for levothyroxine, we performed clinical simulation studies in GPs and compared MIPD to traditional prescription. Fifty-one GPs selected starting dosages for 20 primary hypothyroidism cases without and then with MIPD 2 weeks later. Overdosage and underdosage were defined as higher and lower than 12.5 μg relative to stable maintenance dosage. MIPD decreased overdosage in number (30.5 to 23.9%, P < 0.01) and magnitude (median 50 to 37.5 μg, P < 0.01) and increased optimal starting dosages (18.3 to 30.2%, P < 0.01). GPs considered lab results more often with MIPD and most would use the model frequently. This study demonstrates the clinical relevance, safety, and effectiveness of MIPD for levothyroxine in primary care.

Population pharmacokinetic analysis for dose regimen optimization of vancomycin in Southern Chinese children.

Changes in physiological factors may result in large pharmacokineticvariability of vancomycin in pediatric patients, thereby leading to either supratherapeutic or subtherapeutic exposure and potentially affecting clinical outcomes. This study set out to characterize the disposition of vancomycin, quantify the exposure target and establish an optimal dosage regimen among the Southern Chinese pediatric population. Routine therapeutic drug monitoring data of 453 patients were available. We performed a retrospective population pharmacokineticanalysis of hospitalized children prescribed intravenous vancomycin using NONMEM® software. A one-compartment PPK model of vancomycin with body weight and renal functions as covariates based on a cutoff of 2 years old children was proposed in this study. Both internal and external validation showing acceptable and robust predictive performance of the model to estimate PK parameters. The value of area under the curve over 24 h to minimum inhibitory concentration ratio (AUC0-24/MIC) ≥ 260 was a significant predictor for therapeutic efficacy. Monte Carlo simulations served as a model-informed precision dosing approach and suggested that different optimal dose regimens in various scenarios should be considered rather than flat dosing. The evaluation of vancomycin exposure-efficacy relationship indicated that lower target level of AUC0-24/MIC may be needed to achieve clinical effectiveness in children, which was used to derive the recommended dosing regimen. Further prospective studies will be needed to corroborate and elucidate these results.

Leveraging QSP Models for MIPD: A Case Study for Warfarin/INR.

Warfarin dosing remains challenging due to substantial inter-individual variability, which can lead to unsafe or ineffective therapy with standard dosing. Model-informed precision dosing (MIPD) can help individualize warfarin dosing, requiring the selection of a suitable model. For models developed from clinical data, the dependence on the study design and population raises questions about generalizability. Quantitative system pharmacology (QSP) models promise better extrapolation abilities; however, their complexity and lack of validation on clinical data raise questions about applicability in MIPD. We have previously derived a mechanistic warfarin/international normalized ratio (INR) model from a blood coagulation QSP model. In this article, we evaluated the predictive performance of the warfarin/INR model in the context of MIPD using an external dataset with INR data from patients starting warfarin treatment. We assessed the accuracy and precision of model predictions, benchmarked against an empirically based reference model. Additionally, we evaluated covariate contributions and assessed the predictive performance separately in the more challenging outpatient data. The warfarin/INR model performed comparably to the reference model across various measures despite not being calibrated with warfarin initiation data. Including CYP2C9 and/or VKORC1 genotypes as covariates improved the prediction quality of the warfarin/INR model, even after assimilating 4 days of INR data. The outpatient INR exhibited higher unexplained variability, and predictions slightly exceeded observed values, suggesting that model adjustments might be necessary when transitioning from an inpatient to an outpatient setting. Overall, this research underscores the potential of QSP-derived models for MIPD, offering a complementary approach to empirical model development.

Best Practice for Therapeutic Drug Monitoring of Infliximab: Position Statement from the International Association of Therapeutic Drug Monitoring and Clinical Toxicology.

Infliximab, an anti-tumor necrosis factor monoclonal antibody, has revolutionized the pharmacological management of immune-mediated inflammatory diseases (IMIDs). This position statement critically reviews and examines existing data on therapeutic drug monitoring (TDM) of infliximab in patients with IMIDs. It provides a practical guide on implementing TDM in current clinical practices and outlines priority areas for future research. The endorsing TDM of Biologics and Pharmacometrics Committees of the International Association of TDM and Clinical Toxicology collaborated to create this position statement. Accumulating data support the evidence for TDM of infliximab in the treatment of inflammatory bowel diseases, with limited investigation in other IMIDs. A universal approach to TDM may not fully realize the benefits of improving therapeutic outcomes. Patients at risk for increased infliximab clearance, particularly with a proactive strategy, stand to gain the most from TDM. Personalized exposure targets based on therapeutic goals, patient phenotype, and infliximab administration route are recommended. Rapid assays and home sampling strategies offer flexibility for point-of-care TDM. Ongoing studies on model-informed precision dosing in inflammatory bowel disease will help assess the additional value of precision dosing software tools. Patient education and empowerment, and electronic health record-integrated TDM solutions will facilitate routine TDM implementation. Although optimization of therapeutic effectiveness is a primary focus, the cost-reducing potential of TDM also merits consideration. Successful implementation of TDM for infliximab necessitates interdisciplinary collaboration among clinicians, hospital pharmacists, and (quantitative) clinical pharmacologists to ensure an efficient research trajectory.

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.

Prospective validation of a model-informed precision dosing tool for vancomycin treatment in neonates.

We recruited 48 neonates (50 vancomycin treatment episodes) in a prospective study to validate a model-informed precision dosing (MIPD) software. The initial vancomycin dose was based on a population pharmacokinetic model and adjusted every 36-48 h. Compared with a historical control group of 53 neonates (65 episodes), the achievement of a target trough concentration of 10-15 mg/L improved from 37% in the study to 62% in the MIPD group (P = 0.01), with no difference in side effects.

Personalized Antifungal Therapy Through Model-Informed Precision Dosing of Posaconazole.

Posaconazole is a pharmacotherapeutic pillar for prophylaxis and treatment of invasive fungal diseases. Dose individualization is of utmost importance as achieving adequate antifungal exposure is associated with improved outcome. This study aimed to select and evaluate a model-informed precision dosing strategy for posaconazole. Available population pharmacokinetic models for posaconazole administered as a solid oral tablet were extracted from the literature and evaluated using data from a previously published prospective study combined with data collected during routine clinical practice. External evaluation and selection of the most accurate and precise model was based on graphical goodness-of-fit and predictive performance. Measures for bias and imprecision included mean percentage error (MPE) and normalized relative root mean squared error (NRMSE), respectively. Subsequently, the best-performing model was evaluated for its a posteriori fit-for-purpose and its suitability in a limited sampling strategy. Seven posaconazole models were evaluated using 764 posaconazole plasma concentrations from 143 patients. Multiple models showed adequate predictive performance illustrated by acceptable goodness-of-fit and MPE and NRMSE below ±10% and ±25%, respectively. In the fit-for-purpose analysis, the selected model showed adequate a posteriori predictive performance. Bias and imprecision were lowest in the presence of two prior measurements. Additionally, this model showed to be useful in a limited sampling strategy as it adequately predicted total posaconazole exposure from one (non-)trough concentration. We validated an MIPD strategy for posaconazole for its fit-for-purpose. Thereby, this study is an important first step towards MIPD-supported posaconazole dosage optimization with the goal to improve antifungal treatment in clinical practice.

Performance of Eight Infliximab Population Pharmacokinetic Models in a Cohort of Dutch Children with Inflammatory Bowel Disease.

Efficacy of infliximab in children with inflammatory bowel disease can be enhanced when serum concentrations are measured and further dosing is adjusted to achieve and maintain a target concentration. Use of a population pharmacokinetic model may help to predict an individual's infliximab dose requirement. The aim of this study was to evaluate the predictive performance of available infliximab population pharmacokinetic models in an independent cohort of Dutch children with inflammatory bowel disease. In this retrospective study, we used data of 70 children with inflammatory bowel disease (443 infliximab concentrations) to evaluate eight models that focused on infliximab pharmacokinetic models in individuals with inflammatory bowel disease, preferably aged ≤ 18 years. Predictive performance was evaluated with prior predictions (based solely on patient-specific covariates) and posterior predictions (based on covariates and infliximab trough concentrations). Model accuracy and precision were calculated with relative bias and relative root mean square error and we determined the classification accuracy at the trough concentration target of ≥ 5 mg/L. The population pharmacokinetic model by Fasanmade was identified to be most appropriate for the total dataset (relative bias before/after therapeutic drug monitoring: -20.7%/11.2% and relative root mean square error before/after therapeutic drug monitoring: 84.1%/51.6%), although differences between models were small and several were deemed suitable for clinical use. For the Fasanmade model, sensitivity and specificity for maximum posterior predictions for the next infliximab trough concentration to be ≥ 5 mg/L were respectively 83.5% and 80% with an area under the receiver operating characteristic curve of 0.870. In our paediatric cohort, various models provided acceptable predictive performance, with the Fasanmade model deemed most suitable for clinical use. Model-informed precision dosing can therefore be expected to help to maintain infliximab trough concentrations in the target range.

Physiologically based pharmacokinetic modeling of imatinib and N-desmethyl imatinib for drug-drug interaction predictions.

The first-generation tyrosine kinase inhibitor imatinib has revolutionized the development of targeted cancer therapy and remains among the frontline treatments, for example, against chronic myeloid leukemia. As a substrate of cytochrome P450 (CYP) 2C8, CYP3A4, and various transporters, imatinib is highly susceptible to drug-drug interactions (DDIs) when co-administered with corresponding perpetrator drugs. Additionally, imatinib and its main metabolite N-desmethyl imatinib (NDMI) act as inhibitors of CYP2C8, CYP2D6, and CYP3A4 affecting their own metabolism as well as the exposure of co-medications. This work presents the development of a parent-metabolite whole-body physiologically based pharmacokinetic (PBPK) model for imatinib and NDMI used for the investigation and prediction of different DDI scenarios centered around imatinib as both a victim and perpetrator drug. Model development was performed in PK-Sim® using a total of 60 plasma concentration-time profiles of imatinib and NDMI in healthy subjects and cancer patients. Metabolism of both compounds was integrated via CYP2C8 and CYP3A4, with imatinib additionally transported via P-glycoprotein. The subsequently developed DDI network demonstrated good predictive performance. DDIs involving imatinib and NDMI were simulated with perpetrator drugs rifampicin, ketoconazole, and gemfibrozil as well as victim drugs simvastatin and metoprolol. Overall, 12/12 predicted DDI area under the curve determined between first and last plasma concentration measurements (AUClast) ratios and 12/12 predicted DDI maximum plasma concentration (Cmax) ratios were within twofold of the respective observed ratios. Potential applications of the final model include model-informed drug development or the support of model-informed precision dosing.

Population pharmacokinetics of vancomycin in term neonates with perinatal asphyxia treated with therapeutic hypothermia.

Little is known about the population pharmacokinetics (PPK) of vancomycin in neonates with perinatal asphyxia treated with therapeutic hypothermia (TH). We aimed to describe the PPK of vancomycin and propose an initial dosing regimen for the first 48 h of treatment with pharmacokinetic/pharmacodynamic target attainment. Neonates with perinatal asphyxia treated with TH were included from birth until Day 6 in a multicentre prospective cohort study. A vancomycin PPK model was constructed using nonlinear mixed-effects modelling. The model was used to evaluate published dosing guidelines with regard to pharmacokinetic/pharmacodynamic target attainment. The area under the curve/minimal inhibitory concentration ratio of 400-600 mg*h/L was used as target range. Sixteen patients received vancomycin (median gestational age: 41 [range: 38-42] weeks, postnatal age: 4.4 [2.5-5.5] days, birth weight: 3.5 [2.3-4.7] kg), and 112 vancomycin plasma concentrations were available. Most samples (79%) were collected during the rewarming and normothermic phase, as vancomycin was rarely initiated during the hypothermic phase due to its nonempirical use. An allometrically scaled 1-compartment model showed the best fit. Vancomycin clearance was 0.17 L/h, lower than literature values for term neonates of 3.5kg without perinatal asphyxia (range: 0.20-0.32 L/h). Volume of distribution was similar. Published dosing regimens led to overexposure within 24 h of treatment. A loading dose of 10mg/kg followed by 24 mg/kg/day in 4 doses resulted in target attainment. Results of this study suggest that vancomycin clearance is reduced in term neonates with perinatal asphyxia treated with TH. Lower dosing regimens should be considered followed by model-informed precision dosing.

Modeling Pharmacokinetics in Individual Patients Using Therapeutic Drug Monitoring and Artificial Population Quasi-Models: A Study with Piperacillin.

Population pharmacokinetic (pop-PK) models constructed for model-informed precision dosing often have limited utility due to the low number of patients recruited. To augment such models, an approach is presented for generating fully artificial quasi-models which can be employed to make individual estimates of pharmacokinetic parameters. Based on 72 concentrations obtained in 12 patients, one- and two-compartment pop-PK models with or without creatinine clearance as a covariate were generated for piperacillin using the nonparametric adaptive grid algorithm. Thirty quasi-models were subsequently generated for each model type, and nonparametric maximum a posteriori probability Bayesian estimates were established for each patient. A significant difference in performance was found between one- and two-compartment models. Acceptable agreement was found between predicted and observed piperacillin concentrations, and between the estimates of the random-effect pharmacokinetic variables obtained using the so-called support points of the pop-PK models or the quasi-models as priors. The mean squared errors of the predictions made using the quasi-models were similar to, or even considerably lower than those obtained when employing the pop-PK models. Conclusion: fully artificial nonparametric quasi-models can efficiently augment pop-PK models containing few support points, to make individual pharmacokinetic estimates in the clinical setting.