Vancomycin pharmacokinetic models: informing the clinical management of drug-resistant bacterial infections

Therapeutic Drug Monitoring (TDM) Implementation in Public Hospitals in Greece in 2003 and 2021: A Comparative Analysis of TDM Evolution over the Years

Therapeutic drug monitoring (TDM) is the clinical practice of measuring drug concentrations. TDM can be used to determine treatment efficacy and to prevent the occurrence or reduce the risk of drug-induced side effects, being, thus, a tool of personalized medicine. Drugs for which TDM is applied should have a narrow therapeutic range and exhibit both significant pharmacokinetic variability and a predefined target concentration range. The aim of our study was to assess the current status of TDM in Greek public hospitals and estimate its progress over the last 20 years. All Greek public hospitals were contacted to provide data and details on the clinical uptake of TDM in Greece for the years 2003 and 2021 through a structured questionnaire. Data from 113 out of 132 Greek hospitals were collected in 2003, whereas for 2021, we have collected data from 98 out of 122 hospitals. Among these, in 2003 and 2021, 64 and 51 hospitals, respectively, performed TDM. Antiepileptics and antibiotics were the most common drug categories monitored in both years. The total number of drug measurement assays decreased from 2003 to 2021 (153,313 ± 7794 vs. 90,065 ± 5698; p = 0.043). In direct comparisons between hospitals where TDM was performed both in 2003 and 2021 (n = 35), the mean number of measurements was found to decrease for most drugs, including carbamazepine (198.8 ± 46.6 vs. 46.6 ± 10.1, p < 0.001), phenytoin (253.6 ± 59 vs. 120 ± 34.3; p = 0.001), amikacin (147.3 ± 65.2 vs. 91.1 ± 71.4; p = 0.033), digoxin (783.2 ± 226.70 vs. 165.9 ± 28.9; p < 0.001), and theophylline (71.5 ± 28.7 vs. 11.9 ± 6.4; p = 0.004). Only for vancomycin, a significant increase in measurements was recorded (206.1 ± 96.1 vs. 789.1 ± 282.8; p = 0.012). In conclusion, our findings show that TDM clinical implementation is losing ground in Greek hospitals. Efforts and initiatives to reverse this trend are urgently needed.

Two Innovative Approaches to Optimize Vancomycin Dosing Using Estimated AUC after First Dose: Validation Using Data Generated from Population PK Model Coupled with Monte-Carlo Simulation and Comparison with the First-Order PK Equation Approach

The revised consensus guidelines for optimizing vancomycin doses suggest that maintaining the area under the concentration-time curve to minimal inhibitory concentration ratio (AUC/MIC) of 400–600 mg·h/L is the target pharmacokinetic/pharmacodynamic (PK/PD) index for efficacy. AUC-guided dosing approach uses a first-order pharmacokinetics (PK) equation to estimate AUC using two samples obtained at steady state and one-compartment model, which can cause inaccurate AUC estimation and fail to achieve the effective PK/PD target early in therapy (days 1 and 2). To achieve an efficacy target from the third or fourth dose, two innovative approaches (Method 1 and Method 2) to estimate vancomycin AUC at steady state (AUCSS) using two-compartment model and three or four levels after the first dose are proposed. The feasibility of the proposed methods was evaluated and compared with another published dosing algorithm (Method 3), which uses two samples and a one-compartment approach. Monte Carlo simulation was performed using a well-established population PK model, and concentration-time profiles for virtual patients with various degrees of renal function were generated, with 1000 subjects per group. AUC extrapolated to infinity (AUC0–∞) after the first dose was estimated using the three methods, whereas reference AUC (AUCref) was calculated using the linear-trapezoidal method at steady state after repeated doses. The ratio of AUC0–∞: AUCref and % bias were selected as the indicators to evaluate the accuracy of three methods. Sensitivity analysis was performed to examine the influence of change in each sampling time on the estimated AUC0–∞ using the two proposed approaches. For simulated patients with various creatinine clearance, the mean of AUC0–∞: AUCref obtained from Method 1, Method 2 and Method 3 ranged between 0.98 to 1, 0.96 to 0.99, and 0.44 to 0.69, respectively. The mean bias observed with the three methods was −0.10% to −2.09%, −1.30% to −3.59% and −30.75% to −55.53%, respectively. The largest mean bias observed by changing sampling time while using Method 1 and Method 2 were −4.30% and −10.50%, respectively. Three user-friendly and easy-to-use excel calculators were built based on the two proposed methods. The results showed that our approaches ensured sufficient accuracy and achieved target PK/PD index early and were superior to the published methodologies. Our methodology has the potential to be used for vancomycin dose optimization and can be easily implemented in clinical practice.

Diagnostic and medical needs for therapeutic drug monitoring of antibiotics

Therapeutic drug monitoring (TDM) of antibiotics has been practiced for more than half a century, but it is still not widely applied for infected patients. It has a traditional focus on limiting toxicity of specific classes of antibiotics such as aminoglycosides and vancomycin. With more patients in critical care with higher levels of sickness severity and immunosuppression as well as an increasingly obese and ageing population, an increasing risk of suboptimal antibiotic exposure continues to escalate. As such, the value of TDM continues to expand, especially for beta-lactams which constitute the most frequently used antibiotic class. To date, the minimum inhibitory concentration (MIC) of infectious microbes rather than classification in terms of susceptible and resistant can be reported. In parallel, increasingly sophisticated TDM technology is becoming available ensuring that TDM is feasible and can deliver personalized antibiotic dosing schemes. There is an obvious need for extensive studies that will quantify the improvements in clinical outcome of individual TDM-guided dosing. We suggest that a broad diagnostic and medical investigation of the TDM arena, including market analyses and analytical technology assessment, is a current priority.

An Evaluation of Vancomycin Area Under the Curve Estimation Methods for Children Treated for Acute Pulmonary Exacerbations of Cystic Fibrosis Due to Methicillin-Resistant Staphylococcus aureus

The prevalence of pulmonary methicillin-resistant Staphylococcus aureus infections in patients with cystic fibrosis (CF) has increased over the last 2 decades. Two concentrations-a postdistributive and a trough-are currently used to estimate the area under the curve (AUC) of vancomycin, an antibiotic routinely used to treat these infections, to achieve the target AUC/minimum inhibitory concentration of ≥400 mg·h/L in ensuring optimal dosing of this drug. This study evaluated precision and bias in estimating vancomycin AUCs obtained either from a population pharmacokinetic (PK) model by using a single trough concentration or from standard PK equation-based 2-point monitoring approach. AUCs were either obtained from a single trough concentration-fitted model or derived from a model fitted by 2 concentration points. Children ≥2 years of age with CF received intravenous vancomycin at 2 centers from June 2012 to December 2014. A population PK model was developed in Pmetrics to quantify the between-subject variability in vancomycin PK parameters, define the sources of PK variability, and leverage information from the population to improve individual AUC estimates. Twenty-three children with CF received 27 courses of vancomycin. The median age was 12.3 (interquartile range [IQR] 8.5-16.6) years. From the individual vancomycin PK parameter estimates from the population PK model, median AUC was 622 (IQR 529-680) mg·h/L. Values were not significantly different from the AUC calculated using the standard PK equation-based approach (median 616 [IQR 540-663] mg·h/L) (P = .89). A standard PK equation-based approach using 2 concentrations and a population PK model-based approach using a single trough concentration yielded unbiased and precise AUC estimates. Findings suggest that options exist to implement AUC-based pediatric vancomycin dosing in patients with CF. The findings of this study reveal that several excellent options exist for centers to implement AUC-based pediatric vancomycin dosing for patients with CF.

Population pharmacokinetics of vancomycin and AUC-guided dosing in Chinese neonates and young infants

OBJECTIVES

To develop a population pharmacokinetic (PK) model for vancomycin in Chinese neonates and infants less than 2 months of age (young infants) with a wide gestational age range, in order to determine the appropriate dosing regimen for this population.

METHODS

We performed a retrospective chart review of patients from the neonatal intensive care unit (NICU) at Children's Hospital of Fudan University to identify neonates and young infants treated with vancomycin from May 2014 to May 2017. Vancomycin concentrations and covariates were utilized to develop a one-compartment model with first-order elimination. The predictive performance of the final model was assessed by both internal and external evaluation, and the relationship between trough concentration and AUC_0-24 was investigated. Monte Carlo simulations were performed to design an initial dosing schedule targeting an AUC_0-24 ≥ 400.

RESULTS

The analysis included a total of 330 concentration-time data points from 213 neonates and young infants with gestational age (GA) and body weight of 25-42 weeks and 0.88-5.1 kg, respectively. Body weight, postmenstrual age (PMA) and serum creatinine level were found to be important factors explaining the between-subject variability in vancomycin PK parameters for this population. Both internal and external evaluation supported the prediction of the final vancomycin PK model. The typical population parameter estimates of clearance and distribution volume for an infant weighing 2.73 kg with a PMA of 39.8 weeks and serum creatinine of 0.28 mg/dL were 0.103 L/h/kg and 0.58 L/kg, respectively. Although vancomycin serum trough concentrations were predictive of the AUC, considerable variability was observed in the achievement of an AUC_0-24/MIC of ≥400. For MIC values of ≤0.5 mg/L, AUC_0-24/MIC ≥400 was achieved for 95% of the newborn infants with vancomycin troughs of 5-10 mg/L. When the MIC increased to 1 mg/L, only 15% of the patients with troughs of 5-10 mg/L achieved AUC_0-24/MIC ≥400. For MIC values of 2 mg/L, no infants achieved the target. Simulations predicted that a dose of at least 14 and 15 mg/kg every 12 h was required to attain the target AUC_0-24 ≥ 400 in 90% of infants with a PMA of 30-32 and 32-34 weeks, respectively. This target was also achieved in 93% of simulated infants in the oldest PMA groups (36-38 and 38-40 weeks, respectively) when the dosing interval was extended to 8 h. For infants with a PMA ≥44 weeks, a dose increase to 18 mg/kg every 8 h was needed. The trough concentrations of 5-15 mg/L were highly predictive of an AUC_0-24 of ≥400 when treating invasive MRSA infections with an MIC of ≤1 mg/L.

CONCLUSIONS

The PK parameters for vancomycin in Chinese infants younger than 2 months of age were estimated using the model developed herein. This model has been used to predict individualized dosing regimens in this vulnerable population in our hospital. A large external evaluation of our model will be conducted in future studies.

Using a Vancomycin PBPK Model in Special Populations to Elucidate Case-Based Clinical PK Observations

Simultaneous changes in several physiological factors may contribute to the large pharmacokinetic (PK) variability of vancomycin. This study was designed to systematically characterize the effects of multiple physiological factors to the altered PK of vancomycin observed in special populations. A vancomycin physiologically based pharmacokinetic (PBPK) model was developed as a PK simulation platform to quantitatively assess the effects of changes in physiologies to the PK profiles. The developed model predicted the concentration-time profiles in healthy adults and diseased patients. The implementation of developmental changes in both renal and non-renal elimination pathways to the pediatric model improved the predictability of vancomycin clearance. Simulated PK profiles with a 50% decrease in cardiac output (peak plasma concentration (C_max ), 59.9 ng/mL) were similar to those observed in patients before bypass surgery (C_max , 55.1 ng/mL). The PBPK modeling of vancomycin demonstrated its potential to provide mechanistic insights into the altered disposition observed in patients who have changes in multiple physiological factors.

Population Pharmacokinetics and Pharmacodynamic Target Attainment of Vancomycin in Neonates on Extracorporeal Life Support

OBJECTIVES

To evaluate the population pharmacokinetics and pharmacodynamic target attainment of vancomycin in neonates with a contemporary ¼-inch extracorporeal life support circuit with a Quadrox-iD Pediatric oxygenator (Maquet Cardiovascular, LLC, Wayne, NJ).

DESIGN

Retrospective medical record review.

SETTING

Two free-standing tertiary/quaternary pediatric children's hospitals.

PATIENTS

Neonates receiving either veno-arterial or veno-venous extracorporeal life support and vancomycin for empiric or definitive therapy with resulting serum concentrations.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

Twelve patients with a median gestations age of 39 weeks (range 36-41 wk) and a median postnatal age of 9.5 days (range 0-28 d) accounted for 14 courses of vancomycin therapy while on extracorporeal life support and were included in the analysis. The median weight was 3.1 kg (range 2.2-4.41 kg) with five of 12 patients (41.7%) being female. Vancomycin concentrations were best described by an one-compartment model incorporating allometric scaling of estimated glomerular filtration rate on clearance. The mean total body clearance (mL/min/kg) for the population was 3.48 ± 1.31 mL/min/kg, and the mean total volume of distribution (L/kg) for the population was 1.2 ± 0.4 L/kg. The intermittent and continuous infusion dosing regimens that provided for the highest percentage of trough concentrations in the range of 10-20 mg/L were the 10 mg/kg/dose IV q8h, 12.5 mg/kg/dose IV q8-12h, 15 mg/kg/dose IV q12h, and 20 mg/kg/dose IV q12h, and the 20, 25, and 30 mg/kg/d continuous infusion regimens, respectively. All regimens allowed for an area under the concentration:minimum inhibitory concentration ratio of 400:1 for minimum inhibitory concentrations of less than or equal to 0.5 mg/L for a 90% PTA. None of the simulated regimens had a greater than 90% probability of achieving an area under the concentration:minimum inhibitory concentration ratio of 400:1 for vancomycin minimum inhibitory concentrations greater than or equal to 1 mg/L while maintaining trough concentrations in the range of 10-20 mg/L.

CONCLUSIONS

To our knowledge, this is the first pharmacokinetic and pharmacodynamic study of neonates receiving vancomycin with a contemporary ¼-inch extracorporeal life support circuit including the Quadrox-iD Pediatric oxygenator (Maquet Cardiovascular, LLC). The data suggest differences in vancomycin pharmacokinetics compared with previous extracorporeal life support data, notably a more rapid clearance, which could result in lower vancomycin concentrations. Considering this, a more aggressive initial dosing regimen may need to be employed in infants on extracorporeal life support.

Coagulase negative staphylococcal sepsis in neonates: do we need to adapt vancomycin dose or target?

Background

Despite differences in types of infection and causative organisms, pharmacokinetic-pharmacodynamic (PKPD) targets of vancomycin therapy derived from adult studies are suggested for neonates. We aimed to identify doses needed for the attainment of AUC/MIC > 400 and AUC/MIC > 300 in neonates with sepsis and correlate these targets with recommended doses and treatment outcome.

Methods

Neonates who had Vancomycin therapeutic drug monitoring (TDM) performed between January 1, 2010 and December 31, 2012 were studied. Clinical characteristics, episodes of Gram-positive sepsis with outcomes and all neonatal blood culture isolates in hospital were collected from medical records. To estimate probability of target attainment of AUC/MIC >400 and AUC/MIC >300 a 1000-subject Monte Carlo simulation was performed by calculating AUC using Anderson’s (Anderson et al. 2006) and TDM trough concentrations (C_trough) based population PK models.

Results

Final dataset included 76 patients; 57 with confirmed Gram-positive sepsis. TDM was taken after the 1^st to 44^th dose. 84.1% of C_trough were within the range 5–15 mg/L. Currently recommended doses achieved probability of the targets (PTA) of AUC/MIC >400 and AUC/MIC >300 in less than 25% and 40% of cases, respectively. Doses required for 80% PTA of AUC/MIC > 400 for MIC ≥2 mg/L resulted in C_trough values ≥14 mg/L. Mean AUC/MIC values were similar in treatment failure and success groups.

Conclusion

With currently recommended vancomycin dosing the therapeutic target of AUC/MIC > 400 is achieved only by 25% of neonates. Appropriate PKPD targets and respective dosing regimens need to be defined in prospective clinical studies in this population.

A propensity-matched cohort study of vancomycin-associated nephrotoxicity in neonates

BACKGROUND

The incidence of nephrotoxicity among vancomycin-treated neonates has been reported to range from 2% to 20%. These widely varying estimates have led to confusion and controversy regarding the safety of vancomycin among neonates.

OBJECTIVE

Evaluate the incidence of nephrotoxicity among neonates receiving vancomycin concomitantly with gentamicin.

DESIGN

Retrospective observational cohort study using propensity score matching to provide covariate balance between neonates who did or did not receive vancomycin based on factors known to be related to the development of renal dysfunction.

SETTING

Hospitals (n=22) throughout the Intermountain West, including a quaternary care children's hospital.

PATIENTS

Neonates ≤44 postmenstrual weeks (median gestational age: 31 (IQR 28-36) weeks) receiving intravenous gentamicin with or without exposure to vancomycin from January 2006 to December 2012.

MAIN OUTCOME MEASURES

Nephrotoxicity based on the modified Acute Kidney Injury Network criteria for acute kidney injury (AKI) or serum creatinine concentration ≥1.5 mg/dL persisting for ≥48 h.

RESULTS

The final cohort was comprised of 1066 neonates (533 receiving vancomycin and gentamicin vs 533 receiving gentamicin). In a propensity score-matched cohort that was well balanced across 16 covariates, AKI was not associated with vancomycin use (16 neonates receiving vancomycin vs 7 controls experienced AKI; OR 1.5; 95% CI 0.6 to 4.0). However, the presence of a patent ductus arteriosus, concomitant non-steroidal anti-inflammatory drug use, ≥1 positive blood cultures, low birth weight and higher severity of illness and risk of mortality scores were associated with an increased risk of nephrotoxicity.

CONCLUSIONS

These results corroborate several earlier reports and much anecdotal evidence describing the infrequent occurrence of nephrotoxicity in neonates receiving concomitant vancomycin and gentamicin.

Towards Rational Dosing Algorithms for Vancomycin in Neonates and Infants Based on Population Pharmacokinetic Modeling

Because of the recent awareness that vancomycin doses should aim to meet a target area under the concentration-time curve (AUC) instead of trough concentrations, more aggressive dosing regimens are warranted also in the pediatric population. In this study, both neonatal and pediatric pharmacokinetic models for vancomycin were externally evaluated and subsequently used to derive model-based dosing algorithms for neonates, infants, and children. For the external validation, predictions from previously published pharmacokinetic models were compared to new data. Simulations were performed in order to evaluate current dosing regimens and to propose a model-based dosing algorithm. The AUC/MIC over 24 h (AUC24/MIC) was evaluated for all investigated dosing schedules (target of >400), without any concentration exceeding 40 mg/liter. Both the neonatal and pediatric models of vancomycin performed well in the external data sets, resulting in concentrations that were predicted correctly and without bias. For neonates, a dosing algorithm based on body weight at birth and postnatal age is proposed, with daily doses divided over three to four doses. For infants aged <1 year, doses between 32 and 60 mg/kg/day over four doses are proposed, while above 1 year of age, 60 mg/kg/day seems appropriate. As the time to reach steady-state concentrations varies from 155 h in preterm infants to 36 h in children aged >1 year, an initial loading dose is proposed. Based on the externally validated neonatal and pediatric vancomycin models, novel dosing algorithms are proposed for neonates and children aged <1 year. For children aged 1 year and older, the currently advised maintenance dose of 60 mg/kg/day seems appropriate.

Predictive Performance of a Vancomycin Population Pharmacokinetic Model in Neonates

INTRODUCTION

The pharmacokinetics of vancomycin are highly variable among neonates, which makes dosing challenging in this population. However, adequate drug exposure is critical, especially when treating methicillin-resistant Staphylococcus aureus (MRSA) infections. Utilization of population pharmacokinetic models and Bayesian methods offers the potential for developing individualized therapeutic approaches. To meet this need, a neonatal vancomycin population pharmacokinetic model was recently published. The current study sought to externally evaluate the predictive performance and generalizability of this model.

METHODS

A retrospective chart review of neonates who received vancomycin and had ≥1 peak and ≥1 trough concentrations at five Intermountain Healthcare neonatal intensive care units from 2006 to 2013 was performed and served as the external validation cohort. The published population pharmacokinetic model was implemented in NONMEM 7.2 with the structural and variance parameter values set equal to the estimates reported previously. The model was then used to predict the first peak and trough concentration for each neonate in the validation cohort and the model prediction error and absolute prediction error were calculated. Normalized prediction distribution errors (NPDE) were also evaluated.

RESULTS

A total of 243 neonates were studied with a median postmenstrual age of 33 (range: 23-54) weeks and a median weight of 1.6 (range: 0.4-6.8) kg. The model predicted the observed vancomycin concentrations with reasonable precision. For all vancomycin concentrations, the median prediction error was -0.8 (95% CI: -1.4 to -0.4) mg/L and the median absolute prediction error was 3.0 (95% CI: 2.7-3.5) mg/L. No trends in NPDE across weight, postmenstrual age, serum creatinine, or time after dose were observed.

CONCLUSION

An evaluation of a recently published neonatal vancomycin population pharmacokinetic model in a large external dataset supported the predictive performance and generalizability of the model. This model may be useful in evaluating neonatal vancomycin dosing regimens and estimating the extent of drug exposure.