Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, August 2004, p. 2799-2807, Vol. 48, No. 8
0066-4804/04/$08.00+0 DOI: 10.1128/AAC.48.8.2799-2807.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Julia Chung,2
Susan Liu,2
William Knebel,2 and
Helen Kastrissios2
Cubist Pharmaceuticals, Inc., Lexington, Massachusetts 02421,1 GloboMax Holdings LLC, Hanover, Maryland 210762
Received 9 November 2003/ Returned for modification 26 January 2004/ Accepted 4 April 2004
|
|
|---|
|
|
|---|
1-lactone) is a novel cyclic lipopeptide antibiotic derived from the fermentation of Streptomyces roseosporus. Daptomycin was recently approved for the treatment of complicated skin and skin structure infections (cSSSI) caused by aerobic gram-positive bacteria, including those caused by methicillin-resistant Staphylococcus aureus and methicillin-susceptible S. aureus. In vitro, daptomycin demonstrates a rapid concentration-dependent bactericidal activity against most clinically relevant gram-positive pathogenic bacteria, including bacterial isolates that are resistant to methicillin, vancomycin, and linezolid (9). The MICs for daptomycin at which 90% of isolates tested are inhibited are typically
1 µg/ml for staphylococci and streptococci and 2 to 4 µg/ml for enterococci, including vancomycin-resistant isolates (7). Although the mechanism of action has not been fully defined, it is distinct from those of other antibiotics and appears to be mediated by the disruption of multiple aspects of membrane function (5, 17). In phase 3 trials for the treatment of cSSSI caused by susceptible gram-positive bacteria, clinical and microbiological outcomes of patients treated with daptomycin were comparable to those for patients receiving conventional antibiotic therapy, such as pencillinase-resistant penicillins or vancomycin (1, 18).
Studies of healthy human subjects have demonstrated linear pharmacokinetics after single (20) and multiple (8) intravenous daptomycin doses up to 6 mg/kg of body weight over 14 days. After once-daily doses of 4 mg/kg, the average steady-state trough daptomycin concentration in plasma was 5.89 µg/ml and varied by 27% among individuals (8). Daptomycin is >90% bound to plasma proteins and has a low steady-state volume of distribution, averaging 0.06 to 0.15 liter/kg (8, 19, 20), consistent with distribution into extracellular fluid. Elimination is primarily achieved by renal excretion of unchanged drug. In healthy adult subjects, the mean urinary recovery over a 24-h period is 50 to 60% of the administered dose (8, 19). In phase 1 studies involving healthy adult subjects and subjects with graded renal insufficiencies, including end-stage renal disease, daptomycin plasma clearance (CL) was significantly reduced among subjects with creatinine clearances (CLCR) of
40 ml/min or who were on dialysis (D. A. Sica, T. Gehr, and B. Dvorchik, Abstr. 42nd Intersci. Conf. Antimicrob. Agents Chemother., abstr. 2257, 2002).
For the present analysis, data from 15 clinical trials were combined and a population analysis approach was used to evaluate possible sources of interindividual variability in daptomycin pharmacokinetics. The specific objectives were to develop a model to describe the pharmacokinetics of daptomycin in both healthy volunteers and subjects with acute bacterial infections who were representative of the target patient population and to identify clinical characteristics that impact daptomycin pharmacokinetics. The factors examined included age, sex, weight, and the presence of bacterial infection, end-organ dysfunction, comorbidities, and concomitant medications.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Study designs
|
The covariates explored as possible sources of interindividual variability in daptomycin pharmacokinetics are listed in Table 2. Body surface area was calculated by the following formula (R. D. Mosteller, Letter, N. Engl. J. Med. 317:1098, 1987):
![]() |
![]() |
|
View this table: [in a new window] |
TABLE 2. Subject characteristics
|
(i) Base model selection.
One-, two-, and three-compartment structural models were fit to the data for concentrations in plasma over time; graphical displays of the data were also evaluated. Hypothesis testing to discriminate among alternative hierarchical structural models was performed by using the likelihood ratio test (16). For comparisons of alternative models, the difference in the NONMEM objective function was approximately chi-square distributed, with n degrees of freedom, where n was the difference in the number of parameters between the hierarchical models. A decrease of
3.84 in the value of the NONMEM objective function, which is less than twice the maximum logarithm of the likelihood of the data, is significant in the likelihood ratio test (n = 1; P < 0.05). The goodness of fit was evaluated by using diagnostic scatter plots (not shown).
The duration of infusion (D1) was estimated for a subset of 108 subjects included in the phase 3 clinical trials for whom the date and time of the start of daptomycin infusion and the time of the first pharmacokinetic blood draw for concentration measurements, but not the time of cessation of the infusion, were recorded. All estimates were consistent with the protocols, which specified an infusion time of 30 min. Interindividual variability in D1 could not be estimated and was not included in the model.
All pharmacokinetic parameters were assumed to be logarithmically normally distributed, and exponential interindividual variability terms were included in the pharmacokinetic parameters in the model. Various residual error models were tested, including an evaluation of possible systematic differences between phase 1 and phase 2/3 studies and among studies that used different assay methods.
(ii) Population pharmacokinetic model building. Exploratory analyses were used to guide the model building process. Relationships between individual covariates and Bayesian estimates of the pharmacokinetic parameters were explored graphically. Generalized additive models were used to evaluate both linear and nonlinear relationships between parameters and covariates (14). In addition, measures of body size and renal function markers were tested as possible sources of interindividual variability for each pharmacokinetic parameter. All possible covariate-parameter relationships thus selected were tested, with the exception that possible drug interactions and the effect of the daptomycin dose were examined only for CL and the volume of the central compartment (V1), as appropriate.
Continuous covariates were entered into the population pharmacokinetic model according to the following equation:
![]() |
V is the median value of the covariate in the study population.
1 is the typical value of the parameter (when COV = C
V) and
2 is the slope of the effect of the covariate on the parameter.
Categorical covariates were included in the model by using indicator variables, as shown in the following equation:
![]() |
1 is the typical value of the parameter when the covariate is not present (IND = 0), and
2 is the fractional change in the value of P when the covariate is present (IND = 1). The statistical significance of each covariate-parameter relationship was screened individually in NONMEM, and the model was built by stepwise additions to obtain a full model. Stepwise deletions were used to obtain the final (reduced) model. The likelihood ratio test was used for hypothesis testing to discriminate among alternative hierarchical models. A strict inclusion criterion (P < 0.001) corresponding to a change in the value of the NONMEM objective function of 10.83 (n = 1 degrees of freedom) was used to account for multiple hypothesis testing. At each stage of the analysis, the goodness of fit was evaluated by using diagnostic scatterplots.
(iii) Pharmacokinetic parameter calculations.
The terminal half-life (t1/2), volume of distribution at a steady state (Vss), and area under the curve from time zero to infinity [AUC(0-
)] were calculated from individual pharmacokinetic parameter estimates obtained by Bayesian estimation from the final population pharmacokinetic model. For a two-compartment model (10), the following equations were used:
![]() |
![]() |
![]() |
![]() |
(iv) Statistics.
Individual Bayes estimates and calculated pharmacokinetic parameter values were grouped according to four renal function categories. Three categories were defined by using the estimated CLCR values: the groups were values of
80 ml/min, >40 to <80 ml/min, and
40 ml/min. These ranges were chosen based on an analysis of phase 1 studies with renally impaired subjects (Sica et al., 42nd ICAAC). Subjects on dialysis comprised a fourth category (thus, the categories were CLCR values of
80 ml/min, <80 to >40 ml/min, and
40 ml/min and subjects on dialysis). Differences between groups were evaluated by analysis of variance with Scheffe's test (S-PLUS Professional; Insightful Corp., Seattle, Wash.). P values of <0.05 were considered significant.
|
|
|---|
Imputed covariate values were generated for a total of 37 subjects; laboratory tests of hepatic function were the most frequently missing values. Height, used in the calculation of body surface area, was imputed for 3 subjects; baseline creatinine values in serum were imputed for 13 subjects. Albumin levels in serum were not recorded in one phase 1 study and three phase 2/3 studies and were considered to be missing for all subjects in those studies. Missing body temperature values on the day of pharmacokinetic sampling were not imputed for 29 subjects from the phase 2/3 studies. The population pharmacokinetic model was coded to remove the effect of missing covariates in the model.
Descriptive statistics for all covariates are presented in Table 2.
Pharmacokinetic analysis. (i) Base model. Plots of data for concentrations in plasma versus time (not shown) showed a biphasic disposition of daptomycin. A review of the minimum objective function and diagnostic plots showed that the data for daptomycin concentrations in plasma over time were best described by using a two-compartment open model with first-order elimination. The structural pharmacokinetic model used the following parameters: clearance (CL), the volume of the central compartment (V1), intercompartmental clearance (Q), and the volume of the peripheral compartment (V2). In addition, the duration of infusion (D1) was estimated for several phase 3 subjects. Estimates of D1 were consistent with the duration of infusion specified in the clinical protocols.
The median daptomycin clearance for the study population was estimated to be 0.688 liter/h (11.5 ml/min) and the volume of the central compartment was 4.8 liters. Median estimates for the intercompartmental clearance and volume of distribution of the peripheral compartment were 3.6 liters/h and 3.6 liters, respectively. All pharmacokinetic parameters were precisely estimated, with relative standard errors (RSEs) of <3%. The estimated median duration of infusion for 108 subjects enrolled in phase 2/3 trials was 0.402 h (24 min), with an RSE of 23.8%.
Interindividual variabilities were estimated to be 52.1% for CL, 60.6% for the volume of the central compartment, 31.9% for the volume of the peripheral compartment, and 74.4% for intercompartmental clearance. A simple additive residual error model based on diagnostic plots provided the best fit for the data. A further evaluation of diagnostic plots indicated that there was a larger degree of misfit of model predictions for observations collected in phase 2/3 studies; therefore, different error structures were evaluated for phase 1 versus phase 2/3 studies. On the basis of the likelihood ratio test, the final residual error model was described by a combination of additive errors, reflecting the different assay methods used to determine daptomycin concentrations in plasma. The residual error was slightly lower for the study in which daptomycin concentrations in plasma were assayed by LC/MS/MS than for studies assayed by HPLC (2.08 versus 4.72 µg/ml, respectively), consistent with the higher sensitivity of the former method.
(ii) Population model. Exploratory graphical analyses revealed a direct correlation between daptomycin clearance and various markers of renal function, including estimated creatinine clearance, renal function category, and laboratory markers of renal function. Intercompartmental clearance (in liters per hour) and the volume of the peripheral compartment (in liters) were correlated with body weight. There were no obvious relationships between V1 and any of the tested covariates and no significant association between daptomycin pharmacokinetics and either the concomitant medications or the concomitant diseases evaluated in this patient population.
The final model for daptomycin clearance was determined to be the following:
![]() |
![]() |
Both intercompartmental clearance and the volume of the peripheral compartment were determined to be functions of body weight, as follows:
![]() |
![]() |
The median value for V1 was estimated to be 4.80 liters. Although the interindividual variability in V1 was 57%, none of the covariates investigated, including body weight, was identified as a significant source of variability in V1.
Additional exploratory analyses were performed to evaluate whether any other covariate could explain the effect of sex on daptomycin clearance or the effect of the presence of infection on V2. A review of the covariate graphics indicated that body weight, body surface area, and race differed by sex. Each of these covariates was substituted into the clearance model to determine if it could be substituted for sex in the model, but none produced a significant change in the objective function value. Consequently, the clearance model including sex represented the final model.
Infections were only present in subjects in the phase 2 and/or 3 clinical trials. Therefore, in the model the presence of infection could be a marker for an unmonitored covariate that differed between the phase 1 and phase 2/3 clinical trials. In a graphical evaluation, age and serum albumin were determined to differ between the two groups. Each of these covariates, as well as body temperature, was substituted into the V2 model, and none produced a significant change in the objective function value. The V2 model including infection represented the final model.
Parameter estimates for the final population pharmacokinetic model are presented in Table 3. Pharmacokinetic parameters were precisely estimated and diagnostic plots showed a good fit of the final model to the observed daptomycin concentrations in plasma (Fig. 1).
|
View this table: [in a new window] |
TABLE 3. Population pharmacokinetic parameter estimates for daptomycinb
|
![]() View larger version (41K): [in a new window] |
FIG. 1. Diagnostic plots for daptomycin population pharmacokinetic model. Observed versus predicted daptomycin concentrations in plasma (left), observed versus individual predicted daptomycin concentrations in plasma (middle), and weighted residuals versus predicted daptomycin concentrations in plasma (right panel) are shown. Circles represent individual data points. Dashed lines represent regression lines. Solid lines represent unity.
|
Individual estimates of CL and V1 were obtained from the final population pharmacokinetic model by Bayesian estimation. These were used to calculate individual estimates of the t1/2, Vss, and AUC(0-
) for a single 4-mg/kg intravenous dose from the individual parameter estimates and were summarized by renal function category. Summary statistics for these estimates are presented for all subjects (Table 4) and separately for phase 1 and phase 2/3 subjects (Table 5).
|
View this table: [in a new window] |
TABLE 4. Summary of pharmacokinetic parameters sorted by estimated CLCR and obtained by Bayesian estimation from the final model
|
|
View this table: [in a new window] |
TABLE 5. Summary of pharmacokinetic parameters by study phase and estimated CLCR of individual parameters obtained by Bayesian estimation from the final model
|
) for a single 4-mg/kg intravenous dose were dependent on renal function. An analysis of variance indicated that compared with subjects with CLCR values of >40 ml/min, subjects whose CLCR values were
40 ml/min or who were on dialysis had significantly larger volumes of the central compartment. This factor was not significant in the population pharmacokinetic analysis, most likely because of the large variation in V1 and the relatively small number of subjects on dialysis. Vss was not dependent on renal function.
Relative to that in subjects with normal renal function (CLCR values of
80 ml/min), the daptomycin half-life was increased 2.3-fold in subjects with CLCR values of
40 ml/min and 3.5-fold in subjects who were on dialysis; changes in dose-normalized AUC(0-
) values were 1.8-fold and 3-fold, respectively (Table 5). In comparison, median half-life and dose-normalized AUC(0-
) values in subjects with CLCR values of
80 ml/min and in subjects with CLCR values of <80 and >40 ml/min differed <10%. These differences, although statistically significant, were not considered clinically meaningful.
|
|
|---|
This report represents the first population pharmacokinetic analysis of daptomycin and includes subjects from all three phases of the clinical development program. The increased use of population pharmacokinetic analysis has generated a number of newer software programs that, like NONMEM, have advantages and disadvantages. One topic of discussion has been the ability and ease of use of NONMEM to detect the presence of a nonnormal distribution, especially within a subpopulation. The intelligent use of any pharmacokinetic program is a prerequisite for a meaningful analysis. NONMEM, when used by trained personnel, does allow one to detect a distribution that is substantially nonnormal.
Among healthy subjects, the estimated pharmacokinetics were consistent with those previously reported for a phase 1 study in which single doses of 0.5 to 6 mg of daptomycin/kg were administered intravenously to healthy volunteers (20). The population analysis defined quantitatively the decrease in daptomycin CL associated with reduced renal function, a relationship that was suggested by earlier phase 1 studies. New findings included the increase in the volume of the peripheral compartment (V2) in subjects with bacterial infections relative to healthy subjects as well as the associations between weight and both intercompartmental clearance (Q) and V2.
Renal function, sex, and body temperature accounted for 21.5% of the interindividual variability in daptomycin clearance, with renal function being the single most significant explanatory variable. During the screening of covariates, adding just renal function (i.e., CLCR in nondialysis subjects and a flag for subjects on dialysis) to the clearance model reduced the interindividual variability by 18.9%, from 52.1% in the base model (no covariates) to 33.2% (with the addition of renal function markers) (data not shown). This finding is consistent with the fact that daptomycin, like other hydrophilic antibiotics, is cleared primarily by renal excretion (20).
The median estimated daptomycin clearance for a normothermic male with an estimated CLCR of 91.2 ml/min was 0.807 liter/h (13.5 ml/min). Among subjects on dialysis, the median daptomycin clearance was estimated to be 0.269 liter/h (4.5 ml/min), or approximately one-third that of nondialysis subjects. Among subjects who were not on dialysis, daptomycin clearance was a linear function of CLCR. For example, for an increase or decrease in the estimated CLCR of 10 ml/min, the daptomycin clearance increased or decreased by 0.05 liter/h (0.8 ml/min).
The dose of daptomycin recommended for the treatment of cSSSI is 4 mg/kg administered by intravenous infusion once every 24 h for subjects with CLCR of
30 ml/min and once every 48 h for subjects with lower CLCR values, including those who are on dialysis. These recommendations are based on several observations in addition to the data presented in this report. Sica et al. (42nd ICAAC) determined mean Cmax and AUCss values among 44 subjects with graded renal impairment or who were undergoing dialysis. The Cmax was consistent for all subjects and the AUCss was similar in all subjects who had an estimated CLCR of >40 ml/min. For subjects with an estimated CLCR of
40 ml/min, the AUC was increased 2.33-fold compared to that for subjects with a CLCR of >80 ml/min. The two phase 3 trials for the treatment of cSSSI included a limited number of subjects with estimated CLCR values between 30 and 40 ml/min. There was no increase in adverse events attributed to daptomycin among these subjects; none participated in the pharmacokinetic studies reported here. Additional studies of the pharmacokinetics and safety of daptomycin in renally impaired subjects and in those undergoing dialysis are in progress.
Daptomycin clearance was influenced to a lesser extent by sex and body temperature. Clearance in females was estimated to be approximately 80% that of male subjects with similar renal function. Among subjects with cSSSI treated with daptomycin in two recent large phase 3 trials, there were no clinically or statistically significant differences between the success rates for males (n = 230) and females (n = 192) (74.8 versus 77.1%, respectively; 95% confidence intervals, 7.2 and 8.3) (data on file, Cubist Pharmaceuticals). Thus, although the difference in daptomycin clearance related to sex was statistically significant, it does not appear to be clinically meaningful.
The observation that daptomycin clearance increased with elevated body temperatures (>37.2°C) should be interpreted cautiously since the analysis was limited to data obtained from 100 subjects in the phase 2/3 clinical studies, of whom only 14% were hyperthermic (body temperature of
38°C).
Comorbidities, including diseases producing fluid accumulation (e.g., ascites and edema), diabetes, hypertension, and congestive heart failure, were not significantly correlated with daptomycin clearance. Medications that were tested for possible pharmacokinetic interactions with daptomycin included acidic drugs that are actively secreted in the renal tubule and drugs that are highly (>95%) bound to albumin (3, 13, 15). These had no effect on daptomycin pharmacokinetics.
The estimated increases in Q and V2 for daptomycin with increased body weights were consistent with the physicochemical properties of daptomycin and the physiologic effects of weight. Daptomycin appears to be restricted to the extracellular space which increases with body weight (15). Similarly, the extravascular distribution of daptomycin occurs via diffusion (15), which would also be facilitated by the increased fluid (water) associated with an increased body weight.
V2 was estimated to be approximately twofold larger in subjects with acute bacterial infections than in uninfected subjects. This is consistent with the pathophysiology of acute bacterial infections, which is characterized by an inflammatory response associated with increased vascular permeability and the collection of extracellular fluid at the site of infection. However, since bacterial infections were only present among subjects in the phase 2/3 clinical trials, it is also possible that this factor was a surrogate for another, possibly unmonitored, covariate or an unidentified systematic difference between the phase 1 and phase 2/3 clinical trials. Currently, there is no recommendation regarding increased doses of daptomycin for patients with exceptionally severe infections or impaired host defenses. A trial of daptomycin at 6 mg/kg intravenously once a day for the treatment of infective endocarditis due to S. aureus is in progress.
Although there was an appreciable variability in the estimates of V1, none of the covariates investigated was identified as a significant source of this variability. In a recently reported phase 1 study of daptomycin pharmacokinetics using subjects with graded renal insufficiencies and end-stage renal disease (Sica et al., 42nd ICAAC), the total volume of distribution for daptomycin was increased among subjects on hemodialysis. In the present larger study, this relationship was extended and further defined as an increase in V1 in subjects with CLCR values of
40 ml/min as well as in subjects on dialysis, but only in a supplemental analysis of variance (Table 5). This may be because of the relatively small proportion of subjects who were undergoing dialysis or perhaps because the effect represents another, possibly unmonitored, covariate. Additional studies of daptomycin pharmacokinetics in subjects on hemodialysis are in progress.
In conclusion, this population analysis of daptomycin pharmacokinetics indicates that renal function is the single most significant factor contributing to interindividual variabilities in daptomycin clearance. Because of their reduced daptomycin clearance, patients on dialysis and those with severe renal disease (CLCR of <30 ml/min) will require adjusted dosage regimens to achieve systemic exposures that are clinically and pharmacologically comparable to those seen in subjects with higher levels of renal function. Daptomycin clearance was also impacted by sex and body temperature. However, an analysis of clinical outcomes suggested that the variation associated with sex is not clinically meaningful. The relationship with body temperature should be interpreted cautiously since the analysis was limited to the subset of subjects from phase 2/3 clinical studies, of which only 14% were hyperthermic. The relationships between body weight and the rate and extent of extravascular distribution support the dosing of daptomycin on the basis of milligrams per kilogram of body weight.
We acknowledge the cooperation and assistance of the subjects, investigators, and study personnel who participated in these trials and the support of our colleagues in the Cubist Clinical Department.
Present address: Paratek Pharmaceuticals, Boston, MA 02111. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»