Development of a Population Pharmacokinetic Model of Vancomycin and its Application in Chinese Geriatric Patients with Pulmonary Infections.

BACKGROUND: The optimal use of vancomycin in the elderly requires information about the drug's pharmacokinetics and the influence of various factors on the drug's disposition. However, because of sampling restrictions, it is often difficult to perform traditional pharmacokinetic studies in elderly patients.
OBJECTIVE: This study was conducted to estimate the population pharmacokinetics of vancomycin in Chinese geriatric patients (age ≥ 65 years) with pulmonary infections and to explore the clinical application of this information for vancomycin dose individualization.
METHODS: The steady-state trough concentrations were retrospectively collected from January 2011 to December 2016 and were analyzed using the nonlinear mixed-effect model software. The final model was evaluated using the bootstrap method, goodness-of-fit plots and the normalized prediction distribution error method. MAIN OUTCOME MEASURE: Model parameters and prediction error.
RESULTS: A total of 125 steady-state trough concentrations from 70 patients were retrospectively collected. A one-compartment model was established. The final model was depicted as clearance (CL) [L/h] = 2.45 × (CLCR/56.28) × 0.542; volume of distribution (Vd) [L] = 154. The creatinine clearance (CLCR) was identified as the most significant covariate in the final model. The typical values of CL and Vd in the final model were 2.45 L/h and 154 L, respectively. Model validation outcomes showed that the final model was stable and had satisfactory prediction performance.
CONCLUSION: A population pharmacokinetic model was established to estimate the pharmacokinetics characteristics of Chinese geriatric patients with pulmonary infections, and this model can be used to develop an initial vancomycin dosing regimen for geriatric patients.

Key Points

Introduction

Pulmonary infections are common in hospitalized patients. Vancomycin is the first choice for the treatment of pulmonary infections caused by Gram-positive staphylococcus and is considered the last defense against Gram-positive bacterial infections [1]. These infections are usually more severe in geriatric patients than in younger adult patients. Bacteria frequently develop resistance to antibiotics [2]. Drug-resistant bacteria as well as adverse drug reactions are more likely occurrences in the elderly, which means that drug therapy in this population fails or results in death at a higher rate [3].

Vancomycin, a bacteriostatic glycopeptide antibiotic, remains the most frequently used antibiotic in critically ill patients for the treatment of bacterial infections due to methicillin-resistant Staphylococcus aureus. One of the common adverse reactions of vancomycin is acute kidney injury (AKI) [4, 5], which has a high incidence and a short incubation period (it may occur in 48 h). Presently, the conventional guidance is that therapeutic drug monitoring (TDM) should be carried out within 48 h after stability, but the warning signs of injury may be delayed past this timeframe. However, if the pharmacokinetics of a patient can be predicted by using a population pharmacokinetic (PopPK) model based on renal function, the dose may be adjusted prior to dispensing the medication so that it is safer and more effective, which may effectively reduce the incidence of vancomycin AKI. Vancomycin-induced acute kidney injury (VI-AKI) occurs more readily in geriatric patients with pulmonary infections, so there is a more urgent need in these patients for an accurate dosage regimen.

The renal function of the elderly gradually decreases with age, which has a significant impact on the pharmacokinetics of vancomycin. A larger volume of distribution and a longer half-life of vancomycin have been observed in elderly patients as compared to younger adult patients [6]. Geriatric patients with pulmonary infections usually have poor nutrition, are associated with hypoalbuminemia, and suffer from internal environmental disorders such as hypokalemia, hyponatremia and metabolic acidosis [4, 6]. Thus, the elderly and the young are distinguished in terms of their pharmacokinetic characteristics. However, the pharmacokinetics of vancomycin in the elderly population are less known [6].

Currently, there are limited PopPK models for vancomycin studies in geriatric patients. It is necessary to establish a PopPK model and individual treatment plans for elderly patients to improve the safety and effectiveness of vancomycin therapy in geriatric patients.

The objective of the present study was to estimate the PopPK of vancomycin in Chinese geriatric patients (age ≥ 65 years) with pulmonary infections.

Ethics Approval

This study was approved by the Clinical Research Ethics Committee of Peking University First Hospital on December 15, 2015. The approval number is 2015[998].

Methods

Patients and Data Collection

All patients of Peking University First Hospital between January 2012 and December 2016 who satisfied the screening criteria were included in the study. The clinical data of these patients were retrospectively collected.

The inclusion criteria were as follows: (1) age ≥ 65 years, (2) patient had received intravenous vancomycin infusion for more than 2 days during hospitalization due to pulmonary infection [hospital-acquired pneumonia (HAP) or community-acquired pneumonia (CAP)], and (3) patients had received TDM more than once.

The exclusion criteria were as follows: (1) patient had multiple organ failure, renal replacement therapy or low volume shock, or (2) patient’s clinical data were lacking.

For those patients who qualified for the study, we collected basic information, vancomycin treatment information, TDM information, laboratory test results and information about concomitant drugs. The basic information consisted of sex, age, height and weight; the vancomycin treatment information consisted of dosage, interval, infusion time, and the number of infusions; the TDM information consisted of monitoring time and results; the laboratory test results consisted of serum creatinine, alanine aminotransferase, aspartate aminotransferase, total protein, albumin, total bilirubin, and urea nitrogen; and concomitant drugs consisted of amphotericin B, amikacin, diuretics (furosemide, torsemide, hydrochlorothiazide, spironolactone), and vasoactive drugs (catecholamines, epinephrine, norepinephrine).

Blood Sampling and Vancomycin Assays

The vancomycin dosing regimen consisted of 500 mg (every 6 h, 8 h, 12 h, 24 h, 48 h) and 1000 mg (every 8 h, 12 h), for 1.5–2 h intravenous infusion. The TDM usually occurred 0.5–2 h before and after the fourth or fifth dose of vancomycin, which were regarded as steady-state concentrations. All the TDM data were collected. All serum concentrations of vancomycin were determined via the chemiluminescence microparticle immunoassay method with the ARCHITECT i1000 system (Abbott Laboratories, Chicago, IL, USA). The coefficient of variation was less than 10%, and the linear range was between 3 and 100 mg/L.

PopPK Analysis

The vancomycin PopPK model was established and simulated with nonlinear mixed effect model software NONMEM® (v.7.3.0; ICON Development Solution, Ellicott City, MD, USA), and the Rstudio environment (v.98.1103, https://www.rstudio.com/) was used for statistical analysis and plotting. A first-order conditional estimation method with an inter- and intra-subject variability interaction (FOCE-I) option was used to calculate pharmacokinetic parameters.

The serum concentrations in the study were trough concentrations in steady-state; therefore, we established a one-compartment model that with first-order elimination (ADVAN1 TRANS2) was used as a base model to describe the relationship between serum concentration and time. The inter-individual variability of the model was formulated with an exponential random-effects formula, and the formulae were as follows:

CLi is the relative clearance of patient i, θCL the population typical value, and η the individual variation. The results show normal distribution with a mean of 0 and a variance of ω2 (N (0, ω2)).

The residual unexplained variability was fitted with a mixture model. The formula used was as follows:where C is the observed concentration of ith patient, C the simulated concentration of the ith patient, ε1, and ε2 residual variabilities that assume a normal distribution in which the mean is 0 and the variance is (N (0, ) and (N (0, )), ε1 proportional residual variability, and ε2 the additive residual variability.

The covariates examined in this study were sex (SEX), AGE, weight (WT), serum creatinine (SCR), creatinine clearance (CLCR), alanine aminotransferase (ALT) and aspartate aminotransferase (AST), urea nitrogen (BUN), total protein (TP), albumin (ALB) total bilirubin (TBIL), amikacin, amphotericin B, diuretics, and vasoactive drugs. CLCR was calculated using the Cockcroft–Gault formula. First, we selected covariates that could significantly improve the fitness of the model. Then, we added covariates into the base model and subsequently observed a change in the minimum value of the objective function (OFV); the covariate significantly enhanced the fitness of the model (χ2, p < 0.05) once the OFV decreased by more than 3.84, determining which enhancement could be added into the whole model. The selected covariates were added into the base model in the order of the degree of OFV decrease, starting with the covariate that caused the highest degree.

A backward elimination method was applied to determine the final model, and the covariates were eliminated in the reverse order they were used to build the total quantity model. We observed a change in the minimum value of the OFV. If the OFV increased by more than 6.64, the covariate significantly improved the fitness of the model (χ2, p < 0.01), and persisted in the model. The influences of continuous covariates and categorical covariates were validated by formulae as follows:

Model Evaluation

The internal verification of the final model was carried out by goodness-of-fit plots and statistical methods. Scatter plots were drawn of population-predicted (PRED) and individual-predicted (IPRED) concentrations versus observed concentrations (DV), and TIME versus the conditional weighted residuals (CWRES) and weighted residuals (WRES). A Y = X scatter trend line was added to the IPRED-DV and the PRED-DV. The closer the slope was to 1, the better the fit of the final model. The NPDE method was performed with Rstudio (v.0.98.1103, https://www.rstudio.com/). NPDEs were calculated to evaluate the predictive properties of the model. A simulation database was built, and the NPDE parameters were calculated. The mean, variance, skewness and kurtosis of the results were tested with statistical methods to assess whether the distribution approximated the standard normal distribution. We performed statistical tests associated with NPDE that consisted of: (1) a t test to determine whether there was a significant difference between the mean value of NPDE and 0; (2) a Fisher test to determine whether there was a significant difference between NPDE variance and 1; and (3) a Shapiro–Wilks test to determine whether there was a significant difference between the distribution of NPDE and the normal distribution. Histograms of NPDEs, plots of NPDE versus TIME and NPDE versus predicted concentrations were generated and used as additional simulation-based diagnostics. A total of 1000 bootstrap datasets was generated to test the statistical verification, reliability and stability of the PopPK model. The observed dataset was resampled with replacement data to generate a new dataset with the same size and population characteristics as the original dataset. The robust rate, point estimates and 95% confidence interval of the model were evaluated. If the results satisfied the conditions of success ratio of 1000 bootstraps > 90% and each parameter < 95% confidence interval, the model was considered more stable.

Simulation

We used the simulation module of the NONMEM software to simulate the concentration of vancomycin 1000 times.

Results

The development of the model was based on a database of 125 observations from 70. There were 49 males and 21 females in the study. The mean age of patients was 78.3 (± 6.96) years. A total of 57 patients (81.4%) had HAP while 13 patients (18.6%) had CAP. The number of concentration data points derived from each patient ranged from 1 to 5 (mean = 1.79). The basic information about the patients is shown in Table 1. Normally distributed continuous variables were expressed as the mean ± standard deviation, and non-normally distributed continuous variables were presented as the median (interquartile range).

Table 1

Demographic and pathophysiological characteristics of study patients

Characteristic	Statistics results
Male/female n	49/21
AGE/(years) (mean ± SD)	78.3 ± 6.96
Height/cm (mean ± SD)	161 ± 10
WT/kg (mean ± SD)	60.7 ± 10.2
Daily dose/g/day (mean ± SD)	1.55 ± 0.770
Serum concentration/mg/L (mean ± SD)	17 ± 8.03
SCR/μmol/L (mean ± SD)	90.6 ± 31
CLCR/mL/min (mean ± SD)	56.3 ± 22.1
BUN/mmol/L (minimum–maximum)	10.5 (3.18–86)
TP/g/dL (mean ± SD)	59.8 ± 8.92
ALB/g/dL (mean ± SD)	29.3 ± 4.12
HAP n (%)	57 (81.4%)
CAP n (%)	13 (18.6%)
Respiratory failure n (%)	46 (65.7%)
Hypertension n (%)	38 (54.3%)

WT weight, SCR serum creatinine, CL creatinine clearance rate, BUN blood urea nitrogen, TP total protein, ALB albumin, HAP hospital-acquired pneumonia, CAP community-acquired pneumonia

The results showed that CLCR, WT and AGE significantly improved the initial model. Consequently, CLCR, WT, and AGE were added to build out the whole model. However, when we added AGE to the model, the OFV was < 6.64, so AGE was not used in the final model.The final model is shown in Eqs. (9) and (10). The modeling process is shown in Table 2. When the influences of continuous covariates were validated by Eq. (3), the OFV was satisfied better.

Table 2

Model construction process

Model ID	OFV	ΔOFV	P value	Modeling
1	617.755	–	–	Basic model
Covariates on CL
2	606.197	− 11.558	P < 0.05	Add CLCR on CL in model 1
3	608.973	− 8.782	P < 0.05	Add AGE on CL in model 1
4	609.196	− 8.559	P < 0.05	Add WT on CL in model 1
Forward inclusion
5	601.501	− 4.696	P < 0.05	Add AGE on CL in model 2
6	602.416	− 3.781	P > 0.05	Add WT on CL in model 2
Backward elimination
7	608.973	7.472 > 6.64	P < 0.01	Sub CLCR on CL in model 6
8	606.197	4.696 < 6.64	P > 0.01	Sub AGE on CL in model 6
				Model 2 as final model

WT weight, CL creatinine clearance rate, OFV objective function value

We drew a scatter plot of continuous covariates and CL fitted by a basic model, and a linear trend line was added to the scatter plot. The results are shown in Fig. 1. SCR and AGE were noted to be inversely related to CL. CLCR and WT were noted to be positively related to CL.

Fig. 1

The correlation analysis between the clearance rate (CL) value of the vancomycin base model and continuous covariates. a Weight (WT) versus CL, b AGE versus CL, c serum creatinine (SCR) versus CL, d creatinine clearance rate (CLCR) versus CL

We drew scatter plots of PRED and the IPRED versus the DV. We also drew scatter plots of TIME versus WRES or CWRES. The results show that the value of WRES ranges from − 3  to 3, and the value of CWRES ranges from − 3 to 3 (see Fig. 2; Table 3).

Fig. 2

Diagnostic goodness of fit plots of the final model. a Observed concentration (DV) versus population-predicted concentration (PRED), b DV versus individual predicted concentration (IPRED), c conditional weighted residuals (CWRES) versus time, d weighted residuals (WRES) versus time

Table 3

Comparison of the base and final models

Parameter	Final model		Base model
	Value	RSE%	Value	RSE%
θ1CL (L/h)	2.45	6.9%	2.47	7.94%
θ2Vd (L)	154	9.2%	142	17.0%
θ3CLCR (mL/min) on CL	0.542	35.1%	–	–
ω CL	0.174	21.2%	0.235	40.6%
ω V	0.339	37.8%	0.213	46.8%
σ 1	0.0657	34.2%	0.0722	48.8%
σ 2	0 FIX	–	0 FIX	–

θ population parameter of CL, θ population parameter of Vd, θ population parameter of CLCR, ω intra-individual variation of CL, σ proportional residual variation, σ additive residual variation, CL creatinine clearance rate

An internal validation of the NPDE results showed that the mean and variance of the final model are 0.248 and 1.28, respectively; the values of the t test, Fisher test and Shapiro–Wilks test are 0.0158 (p > 0.001), 0.0388 (p > 0.01) and 0.861 (p > 0.05), respectively. The numerical results show that the distribution of the final model was close to normal distribution, and the model fit well (see Fig. 3). The y = x chart, histogram and random scatter plot showed a similar trend with a normal distribution. The results indicate a good fit, and the 1000 bootstrap statistical verification showed that the model robust rate is 92.1%, the parameter values are all within the 95% confidence interval, and the 95% confidence interval estimate did not consist of 0, indicating that the model was stable. The results are shown in Table 4.

Fig. 3

Normalized predictive distribution error (NPDE) of the population pharmacokinetics final model for vancomycin. a Quantile–quantile plot of NPDE versus normalized distribution, b the distribution chart of predictive distribution error, c NPDE versus time, d NPDE versus PRED

Table 4

Bootstrap results for the final model; 1000 iterations; success rate 92.1%

Parameter	NONMEM parameter	Bootstrap 95% CI	Bootstrap median	Bootstrap 95% CI
θ1CL(L/h)	2.45	(2.02, 2.88)	2.43	(2.09, 2.81)
θ2Vd (L)	154	(110, 197)	154	(117, 191)
θ3CLCR (mL/min) on CL	0.542	(0.141, 0.942)	0.538	(0.206, 0.878)
ω CL	0.174	(0.076, 0.272)	0.162	(0.092, 0.256)
ω V	0.339	(0.079, 0.598)	0.289	(0.121, 0.557)
σ 1	0.0657	(0.0177, 0.114)	0.0691	(0.0254, 0.106)
σ 2	0 FIX	–	0 FIX	–

CL creatinine clearance rate, CL clearance, V volume of distribution

The simulation results show that the initial drug regimen should be established according to the CLCR of patients. The results are shown in Table 5. For patients whose CLCR is > 50 mL/min, the initial vancomycin dosing regimen can be the full dose of 1000 mg administered every 8 h (18.41 ± 8.95 mg/L) or q 12 h (15.10 ± 7.02 mg/L). The concentration, however, may exceed 20 mg/L if the CLCR of the patient is not stable at a dosing regimen of 1000 mg, every 8 h. For patients whose CLCR is ≤ 50 mL/min, the initial dosing regimen can be modified to 1000 mg, q 12 h (17.63 ± 8.44 mg/L) since q 8 h may lead to a high concentration.

Table 5

Simulation data (simulation times = 1000)

Dosage regimen	Concentration (mg/L)
	Total (average plasma concentration)	CLCR > 50 mL/min	CLCR < 50 mL/min
1000 mg, q 8 h	19.26 ± 9.50	18.41 ± 8.95	20.72 ± 9.53
1000 mg, q 12 h	16.02 ± 7.51	15.10 ± 7.02	17.63 ± 8.44
1000 mg, q 24 h	10.29 ± 4.69	9.40 ± 4.38	11.96 ± 5.30
500 mg, q 6 h	11.22 ± 5.85	10.83 ± 4.57	11.88 ± 4.36
500 mg, q 8 h	9.82 ± 4.97	9.38 ± 3.65	10.57 ± 3.57
500 mg, q 12 h	7.98 ± 3.94	7.52 ± 2.68	8.79 ± 3.43
500 mg, q 24 h	5.06 ± 2.53	4.62 ± 1.35	5.89 ± 2.90

CL creatinine clearance rate, q × h every × hours

Discussion

Our study aims to estimate the PopPK model of vancomycin in Chinese geriatric patients (age ≥ 65 years) with pulmonary infections because there is no suitable PopPK model for this population. After retrospectively collecting 125 observations from 70 patients of Peking University First Hospital, we built the following PopPK model: CL(L/h) = 2.45 × (CLCR/56.28)0.542, (9) Vd(L) = 154 (10).

We performed an internal validation and simulation to assess the final model. The distribution of the final model was close to normal distribution, and the model had good fit.

Our study found that CLCR is the only covariate that affects the clearance parameter in geriatric patients with pulmonary infections. Overall, 90% of vancomycin is excreted though the kidney. Renal function is undoubtedly one of the important factors that affects the pharmacokinetics of vancomycin, and the CLCR is a substantial factor related to age, weight and creatinine. Additionally, many studies of PopPK models of vancomycin in adults have demonstrated that CLCR is one of the most important covariates that affect vancomycin pharmacokinetics. In these studies, the typical value of CL in the final model was 2.45 L/h. We placed the average values of CL from these studies into our final model and obtained the typical values of 3.58 L/h and 2.95 L/h [8, 9]. In the study that found a typical value of 3.58 L/h, the subjects were adult patients. However, research has shown that the CL values of geriatric patients are smaller than those of adult patients [3]. In the study that found a typical value of 2.95 L/h, the subjects were geriatric patients and the typical values were similar. The models from all the studies we reference are shown in Table 6.

Table 6

Population pharmacokinetic model of vancomycin in adult patients

Author		Population	Sample size	Age (years)	No. of model	PopPK model
Lin [10]		Adult Chinese patients	179	51.6 ± 16.9	1	CL (L/h) = 7.56 × (CLCR/104.71)0.886; Vd (L) = 101
Deng [7]		Adult Chinese patients	72	NA	1	CL (L/h) = 4.90 (CLCR ≥ 80 ml/min);CL (L/h) = 0.0654 × CLCR if CLCR < 80 ml/min.; Vd (L) = 47.76
Leu [11]		MRSA infectors	76	NA	1	CL (mL/min) = CLCR;Vd (L) = (0.17 × AGE) + (0.22 × ABW) + 15
Revilla [12]		ICU patients	191	61.1 ± 16.3	1	CL (mL/min/kg) = 0.67 × CLCR + AGE-0.24;Vd (L/kg) = 0.82 × 2.49 × AA = 0 or 1 if SCR ≤ or SCR > 1 mg/dl
Sanchez [6]		Adult patients	141	55 ± 14.58	2	CL (L/h) = 0.157 + 0.563 × CLCR;V1 (L) = 0.283 × TBW; V2 (L) = 32.2 × AGE/53.5; Q = 0.563 × TBW
Tanaka [13]		Adult patients	164	74 (17–95)	1	CL (L/h) = 0.875 × GFR; Vd (L) = 0.864 (L/kg)
Dolton [14]		Burn patients	70	34 (15–88)	2	CL (L/h) = 4.7 × (CLCR/6.53);
		Control patients		72 (38 ~ 95)		V1 (L) = 68.4 × WT/70-33.1 × BURN;V2 (L) = 73(L/h) × WT/70; Q (L/h) = 4.54
Yamamoto [15]		Healthy people	106	21.7 (20 ~ 25)	2	CL (L/h) = 3.83, if CLCR ≤ 85 mL/minCL (L/h) = 0.0322 × CLCR + 0.32, if CLCR < 85 mL/min; Vd (L) = 0.206 × WT
Llopis-Salvia [16]		Adult patients	50	60 (18–81)	2	CL(L/h) = 0.034 × CLCR + 0.015 × TBW;V1 (L/h) = 0.414 × TBW; V2 (L/kg) = 1.32 × TBW; Q (L/h) = 7.48
Staatz [17]		Adult patients	102	66 (17–7)	1	CL(L/h) = 2.97 × (1 + 0.0205 × (CLCR − CLCRmedian));Vd (L/h) = 1.24
Buelga [18]		Blood cancer patients	215	64.7 ± 11.3	1	CL (L/h) = 1.08 × CLCR; Vd (L) = 0.98 × WT(L/kg)
Yasuhara [19]		Adult patients	190	64.3 (19.3–89.6)	2	CL (L/h) = 0.0478 × CLCR (CLce ≤ 85 mL/min)CL (L/h) = 3.51 × CLCR (CLCR > 85 mL/min);V1 = 60.7; V2 = 24.63; Q = 60.7
Lin Zhong [20]		Adult infected patients	98	52.9 ± 17.2	1	CL (L/h) = 0.18 × (Scr)0.19 × (Sex)0.31;V1 (L) = = 20.85 × (Scr)0.52 × (Age)0.23 × (WT)0.98
Liu Chang [21]		Adult patients	276	40.83 ± 21.23	1	CL (L/h) = (3.35 × Infection + (1-Infection) × 6.74) × e0.0829;Vd (L) = (42.8 × Infection + (1-infection) × 26.6) × e0.147Infected population Infection = 1, no-infected population Infection = 0
Wang Rong [22]		Severe patients in neurosurgery	42	53 ± 14	2	CL (L/h) = 3.55 × (CLCR/94.1)0.84;V1 (L) = (62.8 − ALB) × (AGE/57.9)0.37;V2 (L) = 35.4 × (1 + (0.81 × GEND) × e(AGE−57.9) × 0.03; Q = 7.04
Zhang Jin [8]		Elderly patients	64	73.99 ± 6.90	1	CL (L/h) = 4.23 × (CLCR/76.82)0.774; V (L) = 71.4
Chen Bing [23]		Adult patients	169	47.3 ± 19.6	2	CL (L/h) = 4.02 × (CLCR/90)0.58; V1 (L) = 21.8; V2 (L) = 67.1; Q (L/h) = 6.01
He Xiaorong [24]		Elderly patients	260	76.79 ± 11.18	1	CL (L/h) = (1.71 + 8.31 × (1 − e(−0.0113(L/h)×CLCR))) × 0.475(AGE/72)Vd (L) = 26.2 × 0.475[AGE/72]
Weng Fangjuan [25]		Adult Gram-positive bacterial infection in ICU	103	NA	1	CL (L/h) = 0.044 × CLCR; V1 (L) = 0.542 × AGE; V2 (L) = 44.2;Q (L/h) = 6.95
Wu Wei [26]		Adult patients	100	51.75 ± 16.54	1	CL (L/h) = 7.56 × (CLCR/104.7133)0.886; Vd (L) = 101
Meng Long [27]		NA	129	NA	2	CL (L/h) = 3.56 × (CLCR/94.1)0.85;V1 (L) = (31 − (ALB-33.1) × 108) × (AGE/57.9)0.34;V2 (L) = 21.6 × (1 + 1.27 × SEX); Q (L/h) = 5.94
Wang Yang [28]		Severe lower respiratory tract infection patients	70	NA	2	CL (L/h) = 1.67 × e(η1),V1 (L) = 33.04 × [1 − 0.199 × (60/SCR)] × e(η2)V2 (L) = 19.29 × e(η4)Q (L/h) = 7.08 × 2.053(ALB) × e(η3),

CL creatinine clearance rate, Vd volume of distribution, V1 volume of distribution of the central compartment; V2 volume of distribution of the peripheral compartment, Q intercompartmental clearance, GEND gender, WT bodyweight, SCR serum creatinine

The simulation results show that the initial dosing regimen should be made according to the CLCR of patients. For patients whose CLCR is > 50 mL/min, the initial regimen can be 1000 mg, q 8 h (18.41 ± 8.95 mg/L) or q 12 h (15.10 ± 7.02 mg/L). However, the concentration may exceed 20 mg/L with dosing of 1000 mg, q 8 h if the CLCR of the patient is not stable. For patients whose CLCR is ≤ 50 mL/min, the initial dosing regimen can be 1000 mg, q 12 h (17.63 ± 8.44 mg/L) since q 8 h may lead to a high concentration.

There are some limitations to our study. The concentration values reflect only steady-state trough concentrations since the concentration data were obtained from routine TDM. Furthermore, our study was a retrospective study, which may affect the accuracy of some of the data and influence the final model. Further research, including a prospective study, should be carried out to build a more precise model.

Conclusion

A one-compartment PopPK model of vancomycin in Chinese geriatric patients with pulmonary infections was established. CLCR was found to be the only covariate that affected vancomycin pharmacokinetic parameters. The predictive performance of the PopPK model in geriatric patients differed significantly from that of the PopPK model in adult patients.

A population pharmacokinetic model was established to estimate the pharmacokinetics characteristics in Chinese geriatric patients with pulmonary infections.

The model could be used to develop an initial assessment of vancomycin dosing in geriatric patients.