METHODS: A 'meta-model' with 4894 concentrations from 1631 neonates was built using NONMEM, and Monte Carlo simulations were performed to design an optimal intermittent infusion, aiming to reach a target AUC0-24 of 400 mg·h/L at steady-state in at least 80% of neonates.

RESULTS: A two-compartment model best fitted the data. Current weight, postmenstrual age (PMA) and serum creatinine were the significant covariates for CL. After model validation, simulations showed that a loading dose (25 mg/kg) and a maintenance dose (15 mg/kg q12h if <35 weeks PMA and 15 mg/kg q8h if ≥35 weeks PMA) achieved the AUC0-24 target earlier than a standard 'Blue Book' dosage regimen in >89% of the treated patients.

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%.

OBJECTIVE: The aim of this proof-of-concept study was to evaluate whether combining population pharmacokinetic and machine learning approaches could provide a more accurate prediction of the clearance of renally eliminated drugs in individual neonates.

METHODS: Six drugs that are primarily eliminated by the kidneys were selected (vancomycin, latamoxef, cefepime, azlocillin, ceftazidime, and amoxicillin) as 'proof of concept' compounds. Individual estimates of clearance obtained from population pharmacokinetic models were used as reference clearances, and diverse machine learning methods and nested cross-validation were adopted and evaluated against these reference clearances. The predictive performance of these combined methods was compared with the performance of two other predictive methods: a covariate-based maturation model and a postmenstrual age and body weight scaling model. Relative error was used to evaluate the different methods.

RESULTS: The extra tree regressor was selected as the best-fit machine learning method. Using the combined method, more than 95% of predictions for all six drugs had a relative error of < 50% and the mean relative error was reduced by an average of 44.3% and 71.3% compared with the other two predictive methods.

METHODS AND FINDINGS: We reviewed all GenBank submissions of HIV-1 reverse transcriptase sequences with or without protease and identified 287 studies published between March 1, 2000, and December 31, 2013, with more than 25 recently or chronically infected ARV-naïve individuals. These studies comprised 50,870 individuals from 111 countries. Each set of study sequences was analyzed for phylogenetic clustering and the presence of 93 surveillance drug-resistance mutations (SDRMs). The median overall TDR prevalence in sub-Saharan Africa (SSA), south/southeast Asia (SSEA), upper-income Asian countries, Latin America/Caribbean, Europe, and North America was 2.8%, 2.9%, 5.6%, 7.6%, 9.4%, and 11.5%, respectively. In SSA, there was a yearly 1.09-fold (95% CI: 1.05-1.14) increase in odds of TDR since national ARV scale-up attributable to an increase in non-nucleoside reverse transcriptase inhibitor (NNRTI) resistance. The odds of NNRTI-associated TDR also increased in Latin America/Caribbean (odds ratio [OR] = 1.16; 95% CI: 1.06-1.25), North America (OR = 1.19; 95% CI: 1.12-1.26), Europe (OR = 1.07; 95% CI: 1.01-1.13), and upper-income Asian countries (OR = 1.33; 95% CI: 1.12-1.55). In SSEA, there was no significant change in the odds of TDR since national ARV scale-up (OR = 0.97; 95% CI: 0.92-1.02). An analysis limited to sequences with mixtures at less than 0.5% of their nucleotide positions—a proxy for recent infection—yielded trends comparable to those obtained using the complete dataset. Four NNRTI SDRMs—K101E, K103N, Y181C, and G190A—accounted for >80% of NNRTI-associated TDR in all regions and subtypes. Sixteen nucleoside reverse transcriptase inhibitor (NRTI) SDRMs accounted for >69% of NRTI-associated TDR in all regions and subtypes. In SSA and SSEA, 89% of NNRTI SDRMs were associated with high-level resistance to nevirapine or efavirenz, whereas only 27% of NRTI SDRMs were associated with high-level resistance to zidovudine, lamivudine, tenofovir, or abacavir. Of 763 viruses with TDR in SSA and SSEA, 725 (95%) were genetically dissimilar; 38 (5%) formed 19 sequence pairs. Inherent limitations of this study are that some cohorts may not represent the broader regional population and that studies were heterogeneous with respect to duration of infection prior to sampling.