IntroductionGuidelines for vancomycin therapeutic monitoring recommend using a Bayesian approach with a population pharmacokinetic model to estimate the 24 h area under the concentration‐time curve over first‐order equations. Thus, we performed an external evaluation of population pharmacokinetic models of vancomycin in neonates and compared Bayesian results with those observed in clinical practice via pharmacokinetic equations to improve therapeutic monitoring by proposing optimized initial dosing nomograms and assessing the feasibility of reduced blood sampling strategies using the most predictive models.MethodsModels were identified from the literature and evaluated via an external neonatal population. A priori predictive performance was first assessed by prediction‐based diagnostics, then by simulation‐based diagnostics and a posteriori analyses only if deemed satisfactory; model‐informed vancomycin exposure was also compared with reference first‐order pharmacokinetic equations. The best‐performing models were ultimately subjected to Monte Carlo simulations to develop new initial dosing nomograms offering the highest probability of achieving therapeutic target.ResultsA total of 28 population pharmacokinetic models were evaluated in the external dataset, which includes 72 neonates and 380 vancomycin concentrations. Eleven models had an adequate predictive performance with bias ≤ ± 15% and imprecision 30%, while the Bayesian approach yielded over 75% agreement with reference exposure values in most cases. Nonetheless, Capparelli et al. and Mehrotra et al. models performed the best overall, showing the lowest imprecisions of 16.8% and 16.9%, respectively; both models recommended higher dosage regimens than the theoretical nomogram currently applied to favor therapeutic target attainment.DiscussionWe externally evaluated numerous neonatal population pharmacokinetic models of vancomycin and used the most predictive ones to advocate new initial dosing nomograms. Clinical implementation of the Bayesian approach could reduce the time needed to reach therapeutic target and limit the number of blood samples in newborns compared with traditional pharmacokinetic equations.