Skip to main content

Establishment of a population pharmacokinetics model of vancomycin in 94 infants with septicemia and its application in individualized therapy



We aim to develop a population pharmacokinetics (PopPK) model of vancomycin for the treatment of septicemia in infants younger than one year. Factors influence of the PK was investigated to optimize vancomycin dosing regimen.


The nonlinear mixed effects modelling software (NONMEM) was used to develop the PopPK model of vancomycin. The stability and predictive ability of the final model were assessed by using normalized prediction distribution errors (NPDE) and bootstrap methods. The final model was subjected to Monte Carlo simulation in order to determine the optimal dose.


A total of 205 trough and peak concentrations in 94 infants (0–1 year of age) with septicemia were analyzed. The interindividual variability of the PK parameter was described by the exponential model. Residual error was better described by the proportional model than the mixed proportional and addition models. Serum creatinine concentration and body weight are the major factors that affect the PK parameters of vancomycin. The clearance was shown to be higher when ceftriaxone was co-treated. More than two model evaluation methods showed better stability than the base model, with superior predictive performance, which can develop individualized dosing regimens for clinical reference. Through prediction of final model, the trough concentration was more likely < 5 mg/L when a routine dose of 10 mg/kg is administered every 6 h to 3–9-month-old infants. Therefore, the dose should be increased in the treatment of infant septicemia.


The stable and effective PopPK model of vancomycin in Chinese infants with septicemia was established. This model has satisfactory predictive ability for clinically individualized dosing regimens in this vulnerable population.

Peer Review reports


Over the last decades, the most frequent causative agents of septicemia in neonates have been reported to be the Gram-positive cocci infections, especially methicillin-resistant Staphylococcus aureus (MRSA) [1]. Vancomycin is widely used in clinical scenarios as the antimicrobial therapy for targeted or empiric treatment of neonatal sepsis [2]. It is the first glycopeptide antibiotic with triple mechanisms, inclusive of inhibiting the synthesis of cell walls and RNA in the cytoplasm, and altering membrane permeability in bacteria [3]. However, there remain numerous disputes regarding the current individualization of vancomycin dose in clinical practice. Even with the routine dose of vancomycin in newborns admitted in the intensive care unit, we are still unable to easily obtain trough and peak plasma concentrations of 5 to 10 mg/liter and 20 to 50 mg/liter in neonate patients [4].

Vancomycin-associated nephrotoxicity is a critical appraisal of risk with high-dose therapy for infants and adults [5]. Nevertheless, the immaturity of physiological processes and organ functions, such as liver’s and kidney’s, predisposes infants to disparate drug disposition and responses compared with that in adults [5, 6]. It is especially important to avoid exposure to unnecessarily high peak or trough concentrations of vancomycin in the premature newborn. The pharmacokinetic (PK) and pharmacodynamic (PD) parameters constantly change across the pediatric groups. The disposition and response of infants are not only different from that of other aged children but also widely vary within themselves [6, 7]. These issues along with the lack of pharmacokinetics data on infants lead to remarkable difficulty in establishing the regimen for infants. Based on the assessment of population PK models, population PK parameter estimates are less biased than those obtained using the naive and standard two-stage approaches. Thus, the NONMEM software is the most widely used software for the characterization of population PK in clinical practice in order to solve the problem [8] . To date, vancomycin has been one of the most studied antibiotics using the PopPK in neonates and adults [9]. Approaches to vancomycin clinical dosing were developed on basis of the PK parameters of individuals [10]. The method does not only conduct statistical analysis using the sparse uneven data for each patient with a few samplings, but also integrate gestational age, gestational weeks, weight, and other important data of infants for PK data processing, which suitable for infants [10, 11]. Our previous investigation characterized the PK parameters of neonates [12], but the optimal use of vancomycin in infants need to be conducted continually. It is challenging to develop and validate a model that prospectively predicts individual exposure in infants and further optimize dosing regimens. Vancomycin is always used to combine with other antibiotics such as ceftriaxone and meropenem. In order to know whether drug interactions influence the blood concentration of vancomycin, this study aims to identify individual factors that affect drug interactions of ceftriaxone that subsequently influences variability and to establish dosing regimens for infants in China.

In this study, the PK parameters that affect the PK of vancomycin were completely elucidated, and the PopPK of vancomycin in infants (0–1 year of age) with septicemia in a Chinese population was developed. Based on the PopPK model, we screened for factors that might significantly influence the PK of vancomycin and provide a reference for the individualization of treatment with vancomycin in clinical practice. Our work would potentially reduce the incidence of adverse reactions originated from the use of vancomycin and improve efficacy and safety.


Patients and data collection

A total of 94 Chinese infants who received vancomycin for the treatment of septicemia at a grade A special hospital between January 2009 and December 2015 were enrolled. The following data were collected: blood concentration of vancomycin, daily vancomycin dose, demographic characteristics of patients, blood examination results, hepatic and renal function data, and co-administered medicines. Serum vancomycin concentration was determined by chemiluminescence analysis. Demographic and clinical data of infants are summarized and presented in Table 1.

Table 1 Demographic and Clinical Data of Infants

Using the fluorescence polarization immunoassay method with the AxSYM system, serum vancomycin concentrations were measured as part of the therapeutic drug monitoring activity in the hospital. The assay sensitivity limit was 2.00 mg/L. Coefficients of variation were 4.26, 2.94, and 4.06% at concentrations of 7.0, 35.0, and 75.0 mg/L, respectively (package insert, Axsym system, Abbot Laboratories, AbbotPark, IL, USA).

PopPK modeling

The PopPK model was developed using the NONMEM version 7.2 (Icon Inc., PA, the USA), Fortran compiler: Intel Fortran XE 2011 Update 13 (Intel Corp, CA, the USA),Wings for NONMEM (version 7.2, Nick Holford, University of Auckland, Auckland, New Zealand,, Xpose (version 4.2.1, Department of Pharmaceutical Biosciences at Uppsala University), R package (version 2.15.1,, Perl-speaks-NONMEM version 3.5.3 (, Perl (ActivePerl-5.14.2, 64 bit), and Pirana (version 2.7.0, NONMEM was used to employed pharmacokinetic analysis. Based on individual dataset, population values of the PK parameters were considered as fixed-effect parameters. Inter- and intra-individual variabilities were estimated as random effects.

205 trough and peak about vancomycin concentrations versus time were fitted to both one- and two-compartmental models with first-order elimination. The most suitable compartmental model was determined to be specified to NONMEM by ADVAN1-TRANS2 or ADVAN3-TRANS4 subroutines. By using these model specifications, the fixed-effect PK parameters were directly estimated. Total body clearance (CL) and volume of distribution (V) are for the one-compartmental model, while CL, volume of distribution of the central compartment (V1), intercompartment clearance (Q), and volume of distribution of the peripheral compartment (V2) for the two-compartmental model. Additive, proportional, and exponential error models were tested to describe inter- and intraindividual variabilities.

In preliminary screening phase, each of covariates that significantly improved the predictive ability of the basic model would be included. Observed concentration–time profile was used to compare. The influence of covariates on the PK parameters of vancomycin as follows: weight, age, sex, serum creatinine (SCR) concentration, blood urea nitrogen, creatinine clearance rate [CLCR, using the Modification of Diet in Renal Disease/MDRD4 equation [10] and Cockcroft–Gault (C–G) equation [13]], serum albumin concentration, aspartate transaminase and alanine transaminase levels, and concomitant drugs (ceftriaxone, meropenem, gentamicin, furosemide, ibuprofen, and dexamethasone). The significant covariates were then cumulatively added (forward stepwise fashion) to the model in the order of their contribution in reducing the objective function value (OFV, − 2 log likelihood difference) in the preliminary analysis until no further reduction in OFV is observed. Finally, back elimination was conducted to eliminate any unnecessary covariate from the full regression model in the descending order of their contribution to the change in OFV.

The best compartmental model (one vs. two compartments), error model (additive vs. proportional vs. combined error models), and the retention of covariate (s) in the model were determined with the statistical significance of the model and were evaluated via the likelihood ratio test using the minimum value of the OFV, as produced by the NONMEM program. Changes in OFV of > 6.63 and > 9.21 were considered significant based on the χ2 distribution with degree of freedom (df) = 1 and 2 (both P < 0.01), respectively [14].

Other diagnostic criteria were a reduction in the unexplained interindividual variability for the associated PK parameters and an improvement in the graphic diagnostic model. Graphics were obtained using plotted observed vs. predictive concentrations and observed vs. weighted residuals (predicted minus observed concentrations and weighted by standard deviation). First-order conditional estimation interaction was utilized to estimate the PopPK parameters. Covariates were screened according to the stepwise method, and the final PopPK model of vancomycin was then established by backward elimination. The exponential model was used to describe interindividual variability of the PK parameters. Residual error was better described by the exponential model than the mixed exponential and addition models.

Covariate screening

Before covariate analyses, any correlation between the covariates (correlation coefficient > 0.5) should be examined, and the study should only choose one for the analysis. Using the Xpose 4 program package of the R language to draw each co-distribution map, it was examined whether the distribution of data is normal [15]. Parameter values of the final model is shown in Table 2 and screening process was summarized in Table 3.

Table 2 Parameter Value of the Final Model
Table 3 Model Selection Processa

Model evaluation

Internal evaluation is the method in which original modeling data are used in resampling technology and in the establishment of a validation dataset to validate the model. The normalized prediction distribution errors (NPDE), and bootstrap methods were used to test the model [16]. The results of boostrap in final model were summarized in Table 4. The stability and predictive ability of the final model were evaluated using NPDE, bootstrap method, and external data. The results of the model average prediction error (MPE) and average absolute error (MAE) calculated by the external evaluation method are shown in Table 5.

Table 4 Bootstrap Results of Final Model
Table 5 Prediction Ability of the Basic and Final Models

Monte Carlo simulation

The final model was subjected to Monte Carlo simulation to determine the optimal dose under which the drug concentration of a certain patient can be controlled within the reference concentration. According to the range of the commonly used dose in infants and factors that set the scene, data files of each scene were written. Based on the PopPK model of vancomycin, data files of each scenario were simulated using the Monte Carlo method, and data of 5000 were generated, which were completed by the ONLYSIMULATION and SUBPROBLEMS modules in the NONMEM software SIMULATION.


Basic model

This study established the PopPK model according to results of the conventional therapeutic drug monitoring (TDM), which only had the fixed sampling point trough and peak concentrations. Data on vancomycin concentrations versus time were fitted to both one- and two-compartmental models, respectively. Because sparse sampling method, there is not available data enough to establish a two-compartmental (Table 3). The results preferred the one-compartmental model as more suitable model to describe this set of data.

The additive, proportional, and exponential error models were tested to describe inter- and intra-individual variabilities. Interindividual variability was estimated using an exponential error model. Because of the extremely large relative standard deviation  in the interindividual variation  of V, which was not allowed, the ETA test was not conducted. Bootstrap is used to estimate statistical variance and interval estimation of the statistical method, which could be used to estimate the confidence interval of parameter values of 2.5–97.5%. If parameter values of the 95% confidence interval include 0, then the parameter estimates are not reliable.

For the residual variance model, the mixed model of proportional and additive errors was examined, and the proportional residual model is adopted.

Covariate Models

As depicted in Fig. 1, data obtained in this study had significant differences between the newborns and infants. The distribution of some covariates was not normal, including body weight, age, and SCR. Because the relationship between the clearance rate and body weight was obviously different at different phages, the method of subsection and different velocities should be adopted. Because the OFV value after adding weight generally decreased by > 100, body weight was considered an important covariate that influences the clearance rate and should be retained. Results of the second method showed that the correction of the prediction error was > 50%, which did not meet the requirements, and results of the first method had good stability. Thus, model 2 should be chosen for the next step of modeling.

Fig. 1

Correlation Analysis of Covariates in Infants

First, based on the first method model, SCR concentration and CLCR were added to the CL as an impact factor; then, the two covariates were compared to identify the one that could better describe the CL. Results showed that OFV values decrease after the addition of SCR concentration and CLCR, and the OFV of the multiplication method would decrease more than the addition method. Because the OFV of the two models were similar and creatinine value was more convenient to obtain, which is much useful in future clinical development, the model for SCR as the covariate should be chosen.

The effects of concomitant drugs on clearance were investigated. The OFV of meropenem and ceftriaxone covariates decreased by > 3.84. Thus, meropenem and ceftriaxone should both be incorporated. The OFV of meropenem elimination decreased < 6.635 (P > 0.01), which was not statistically significant. Meanwhile, combination with meropenem neither had an effect on CL (< 20%) nor had clinical significance. Thus, meropenem should be removed and other covariates should be retained.

Final model

The final regression model was as follows:

$$ CL\left(\frac{L}{hr}\right)=10.3\times {\left(\frac{\mathrm{WT}\left(\mathrm{kg}\right)}{70}\right)}^{1.06}\times {\left(\frac{\mathrm{SCR}\left(\mu \mathrm{mol}/\mathrm{L}\right)}{20}\right)}^{-0.315}\times {1.46}^{DC} $$
$$ V(L)=50.6\times \left(\frac{\mathrm{WT}\left(\mathrm{kg}\right)}{70}\right) $$

DC = 1 when incorporated with ceftriaxone, then DC = 0.

SCR and body weight were the major factors influencing the PK parameters of vancomycin. The clearance would be larger when ceftriaxone is incorporated.

Compared with the basic model, the external validation of the final model showed that the fitting degree of the predicted and measured values obviously improved, and the individual value average prediction errors of the final model was smaller but higher in the forecast precision. Then, the final model showed better fitting effect after adding the fixed-effect factors (Table 2).

Model evaluation

In order to calculate the predictive error of each sampling point, we made MPE and MAE diagram for the basic model and the final model. It can be seen from Table 5 that the inter-individual variation of the pharmacokinetic parameters in the final vancomycin PopPK model is significantly reduced compared with the basic model. The NPDE was used to validate the prediction performance of the model. The NPDE was homogeneity of variance and normal distribution according to the quantile plots (Fig. 2), NPDE distribution, and statistical test results. The t-test result was 0.566, and the Fisher’s variance test result was 0.859. Such results indicated that the NPDE-preferred model can be used to generate analog data.

Fig. 2

Normalized prediction distribution error (NPDE) for the final model. a Quantile-quantile plots of NPDE vs. the expected standard normal distribution; b Histogram of NPDE values with the standard normal distribution overlay; c Scatter plot of the time vs. NPDE; d. Scatterplot of predictions vs. NPDE

At the same time, 1000 data were created for bootstrap validation based on a non-parametric method to obtain 1000 set of model parameters for statistical calculations. If the parameter value was not significantly different, then the model was stable. In this study, more than two model validation methods showed that the final internal model had good stability with superior predictive performance, which can be used to develop individualized dosing regimens for clinical reference.

The performance of the final covariate model was evaluated by visual inspection of diagnostic scatter plots. The robustness of the model was assessed using a nonparametric bootstrap, with replacement, of 1000 NONMEM runs of the final model; the bootstrap median parameter values and the percentile bootstrap 95% intervals were compared with the respective values estimated from the final model. The final model’s 25 estimates were terminated. In addition, the 2.5 to 97.5% confidence interval of CL, V and the estimated results of the covariable parameters were not covered to 0, which met the requirements. All of results are summarized in Table 4.

Monte Carlo simulation

The steady state of infants administered with vancomycin was designed, and an analysis was conducted using Monte Carlo simulation to generate simulated data. In children with infection who were taking vancomycin at a steady state, the Monte Carlo simulation generated analog data analysis results showing that trough concentrations of the dosing regimen (15 mg/kg) every 24 h in infants who are 29 weeks premature would be < 5 mg/L, and this did not even meet the dose requirement of 10–15 mg/L [4]. The dosing interval should be shortened; although trough concentrations in 28-day-old neonates were only 4.3 mg/L, the peak concentrations were up to 72 mg/L. Trough concentrations are more likely to be < 5 mg/L when a dose of 10 mg/kg is administered every 6 h in 3–9-month-old infants; thus, the dose should be increased. The clinical application mode of vancomycin has been established and used in individual drug administration in clinical practice, as shown in Fig. 3.

Fig. 3

Distribution of the concentration range of vancomycin. Vancomycin concentration–time profile from the Monte Carlo simulation with five typical cases. Deep and light pink colors represent the 5th and 95th percentiles of the simulated data, respectively. Red solid line represents the lower and upper boundaries of the therapeutic range of vancomycin. Horizontal axis: Time after achieving steady state (h); vertical axis: vancomycin concentration (mg/L). a. Neonates with 0.95 kg body weight (gestational age: 29 weeks, SCR: 100umol/L, 19 mg every 24 h); b. Neonates with 2.4 kg body weight (gestational age: 39 weeks, SCR: 70 μmol/L, 36 mg every 8 h: 0:00, 8:00, and 16:00); c. Neonates with 4 kg body weight (gestational age: 39 weeks, 28 days, SCR: 60 μmol/L, 60 mg every 8 h: 0:00,8:00, and 16:00); d. Infants with 5 kg body weight (gestational age: 39 weeks, 3 months, SCR: 32 μmol/L, 50 mg every 6 h: 0:00,6:00,12:00, and 18:00); e. Infants with 8 kg body weight (gestational age: 39 weeks, 9 months, SCR: 28 μmol/L, 80 mg every 6 h: 0:00, 6:00, 12:00, and 18:00)


The PK of vancomycin in compartmental model studies, including one-, two-, and three-compartmental models, is controversial. In recent years, one-compartmental model has been successfully used in the study of neonates [17], children [18], adults [19], and different disease statuses of patients [20], which can better meet the needs of the fitting model and estimation of the PK parameters. Because of limited information on drugs, the results preferred two-compartment model as better model to describe the present population. In this study, the PopPK model was successfully established using the one-compartment model as the basis of the structural model for PopPK analysis.

The PopPK model can fully consider the intra- and inter-individual variability; thus, it is advantageous for individual drug delivery. The mixture model was fitted with the simultaneous presence of proportional and additive errors. Factors that influence the PK of vancomycin were selected and investigated, and the effects of weight, SCR, and combined use of drugs were investigated. When screening weight, two methods were used, and meropenem and ceftriaxone as concomitant drugs were also assessed.

Most of the vancomycin in the body were eliminated by the prototype, and the effect of renal function [CLCR or SCR concentrations] on vancomycin clearance rate was investigated. The CLCR data were estimated using the Schwartz formula, and the final model indicated that SCR was an important factor affecting the vancomycin clearance rate. Moreover, the results were consistent with the those reported by Reveilla [21] and Llopis-Salvi [22]. This study first used LZL12B1 as the basic model, adding SCR concentration and CLCR as factors influencing CL, and compared the two covariates to determine which among the two can better describe the CL. Results showed that the addition of SCR concentration and CLCF after the model OFV decreased, and there was a greater reduction in multiplication and addition after the addition of OFV.

Age, weight, and height had an extremely high linear correlation, and height and TP and TP and ALB also had correlation coefficient of > 0.5; thus, attention should be paid to the association between these factors, and such factors must not be added at the same time. In addition, after adding weight, the OFV value decreased by > 100; thus, weight should be retained because it has an important influence on the clearance rate of covariates. Stability and accuracy were in accordance with the requirements.

This study examined the effect on vancomycin elimination rate while administered in combination of other antimicrobial agents (ceftriaxone, meropenem, and gentamicin), furosemide, ibuprofen, and dexamethasone. Based on the basic model, factors of drug combination were investigated, and those that accounted for > 10% of cases were internalized. The effect of drug combination on the clearance rate was investigated using LZL12WS3, respectively in additional method with meropenem and ceftriaxone covariates. Finally, after drawing the Xpose program [15], the distribution of covariates showed a straight line, indicating that the difference between values of covariates and the final model was not statistically different. Thus, the final model only included ceftriaxone.

Because the renal pathway is primarily involved in vancomycin elimination, it is not surprising that the creatinine level is a significant factor of vancomycin clearance. More that 80% of vancomycin intravenously administered are excreted via the urine, and the total body clearance rate is related to renal function. Thus, ceftriaxone may competitively inhibit vancomycin. This study showed that the use of ceftriaxone increases vancomycin clearance rate. Therefore, the protein binding rate of ceftriaxone may be high and have competitive inhibition to vancomycin, and the uncombined vancomycin levels increasingly cleared up [23]. Besides, ceftriaxone may affect the glomerular filtration of vancomycin because the organs of infants were not mature. Similar reports have not yet been published, and the detailed mechanism must be explored in future studies.

Based on the recent guidelines, vancomycin trough concentration of > 10 mg/L was recommended to prevent the development of resistance [24]. Vancomycin trough concentrations of 15–20 mg/L were used to treat serious infections such as endocarditis, osteomyelitis, meningitis, and hospital acquired pneumonia. Predictive performance in the validation step was evaluated by applying the final model to the validation group. The results of this study can be directly applied in clinical practice, and the model can be used to obtain specific PK parameters of a patient to establish vancomycin dosage regimen in patients similar to those in the present study.

The PK parameters of vancomycin in infants with septicemia younger than one year in China were established, and factors influencing the variability of these PK parameters were identified. Inter- and intra-individual variabilities of the PK parameters were determined. The established PK model in this study was stable and had good prediction ability, which can be used to obtain specific PK parameters to establish vancomycin dosage regimen in patients similar to those in the present study and to develop individualized clinical dosing.


The commended vancomycin dose for infants according to the latest guidelines and instructions is more likely not to obtain target trough concentrations [4, 25]. The median concentration during 29 weeks of gestation was only 3.66 mg/L. We found the dose should be increased and the dosing interval should be decreased by applying final model. As a result, the clinical application of the PopPK model of vancomycin established in the present study can promote the rational use of drugs in clinical practice.

Availability of data and materials

All data generated or analysed during this study are included in this published article.



Population pharmacokinetics


Methicillin-resistant Staphylococcus aureus






Serum creatinine


Objective function value


Normalized prediction distribution errors


Therapeutic drug monitoring


Aspartate aminotransferase


Alanine aminotransferase


Blood Urea Nitrogen


Total protein


  1. 1.

    Kriebs JM. Staphylococcus infections in pregnancy: maternal and neonatal risks. J Perinat Neonatal Nurs. 2016;30(2):115–23.

    Article  PubMed  Google Scholar 

  2. 2.

    Peixoto PB, Massinhani FH, Netto Dos Santos KR, Chamon RC, Silva RB, Lopes Correa FE, et al. Methicillin-resistant Staphylococcus epidermidis isolates with reduced vancomycin susceptibility from bloodstream infections in a neonatal intensive care unit. J Med Microbiol. 2020;69(1):41–5.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Liu C, Bayer A, Cosgrove SE, Daum RS, Fridkin SK, Gorwitz RJ, et al. Clinical practice guidelines by the infectious diseases society of america for the treatment of methicillin-resistant Staphylococcus aureus infections in adults and children: executive summary. Clin Infect Dis. 2011;52(3):285–92.

    Article  PubMed  Google Scholar 

  4. 4.

    Tobin CM, Darville JM, Thomson AH, Sweeney G, Wilson JF, MacGowan AP, et al. Vancomycin therapeutic drug monitoring: is there a consensus view? The results of a UK National External Quality Assessment Scheme (UK NEQAS) for antibiotic assays questionnaire. J Antimicrob Chemother. 2002;50(5):713–8.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Wong-Beringer A, Joo J, Tse E, Beringer P. Vancomycin-associated nephrotoxicity: a critical appraisal of risk with high-dose therapy. Int J Antimicrob Agents. 2011;37(2):95–101.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Gonzalez D, Melloni C, Yogev R, Poindexter BB, Mendley SR, Delmore P, et al. Use of opportunistic clinical data and a population pharmacokinetic model to support dosing of clindamycin for premature infants to adolescents. Clin Pharmacol Ther. 2014;96(4):429–37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Kimura T, Sunakawa K, Matsuura N, Kubo H, Shimada S, Yago K. Population pharmacokinetics of arbekacin, vancomycin, and panipenem in neonates. Antimicrob Agents Chemother. 2004;48(4):1159–67.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Ette EI, Williams PJ. Population pharmacokinetics II: estimation methods. Ann Pharmacother. 2004;38(11):1907–15.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Spruill WJ, Wade WE, Cobb HH 3rd. Estimating glomerular filtration rate with a modification of diet in renal disease equation: implications for pharmacy. Am J Health Syst Pharm. 2007;64(6):652–60.

    Article  PubMed  Google Scholar 

  10. 10.

    Brendel K, Dartois C, Comets E, Lemenuel-Diot A, Laveille C, Tranchand B, et al. Are population pharmacokinetic and/or pharmacodynamic models adequately evaluated? A survey of the literature from 2002 to 2004. Clin Pharmacokinet. 2007;46(3):221–34.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Maitre PO, Buhrer M, Thomson D, Stanski DR. A three-step approach combining Bayesian regression and NONMEM population analysis: application to midazolam. J Pharmacokinet Biopharm. 1991;19(4):377–84.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Li ZL, Liu YX, Jiao Z, Qiu G, Huang JQ, Xiao YB, et al. Population pharmacokinetics of Vancomycin in Chinese ICU neonates: initial dosage recommendations. Front Pharmacol. 2018;9:603.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Cockcroft DW, Gault MH. Prediction of Creatinine clearance from serum Creatinine. Nephron. 1976;16(1):31–41.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Purwonugroho TA, Chulavatnatol S, Preechagoon Y, Chindavijak B, Malathum K, Bunuparadah P. Population pharmacokinetics of vancomycin in Thai patients. ScientificWorldJournal. 2012;2012:762649–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Keizer RJ, Karlsson MO, Hooker A. Modeling and simulation workbench for NONMEM: tutorial on Pirana, PsN, and Xpose. CPT Pharmacometrics Syst Pharmacol. 2013;2(6):e50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Deng C, Liu T, Wu K, Wang S, Li L, Lu H, et al. Predictive performance of reported population pharmacokinetic models of vancomycin in Chinese adult patients. J Clin Pharm Ther. 2013;38(6):480–9.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Balch AH, Constance JE, Thorell EA, Stockmann C, Korgenski EK, Campbell SC, et al. Pediatric vancomycin dosing: trends over time and the impact of therapeutic drug monitoring. J Clin Pharmacol. 2015;55(2):212–20.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Anderson BJ, Allegaert K, Van den Anker JN, Cossey V, Holford NH. Vancomycin pharmacokinetics in preterm neonates and the prediction of adult clearance. Br J Clin Pharmacol. 2007;63(1):75–84.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Sanchez JL, Dominguez AR, Lane JR, Anderson PO, Capparelli EV, Cornejo-Bravo JM. Population pharmacokinetics of vancomycin in adult and geriatric patients: comparison of eleven approaches. Int J Clin Pharmacol Ther. 2010;48(8):525–33.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    de Gatta Garcia M d MF, Revilla N, Calvo MV, Dominguez-Gil A, Sanchez Navarro A. Pharmacokinetic/pharmacodynamic analysis of vancomycin in ICU patients. Intensive Care Med. 2007;33(2):279–85.

    CAS  Article  Google Scholar 

  21. 21.

    Revilla N, Martin-Suarez A, Perez MP, Gonzalez FM, Fernandez de Gatta Mdel M. Vancomycin dosing assessment in intensive care unit patients based on a population pharmacokinetic/pharmacodynamic simulation. Br J Clin Pharmacol. 2010;70(2):201–12.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Llopis-Salvia P, Jimenez-Torres NV. Population pharmacokinetic parameters of vancomycin in critically ill patients. J Clin Pharm Ther. 2006;31(5):447–54.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Schleibinger M, Steinbach CL, Topper C, Kratzer A, Liebchen U, Kees F, et al. Protein binding characteristics and pharmacokinetics of ceftriaxone in intensive care unit patients. Br J Clin Pharmacol. 2015;80(3):525–33.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Rybak M, Lomaestro B, Rotschafer JC, Moellering R Jr, Craig W, Billeter M, et al. Therapeutic monitoring of vancomycin in adult patients: a consensus review of the American Society of Health-System Pharmacists, the Infectious Diseases Society of America, and the Society of Infectious Diseases Pharmacists. Am J Health Syst Pharm. 2009;66(1):82–98.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Gwee A, Cranswick N, McMullan B, Perkins E, Bolisetty S, Gardiner K, et al. Continuous Versus Intermittent Vancomycin Infusions in Infants: A Randomized Controlled Trial. Pediatrics. 2019;143(2).

Download references


We would like to express my gratitude to all those who helped me during the writing of this paper. Our deepest gratitude goes first and foremost to Professor Xiaohu Gong for assisting in preparation of this manuscript. Second, we would like to express my heartfelt gratitude to Professor Cheng Cai, Yong Hu and Yihuan Chen, who have instructed and helped me a lot in collecting the patient’s information.


The study design, data acquisition and language editing of this work was supported by the Wu Jieping Foundation of Clinical Pharmacy Branch of Chinese Medical Association (No. 320.6750.19090–1).

Author information




ZJ, FX and HJS conceived and designed the study. ZLL, HJL, and CYW collected the information of the neonates. ZLL, HJL, and CYW performed the experiments. ZLL, HJL, and ZJ wrote the paper. ZLL, HJL, and CYW built and evaluated the model. ZJ, FX and HJS reviewed and edited the manuscript. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Zheng Jiao or Feng Xu or Huajun Sun.

Ethics declarations

Ethics approval and consent to participate

This study was approved by Ethic Committee of Shanghai Children’s Hospital (2015R003-F01). Written informed consent was obtained from guardian.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, Z., Li, H., Wang, C. et al. Establishment of a population pharmacokinetics model of vancomycin in 94 infants with septicemia and its application in individualized therapy. BMC Pharmacol Toxicol 22, 26 (2021).

Download citation


  • Infants septicemia
  • Vancomycin
  • Population pharmacokinetic
  • Monte Carlo
  • Individualized administration