Prediction of Tacrolimus Dose/Weight-Adjusted Trough Concentration in Pediatric Refractory Nephrotic Syndrome: A Machine Learning Approach

Tacrolimus pharmacokinetics in pediatric nephrotic syndrome: A combination of population pharmacokinetic modelling and machine learning approaches to improve individual prediction

Background and Aim: Tacrolimus (TAC) is a first-line immunosuppressant for the treatment of refractory nephrotic syndrome (RNS), but the pharmacokinetics of TAC varies widely among individuals, and there is still no accurate model to predict the pharmacokinetics of TAC in RNS. Therefore, this study aimed to combine population pharmacokinetic (PPK) model and machine learning algorithms to develop a simple and accurate prediction model for TAC.

Methods: 139 children with RNS from August 2013 to December 2018 were included, and blood samples of TAC trough and partial peak concentrations were collected. The blood concentration of TAC was determined by enzyme immunoassay; CYP3A5 was genotyped by polymerase chain reaction-restriction fragment length polymorphism method; MYH9, LAMB2, ACTN4 and other genotypes were determined by MALDI-TOF MS method; PPK model was established by nonlinear mixed-effects method. Based on this, six machine learning algorithms, including eXtreme Gradient Boosting (XGBoost), Random Forest (RF), Extra-Trees, Gradient Boosting Decision Tree (GBDT), Adaptive boosting (AdaBoost) and Lasso, were used to establish the machine learning model of TAC clearance.

Results: A one-compartment model of first-order absorption and elimination adequately described the pharmacokinetics of TAC. Age, co-administration of Wuzhi capsules, CYP3A5 *3/*3 genotype and CTLA4 rs4553808 genotype were significantly affecting the clearance of TAC. Among the six machine learning models, the Lasso algorithm model performed the best (R2 = 0.42).

Conclusion: For the first time, a clearance prediction model of TAC in pediatric patients with RNS was established using PPK combined with machine learning, by which the individual clearance of TAC can be predicted more accurately, and the initial dose of administration can be optimized to achieve the goal of individualized treatment.

Development and evaluation of uncertainty quantifying machine learning models to predict piperacillin plasma concentrations in critically ill patients

BACKGROUND

Beta-lactam antimicrobial concentrations are frequently suboptimal in critically ill patients. Population pharmacokinetic (PopPK) modeling is the golden standard to predict drug concentrations. However, currently available PopPK models often lack predictive accuracy, making them less suited to guide dosing regimen adaptations. Furthermore, many currently developed models for clinical applications often lack uncertainty quantification. We, therefore, aimed to develop machine learning (ML) models for the prediction of piperacillin plasma concentrations while also providing uncertainty quantification with the aim of clinical practice.

METHODS

Blood samples for piperacillin analysis were prospectively collected from critically ill patients receiving continuous infusion of piperacillin/tazobactam. Interpretable ML models for the prediction of piperacillin concentrations were designed using CatBoost and Gaussian processes. Distribution-based Uncertainty Quantification was added to the CatBoost model using a proposed Quantile Ensemble method, useable for any model optimizing a quantile function. These models are subsequently evaluated using the distribution coverage error, a proposed interpretable uncertainty quantification calibration metric. Development and internal evaluation of the ML models were performed on the Ghent University Hospital database (752 piperacillin concentrations from 282 patients). Ensuing, ML models were compared with a published PopPK model on a database from the University Medical Centre of Groningen where a different dosing regimen is used (46 piperacillin concentrations from 15 patients.).

RESULTS

The best performing model was the Catboost model with an RMSE and [Formula: see text] of 31.94-0.64 and 33.53-0.60 for internal evaluation with and without previous concentration. Furthermore, the results prove the added value of the proposed Quantile Ensemble model in providing clinically useful individualized uncertainty predictions and show the limits of homoscedastic methods like Gaussian Processes in clinical applications.

CONCLUSIONS

Our results show that ML models can consistently estimate piperacillin concentrations with acceptable and high predictive accuracy when identical dosing regimens as in the training data are used while providing highly relevant uncertainty predictions. However, generalization capabilities to other dosing schemes are limited. Notwithstanding, incorporating ML models in therapeutic drug monitoring programs seems definitely promising and the current work provides a basis for validating the model in clinical practice.