This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hurwitz, S. J.
Right arrow Articles by Schinazi, R. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hurwitz, S. J.
Right arrow Articles by Schinazi, R. F.

 Previous Article  |  Next Article 

Antimicrobial Agents and Chemotherapy, December 2008, p. 4241-4250, Vol. 52, No. 12
0066-4804/08/$08.00+0     doi:10.1128/AAC.00054-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Development of an Optimized Dose for Coformulation of Zidovudine with Drugs That Select for the K65R Mutation Using a Population Pharmacokinetic and Enzyme Kinetic Simulation Model{triangledown}

Selwyn J. Hurwitz,1,2 Ghazia Asif,1,2 Nancy M. Kivel,3 and Raymond F. Schinazi1,2*

Center for AIDS Research, Laboratory of Biochemical Pharmacology, Department of Pediatrics, Emory University School of Medicine, Atlanta, Georgia 30322,1 Veterans Affairs Medical Center, Decatur, Georgia 30033,2 RFS Pharma, LLC, Tucker, GA 300863

Received 14 January 2008/ Returned for modification 5 April 2008/ Accepted 28 September 2008


arrow
ABSTRACT
 
In vitro selection studies and data from large genotype databases from clinical studies have demonstrated that tenofovir disoproxil fumarate and abacavir sulfate select for the K65R mutation in the human immunodeficiency virus type 1 polymerase region. Furthermore, other novel non-thymine nucleoside reverse transcriptase (RT) inhibitors also select for this mutation in vitro. Studies performed in vitro and in humans suggest that viruses containing the K65R mutation remained susceptible to zidovudine (ZDV) and other thymine nucleoside antiretroviral agents. Therefore, ZDV could be coformulated with these agents as a "resistance repellent" agent for the K65R mutation. The approved ZDV oral dose is 300 mg twice a day (b.i.d.) and is commonly associated with bone marrow toxicity thought to be secondary to ZDV-5'-monophosphate (ZDV-MP) accumulation. A simulation study was performed in silico to optimize the ZDV dose for b.i.d. administration with K65R-selecting antiretroviral agents in virtual subjects using the population pharmacokinetic and cellular enzyme kinetic parameters of ZDV. These simulations predicted that a reduction in the ZDV dose from 300 to 200 mg b.i.d. should produce similar amounts of ZDV-5'-triphosphate (ZDV-TP) associated with antiviral efficacy (>97% overlap) and reduced plasma ZDV and cellular amounts of ZDV-MP associated with toxicity. The simulations also predicted reduced peak and trough amounts of cellular ZDV-TP after treatment with 600 mg ZDV once a day (q.d.) rather than 300 or 200 mg ZDV b.i.d., indicating that q.d. dosing with ZDV should be avoided. These in silico predictions suggest that 200 mg ZDV b.i.d. is an efficacious and safe dose that could delay the emergence of the K65R mutation.


arrow
INTRODUCTION
 
Current first-line highly active antiretroviral therapy (HAART) for the treatment of human immunodeficiency virus (HIV-1) infections combines two nucleoside reverse transcriptase inhibitors (NRTI) together with either a protease inhibitor or a non-NRTI (18, 19, 54). These drug combinations have markedly decreased mortality and morbidity from HIV-1 infections in the developed world (11). Existing therapeutic modalities cannot eradicate HIV-1 infection because of the compartmentalization of the virus and its latent properties (80, 81). Therefore, chronic therapy remains the standard of care for the foreseeable future. HAART regimens are selected in part to minimize cross-resistance and thereby delay the emergence of resistant viruses. However, all regimens eventually fail, due primarily to a lack of adherence to strict regimens, delayed toxicities, and/or the emergence of drug-resistant HIV-1 strains (68). Thus, it is a major imperative to develop regimens that delay, prevent, or attenuate the onset of resistance for second-line treatments for infected individuals who have already demonstrated mutations in the systemic circulation. The occurrence of common resistance mutations, including thymine analog mutations (specifically, D67N, K70R, T215Y, and T219Q), K65R, and M184V, needs to be the continued focus in the drug development of HIV-1 NRTI (79).

Data from large genotype databases demonstrated an increased incidence of the K65R mutation from 0.8% in 1998 to 3.8% in 2003, presumably as a result of the increased use of K65R-selecting drugs (33, 76). This mutation produces a single amino acid change from lysine to arginine in the HIV-1 reverse transcriptase (RT) gene. The in vitro selection of K65R, accompanied with moderate resistance, has been demonstrated for nonthymine NRTI, including abacavir sulfate (ABC), tenofovir disoproxil fumarate (TDF), zalcitabine, didanosine, adefovir dipivoxil and lamivudine (3TC), β-D-2',3'-didehydro-2',3'-dideoxy-5-fluorocytidine (D-d4FC, dexelvucitabine; Reverset), and β-D-(2R,4R)-1,3-dioxolane guanosine (DXG) (39, 76, 88, 89). The guanosine nucleoside prodrug of DXG, amdoxovir [AMDX; (–)-β-D-2,6-diaminopurine dioxolane (DAPD)] (12, 25, 31), is being developed by RFS Pharma, LLC, primarily for the second-line treatment of HIV-1 infections (35, 52). The advantages of DXG include an increased sensitivity to M184V/I strains in vitro and activity against thymine analog mutations that may have been selected during previous antiretroviral therapy and 69SS double insert (36, 37, 55). DXG is synergistic with several NRTI, including zidovudine (ZDV), 3TC, and nevirapine (36). In vitro studies with HIV-1 in culture with MT-2 cells demonstrated a slow onset of resistance to DXG that was associated with the K65R mutation (66, 67, 89). An in vitro study demonstrated that ZDV alone selected for a mixture of K70K/R mutations at week 25, and DAPD alone selected for a mixture of the K65R and L74V mutations at week 20. However, when DAPD and ZDV were used in combination in HIV-infected primary human lymphocytes, no drug-resistant mutations were detected through week 28 (71). Furthermore, mechanistic studies demonstrated that K65R mutants remained susceptible to thymine NRTI, including ZDV and stavudine (8, 29, 38, 39, 66). Regimens containing TDF in combination with lamivudine and abacavir have demonstrated high failure rates due to the emergence of drug-resistance mutations, including K65R. A composite analysis of data from those regimens revealed a success rate of 86% when the regimen contained ZDV compared to 62% when it did not. Furthermore, no K65R mutations were observed in subjects on regimens containing ZDV, suggesting that ZDV has value in preventing selection for this mutation (33). A previous study has also demonstrated a reduced selection of the K65R mutation when ZDV was added to a regimen containing ABC (50). Therefore, the coformulation of ZDV with K65R-mutation-selecting drugs is warranted (33, 53, 69), since ZDV has the potential to serve as a "resistance repellent" agent for the K65R mutation. The addition of ZDV may not be appropriate if it competes for rate-limiting enzyme phosphorylation with other NRTI in the HAART regimen. However, the enzyme used for the rate-limiting phosphorylation step of ZDV differs from those of TDF, ABC, and DXG (3, 4, 20, 22-25, 30, 41, 55, 60, 85).

ZDV was the first antiretroviral drug tested in the clinic initially as a monotherapy drug and later as a component of HAART regimens (11, 17, 26) and was approved as a generic formulation in September 2005 by the United States Food and Drug Administration. Like other NRTI, ZDV undergoes three intracellular phosphorylation steps to form the active ZDV-5'-triphosphate (ZDV-TP), which inhibits wild-type HIV-1 RT with a median inhibitory concentration of about 0.035 µM (72). The single-dose plasma pharmacokinetics of ZDV have been well described in HIV-1-infected individuals following intravenous and oral administration (1, 10, 20, 34, 65, 78, 90).

ZDV treatment is limited by toxic side effects, including nausea and malaise, as well as serious bone marrow cytotoxicities, including anemia and neutropenia (15, 73, 82). The bone marrow cytotoxicities of ZDV are believed to be associated with mitochondrial damage and correlate with the ZDV-5'-monophosphate (ZDV-MP) content (87). The current approved oral dose of ZDV in most of the world is 300 mg twice a day (b.i.d.). A double-blind, parallel group multicenter study involving 474 HIV-infected patients, comparing ZDV monotherapy daily doses of 400, 800, and 1,600 mg, was published in 1992 (62). A dose-dependent statistically significant increase in the incidence of anemia and leucopenia (P = 0.008), neutropenia (P = 0.0005), and neuropathy (P = 0.03) was observed, and the percentage of subjects who failed to complete the study due to side effects was dose related (21%, 31%, and 32% for the 400-, 800-, and 1,600-mg daily doses of ZDV, respectively). Although there was a trend toward fewer cases of AIDS dementia complex, it was not statistically significant. It was concluded that lower ZDV doses reduced toxicity, and doses >400 mg to 600 mg a day offered no clinical advantage (62). Furthermore, a study by Barry et al. (7) reported that a reduced dose of 100 mg ZDV three times a day (t.i.d.) produced similar amounts of cellular ZDV-TP, which mediates its antiviral effect, with significantly decreased ZDV plasma concentrations and intracellular amounts of ZDV-MP, lending support to enzymatic studies that suggest thymidylate kinase (TMPK) may be oversaturated at clinical doses (30). There are conflicting reports regarding the clinical relevance of the saturation of TMPK at clinical ZDV doses. Fletcher et al. reported higher average amounts of cellular ZDV-TP and decreased variance (0.62 nM and 32% coefficient of variation [CV] versus 0.76 nM and 16% CV, respectively) when ZDV doses were adjusted to maintain a target plasma concentration (27), suggesting that variability between subjects in pharmacokinetics is important. However, the accumulation of ZDV-TP is further complicated, since phosphorylation is cell cycle dependent and the fraction of dividing cells may vary between HIV-infected individuals (48, 49, 59, 86). Once activated, cells from different hosts may have various TMPK activities (49). Given the large variation in cellular ZDV-TP concentrations measured in peripheral blood mononuclear (PBM) cells (a CV of approximately 100%) (2, 28), large clinical trials would be needed to statistically demonstrate whether ZDV phosphorylation is saturable. A potentially useful approach could be to merge the enzyme kinetic data, which are best characterized in vitro, with the in vivo-derived population pharmacokinetic parameters of ZDV in HIV-infected individuals, in an in silico simulation study. This approach would allow large populations of virtual subjects to be simulated and compared with actual clinical data.

The results of Barry et al. (7) suggest the potential for dose saturation of ZDV phosphorylation; moreover, in vitro data indicate that virus containing the K65R mutation is more sensitive to ZDV (8, 50, 66, 67). Furthermore, a regimen of 600 mg ZDV once a day (q.d.) was recently shown to produce a less-pronounced and slower onset of viral depletion than was observed for the standard regimen of 300 mg b.i.d. (77). Therefore, the objectives of this study were to develop a population pharmacokinetic and pharmacodynamic model that combines population pharmacokinetic parameters and population variance in cellular enzyme levels in HIV-1-infected individuals, to determine whether dose saturation of ZDV-TP is supported mechanistically, to develop an optimal dosage regimen of ZDV for coformulation as a K65R-resistant repellent with DAPD and other NRTI when the initial load of K65R is low, and to assess whether the outcomes of the 600 mg ZDV q.d. trial were predictable based on plasma pharmacokinetics and enzyme dynamic data derived in vitro. Since ZDV reduces the overall viral load by <1 log, it was considered prudent to target the cellular contents of ZDV-TP similar to those observed presently in the clinic.


arrow
MATERIALS AND METHODS
 
The structure of the population pharmacokinetics and the cellular pharmacology models used in this study are summarized in Fig. 1.


Figure 1
View larger version (42K):
[in this window]
[in a new window]

 
FIG. 1. Overall schematic representation of the pharmacokinetic and viral dynamic model of ZDV phosphorylation.

Population pharmacokinetics of ZDV. Population characteristics and pharmacokinetic parameters are summarized in Table 1. The two-compartment model of Zhou et al. (90) was fitted using data from 175 individuals who received 200 mg ZDV t.i.d. However, only 93 subjects were used to model ZDV in the absence of nevirapine, which produced a pharmacokinetic interaction with ZDV in that study. Drug absorption was approximated as zero order, since insufficient data were available to characterize the absorption phase, and subjects were divided into groups of fast and slow absorbers. Fast absorbers (41.7% of individuals) absorbed ZDV over 0.25 h, while slow absorbers (the remainder) absorbed ZDV over 1.57 h. The systemic clearance (CL/F; liters/h) was related to the covariate's age (years) and body weight (kg) using the following equations. For individuals <30 years old, the equation CL/F = 127 + 0.93 x (body weight – 70) was used; otherwise, the equation CL/F = 127 + 0.93 x (body weight – 70) + 6.52 x (age – 25) was used. The steady-state volume of distribution (Vss/F; liters) was related to body weight using the equation Vss/F = 464 + 9.83 x (body weight – 70). Similarly, the volume of the central compartment (V1/F; liters) is calculated using a constant fitted ratio of Formula using the equation Formula. The volume of the peripheral tissue compartment (V2/F) was calculated by using the ratio RV1/V2 (Table 1).


View this table:
[in this window]
[in a new window]

 
TABLE 1. Population demographics and pharmacokinetic parameters of oral ZDV used for the simulation study

Both CL/F and Vss/F were assumed to follow log-normal distributions, with variance equal to {omega}2CL/F and {phi} x {omega}2CL/F, respectively, where {phi} is the fixed ratio between variances. The residual/interindividual variance was modeled as {varepsilon}ij x b x Cp, where {varepsilon}ij (the error term) is normally distributed with a mean of 0 and a variance of {sigma}2, the factor b is a constant parameter that is assumed equal to 1 after the absorption phase (>2 h) and allows for a greater variability during the drug absorption phase, and Cp is the plasma concentration of ZDV. Given the limited number of subjects, large number of covariates, and lack of a reported variance-covariance matrix, the available data might not have supported the successful estimation of all potentially important covariate effects and covariances. It is also possible that linking intersubject variances in CL/F and V/F could result in a reduced overall variability in the simulation. The covariate distributions from Zhou et al. (90) were used unchanged in the simulation.

The plasma pharmacokinetics of ZDV are characterized by a biexponential decay after intravenous injection (9). However, some pharmacokinetic studies have utilized a one-compartment model to describe ZDV disposition (65), in part due to some obscuring of the distribution phase by the absorption phase. Panhard et al. fitted a one-compartment model in 10 infected individuals receiving ZDV alone (65). The residual/interindividual variance was described using the combined proportional and additive error model {sigma}2total = {sigma}2 x (a + Cij), where a is a fitted constant and Cij is the predicted plasma concentration. Since much of the variance occurs during drug absorption, this approach may tend to overestimate variance in the postabsorption phase, resulting in an overpredicted variability of ZDV-MP and -TP. The parameters for both pharmacokinetic models are presented in Table 1.

Accumulation of ZDV and its nucleotides in human PBM cells. Plasma and cytosolic concentrations of ZDV were assumed to equilibrate instantaneously, since ZDV is not appreciably protein bound and equilibration between the plasma and cytoplasm is achieved rapidly due to the action of equilibrative nucleoside transporters present on the cell membranes of lymphocytes (5, 70). The initial intracellular phosphorylation step of ZDV is catalyzed by cellular thymidine kinase (TK). The primary enzyme for 5'-monophosphorylation to ZDV-MP is TK1, which is located primarily in the cytosol of cells in the S phase. However, mitochondrial TK (TK2) has also been shown to phosphorylate ZDV in cultured monocytes/macrophages that do not express TK1 but to a much lesser degree (3, 4, 22, 23, 30, 60). The Km (concentration at 50% of the maximal metabolism rate) of ZDV versus TK1 is 3 µM (23, 30), and the activity of TK1 versus ZDV is 0.6 of the activity versus thymidine. TMPK catalyzes the subsequent phosphorylation to ZDV-5'-diphosphate (ZDV-DP) and is rate limiting with a Km of 7.6 to 8 µM versus that of ZDV-MP (30, 45). Distributions describing the maximal rates of conversion of the phosphorylation of ZDV (Formula) and ZDV-MP (Vmax,TMP), respectively, were estimated using data from previous enzymatic studies conducted on PBM cells isolated from cohorts of infected individuals (48, 49). These and other studies have reported decreased average levels of TK1 and TMPK in the in vitro-activated cells of individuals treated with ZDV for more than 6 months. Therefore, the distribution of TK1 and TMPK activities was calculated using the TK1 activities of 24 HIV-positive individuals not previously exposed to AZT (48) and 27 HIV-positive individuals with various exposures to ZDV (49), respectively. Since TK1 and TMPK are cell cycle dependent (32), the PBM cells of infected individuals were stimulated ex vivo using phytohemagglutinin (PHA) (48, 49), followed by lysing and measurement of the enzyme activities. The activities of TK1 were measured using ZDV (20 µM), which is a higher concentration than that used in vivo but is sufficient to saturate the enzyme, allowing the capacity of the enzyme (Formula) to be measured. TMPK activities were measured using thymidine-MP as a substrate (50 µM). Vmax,TMP was calculated by assuming a TMPK efficiency of 1% of that of ZDV-MP relative to thymidine-MP (45). The relatively slow Vmax of ZDV-MP is believed to be related to steric hindrance in the binding of ZDV-MP to TMPK (51). Vmax,TMP and Formula were assumed to be distributed in a log-normal manner. To ensure that Vmax values remained in the physiological range, maximal values were constrained to be less than the largest values reported for noninfected individuals in the study, since the values were higher in noninfected than in HIV-infected individuals (48, 49). Units of the apparent Vmax were converted to µmol of ZDV-DP/h/liter, using an average of two reported measurements of the protein content of PHA-stimulated PBM cells (0.053 mg/106 cells) (4, 43). The aqueous volume of activated lymphocytes was calculated using a mean projected cell surface area of 80 µm2 (84), assuming a specific gravity of 1.06, 70% water, and spherical geometry. This volume was then used to convert pmol/106 cells to µM. The parameters used for modeling the cellular phosphorylation and dephosphorylation of ZDV are presented in Table 2. It was assumed that PHA stimulates about 40% of the lymphocyte population in vitro (86).


View this table:
[in this window]
[in a new window]

 
TABLE 2. Population enzyme kinetic parameters for cellular metabolism of ZDV used for the simulation study

The cytoplasmic accumulation of ZDV-MP is described by the equation Formula, where C1 is the concentration of ZDV in cells (at equilibrium, C1 = Cp), kDP,MP is the rate of dephosphorylation of ZDV-DP to ZDV-MP, kMP is the rate of dephosphorylation of ZDV-MP to ZDV, and Vmax,TMP is the maximal rate of conversion of the phosphorylation of ZDV-MP to DP.

The final phosphorylation step to active ZDV-TP is catalyzed by nucleoside diphosphate kinase, is not rate limiting under physiological conditions, and takes place too rapidly to be easily characterized in vitro (86). Considering the area-under-the-curve ratio of ZDV-TP to ZDV-MP in vivo (0.42 ± 0.42:0.52 ± 0.32 [pmol/106 cells/h ± the standard error] for a 100-mg dose and 0.61 ± 0.81: 0.56 ± 0.57 for a 300-mg dose, respectively [7]), a constant 1:1 ratio of ZDV-TP to -DP was assumed. kDP,MP and kMP were calculated using steady-state data from an in-vitro study of ZDV nucleotides measured in PHA-stimulated PBM cells following a 4-h incubation with various concentrations of ZDV (86). Thus, at steady state, the rate of accumulation of ZDV-MP [Formula] equals the rate of decomposition of ZDV-MP to ZDV (ZDV-MP x kMP), allowing kMP to be calculated as the regression slope of Formula versus ZDV. Likewise, the first-order decomposition rate of ZDV-DP to ZDV-MP was calculated as the regression slope of Vmax,TMP x ZDV-MP/(ZDV-MP + Km,TMP) versus ZDV-DP (86).

The accumulation of ZDV-DP and ZDV-TP (ZDV-DP,TP) was described by the equation dZDV-DP,TP/dt = Vmax,TMP x ZDV-MP/(ZDV-MP + Km,TMP) – ZDV-DP x kDP,MP, where d(ZDV-DP,TP)/dt and is the rate of change in the total cellular concentration of ZDV-DP and -TP.

The cellular contents of the various nucleotides of ZDV in activated cells were used to estimate the ZDV nucleotide contents of cells in vivo, assuming an in-vivo-activated cell fraction of 8% (40) and that phosphorylation in active cells is minimal in nonactivated cells (86). It was previously shown that the ratio between ZDV-MP, -DP, and -TP was similar in resting and PHA-stimulated PBM cells, while the absolute concentrations were 60 to 150 times higher in stimulated cells (86). Thus, assuming an 8% stimulation of PBM cells in vitro, activated PBM cells would contribute between 84 and 92% of the total ZDV nucleotides in humans. The average percentage (88%) was used to calculate the total amounts of ZDV-MP, -DP, and -TP using the simulated amounts in activated cells. Plasma concentrations of ZDV and cellular amounts of ZDV-MP were compared with data from existing clinical studies (7, 28, 65, 90).

Computer simulations. Monte Carlo population pharmacokinetic and enzyme dynamic simulations were conducted for 5,000 individuals per simulated dose regimen, using Trial Simulator version 2.2.1 (Pharsight Corp., Mountain View, CA), which utilizes a fifth-order Runga-Cutta algorithm for numerical integration. This program allows customized differential equations, together with probability distributions for each parameter in the equations, to be entered. The simulated results were then analyzed using routines in the S-Plus computer program (version 6.0 professional; Insightful Corp., Seattle, WA) embedded in the Trial Simulator software. The medians, 25th and 75th percentiles versus time of ZDV concentrations in plasma, and cytoplasmic content of ZDV-MP and ZDV-TP were calculated. Three simulations of 5,000 individuals each were performed to ensure reproducibility of the output (<5% differences in medium and interquartile ranges).


arrow
RESULTS
 
The outline of the model is shown in graphic form in Fig. 1. Simulations were performed using both pharmacokinetic models at ZDV doses of 50, 100, 200, and 300 mg b.i.d. and 600 mg q.d. Median and interquartile ranges of plasma ZDV and cellular ZDV-MP and -TP were plotted for the 200 and 300 mg b.i.d. doses (Fig. 2 to 4). Box plots of between-dose peak and trough amounts of ZDV-TP (pmol/106 cells) (median and interquartile ranges) for all doses using the one- and two-compartment population pharmacokinetic models are shown in Fig. 5A and B, respectively. The levels of ZDV and ZDV-MP, -DP, and -TP essentially achieved steady state within 48 h. Plasma concentrations were assumed to increase linearly with dosage in both models. The two-compartment model (90) produced higher maximum concentration (Cmax) values, a more gradual terminal slope, and a smaller variance after the absorption phase (t > 2 h) than the one-compartment model (Fig. 2). Likewise, peak amounts of cellular ZDV-MP were higher and terminal decline slopes were slower for the two-compartment model than the one-compartment model (Fig. 3). Both models predicted a noticeable increase in predicted ZDV plasma concentrations and amounts of ZDV-MP when doses were increased from 200 to 300 mg b.i.d.


Figure 2
View larger version (28K):
[in this window]
[in a new window]

 
FIG. 2. Simulated plasma concentrations of ZDV (n = 5,000) for the 200- and 300-mg doses (gray and black, respectively) according to the models of Zhou et al. (median, solid line; 25th and 75th percentiles [p25 and p75], dashed lines) and Panhard et al. (median, dotted and dashed line; p25 and p75, dotted lines) over 60 h (65, 90).


Figure 4
View larger version (26K):
[in this window]
[in a new window]

 
FIG. 4. Predicted levels of ZDV-TP per 106 cells (n = 5,000) for the 200- and 300-mg doses (gray and black, respectively) according to the models of Zhou et al. (median, solid line; 25th and 75th percentiles [p25 and p75], dashed lines) and Panhard et al. (median, dotted and dashed line; p25 and p75, dotted lines) over 60 h (65, 90).


Figure 5
View larger version (17K):
[in this window]
[in a new window]

 
FIG. 5. Box plots (medians and interquartile ranges; whisker plots are p10 to p90 ranges) of peak (checked) and trough (shaded) levels of ZDV-TP in lymphocytes in vivo following ZDV doses of 50, 100, 150, 200, and 300 mg b.i.d. and 600 mg q.d., using the pharmacokinetic model of Zhou et al. (A) and Panhard et al. (B), respectively. Prediction lines were derived by a regression of median values for the 50- and 100-mg b.i.d. doses.


Figure 3
View larger version (34K):
[in this window]
[in a new window]

 
FIG. 3. Simulated cellular levels of ZDV-MP per 106 cells (n = 5,000) for the 200- and 300-mg doses (gray and black, respectively) according to the models of Zhou et al. (median, solid line; 25th and 75th percentiles [p25 and p75], dashed lines) and Panhard et al. (median, dotted and dashed line; p25 and p75, dotted lines) over 60 h (65, 90).

A high degree of overlap of the peak amounts of ZDV-TP was predicted using the two-compartment model of Zhou et al. as well as the one-compartment model of Panhard et al. for ZDV-TP at the 200- and 300-mg b.i.d. doses (98% and 97%, respectively; Fig. 4 and Fig. 5A and B). The respective peak amounts of ZDV-TP at 300 mg b.i.d. were 0.061 pmol/106 cells (interquartile range, 0.031 to 0.12) and 0.062 pmol/106 cells (interquartile range, 0.032 to 0.12) for one- and two-compartment models, respectively. However, the one-compartment model predicted lower trough amounts of ZDV-TP between doses (Fig. 4). Comparison of the median peak amounts of ZDV-TP at steady state in the dose range from 50 to 300 mg b.i.d. with linear predictions using median peak values of the 50- and 100-mg b.i.d. doses suggests some saturation of the phosphorylation to ZDV-TP, demonstrated by the lower-than-expected median peak ZDV-TP amount at 300 mg b.i.d. Furthermore, the simulation based on the two-compartment model suggested higher amounts of cellular ZDV-TP (Fig. 5A) than the one-compartment model (Fig. 5B). A less-than-linear increase in between-dose minimal (trough) amounts of ZDV-TP was suggested by the simulation using the two-compartment model (Fig. 5A) that was not evident for the one-compartment simulation (Fig. 5B). Both simulations predicted reduced peak and trough amounts of cellular ZDV-TP after treatment with 600 mg ZDV q.d. than with 300 or 200 mg ZDV b.i.d.


arrow
DISCUSSION
 
Pharmacokinetic and pharmacodynamic simulations are useful tools for consolidating all available drug information into a usable form and are gaining favor in the pharmaceutical industry to design clinical trials, since they allow detailed analyses of dosage regimens in silico before the actual studies are conducted (13, 21, 56, 57, 63, 64). Rosario et al. recently utilized clinical trial simulations to streamline the phase 2a development of the CCR5 receptor blocking agent maraviroc (75). Furthermore, Hurwitz et al. recently developed a pharmacodynamic model of viral depletion for use with lamivudine (46).

The objective of this study was to superimpose previously reported population pharmacokinetic parameters together with distribution summaries of the cellular enzyme kinetics of ZDV metabolism derived from HIV-1-infected individuals, who were not previously treated with ZDV, using a simulation model. Activated CD4+ lymphocytes are the dominant substrates for HIV-1 infection and could be a significant site for the selection of the K65R mutant virus. The model was then used to predict the accumulation of ZDV nucleotides in activated CD4+ lymphocytes versus dose regimen. Predicted concentrations of ZDV-TP in activated CD4+ PBM cells may be useful for possible incorporation into a virus pharmacokinetic-pharmacodynamic model that relates virus depletion profiles versus time and dose of administration (46, 47).

Simulations were repeated using a one- and a two-compartment population pharmacokinetic model, with different error structures (65, 90). The higher Cmax values and smaller variance after the absorption phase predicted by the two-compartment model of Zhou et al. (90) compared to the one-compartment model may result in part from the different variance structures with a boosted variance during the absorption phase in the two-compartment model. Furthermore, the more gradual terminal slope was expected from the two-compartment model, which predicts a more gradual elimination rate in the latter part of the dosing interval.

The peak cellular levels of ZDV-MP (associated with bone marrow toxicity) predicted by the two-compartment model were higher and the terminal elimination slightly slower than those predicted by the one-compartment model in accordance with the plasma concentration-versus-time profiles of ZDV (Fig. 2). Peak amounts of ZDV-MP were higher in the simulation compared to those reported by Barry et al. in 10 HIV-infected subjects given 300 mg ZDV b.i.d. (2.25 ± 1.5 [mean ± standard deviation] pmol/106 cells; n = 10). Calculation of the ZDV-MP levels assumed a stimulated fraction of PBM cells in humans of 8%. Therefore, ZDV-MP levels in activated cells are expected to be much larger and would exceed the Km value versus TMPK (7.6 µM) for much of the dose interval.

Simulations using the one- and two-compartment pharmacokinetic models predicted a high degree of overlap in cellular ZDV-TP between the 200- and 300-mg b.i.d. doses (97% and 98%, respectively) (Fig. 4). Peak amounts of predicted ZDV-TP at the 300-mg b.i.d. dose were 0.061 pmol/106 cells (interquartile range, 0.031 to 0.12) and 0.062 pmol/106 cells (interquartile range, 0.032 to 0.12), respectively, and were comparable to clinical observations from the study by Barry et al. of 10 HIV-infected individuals given 300 mg ZDV b.i.d. (0.07 ± 0.09 [mean ± standard deviation] pmol/106 cells; n = 10) (7). The corresponding terminal half-lives for ZDV-TP estimated using the median simulated curves were 9.8 and 5.4 h, respectively, compared to 6.5-h half-life estimated using the mean curve for the same dose regimen (7). Interestingly, similar peak amounts of ZDV-TP were reported for the 100- and 300-mg b.i.d. doses in the same study. The more rapid terminal decline in cellular ZDV-TP observed for the lower dose could result from a declining replacement of decaying cellular ZDV-TP from the phosphorylation of ZDV entering cells when concentrations of ZDV are low. Similarly, the more rapid decline in ZDV-TP is expected for the one-compartment model, since it predicts a more rapid terminal elimination of ZDV from the plasma. A decreased replacement in ZDV-TP could also explain the more rapid cellular half-life for ZDV-TP reported in vitro (<2.5 h) (86) than in vivo observations (7 to 10 h) (2, 28, 74, 83), since cellular half-lives are typically measured in vitro after cells are resuspended in medium in the absence of extracellular ZDV.

The higher median ZDV-TP levels predicted by the two-compartment model (Fig. 5A) compared to those of the one-compartment model (Fig. 5B) could be related to the higher predicted plasma Cmax values of ZDV. Comparisons of the median peak amounts of ZDV-TP at steady state, in the dose range of 200 to 300 mg b.i.d., with predictions using a linear regression of median peak values of the 50- and 100-mg b.i.d. doses suggested at least partial saturation of the phosphorylation to ZDV-TP. The less-than-linear increase in the between-dose trough ZDV-TP concentration predicted using the two-compartment model (Fig. 5A) compared with the one-compartment simulation (Fig. 5B) may be explained in part by the more rapid decline in ZDV-TP during the terminal phase of the dosing interval. ZDV-TP acts as a competitive inhibitor and chain terminator during reverse transcription. Therefore, the trough concentrations of ZDV-TP are of critical importance, since CD4+ cells are maximally sensitive to HIV infection when ZDV-TP concentrations are lowest. Both models suggested that similar peak amounts of ZDV-TP follow 600-mg and 300-mg b.i.d. doses, which could result from enzyme saturation and/or a lack of accumulation of ZDV-TP resulting from the prolonged dose interval. The median peak amounts of ZDV-TP predicted for the 600-mg q.d. dose were 0.054 and 0.057 pmol/106 cells, using the two- and one-compartment pharmacokinetic models, respectively, and the corresponding median trough amounts of ZDV-TP were 0.03 and 0.0033 pmol/106 cells, respectively. The trough amounts of ZDV-TP were comparable with those reported by Flynn et al. at the same dose when coadministered with 3TC in 27 HIV-infected adolescents: 0.045 pmol/106 cells (95% confidence interval, 0.035 to 0.065; n = 27) and 0.01 pmol/106 cells (95% confidence interval, 0.01 to 0.02), respectively (28).

Both pharmacokinetic models produced similar overall conclusions for ZDV phosphorylation. However, the population enzymology component of the model was limited by the availability of only one data set each, describing the phosphorylation potential (Vmax) for TK1 and TMPK, respectively, in the PBM cells of infected subjects following ex vivo stimulation (48, 49). These studies indicated substantial variances in the respective Vmax values between HIV-infected individuals. A decrease in TK1 and TMPK activities in the PBM cells of infected individuals was reported in individuals treated with ZDV for more than 6 months. However, the significance of this finding is not clear, since statistically different cellular amounts of ZDV-TP were not demonstrated in a study of subjects exposed to ZDV over a similar time interval (6, 44). ZDV-MP levels tended to be high, while amounts of ZDV-TP agreed closely with in vivo amounts. However, in silico predictions depend on the model structure together with parameter estimates and their associated variance structures. An average estimate for the aqueous cell volume in CD4+ cells was used in these simulations. In reality, cell volumes in an activated lymphocyte population vary according to stage in the cell cycle. The primary substrates for HIV infection are activated T-CD4+ lymphocytes. In most studies, the ZDV-TP content is measured in PBM cells obtained by centrifugation using a Ficoll gradient without further enriching for activated CD4+ cells (7, 28, 83). Since TK1 and TMPK expressions are cell cycle dependent, consideration of data from these measurements underestimates the ZDV-TP content in activated lymphocytes. In this study, the ratio of activated PBM cells was assumed to be constant. However, the activated lymphocyte fraction tends to be elevated in untreated HIV-infected individuals and may decrease toward more normal levels once the infection is stabilized using HAART (59). Therefore, it may be useful for future studies to measure the amounts in activated cells if feasible or to provide a cell cycle distribution together with the ZDV-TP content. The parameter distributions were from different studies, so statistical correlations between parameter covariates could not be assessed. The large variances in enzyme levels between infected individuals make it difficult to estimate the ratios of activities between TK1 (producing ZDV-MP) and TMPK (producing ZDV-DP) (48, 49). This simulation relied on the interplay between the various parameters, and it is possible that there could be a compounding or compensation in the overall error associated with the variance parameters.

The Thai national guidelines for the management of HIV recommend that ZDV doses be reduced from 300 to 200 mg b.i.d. in patients weighing less than 60 kg, which has resulted in fewer side effects and improved long-term tolerability without evidence of reduced efficacy (58). A pharmacokinetic study conducted in Thailand demonstrated similar plasma concentrations in subjects weighing less than 60 kg treated at the 200-mg b.i.d. dose, compared with a historical data set from heavier subjects administered 300 mg b.i.d. (16). Although the oral clearance of ZDV may be reduced in individuals with body weights <60 kg (10, 90), other studies have failed to demonstrate a linear relationship between body weight and ZDV clearance (78). Furthermore, the wide interindividual variance in cellular TK1 and TMPK (45, 48) suggests that plasma concentrations of ZDV alone would be insufficient to predict levels of the active ZDV-TP. The in silico findings of this study suggest that the current ZDV dose could also be lowered from 300 to 200 mg b.i.d. in subjects with body weights more typical of Western populations to reduce toxicities and still maintain adequate ZDV-TP concentrations. However, it would not be prudent to decrease the dose to levels much lower than 200 mg b.i.d., as depicted in Fig. 5A and B, since a more linear decrease was predicted for median peak and trough amounts of ZDV-TP between 50 to 200 mg b.i.d. but not for 300 mg b.i.d. However, the simulations assumed ZDV was not coadministered with agents that alter its metabolism. Therefore, ZDV doses may need modification when coadministered with protease inhibitors, including nevirapine which may induce glucuronidation and/or reduce absorption (65, 90), resulting in lower plasma concentrations. The NRTI TDF, 3TC, and DAPD are not known to affect ZDV metabolism (42, 61).

TDF is administered q.d. in accordance with the long cellular half-life of tenofovir-DP (>60 h). Although the addition of ZDV to q.d. TDF should be beneficial, a coformulation with ZDV would result in a less-effective ZDV response. The half-lives of DXG-TP and carbovir-TP (the active metabolite of ABC) are ~16 h and 12 to 24 h, respectively (42, 79). DAPD is administered b.i.d., while ABC demonstrated similar efficacy when administered q.d. or b.i.d. (14). Therefore, ZDV may be a candidate for coformulation at an optimal b.i.d. dose regimen with an NRTI such as DAPD, the prodrug of DXG or ABC.

Based on these in silico results, an intensive pharmacokinetic phase 2 clinical study of 24 HIV-1-infected subjects receiving 200 or 300 mg ZDV b.i.d. in combination with 500 mg DAPD b.i.d. for 10 days has recently been completed and supports our in silico findings that 200 mg ZDV b.i.d. produces a similar antiviral response as the approved ZDV dose of 300 mg b.i.d. without any untoward effects (61). The data generated from this clinical study will be important for positioning DAPD coformulated with 200 mg ZDV b.i.d. in advanced phase 2/3 studies. Furthermore, the utility of a lower effective dose of ZDV without its associated toxicity, particularly bone marrow toxicity, would be beneficial in future HAART combinations when used with drugs that select for the K65R mutation.


arrow
ACKNOWLEDGMENTS
 
This work was supported in part by the Emory Center for AIDS Research (CFAR) NIH grants 5P30-AI-41980 and 5R37-AI-041980 and the Department of Veterans Affairs.

We thank Monica Hurwitz for graphic-design assistance. R.F.S. is the founder, director, and major shareholder of RFS Pharma, LLC, and as the inventor of DAPD, he may receive future royalties from the sale of this drug. His financial interest in RFS Pharma and DAPD has been reviewed and approved by Emory University in compliance with its conflict-of-interest policies.


arrow
FOOTNOTES
 
* Corresponding author. Mailing address: Veterans Affairs Medical Center, Medical Research 151H, 1670 Clairmont Road, Decatur, GA 30033. Phone: (404) 728-7711. Fax: (404) 728-7726. E-mail: rschina{at}emory.edu Back

{triangledown} Published ahead of print on 6 October 2008. Back


arrow
REFERENCES
 
    1
  1. Acosta, E. P., L. M. Page, and C. V. Fletcher. 1996. Clinical pharmacokinetics of zidovudine. An update. Clin. Pharmacokinet. 30:251-262.[Medline]
  2. 2
  3. Anderson, P. L., T. N. Kakuda, S. Kawle, and C. V. Fletcher. 2003. Antiviral dynamics and sex differences of zidovudine and lamivudine triphosphate concentrations in HIV-infected individuals. AIDS 17:2159-2168.[CrossRef][Medline]
  4. 3
  5. Arner, E. S., T. Spasokoukotskaja, and S. Eriksson. 1992. Selective assays for thymidine kinase 1 and 2 and deoxycytidine kinase and their activities in extracts from human cells and tissues. Biochem. Biophys. Res. Commun. 188:712-718.[CrossRef][Medline]
  6. 4
  7. Arner, E. S., A. Valentin, and S. Eriksson. 1992. Thymidine and 3'-azido-3'-deoxythymidine metabolism in human peripheral blood lymphocytes and monocyte-derived macrophages. A study of both anabolic and catabolic pathways. J. Biol. Chem. 267:10968-10975.[Abstract/Free Full Text]
  8. 5
  9. Balimane, P. V., and P. J. Sinko. 1999. Involvement of multiple transporters in the oral absorption of nucleoside analogues. Adv. Drug Deliv. Rev. 39:183-209.[CrossRef][Medline]
  10. 6
  11. Barry, M., M. Wild, G. Veal, D. Back, A. Breckenridge, R. Fox, N. Beeching, F. Nye, P. Carey, and D. Timmins. 1994. Zidovudine phosphorylation in HIV-infected patients and seronegative volunteers. AIDS 8(8):F1-F5.[Medline]
  12. 7
  13. Barry, M. G., S. H. Khoo, G. J. Veal, P. G. Hoggard, S. E. Gibbons, E. G. Wilkins, O. Williams, A. M. Breckenridge, and D. J. Back. 1996. The effect of zidovudine dose on the formation of intracellular phosphorylated metabolites. AIDS 10:1361-1367.[Medline]
  14. 8
  15. Bazmi, H. Z., J. L. Hammond, S. C. H. Cavalcanti, C. K. Chu, R. F. Schinazi, and J. W. Mellors. 2000. In vitro selection of mutations in the human immunodeficiency virus type 1 reverse transcriptase that decrease susceptibility to (–)-β-D-dioxolane-guanosine and suppress resistance to 3'-azido-3'-deoxythymidine. Antimicrob. Agents Chemother. 44:1783-1788.[Abstract/Free Full Text]
  16. 9
  17. Blum, M. R., S. H. Liao, S. S. Good, and P. de Miranda. 1988. Pharmacokinetics and bioavailability of zidovudine in humans. Am. J. Med. 85:189-194.[Medline]
  18. 10
  19. Burger, D. M., P. L. Meenhorst, C. H. H. ten Napel, J. W. Mulder, C. Neef, C. H. Koks, A. Bult, and J. H. Beijnen. 1994. Pharmacokinetic variability of zidovudine in HIV-infected individuals: subgroup analysis and drug interactions. AIDS 8:1683-1689.[Medline]
  20. 11
  21. Carpenter, C. C., D. A. Cooper, M. A. Fischl, J. M. Gatell, B. G. Gazzard, S. M. Hammer, M. S. Hirsch, D. M. Jacobsen, D. A. Katzenstein, J. S. Montaner, D. D. Richman, M. S. Saag, M. Schechter, R. T. Schooley, M. A. Thompson, S. Vella, P. G. Yeni, and P. A. Volberding. 2000. Antiretroviral therapy in adults: updated recommendations of the International AIDS Society-USA panel. JAMA 283:381-390.[Abstract/Free Full Text]
  22. 12
  23. Chen, H., R. F. Schinazi, P. Rajagopalan, Z. Gao, C. K. Chu, H. M. McClure, and F. D. Boudinot. 1999. Pharmacokinetics of (-)-beta-D-dioxolane guanine and prodrug (-)-beta-D-2,6-diaminopurine dioxolane in rats and monkeys. AIDS Res. Hum. Retrovir. 15:1625-1630.[CrossRef][Medline]
  24. 13
  25. Chien, J. Y., S. Friedrich, M. A. Heathman, D. P. de Alwis, and V. Sinha. 2005. Pharmacokinetics/pharmacodynamics and the stages of drug development: role of modeling and simulation. AAPS J. 7:E544-E559.[CrossRef][Medline]
  26. 14
  27. Cohen, C. J., M. Kubota, P. S. Brachman, W. B. Harley, S. Schneider, V. C. Williams, D. H. Sutherland-Phillips, M. L. Lim, R. B. Balu, and M. S. Shaefer. 2008. Short-term safety and tolerability of a once-daily fixed-dose abacavir-lamivudine combination versus twice-daily dosing of abacavir and lamivudine as separate components: findings from the ALOHA study. Pharmacotherapy 28:314-322.[Medline]
  28. 15
  29. Collier, A. C., S. Bozzette, R. W. Coombs, D. M. Causey, D. A. Schoenfeld, S. A. Spector, C. B. Pettinelli, G. Davies, D. D. Richman, J. M. Leedom, et al. 1990. A pilot study of low-dose zidovudine in human immunodeficiency virus infection. N. Engl. J. Med. 323:1015-1021.[Abstract]
  30. 16
  31. Cressey, T. R., P. Leenasirimakul, G. Jourdain, Y. Tawon, P. O. Sukrakanchana, and M. Lallemant. 2006. Intensive pharmacokinetics of zidovudine 200 mg twice daily in HIV-1-infected patients weighing less than 60 kg on highly active antiretroviral therapy. J. Acquir. Immune Defic. Syndr. 42:387-389.[Medline]
  32. 17
  33. Darbyshire, J., M. Foulkes, R. Peto, W. Duncan, A. Babiker, R. Collins, M. Hughes, T. Peto, and A. Walker. 2000. Immediate versus deferred zidovudine (AZT) in asymptomatic or mildly symptomatic HIV infected adults. Cochrane Database Syst. Rev. CD002039. doi:10.1002/14651858.CD002039.
  34. 18
  35. de Clercq, E. 1996. Non-nucleoside reverse transcriptase inhibitors (NNRTIs) for the treatment of human immunodeficiency virus type 1 (HIV-1) infections: strategies to overcome drug resistance development. Med. Res. Rev. 16:125-157.[CrossRef][Medline]
  36. 19
  37. Deeks, S. G., M. Smith, M. Holodniy, and J. O. Kahn. 1997. HIV-1 protease inhibitors. A review for clinicians. JAMA 277:145-153.[Abstract/Free Full Text]
  38. 20
  39. Dudley, M. N. 1995. Clinical pharmacokinetics of nucleoside antiretroviral agents. J. Infect. Dis. 171(Suppl. 2):S99-S112.
  40. 21
  41. el-Tahtawy, A. A., T. N. Tozer, F. Harrison, L. Lesko, and R. Williams. 1998. Evaluation of bioequivalence of highly variable drugs using clinical trial simulations. II: Comparison of single and multiple-dose trials using AUC and Cmax. Pharm. Res. 15:98-104.[CrossRef][Medline]
  42. 22
  43. Eriksson, B. F., C. K. Chu, and R. F. Schinazi. 1989. Phosphorylation of 3'-azido-2',3'-dideoxyuridine and preferential inhibition of human and simian immunodeficiency virus reverse transcriptases by its 5'-triphosphate. Antimicrob. Agents Chemother. 33:1729-1734.[Abstract/Free Full Text]
  44. 23
  45. Eriksson, S., B. Kierdaszuk, B. Munch-Petersen, B. Oberg, and N. G. Johansson. 1991. Comparison of the substrate specificities of human thymidine kinase 1 and 2 and deoxycytidine kinase toward antiviral and cytostatic nucleoside analogs. Biochem. Biophys. Res. Commun. 176:586-592.[CrossRef][Medline]
  46. 24
  47. Faletto, M. B., W. H. Miller, E. P. Garvey, M. H. St. Clair, S. M. Daluge, and S. S. Good. 1997. Unique intracellular activation of the potent anti-human immunodeficiency virus agent 1592U89. Antimicrob. Agents Chemother. 41:1099-1107.[Abstract]
  48. 25
  49. Feng, J. Y., W. B. Parker, M. L. Krajewski, D. Deville-Bonne, M. Veron, P. Krishnan, Y. C. Cheng, and K. Borroto-Esoda. 2004. Anabolism of amdoxovir: phosphorylation of dioxolane guanosine and its 5'-phosphates by mammalian phosphotransferases. Biochem. Pharmacol. 68:1879-1888.[CrossRef][Medline]
  50. 26
  51. Fischl, M. A., D. D. Richman, D. M. Causey, M. H. Grieco, Y. Bryson, D. Mildvan, O. L. Laskin, J. E. Groopman, P. A. Volberding, R. T. Schooley, et al. 1989. Prolonged zidovudine therapy in patients with AIDS and advanced AIDS-related complex. JAMA 262:2405-2410.[Abstract/Free Full Text]
  52. 27
  53. Fletcher, C. V., E. P. Acosta, K. Henry, L. M. Page, C. R. Gross, S. P. Kawle, R. P. Remmel, A. Erice, and H. H. Balfour, Jr. 1998. Concentration-controlled zidovudine therapy. Clin. Pharmacol. Ther. 64:331-338.[CrossRef][Medline]
  54. 28
  55. Flynn, P. M., J. Rodman, J. C. Lindsey, B. Robbins, E. Capparelli, K. M. Knapp, J. F. Rodriguez, J. McNamara, L. Serchuck, B. Heckman, and J. Martinez. 2007. Intracellular pharmacokinetics of once versus twice daily zidovudine and lamivudine in adolescents. Antimicrob. Agents Chemother. 51:3516-3522.[Abstract/Free Full Text]
  56. 29
  57. Foli, A., K. M. Sogocio, B. Anderson, M. Kavlick, M. W. Saville, M. A. Wainberg, Z. Gu, J. M. Cherrington, H. Mitsuya, and R. Yarchoan. 1996. In vitro selection and molecular characterization of human immunodeficiency virus type 1 with reduced sensitivity to 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA). Antivir. Res. 32:91-98.[CrossRef][Medline]
  58. 30
  59. Furman, P. A., J. A. Fyfe, M. H. St. Clair, K. Weinhold, J. L. Rideout, G. A. Freeman, S. N. Lehrman, D. P. Bolognesi, S. Broder, H. Mitsuya, et al. 1986. Phosphorylation of 3'-azido-3'-deoxythymidine and selective interaction of the 5'-triphosphate with human immunodeficiency virus reverse transcriptase. Proc. Natl. Acad. Sci. USA 83:8333-8337.[Abstract/Free Full Text]
  60. 31
  61. Furman, P. A., J. Jeffrey, L. L. Kiefer, J. Y. Feng, K. S. Anderson, K. Borroto-Esoda, E. Hill, W. C. Copeland, C. K. Chu, J.-P. Sommadossi, I. Liberman, R. F. Schinazi, and G. R. Painter. 2001. Mechanism of action of 1-β-D-2,6-diaminopurine dioxolane, a prodrug of the human immunodeficiency virus type 1 inhibitor 1-β-D-dioxolane guanosine. Antimicrob. Agents Chemother. 45:158-165.[Abstract/Free Full Text]
  62. 32
  63. Gao, W. Y., R. Agbaria, J. S. Driscoll, and H. Mitsuya. 1994. Divergent anti-human immunodeficiency virus activity and anabolic phosphorylation of 2',3'-dideoxynucleoside analogs in resting and activated human cells. J. Biol. Chem. 269:12633-12638.[Abstract/Free Full Text]
  64. 33
  65. Gilliam, B. L., M. M. Sajadi, A. Amoroso, C. E. Davis, F. R. Cleghorn, and R. R. Redfield. 2007. Tenofovir and abacavir combination therapy: lessons learned from an urban clinic population. AIDS Patient Care STDS 21:240-246.[Medline]
  66. 34
  67. Gitterman, S. R., G. L. Drusano, M. J. Egorin, H. C. Standiford, et al. 1990. Population pharmacokinetics of zidovudine. Clin. Pharmacol. Ther. 48:161-167.[Medline]
  68. 35
  69. Gripshover, B. M., H. Ribaudo, J. Santana, J. G. Gerber, T. B. Campbell, E. Hogg, B. Jarocki, S. M. Hammer, and D. R. Kuritzkes. 2006. Amdoxovir versus placebo with enfuvirtide plus optimized background therapy for HIV-1-infected subjects failing current therapy (AACTG A5118). Antivir. Ther. 11:619-623.[Medline]
  70. 36
  71. Gu, Z., M. A. Wainberg, N. Nguyen-Ba, L. L'Heureux, J. M. de Muys, T. L. Bowlin, and R. F. Rando. 1999. Mechanism of action and in vitro activity of 1',3'-dioxolanylpurine nucleoside analogues against sensitive and drug-resistant human immunodeficiency virus type 1 variants. Antimicrob. Agents Chemother. 43:2376-2382.[Abstract/Free Full Text]
  72. 37
  73. Gu, Z., M. A. Wainberg, P. Nguyen-Ba, L. L'Heureux, J. M. de Muys, and R. F. Rando. 1999. Anti-HIV-1 activities of 1,3-dioxolane guanine and 2,6-diaminopurine dioxolane. Nucleosides Nucleotides 18:891-892.[Medline]
  74. 38
  75. Hammond, J. L., D. L. Koontz, H. Z. Bazmi, J. R. Beadle, S. E. Hostetler, G. D. Kini, K. A. Aldern, D. D. Richman, K. Y. Hostetler, and J. W. Mellors. 2001. Alkylglycerol prodrugs of phosphonoformate are potent in vitro inhibitors of nucleoside-resistant human immunodeficiency virus type 1 and select for resistance mutations that suppress zidovudine resistance. Antimicrob. Agents Chemother. 45:1621-1628.[Abstract/Free Full Text]
  76. 39
  77. Hammond, J. L., U. M. Parikh, D. L. Koontz, S. Schlueter-Wirtz, C. K. Chu, H. Z. Bazmi, R. F. Schinazi, and J. W. Mellors. 2005. In vitro selection and analysis of human immunodeficiency virus type 1 resistant to derivatives of β-2',3'-didehydro-2',3'-dideoxy-5-fluorocytidine. Antimicrob. Agents Chemother. 49:3930-3932.[Abstract/Free Full Text]
  78. 40
  79. Hellerstein, M., M. B. Hanley, D. Cesar, S. Siler, C. Papageorgopoulos, E. Wieder, D. Schmidt, R. Hoh, R. Neese, D. Macallan, S. Deeks, and J. M. McCune. 1999. Directly measured kinetics of circulating T lymphocytes in normal and HIV-1-infected humans. Nat. Med. 5:83-89.[CrossRef][Medline]
  80. 41
  81. Hernandez-Santiago, B. I., H. Chen, G. Asif, T. Beltran, S. Mao, S. J. Hurwitz, J. Grier, H. M. McClure, C. K. Chu, D. C. Liotta, and R. F. Schinazi. 2005. Pharmacology and pharmacokinetics of the antiviral agent β-D-2',3'-dideoxy-3'-oxa-5-fluorocytidine in cells and rhesus monkeys. Antimicrob. Agents Chemother. 49:2589-2597.[Abstract/Free Full Text]
  82. 42
  83. Hernandez-Santiago, B. I., A. Obikhod, E. Fromentin, S. J. Hurwitz, and R. F. Schinazi. 2007. Cellular pharmacology of 9-(beta-D-1,3-dioxolan-4-yl) guanine and its lack of drug interactions with zidovudine in primary human lymphocytes. Antivir. Chem. Chemother. 18:343-346.[Medline]
  84. 43
  85. Hirschhorn, R., G. Brittinger, K. Hirschhorn, and G. Weissmann. 1968. Studies on lysosomes. XII. Redistribution of acid hydrolases in human lymphocytes stimulated by phytohemagglutinin. J. Cell Biol. 37:412-423.[Abstract/Free Full Text]
  86. 44
  87. Hoggard, P. G., J. Lloyd, S. H. Khoo, M. G. Barry, L. Dann, S. E. Gibbons, E. G. Wilkins, C. Loveday, and D. J. Back. 2001. Zidovudine phosphorylation determined sequentially over 12 months in human immunodeficiency virus-infected patients with or without previous exposure to antiretroviral agents. Antimicrob. Agents Chemother. 45:976-980.[Abstract/Free Full Text]
  88. 45
  89. Hsu, C.-H., R. Hu, G. E. Dutschman, G. Yang, P. Krishnan, H. Tanaka, M. Baba, and Y.-C. Cheng. 2007. Comparison of the phosphorylation of 4'-ethynyl 2',3'-dihydro-3'-deoxythymidine with that of other anti-human immunodeficiency virus thymidine analogs. Antimicrob. Agents Chemother. 51:1687-1693.[Abstract/Free Full Text]
  90. 46
  91. Hurwitz, S. J., G. Asif, and R. F. Schinazi. 2007. Development of a population simulation model for HIV monotherapy virological outcomes using lamivudine. Antivir. Chem. Chemother. 18:329-341.[Medline]
  92. 47
  93. Hurwitz, S. J., and R. F. Schinazi. 2002. Development of a pharmacodynamic model for HIV treatment with nucleoside reverse transcriptase and protease inhibitors. Antivir. Res. 56:115-127.[Medline]
  94. 48
  95. Jacobsson, B., S. Britton, Q. He, A. Karlsson, and S. Eriksson. 1995. Decreased thymidine kinase levels in peripheral blood cells from HIV-seropositive individuals: implications for zidovudine metabolism. AIDS Res. Hum. Retrovir. 11:805-811.[Medline]
  96. 49
  97. Jacobsson, B., S. Britton, Y. Tornevik, and S. Eriksson. 1998. Decrease in thymidylate kinase activity in peripheral blood mononuclear cells from HIV-infected individuals. Biochem. Pharmacol. 56:389-395.[CrossRef][Medline]
  98. 50
  99. Lanier, E. R., N. Givens, C. Stone, P. Griffin, D. Gibb, S. Walker, M. Tisdale, D. Irlbeck, M. Underwood, M. St. Clair, and M. Ait-Khaled. 2004. Effect of concurrent zidovudine use on the resistance pathway selected by abacavir-containing regimens. HIV Med. 5:394-399.[CrossRef][Medline]
  100. 51
  101. Lavie, A., I. Schlichting, I. R. Vetter, M. Konrad, J. Reinstein, and R. S. Goody. 1997. The bottleneck in AZT activation. Nat. Med. 3:922-924.[CrossRef][Medline]
  102. 52
  103. Margolis, D. M., A. L. Mukherjee, C. V. Fletcher, E. Hogg, D. Ogata-Arakaki, T. Petersen, D. Rusin, A. Martinez, and J. W. Mellors. 2007. The use of beta-D-2,6-diaminopurine dioxolane with or without mycophenolate mofetil in drug-resistant HIV infection. AIDS 21:2025-2032.[Medline]
  104. 53
  105. Martin-Carbonero, L., P. Gil, T. Garcia-Benayas, P. Barreiro, F. Blanco, C. de Mendoza, I. Maida, J. Gonzalez-Lahoz, and V. Soriano. 2006. Rate of virologic failure and selection of drug resistance mutations using different triple nucleos(t)ide analogue combinations in HIV-infected patients. AIDS Res. Hum. Retrovir. 22:1231-1235.[Medline]
  106. 54
  107. Merrill, D. P., M. Moonis, T. C. Chou, and M. S. Hirsch. 1996. Lamivudine or stavudine in two- and three-drug combinations against human immunodeficiency virus type 1 replication in vitro. J. Infect. Dis. 173:355-364.[Medline]
  108. 55
  109. Mewshaw, J. P., F. T. Myrick, D. A. Wakefield, B. J. Hooper, J. L. Harris, B. McCreedy, and K. Borroto-Esoda. 2002. Dioxolane guanosine, the active form of the prodrug diaminopurine dioxolane, is a potent inhibitor of drug-resistant HIV-1 isolates from patients for whom standard nucleoside therapy fails. J. Acquir. Immune Defic. Syndr. 29:11-20.[Medline]
  110. 56
  111. Michelson, S., A. Sehgal, and C. Friedrich. 2006. In silico prediction of clinical efficacy. Curr. Opin. Biotechnol. 17:666-670.[Medline]
  112. 57
  113. Miller, R., W. Ewy, B. W. Corrigan, D. Ouellet, D. Hermann, K. G. Kowalski, P. Lockwood, J. R. Koup, S. Donevan, A. El-Kattan, C. S. Li, J. L. Werth, D. E. Feltner, and R. L. Lalonde. 2005. How modeling and simulation have enhanced decision making in new drug development. J. Pharmacokinet. Pharmacodyn. 32:185-197.[CrossRef][Medline]
  114. 58
  115. Ministry of Public Health, Thailand. 2005. National guidelines for the clinical management of HIV infection in children and adults. Ministry of Public Health, AIDS Division, Department of Communicable Disease Control, Nonthaburi, Thailand.
  116. 59
  117. Mohri, H., A. S. Perelson, K. Tung, R. M. Ribeiro, B. Ramratnam, M. Markowitz, R. Kost, A. Hurley, L. Weinberger, D. Cesar, M. K. Hellerstein, and D. D. Ho. 2001. Increased turnover of T lymphocytes in HIV-1 infection and its reduction by antiretroviral therapy. J. Exp. Med. 194:1277-1287.[Abstract/Free Full Text]
  118. 60
  119. Munch-Petersen, B., L. Cloos, G. Tyrsted, and S. Eriksson. 1991. Diverging substrate specificity of pure human thymidine kinases 1 and 2 against antiviral dideoxynucleosides. J. Biol. Chem. 266:9032-9038.[Abstract/Free Full Text]
  120. 61
  121. Murphy, R., C. Zala, C. Ochoa, P. Tharnish, J. Mathew, E. Fromentin, G. Asif, S. J. Hurwitz, N. M. Kivel, and R. F. Schinazi. 2008. Pharmacokinetics and potent anti-HIV-1 activity of amdoxovir plus zidovudine in a randomized double-blind placebo-controlled study, abstr. 794. Abstr. 14th Conf. Retrovir. Opportunistic Infect. Boston, MA.
  122. 62
  123. Nordic Medical Research Council HIV Therapy Group. 1992. Double blind dose-response study of zidovudine in AIDS and advanced HIV infection. BMJ 304:13-17.[Abstract/Free Full Text]
  124. 63
  125. Ogungbenro, K., and L. Aarons. 2007. Design of population pharmacokinetic experiments using prior information. Xenobiotica 37:1311-1330.[Medline]
  126. 64
  127. Ogungbenro, K., I. Gueorguieva, O. Majid, G. Graham, and L. Aarons. 2007. Optimal design for multiresponse pharmacokinetic-pharmacodynamic models—dealing with unbalanced designs. J. Pharmacokinet. Pharmacodyn. 34:313-331.[Medline]
  128. 65
  129. Panhard, X., M. Legrand, A. M. Taburet, B. Diquet, C. Goujard, and F. Mentre. 2007. Population pharmacokinetic analysis of lamivudine, stavudine and zidovudine in controlled HIV-infected patients on HAART. Eur. J. Clin. Pharmacol. 63:1019-1029.[CrossRef][Medline]
  130. 66
  131. Parikh, U. M., L. Bacheler, D. Koontz, and J. W. Mellors. 2006. The K65R mutation in human immunodeficiency virus type 1 reverse transcriptase exhibits bidirectional phenotypic antagonism with thymidine analog mutations. J. Virol. 80:4971-4977.[Abstract/Free Full Text]
  132. 67
  133. Parikh, U. M., D. L. Koontz, C. K. Chu, R. F. Schinazi, and J. W. Mellors. 2005. In vitro activity of structurally diverse nucleoside analogs against human immunodeficiency virus type 1 with the K65R mutation in reverse transcriptase. Antimicrob. Agents Chemother. 49:1139-1144.[Abstract/Free Full Text]
  134. 68
  135. Pereira, C. F., and J. T. Paridaen. 2004. Anti-HIV drug development—an overview. Curr. Pharm. Des. 10:4005-4037.[CrossRef][Medline]
  136. 69
  137. Perez-Elias, M. J., S. Moreno, C. Gutierrez, D. Lopez, V. Abraira, A. Moreno, F. Dronda, J. L. Casado, A. Antela, and M. A. Rodriguez. 2005. High virological failure rate in HIV patients after switching to a regimen with two nucleoside reverse transcriptase inhibitors plus tenofovir. AIDS 19:695-698.[Medline]
  138. 70
  139. Plagemann, P. G., R. M. Wohlhueter, and C. Woffendin. 1988. Nucleoside and nucleobase transport in animal cells. Biochim. Biophys. Acta 947:405-443.[Medline]
  140. 71
  141. Rapp, K. L., M. Ruckstuhl, and R. F. Schinazi. 2007. The combination of zidovudine and amdoxovir prevents the selection of thymidine analogue mutations in primary human lymphocytes. Antivir. Ther. 12:S130.
  142. 72
  143. Reardon, J. E., and W. H. Miller. 1990. Human immunodeficiency virus reverse transcriptase. Substrate and inhibitor kinetics with thymidine 5'-triphosphate and 3'-azido-3'-deoxythymidine 5'-triphosphate. J. Biol. Chem. 265:20302-20307.[Abstract/Free Full Text]
  144. 73
  145. Richman, D. D., M. A. Fischl, M. H. Grieco, M. S. Gottlieb, P. A. Volberding, O. L. Laskin, J. M. Leedom, J. E. Groopman, D. Mildvan, M. S. Hirsch, et al. 1987. The toxicity of azidothymidine (AZT) in the treatment of patients with AIDS and AIDS-related complex. A double-blind, placebo-controlled trial. N. Engl. J. Med. 317:192-197.[Abstract]
  146. 74
  147. Rodriguez, J. F., J. L. Rodriguez, J. Santana, H. Garcia, and O. Rosario. 2000. Simultaneous quantitation of intracellular zidovudine and lamivudine triphosphates in human immunodeficiency virus-infected individuals. Antimicrob. Agents Chemother. 44:3097-3100.[Abstract/Free Full Text]
  148. 75
  149. Rosario, M. C., B. Poland, J. Sullivan, M. Westby, and E. van der Ryst. 2006. A pharmacokinetic-pharmacodynamic model to optimize the phase IIa development program of maraviroc. J. Acquir. Immune Defic. Syndr. 42:183-191.[Medline]
  150. 76
  151. Ruane, P. J., and A. D. Luber. 2004. K65R-associated virologic failure in HIV-infected patients receiving tenofovir-containing triple nucleoside/nucleotide reverse transcriptase inhibitor regimens. MedGenMed 6:31.[Medline]
  152. 77
  153. Ruane, P. J., G. J. Richmond, E. DeJesus, C. E. Hill-Zabala, S. C. Danehower, Q. Liao, J. Johnson, and M. S. Shaefer. 2004. Pharmacodynamic effects of zidovudine 600 mg once/day versus 300 mg twice/day in therapy-naive patients infected with human immunodeficiency virus. Pharmacotherapy 24:307-312.[CrossRef][Medline]
  154. 78
  155. Sahai, J., K. Gallicano, E. Ormsby, G. Garber, and D. W. Cameron. 1994. Relationship between body weight, body surface area and serum zidovudine pharmacokinetic parameters in adult, male HIV-infected patients. AIDS 8:793-796.[Medline]
  156. 79
  157. Schinazi, R. F., B. I. Hernandez-Santiago, and S. J. Hurwitz. 2006. Pharmacology of current and promising nucleosides for the treatment of human immunodeficiency viruses. Antivir. Res. 71:322-334.[CrossRef][Medline]
  158. 80
  159. Siliciano, J. D., and R. F. Siliciano. 2004. A long-term latent reservoir for HIV-1: discovery and clinical implications. J. Antimicrob. Chemother. 54:6-9.[Abstract/Free Full Text]
  160. 81
  161. Siliciano, R. F. 2005. Scientific rationale for antiretroviral therapy in 2005: viral reservoirs and resistance evolution. Top. HIV Med. 13:96-100.[Medline]
  162. 82
  163. Spiga, M. G., D. A. Weidner, C. Trentesaux, R. D. LeBoeuf, and J. P. Sommadossi. 1999. Inhibition of beta-globin gene expression by 3'-azido-3'-deoxythymidine in human erythroid progenitor cells. Antivir. Res. 44:167-177.[Medline]
  164. 83
  165. Stretcher, B. N., A. J. Pesce, P. T. Frame, and D. S. Stein. 1994. Pharmacokinetics of zidovudine phosphorylation in peripheral blood mononuclear cells from patients infected with human immunodeficiency virus. Antimicrob. Agents Chemother. 38:1541-1547.[Abstract/Free Full Text]
  166. 84
  167. Teague, T. K., L. Munn, K. Zygourakis, and B. W. McIntyre. 1993. Analysis of lymphocyte activation and proliferation by video microscopy and digital imaging. Cytometry 14:772-782.[Medline]
  168. 85
  169. Thompson, M. A., H. A. Kessler, J. J. Eron, Jr., J. M. Jacobson, N. Adda, G. Shen, J. Zong, J. Harris, C. Moxham, and F. S. Rousseau. 2005. Short-term safety and pharmacodynamics of amdoxovir in HIV-infected patients. AIDS 19:1607-1615.[Medline]
  170. 86
  171. Tornevik, Y., B. Jacobsson, S. Britton, and S. Eriksson. 1991. Intracellular metabolism of 3'-azidothymidine in isolated human peripheral blood mononuclear cells. AIDS Res. Hum. Retrovir. 7:751-759.[Medline]
  172. 87
  173. Tornevik, Y., B. Ullman, J. Balzarini, B. Wahren, and S. Eriksson. 1995. Cytotoxicity of 3'-azido-3'-deoxythymidine correlates with 3'-azidothymidine-5'-monophosphate (AZTMP) levels, whereas anti-human immunodeficiency virus (HIV) activity correlates with 3'-azidothymidine-5'-triphosphate (AZTTP) levels in cultured CEM T-lymphoblastoid cells. Biochem. Pharmacol. 49:829-837.[CrossRef][Medline]
  174. 88
  175. Wainberg, M. A., M. D. Miller, Y. Quan, H. Salomon, A. S. Mulato, P. D. Lamy, N. A. Margot, K. E. Anton, and J. M. Cherrington. 1999. In vitro selection and characterization of HIV-1 with reduced susceptibility to PMPA. Antivir. Ther. 4:87-94.[Medline]
  176. 89
  177. White, K. L., N. A. Margot, T. Wrin, C. J. Petropoulos, M. D. Miller, and L. K. Naeger. 2002. Molecular mechanisms of resistance to human immunodeficiency virus type 1 with reverse transcriptase mutations K65R and K65R+M184V and their effects on enzyme function and viral replication capacity. Antimicrob. Agents Chemother. 46:3437-3446.[Abstract/Free Full Text]
  178. 90
  179. Zhou, X.-J., L. B. Sheiner, R. T. D'Aquila, M. D. Hughes, M. S. Hirsch, M. A. Fischl, V. A. Johnson, M. Myers, J.-P. Sommadossi, and the National Institute of Allergy and Infectious Diseases AIDS Clinical Trials Group Protocol 241 Investigators. 1999. Population pharmacokinetics of nevirapine, zidovudine, and didanosine in human immunodeficiency virus-infected patients. Antimicrob. Agents Chemother. 43:121-128.[Abstract/Free Full Text]


Antimicrobial Agents and Chemotherapy, December 2008, p. 4241-4250, Vol. 52, No. 12
0066-4804/08/$08.00+0     doi:10.1128/AAC.00054-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.




This article has been cited by other articles:

  • Behrens, G, Gatell, J M, Geretti, A M, Waters, L (2009). Panel discussion. Int J STD AIDS 20: 18-19 [Full Text]  

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Hurwitz, S. J.
Right arrow Articles by Schinazi, R. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Hurwitz, S. J.
Right arrow Articles by Schinazi, R. F.