Previous Article | Next Article ![]()
Journal of Virology, July 2006, p. 7020-7027, Vol. 80, No. 14
0022-538X/06/$08.00+0 doi:10.1128/JVI.02747-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Section of Retroviral Therapeutics, Brigham and Women's Hospital, and Division of AIDS, Harvard Medical School, Boston, Massachusetts
Received 31 December 2005/ Accepted 11 May 2006
|
|
|---|
|
|
|---|
Resistance to ZDV results from the sequential accumulation of thymidine analog resistance mutations (TAMs) at reverse transcriptase (RT) codons 41, 67, 70, 210, 215, and 219 (2, 13, 22). The K70R mutation is usually the first mutational change in HIV-1 RT to emerge during ZDV therapy (2). Although this mutation produces only a modest (eightfold) decrease in ZDV susceptibility in molecular clones of HIV-1 (21), variants of HIV-1 carrying the K70R mutation are selected rapidly in vivo after initiation of ZDV monotherapy (7, 20). Subsequently, variants with mutations at codon 215 and 41 emerge and replace the K70R mutants. Two mutations are possible at codon 215 (Y or F), both of which involve double-nucleotide mutations (ACC to TAC [Y] or TTC [F]). The combined presence of M41L and T215Y confers a 60-fold increase in the 50% inhibitory concentration for ZDV and a 1.8-fold increase in the risk of disease progression and/or death (2, 16, 18). Continued evolution leads to the accumulation of mutations at codons 67 and 210 (13). The combined presence of three to six TAMs results in high-level (>500-fold) ZDV resistance and confers cross-resistance to other nucleoside RT inhibitors (35).
Data from several studies suggest that TAMs are found in two distinct clusters. Mutations that occur together with T215Y (including M41L, L210W, and sometimes D67N) constitute the TAM-1 cluster; mutations that occur together with K70R (including D67N, T215F, and K219Q) constitute the TAM-2 cluster (10, 26, 36). The division of TAMs into two distinct clusters has important clinical significance: ZDV-resistant viruses carrying TAM-1 mutations usually are cross-resistant to didanosine and tenofovir, whereas viruses carrying TAM-2 mutations usually remain susceptible to those drugs (27, 30). The T215Y mutation may be found by itself or in combination with M41L and L210W, but T215F rarely occurs alone or with the M41L and L210W mutations. The L210W mutation, which generally occurs in combination with M41L and T215Y, rarely occurs with the T215F or other TAM-2 mutations (37). To explore the virologic basis for this clustering, we compared the relative replicative fitness and infectivity of HIV-1 recombinants carrying various combinations of TAMs in the absence and presence of ZDV.
(These data were presented in part at the following meetings: (i) 11th Conference on Retroviruses and Opportunistic Infections, 10 to 14 February 2004, San Francisco, Calif. [abstr. 638]; and (ii) 13th International HIV-1 Drug Resistance Workshop, 9 to 11 June 2004, Tenerife, Canary Islands, Spain [abstr. 59].)
|
|
|---|
Construction of recombinant marker viruses.
A 3,314-bp fragment encompassing the entire coding sequence of HIV-1 protease and reverse transcriptase (PR-RT) was amplified by PCR from proviral clone pNL4-3 (corresponding to nucleotides 1998 to 5312 of the HIV-1 NL4-3 sequence [http://hiv-web.lanl.gov]) and cloned into vector pGEM-T Easy vector (Promega, Madison, Wisconsin) to generate plasmid pGEM-PR-RT. Subsequently, 1266 nucleotides of the PR-RT-coding sequence (corresponding to nucleotides 2243 to 3509 in the HIV-1 NL4-3 sequence) were deleted, and a unique BstEII restriction enzyme site was introduced at the deletion junction by PCR-based site-directed mutagenesis using the QuikChange site-directed mutagenesis kit (Stratagene, La Jolla, California). The PR-RT gene carrying this deletion was then cloned into pNL4-3 to yield pHIV
PR.RTBstEII. A segment of the Salmonella enterica serovar Typhimurium histidinol dehydrogenase gene (hisD) or the green fluorescent protein gene (GFP) was then introduced into the XhoI site in nef to serve as a sequence tag, yielding plasmid pHIV
PR.RTBstEIInef-hisD or pHIV
PR.RTBstEIInef-GFP, respectively. The hisD gene fragment was amplified with primers 5'-CCGCTCGAGCGATATTCTGGAAAGCAATGCCAG-3' and 5'-CTGGTCTCGAGCAGGTCAGAAAAAATACGTTC-3', and the GFP fragment was amplified with primers 5'-CCGCTCGAGATGGTGAGCAAGGGCGAGGAGC-3' and 5'-CCGCACGTGTCTTCTGCTTGTCGGCCATG-3'.
Site-directed mutagenesis. Mutations at RT codons 41, 67, 70, 210, 215, and 219 were introduced into the pGEM-PR-RT plasmid using the QuikChange site-directed mutagenesis kit (Stratagene). The presence of mutant sequences was confirmed by automated sequencing of the final plasmid clone with an ABI 377 automated sequencer.
Generation of recombinant marker viruses. Infectious recombinant marker viruses carrying the desired mutations in RT were generated as described previously (24). The RT-coding sequence was amplified from proviral DNA of infected cells at the end of virus culture and analyzed by automated DNA sequencing to verify the presence of the correct alleles at codons 41, 67, 70, 210, 215, and 219. Susceptibility of HIV-1 recombinants to ZDV was determined by a standardized assay modified for use with MT-2 cells (15).
Viral replication kinetics. Inoculation of 500 tissue culture infectious doses was used to infect 5 x 105 MT-2 cells (multiplicity of infection [MOI] of 0.001). After incubation for 2 h at 37°C, cells were washed twice with phosphate-buffered saline (PBS) and resuspended in R-10 medium at 7.5 x 104 cells/ml. Triplicate 2-ml cultures were tested in the absence and presence of ZDV (0.20 µM). Viral replication was quantified by measuring HIV-1 p24 antigen production in cell-free culture supernatants. Experiments were performed three or four times, at least once with recombinants carrying the hisD marker and at least once with the GFP marker, and the results were averaged.
Growth competition assays. Recombinant marker viruses of interest carrying the hisD or GFP sequence tag were mixed together at a ratio of 50:50, 80:20, or 20:80 and inoculated onto 1.5 x 106 MT-2 cells suspended in 300 µl of R-10 medium to yield an MOI of 0.001. After incubation at 37°C for 2 h, cells were washed twice with phosphate-buffered saline (PBS), resuspended in 10 ml of R-10 medium at a concentration of 0.15 x 106 cells/ml in 25-cm2 tissue culture flasks, and reincubated (day 0). Cultures were passaged by inoculating 200 µl of supernatant onto 10 x 106 fresh MT-2 cells every 3 or 4 days. The proportions of the two competing variants were estimated by quantifying GFP and hisD sequences present in culture supernatants on days 1, 4, 7, 11, and 14 using real-time reverse transcriptase PCR with an ABI PRISM 7000 sequence detection system (Applied Biosystems, Foster City, California). Primers GFP86F (5'-GACGTAAACGGCCACAAGTTC-3') and GFP180R (5'-TTGCCGGTGGTGCAGAT-3') and the probe 5'-FAM-CAAGCTGACCCTGAAGT-TAMRA-3' (FAM is 6-carboxyfluorescein; TAMRA is tetramethylrhodamine) were used to detect the GFP fragment, and primers HisD13F (5'-TCGAGCGATATTCTGGCAAA-3') and HisD71R (5' CGGAAAGTACGGTCGCTACCT-3') and the probe sequence 5'-VIC-CCTTGCCTGGCATTG-TAMRA-3' were used to detect the hisD fragment. Viral RNA was extracted from culture supernatants using the QIAGEN kit and treated with RNase-free DNase (QIAGEN, Valencia, California). Parameters for the real-time PCR were described previously (24), except that the initial reverse transcription reaction was performed at 50°C for 30 min. Quantitative real-time reverse transcriptase PCR was performed in triplicate for each sample. Control experiments in which two wild-type recombinants tagged with the hisD and GFP markers were competed against each other showed no change in relative proportions of the two viruses over time (data not shown). Similar results were obtained in control experiments using T215Y mutants tagged with these two markers, indicating that the sequence tags had no significant effect on the relative fitness of the recombinants. For each pair of viruses tested, reciprocal growth competition assays were conducted in which the RT genes of interest were linked to hisD and GFP or vice versa; in each case, similar results were obtained. Therefore, data shown represent the means ± standard deviations (SD) of reciprocal experiments.
Estimation of viral fitness. Quantitative estimates of relative viral fitness were calculated as described previously (33). For each pair of recombinant viruses tested, the final ratio of the two viruses was determined by quantitative real-time PCR at day 11 as described above. The relative fitness (w) of each virus was obtained from the average of the results of three independent growth competition assays (inoculated at ratios of 80:20, 50:50, and 20:80) (33). The fitness difference (WD) was estimated by the ratio of the relative fitness values (WD = WM/WL), where WM is the more-fit virus and WL the less-fit virus in the growth competition assay. Growth competition assays that compared two wild-type recombinants carrying the hisD and GFP markers gave a mean fitness difference (±SD) of 1.0 ± 0.02-fold. Therefore, we considered fitness differences greater than 1.1 to be significant.
Infectivity and fitness profile assays. Viral infectivity was determined on MAGI cells (19) (obtained from the NIH AIDS Research and Reference Reagent Program) by a modification of the method of Mammano et al. (25). Cells were plated in DMEM in 24-well plates at 4 x 104 cells/well the day before infection. After the medium was removed, virus stocks at an MOI of 0.05 were adsorbed onto cells in triplicate wells for 2 h. Subsequently, DMEM was added to a final volume of 1 ml/well, and cultures were incubated at 37°C. After 3 days, the medium was removed, cells were washed with PBS and lysed in 500 µl of lysis buffer (5 mM MgCl2, 0.1% NP-40 in PBS), and ß-galactosidase activity was assayed by hydrolysis of chlorophenol red-ß-D-galactopyranoside (CPRG) as described previously (25). Infectivity was determined from the mean CPRG units (optical density at 570 nm) in the triplicate wells and expressed as a percentage of wild-type infectivity. To determine the replicative advantage of mutant viruses as a function of ZDV concentration, we calculated the ratio of mutant to wild-type infectivity (in CPRG units) for each mutant at ZDV concentrations ranging from 0.001 to 10 µM. Three independent experiments were performed for each mutant; the values plotted represent the mean infectivity ratios ± 95% confidence intervals obtained at each drug concentration.
|
|
|---|
![]() View larger version (12K): [in a new window] |
FIG. 1. Replication kinetics of wild-type (WT) or T215Y (215Y) or T215F (215F) recombinant HIV-1 in the absence (A) or presence (B) of 0.20 µM ZDV. The data shown are the means ± SD (error bars) of three or four independent experiments.
|
![]() View larger version (18K): [in a new window] |
FIG. 2. Growth competition assays between recombinant virus carrying T215Y (215Y) or T215F (215F) in the absence (A to C) and presence (D to F) of 2.0 µM ZDV. The T215Y and T215F recombinants were inoculated at a ratio of 80:20 (A and D), 50:50 (B and E), and 20:80 (C and F), respectively. The data shown are the means ± SD (error bars) of two independent experiments.
|
|
View this table: [in a new window] |
TABLE 1. Relative fitness values of HIV-1 recombinants carrying thymidine analog resistance mutations from growth competition assays in the absence or presence of zidovudine
|
![]() View larger version (8K): [in a new window] |
FIG. 3. Fitness profiles of recombinants carrying ZDV resistance mutations at codon 215. The infectivity ratio of the T215Y (A) or T215F (B) mutant to the wild-type virus or the ratio of the T215Y mutant to the T215F mutant (C) (measured as relative ß-galactosidase activity) was determined in the presence of 0.001, 0.01, 0.05, 0.1, 0.5, 1.0, 5.0, and 10.0 µM ZDV. The data shown are the means ± 95% confidence limits (error bars) of three independent experiments.
|
![]() View larger version (11K): [in a new window] |
FIG. 4. Growth competition assays comparing the fitness of recombinants carrying M41L/T215Y (41L/215Y) and M41L/L210W/T215Y (41L/210W/215Y) mutations inoculated at a ratio of 50:50 in the absence (A) or presence (B) of 2 µM ZDV. The data shown are the means ± SD (error bars) of two independent experiments.
|
![]() View larger version (28K): [in a new window] |
FIG. 5. Growth competition assays comparing the fitness of recombinants carrying the L210W or T215F mutation in a TAM-2 background. Panels show the results of competition assays between D67N/K70R/K219Q versus D67N/K70R/T215F/K219Q (A and B), D67N/K70R/K219Q versus D67N/K70R/L210W/K219Q (C and D), and D67N/K70R/T215F/K219Q versus D67N/K70R/L210W/K219Q (E and F) mutants inoculated at a ratio of 50:50 in the absence (A, C, and E) or presence (B, D, and F) of 2 µM ZDV. The data shown are the means ± SD (error bars) of at least two independent experiments.
|
|
|
|---|
Our findings extend results of a previous study in which recombinant viruses carrying ZDV resistance mutations were competed against wild-type virus (12). In that study, mutants carrying the T215Y mutation with or without M41L were less fit than wild-type virus in the absence of drug but had a fitness advantage at ZDV concentrations of 0.01 µM or greater. However, that study did not explore the effect of the T215F mutation on viral fitness. Our results also agree with an earlier study that demonstrated adverse effects of the L210W mutation on viral fitness in the absence of ZDV (13). However, in contrast to that study, we found that the M41L/L210W/T215Y mutant was fitter than a M41L/T215Y mutant in the presence of ZDV. A possible explanation for these discrepant results could be differences in the ZDV concentration testedwhereas the earlier study used a concentration of 5 µM ZDV, we performed the growth competition assays in the presence of 2 µM ZDV. The choice of a lower concentration of ZDV for use in the current study was based on the results of 50% inhibitory concentration determinations (data not shown). Another possible explanation is that we expressed the mutant Hxb2 RTs in NL4-3 recombinants, whereas in the earlier study the RTs were expressed in an Hxb2 background. Given the close sequence homology of Hxb2 and NL4-3 pol genes, it seems unlikely that this factor accounts for the differences in our results.
Initial studies reporting the occurrence of L210W as a ZDV resistance mutation suggested that it was of marginal significance and observed relatively infrequently in clinical samples (13). Subsequent studies have shown that the L210W mutation commonly occurs together with M41L and T215Y and contributes significantly to ZDV resistance (14, 37). Moreover, didanosine and tenofovir are significantly less active against viruses that carry the 41L/210W/215Y combination (27, 30). Likewise, we have found L210W to be present in HIV-1 sequences from approximately 50% of highly treatment-experienced patients tested in our laboratory at the time of screening for entry into salvage therapy trials (M. Marcial and D. R. Kuritzkes, unpublished observations). Results of the current study showing a fitness advantage of the M41L/L210W/T215Y mutant in the presence of ZDV are consistent with the results of clinical studies that demonstrate the importance of the TAM-1 cluster. It would be interesting to explore the fitness of M41L/L210W/T215Y mutants in the presence of other antiretroviral drugs to which this cluster of mutations confers cross-resistance.
Our results provide additional insights that help explain the ordered appearance of ZDV resistance mutations. Previous work showed that the single point mutation at codon 70 (K
R) emerges first but confers relatively low-level ZDV resistance (2). Subsequently, the T215Y mutation emerges in different viral genomes and these viruses overgrow the K70R mutants. Because the T215Y substitution requires a double-nucleotide mutation (ACC to TAC), it would be expected to occur at a much lower frequency than that of the K70R point mutation, thus explaining its later emergence. Continued replication of partially resistant T215Y mutants in the setting of nonsuppressive regimens eventually results in the accumulation of the M41L, L210W, and D67N mutations. Alternatively, the D67N mutation arises in viruses already carrying a K70R mutation, with eventual accumulation of K219Q, K219E, or K219R substitution and the T215F substitution (which, like T215Y, requires a double mutation [ACC to TTC]). The poor replication kinetics and lower fitness of T215F compared to T215Y along with the modest advantage over the wild-type virus in the presence of ZDV most likely explain why this mutation rarely occurs as an initial ZDV resistance mutation or in association with M41L.
Biochemical studies show that the T215Y and T215F mutations confer resistance to ZDV by accelerating the rate at which the terminal azidothymidine monophosphate is removed through phosphorolytic cleavage, thereby allowing DNA polymerization to proceed (1, 28, 29). The aromatic side chains of Y or F at position 215 are proposed on the basis of modeling studies to enhance ATP binding through stabilizing interactions with the adenine ring of ATP, the physiologically relevant pyrophosphate donor for intracellular nucleotide excision (3). Given the structural similarities of the two amino acid side chains, the fitness difference between T215Y and T215F mutants is unexpected. A possible explanation is suggested by the RT crystal structure. Molecular modeling based on published coordinates (4) suggests that the side chain of the amino acid at position 215 is exposed to the aqueous environment at the enzyme surface (Z. Hu and D. R. Kuritzkes, unpublished observations). Whereas the hydroxyl group on a tyrosine side chain at position 215 could still interact favorably with the aqueous solvent, substitution of the hydrophobic aromatic side chain of a phenylalanine at this position potentially could have a destabilizing effect on the enzyme. Additional biochemical and structural studies are needed to explore this hypothesis.
Although our results are consistent with ZDV resistance patterns observed in the clinic, in vitro conditions cannot reproduce the full array of selective pressures to which HIV-1 is subjected in a patient receiving ZDV-containing antiretroviral therapy. For example, the host HLA haplotype can affect the rate and pattern of emergence of drug resistance mutations because of overlap between key drug resistance loci and epitopes recognized by cytotoxic T lymphocytes (11, 17). It is also possible that differences in the intracellular concentrations of nucleoside analog triphosphates and the deoxynucleoside triphosphate pool result in different selective pressures on RT (6). In addition, viral subtype and polymorphisms in RT that are not associated with drug resistance could influence the emergence of specific TAMs through mutational interactions (32).
In conclusion, HIV-1 recombinants carrying the T215Y mutation for thymidine analog resistance have a clear advantage over T215F mutants with respect to relative replicative fitness and infectivity. Differences in the fitness profiles of these two mutants relative to the wild-type virus in the presence of ZDV likely explain the relative frequency with which T215Y and T215F mutants are observed early in the evolution of ZDV-resistant viruses. Likewise, the differential effects of the L210W and T215F mutations on viral fitness in different TAM backbones help explain the uncommon occurrence of the L210W mutation in combination with T215F in clinical isolates. It will be interesting to determine whether similar fitness differences are observed in RT backbones from other (non-B) viral subtypes and to determine the influence of the RT backbone on the pattern of emergence of TAMs.
We thank Stephen Hughes, Paul Boyer, Stefan Sarafianos, Edward Arnold, and Wenhui Li for helpful discussions; Kelly Hartman for expert technical assistance; and Janet Steele and Lindsay Ware for administrative support.
|
|
|---|
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»