Next Article in Journal
Staphylococcal Infections and Neonatal Skin: Data from Literature and Suggestions for the Clinical Management from Four Challenging Patients
Next Article in Special Issue
All Roads Lead to Rome: Enhancing the Probability of Target Attainment with Different Pharmacokinetic/Pharmacodynamic Modelling Approaches
Previous Article in Journal
Screening for Tuberculosis Infection among Migrants: A Cost-Effectiveness Analysis in the Italian Context
Previous Article in Special Issue
Pulmonary Pharmacokinetic and Pharmacodynamic Evaluation of Ampicillin/Sulbactam Regimens for Pneumonia Caused by Various Bacteria, including Acinetobacter baumannii
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting the Area under the Plasma Concentration-Time Curve (AUC) for First Dose Vancomycin Using First-Order Pharmacokinetic Equations

by
Kritsaporn Sujjavorakul
1,2,
Wasan Katip
3,4,
Stephen J. Kerr
5,6,7,
Noppadol Wacharachaisurapol
8,9 and
Thanyawee Puthanakit
1,9,*
1
Department of Pediatrics, Faculty of Medicine, Chulalongkorn University, Bangkok 10330, Thailand
2
Critical Care Excellence Center, King Chulalongkorn Memorial Hospital, Bangkok 10330, Thailand
3
Department of Pharmaceutical Care, Faculty of Pharmacy, Chiang Mai University, Chiang Mai 50200, Thailand
4
Epidemiology Research Group of Infectious Disease (ERGID), Chiang Mai University, Chiang Mai 50200, Thailand
5
Biostatistics Centre, Faculty of Medicine, Chulalongkorn University, Bangkok 10330, Thailand
6
The Kirby Institute, The University of New South Wales, Kensington, NSW 2052, Australia
7
HIV-NAT, The Thai Red Cross AIDS Research Centre, Bangkok 10330, Thailand
8
Center of Excellence in Clinical Pharmacokinetics and Pharmacogenomics, Department of Pharmacology, Faculty of Medicine, Chulalongkorn University, Bangkok 10330, Thailand
9
Center of Excellence for Pediatric Infectious Diseases and Vaccines, Faculty of Medicine, Chulalongkorn University, Bangkok 10330, Thailand
*
Author to whom correspondence should be addressed.
Antibiotics 2023, 12(4), 630; https://doi.org/10.3390/antibiotics12040630
Submission received: 20 February 2023 / Revised: 15 March 2023 / Accepted: 21 March 2023 / Published: 23 March 2023

Abstract

:
To treat critically ill patients, early achievement of the target area under the plasma concentration-time curve/minimum inhibitory concentration (AUC/MIC) in the first 24 h is recommended. However, accurately calculating the AUC before steady state is an obstacle to this goal. A first-order pharmacokinetic equation to calculate vancomycin AUC after a first dose of vancomycin has never been studied. We sought to estimate AUC using two first-order pharmacokinetic equations, with different paired concentration time points, and to compare these to the actual first dose vancomycin AUC calculated by the linear-log trapezoid rule as a reference. The equations were validated using two independent intensive first dose vancomycin concentration time data sets, one from 10 adults and another from 14 children with severe infection. The equation with compensation for the alpha distribution phase using a first vancomycin serum concentration from 60 to 90 min and the second concentration from 240 to 300 min after the completed infusion showed good agreement and low bias of calculated AUC, with mean differences <5% and Lin’s correlation coefficient >0.96. Moreover, it gave an excellent correlation with Pearson’s r > 0.96. Estimating the first dose vancomycin AUC calculated using this first-order pharmacokinetic equation is both reliable and reproducible in clinical practice settings.

1. Introduction

Vancomycin is a commonly used broad-spectrum antibiotic which covers multidrug-resistant Gram-positive bacteria, e.g., methicillin-resistant staphylococci, penicillin-resistant Streptococcus pneumoniae (PRSP), and ampicillin-resistant enterococci [1]. The current guidelines on therapeutic drug monitoring (TDM) for vancomycin recommend achieving the target area under the plasma concentration-time curve (AUC) within the first 24–48 h. This approach maximizes the therapeutic effect, while simultaneously minimizing adverse effects, especially vancomycin-associated acute kidney injury (vAKI) [2]. Early TDM and early dose adjustment are especially important in critically ill patients whose mortality rates are high during the initial phase of illness. For example, between one-third and half of deaths in patients with severe sepsis and septic shock occur within the first 48 h after diagnosis [3,4,5]. However, in clinical practice, routine TDM for vancomycin is measured at steady state at least 24 h after the initial dose. This limits the opportunity to optimize the AUC during critical periods of illness.
The largest pooled population-pharmacokinetic modeling study characterized vancomycin pharmacokinetics as a two-compartment distribution model [6]. There is no apparent metabolism [7], and the absolute elimination rate is a linear function of its concentration in plasma following first-order kinetics [8].
Three methods are commonly used to estimate the AUC: the linear-trapezoid rule, Bayesian-derived AUC monitoring, and first-order pharmacokinetic equations. An accurate collection of multiple concentrations over the same dosing interval is required for the linear-trapezoid rule, making this method unsuitable outside the research setting. Bayesian-derived AUC monitoring is recommended in vancomycin TDM. However, a well-developed population-specific pharmacokinetic model of vancomycin as the Bayesian prior is required, and thus Bayesian dose-optimizing software is not available in some special populations, for example obese patients, critically ill patients, pediatric patients, and patients with unstable renal function [2]. The first-order pharmacokinetic equation method makes fewer assumptions than the Bayesian approach and is simple enough for routine clinical practice use. The advantage of this method is that it provides a snapshot AUC for the sampling period, which is beneficial in groups of patients with high variability in drug pharmacokinetics [2].
First-order pharmacokinetic equations for vancomycin TDM at steady state have been validated [9,10] and are recommended in guidelines [2]. However, the use of this method for vancomycin TDM after the initial dose has not been studied. The aim of this study was to evaluate the usefulness of first-order pharmacokinetic equations for predicting the first-dose vancomycin AUC, and to identify the optimal sampling times to be used in the calculation. We hypothesized that an equation using two appropriate concentration time points could approximate the AUC of the initial vancomycin dose calculated by intensive sampling with reasonable precision, low bias, and high correlation.

2. Results

We used two previously published full pharmacokinetic data sets to evaluate the pharmacokinetics of initial dose vancomycin. These were an adult data set from 10 septic shock patients [11], and a pediatric data set from 14 children with severe infection [12]. In the adult data set, the median age was 59 (IQR 45.8–79.5) years and median creatinine clearance was 43.5 (IQR 28.8–92.5) mL/min. An initial vancomycin dose of 30 mg/kg was infused over 120 min. In the pediatric data set, the median age was 6.4 (IQR 3.3–10.7) years and the median creatinine clearance was 183.1 (IQR 148.2–219.5) mL/min/1.73 m2; the initial vancomycin dose of 15 mg/kg was infused over 60 min. Detailed characteristics of these studies are described in Table S1. Details about the building of pharmacokinetic models for the reference standard and two first-order pharmacokinetic equations are provided in the Materials and Methods section.
The analysis results from the adult data set using Model 1 are shown in Table S2. In this analysis, there were no equation-paired concentration time points with an acceptable mean difference from the Bland–Altman analysis. In contrast, applying Model 2 to the adult data set (Table S3), C1 sampled at 0, 40, 60, and 90 min after the completed infusion combined with C2 at 240 min after the completed infusion gave an acceptable −4.6 to 4.1 percent mean difference. The correlation and agreement were better using C1 obtained 60 and 90 min after the end of the infusion (Pearson’s r = 0.976 and Lin’s correlation coefficient = 0.971 and 0.967, respectively) than those obtained at 0 and 40 min (Pearson’s r = 0.899, 0.929 and Lin’s correlation coefficient = 0.869 and 0.923, respectively).
In the pediatric data set, Model 1 (Table S4) using C2 at 240 min after the completed infusion, only C1 sampled immediately after completing the infusion showed an acceptable −2.7 percent mean difference, with Pearson’s r = 0.974 and Lin’s correlation coefficient = 0.940. Model 1 using C2 sampled 300 min after the completed infusion and C1 at 0 and 15 min gave 3.4 and −4.8 percent mean difference with Pearson’s r = 0.977, 0.98 and Lin’s correlation coefficient = 0.965, 0.960, respectively.
Table S5 shows the results from the pediatric data set using Model 2. C1 between samples 60 and 180 min after the infusion ended resulted in a −4.2 to 2.8 percent mean difference when C2 was taken either at 240 or 300 min after the completed infusion. Alternatively, C1 sampled 30 min after the completed infusion gave an acceptable percent mean difference only with C2 at 240 min after the infusion. In terms of correlation and agreement, all mentioned paired concentration time points had Pearson’s r > 0.96 and Lin’s correlation coefficient >0.95, with the exception of C1 at 180 min paired with C2 at 240 min after the completed infusion and C1 at 180 min paired with C2 at 300 min after the completed infusion, both of which had lower correlation and agreement (Pearson’s r = 0.883 and Lin’s correlation coefficient = 0.881; Pearson’s r = 0.939 and Lin’s correlation coefficient = 0.933, respectively).
Selected equation-paired concentration time points with an acceptable mean difference <5% from the Bland–Altman analysis from both the adult and pediatric data sets are shown in Table 1.
Considering all results and looking for intersected equation-paired concentration time points that provided an acceptable percent mean difference, Model 2 consistently provided good agreement, low bias, and high correlation at C1 sampled in a time window from 60 to 90 min after the completed infusion, and C2 sampled in a time window from 240 to 300 min after the completed infusion (Figure 1). Bland–Altman plots for the selected equation-paired concentration time points are shown in Figure 2.
The pooled data available from Model 2 with C1 sampled from 60 to 90 min, and C2 from 240 to 300 min after the completed infusion, showed a strong correlation and goodness of fit between the reference standard and the proposed selected equation-paired concentration time point combinations. Figure 3 shows the regression line of best fit between the two methods. The regression line had a slope of 1.03 and an intercept of 0.15, with a correlation coefficient between the two methods of r = 0.994 (p < 0.001).
Abbreviations: AUCf, full area under the plasma concentration-time curve.
In addition to the full area under the plasma concentration time curve (AUCf) calculation, the first-order pharmacokinetic equations were also used for calculating other pharmacokinetic parameters including vancomycin clearance (Cl), volume of distribution (Vd), and half-life (t1/2). The corresponding mean and standard deviation for these parameters, as well as the reference standards, are shown in Table 2 and Table 3.

3. Discussion

We used two previously published intensive pharmacokinetic data sets and two first-order pharmacokinetic equations to assess which paired post-infusion concentration time point combinations could best predict the actual vancomycin AUCf. We found that Model 2 with a compensation area for the alpha-phase, vancomycin serum concentrations collected from 60 to 90 min after the completed infusion as C1, and those collected 240–300 min after the completed infusion as C2 produced a calculated AUCf with less than 5% mean difference compared to the full linear-log trapezoid method. Pearson’s correlation coefficient and Lin’s agreement coefficient for these combinations were >0.96. The mean difference and 95% limit of agreement for this equation-paired concentration time points combination were sufficiently narrow for clinical purposes. These findings were consistent with whether the infusion period was one or two hours, and regardless of the difference in kidney function in the two data sets. As shown in Table 2 and Table 3, other pharmacokinetic parameters, including vancomycin Cl, Vd, and t1/2, showed comparable results with the log-linear trapezoid method. These findings suggested that Model 2 is highly relevant for determining pharmacokinetic parameters after the first dose of vancomycin.
We focused on using two time points within 6 h after the start of the infusion for two reasons. First, this would allow for early dose adjustment, if necessary, in critically ill patients who need to reach the target AUC as soon as possible because achieving the appropriate target of AUC values can reduce mortality in critically ill patients [13]. Second, using two time points within the first six hours makes this approach applicable to both adult and pediatric patients, since the dosing interval in pediatric patients is usually shorter at 6–8 h.
A previous study found that first-order pharmacokinetic equations accurately predicted aminoglycoside AUC [14], daptomycin AUC [15], and steady state vancomycin AUC [9]. The consistency of our results with these other antimicrobial agents is likely due to the similarity in concentration time profile shape.
Pai and Rodvold reported the AUC for aminoglycosides calculation using a first order-pharmacokinetic equation analogous to Model 1 in our study, which had a percent bias of −2.80 (95% CI = −3.06 to −2.54) [14]. In addition, they found that the same method could estimate the AUC of daptomycin with a bias of no more than ±10% [15]. These levels of bias were not expected to be clinically significant; as a result, these methods are widely used for optimizing the administration of aminoglycoside and daptomycin in clinical practice.
Pai et al. also conducted an approach using first-order pharmacokinetic equations to calculate the steady-state vancomycin AUC [9]. The authors proposed two equations analogous to those of Model 1 and Model 2 in our study. In steady state both models showed good agreement with median error < 2%, but their Model 1 tended to underestimate the prediction of the AUC versus the reference standard, consistent with the results of our study.
The most likely reason that Model 1 underestimates AUC versus the reference standard is because this model fails to capture the alpha-phase of the concentration-time curve (Figure S1). Vancomycin is infused over one or two hours, and this equation uses concentration time values that occur after the infusion period; therefore, data from the alpha distribution phase are excluded, and this has a significant impact on the AUCf calculation. As shown in Figure S2, Model 2 captures additional area that at least partially compensates for the alpha-phase AUC, which is not captured by Model 1, resulting in lower bias and improved prediction.
The main advantages of using first-order pharmacokinetic equations are their simplicity and generalizability. Calculations are easy to make using basic computer calculation software, so the method can be easily applied in routine clinical practice. The first-order pharmacokinetic equation method relies primarily on accurate sampling times and concentration measurements, not on the patient’s clinical status, as the absolute rate of vancomycin elimination is a linear function, as previously mentioned [8]. Snapshot measurements with minimal assumptions are especially useful in patients with dynamic clinical states, for example septic patients. As demonstrated in this study, equation-paired concentration time point combinations could accurately predict pharmacokinetic parameters in both intensive pharmacokinetic data sets, despite the disparate demographic characteristics, kidney function, initial infusion dose, and infusion period.
There are a number of limitations to our research. First, although the target AUC for the first dose of vancomycin is unclear, DeRyke and Alexander [16] demonstrate that pharmacokinetic parameters derived from pharmacokinetic equations can be used to derive an area under the plasma concentration-time curve in the first 24 h of treatment (AUC0–24) for each individual patient. Calculating pharmacokinetic parameters based on the first vancomycin dose can be used for optimizing the subsequent doses using the method described by Sawchuk-Zaske [17]. Several studies using this method to individualize vancomycin regimens showed that first-dose TDM resulted in faster target AUC attainment than routine steady-state TDM in various group of patients, including critically ill adults and neonates [18,19,20]. Second, our study aimed to identify the best-performing first-order pharmacokinetic equation and optimal sampling times that reliably reproducibly estimated first dose vancomycin AUC and other pharmacokinetic parameters using two independent data sets. The data sets we used to compare the actual AUC against the AUC predicted by each equation were from a group of adults with septic shock and a group of children with severe infection (Table S1). Our findings may therefore need further evaluation before generalization to other groups of patients, such as those who are overweight/obese, and those with renal impairment. Lastly, the data sets used in our study were comparatively small, and resulted in wide standard deviations around the mean difference in Bland-Altman analysis. Although conducting a new similar study with a larger sample size would likely lower this variability, accessing intensive pharmacokinetic data sets to make these calculations is difficult and designing a new study is an ethically challenging because the subjects are vulnerable. Furthermore, the 95% limit of agreements of selected equation-paired concentration time point combinations in our study was not expected to be clinically significant. Therefore, despite the limitations mentioned and the small sample size, our current study results provide information which can guide clinical practice in some groups of patients, and pharmacokinetic parameters based on the first vancomycin dose are useful for further AUC0–24 research study (TCTR20210617001).

4. Materials and Methods

The study flow diagram of this study was presented in Figure 4.

4.1. Pharmacokinetic Model Building

4.1.1. The Reference Standard

The reference AUCf of the individuals in both studies were calculated using all plasma concentrations available (10 time points for subjects in the adult data set and 9 time points for subjects in the pediatric data set, Table S1) and the linear-log trapezoid rule. The pharmacokinetic analysis of the serum concentration-time data was performed using PKSolver (version 2.0, China Pharmaceutical University, Nanjing, China) [21].

4.1.2. The First-Order Pharmacokinetic Equations [9]

Using two available serum vancomycin concentrations, the elimination rate constant (Ke) can be calculated using the following equation [22]:
Ke = Ln C 1 C 2 t ,
where C1 is the first concentration measured after the infusion has been completed, C2 is the second concentration collected toward the end of the dosing interval, and t is the difference in time between C1 and C2. Once the Ke is computed, it can be used to compute theoretical concentrations through forward- and backward-extrapolation.
In the adult data set, vancomycin was infused over 120 min and concentration levels were sampled at 0, 10, 20, 40, 60, 90, 120, and 240 min after the completed infusion. C1 was sequentially chosen from 0 to 120 min after the completed infusion and C2 was the last sample, collected 240 min after the completed infusion. In the pediatric data set, vancomycin was infused over 60 min and levels were sampled at 0, 15, 30, 60, 120, 180, 240, and 300 min after the completed infusion. C1 was sequentially chosen from 0 to 240 min after the completed infusion and C2 was chosen from 240 or 300 min after the completed infusion.
In this study, we proposed and compared two first-order pharmacokinetic equation models analogous to those used by Pai et al. in their study on optimizing vancomycin delivery at steady-state [9]. The models can be described as follows:
For the first dose, the vancomycin concentration at the start of infusion is zero. The end of infusion concentration is estimated (Ceoi′) regardless effect of the alpha-phase of the 2-compartment distribution model. The area between the start and the end of the infusion time (teoi) can be measured as the area of the triangle:
AUC 0 t eoi = t eoi × C eoi 2
The area under the exponential part from teoi to the end of dose (t) is:
AUC t eoi t = C eoi Ke
Therefore, the AUCf for first dose vancomycin in Model 1 can be simplified to:
AUC f = t eoi × C eoi 2 + C eoi Ke
This model captures additional area by backward extrapolating peak concentration to the start of infusion (C0′), aiming to compensate for the unmeasured alpha-phase. Under this model, the equation can be simplified to:
AUC f = C 0 Ke
First-order pharmacokinetic equation calculations were performed using Microsoft Office Excel 365 (Redmond, Washington, DC, USA).

4.2. Statistical Analysis

To evaluate Equations (4) and (5), the linear-log trapezoid rule was used to calculate the AUCf using all available samples from the full pharmacokinetic studies, and served as the reference standard. Bland-Altman analysis was used to assess agreement and bias [23] between the reference standard and each combination of the first-order pharmacokinetic equations and selected paired concentration time points. To facilitate the comparison of the equation-derived and gold standard AUC from the adult and pediatric studies, results were expressed as percent change. An a priori acceptable mean difference was set at 5%. Pearson’s correlation and linear regression were used to assess the linear correlation of estimates. Lin’s correlation coefficient was also used to assess agreement. Statistical significance was defined as p < 0.05. Furthermore, a summary of pharmacokinetic parameters derived from selected equation-paired concentration time points combination including vancomycin Cl, Vd, and t1/2 were calculated and compared with the parameters derived from the linear-log trapezoid rule, which served as the reference standard. All statistical analyses were performed using Stata version 15.1 (Stata Corp LCC, College Station, TX, USA).

5. Conclusions

Of the two first-order pharmacokinetic equations tested, the one that used plasma samples obtained from 60 to 90 and 240 to 300 min after the completed infusion reliably and reproducibly estimated first dose vancomycin AUC. If physicians used this equation in clinical practice settings, doses could be adjusted earlier than in conventional practice and lead to a more rapid achievement of the target AUC. However, whether achieving the target AUC earlier would definitively improve patient outcomes by enhancing the therapeutic effect and reducing adverse drug effects needs further investigation.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antibiotics12040630/s1, Table S1: Data characteristics of adult and pediatric data sets; Table S2: Agreement, bias, and correlation results from the adult data set using Model 1; Table S3: Agreement, bias, and correlation results from the adult data set using Model 2; Table S4: Agreement, bias, and correlation results from the pediatric data set using Model 1; Table S5: Agreement, bias, and correlation results from the pediatric data set using Model 2; Figure S1: Expected area under the concentration curve calculated using Model 1 compared to the expected first dose vancomycin concentration time profile; Figure S2: Expected area under the concentration curve calculated using Model 2 compared to the expected first dose vancomycin concentration time profile.

Author Contributions

Conceptualization, K.S. and T.P.; Data curation, K.S. and W.K.; Formal analysis, K.S. and S.J.K.; Funding acquisition, K.S. and T.P.; Investigation, K.S. and W.K.; Methodology, K.S., W.K., S.J.K., N.W. and T.P.; Project administration, K.S. and T.P.; Resources, W.K. and N.W.; Software, K.S., W.K., S.J.K. and N.W.; Supervision, T.P.; Validation, K.S., W.K., S.J.K., N.W. and T.P.; Visualization, K.S., W.K., S.J.K., N.W. and T.P.; Writing—original draft, K.S.; Writing—review and editing, K.S., W.K., S.J.K., N.W. and T.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research project was supported by the Second Century Fund (C2F), Chulalongkorn University.

Institutional Review Board Statement

Ethical review and approval were waived by the Institutional Review Board, Faculty of Medicine, Chulalongkorn University for this study because this manuscript is based on already collected and published data.

Informed Consent Statement

Patient consent was waived because this manuscript is based on already collected and published data.

Data Availability Statement

No new data were created in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors thank Satja Issaranggoon Na Ayuthaya and Peninnah Oberdorfer, the principal investigators of the pediatric data study, for kindly providing the original data set used for this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Weiss, S.L.; Peters, M.J.; Alhazzani, W.; Agus, M.S.D.; Flori, H.R.; Inwald, D.P.; Nadel, S.; Schlapbach, L.J.; Tasker, R.C.; Argent, A.C.; et al. Surviving Sepsis Campaign International Guidelines for the Management of Septic Shock and Sepsis-Associated Organ Dysfunction in Children. Pediatr. Crit. Care Med. 2020, 21, e52–e106. [Google Scholar] [CrossRef] [PubMed]
  2. Rybak, M.J.; Le, J.; Lodise, T.P.; Levine, D.P.; Bradley, J.S.; Liu, C.; Mueller, B.A.; Pai, M.P.; Wong-Beringer, A.; Rotschafer, J.C.; et al. Therapeutic monitoring of vancomycin for serious methicillin-resistant Staphylococcus aureus infections: A revised consensus guideline and review by the American Society of Health-System Pharmacists, the Infectious Diseases Society of America, the Pediatric Infectious Diseases Society, and the Society of Infectious Diseases Pharmacists. Am. J. Health Syst. Pharm. 2020, 77, 835–864. [Google Scholar] [PubMed] [Green Version]
  3. Weiss, S.L.; Balamuth, F.; Hensley, J.; Fitzgerald, J.C.; Bush, J.; Nadkarni, V.M.; Thomas, N.J.; Hall, M.; Muszynski, J. The Epidemiology of Hospital Death Following Pediatric Severe Sepsis: When, Why, and How Children with Sepsis Die. Pediatr. Crit. Care Med. 2017, 18, 823–830. [Google Scholar] [CrossRef] [PubMed]
  4. Cvetkovic, M.; Lutman, D.; Ramnarayan, P.; Pathan, N.; Inwald, D.P.; Peters, M.J. Timing of death in children referred for intensive care with severe sepsis: Implications for interventional studies. Pediatr. Crit. Care Med. 2015, 16, 410–417. [Google Scholar] [CrossRef] [PubMed]
  5. Law, A.C.; Stevens, J.P.; Walkey, A.J. National Trends in Timing of Death among Patients with Septic Shock, 1994–2014. Crit. Care Med. 2019, 47, 1493–1496. [Google Scholar] [CrossRef]
  6. Colin, P.J.; Allegaert, K.; Thomson, A.H.; Touw, D.J.; Dolton, M.; de Hoog, M.; Roberts, J.A.; Adane, E.D.; Yamamoto, M.; Santos-Buelga, D.; et al. Vancomycin Pharmacokinetics Throughout Life: Results from a Pooled Population Analysis and Evaluation of Current Dosing Recommendations. Clin. Pharm. 2019, 58, 767–780. [Google Scholar] [CrossRef] [Green Version]
  7. National Center for Biotechnology Information. PubChem Compound Summary for CID 14969, Vancomycin. Available online: https://pubchem.ncbi.nlm.nih.gov/compound/Vancomycin (accessed on 17 May 2020).
  8. Buxton, I.L.O. Pharmacokinetics: The Dynamics of Drug Absorption, Distribution, Metabolism, and Elimination. In Goodman & Gilman’s: The Pharmacological Basis of Therapeutics, 13e; Brunton, L.L., Hilal-Dandan, R., Knollmann, B.C., Eds.; McGraw-Hill Education: New York, NY, USA, 2017. [Google Scholar]
  9. Pai, M.P.; Neely, M.; Rodvold, K.A.; Lodise, T.P. Innovative approaches to optimizing the delivery of vancomycin in individual patients. Adv. Drug Deliv. Rev. 2014, 77, 50–57. [Google Scholar] [CrossRef] [Green Version]
  10. Suchartlikitwong, P.; Anugulruengkitt, S.; Wacharachaisurapol, N.; Jantarabenjakul, W.; Sophonphan, J.; Theerawit, T.; Chatsuwan, T.; Wattanavijitkul, T.; Puthanakit, T. Optimizing Vancomycin Use Through 2-Point AUC-Based Therapeutic Drug Monitoring in Pediatric Patients. J. Clin. Pharm. 2019, 59, 1597–1605. [Google Scholar] [CrossRef]
  11. Katip, W.; Jaruratanasirikul, S.; Pattharachayakul, S.; Wongpoowarak, W.; Jitsurong, A.; Lucksiri, A. The pharmacokinetics of vancomycin during the initial loading dose in patients with septic shock. Infect. Drug Resist. 2016, 9, 253–260. [Google Scholar] [CrossRef] [Green Version]
  12. Issaranggoon Na Ayuthaya, S.; Katip, W.; Oberdorfer, P.; Lucksiri, A. Correlation of the vancomycin 24-h area under the concentration-time curve (AUC24) and trough serum concentration in children with severe infection: A clinical pharmacokinetic study. Int. J. Infect. Dis. 2020, 92, 151–159. [Google Scholar] [CrossRef] [Green Version]
  13. Zelenitsky, S.; Rubinstein, E.; Ariano, R.; Iacovides, H.; Dodek, P.; Mirzanejad, Y.; Kumar, A.; Cooperative Antimicrobial Therapy of Septic Shock-CATSS Database Research Group. Vancomycin pharmacodynamics and survival in patients with methicillin-resistant Staphylococcus aureus-associated septic shock. Int. J. Antimicrob. Agents 2013, 41, 255–260. [Google Scholar] [CrossRef]
  14. Pai, M.P.; Rodvold, K.A. Aminoglycoside dosing in patients by kidney function and area under the curve: The Sawchuk-Zaske dosing method revisited in the era of obesity. Diagn. Microbiol. Infect. Dis. 2014, 78, 178–187. [Google Scholar] [CrossRef]
  15. Pai, M.P.; Russo, A.; Novelli, A.; Venditti, M.; Falcone, M. Simplified equations using two concentrations to calculate area under the curve for antimicrobials with concentration-dependent pharmacodynamics: Daptomycin as a motivating example. Antimicrob. Agents Chemother. 2014, 58, 3162–3167. [Google Scholar] [CrossRef] [Green Version]
  16. Deryke, C.A.; Alexander, D.P. Optimizing Vancomycin Dosing through Pharmacodynamic Assessment Targeting Area under the Concentration-Time Curve/Minimum Inhibitory Concentration. Hosp. Pharm. 2009, 44, 751–765. [Google Scholar] [CrossRef]
  17. Sawchuk, R.J.; Zaske, D.E.; Cipolle, R.J.; Wargin, W.A.; Strate, R.G. Kinetic model for gentamicin dosing with the use of individual patient parameters. Clin. Pharm. 1977, 21, 362–369. [Google Scholar] [CrossRef]
  18. Shahrami, B.; Najmeddin, F.; Mousavi, S.; Ahmadi, A.; Rouini, M.R.; Sadeghi, K.; Mojtahedzadeh, M. Achievement of Vancomycin Therapeutic Goals in Critically Ill Patients: Early Individualization May Be Beneficial. Crit. Care Res. Pract. 2016, 2016, 1245815. [Google Scholar] [CrossRef] [Green Version]
  19. Truong, J.; Smith, S.R.; Veillette, J.J.; Forland, S.C. Individualized Pharmacokinetic Dosing of Vancomycin Reduces Time to Therapeutic Trough Concentrations in Critically Ill Patients. J. Clin. Pharm. 2018, 58, 1123–1130. [Google Scholar] [CrossRef]
  20. Crumby, T.; Rinehart, E.; Carby, M.C.; Kuhl, D.; Talati, A.J. Pharmacokinetic comparison of nomogram-based and individualized vancomycin regimens in neonates. Am. J. Health Syst. Pharm. 2009, 66, 149–153. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Huo, M.; Zhou, J.; Xie, S. PKSolver: An add-in program for pharmacokinetic and pharmacodynamic data analysis in Microsoft Excel. Comput. Methods Programs Biomed. 2010, 99, 306–314. [Google Scholar] [CrossRef]
  22. Sawchuk, R.J.; Zaske, D.E. Pharmacokinetics of dosing regimens which utilize multiple intravenous infusions: Gentamicin in burn patients. J. Pharm. Biopharm. 1976, 4, 183–195. [Google Scholar] [CrossRef]
  23. Giavarina, D. Understanding Bland Altman analysis. Biochem. Med. 2015, 25, 141–151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Summary of intersected C1 and C2 windows as inputs of the first-order pharmacokinetic equation using Model 2.
Figure 1. Summary of intersected C1 and C2 windows as inputs of the first-order pharmacokinetic equation using Model 2.
Antibiotics 12 00630 g001
Figure 2. Bland-Altman plots of the reference standard and each selected equation-paired concentration time points combination. (A) Adult data set using Model 2 with C1 at 60 min and C2 at 240 after the completed infusion had mean bias of 13.0 (3.0%) and 95% limits of agreement of −45.0 to 70.9. (B) Adult data set using Model 2 with C1 at 90 min and C2 at 240 after the completed infusion had mean bias of 18.0 (4.1%) and 95% limits of agreement of −40.8 to 76.7. (C) Pediatric data set using Model 2 with C1 at 60 min and C2 at 240 after the completed infusion had bias of −0.1 (−0.1%) and 95% limits of agreement of −22.6 to 22.3. (D) Pediatric data set using Model 2 with C1 at 60 min and C2 at 300 after the completed infusion had bias of 2.1 (1.9%) and 95% limits of agreement of −14.4 to 18.7.
Figure 2. Bland-Altman plots of the reference standard and each selected equation-paired concentration time points combination. (A) Adult data set using Model 2 with C1 at 60 min and C2 at 240 after the completed infusion had mean bias of 13.0 (3.0%) and 95% limits of agreement of −45.0 to 70.9. (B) Adult data set using Model 2 with C1 at 90 min and C2 at 240 after the completed infusion had mean bias of 18.0 (4.1%) and 95% limits of agreement of −40.8 to 76.7. (C) Pediatric data set using Model 2 with C1 at 60 min and C2 at 240 after the completed infusion had bias of −0.1 (−0.1%) and 95% limits of agreement of −22.6 to 22.3. (D) Pediatric data set using Model 2 with C1 at 60 min and C2 at 300 after the completed infusion had bias of 2.1 (1.9%) and 95% limits of agreement of −14.4 to 18.7.
Antibiotics 12 00630 g002
Figure 3. The regression line between the reference standard AUCf and pooled selected equation-paired concentration time point combinations. The regression equation was y = 0.15 + 1.03x (R2 = 0.988, p < 0.001). Pearson’s correlation coefficient between the two methods was r = 0.994 (p < 0.001).
Figure 3. The regression line between the reference standard AUCf and pooled selected equation-paired concentration time point combinations. The regression equation was y = 0.15 + 1.03x (R2 = 0.988, p < 0.001). Pearson’s correlation coefficient between the two methods was r = 0.994 (p < 0.001).
Antibiotics 12 00630 g003
Figure 4. Study flow diagram. Abbreviations: AUCf, full area under the plasma concentration-time curve.
Figure 4. Study flow diagram. Abbreviations: AUCf, full area under the plasma concentration-time curve.
Antibiotics 12 00630 g004
Table 1. Selected equation-paired concentration time points with an acceptable mean difference from the Bland–Altman analysis from adult and pediatric data sets.
Table 1. Selected equation-paired concentration time points with an acceptable mean difference from the Bland–Altman analysis from adult and pediatric data sets.
Data SetTime PointsBland–Altman AnalysisCorrelationLin’s
Coefficients
C1C2MeanDifference95% Limits of AgreementPearson’s rp-ValueRho_cp-Value
Mean (%)SD
Model 1
Pediatric0240111.9−3.0 (−2.7)13.4(−29.1, 23.2)0.974<0.0010.940<0.001
0300115.33.9 (3.4)10.2(−16.2, 24.0)0.977<0.0010.965<0.001
15300109.3−5.3 (−4.8)9.9(−24.7, 14.2)0.98<0.0010.960<0.001
Model 2
Adult0240417.5−19.3 (−4.6)60.2(−137.3, 98.8)0.899<0.0010.869<0.001
40240422.2−9.9 (−2.3)50.0(−107.9, 88.0)0.929<0.0010.923<0.001
60240433.613.0 (3.0)30.0(−45.0, 70.9)0.976<0.0010.971<0.001
90240436.118.0 (4.1)30.0(−40.8, 76.7)0.976<0.0010.967<0.001
Pediatric30240114.14.3 (3.8)11.3(−17.8, 26.4)0.969<0.0010.955<0.001
60240111.8−0.1 (−0.1)11.4(−22.6, 22.3)0.964<0.0010.962<0.001
120240109.6−4.6 (−4.2)11.6(−27.4, 18.2)0.963<0.0010.957<0.001
180240113.32.8 (2.5)20.5(−37.4, 43.0)0.883<0.0010.881<0.001
60300113.02.1 (1.9)8.4(−14.4, 18.7)0.981<0.0010.979<0.001
120300109.7−4.4 (−4.01)10.6(−25.1, 16.2)0.974<0.0010.966<0.001
180300113.53.2 (2.8)16.2(−28.5, 34.9)0.939<0.0010.933<0.001
Table 2. Summary of pharmacokinetic parameters of the adult data set derived from selected equation-paired concentration time point combinations compared to the reference standard.
Table 2. Summary of pharmacokinetic parameters of the adult data set derived from selected equation-paired concentration time point combinations compared to the reference standard.
Adult Data SetReference StandardModel 2
C1 at 60 and C2 at 240 min after the Completed InfusionC1 at 90 and C2 at 240 min after the Completed Infusion
AUC, mean ± SD (mg/L × h)427.14 ± 135.26440.11 ± 132.85445.09 ± 138.21
Vancomycin Cl, mean ± SD (L/h)4.62 ± 1.454.47 ± 1.384.43 ± 1.40
Vd, mean ± SD (L)39.35 ± 8.9539.88 ± 8.9241.33 ± 7.76
Half-life ± SD (h)6.28 ± 1.976.56 ± 2.026.83 ± 1.69
Abbreviations: AUC, area under the concentration curve; Cl, clearance; Vd, volume of distribution.
Table 3. Summary of pharmacokinetic parameters of the pediatric data set derived from selected equation-paired concentration time point combinations compared to the reference standard.
Table 3. Summary of pharmacokinetic parameters of the pediatric data set derived from selected equation-paired concentration time point combinations compared to the reference standard.
Pediatric Data SetReference StandardModel 2
C1 at 60 and C2 at 240 min after the Completed InfusionC1 at 60 and C2 at 300 min after the Completed Infusion
AUC, mean ± SD (mg/L × h)111.91 ± 42.73111.78 ± 39.83114.03 ± 43.39
Vancomycin Cl, mean ± SD (L/kg/h)0.16 ± 0.040.16 ± 0.040.15 ± 0.04
Vd, mean ± SD (L/kg)0.55 ± 0.100.51 ± 0.060.52 ± 0.07
Half-life ± SD (h)2.58 ± 0.912.38 ± 0.692.47 ± 0.84
Abbreviations: AUC, area under the concentration curve; Cl, clearance; Vd, volume of distribution.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sujjavorakul, K.; Katip, W.; Kerr, S.J.; Wacharachaisurapol, N.; Puthanakit, T. Predicting the Area under the Plasma Concentration-Time Curve (AUC) for First Dose Vancomycin Using First-Order Pharmacokinetic Equations. Antibiotics 2023, 12, 630. https://doi.org/10.3390/antibiotics12040630

AMA Style

Sujjavorakul K, Katip W, Kerr SJ, Wacharachaisurapol N, Puthanakit T. Predicting the Area under the Plasma Concentration-Time Curve (AUC) for First Dose Vancomycin Using First-Order Pharmacokinetic Equations. Antibiotics. 2023; 12(4):630. https://doi.org/10.3390/antibiotics12040630

Chicago/Turabian Style

Sujjavorakul, Kritsaporn, Wasan Katip, Stephen J. Kerr, Noppadol Wacharachaisurapol, and Thanyawee Puthanakit. 2023. "Predicting the Area under the Plasma Concentration-Time Curve (AUC) for First Dose Vancomycin Using First-Order Pharmacokinetic Equations" Antibiotics 12, no. 4: 630. https://doi.org/10.3390/antibiotics12040630

APA Style

Sujjavorakul, K., Katip, W., Kerr, S. J., Wacharachaisurapol, N., & Puthanakit, T. (2023). Predicting the Area under the Plasma Concentration-Time Curve (AUC) for First Dose Vancomycin Using First-Order Pharmacokinetic Equations. Antibiotics, 12(4), 630. https://doi.org/10.3390/antibiotics12040630

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop