Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Pediatr., 03 December 2025

Sec. Obstetric and Pediatric Pharmacology

Volume 13 - 2025 | https://doi.org/10.3389/fped.2025.1683683

This article is part of the Research TopicPrecision Medicine in Pediatrics - Volume IIView all 25 articles

Evaluation of a Bayesian dosing calculator for vancomycin in pediatric patients with augmented renal clearance


Abdullah Alsultan,
Abdullah Alsultan1,2*Muneera Al-JelaifyMuneera Al-Jelaify3Huda AlshahraniHuda Alshahrani4Noura M. AlajmiNoura M. Alajmi5Ghadeer AlfuhaydiGhadeer Alfuhaydi6Saeed Alqahtani,Saeed Alqahtani1,2Ali SomilyAli Somily7Mashal M. AlmutairiMashal M. Almutairi8Manal Abouelkheir,

Manal Abouelkheir3,9*
  • 1Department of Clinical Pharmacy, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
  • 2Clinical Pharmacokinetics and Pharmacodynamics Unit, King Saud University Medical City, Riyadh, Saudi Arabia
  • 3Pharmacy Services, King Saud University Medical City, Riyadh, Saudi Arabia
  • 4College of Pharmacy, King Khalid University, Abha, Saudi Arabia
  • 5Department of Pharmaceutical Care Services, King Abdulaziz Medical City, Jeddah, Saudi Arabia
  • 6College of Pharmacy, Qassim University, Qassim, Saudi Arabia
  • 7Department of Pathology, College of Medicine, King Saud University and King Saud University Medical City, Riyadh, Saudi Arabia
  • 8Department of Pharmacology and Toxicology, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
  • 9Department of Clinical Pharmacy, Faculty of Pharmacy, Misr International University, Cairo, Egypt

Background: Augmented renal clearance (ARC) is increasingly recognized among pediatric oncology and intensive care patients. It can result in subtherapeutic concentrations of renally eliminated drugs like vancomycin. Bayesian dosing tools are recommended to individualize therapy, yet their performance in pediatric ARC remains underexplored. This study evaluated the predictive accuracy of the freely available Bayesian dosing calculator, NextDose, for vancomycin in pediatric patients with and without ARC. NextDose is a web-based application built on a large-population pharmacokinetic model encompassing neonates to adults with varying renal function.

Methods: A retrospective observational study included pediatric patients (1–12 years) who received vancomycin and had at least one steady-state serum concentration. ARC was defined as estimated glomerular filtration rate (eGFR) > 130 mL/min/1.73 m2. Predictive performance was assessed using relative median prediction error (rMPE, bias) and relative median absolute prediction error (rMAPE, precision; lower values indicate higher precision). Both a priori (using only clinical/demographic data) and a posteriori (drug-level informed) predictions were evaluated.

Results: A total of 112 pediatric patients were included, of whom 47 (42%) met ARC criteria. The mean age was 5.9 ± 3.4 years; 10.7% were younger than 2 years, 57.2% were aged 2–7 years, and 32.1% were 7–12 years old. a priori predictions showed high bias (rMPE 27%) and moderate precision (rMAPE 31%), with no significant differences between the ARC and non-ARC groups. In contrast, a posteriori predictions demonstrated marked improvement (rMPE −3.9%, rMAPE 13.5%), with 86% of predictions meeting the <30% prediction-error threshold. Patients with ARC exhibited superior predictive accuracy than non-ARC counterparts (rMAPE 12% vs. 17.5%, p = 0.03).

Conclusion: NextDose overestimated vancomycin concentrations in a priori predictions, suggesting it may not be suitable for initial dose calculations. Incorporating one or two measured concentrations significantly improved predictive accuracy, particularly in patients with ARC, supporting its use alongside therapeutic drug monitoring to personalize vancomycin monitoring.

Background

Augmented renal clearance (ARC) is well-recognized but often underappreciated in pediatric patients. ARC in pediatric patients is typically defined as an estimated glomerular filtration rate (eGFR) exceeding 130–160 mL/min/1.73 m2, most commonly calculated using the Schwartz formula (13). ARC is particularly prevalent among patients in pediatric intensive care units (PICUs) and pediatric oncology settings (4, 5). Studies have reported that ARC may affect up to 8%–78% in PICU patients (5). ARC results from hemodynamic alterations, including increased cardiac output and renal blood flow, as well as systemic vasodilation driven by the release of inflammatory mediators such as cytokines and prostaglandins. In addition, aggressive fluid resuscitation and the use of vasoactive agents can enhance cardiac preload, further contributing to elevated renal perfusion and augmented clearance (1). Risk factors associated with ARC in pediatric patients include low serum creatinine, febrile neutropenia, male sex, older pediatric age (e.g., adolescents), septic shock, and the use of vasoactive medications (5).

In patients with ARC, standard dosing regimens of renally eliminated drugs often result in subtherapeutic concentrations, increasing the risk of treatment failure and the emergence of antimicrobial resistance (1, 58). Despite this clinical concern, most dosing guidelines focus on patients with normal or impaired renal function, providing limited guidance for those with ARC. Patients with ARC have high eGFR and might require larger doses of renally eliminated drugs (1, 7, 9).

One such drug is vancomycin, a glycopeptide antibiotic widely used to treat serious Gram-positive infections in the pediatric population. These infections are associated with high mortality and morbidity, particularly in vulnerable populations such as those in the PICU and oncology wards. The drug has high inter- and intra-patient variability and a narrow therapeutic index. Dosing is especially challenging in patients with ARC, with multiple studies identifying it as a risk factor contributing to subtherapeutic vancomycin concentrations (8, 10). For example, a previous study has reported vancomycin trough concentrations below 5 mg/L in 35.2% of ARC cases (4). These studies highlight the importance of individualized dosing and close monitoring of vancomycin.

One approach to improve the dosing and therapeutic drug monitoring (TDM) of vancomycin is to use model-informed precision dosing (MIPD) or Bayesian dose calculators. These tools have emerged as a promising strategy to personalize antimicrobial therapy in pediatric patients. They are increasingly used to tailor dosing based on individual pharmacokinetic profiles, potentially improving target attainment (1116). Historically, vancomycin dosing was guided by trough concentrations; however, recent consensus guidelines recommend targeting an area under the concentration-time curve over the minimum inhibitory concentration (AUC24/MIC) of 400–600 to optimize efficacy and minimize nephrotoxicity (17). Bayesian forecasting facilitates estimation of AUC values using one or two measured concentrations, allowing precise, patient-specific dosing adjustments without extensive sampling. The 2020 IDSA/ASHP/SIPD/PID vancomycin dosing and monitoring guidelines specifically recommend the use of Bayesian software to achieve these AUC-based targets (17). However, the predictive performance of these calculators in pediatric patients with ARC remains largely unknown.

Given the high prevalence and clinical implications of ARC, especially in vulnerable subpopulations such as PICU and oncology patients, it is important to evaluate whether current Bayesian dosing tools can reliably predict vancomycin concentrations in these populations. This study aims to evaluate the predictive accuracy of a freely available online Bayesian dose calculator in pediatric patients with ARC and to compare its performance to that in patients without ARC.

Method

Study design

This was a single-center retrospective observational study conducted at King Saud University Medical City, Riyadh, Saudi Arabia. The inclusion criteria were pediatric patients aged 1 to 12 years who received vancomycin for a suspected or confirmed Gram-positive infection with at least one measured drug concentration available at the steady state. Patients were classified as having ARC if their eGFR exceeded 130 mL/min/1.73 m2 (5). eGFR was calculated using the bedside Schwartz equation (18). Collected data included demographic and clinical variables such as age, weight, height, sex, hospital ward, oncology status, vancomycin indication and dosing information (dose, frequency and administration time), as well as the timing and values of serum drug concentrations. Vancomycin concentrations were determined at the pharmacokinetic laboratory within King Saud University Medical City using the Architect i4000SR immunoassay analyzer (Abbott), with a standard calibration curve ranging from 3 to 100 mg/L. As per the hospital protocol, vancomycin peaks are to be collected one hour after the end of vancomycin infusion and trough levels are to be collected 30 min prior to the next dose. Serum creatinine was also measured as part of the routine clinical practice in the hospital laboratory using the Jaffe method.

Prediction accuracy of the Bayesian model

To assess predictive accuracy, we used the freely available online dose calculator NextDose (https://www.nextdose.org/) to predict vancomycin concentrations. Up to date, it is the only available open-access Bayesian dose calculator that incorporates a pediatric-specific vancomycin model. For each patient, the accuracy of the calculator was evaluated for both a priori and posteriori predictions. a priori is a term that refers to the use of a population pharmacokinetic model that is based solely on the patient's baseline demographic and clinical characteristics (e.g., age, weight, renal function) to calculate an initial dosing regimen. In contrast, a posteriori involves also integrating drug concentrations to refine individual PK estimates and perform dose adjustments accordingly. The vancomycin model implemented in NextDose is based on the population PK model described in a recent publication by Holford et al. (19). This joint population PK model for gentamicin, amikacin, and vancomycin was developed using the extensive “GAVamycin” dataset, which includes 9,901 patients ranging from premature neonates to adults across 18 studies. Key covariates incorporated in the model include: total bodyweight, postnatal age, postmenstrual age and GFR. These covariates influence both clearance and volume of distribution. Weight and age primarily affect size and maturation, whereas GFR is the principal physiological determinant of renal clearance. The model aims to utilize GFR as the key physiological parameter to distinguish the components of renal clearance for these agents, with the ultimate goal of supporting individualized dosing.

Prediction-based diagnostics and statistical analysis

We evaluated prediction accuracy using established metrics of bias and precision. Bias was calculated as the relative median prediction error (rMPE), and precision was assessed using the relative median absolute prediction error (rMAPE), where lower values indicate higher precision, using the following formulas:

rMPE(%)=1Ni=1N(Ci_PredictedCi_ObservedCi_Observed)×100
rMAPE(%)=1Ni=1N|Ci_PredictedCi_ObservedCi_Observed|×100

The criteria for acceptance were set as follows:

For both a priori and a posteriori predictions, bias (rMPE) was considered acceptable if <10% (20). While no universally accepted threshold exists for precision (rMAPE) has been established in the literature, values <30% are commonly regarded as indicative of acceptable predictive performance for a posteriori prediction (21, 22). We also calculated the percentage of observations that met the predefined precision criterion (<30% prediction error).

While the primary objective of the analysis was to evaluate the predictive accuracy across the entire cohort of pediatrics and to compare accuracy for ARC vs. non-ARC patients, we conducted additional analyses to explore potential sources of variability in predictive accuracy. The categorical variables that we tested included:

• Oncology vs. non-oncology patients;

• PICU vs. general ward admissions;

• Male vs. female

• Patients with one vs. two measured concentrations during the initial dosing occasion.

Continuous variables we tested included age, eGFR and body weight.


Continuous data are presented as mean ± standard deviations (SD), while categorical variables are presented as frequencies and percentages. Comparisons of predictive accuracy between categorical groups were performed using independent t-tests, while the influence of continuous variables on predictive accuracy was evaluated using simple linear regression. Categorical variables were compared using the Chi-square test. P-value <0.05 was considered statistically significant. All statistical analyses were conducted using the R statistical software (version 4.3.1).

Results

Baseline demographics

A total of 112 pediatric patients were included in the analysis. The mean age was 5.9 ± 3.4 years, with 10.7% of patients under 2 years, 57.2% between 2 and 7 years, and 32.1% between 7 and 12 years. The average weight was 19.8 ± 12.9 kg, and the mean height was 108.5 ± 24.4 cm. Males accounted for 55.4% of the cohort. The mean serum creatinine was 35.4 ± 17.8 µmol/L, and the average eGFR, calculated using the Schwartz bedside equation, was 130.2 ± 24.3 mL/min/1.73 m2. Oncology patients made up 15.2% of the cohort, while 42% were PICU patients. The most common indications for vancomycin therapy were sepsis or septic shock (43%), followed by pneumonia (12.5%), meningitis (11.6%), and febrile neutropenia (10.7%). The mean total daily vancomycin dose was 59.6 ± 10.5 mg/kg/day. At the first sampling occasion, a total of 185 vancomycin concentrations were measured: 73 peak levels (mean 20.4 ± 5.3 mg/L) and 112 trough levels (mean 12 ± 6.6 mg/L). On the second sampling occasion, a total of 69 observations for 53 patients were available, comprising 24 peaks and 45 trough samples.

As shown in Table 1, among the 112 patients, 47 (42%) met the ARC definition criteria. Compared to non-ARC patients, those with ARC had significantly lower serum creatinine levels (24.8 ± 7.6 vs. 43.2 ± 19.1 µmol/L, p < 0.01) and higher eGFR values (177.5 ± 48.1 vs. 96.1 ± 25.0 mL/min/1.73 m2, p < 0.01), consistent with enhanced renal function. Patients with ARC were also generally older (6.4 ± 3.7 vs. 5.9 ± 4.0 years) and had slightly higher body weight (21.2 ± 13.1 vs. 18.8 ± 12.7) than the non-ARC group, although these differences did not reach statistical significance (p > 0.05). Although oncology diagnoses were more common in the ARC group (23.4% vs. 9.2%, p = 0.07), this did not reach statistical significance. No major differences were observed in sex distribution, infection types, or PICU admission rates between groups. The overall total vancomycin daily dose was slightly higher in ARC patients; they had slightly lower average peak and trough levels, although not significant.

Table 1
www.frontiersin.org

Table 1. Baseline demographics and clinical characteristics.

Prediction accuracy of the Bayesian model

A priori predictions

For a priori predictions (n = 185), the overall bias (rMPE) for the full cohort was 27% and precision (rMAPE) was 31%, with 49% of predictions achieving the predefined precision criterion of <30% prediction error. As shown in Table 2, subgroup analyses showed comparable performance between patients with and without ARC, with bias values of 25% vs. 29% in non-ARC, p = 0.43, and precision values of 32% vs. 31% in non-ARC, p = 0.65. Similarly, 47% of predictions in the ARC group and 50% in the non-ARC group achieved precision within the <30% error threshold. Figure 1 displays observed vs. predicted concentrations for the entire cohort, as well as stratified by ARC status. Oncology patients exhibited slightly higher bias (34.5%) and lower precision (37.5%) compared to non-oncology patients (bias 26.5%, rMAPE 30.5%), though these differences were not statistically significant (p = 0.54, p = 0.38). Across clinical setting, the model's predictions in PICU patients showed greater bias (30%) and lower precision (44%) relative to those in general ward patients (bias 27%, precision 29%), with only 36% of PICU predictions falling within the <30% prediction error threshold compared to 55% in the general ward group (p = 0.3, p = 0.21). Regarding the effect of sex, predictive performance was comparable between males and females, with no significant differences observed in either bias or precision. The effects of continuous variables, age and body weight, were not statistically significant (R2 = 0.014 and 0.010 for age and weight on bias, respectively; R2 = 0.012 and 0.010 for age and weight on precision, respectively.

Table 2
www.frontiersin.org

Table 2. Bias and precision for a priori predictions.

Figure 1
Three scatter plots comparing observed and predicted concentrations. Each plot shows data points clustered around a diagonal line indicating a positive correlation. The axes range from zero to sixty for both the predicted and observed concentrations.

Figure 1. Observed vs. predicted concentration for a priori predictions. Left is the full cohort, Middle is the ARC and right is the non-ARC.

A posteriori predictions

In contrast, posteriori predictions (n = 69), which incorporate measured drug concentrations, improved due to Bayesian estimation by the model. The overall bias for the full cohort was −3.9%, and the precision was 13.5%, with 86% of predictions meeting the <30% prediction error threshold.

All subgroups showed acceptable levels of bias and precision, but there were differences in accuracy between the subgroups as shown in Table 3. ARC patients benefited significantly from having earlier measurements, with bias and precision values of −4.5% and 12%, respectively, and 92% of predictions within acceptable prediction error limits (<30% prediction error). In contrast, non-ARC patients had significantly higher bias (6.9%) and lower precision (17.5%), with only 67% of predictions achieving acceptable precision compared to ARC patients (p = 0.04 for bias, p = 0.03 for precision). Predictive performance was comparable between oncology and non-oncology patients (bias: −5.5% vs. −3%, p = 0.64; precision: 12% vs. 15%, p = 0.86). Patients in the general ward had better predictive accuracy (bias: −4.9%, precision: 9.4%) than those in the PICU (bias: 0%, precision: 17%; p = 0.045 for precision). Additionally, patients with two concentration samples on the first occasion had significantly higher precision (7.5%) compared to those with only one sample (17.5%), p = 0.02; with 97% of these predictions meeting the predefined precision criterion (<30% prediction error). Observed vs. predicted concentrations for both a priori and a posteriori predictions are presented in Figures 1, 2. These plots visually confirm improved alignment with tighter dispersion around the line of identity with a posteriori predictions, particularly among patients with ARC. Regarding the effects of continuous variables, age and body weight, none showed a significant correlation with bias or precision.

Table 3
www.frontiersin.org

Table 3. Bias and precision for posteriori predictions.

Figure 2
Three side-by-side scatter plots compare observed versus predicted concentrations. Each plot includes black dots representing data points and a diagonal line indicating perfect prediction. The plots show a moderate positive correlation, with data points scattered around the diagonal line. The axes range from zero to sixty on both scales.

Figure 2. Observed VS predicted concentration for a posterioiri predictions. Left is the full cohort, Middle is the ARC and right is the non-ARC.

Discussion

In this study, we evaluated the predictive performance of the freely accessible online Bayesian dose calculator (NextDose) for pediatric patients receiving vancomycin, with a particular focus on those with ARC. For a priori predictions (predictions based solely on demographic and clinical data), the calculator showed moderate precision but with a high level of bias, with a consistent tendency to overpredict vancomycin concentrations in all subgroups. None of the subgroups met the predefined threshold for acceptable bias, and no statistically significant differences in predictive accuracy were observed between groups at this stage. However, incorporating one or two measured concentrations markedly improved predictive performance across all subgroups, which demonstrated an acceptable level of bias (rMPE <7%) and high precision (rMAPE < 18%). These findings underscore the clinical value of Bayesian forecasting combined with therapeutic drug monitoring (TDM) to individualize dosing and suggest that, owing to the accuracy of a posteriori predictions, the need for repeated sampling may be limited.

Improved accuracy was particularly evident in specific subgroups, including patients with ARC, those managed in general wards, and those with two concentration samples obtained during the initial dosing interval. ARC patients demonstrated both lower bias and higher precision, whereas patients in the general ward and those with two samples collected during the initial dosing occasion showed significantly improved precision. The improved performance in general ward patients may be attributed to their relatively stable physiology, in contrast to the greater pharmacokinetic variability often seen in PICU patients. The improved accuracy with two samples is expected, as additional data points enhance model fitting and reduce predictive uncertainty. However, the superior accuracy observed in ARC patients is less straightforward. The model implemented in NextDose is based on a large and diverse dataset comprising 3,232 patients with 8,877 vancomycin observations across a wide age range from premature neonates to elderly adults with varying serum creatinine levels. This extensive and heterogeneous dataset may contribute to the model's robustness and its capacity to adjust for clearance variability, including the increased clearance observed in patients with ARC.

It is challenging to pinpoint the exact reason for the high bias observed in the a priori predictions, as none of the subgroups achieved acceptable levels of bias or precision. One possible explanation is that the underlying model was developed using a population with different demographic or clinical characteristics, suggesting potential interethnic variability. Additionally, reliance on retrospectively collected routine TDM data may have introduced variability that contributed to the observed error. Nonetheless, it is well-documented that a priori predictions often exhibit substantial bias, consistent with findings from previous studies. These results align with earlier research demonstrating the limitations of pharmacokinetic models when relying solely on demographic and clinical variables without incorporating drug concentration data. For instance, a study by Lv et al. evaluated the predictive performance of 10 population pharmacokinetic models in pediatric patients (23). None of the models demonstrated acceptable levels of bias or precision in the absence of drug concentration data. Predictive accuracy improved substantially when one or two observed drug concentrations were incorporated to inform subsequent predictions (a posteriori). Similarly, Frymoyer et al. demonstrated that Bayesian forecasting significantly improves target attainment for vancomycin in children, particularly when one or more drug levels are used for dose adjustment. This highlights the value of a posteriori predictions (24). On the other hand, Le J. et al. raised concerns that not all Bayesian tools are equally effective across different age groups and that their predictive performance may vary depending on the underlying population PK models they use (25). Zhao et al. found that while Bayesian models improve dose individualization, they may still underperform in neonates or critically ill children without prior model calibration using similar populations (26). In the study by Hughes et al., they evaluated the predictive performance of nine published pharmacokinetic models for vancomycin in a cohort of infants. Only two models showed an acceptable level of precision (27). Likewise, in the study by Hahn et al., they evaluated the predictive performance of two pharmacokinetic models in a cohort of pediatric patients; both models tested showed an acceptable level of bias and precision (28).

Several limitations should be acknowledged in this study. First, this was a single-center retrospective study, which may limit the generalizability of the findings to other settings. Second, the sample size for the a posteriori analyses was relatively small. Third, we only evaluated one dosing calculator; although other Bayesian dose calculators exist, they are not freely accessible and were therefore not included in this analysis. In addition, because this study used retrospective routine TDM records, factors such as acute illness severity, fluid resuscitation, or possible sampling or documentation errors were not consistently available. Prospective studies with detailed clinical documentation are warranted to better identify determinants of prediction accuracy. Finally, ARC was defined using the Schwartz equation based on serum creatinine, which may not always reflect true renal function in children. Emerging data suggest that cystatin C may provide a more accurate and dynamic marker of renal function in this population, and its use may further refine ARC diagnosis and dosing strategies (29). However, cystatin C measurement is not yet routinely implemented in clinical practice. Future studies should include prospective, multicenter studies to validate findings and improve generalizability. Studies evaluating different Bayesian dosing platforms in pediatric populations are also needed. Furthermore, integrating alternative biomarkers such as cystatin C into pharmacokinetic models could improve renal function assessment. Lastly, future investigations should assess the impact of Bayesian-guided dosing on outcomes such as target attainment, therapeutic efficacy, cost-effectiveness, and toxicity reduction.

In conclusion, the Bayesian dose calculator exhibited high bias in a priori predictions, often overestimating vancomycin concentrations, suggesting that it is not appropriate for determining initial dosing regimens. However, the incorporation of even a single measured drug concentration substantially improved the predictive accuracy across all subgroups. These findings support the use of the NextDose application with TDM to personalize vancomycin dose adjustments in pediatric patients, including those with ARC.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by The study was approved by the IRB at King Saud University Medical City (No. E-20-5333). The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants' legal guardians/next of kin because It was a non interventional retrospective study.

Author contributions

AA: Writing – original draft, Visualization, Funding acquisition, Formal analysis, Methodology, Writing – review & editing, Supervision, Conceptualization. MA-J: Visualization, Data curation, Methodology, Conceptualization, Writing – review & editing. HA: Formal analysis, Data curation, Writing – review & editing, Software. NA: Writing – review & editing, Formal analysis, Software, Writing – original draft. GA: Formal analysis, Data curation, Software, Writing – review & editing. SA: Conceptualization, Funding acquisition, Writing – review & editing, Formal analysis. AS: Funding acquisition, Data curation, Conceptualization, Writing – review & editing. MasA: Formal analysis, Data curation, Writing – review & editing. ManA: Conceptualization, Investigation, Writing – original draft, Writing – review & editing, Data curation, Supervision, Formal analysis, Methodology.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was funded by King Abdulaziz City for Science and Technology (KACST), Riyadh, Saudi Arabia, grant number: KACST 2-17-03-001-0014.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Dhont E, Van Der Heggen T, De Jaeger A, Vande Walle J, De Paepe P, De Cock PA. Augmented renal clearance in pediatric intensive care: are we undertreating our sickest patients? Pediatr Nephrol. (2020) 35(1):25–39. doi: 10.1007/s00467-018-4120-2

PubMed Abstract | Crossref Full Text | Google Scholar

2. Avedissian SN, Bradley E, Zhang D, Bradley JS, Nazer LH, Tran TM, et al. Augmented renal clearance using population-based pharmacokinetic modeling in critically ill pediatric patients. Pediatr Crit Care Med. (2017) 18(9):e388–94. doi: 10.1097/PCC.0000000000001228

PubMed Abstract | Crossref Full Text | Google Scholar

3. Hirai K, Ihara S, Kinae A, Ikegaya K, Suzuki M, Hirano K, et al. Augmented renal clearance in pediatric patients with febrile neutropenia associated with vancomycin clearance. Ther Drug Monit. (2016) 38(3):393–7. doi: 10.1097/FTD.0000000000000270

PubMed Abstract | Crossref Full Text | Google Scholar

4. Wang F, Zhou M, Wang W, Zhu Z, Yan Y. Age-associated augmented renal clearance and low BMI trigger suboptimal vancomycin trough concentrations in children with haematologic diseases: data of 1453 paediatric patients from 2017 to 2022. BMC Pediatr. (2023) 23(1):528. doi: 10.1186/s12887-023-04288-4

PubMed Abstract | Crossref Full Text | Google Scholar

5. Rhoney DH, Metzger SA, Nelson NR. Scoping review of augmented renal clearance in critically ill pediatric patients. Pharmacotherapy. (2021) 41(10):851–63. doi: 10.1002/phar.2617

PubMed Abstract | Crossref Full Text | Google Scholar

6. Alsultan A, Aldawsari MR, Alturaiq NK, Syed SA, Alsubai A, Kurdee Z, et al. Evaluation of pharmacokinetic pharmacodynamic target attainment of meropenem in pediatric patients. Pediatrics & Neonatology. (2024) 65(4):386–90. doi: 10.1016/j.pedneo.2023.09.008

PubMed Abstract | Crossref Full Text | Google Scholar

7. Gijsen M, Vlasselaers D, Spriet I, Allegaert K. Pharmacokinetics of antibiotics in pediatric intensive care: fostering variability to attain precision medicine. Antibiotics. (2021) 10(10):1182. doi: 10.3390/antibiotics10101182

PubMed Abstract | Crossref Full Text | Google Scholar

8. Scully PT, Lam WM, Coronado Munoz AJ, Modem VM. Augmented renal clearance of vancomycin in suspected sepsis: single-center, retrospective pediatric cohort. Pediatr Crit Care Med. (2022) 23(6):444–52. doi: 10.1097/PCC.0000000000002918

PubMed Abstract | Crossref Full Text | Google Scholar

9. André P, Diezi L, Dao K, Crisinel PA, Rothuizen LE, Chtioui H, et al. Ensuring sufficient trough plasma concentrations for broad-spectrum beta-lactam antibiotics in children with malignancies: beware of augmented renal clearance!. Front Pediatr. (2021) 9:768438. doi: 10.3389/fped.2021.768438

PubMed Abstract | Crossref Full Text | Google Scholar

10. Bakke V, Sporsem H, Von der Lippe E, Nordøy I, Lao Y, Nyrerød HC, et al. Vancomycin levels are frequently subtherapeutic in critically ill patients: a prospective observational study. Acta Anaesthesiol Scand. (2017) 61(6):627–35. doi: 10.1111/aas.12897

PubMed Abstract | Crossref Full Text | Google Scholar

11. Alsultan A, Alghamdi WA, Alghamdi J, Alharbi AF, Aljutayli A, Albassam A, et al. Clinical pharmacology applications in clinical drug development and clinical care: a focus on Saudi Arabia. Saudi Pharm J. (2020) 28(10):1217–27. doi: 10.1016/j.jsps.2020.08.012

PubMed Abstract | Crossref Full Text | Google Scholar

12. Alsultan A, Almofada R, Alomair S, Egelund EF, Albassam AA, Ali M, et al. Evaluation of the predictive performance of an online voriconazole dose calculator in children. Eur J Clin Pharmacol. (2024) 80:1989–93. doi: 10.1007/s00228-024-03762-x

PubMed Abstract | Crossref Full Text | Google Scholar

13. Del Valle-Moreno P, Suarez-Casillas P, Mejías-Trueba M, Ciudad-Gutiérrez P, Guisado-Gil AB, Gil-Navarro MV, et al. Model-informed precision dosing software tools for dosage regimen individualization: a scoping review. Pharmaceutics. (2023) 15(7):1859. doi: 10.3390/pharmaceutics15071859

PubMed Abstract | Crossref Full Text | Google Scholar

14. Frymoyer A, Schwenk HT, Brockmeyer JM, Bio L. Impact of model-informed precision dosing on achievement of vancomycin exposure targets in pediatric patients with cystic fibrosis. Pharmacotherapy. (2023) 43(10):1007–14. doi: 10.1002/phar.2845

PubMed Abstract | Crossref Full Text | Google Scholar

15. Darwich AS, Polasek TM, Aronson JK, Ogungbenro K, Wright DFB, Achour B, et al. Model-Informed precision dosing: background, requirements, validation, implementation, and forward trajectory of individualizing drug therapy. Annu Rev Pharmacol Toxicol. (2021) 61:225–45. doi: 10.1146/annurev-pharmtox-033020-113257

PubMed Abstract | Crossref Full Text | Google Scholar

16. Berezowska M, Hayden IS, Brandon AM, Zats A, Patel M, Barnett S, et al. Recommended approaches for integration of population pharmacokinetic modelling with precision dosing in clinical practice. Br J Clin Pharmacol. (2024) 91(4):1064–79. doi: 10.1111/bcp.16335

PubMed Abstract | Crossref Full Text | Google Scholar

17. Rybak MJ, Le J, Lodise TP, Levine DP, Bradley JS, Liu C, et al. Therapeutic monitoring of vancomycin for serious methicillin-resistant Staphylococcus aureus infections: a revised consensus guideline and review by the American society of health-system pharmacists, the infectious diseases society of America, the pediatric infectious diseases society, and the society of infectious diseases pharmacists. Clin Infect Dis. (2020) 71(6):1361–4. doi: 10.1093/cid/ciaa303

PubMed Abstract | Crossref Full Text | Google Scholar

18. Schwartz GJ, Muñoz A, Schneider MF, Mak RH, Kaskel F, Warady BA, et al. New equations to estimate GFR in children with CKD. J Am Soc Nephrol. (2009) 20(3):629–37. doi: 10.1681/ASN.2008030287

PubMed Abstract | Crossref Full Text | Google Scholar

19. Holford N, O'Hanlon CJ, Allegaert K, Anderson B, Falcão A, Simon N, et al. A physiological approach to renal clearance: from premature neonates to adults. Br J Clin Pharmacol. (2024) 90(4):1066–80. doi: 10.1111/bcp.15978

PubMed Abstract | Crossref Full Text | Google Scholar

20. Johansson AM, Ueckert S, Plan EL, Hooker AC, Karlsson MO. Evaluation of bias, precision, robustness and runtime for estimation methods in NONMEM 7. J Pharmacokinet Pharmacodyn. (2014) 41(3):223–38. doi: 10.1007/s10928-014-9359-z

PubMed Abstract | Crossref Full Text | Google Scholar

21. Hara M, Masui K, Eleveld DJ, Struys MMRF, Uchida O. Predictive performance of eleven pharmacokinetic models for propofol infusion in children for long-duration anaesthesia. Br J Anaesth. (2017) 118(3):415–23. doi: 10.1093/bja/aex007

PubMed Abstract | Crossref Full Text | Google Scholar

22. Cunio CB, Uster DW, Carland JE, Buscher H, Liu Z, Brett J, et al. Towards precision dosing of vancomycin in critically ill patients: an evaluation of the predictive performance of pharmacometric models in ICU patients. Clin Microbiol Infect. (2021) 27(5):783.e7–e14. doi: 10.1016/j.cmi.2020.07.005

Crossref Full Text | Google Scholar

23. Lv C, Lu J, Jing L, Liu TT, Chen M, Zhang R, et al. Systematic external evaluation of reported population pharmacokinetic models of vancomycin in Chinese children and adolescents. J Clin Pharm Ther. (2021) 46(3):820–31. doi: 10.1111/jcpt.13363

PubMed Abstract | Crossref Full Text | Google Scholar

24. Frymoyer A, Guglielmo BJ, Hersh AL. Desired vancomycin trough serum concentration for treating invasive methicillin-resistant staphylococcal infections. Pediatr Infect Dis J. (2013) 32(10):1077–9. doi: 10.1097/INF.0b013e318299f75c

PubMed Abstract | Crossref Full Text | Google Scholar

25. Le J, Bradley JS, Murray W, Romanowski GL, Tran TT, Nguyen N, et al. Improved vancomycin dosing in children using area under the curve exposure. Pediatr Infect Dis J. (2013) 32(4):e155–63. doi: 10.1097/INF.0b013e318286378e

PubMed Abstract | Crossref Full Text | Google Scholar

26. Zhao W, Kaguelidou F, Biran V, Zhang D, Allegaert K, Capparelli EV, et al. External evaluation of population pharmacokinetic models of vancomycin in neonates: the transferability of published models to different clinical settings. Br J Clin Pharmacol. (2013) 75(4):1068–80. doi: 10.1111/j.1365-2125.2012.04406.x

PubMed Abstract | Crossref Full Text | Google Scholar

27. Hughes JH, Tong DMH, Faldasz JD, Frymoyer A, Keizer RJ. Evaluation of neonatal and paediatric vancomycin pharmacokinetic models and the impact of maturation and serum creatinine covariates in a large multicentre data set. Clin Pharmacokinet. (2023) 62(1):67–76. doi: 10.1007/s40262-022-01185-4

PubMed Abstract | Crossref Full Text | Google Scholar

28. Hahn A, Frenck RW, Zou Y, Vinks AA. Validation of a pediatric population pharmacokinetic model for vancomycin. Ther Drug Monit. (2015) 37(3):413–6. doi: 10.1097/FTD.0000000000000153

PubMed Abstract | Crossref Full Text | Google Scholar

29. Downes KJ, Zuppa AF, Sharova A, Neely MN. Optimizing vancomycin therapy in critically ill children: a population pharmacokinetics study to inform vancomycin area under the curve estimation using novel biomarkers. Pharmaceutics. (2023) 15(5):1336. doi: 10.3390/pharmaceutics15051336

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: pediatrics, augmented renal clearance, vancomycin, Bayesian dose calculators, model-informed precision dosing

Citation: Alsultan A, Al-Jelaify M, Alshahrani H, Alajmi NM, Alfuhaydi G, Alqahtani S, Somily A, Almutairi MM and Abouelkheir M (2025) Evaluation of a Bayesian dosing calculator for vancomycin in pediatric patients with augmented renal clearance. Front. Pediatr. 13:1683683. doi: 10.3389/fped.2025.1683683

Received: 11 August 2025; Revised: 30 October 2025;
Accepted: 19 November 2025;
Published: 3 December 2025.

Edited by:

Yang Zhou, Brown University, United States

Reviewed by:

Georg Hempel, University of Münster, Germany
Jumpei Saito, National Center for Child Health and Development (NCCHD), Japan

Copyright: © 2025 Alsultan, Al-Jelaify, Alshahrani, Alajmi, Alfuhaydi, Alqahtani, Somily, Almutairi and Abouelkheir. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Abdullah Alsultan, YWJzdWx0YW5Aa3N1LmVkdS5zYQ==; Manal Abouelkheir, bWFuYWwuYWJvZWxraGVpckBtaXVlZ3lwdC5lZHUuZWc=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.