Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, November 2006, p. 3680-3688, Vol. 50, No. 11
0066-4804/06/$08.00+0 doi:10.1128/AAC.00369-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
School of Medicine, University of Manchester, Stopford Building, Oxford Road, Manchester M13 9PT, United Kingdom,1 Pediatric Oncology Branch, National Cancer Institute, National Institutes of Health, 10 Center Dr., Bethesda, Maryland 20892,2 Ordway Research Institute, 150 New Scotland Avenue, Albany, New York 12208,3 Wythenshawe Hospital, Southmoor Rd., Manchester, M23 9LT, United Kingdom4
Received 25 March 2006/ Returned for modification 28 May 2006/ Accepted 27 August 2006
|
|
|---|
|
|
|---|
Flucytosine (5FC) is a fluorinated pyrimidine which acts as a subversive substrate within the pyrimidine salvage pathway (16, 17). Flucytosine exhibits antifungal activity against Candida spp. and Cryptococcus neoformans. In clinical contexts, 5FC has been most frequently used in combination with amphotericin B to promote an additive to synergistic interaction (5). Furthermore, the addition of 5FC to amphotericin B has been used as an amphotericin B-sparing strategy and may minimize the probability of the emergence of 5FC resistance. Andes and van Ogtrop characterized the pharmacodynamics of 5FC in a neutropenic murine model of disseminated candidiasis and demonstrated the fraction of the dosing interval that serum levels of 5FC remained above the MIC (T>MIC) was the dynamic variable which best linked drug exposure to the observed drug effect (3).
The primary objective of the current study was the derivation of an in vivo drug exposure breakpoint for 5FC against Candida albicans. We sought to further explore the clinical implications of these findings by using Monte Carlo simulation. The first key step was achieved using a well-characterized murine model of disseminated candidiasis and a single strain of C. albicans. At this point, we explored the extent to which the MIC and the T>MIC accounted for the outcome of 5FC exposure. Subsequently, the exposure-response relationships for an additional five strains of C. albicans were also studied. These strains were chosen both to encompass the distribution of MICs of 5FC against C. albicans and on the basis of a previously defined putative 5FC resistance mechanism (17). This strategy provided an opportunity to investigate the relationship between the 5FC resistance genotype and the outcome following 5FC exposure.
|
|
|---|
64 (
64) mg/liter, respectively. The 5FC MIC was also determined using CLSI (formerly NCCLS) methodology (20), which yielded concordant results, except for strains F/9651 and F/9464, both of which showed an MIC of
64 mg/liter. The 5FC resistance genotype had been previously defined for five of the six strains (17). The genotype of the remaining isolate, F/6862, was defined according to the methodology outlined in a previous publication (17). The relevant details are summarized in Table 1. |
View this table: [in a new window] |
TABLE 1. Study isolates
|
Murine pharmacokinetics of 5FC and determination of the exposure-response relationships for 5FC. The serum pharmacokinetics of 5FC in infected CD1 mice were defined in a previous study (18). The parameter estimates (means ± standard deviations [SD]) were as follows: volume of the central compartment, 0.022 ± 0.017 liter; clearance, 0.021 ± 0.013 liter/h; the intercompartmental rate constants Kcp, 13.440 ± 7.330 h1, and Kpc, 24.070 ± 6.020 h1, and the term describing the first-order absorption of drug from the peritoneal cavity into the central compartment, Ka, 19.580 ± 9.660 h1. The mean parameter estimates from the in vivo pharmacokinetic model (above) were used to transform the dose-response relationships to T>MIC versus effect.
Exposure-response relationships for strain F/6862.
The dose exposure-response relationship for 5FC against F/6862 was determined using 5FC at dosages of 0, 0.5, 1.0, 1.5, 3, 6.25, and 10 mg/kg of body weight administered i.p. every 8 hours for 24 h. A sigmoid Emax model (see equation 1 below) was fitted to the mean fungal densities, and the data were weighted by the inverse of the observed variance at each point. The model was implemented in the ADAPT II package software (11). The sigmoid Emax model was as follows:
![]() | (1) |
Establishing the target to define therapeutic success and failure in mice. The duration of the pharmacodynamic model was extended to 10 days to examine the relationship between the fungal tissue burden induced after 24 h of therapy and the ultimate survival of mice. Cyclophosphamide was readministered every 5 days to reinduce neutropenia (the period of neutropenia extends from day 2 to 6 after cyclophosphamide administration). Four cohorts of mice, consisting of a control group (n = 10) and three treatment groups (n = 10 per group) were infected with strain F/6862. Mice in the treatment groups received 5FC i.p. every 8 h for 24 h only (i.e., only three doses of drug were administered), and animals were then observed. Note that because of the abbreviated therapy, this experiment does not represent a traditional survival model. The control group (cohort 1) received saline (0.4 ml i.p., three doses only), and cohorts 2, 3, and 4 were administered 5FC at 0.18, 0.37, and 1 mg/kg in 0.4-ml injections i.p. every 8 h for 24 h only; these dosage schedules were designed to induce a mean decline in the Candida kidney density of 0.5 log, 1 log, and 2 logs, respectively (Fig. 1). Survival was plotted following a Kaplan-Meier survival analysis.
![]() View larger version (15K): [in a new window] |
FIG. 1. (A) Sigmoid-Emax relationship for isolate F/6862 (MIC, 0.125 mg/liter). Data are means ± SD. The doses producing 0.5-, 1-, and 2-log falls in log10CFU/g were 0.18, 0.37, and 1 mg/kg administered every 8 h, respectively. (B) Kaplan-Meier survival curves for groups of mice infected with strain F/6862 (MIC, 0.125 mg/liter) following 24 h of 5FC therapy, to further examine the consequences of various degrees of antifungal killing. , 2-log drop versus 1-log drop, P = 0.039; , 1-log drop versus control, P = 0.002; , 0.5-log drop versus control, P = 0.854.
|
2-log fall in tissue burden) as a function of drug exposure (fraction of the dosing interval that serum concentrations exceeded the MIC) which developed in each mouse. This analysis was performed with the use of SYSTAT version 11 (Richmond, CA). Population pharmacokinetic data for 5FC in humans and Monte Carlo simulation. The available literature describing the population pharmacokinetics of 5FC was retrieved and reviewed (10, 15, 26); ultimately, data from Ette et al. (15) were chosen. The population parameter mean values and their respective dispersions based upon a 70-kg male with the mean renal function of the population described in the study were used for the simulations. The mean population parameter vectors were embedded in subroutine PRIOR of ADAPT II (11), and a 9,999-subject Monte Carlo simulation was performed. The proportion of simulated patients receiving 25, 50, and 100 mg/kg/day of 5FC in which the serum levels of 5FC exceeded 45% of the dosing interval was determined. The MIC distribution of C. albicans was defined using 1,941 C. albicans isolates, obtained from 10 European data sources and tested using EUCAST methodology (data kindly supplied by Gunnar Kahlmeter, EUCAST).
Effect of genotype, MIC, and T>MIC on the exposure-response relationships. The dose effect, T>MIC-versus-effect relationships, and influence of 5FC resistance genotype on the antifungal response were determined. This was achieved using the same murine model of disseminated candidiasis and the optimized inocula for each strain.
Relationship between in vivo growth of Candida isolates and Emax.
The T>MIC-versus-effect relationships for the six C. albicans isolates did not overlie one another; consequently, it was clear that the T>MIC did not account for all of the residual system variance within the dose-response relationships. Previous studies have suggested that the rate of bacterial growth may be a determinant of the magnitude of the observed effect following exposure to antimicrobial agents. Consequently, we studied the rate of in vivo growth and examined the relationship between the rate of growth and the extent of 5FC-induced killing. These experiments were restricted to four of the six isolates in which a dose-response relationship was evident (F/10141, F/6862, F/7421, and F/9651, with MICs of 0.06, 0.125, 0.5 and 4.0 mg/liter, respectively). A total of nine mice per strain were used. The same inoculum was used for each strain; each mouse received 2 x 104 Candida organisms intravenously in 0.2 ml phosphate-buffered saline. For each isolate, groups of three mice were sequentially sacrificed at 5, 19, and 29 h following the initiation of infection; the fungal density within the kidneys was then determined using the methodology described above. The in vivo growth of each isolate was modeled using equation 2
![]() | (2) |
Incorporation of the rate of in vivo growth into the final model of drug effect.
Since a linear relationship was apparent between Kg and the maximum effect induced by 5FC (see Fig. 6), two competing models to account for drug action were compared. In the first case, drug effect was modeled using a conventional sigmoid-Emax relationship as defined in equation 1. In a second, expanded model represented by equation 3, we allowed the growth rate to independently affect the model estimates of the Candida kidney density following exposure to 5FC. Since Kg appeared to be directly proportional to the maximum effect, the second term in equation 3, which incorporates Kg, is a linear function which contributes a fraction, F, to the overall estimate of drug effect:
![]() | (3) |
![]() View larger version (11K): [in a new window] |
FIG. 6. Correlation between the growth constant (Kg) and the maximum effect induced by 5FC. The maximum effect was taken from the mean parameter estimates shown in Fig. 1.
|
|
View this table: [in a new window] |
TABLE 2. Parameter estimates from the population analysis of the rate of in vivo growth, including the Bayesian estimates for the parameter means and their dispersions describing the in vivo growth of four isolates
|
|
|
|---|
Exposure-response relationships for 5FC. The dose-response relationship for 5FC against strain F/6862 is shown in Fig. 1. The data were accounted for by a sigmoid Emax model with r2 = 98.7%. The doses which resulted in 0.5-, 1-, and 2-log declines in the Candida kidney density at 24 h after treatment initiation were 0.18, 0.37, and 1 mg/kg administered every 8 hours, respectively (Fig. 1 and Table 3).
|
View this table: [in a new window] |
TABLE 3. Magnitudes of drug exposure required to produce various degrees of anti-Candida effect
|
2 = 0.034; 1 df; P = 0.85). In contrast, there was a statistically significant prolongation in the survival of mice in which a 1-log fall had been induced in comparison with controls (
2 = 9.862; 1 df; P = 0.002) as well as in mice in which a 2-log fall had been induced in comparison with 1-log fall (
2 = 4.257; 1 df; P = 0.04). On the basis of these results, a 2-log fall was employed as the microbiological target within the murine model.
Defining the probability of target attainment in vivo.
The logistic regression analysis revealed a statistically significant relationship between the magnitude of drug exposure, quantified in terms of T>MIC, and attainment of the microbiological target of a 2-log fall of CFU/g within the kidney (
2 = 40.13; df = 1; P < 0.001) (Fig. 2). Classification and regression tree analysis revealed that a T>MIC for at least 45% of the dosing interval optimally separated mice attaining the target into high- and low-probability groups; this magnitude of drug exposure corresponded to a 94% probability of attaining the target (Fig. 2) and was employed as the drug exposure breakpoint which was bridged to humans.
![]() View larger version (12K): [in a new window] |
FIG. 2. Relationship between the continuous variable T>MIC and the probability of achieving a 2-log fall in log10CFU/g in the kidneys of mice infected with strain F/6862 (MIC, 0.125 mg/liter), as defined using logistic regression. The regression coefficients for the constant and T>MIC were 7.506 (P = 0.008) and 22.08 (P = 0.004), respectively. The breakpoint value of T>MIC for 45% of the dosing interval was defined using classification and regression tree analysis. A serum level above the MIC for at least 45% of the dosing interval optimally separated the population into groups with high and low probabilities of a 2-log fall within the kidney. Circles represent raw data obtained from eight mice (except for controls and the group in which the T>MIC was 0.83, which consisted of four mice).
|
45%) was attained in at least 95% of simulated patients who were infected with an isolate (of C. albicans) with an MIC of 16 mg/liter or less and who received at least 100 mg/kg/day in four divided doses. This finding suggests that the MIC which defines 5FC resistance, using EUCAST methodology, is 32 mg/liter. Importantly, however, such a breakpoint is tenable only if a human dose of 100 mg/kg/day is used; a downward revision in dose to 50 or 25 mg/kg/day would require a simultaneous reduction in the susceptibility breakpoint to 16 or 8 mg/liter, respectively. Thus, Fig. 3 highlights that dose and susceptibility breakpoints are inextricably linked.
![]() View larger version (18K): [in a new window] |
FIG. 3. Fraction of the 9,999 simulated patients which achieved the drug exposure target (i.e., serum 5FC concentrations greater than the MIC for at least 45% of the dosing interval) following the administration of 5FC at 25, 50 and 100 mg/kg/day, depicted by open circles, triangles, and squares, respectively. The distribution of the MICs for 1,941 C. albicans isolates, obtained using the EUCAST method, is shown by the solid circles.
|
![]() View larger version (26K): [in a new window] |
FIG. 4. Dose-response relationships for each of the study strains. Data are means ± SD. (A) F/10141 (MIC, 0.06 mg/liter; homozygous susceptible); (B) F/6862 (MIC, 0.125 mg/liter; homozygous susceptible); (C) F/7421 (MIC, 0.5 mg/liter; heterozygous); (D) F/9651 (MIC, 4.0 mg/liter; resistance mechanism unknown); (E) F/9464 (MIC, 8.0 mg/liter; homozygous resistant); (F) F/8341 (MIC, 64 mg/liter; homozygous resistant) (for this strain only the maximum possible dose was studied and growth was observed). The model fits from the population analysis and mean parameter estimates for each strain are shown. The EC50s vary 25-fold.
|
![]() View larger version (28K): [in a new window] |
FIG. 5. Relationship between the fraction of the dosing interval that serum 5FC concentrations are above the MIC (T>MIC) and the observed response for each of the study strains. The model was fitted to the data using a population methodology; r2 = 99.2%. The importance of the T>MIC in accounting for the exposure-response relationships is highlighted by the fact that the estimates for the EC50s now vary approximately twofold. The decline in log10CFU/g varies between strains and is independent of the 5FC MIC.
|
The fits of a standard sigmoid Emax model and an expanded model (Fig. 7) which incorporated the growth constant, Kg, to the data were both highly acceptable; in the case of the former, the coefficient of determination for the observed versus predicted values after the Bayesian step was 99.2% and the log-likelihood value was 9.40. The fit of the expanded model, in which the in vivo growth rate was allowed to influence the estimate of drug effect, yielded the same coefficient of determination (r2 = 99.2%), but the model produced a more likely solution (log-likelihood values of 26.92 versus 9.40 for the expanded versus standard models, respectively (twice difference in likelihood value
2 = 35; df = 3; P < 0.001), indicating that the rate of in vivo growth exerted an influence on the observed antifungal effect.
![]() View larger version (19K): [in a new window] |
FIG. 7. Expanded model describing the effect of 5FC. The exposure-response relationships are influenced by both the fraction of the dosing interval that the serum concentrations of 5FC are above the MIC (T>MIC) and the growth rate (Kg).
|
|
|
|---|
In terms of the end point, the most straightforward measure of outcome is a dichotomous variable, such as survival, microbiological clearance, or the resolution of an infectious syndrome (see, for example, methodology employed by Drusano et al. (12, 13) and Andes and Craig (1). In the majority of pharmacodynamic models, however, the biomarker is a continuous variable, and the experimental design is such that sterilization of the affected site is rarely achieved. Consequently, if pharmacodynamic models are to be used for breakpoint definition, a cutoff value must be sought and applied. In the pharmacodynamic literature, two end points have been used: (i) the attainment of stasis (i.e., the reduction of the infectious burden to the level observed at the initiation of therapy) and (ii) some proportional reduction in the infectious burden (e.g., E50 or E90) (see for example references 2 and 6). There is no universally applicable end point; the choice depends on the context. For severe infections, orders of logarithmic killing are probably required; for other infections, a less stringent end point may suffice. In the current study, we attempted to link various orders of logarithmic killing with survival; a 2-log fall induced after 24 h of therapy was associated with prolongation in survival. As shown in Table 3, we could equally have used a 50% reduction in infectious burden and obtained the same result.
The probability with which one wishes to attain the target must also be explicitly defined. If the drug exposure required to achieve a mean 2-log fall was used as the pharmacodynamic target, then only approximately 50% of mice would attain that target; such a response rate would be clinically unacceptable. Ideally one would prefer that a high proportion of the population attain the microbiological target of interest. Using logistic regression, we demonstrated that if serum levels were greater than the MIC for 45% of the dosing interval, there was a 94% probability within the murine population of achieving a 2-log fall in fungal tissue burden. As demonstrated in Fig. 2, progressively higher drug exposures are required to further increase the probability of attaining the target, but the gains are small. Conversely, in circumstances in which drug exposure is suboptimal, a proportion of drug recipients may still respond; i.e., the response rate is not necessarily zero.
The British Society for Medical Mycology (7) and the CLSI (20) have published susceptibility breakpoints for 5FC against Candida spp. According to the British Society for Medical Mycology, the breakpoint for sensitive is
1.0 mg/liter, that for intermediate susceptibility is 2.0 to 8.0 mg/liter, and that for resistant
16 mg/liter (7). The CLSI breakpoints are as follows: susceptible,
4.0 mg/liter; intermediate susceptibility, 8.0 to 16.0 mg/liter; and resistant,
32 mg/liter (20). EUCAST has yet to publish susceptibility breakpoints for 5FC; in this regard, the current study provides the following perspectives that may be useful in this process. First, 5FC resistance could be defined by an MIC of
32 mg/liter, although this is strictly contingent upon a dose of 100 mg/kg/day being used; the use of a smaller dose would require a simultaneous downward revision of the susceptibility breakpoint. Second, if 100 mg/kg/day is to be used, there is no requirement to define an intermediate resistant category; such a construct would be useful only for dosages of 25 to 50 mg/kg/day to indicate that a higher dose would be required to treat certain strains for which therapy would likely fail at the lower dose (e.g., if 25 mg/kg/day was used, isolates with MICs of 8 and 16 could be designated "susceptible dose dependent" to indicate that 100 mg/kg/day would be required to maximize the probability of therapeutic success; there are parallels to draw here with susceptibility breakpoints for fluconazole as defined by the CLSI [20]). Finally, there is an interplay between the effect induced by various simulated doses, as assessed by fractional target attainment rate. As demonstrated in Fig. 4, the fractional population attainment rates for 25, 50, and 100 mg/kg/day are comparable simply because the vast majority of isolates have an MIC of
0.25 mg/liter. These data suggest that the dose of 5FC for the treatment of disseminated candidiasis due to C. albicans could be lowered (with concomitant reduction in the 5FC susceptibility breakpoint, as discussed above). Importantly, however, these results may not apply to non-C. albicans species of Candida, to other yeasts such as Cryptococcus neoformans, or to the treatment of infections within sanctuary sites such as the central nervous system.
A pharmacodynamic approach also provides an opportunity to reflect upon the predictive value of an MIC or other in vitro measures of drug potency, such as the drug resistance genotype. The MICs of 5FC for C. albicans may be difficult to read because of trailing, thus raising the question as to which is the "right" MIC. Clearly, the MIC of strain F/9464, which was consistently estimated at 8.0 mg/liter using EUCAST methodology, is likely to be falsely low given the lack of any observed exposure-response relationship, despite seemingly adequate levels of drug exposure. For this strain, an MIC of
64 as determined using CLSI methodology better predicts the lack of an in vivo response. Conversely, however, an MIC of 4.0 mg/liter for isolate F/9651, as determined using the EUCAST method, does predict the response to 5FC, whereas the estimate of
64 mg/liter produced by the CLSI method does not (falsely high). These results suggest that both methodologies misclassify certain strains, and consequently, further technical refinements are warranted (see, for example, reference 24). The determination of the 5FC resistance genotype may be a useful adjunct to more conventional techniques, especially for strains with reduced susceptibility without frank resistance.
One of the most significant findings of the current analyses is that the rate of in vivo growth appears to be an important and additional determinant of the drug effect. Decreased antimicrobial susceptibility has been described in the setting of reduced growth rates induced by nutrient deficiencies, noxious insults, organisms in a stationary phase of growth, or organisms growing within biofilms (4, 8, 25). The mechanism by which the rate of in vivo growth is related to the antifungal effect is unclear. An explanation which specifically relates to 5FC involves its mechanism of action. Since the period of rapid DNA synthesis occupies a finite time within the cell cycle (9) and is independent of the growth rate, organisms that are rapidly growing spend a higher proportion of their time synthesizing DNA; this may render these strains more susceptible to the action of 5FC.
In summary, experimental models of disseminated candidiasis enable a further understanding of exposure-response relationships and therefore can play a central role in the provision of decision support in the process of setting susceptibility breakpoints. Further research is required to provide a better understanding of the factors related to both the host and fungus which determine the outcome of antimicrobial therapy.
William Hope is supported by an unrestricted educational grant from Merck & Co. This study was supported by a grant from Valeant Pharmaceuticals and the Fungal Research Trust.
Published ahead of print on 5 September 2006. ![]()
|
|
|---|
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»