key: cord-010903-kuwy7pbo
authors: Liu, Jiajun; Neely, Michael; Lipman, Jeffrey; Sime, Fekade; Roberts, Jason A.; Kiel, Patrick J.; Avedissian, Sean N.; Rhodes, Nathaniel J.; Scheetz, Marc H.
title: Development of Population and Bayesian Models for Applied Use in Patients Receiving Cefepime
date: 2020-03-05
journal: Clin Pharmacokinet
DOI: 10.1007/s40262-020-00873-3
sha: 
doc_id: 10903
cord_uid: kuwy7pbo

BACKGROUND AND OBJECTIVE: Understanding pharmacokinetic disposition of cefepime, a β-lactam antibiotic, is crucial for developing regimens to achieve optimal exposure and improved clinical outcomes. This study sought to develop and evaluate a unified population pharmacokinetic model in both pediatric and adult patients receiving cefepime treatment. METHODS: Multiple physiologically relevant models were fit to pediatric and adult subject data. To evaluate the final model performance, a withheld group of 12 pediatric patients and two separate adult populations were assessed. RESULTS: Seventy subjects with a total of 604 cefepime concentrations were included in this study. All adults (n = 34) on average weighed 82.7 kg and displayed a mean creatinine clearance of 106.7 mL/min. All pediatric subjects (n = 36) had mean weight and creatinine clearance of 16.0 kg and 195.6 mL/min, respectively. A covariate-adjusted two-compartment model described the observed concentrations well (population model R(2), 87.0%; Bayesian model R(2), 96.5%). In the evaluation subsets, the model performed similarly well (population R(2), 84.0%; Bayesian R(2), 90.2%). CONCLUSION: The identified model serves well for population dosing and as a Bayesian prior for precision dosing. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s40262-020-00873-3) contains supplementary material, which is available to authorized users.

Cefepime is a commonly utilized antibiotic for nosocomial infections. Rising resistance, manifesting as increased cefepime minimum inhibitory concentrations (MICs), has led to more frequent clinical failures [1, 2] . To advise clinical outcomes according to MICs, the Clinical and Laboratory Standards Institute updated the susceptibility breakpoints and then created a category of susceptibledose dependent for MICs of 4 and 8 mg/L for Enterobacteriaceae spp. [3] . Achieving goal pharmacokinetic exposures to effectively treat these higher MICs can require a precision dosing approach.

The online version of this article (https ://doi.org/10.1007/s4026 2-020-00873 -3) contains supplementary material, which is available to authorized users.

Cefepime, like other β-lactams, has pharmacodynamic activity governed by 'time-dependent' activity. The fraction of time that the unbound drug concentration exceeds the MIC (fT >MIC ) for the dosing interval is the pharmacokinetic/pharmacodynamic (PK/PD) efficacy target for cefepime [4] , and a target of 60-74% has been previously proposed [5] [6] [7] [8] . For the currently approved cefepime product and combination agents in the pipeline [9, 10] , understanding cefepime disposition and variability is crucial to for optimal treatment of patients. As inter-and intra-patient PK variability can impact the achievement of PD goals, understanding the precision of population dosing is important. Further, to fully realize precision dosing, individualized models (e.g., Bayesian models) are needed. Once developed, these models will form the basis for adaptive feedback and control strategies when paired with real-time drug assays. The purpose of this study was to: (1) develop and evaluate a unified cefepime population PK model for adult and pediatric patients, and (2) construct an individualized model that can be utilized to deliver precision cefepime dosing.

Data from four clinical cefepime PK studies representing unique groups of patients were compiled. Subject demographics and study methodologies have been previously described [11] [12] [13] [14] . In brief, populations represented were febrile neutropenic adults with hematologic malignancies [13, 14] , those with critical illness [12] , and children with presumed or documented bacterial infections [11] . For the two studies that evaluated adults with neutropenic fever, Sime et al. prospectively enrolled 12 patients receiving chemotherapy and/or stem cell transplant who subsequently developed febrile neutropenia and were administered maximum doses of cefepime [13] . A total of 53 cefepime plasma concentrations in presumably steady-state dosing intervals (third, sixth, and ninth) were analyzed for PK target attainment. Whited et al. prospectively studied similar patients (n = 9) who were admitted to hematology-oncology services and were receiving cefepime at a maximum dosage for febrile neutropenia [14] . Cefepime PK samples were obtained during steady state and analyzed for population parameters. Critically ill adults were studied by Roberts et al. as a prospective multinational PK study and included 14 patients who received cefepime (only n = 13 were included for model evaluation) [12] . Last, Reed et al. characterized cefepime pharmacokinetics in hospitalized pediatric patients (above 2 months of age) who received cefepime as monotherapy for bacterial infections [11] . For our study, only those who received intravenous cefepime were included for model development.

Adult (n = 9) and partial pediatric (n = 24) datasets were utilized for PK model building ( Fig. 1) [11, 14] . Model evaluation was performed with other datasets consisting of independent adult (n = 12, n = 13) and pediatric (n = 12) patients [11] [12] [13] 15] . Pediatric patients from Reed et al. [11] were randomized into the model building or the evaluation dataset. All clinical patient-level data included age, weight, and serum creatinine. An estimated creatinine clearance (CrCL) was calculated for each patient [16] . The Cockcroft-Gault formula (applied to all subjects) served as a standardized descriptor for the elimination rate constant. This study was exempted by the Institutional Review Board at Midwestern University Chicago College of Pharmacy.

To construct the base PK models, the nonparametric adaptive grid (NPAG) algorithm [17, 18] within the Pmetrics (Version 1.5.2) package [18] for R [19] was utilized. Multiple physiologically relevant, one-and two-compartmental PK models were built and assessed. The one-compartment structural model included an intravenous cefepime dose into and parameterized total cefepime elimination rate constant (K e ) from the central compartment. The two-compartment model included additional parameterizations of intercompartmental transfer constants between central and peripheral compartments (K CP and K PC ). In candidate models, total cefepime elimination was explored according to full renal and partial renal clearance (CL) models [i.e., nonrenal elimination (K e Intcpt) and renal elimination descriptor (K e 0 vectorized as a function of glomerular filtration estimates)] [9, 20] .

Assay error was included into the model using a polynomial equation in the form of standard deviation (SD)

A unified cefepime population pharmacokinetic model has been developed from adult and pediatric patients and evaluates well in independent populations.

When paired with real-time β-lactam assays, a precision dosing approach will optimize drug exposure and improve clinical outcomes. as a function of each observed concentration, Y (i.e., SD = C 0 + C 1 · Y). Observation weighting was performed using gamma (i.e., error = SD · gamma), a multiplicative variance model to account for extra process noise. Gamma was initially set at 4 with C 0 and C 1 equal to 0.5 and 0.15, respectively.

Covariate relationships were assessed using the 'PMStep' function in Pmetrics by applying stepwise linear regressions (forward selection and backwards elimination) of all covariates on PK parameters. Additionally, a priori analyses examined the effect of covariates on cefepime K e , and both weight and CrCL were variables considered a priori to have a high potential likelihood to impact cefepime pharmacokinetics ( [9, 21, 22] . Weight and CrCL were standardized to 70 kg and 120 mL/min, respectively. Further, an allometric scaler was applied to standardized weight (i.e., quotient of weight in kg divided by 70 kg raised to the negative 0.25th power) as a covariate adjustment to K e (ESM). Ultimate model retention was governed according to criteria described below.

The best-fit PK and error model was identified by the change in objective function value (OFV) calculated as differences in − 2 log-likelihood, with a reduction of 3.84 in OFV corresponding to p < 0.05 based on Chi-square distribution and one degree of freedom. Further, the best-fit model was selected based on the rule of parsimony and the lowest Akaike's information criterion scores. Goodness of fit of the competing models were evaluated by regression on observed vs. predicted plots, coefficients of determination, and visual predictive checks. Predictive performance was assessed using bias and imprecision in both population and individual prediction models. Bias was defined as mean weighted prediction error; imprecision was defined as bias-adjusted mean weighted squared prediction error. Posterior-predicted cefepime concentrations for each study subject were calculated using individual median Bayesian posterior parameter estimates.

To evaluate the final adjusted model, the NPAG algorithm [17, 18] was employed to assess the performance with separate data sets (Fig. 1 ). The population joint density from the best-fit covariate adjusted model was employed as a Bayesian prior for the randomly withheld pediatric data 

Adult Data 3 Adult Data 4 Adult Data 1 Fig. 1 Schematic for data sources in model development and evaluation and separate adult data. In the evaluation process, structural model, model parameters, assay error, and observation weighting were unchanged. Goodness of fit of the competing models were determined as described above.

Simulation was performed to examine the exposures predicted by the final model, employing all support points from the population parameter joint density in the final NPAG analysis [18, 23] . Each support point was treated as a mean vector surrounded by the population variance-covariance matrix (i.e., covariance equal to the population covariance divided by the total number of support points). For each subject, 1000 simulated profiles were created with predicted outputs at 0.1-h intervals. Covariate values for each simulated subject were fixed based on arithmetic means of observed weight and CrCL for corresponding adult and pediatric populations. Semi-parametric Monte Carlo sampling was performed from the multimodal multivariate distribution of parameters with the parameter space concordant with the NPAG population analysis results (i.e., best-fit model) [ Table S1 of the ESM] [23] . Maximum dosing regimens were simulated for adult and pediatric populations (total n = 33): 2 g every 8 h infused over 0.5 h and 50 mg/kg every 8 h infused over 0.5 h, respectively. Protein binding of 20% (i.e., 80% free fraction of total cefepime dose) was accounted for in predicting cefepime concentrations [9] . The PK/PD target of fT >MIC ≥ 68% was utilized across doubling MICs of 0.25-32 mg/L over the first 24 h of cefepime therapy [5] .

Estimates are provided from the first 24 h of simulations as timely administration of effective antimicrobial agents is associated with increased survival [24] .

A total of 70 clinically diverse subjects, contributing 683 cefepime concentrations, were included in this study (n = 33 subjects for model development; n = 37 subjects for evaluation) (Fig. 1 

A total of 428 cefepime observations were available for model development. Cefepime concentrations ranged from 0.5 to 249.7 μg/mL. The base one-and two-compartment models (without covariate adjustment) produced reasonable fits for observed and Bayesian posterior-predicted cefepime concentrations (R 2 = 84.7% and 85.2%, respectively), but population estimates were unsatisfactory (R 2 = 22.7% and 27.8%, respectively) ( Table 1) . Weight and CrCL displayed relationships with the standard two-compartment model (i.e., base two-compartment model). Volume of distribution was associated with weight (p < 0.2) and Ke (total) was associated with CrCL (p < 0.2). After standardizing weight (to 70 kg) without an allometric scaler in the base two-compartment model, fits for both population and Bayesian posterior estimates against the observed data improved (R 2 = 60.7% and 96.5%, respectively; OFV change, 4). Bias and imprecision for Bayesian posterior fits were − 0.18 and 1.12, respectively. When covariates (i.e., weight to volume of distribution and K e ; CrCL to K e ) and the allometric scaler were applied in the two-compartment model, Bayesian posteriors fit well (R 2 = 96.5%; Fig. 2 right) with low bias and imprecision (− 0.15 and 1.07, respectively), and the population PK model produced good fits of the observed cefepime concentrations (R 2 = 87.0%, bias = 0.53, imprecision = 7.75; Fig. 2 left) . The OFV change from the weightadjusted two-compartment model to the final model was significant at − 34 (p < 0.05) [ Table 1 ]. The final model also produced acceptable predicted checks (Fig. 3) . Thus, a two-compartment model with weight and CrCL as covariate adjustments and allometric scaling was selected as the final PK model. The population parameter values from the final PK model are summarized in Table 2 . Structural model and differential equations that define the population PK are listed in the ESM. The population parameter value covariance matrix can be found in Table 3 . Additionally, weighted residual error plots for the best-fit model (Fig. S1 ) and scatter plots for covariates for the base structural model (Fig. S2) can be found included in the ESM. For the evaluation subset, Bayesian priors resulted in reasonably accurate and precise predictions (population R 2 = 84.0%, Bayesian R 2 = 90.2%; Fig. 4 ). Fig. 2 Goodness-of-fit plots for best-fit population cefepime PK model (model development)

Results of the probability of target attainment (PTA) analysis are shown in Table 4 

This study created a population and individual PK model for adult and pediatric patients and can serve as a Bayesian prior for precision dosing. When paired with a real-time assay for cefepime, this model allows for precise and accurate predictions of cefepime disposition via adaptive feedback control. In the absence of real-time assays, these cefepime PK parameters facilitate more accurate population-based dosing Table 3 Population parameter value covariance matrix for the best-fit model strategies. Previous work by Rhodes et al. has shown an absolute difference of approximately 20% in survival probability across the continuum of achieving 0-100% fT >MIC in adult patients with Gram-negative bloodstream infections, thus understanding the dose and re-dosing interval necessary to achieve optimal PK exposures should greatly improve clinical outcomes for patients treated with cefepime [5] . Individualized dosing and therapeutic drug monitoring of β-lactam antibiotics (e.g., cefepime) are critically important to achieving optimal drug exposure (i.e., optimal fT >MIC as the PK/PD target) and improving clinical outcomes [4, 25, 26] . Precision medicine has been named as a major focus for the National Health Institute with $215 million invested [27] , yet precision medicine has mostly focused on genomic differences [28, 29] . Precision dosing is an important facet of precision medicine, and renewed efforts in precision dosing in the real-world setting are being pursued [30] . Cefepime is a highly relevant example. While rigorous reviews and analyses are conducted during the development phase of an antibiotic, dose optimization is far less ideal for the types of patients who ultimately receive the drug. This is highlighted by the fact that although cefepime-associated neurotoxicity is rare, this serious and potentially life-threatening adverse event has been increasingly reported and few strategies exist for optimizing and delivering precision exposures [31, 32] . Lamoth et al. found that a cefepime trough concentration of ≥ 22 mg/L has a 50% probability of predicting neurotoxicity [33] . Huwyler et al. identified a similar predictive threshold of > 20 mg/dL (five-fold increased risk for neurologic events) [34] . In contrast, Rhodes et al. found the cut-off of 22 mg/L to be suboptimal [35] . Furthermore, Rhodes et al. performed simulations from literature cefepime data and observed a high intercorrelation amongst all PK parameters (i.e., area under the curve at steady state, maximum plasma concentration, and minimum plasma concentration), suggesting that more work is needed to establish the pharmacokinetic/toxicodynamic (PK/TD) profile for cefepime. In addition to complications by these less-than-ideal PK/ TD data, clinicians are left to treat patients with extreme age differences, organ dysfunction, and comorbid conditions affecting antibiotic pharmacokinetics/pharmacodynamics [26] . Further, a contemporary dose reduction strategy based on estimated renal function (e.g., estimated CrCL using the Cockcroft-Gault formula) is also likely to be confounded in these patients by intrinsic PK variability, such as changes in volume of distribution, and the challenges of accurately estimating the glomerular filtration rate at any point in time, leading to more 'uncertainties' in balancing dose optimization and adverse events [9, 36] . These 'real-world' patients are often under-represented, and thus not well understood, from a PK/PD and PK/TD standpoint during the drug approval process. Bridging to the more typical patients that are clinically treated is important and central to the mission of precision medicine. The findings of this study can be used to guide cefepime dosing in these 'real-world' patients.

Several other studies have reviewed population cefepime pharmacokinetics. Sime [38] . In our pediatric population, means of CL and elimination half-life were 3.1 L/h and 3.0 h, respectively. Our simulation findings are similar to those of Shoji et al. that the maximum pediatric cefepime dosing did not adequately achieve optimal exposure to target higher MICs. While the cefepime package labeling recommends maximum dosages of 2 g every 8 h for adult patients with neutropenic fever and 50 mg/kg every 8 h for pediatric patients with pneumonia and/or neutropenic fever, there may be a need to extend these dosing regimens to other populations (in the absence of aforementioned indications) to achieve the best clinical outcomes by optimizing the PK/PD attainment goals [9] .

Other studies also performed a simulation for PTA with different cefepime regimens and renal functions. Tam et al. found that with a PD target of 67% f T>MIC, 2 g every 8 h (30-minute infusion) achieved approximately 90% PTA for MIC of 8 mg/L in patients with CrCL of 120 mL/min while 2 g every 12 h achieved barely above 80% PTA for an MIC of 4 mg/L in the same population [39] . Nicasio et al. also conducted a simulation using a PD target of 50% f T>MIC in the critically ill with varying renal function. The maximum recommended dosage (2 g every 8 h) in patients with CrCL between 50 and 120 mL/min achieved a PTA of 78.1% at an MIC of 16 mg/L; however, when the same regimen was infused over 0.5 h, the PTA achieved was significantly lower [37] . Collectively, these findings suggest that cefepime exposure is highly variable and may be clinically suboptimal in a large number of patients commonly treated with cefepime. These findings support the need for precision dosing and therapeutic drug monitoring for β-lactam antibiotics to reach optimal PK/PD targets given the high variability in drug exposures.

Our study is not without limitations. Although a relatively large and diverse cohort was included in model development and evaluation, we did not specifically assess certain subgroups such as patients with morbid obesity and severe renal dysfunction. These conditions may require patient-specific models. Second, many studies to date included 'real-world' patients with various disease sates (e.g., neutropenic fever, renal failure, sepsis); however, all studies were conducted under the research protocol where doses, and administration times were all carefully confirmed. Additional efforts will be needed to evaluate model performance in clinical contexts.

A unified population model for cefepime in adult and pediatric populations was developed and demonstrated excellent performance on evaluation. Current cefepime dosages are often suboptimal, and population variability is high. Precision dosing approaches and real-time assays are needed for cefepime to optimize drug exposure and improve clinical outcomes.

Failure of current cefepime breakpoints to predict clinical outcomes of bacteremia caused by Gram-negative organisms

Evaluation of clinical outcomes in patients with Gram-negative bloodstream infections according to cefepime MIC

Performance standards for antimicrobial susceptibility testing

Pharmacokinetic/pharmacodynamic parameters: rationale for antibacterial dosing of mice and men

Defining clinical exposures of cefepime for Gram-negative bloodstream infections that are associated with improved survival

Clinical pharmacodynamics of cefepime in patients infected with Pseudomonas aeruginosa

Relationship between PK/PD of cefepime and clinical outcome in febrile neutropenic patients with normal renal function

Interrelationship between pharmacokinetics and pharmacodynamics in determining dosage regimens for broad-spectrum cephalosporins

Cefepime/VNRX-5133 broad-spectrum activity is maintained against emerging KPC-and PDC-variants in multidrug-resistant K. pneumoniae and P. aeruginosa

Pharmacokinetics of intravenously and intramuscularly administered cefepime in infants and children

DALI: defining antibiotic levels in intensive care unit patients: are current beta-lactam antibiotic doses sufficient for critically ill patients?

Adequacy of high-dose cefepime regimen in febrile neutropenic patients with hematological malignancies

Pharmacokinetics of cefepime in patients with cancer and febrile neutropenia in the setting of hematologic malignancies or hematopoeitic cell transplantation. Pharmacotherapy

Guidance Document: Population Pharmacokinetics Guidance for Industry

Prediction of creatinine clearance from serum creatinine

An adaptive grid non-parametric approach to pharmacokinetic and dynamic (PK/PD) population models

Accurate detection of outliers and subpopulations with Pmetrics, a nonparametric and parametric pharmacometric modeling and simulation package for R

R: a language and environment for statistical computing

Cefepime clinical pharmacokinetics

Pharmacokinetics of cefepime in patients with respiratory tract infections

Cefepime in intensive care unit patients: validation of a population pharmacokinetic approach and influence of covariables

Population modeling and Monte Carlo simulation study of the pharmacokinetics and antituberculosis pharmacodynamics of rifampin in lungs

Duration of hypotension before initiation of effective antimicrobial therapy is the critical determinant of survival in human septic shock

Pharmacokinetics-pharmacodynamics of antimicrobial therapy: it's not just for mice anymore

Therapeutic drug monitoring of the beta-lactam antibiotics: what is the evidence and which patients should we be using it for?

The White House. FACT SHEET: President Obama's Precision Medicine Initiative

Precision medicine: from science to value. Health Aff (Millwood)

Precision medicine: changing the way we think about healthcare

Precision dosing: defining the need and approaches to deliver individualized drug dosing in the real-world setting

Cefepime and risk of seizure in patients not receiving dosage adjustments for kidney impairment

Characterizing cefepime neurotoxicity: a systematic review

High cefepime plasma concentrations and neurological toxicity in febrile neutropenic patients with mild impairment of renal function

Cefepime plasma concentrations and clinical toxicity: a retrospective cohort study

An exploratory analysis of the ability of a cefepime trough concentration greater than 22 mg/L to predict neurotoxicity

Performance of the Cockcroft-Gault, MDRD, and new CKD-EPI formulas in relation to GFR, age, and body size

Population pharmacokinetics of high-dose, prolonged-infusion cefepime in adult critically ill patients with ventilator-associated pneumonia

Population pharmacokinetic assessment and pharmacodynamic implications of pediatric cefepime dosing for susceptible-dose-dependent organisms

Pharmacokinetics and pharmacodynamics of cefepime in patients with various degrees of renal function

Acknowledgements J.A. Roberts would like to acknowledge funding