Optimization of Vancomycin Initial Dose in Term and Preterm Neonates by Machine Learning

Machine Learning: A Potential Therapeutic Tool to Facilitate Neonatal Therapeutic Decision Making

Bacterial infection is one of the major causes of neonatal morbidity and mortality worldwide. Finding rapid and reliable methods for early recognition and diagnosis of bacterial infections and early individualization of antibacterial drug administration are essential to eradicate these infections and prevent serious complications. However, this is often difficult to perform due to non-specific clinical presentations, low accuracy of current diagnostic methods, and limited knowledge of neonatal pharmacokinetics. Although neonatal medicine has been relatively late to embrace the benefits of machine learning (ML), there have been some initial applications of ML for the early prediction of neonatal sepsis and individualization of antibiotics. This article provides a brief introduction to ML and discusses the current state of the art in diagnosing and treating neonatal bacterial infections, gaps, potential uses of ML, and future directions to address the limitations of current studies. Neonatal bacterial infections involve a combination of physiologic development, disease expression, and treatment response outcomes. To address this complex relationship, future models could consider appropriate ML algorithms to capture time series features while integrating influences from the host, microbes, and drugs to optimize antimicrobial drug use in neonates. All models require prospective clinical trials to validate their clinical utility before clinical use.

Optimizing vancomycin dosing in pediatrics: a machine learning approach to predict trough concentrations in children under four years of age

BACKGROUND

Vancomycin trough concentration is closely associated with clinical efficacy and toxicity. Predicting vancomycin trough concentrations in pediatric patients is challenging due to significant inter-individual variability and rapid physiological changes during maturation.

AIM

This study aimed to develop a machine learning model to predict vancomycin trough concentrations and determine optimal dosing regimens for pediatric patients < 4 years of age using ML algorithms.

METHOD

A single-center retrospective observational study was conducted from January 2017 to March 2020. Pediatric patients who received intravenous vancomycin and underwent therapeutic drug monitoring were enrolled. Seven ML models [linear regression, gradient boosted decision trees, support vector machine, decision tree, random forest, Bagging, and extreme gradient boosting (XGBoost)] were developed using 31 variables. Performance metrics including R-squared (R2), mean square error (MSE), root mean square error (RMSE), and mean absolute error (MAE) were compared, and important features were ranked.

RESULTS

The study included 120 eligible trough concentration measurements from 112 patients. Of these, 84 measurements were used for training and 36 for testing. Among the seven algorithms tested, XGBoost showed the best performance, with a low prediction error and high goodness of fit (MAE = 2.55, RMSE = 4.13, MSE = 17.12, and R2 = 0.59). Blood urea nitrogen, serum creatinine, and creatinine clearance rate were identified as the most important predictors of vancomycin trough concentration.

CONCLUSION

An XGBoost ML model was developed to predict vancomycin trough concentrations and aid in drug treatment predictions as a decision-support technology.

Improving prediction of tacrolimus concentration using a combination of population pharmacokinetic modeling and machine learning in chinese renal transplant recipients

Aims

The population pharmacokinetic (PPK) model-based machine learning (ML) approach offers a novel perspective on individual concentration prediction. This study aimed to establish a PPK-based ML model for predicting tacrolimus (TAC) concentrations in Chinese renal transplant recipients.

Methods

Conventional TAC monitoring data from 127 Chinese renal transplant patients were divided into training (80%) and testing (20%) datasets. A PPK model was developed using the training group data. ML models were then established based on individual pharmacokinetic data derived from the PPK basic model. The prediction performances of the PPK-based ML model and Bayesian forecasting approach were compared using data from the test group.

Results

The final PPK model, incorporating hematocrit and CYP3A5 genotypes as covariates, was successfully established. Individual predictions of TAC using the PPK basic model, postoperative date, CYP3A5 genotype, and hematocrit showed improved rankings in ML model construction. XGBoost, based on the TAC PPK, exhibited the best prediction performance.

Conclusion

The PPK-based machine learning approach emerges as a superior option for predicting TAC concentrations in Chinese renal transplant recipients.

A machine learning approach to predict daptomycin exposure from two concentrations based on Monte Carlo simulations

Daptomycin is a concentration-dependent lipopeptide antibiotic for which exposure/effect relationships have been shown. Machine learning (ML) algorithms, developed to predict the individual exposure to drugs, have shown very good performances in comparison to maximum a posteriori Bayesian estimation (MAP-BE). The aim of this work was to predict the area under the blood concentration curve (AUC) of daptomycin from two samples and a few covariates using XGBoost ML algorithm trained on Monte Carlo simulations. Five thousand one hundred fifty patients were simulated from two literature population pharmacokinetics models. Data from the first model were split into a training set (75%) and a testing set (25%). Four ML algorithms were built to learn AUC based on daptomycin blood concentration samples at pre-dose and 1 h post-dose. The XGBoost model (best ML algorithm) with the lowest root mean square error (RMSE) in a 10-fold cross-validation experiment was evaluated in both the test set and the simulations from the second population pharmacokinetic model (validation). The ML model based on the two concentrations, the differences between these concentrations, and five other covariates (sex, weight, daptomycin dose, creatinine clearance, and body temperature) yielded very good AUC estimation in the test (relative bias/RMSE = 0.43/7.69%) and validation sets (relative bias/RMSE = 4.61/6.63%). The XGBoost ML model developed allowed accurate estimation of daptomycin AUC using C0, C1h, and a few covariates and could be used for exposure estimation and dose adjustment. This ML approach can facilitate the conduct of future therapeutic drug monitoring (TDM) studies.

Machine Learning: A New Approach for Dose Individualization

The application of machine learning (ML) has shown promising results in precision medicine due to its exceptional performance in dealing with complex multidimensional data. However, using ML for individualized dosing of medicines is still in its early stage, meriting further exploration. A systematic review of study designs and modeling details of using ML for individualized dosing of different drugs was performed. We have summarized the status of the study populations, predictive targets, and data sources for ML modeling, the selection of ML algorithms and features, and the evaluation and validation of their predictive performance. We also used the Prediction model Risk of Bias Assessment Tool (PROBAST) to assess the risk of bias of included studies. Currently, ML can be used for both a priori and a posteriori dose selection and optimization, and it can also assist the implementation of therapeutic drug monitoring. However, studies are mainly focused on drugs with narrow therapeutic windows, predominantly immunosuppressants (N = 23, 35.9%) and anti-infectives (N = 21, 32.8%), and there is currently only very limited attention for special populations, such as children (N = 22, 34.4%). Most studies showed poor methodological quality and a high risk of bias. The lack of external validation and clinical utility evaluation currently limits the further clinical implementation of ML for dose individualization. We therefore have proposed several ways to improve the clinical relevance of the studies and facilitate the translation of ML models into clinical practice.

Optimization of Ganciclovir and Valganciclovir Starting Dose in Children by Machine Learning

BACKGROUND AND OBJECTIVES

Ganciclovir (GCV) and valganciclovir (VGCV) show large interindividual pharmacokinetic variability, particularly in children. The objectives of this study were (1) to develop machine learning (ML) algorithms trained on simulated pharmacokinetics profiles obtained by Monte Carlo simulations to estimate the best ganciclovir or valganciclovir starting dose in children and (2) to compare its performances on real-world profiles to previously published equation derived from literature population pharmacokinetic (POPPK) models achieving about 20% of profiles within the target.

MATERIALS AND METHODS

The pharmacokinetic parameters of four literature POPPK models in addition to the World Health Organization (WHO) growth curve for children were used in the mrgsolve R package to simulate 10,800 pharmacokinetic profiles. ML algorithms were developed and benchmarked to predict the probability to reach the steady-state, area-under-the-curve target (AUC0-24 within 40-60 mg × h/L) based on demographic characteristics only. The best ML algorithm was then used to calculate the starting dose maximizing the target attainment. Performances were evaluated for ML and literature formula in a test set and in an external set of 32 and 31 actual patients (GCV and VGCV, respectively).

RESULTS

A combination of Xgboost, neural network, and random forest algorithms yielded the best performances and highest target attainment in the test set (36.8% for GCV and 35.3% for the VGCV). In actual patients, the best GCV ML starting dose yielded the highest target attainment rate (25.8%) and performed equally for VGCV with the Franck model formula (35.3% for both).

CONCLUSION

The ML algorithms exhibit good performances in comparison with previously validated models and should be evaluated prospectively.

Optimization of Rituximab Therapy in Adult Patients With PLA2R1-Associated Membranous Nephropathy With Artificial Intelligence

Introduction

Rituximab is a first-line treatment for membranous nephropathy. Nephrotic syndrome limits rituximab exposure due to urinary drug loss. Rituximab underdosing (serum level <2 μg/ml at month-3) is a risk factor for treatment failure. We developed a machine learning algorithm to predict the risk of underdosing based on patients' characteristics at rituximab infusion. We investigated the relationship between the predicted risk of underdosing and the cumulative dose of rituximab required to achieve remission.

Methods

Rituximab concentrations were measured at month-3 in 92 sera from adult patients with primary membranous nephropathy, split into a training (75%) and a testing set (25%). A forward-backward machine-learning procedure determined the best combination of variables to predict rituximab underdosing in the training data set, which was tested in the test set. The performances were evaluated for accuracy, sensitivity, and specificity in 10-fold cross-validation training and test sets.

Results

The best variables combination to predict rituximab underdosing included age, gender, body surface area (BSA), anti-phospholipase A2 receptor type 1 (anti-PLA2R1) antibody titer on day-0, serum albumin on day-0 and day-15, and serum creatinine on day-0 and day-15. The accuracy, sensitivity, and specificity were respectively 79.4%, 78.7%, and 81.0% (training data set), and 79.2%, 84.6% and 72.7% (testing data set). In both sets, the algorithm performed significantly better than chance (P < 0.05). Patients with an initial high probability of underdosing experienced a longer time to remission with higher rituximab cumulative doses required to achieved remission.

Conclusion

This algorithm could allow for early intensification of rituximab regimen in patients at high estimated risk of underdosing to increase the likelihood of remission.

Bridging the Worlds of Pharmacometrics and Machine Learning

Precision medicine requires individualized modeling of disease and drug dynamics, with machine learning-based computational techniques gaining increasing popularity. The complexity of either field, however, makes current pharmacological problems opaque to machine learning practitioners, and state-of-the-art machine learning methods inaccessible to pharmacometricians. To help bridge the two worlds, we provide an introduction to current problems and techniques in pharmacometrics that ranges from pharmacokinetic and pharmacodynamic modeling to pharmacometric simulations, model-informed precision dosing, and systems pharmacology, and review some of the machine learning approaches to address them. We hope this would facilitate collaboration between experts, with complementary strengths of principled pharmacometric modeling and flexibility of machine learning leading to synergistic effects in pharmacological applications.

Use of Machine Learning for Dosage Individualization of Vancomycin in Neonates

BACKGROUND AND OBJECTIVE

High variability in vancomycin exposure in neonates requires advanced individualized dosing regimens. Achieving steady-state trough concentration (C0) and steady-state area-under-curve (AUC0-24) targets is important to optimize treatment. The objective was to evaluate whether machine learning (ML) can be used to predict these treatment targets to calculate optimal individual dosing regimens under intermittent administration conditions.

METHODS

C0 were retrieved from a large neonatal vancomycin dataset. Individual estimates of AUC0-24 were obtained from Bayesian post hoc estimation. Various ML algorithms were used for model building to C0 and AUC0-24. An external dataset was used for predictive performance evaluation.

RESULTS

Before starting treatment, C0 can be predicted a priori using the Catboost-based C0-ML model combined with dosing regimen and nine covariates. External validation results showed a 42.5% improvement in prediction accuracy by using the ML model compared with the population pharmacokinetic model. The virtual trial showed that using the ML optimized dose; 80.3% of the virtual neonates achieved the pharmacodynamic target (C0 in the range of 10-20 mg/L), much higher than the international standard dose (37.7-61.5%). Once therapeutic drug monitoring (TDM) measurements (C0) in patients have been obtained, AUC0-24 can be further predicted using the Catboost-based AUC-ML model combined with C0 and nine covariates. External validation results showed that the AUC-ML model can achieve an prediction accuracy of 80.3%.

CONCLUSION

C0-based and AUC0-24-based ML models were developed accurately and precisely. These can be used for individual dose recommendations of vancomycin in neonates before treatment and dose revision after the first TDM result is obtained, respectively.

Predicting busulfan exposure in patients undergoing hematopoietic stem cell transplantation using machine learning techniques

PURPOSE

This study aimed to establish an optimal model to predict the busulfan (BU) area under the curve at steady state (AUCss) by using machine learning (ML).

PATIENTS AND METHODS

Seventy-nine adult patients (age ≥18 years) who received BU intravenously and underwent therapeutic drug monitoring from 2013 to 2021 at Fujian Medical University Union Hospital were enrolled in this retrospective study. The whole dataset was divided into a training group and test group at the ratio of 8:2. BU AUCss were considered as the target variable. Nine different ML algorithms and one population pharmacokinetic (pop PK) model were developed and validated, and their predictive performance was compared.

RESULTS

All ML models were superior to the pop PK model (R2 = 0.751, MSE = 0.722, 14 and RMSE = 0.830) in model fitting and had better predictive accuracy. The ML model of BU AUCss established through support vector regression (SVR) and gradient boosted regression trees (GBRT) had the best predictive ability (R2 = 0.953 and 0.953, MSE = 0.323 and 0.326, and RMSE = 0.423 and 0.425).

CONCLUSION

All the ML models can potentially be used to estimate BU AUCss with the aim of facilitating rational use of BU on the individualized level, especially models built by SVR and GBRT algorithms.

A Hybrid Algorithm Combining Population Pharmacokinetic and Machine Learning for Isavuconazole Exposure Prediction

OBJECTIVES

Maximum a posteriori Bayesian estimation (MAP-BE) based on a limited sampling strategy and a population pharmacokinetic (POPPK) model is used to estimate individual pharmacokinetic parameters. Recently, we proposed a methodology that combined population pharmacokinetic and machine learning (ML) to decrease the bias and imprecision in individual iohexol clearance prediction. The aim of this study was to confirm the previous results by developing a hybrid algorithm combining POPPK, MAP-BE and ML that accurately predicts isavuconazole clearance.

METHODS

A total of 1727 isavuconazole rich PK profiles were simulated using a POPPK model from the literature, and MAP-BE was used to estimate the clearance based on: (i) the full PK profiles (refCL); and (ii) C24h only (C24h-CL). Xgboost was trained to correct the error between refCL and C24h-CL in the training dataset (75%). C24h-CL as well as ML-corrected C24h-CL were evaluated in a testing dataset (25%) and then in a set of PK profiles simulated using another published POPPK model.

RESULTS

A strong decrease in mean predictive error (MPE%), imprecision (RMSE%) and the number of profiles outside ± 20% MPE% (n-out20%) was observed with the hybrid algorithm (decreased in MPE% by 95.8% and 85.6%; RMSE% by 69.5% and 69.0%; n-out20% by 97.4% and 100% in the training and testing sets, respectively. In the external validation set, the hybrid algorithm decreased MPE% by 96%, RMSE% by 68% and n-out20% by 100%.

CONCLUSION

The hybrid model proposed significantly improved isavuconazole AUC estimation over MAP-BE based on the sole C24h and may improve dose adjustment.

A HPLC-DAD method to facilitate large-scale therapeutic drug monitoring of dalbavancin

Dalbavancin, a long-acting lipoglycopeptide antibiotic targeting susceptible Gram-positive bacteria, is WHO critically important antibiotic, increasingly used in critical situations such as osteoarticular infections. To ensure its effectiveness and its safety, the therapeutic drug monitoring (TDM) of dalbavancin is strongly recommended. In the absence of an available minimum inhibitory concentration (MIC), the European Committee on Antimicrobial Susceptibility Testing (EUCAST) recommends a breakpoint of 0.125 mg/L for Staphylococcus aureus, corresponding to a trough target concentration of 25 mg/L. Nowadays, the TDM is usually performed using a high-performance liquid chromatography (HPLC) method coupled with a tandem mass spectrometry. However, this expensive and specialized equipment and reagents may be difficult to acquire for non-specialized laboratories. The use of HPLC coupled with diode array detector (DAD) facilitates TDM with a lower cost, while preserving the reliability of the results. Our aim was to provide a sensitive and specific method, relying on HPLC-DAD for extending the TDM of dalbavancin beyond non-specialized labs, therefore maximizing its efficiency and Benefit/risk ratio. Our method complied with the European Medicines Agency guidelines of bioanalytical validation. Irrespective of the concentrations of dalbavancin, the coefficient of variation < 10% confirmed the reliability of this analytical method, with a calibration curve ranging from 5 to 100 mg/L. No interferences nor carryover was observed. Our HPLC-DAD method, combined with a simple extraction, provides a widely usable, inexpensive and easy-to-implement new asset for the TDM of Dalbavancin.