# Population Pharmacokinetics of Vancomycin in Chinese ICU Neonates: Initial Dosage Recommendations

^{1}Department of Pharmacy, Shanghai Children's Hospital, Shanghai Jiao Tong University, Shanghai, China^{2}Department of Pharmacy, Huashan Hospital, Fudan University, Shanghai, China^{3}School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, China^{4}Department of Neonatology, Shanghai Children's Hospital, Shanghai Jiao Tong University, Shanghai, China^{5}Department of Pharmacy, Renmin Hospital of Wuhan University, Wuhan, China^{6}Department of Pharmacy, Gansu Provincial Hospital, Lanzhou, China

The main goal of our study was to characterize the population pharmacokinetics of vancomycin in critically ill Chinese neonates to develop a pharmacokinetic model and investigate factors that have significant influences on the pharmacokinetics of vancomycin in this population. The study population consisted of 80 neonates in the neonatal intensive care unit (ICU) from which 165 trough and peak concentrations of vancomycin were obtained. Nonlinear mixed effect modeling was used to develop a population pharmacokinetic model for vancomycin. The stability and predictive ability of the final model were evaluated based on diagnostic plots, normalized prediction distribution errors and the bootstrap method. Serum creatinine (Scr) and body weight were significant covariates on the clearance of vancomycin. The average clearance was 0.309 L/h for a neonate with Scr of 23.3 μmol/L and body weight of 2.9 kg. No obvious ethnic differences in the clearance of vancomycin were found relative to the earlier studies of Caucasian neonates. Moreover, the established model indicated that in patients with a greater renal clearance status, especially Scr < 15 μmol/L, current guideline recommendations would likely not achieve therapeutic area under the concentration-time curve over 24 h/minimum inhibitory concentration (AUC_{24h}/MIC) ≥ 400. The exceptions to this are British National Formulary (2016–2017), Blue Book (2016) and Neofax (2017). Recommended dose regimens for neonates with different Scr levels and postmenstrual ages were estimated based on Monte Carlo simulations and the established model. These findings will be valuable for developing individualized dosage regimens in the neonatal ICU setting.

## Introduction

After more than 60 years of widespread clinical use, vancomycin remains the gold standard antibiotic prescribed for the treatment of sepsis caused by *coagulase-negative Staphylococci* and methicillin-resistant *Staphylococcus aureus* in neonatal intensive care units (ICU) (Tong et al., 2016). In 2011, guidelines issued by the Infectious Disease Society of America and pediatric-specific guidance recommended targeting trough concentrations of >15 mg/L for critically ill children and >10 mg/L for all other pediatric patients (Liu et al., 2011). Low concentrations may result in less-effective therapy and an increased propensity for the development of bacterial resistance, whereas high concentrations are reported to be associated with nephrotoxicity and ototoxicity, although these toxic effects are less common in neonates (Anderson et al., 2007).

Vancomycin is mainly eliminated by the kidneys and is highly correlated with creatinine clearance. In neonates, the pharmacokinetics of vancomycin are highly variable because of developmental and pathophysiological changes (Stockmann et al., 2015). It is both challenging and imperative to optimize the vancomycin dosage regimen to achieve adequate exposure within a short period of time.

The maximum a *posteriori* Bayesian estimation method has already been used to support vancomycin dosing decisions in adults (Deng et al., 2013; Jacqz-Aigrain et al., 2015) and children (Le et al., 2014). To obtain accurate estimation with this method for individualized therapy, it is crucial that reliable population pharmacokinetic characteristics are known for the target patients. Several pharmacokinetic studies have been conducted in different ethnic groups of adults and showed different clearance (CL) for vancomycin both in Chinese patients [6.05 (2.38–6.06) L/h and 6.06 ± 2.46 L/h] (He et al., 2014; Lin et al., 2015) and in Caucasian patients (4.03 ± 1.7 and 5.83 ± 2.39 L/h) (Guay et al., 1993; Sánchez et al., 2010). Moreover, several population pharmacokinetic studies have been conducted in Caucasian (Seay et al., 1994; Grimsley and Thomson, 1999; Capparelli et al., 2001; Mulla and Pooboni, 2005; Allegaert et al., 2007; Anderson et al., 2007; Marques Minana et al., 2010; Oudin et al., 2011; Mehrotra et al., 2012; Zhao et al., 2013; Frymoyer et al., 2014), Japanese (Kimura, 2002), and Malaysian (Lo et al., 2010) neonates, but a large variability in clearance was observed among groups. The average clearance across all three populations ranged from 0.08 to 0.50 L/h for a neonatal patient with a weight of 3 kg and postmenstrual age (PMA) of 40 weeks.

As little information about vancomycin population pharmacokinetics is known for Chinese patients, the goal of the present study was to establish a population pharmacokinetic model for critically ill Chinese neonates receiving vancomycin therapy and to provide a rational dosage regimen for Chinese neonates.

## Methods

### Patients

Chinese neonates who received vancomycin in the neonatal ICU at Shanghai Children's Hospital (Shanghai, China) between January 2013 and December 2016 were enrolled. Patients included in this study were neonates with PMA ≤ 48 weeks for preterm neonates and with postnatal age (PNA) ≤ 28 days for term neonates. All patients were treated with vancomycin for at least 3 days, and at least one vancomycin concentration was assayed. Patients with extracorporeal membrane oxygenation or continuous renal replacement therapy were excluded from the current study. The following information was retrospectively collected from electronic medical records: gestational age (GA), PMA, PNA, body weight (WT), dosing history and concentration of vancomycin, serum creatinine (Scr) levels, clinical laboratory tests of other renal and hepatic functions, and co-administered medications.

This study was carried out in accordance with the recommendations of the Declaration of Helsinki (2000). The protocol was approved by the Ethics Committee of Shanghai Children's Hospital. Parents or guardians of patients gave their informed consent before enrollment.

The dose of vancomycin (Vancocin, Lilly, S.A, Suzhou, China) was 10–15 mg/kg, which was administrated every 8 h (q8h) or every 12 h (q12h) with a 2-h infusion, in accordance with local protocols. Blood samples were collected 1 h after completion of drug infusion (peak concentration) or half an hour before the start of vancomycin administration (trough concentration) for each patient. Trough and peak concentrations were obtained after at least four repeated doses. Potential outlier data points (observations) were identified by employing conditional weighted residual (CWRES) results outside a range of ± 6 (Byon et al., 2013).

### Bioassay

Concentrations of vancomycin were determined using the fluorescence polarization immunoassay with an ARCHITECT i2000SR (Abbott Laboratories, Chicago, IL, USA). The limit of detection was 1 mg/L, and the calibration range of this assay was 3 to 50 mg/L. The intra-day and inter-day coefficients of variation (CVs) were all < 20%.

Scr assays were performed with the enzymatic method and were analyzed with a 7180 Automatic Analyzer (Hitachi High-Tech Science Systems Corporation, Tokyo, Japan). The calibration range was from 3 to 100 mg/L. The intra-day and inter-day CVs were all < 3.75%. Creatinine can pass through the placenta and many endogenous factors can influence creatinine determination in neonates, which can limit the concentration of creatinine detected.

### Model Building

Population pharmacokinetic modeling was performed using the NONMEM program (Version 7.4; Icon Inc, PA, USA) compiled with gFortran (Version 4.9.2; http://www.gfortran.org). R (Version 3.3.1; http://www.r-project.org) and the Xpose package (Version 4.5.3; http://xpose.sourceforge.net) were used for visual diagnosis. The first-order conditional estimation method with η–ε interaction (FOCE-I) was used throughout the model-building procedure (Wählby et al., 2001).

#### Base Model

The vancomycin concentration data were fitted by the one or two compartment with a first-order elimination model. NONMEM subroutines were specified as ADVAN1-TRANS2 or ADVAN3-TRANS4, respectively.

An exponential model (Equation 1) was used to account for inter-individual variability (IIV):

where TV(P) is the typical value of the pharmacokinetic parameter, and P_{i} refers to the pharmacokinetic parameter of the *i*th patient with random variable η_{i}, which is normally distributed with a mean of zero and variance of ω^{2}.

Residual unexplained variability was tested by an additive model (Equation 2), exponential model (Equation 3), or a combined additive and exponential model (Equation 4):

where Y represents the observation, IPRED is the individual predicted concentration, and ε is a symmetrically distributed variable, with a mean of zero and variance of σ^{2}.

#### Covariate Models

The continuous covariates, including WT, GA, PNA, PMA, Scr, blood urea nitrogen, serum albumin concentration, aspartate transaminase and alanine transaminase, and categorical covariates, including gender, concomitant drug (ceftriaxone, meropenem, gentamicin, furosemide, ibuprofen, and dexamethasone), were screened for their influence on clearance and the volume of distribution.Covariate screening was performed according to the following four steps:

**Step 1** Body weight and age have significant impacts on the pharmacokinetics of vancomycin in neonates (Wallis and Williamson, 2011; Jacqz-Aigrain et al., 2015) and physical maturation is a time-dependent characteristic that must be considered in neonates (Schmidt and Derendorf, 2014). Therefore, body weight (WT) and age (PNA, GA, and PMA) were screened first. Four different models based on allometric scaling were tested using Equation (5):

where COV_{median} is the median of the covariate, MF is the maturation factor that is defined as the process of becoming mature. The model displaying the best fit was selected for further analysis.

*Model I*: In the simplest exponent model, the exponent θ was estimated, and MF was fixed to 1, indicating that maturation was not considered. This model is shown as Equation (6):

*Model II*: For the maturation model (Holford et al., 2013), the exponent θ was assigned a fixed value of 0.75, and MF was calculated according to Equation (7):

where TM_{50} is the age (in terms of GA, PMA, or PNA) at which clearance maturation reaches 50% of that of adults, and Hill is the slope parameter for the sigmoid E_{max} maturation model.

*Model III*: This is referred to as the WT-dependent exponent model (Holford et al., 2013):

*Model IV*: This is referred to as the age-dependent exponent model (Ding et al., 2015):

where θ_{0} is the value of the exponent at a theoretical WT of zero (Equation 8) or at birth (0 years) (Equation 9), *k*_{max} is the maximum decrease of the exponent, *k*_{50} is the WT (Equation 8) or age (Equation 9) at which a 50% decrease relative to the maximum decrease is attained, and the Hill coefficient is used to determine the steepness of the sigmoid decline.

**Step 2** In previous studies, renal function has been identified as an important covariate (Grimsley and Thomson, 1999; Capparelli et al., 2001; Kimura, 2002; Oudin et al., 2011; Mehrotra et al., 2012; Zhao et al., 2013; Derschmills et al., 2014; Frymoyer et al., 2014). We thus investigated and Scr by using the exponential model (Equation 6) and a linear model (Equation 10):

The one displaying the best fit was used for further analysis.

**Step 3**. The remaining covariates were then accessed sequentially by forward inclusion and backward elimination approaches using the exponential model (Equation 8) or linear model (Equation 10) for continuous variables and a proportional model (Equation 11) for categorical variables, such as gender in the P_{i} of vancomycin:

Where θ is the fractional change in TV(P) for males.

**Step 4** Taking into account the rapid variation in the physical status of neonates, the covariates identified above could be defined as time-varying covariates to illustrate IIV. The time-varying covariates model splits the individual covariate effects into baseline and change-from-baseline effects. Two different models based on covariate scaling of the pharmacokinetic parameters were tested using Equations (12) and (13):

where BCOV is the baseline value of the covariate; BCOV_{median} is the median of the baseline value of covariate; θ_{BCOV} describes the effect of IIV. DCOV is equal to the current covariate value minus BCOV at each time point and corresponds to the fractional change in the typical value with each unit difference in BCOV relative to BCOV_{median.} θ_{DCOV} describes the effect of covariate variation within an individual and is the fractional change in the typical value with individual changes in COV, and (η, P_{i}) refers to variable with a mean of zero and variance of ω^{2} that describes the random effect of P_{i}.

If θ_{BCOV} and θ_{DCOV} are different, DCOV is fixed to zero, and an additional variance parameter (η_{COV}, P_{i}) that accounts for IIV to influence the covariates for the population parameter estimates is included (Equation 14).

θ_{COV} is the parameter estimation value of the covariate, COV is the value of the covariate, COV_{median} is the median of the covariate, η_{COV} is random variable (with zero mean and variance ω^{2}).

#### Model Selection Criteria

Structural models were selected through Akaike information criteria (AIC) and Bayesian information criteria (BIC) calculated using Pirana software (version 2.9.0, http://www.pirana-software.com/) (Keizer et al., 2011) Models with lower AIC and BIC values were considered superior (Byon et al., 2013).

Nested models in covariate screening were compared statistically using a likelihood ratio test on the differences in the objective function value (OFV). The change was considered significant if the decrease in OFV was >3.83 (χ^{2}*, df* = 1, *P* < 0.05) for the forward inclusion step, and the increase in OFV was >6.63 (χ^{2}*, df* = 1, *P* < 0.01) for the backward elimination step.

The covariates in the model were selected based on physiological plausibility of parameter estimates, goodness-of-fit plots, and statistical significance.

### Model Evaluation

The performance of the final model was first evaluated by visual inspection of diagnostic goodness-of-fit plots. Goodness-of-fit plots included the following scatterplots: observation (DV) vs. individual prediction (IPRED), DV vs. population prediction (PRED), conditional weighted residual errors (CWRES) vs. time, and CWRES vs. PRED (Hooker et al., 2007).

The robustness of the model was assessed using a nonparametric bootstrap (Ette, 1997), with repetition of 2000 NONMEM runs of the final model. The bootstrap median parameter values and the percentile bootstrap 95% confidence intervals were compared with the respective values estimated from the final model.

Normalized prediction distribution error (NPDE) (Comets et al., 2008) was used to evaluate the predictive performance of the model based on a Monte Carlo simulation with the R package (version 2.0, http://www.npde.biostat.fr/). NPDE results were summarized graphically using (1) quantile-quantile plot of the NPDE, (2) a histogram of the NPDE, (3) scatterplot of NPDE vs. time, and (4) scatterplot of NPDE vs. PRED. If the predictive performance is satisfied, the NPDE will follow a normal distribution (Shapiro-Wilk test) with a mean value of zero (*t*-test) and a variance of one (Fisher's test).

### Dosing Regimen Design

The final established population pharmacokinetic model was used to obtain dosing regimens for vancomycin to reach AUC_{24h}/MIC ≥ 400 which is known to produce an effective therapeutic outcome (Jacqz-Aigrain et al., 2015). When MIC = 1 mg/L, the daily dose can be calculated by the final model for determining the clearance and Equation 14:

Simulations were performed for virtual patients with various levels of renal function and ages, to determine the most appropriate scheme to satisfy the therapeutic criteria. For this purpose, 1,000 replicates of each scenario were simulated by the Monte Carlo method and were completed by the $SIMULATION modules in NONMEM software.

Virtual patients were designated as having PMAs covering a 2-week window between 28 and 44 weeks, and their corresponding WTs were calculated according to the World Health Organization growth chart for infants (Centers for Disease Control Prevention, 2009). Patients were designated a Scr level of 15, 20, or 35 μmol/L and with an age of PNA 7 days (≤1 week) or PNA 15 days (>2 weeks).

Vancomycin dosage recommendations are highly variable as illustrated by differences in various guidelines (Jacqz-Aigrain et al., 2015). The dose regimens were then compared to the guidelines used in different regions or medical centers, including the FDA's labeled dosage regimen, British National Formulary (BNF) for Children 2016–2017 and Blue Book 2016 from the UK, Neonatal Formulary 2015 from the Europe, Red Book 2015, Pediatric & Neonatal Dosage Handbook 2015–2016, and Neofax 2017 (http://neofax.micromedexsolutions.com/neofax/neofax.php) from the US.

## Results

### Patients

Data from 165 vancomycin measurements, with a trough concentration for 75 of the measurements and a peak concentration for 90 of the measurements, from 80 subjects were available to perform the population pharmacokinetics analysis. For each subject, an average of two samples were obtained. The GA range was 25.7 weeks to 41.1 weeks, with a mean WT of 2.87 kg. Of these patients, 59% were preterm infants and 57% had respiratory tract infections. All the observation were included, and no outlier records were identified. Clinical characteristics of the neonates included in the analysis are summarized in Table 1.

### Model Building

A one-compartment model with first order elimination described the pharmacokinetic of vancomycin. The residual unexplained variability was described best by a proportional model. As only peak and trough samples were collected, the relative standard error (RSE) of the IIV of the volume of distribution was poor (>50%) and was not estimated.

For the first step of covariate screening, several models were tested, and the results are shown in Table 2. Among the four models examined, the simple exponent model (Model I) and maturation model (Model II) had lower Akaike information criteria and Bayesian information criteria than the ADE model (Model III) and BDE model (Model IV). The maturation model had a condition value of 23,082, much greater than 1,000, indicating model instability (Byon et al., 2013). Moreover, the RSE of most PK parameters in the maturation model were more than 100%, implying inaccuracy of the model parameter estimates. Therefore, the simple exponent model was employed in further analyses.

Second, Scr and WT were identified as significant covariates and were thus retained in the model. Further incorporation of time-varying covariates did not improve the model performance, which could be attributed to the short treatment duration of vancomycin.

The final model for vancomycin clearance was represented by Equation 15, and WT was added to volume of distribution for physiologic plausibility as shown in Equation 16.

where CL is vancomycin clearance, WT is current body weight in neonates, Scr is serum creatinine, V is volume of distribution for vancomycin. The final model parameter estimates are shown in Table 3.

**Table 3**. Population pharmacokinetic estimates of vancomycin of final model and Bootstrap evaluation.

### Model Evaluation

Goodness-of-fit plots for the base model and final model are shown in Figure 1. Compared with the base model, the final model showed no obvious bias or significant trends within these scatterplots. Moreover, the data fitting for the final model was much improved relative to that of the base model.

**Figure 1**. Diagnostic goodness-of-fit plots for the base model (1) and the final model (2). **(A)** The individual predicted concentration (IPRED) vs. the observed concentration. **(B)** The population predicted concentration (PRED) vs. the observed concentration. **(C)** The PRED vs. the conditional weighted residual errors (CWRES). **(D)** The time after dose vs. CWRES. The black solid lines are the reference lines, and red solid lines are the loess smooth lines.

The results from the bootstrap procedure are shown in Table 3. The median values from the bootstrap procedure were close to the parameter estimates from the NONMEM, with <5% bias. In addition, more than 99% of the bootstrap runs were successful, indicating that the model was stable.

The NPDE distribution and histogram are presented in Figure 2. The assumption of a normal distribution for the differences between predictions and observations was acceptable. The quantile-quantile plots and histogram also confirmed the normality of the NPDE (Figure 2).

**Figure 2**. Normalized prediction distribution error (NPDE) for the final model. Quantile-quantile plots of NPDE vs. the expected standard normal distribution (upper left). Histogram of NPDE values with the standard normal distribution overlayed (upper right). Scatter plot of the time vs. NPDE (lower left). Scatterplot of predictions vs. NPDE (lower right).

### Dosing Regimen Design

The dosage regimens recommended by current guidelines and estimated by the established model are displayed in Figure 3 and Table 4. The guideline schemes from the BNF for children, the Blue Book and the Neonatal Formulary are nearly in accordance with the 15% to 85% dosage interval from the present model that achieves a therapeutic target of 400 ≤ AUC_{24h}/MIC <800. The targeted vancomycin concentration based on the FDA labeled dosage, the Red Book, the Pediatric & Neonatal Dosage Handbook, and Neofax was not likely sufficient, especially for neonates with Scr of 15 μmol/L. Moreover, Neonatal Formulary shows an overdose for a subpopulation of neonates with Scr of 35 μmol/L, which might indicate an increased risk of toxicity.

**Figure 3**. Vancomycin dosage regimen recommended by the six guidelines relative to the regimen recommended by our final model in typical patients when AUC_{24h}/MIC ≥ 400. The six guides correspond to (1) the FDA labeled dosage (2017), (2) the British National Formulary (2016–2017) and the Blue Book (2016), (3) the Neonatal Formulary (2015), (4) the Red Book (2015), (5) the Pediatric and Neonatal Dosage Handbook (2015-2016), and (6) Neofax (2017). The red and blue lines in (1) and (6) refer to dosage guidelines for patients with a PNA of 7 days and 15 days respectively. The blue lines in (2),(3),(4), and (5) refers to dosage for PNA of both 7 and 15 days. The dark blue smooth curve represents the mean dosage for the present study, and the light blue ribbon corresponds to the 15–85% dosage interval.

## Discussion

In the present study, we first developed a population pharmacokinetic model for vancomycin in Chinese neonates in the ICU and found that WT and Scr levels have significant influences on clearance. Little obvious ethnic difference of vancomycin clearance was shown in Asian and Caucasian neonates from our study.

Differences in vancomycin pharmacokinetics have been noted between Asian and Caucasian populations based on two recently published studies (Lin et al., 2015, 2016). To expand this analysis, we first looked at a group of 12 previous studies that examined neonatal populations from various countries. The population within some of these studies were heterogeneous. Patients with extracorporeal membrane oxygenation were included in the model of Mulla et al. (Mulla and Pooboni, 2005), and the majority of neonates used to develop the models of Allegaert et al. (2007) and Capparelli et al. (2001) were preterm. Due to differences in the physiologic development of neonates within these populations, models from these three studies were excluded from the comparison. Additionally, neonates within 6 other studies (Seay et al., 1994; Grimsley and Thomson, 1999; Lo et al., 2010; Marques Minana et al., 2010; Oudin et al., 2011; Zhao et al., 2013) were smaller, with a mean WT <2 kg. Given the relative weight of neonates in our study and the impact of WT on vancomycin pharmacokinetic parameters, these six studies were excluded from comparisons. The three remaining studies (Kimura, 2002; Mehrotra et al., 2012; Frymoyer et al., 2014) were included in our analysis of ethnic differences regarding vancomycin clearance.

The clearance of a standardized patient, as determined by the different models, was calculated to investigate differences in vancomycin clearance relative to ethnicity. The standardized patient had a WT of 2.8 kg, PMA of 37 weeks with different Scr levels.

As shown in Table 5, across varying Scr levels ranging from 20 μmol/L to 50 μmol/L. there was 27~39 % difference in vancomycin clearance between the current study and the study by Mehrotra (Mehrotra et al., 2012), but comparable to another study (<13%) which was also conducted in US (Frymoyer et al., 2014). This information does not support the conclusion that there are obvious differences in vancomycin clearance between Chinese and Caucasians. However, the estimated clearance was much higher than the study conducted in 19 Japanese neonates (Kimura, 2002). The reason was unclear. The present study might be under-power to conclure ethnic impact on vancomycin PK in neonates. Factors in the current study, including analytical methods, and disease progression could affect the assessment of ethnic differences.

Dosage recommendations by the label and reference books are variable as shown in Table 6. The variabilities were attributed to the different covariates considered within these recommendations. Body weight is the most notable covariate for vancomycin dosing found within all references. Age-based (as PMA, PNA, and GA) dosing is also included in most of the references, such as the FDA label recommendation, BNF for children, Blue Book, Neonatal Formulary, Pediatric and Neonatal Dosage Handbook and Neofax. Dosing based on Scr is only included in 2 references, the Red Book, and Pediatric & Neonatal Dosage Handbook and only covered Scr level > 61.9 μmol/L. However, Scr was identified to have large impacts on the CL of vancomycin in all previous population pharmacokinetic studies as well as the present study.

Based off data from our final model, we found that the current recommended doses of vancomycin from FDA labeled dosage, Red Book, Pediatric & Neonatal Dosage Handbook, and Neofax may be inadequate to meet a treatment target of AUC_{24h}/MIC ≥ 400, especially for patients with a greater renal clearance status, especially Scr < 15 μmol/L.

This finding is consistent with several previous studies (Krivoy et al., 1998; Liu et al., 2011; Abdel et al., 2015; Zhao et al., 2015). They found that the usual recommended dose of 60 mg/kg/day did not achieve vancomycin pharmacodynamic targets in most patients. Silva, D C et al (Silva et al., 2012) reported that a vancomycin dose of 81 mg/kg/day was required to achieve an AUC_{24h}/MIC > 400 in 56% of patients. Doses as high as 120 mg/kg/day were also recommended to improve the therapeutic pharmacodynamic targets (Abdel et al., 2015). Higher than usual vancomycin doses may be required to treat patients with severe Gram-positive infections.

There are several limitations to this study. The current guidelines or consensus for therapeutic drug management (TDM) of vancomycin in the United States (Rybak et al., 2009), Japan (Matsumoto et al., 2013) and China (Ye et al., 2016) recommends that trough concentrations should be collected at regular intervals in clinical settings. Therefore, only peak and trough concentrations were collected in our study and a one compartment model was applied for the structural model even though vancomycin is more normally modeled with a two-compartment model. The simplification to a one compartment model may lead to deviation of clearance estimation. However, the bias usually is <20% and does not obviously affect the estimation of the area under the curve (AUC). (Ling et al., 2014) (Kowalski and Hutmacher, 2001). The recommended regimen was based on a study population from our hospital. Therefore, generalization to other ICU neonates treated with vancomycin, especially if their covariate characteristics lie outside the range of our study population, would require additional investigations.

Moreover, several cofactors may affect the clinical outcomes of patients, such as disease progression, baseline weight, gestational age, and medication interactions. Therefore, with a limited number of patients enrolled in the study, we did not compare the clinical outcomes between those who had different exposure levels. Furthermore, as a pharmacodynamics indicator, the MIC value was obtained for only a few patients, such that we could not build a population pharmacokinetic-pharmacodynamics model, which would have been a better predictor of vancomycin's therapeutic effect.

This study offers initial vancomycin dosing regimen with varying degrees of PMA, WT, and serum creatinine for neonates. For patients with complex disease conditions, Bayesian approaches might be used to provide individualized dose recommendations instead of look-up tables or nomograms. Dose calculators and other decision support tools based on population pharmacokinetic (PPK) models (Fuchs et al., 2013) could contribute to simplifying the complex Bayes calculations and making them more intuitive to the user in clinical practice.

Recently, a nonparametric (NP) population modeling approach was reported to have advantages in patient's individual dosing adjustment, which permits development of dosage regimens to hit desired therapeutic targets with maximum precision (Jelliffe et al., 2011; Neely et al., 2018). This will be further investigated.

## Conclusion

In summary, this study has built a population pharmacokinetic model of vancomycin for Chinese neonates. WT and Scr levels were the important covariates, which affect clearance. Moreover, this study found no obvious differences in the clearance of vancomycin comparing Caucasian and Chinese neonates. For patients with normal renal function, the dosing recommendations are likely not sufficient based on the target of AUC_{24h}/MIC ≥400, with the exception of British National Formulary (2016–2017), Blue Book (2016) and Neofax (2017). However, these sources provide little information on dosing adjustments based on patient renal function, thus our model provides a method for adjusting the vancomycin dose accordingly.

## Author Contributions

ZL, ZJ, and HS conceived and designed the study. YL, SW, CW, HS and GQ collected the information of the neonates. YL, JH, YX, and WH performed the experiments. ZL, ZJ, and YL wrote the paper. YL, SW, YX and CW built the model and evaluated. ZL, YL, and ZJ reviewed and edited the manuscript. All authors read and approved the manuscript.

## Funding

This work was supported by the National Natural Science Foundation of China (No. 81603199), Shanghai Municipal Population and Family Planning Commission (No.20154Y0019), Weak Discipline Construction Project of Shanghai Municipal Population and Family Planning Commission (No.2016ZB0305–01& 2016ZB0301–01), and Cross-fund for medical and engineering workers in Shanghai Jiaotong University (YG2015QN29).

## Conflict of Interest Statement

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.

## Acknowledgments

I would like to express my gratitude to all those who helped me during the writing of this paper. My deepest gratitude goes first and foremost to Professor Feng Xu for assisting in preparation of this manuscript. Second, I 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.

## References

Abdel, H. O., Al, O. S., Nazer, L. H., Mubarak, S., and Le, J. (2015). Vancomycin pharmacokinetics and predicted dosage requirements in pediatric cancer patients. *J. Oncol. Pharm. Pract.* 22. 448–453. doi: 10.1177/1078155215591386

Allegaert, K., Anderson, B. J., Jn, V. D. A., Vanhaesebrouck, S., and De, Z. F. (2007). Renal drug clearance in preterm neonates: relation to prenatal growth. *Ther. Drug Monit.* 29, 284–291. doi: 10.1097/FTD.0b013e31806db3f5

Anderson, B. J., Allegaert, K., Jn, V. D. A., Cossey, V., and Holford, N. H. (2007). Vancomycin pharmacokinetics in preterm neonates and the prediction of adult clearance. *Br. J. Clin. Pharmacol.* 63, 75–84. doi: 10.1111/j.1365-2125.2006.02725.x

Byon, W., Smith, M. K., Chan, P., Tortorici, M. A., Riley, S., Dai, H., et al. (2013). Establishing best practices and guidance in population modeling: an experience with an internal population pharmacokinetic analysis guidance. *CPT Pharmacometrics Syst. Pharmacol.* 2:e51. doi: 10.1038/psp.2013.26

Capparelli, E. V., Lane, F. R., Romanowski, G. L., Pharm, M. F., Murray, W., Sousa, P., et al. (2001). The influences of renal function and maturation on vancomycin elimination in newborns and infants. *J. Clin. Pharmacol.* 41, 927–934. doi: 10.1177/009127000104100901

Centers for Disease Control and Prevention (2009). *WHO Child Growth Standards*. Available online at: http://www.who.int/childgrowth/en [EB/OL] (Accessed December 9, 2017).

Comets, E., Brendel, K., and Mentré, F. (2008). Computing normalised prediction distribution errors to evaluate nonlinear mixed-effect models: the npde add-on package for R. *Comput. Methods Programs Biomed.* 90, 154–166. doi: 10.1016/j.cmpb.2007.12.002

Deng, C., Liu, T., Wu, K., Wang, S., Li, L., Lu, H., et al. (2013). Predictive performance of reported population pharmacokinetic models of vancomycin in chinese adult patients. *J. Clin. Pharm. Ther.* 38, 480–489. doi: 10.1111/jcpt.12092

Derschmills, D., Bengry, T., Akierman, A., Alshaikh, B., and Yusuf, K. (2014). Assessment of initial vancomycin dosing in neonates. *Paediatr. Child Health* 19:e30–e34.

Ding, J., Wang, Y., Lin, W., Wang, C., Zhao, L., Li, X., et al. (2015). A population pharmacokinetic model of valproic acid in pediatric patients with epilepsy: a non-linear pharmacokinetic model based on protein-binding saturation. *Clin. Pharmacokinet.* 54, 305–317. doi: 10.1007/s40262-014-0212-8

Ette, E. I. (1997). Stability and performance of a population pharmacokinetic model. *J. Clin. Pharmacol.* 37, 486–495.

Frymoyer, A., Hersh, A. L., Elkomy, M. H., Gaskari, S., Su, F., Drover, D. R., et al. (2014). Association between vancomycin trough concentration and area under the concentration-time curve in neonates. *Antimicrob. Agents Chemother.* 58, 6454–6461. doi: 10.1128/AAC.03620-14

Fuchs, A., Csajka, C., Thoma, Y., Buclin, T., and Widmer, N. (2013). Benchmarking therapeutic drug monitoring software: a review of available computer tools. *Clin. Pharmacokinet*. 52, 9–22. doi: 10.1007/s40262-012-0020-y

Grimsley, C., and Thomson, A. H. (1999). Pharmacokinetics and dose requirements of vancomycin in neonates. *Arch. Dis. Child. Fetal Neonatal Ed.* 81, 221–227.

Guay, D. R., Vancebryan, K., Gilliland, S., Rodvold, K., and Rotschafer, J. (1993). Comparison of vancomycin pharmacokinetics in hospitalized elderly and young patients using a Bayesian forecaster. *J. Clin. Pharmacol.* 33, 918–922. doi: 10.1002/j.1552-4604.1993.tb01922.x

He, X. R., Liu, Z. H., Ji, S. M., Liu, T. T., Li, L., Zhou, T. Y., et al. (2014). Population pharmacokinetics of vancomycin and prediction of pharmacodynamics in the Chinese people. *Acta Pharm. Sin.* 49, 1528–1535.

Holford, N., Heo, Y. A., and Anderson, B. (2013). A pharmacokinetic standard for babies and adults. *J. Pharm. Sci.* 102, 2941–2952. doi: 10.1002/jps.23574

Hooker, A. C., Staatz, C. E., and Karlsson, M. O. (2007). Conditional weighted residuals (cwres): a model diagnostic for the foce method. *Pharm. Res.* 24, 2187–2197. doi: 10.1007/s11095-007-9361-x

Jacqz-Aigrain, E., Leroux, S., Zhao, W., van den Anker, J. N., and Sharland, M. (2015). How to use vancomycin optimally in neonates: remaining questions. *Expert Rev. Clin. Pharmacol.* 8, 635–648. doi: 10.1586/17512433.2015.1060124

Jelliffe, R., Neely, M., Schumitzky, A., Bayard, D., Van Guilder, M., Botnen, A., et al. (2011). Nonparametric population modeling and Bayesian analysis. *Pharmacol. Res*. 64:426. doi: 10.1016/j.phrs.2011.04.008

Keizer, R. J., Van, B. M., Beijnen, J. H., Schellens, J. H., and Huitema, A. D. (2011). Pirana and pcluster: a modeling environment and cluster infrastructure for nonmem. *Comput. Methods Programs Biomed.* 101, 72–79. doi: 10.1016/j.cmpb.2010.04.018

Kimura, T. (2002). Population pharmacokinetics of antibiotics in neonates. *Jpn. J.Antibiot.* 55, 337. doi: 10.1128/AAC.48.4.1159-1167.2004

Krivoy, N., Peleg, S., Postovsky, S., and Arush, M. W. B. (1998). Pharmacokinetic analysis of vancomycin in steady state in pediatric cancer patients. *Eur. Paediatr. Haematol. Oncol.* 15, 333–338.

Kowalski, K. G., and Hutmacher, M. M. (2001). Design evaluation for a population pharmacokinetic study using clinical trial simulations: a case study. *Stat. Med.* 20, 75–91.

Le, J., Ngu, B., Bradley, J. S., Murray, W., Nguyen, A., Nguyen, L., et al. (2014). Vancomycin monitoring in children using bayesian estimation. *Ther. Drug Monit.* 36, 510–518. doi: 10.1097/FTD.0000000000000039

Lin, W. F., Liu, S. S., Yang, T. Y., Win, M. F., Lin, S. W., Huang, C. F., et al. (2015). A larger dose of vancomycin is required in adult neurosurgical intensive care unit patients due to augmented clearance. *Ther. Drug Monit.* 37, 609–618. doi: 10.1097/FTD.0000000000000187

Lin, W. W., Wu, W., Jiao, Z., Lin, R. F., Jiang, C. Z., Huang, P. F., et al. (2016). Population pharmacokinetics of vancomycin in adult chinese patients with post-craniotomy meningitis and its application in individualised dosage regimens. *Eur. J. Clin. Pharmacol.* 72, 29–37. doi: 10.1007/s00228-015-1952-6

Ling, J., Qian, L. X., Ding, J. J., and Jiao, Z. (2014). Effects of multiple-trough sampling design and algorithm on the estimation of population and individual pharmacokinetic parameters. *Acta Pharm. Sin.* 49, 686–694.

Liu, C., Bayer, A., Cosgrove, S. E., Daum, R. S., Fridkin, S. K., Gorwitz, R. J., et al. (2011). 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.* 52, 285–292. doi: 10.1093/cid/cir034

Lo, Y. L., Sc, V. H. J., Lim, C. T., Lee, T. C., and Charles, B. G. (2010). Population pharmacokinetics of vancomycin in premature malaysian neonates: identification of predictors for dosing determination. *Antimicrob. Agents Chemother.* 54, 2626–2632. doi: 10.1128/AAC.01370-09

Marqués-Miñana, M. R., Saadeddin, A., and Peris, J. E. (2010). Population pharmacokinetic analysis of vancomycin in neonates. A new proposal of initial dosage guideline. *Br. J. Clin. Pharmacol.* 70, 713–720. doi: 10.1111/j.1365-2125.2010.03736.x

Matsumoto, K., Takesue, Y., Ohmagari, N., Mochizuki, T., Mikamo, H., Seki, M., et al. (2013). Practice guidelines for therapeutic drug monitoring of vancomycin: a consensus review of the Japanese Society of Chemotherapy and the Japanese Society of Therapeutic Drug Monitoring. *J. Infect.Chemother.* 19, 365–380. doi: 10.1007/s10156-013-0599-4

Mehrotra, N., Tang, L., Phelps, S. J., and Meibohm, B. (2012). Evaluation of vancomycin dosing regimens in preterm and term neonates using montecarlo simulations. *Pharmacotherapy* 32, 408–419. doi: 10.1002/j.1875-9114.2012.01029.x

Mulla, H., and Pooboni, S. (2005). Population pharmacokinetics of vancomycin in patients receiving extracorporeal membrane oxygenation. *Br. J. Clin. Pharmacol.* 60, 265–275. doi: 10.1111/j.1365-2125.2005.02432.x

Neely, M. N., Kato, L., Youn, G., Kraler, L., Bayard, D., van Guilder, M., et al. (2018). Prospective trial on the use of trough concentration versus area under the curve to determine therapeutic vancomycin dosing. *Antimicrob Agents Chemother*. 62:e02042-17. doi: 10.1128/AAC.02042-17

Oudin, C., Vialet, R., Boulamery, A., Martin, C., and Simon, N. (2011). Vancomycin prescription in neonates and young infants: toward a simplified dosage. *Arch. Dis. Child. Fetal Neonatal Ed.* 96, F365–F370. doi: 10.1136/adc.2010

Rybak, M., Lomaestro, B., Rotschafer, J. C., Moellering, R. Jr., Craig, W., Billeter, M., et al. (2009). 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.* 66, 82–98. doi: 10.2146/ajhp080434

Sánchez, J. L., Dominguez, A. R., Lane, J. R., Anderson, P. O., Capparelli, E. V., and Cornejo-Bravo, J. M. (2010). Population pharmacokinetics of vancomycin in adult and geriatric patients: comparison of eleven approaches. *Int. J. Clin. Pharmacol. Ther.* 48, 525–533.

Seay, R. E., Brundage, R. C., Jensen, P. D., Schilling, C. G., and Edgren, B. E. (1994). Population pharmacokinetics of vancomycin in neonates. *Clin. Pharmacol. Ther.* 56, 169–175.

Silva, D. C., Seixas, G. T., Araujo, O. R., Arduini, R. G., Carlesse, F. A., and Petrilli, A. S. (2012). Vancomycin serum concentrations in pediatric oncologic/hematologic intensive care patients. *Braz. J. Infect. Dis.* 16, 361–365. doi: 10.1016/j.bjid.2012.06.011

Stockmann, C., Hersh, A. L., Roberts, J. K., Bhongsatiern, J., Korgenski, E. K., Spigarelli, M. G., et al. (2015). Predictive performance of a vancomycin population pharmacokinetic model in neonates. *Infect. Dis.Ther.* 4, 187–198. doi: 10.1007/s40121-015-0067-9

Tong, S. Y. C., Nelson, J., Paterson, D. L., Fowler, V. G. Jr., Howden, B. P., Cheng, A. C., et al. (2016). Camera2 – combination antibiotic therapy for methicillin-resistant *Staphylococcus aureus* infection: study protocol for a randomised controlled trial. *Trials* 17:170. doi: 10.1186/s13063-016-1295-3

Wählby, U., Jonsson, E. N., and Karlsson, M. O. (2001). Assessment of actual significance levels for covariate effects in nonmem. *J. Pharmacokinet Pharmacodyn.* 28, 231–252. doi: 10.1023/A:1011527125570

Wallis, S., and Williamson, K. A. (2011). Vancomycin in neonates – time to review current dosing recommendations? *Arch. Dis. Child. Fetal Neonatal Ed.* 96:Fa41. doi: 10.1136/archdischild.2011.300164.82

Ye, Z. K., Chen, Y. L., Chen, K., Zhang, X. L., Du, G. H., He, B., et al. (2016). Therapeutic drug monitoring of vancomycin: a guideline of the division of therapeutic drug monitoring, Chinese Pharmacological Society. *J. Antimicrob. Chemother.* 71, 3020–3025. doi: 10.1093/jac/dkw254

Zhao, W., Lopez, E., Biran, V., Durrmeyer, X., Fakhoury, M., and Jacqzaigrain, E. (2013). Vancomycin continuous infusion in neonates: dosing optimisation and therapeutic drug monitoring. *Arch. Dis. Child.* 98, 449–453. doi: 10.1136/archdischild-2012-302765

Keywords: neonate, vancomycin, population pharmacokinetic, Monte Carlo simulation, individualized therapy

Citation: Li Z, Liu Y, Jiao Z, Qiu G, Huang J, Xiao Y, Wu S, Wang C, Hu W and Sun H (2018) Population Pharmacokinetics of Vancomycin in Chinese ICU Neonates: Initial Dosage Recommendations. *Front. Pharmacol*. 9:603. doi: 10.3389/fphar.2018.00603

Received: 31 January 2018; Accepted: 21 May 2018;

Published: 26 June 2018.

Edited by:

Jun Shi, Roche, SwitzerlandReviewed by:

Tamorah Rae Lewis, Children's Mercy Hospital, United StatesWei Zhao, Shandong University, China

Copyright © 2018 Li, Liu, Jiao, Qiu, Huang, Xiao, Wu, Wang, Hu and Sun. 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 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: Zheng Jiao, zjiao@fudan.edu.cn

Hua-jun Sun, sunhj1@shchildren.com.cn

^{†}These authors have contributed equally to this work and co-first authors.