Towards Model-informed Precision Dosing of Piperacillin: Multicenter Systematic External Evaluation of Pharmacokinetic Models in Critically Ill Adults with a Focus on Bayesian Forecasting

Purpose: Inadequate piperacillin (PIP) exposure in intensive care unit (ICU) patients threatens therapeutic success. Model-informed precision dosing (MIPD) might be promising to individualize dosing; however, the transferability of published models to external populations is uncertain. This study aimed to externally evaluate the available PIP population pharmacokinetic (PopPK) models.

Methods: A multicenter dataset of 561 ICU patients (11 centers/3654 concentrations) was used for the evaluation of 24 identified models. Model performance was investigated for a priori (A) predictions, i.e., considering dosing records and patient characteristics only, and for Bayesian forecasting, i.e., additionally including the first (B1) or first and second (B2) therapeutic drug monitoring (TDM) samples per patient. Median relative prediction error (MPE) [%] and median absolute relative prediction error (MAPE) [%] were calculated to quantify accuracy and precision.

Results: The evaluation revealed a large inter-model variability (A: MPE - 135.6-78.3% and MAPE 35.7-135.6%). Integration of TDM data improved all model predictions (B1/B2 relative improvement vs. A: |MPE| 45.1/67.5%; MAPE 29/39%). The model by Kim et al. was identified to be most appropriate for the total dataset (A/B1/B2: MPE - 9.8/- 5.9/- 0.9%; MAPE 37/27.3/23.7%), Udy et al. performed best in patients receiving intermittent infusion, and Klastrup et al. best predicted patients receiving continuous infusion. Additional evaluations stratified by sex and renal replacement therapy revealed further promising models.

Conclusion: The predictive performance of published PIP models in ICU patients varied considerably, highlighting the relevance of appropriate model selection for MIPD. Our differentiated external evaluation identified specific models suitable for clinical use, especially in combination with TDM.