ABSTRACT
Human immunodeficiency virus type 1 (HIV-1) isolates resistant to (−)-β-d-dioxolane-guanosine (DXG), a potent and selective nucleoside analog HIV-1 reverse transcriptase (RT) inhibitor, were selected by serial passage of HIV-1LAI in increasing drug concentrations (maximum concentration, 30 μM). Two independent selection experiments were performed. Viral isolates for which the DXG median effective concentrations (EC50s) increased 7.3- and 12.2-fold were isolated after 13 and 14 passages, respectively. Cloning and DNA sequencing of the RT region from the first resistant isolate identified a K65R mutation (AAA to AGA) in 10 of 10 clones. The role of this mutation in DXG resistance was confirmed by site-specific mutagenesis of HIV-1LAI. The K65R mutation also conferred greater than threefold cross-resistance to 2′,3′-dideoxycytidine, 2′,3′-dideoxyinosine, 2′,3′-dideoxy-3′-thiacytidine, 9-(2-phosphonylmethoxyethyl)adenine, 2-amino-6-chloropurine dioxolane, dioxolanyl-5-fluorocytosine, and diaminopurine dioxolane but had only marginal effects on 3′-azido-3′-deoxthymidine (AZT) susceptibility. However, when introduced into a genetic background for AZT resistance (D67N, K70R, T215Y, T219Q), the K65R mutation reversed the AZT resistance. DNA sequencing of RT clones derived from the second resistant isolate identified the L74V mutation, previously reported to cause ddI resistance. The L74V mutation also decreased the AZT resistance when the mutation was introduced into a genetic background for AZT resistance (D67N, K70R, T215Y, T219Q) but to a lesser degree than the K65R mutation did. These findings indicate that DXG and certain 2′,3′-dideoxy compounds (e.g., ddI) can select for the same resistance mutations and thus may not be optimal for use in combination. However, the combination of AZT with DXG or its orally bioavailable prodrug (−)-β-d-2,6-diaminopurine-dioxolane should be explored because of the suppressive effects of the K65R and L74V mutations on AZT resistance.
Breakthrough replication of drug-resistant viral variants is a major obstacle in the search for effective treatments for human immunodeficiency virus type 1 (HIV-1) infection. Resistant variants have emerged with each Food and Drug Administration-approved class of antiretroviral agents, including nucleoside analogue reverse transcriptase (RT) inhibitors, nonnucleoside RT inhibitors, and protease inhibitors. Additionally, variants resistant to combinations of antiretroviral agents have been reported (2, 5). Problems with drug toxicities and suboptimal pharmacokinetics have also hindered the development of effective treatments for HIV-1 infection.
In response to these shortcomings, a variety of new antiretroviral compounds have been sought and discovered. Testing of each of these in clinical trials is neither economically feasible nor practical, making in vitro analysis a vital aspect of preclinical evaluation of drug candidates. Many previous studies have demonstrated that in vitro selection and characterization of drug-resistant viral variants is a useful means of identifying the genetic and biochemical mechanisms of resistance to specific compounds (24), as well as determining the time course for the development of resistance (14). This knowledge can help identify the most promising antiretroviral agents for clinical development. Additionally, in vitro studies can identify drug candidates for combination therapy by characterizing cross-resistance mutations as well as suppressor mutations that reverse phenotypic resistance (6, 14, 15).
(−)-β-d-Dioxolane-guanosine (DXG) was recently synthesized and licensed for clinical development (10). DXG is a purine nucleoside analog with the natural dconfiguration and is both a potent and selective inhibitor of HIV-1 and hepatitis B virus in vitro (4, 8, 10, 22, 26). An unique feature of DXG is the substitution of a dioxolane ring for the sugar moiety found in natural nucleosides (Fig.1). The triphosphorylated form of DXG is recognized as a substrate by HIV-1 RT and is subsequently incorporated into DNA, resulting in chain termination because DXG lacks a 3′-OH group necessary for chain extension (Fig. 1) (8). The median effective concentration (EC50) of DXG is ∼27 nM in activated peripheral blood mononuclear cells (PBMCs), making it the most potent purine nucleoside analog against HIV-1 in these cells. In addition, DXG has no apparent cytotoxicity up to 100 μM. One potential limitation of DXG is its limited aqueous solubility (5); however, promising orally bioavailable prodrugs of DXG are available. Specifically, (−)-β-d-2,6-diaminopurine-dioxolane (DAPD) and (−)-β-d-2-amino-6-chloropurine dioxolane (ACPD) have favorable pharmacokinetic properties and undergo rapid biotransformation to DXG in vivo (4, 22). DAPD is in clinical development, and thus, this study was undertaken to determine the potential for HIV-1 to develop resistance to DXG.
Chemical structure of DXG.
MATERIALS AND METHODS
Chemicals.DXG, 2′,3′-dideoxycytidine (ddC), 3′-azido-3′-deoxythymidine (AZT), 3′-deoxy-2′,3′-didehydrothymidine (D4T), 2′,3′-dideoxy-3′-thiacytidine (3TC), ACPD, DAPD (10), 9-(2-phosphonylmethoxyethyl)adenine (PMEA), and β-d-dioxolanyl-5-fluorocytosine (d-FDOC) were synthesized by our group and were provided by R. F. Schinazi. 2′,3′-Dideoxyinosine (ddI) and phosphonoformate (PFA) were purchased from Sigma Chemical Company, St. Louis, Mo. The compounds were prepared as 10 or 30 mM stock solutions in sterile water or dimethyl sulfoxide as appropriate and were stored at −20°C. Stock solutions were diluted to the appropriate concentrations with RPMI 1640 medium (Whittaker MA Bioproducts, Walkersville, Md.) immediately before use.
Cells.MT-2 cells (AIDS Research and Reference Reagent Program, National Institute of Allergy and Infectious Diseases, National Institutes of Health; contributed by D. Richman) were cultured in RPMI 1640 medium supplemented with 10% fetal bovine serum (JRH Biosciences, Lenexa, Kans.), 2 μM l-glutamine, 10 mM HEPES buffer, and antibiotics (50 IU of penicillin per ml, 50 μg of streptomycin per ml). PBMCs were isolated from HIV-1-seronegative donors and were activated with phytohemagglutinin (10 μg/ml; Difco Labs, Detroit, Mich.) for 3 days prior to infection with HIV-1. At the time of infection, PBMCs were transferred to and maintained in RPMI 1640 medium supplemented with 10% fetal bovine serum, 2 μMl-glutamine, 10% interleukin-2 (Cellular Products, Buffalo, N.Y.), and antibiotics.
Viruses.Stock preparations of HIV-1LAI were prepared as described previously (21). Briefly, MT-2 cells were electroporated with 10 μg of plasmid DNA that encodes the HIV-1LAI proviral clone. Culture supernatants were harvested 5 to 7 days after infection at the peak viral cytopathic effect (CPE). Prior to selection for DXG-resistant virus, the plasmid-derived virus was passed 10 times as cell-free virus in MT-2 cells. The infectivity of the virus preparation was determined by threefold endpoint dilution in MT-2 cells (six cultures per dilution), and the 50% tissue culture infective dose was calculated by use of the equation of Reed and Muench (23).
Selection of resistant viruses.Resistant virus was selected in two independent experiments by serial passage of HIV-1LAI in MT-2 cells in increasing concentrations of DXG. MT-2 cells (1.0 × 106) were pretreated with DXG for 2 h prior to inoculation with HIV-1LAI. Virus replication was measured by daily monitoring of viral cytopathic effects (CPEs). At a +2 CPE (more than two syncytia per field at a ×100 magnification), the culture fluids were clarified by centrifugation (200 × g for 10 min) and the supernatant was harvested. A 0.1-ml aliquot of the harvested supernatant was subsequently used to initiate a new passage in fresh MT-2 cells. Virus was passaged two to three times for each concentration, depending on how rapidly the virus grew in the presence of the drug at a particular concentration. The selective pressure was increased from an initial DXG concentration of 2.5 μM to a final concentration of 30 μM after 13 (selection 1) or 14 (selection 2) passages. After every two passages, the log10 reduction in viral infectivity in 30 μM DXG was determined to monitor the passaged virus for a reduction in susceptibility to DXG (19).
Antiviral susceptibility determinations.Virus susceptibility to DXG and other compounds was determined in MT-2 cells as described previously (19). Briefly, MT-2 cells were infected at a multiplicity of 0.01 50% tissue culture infective dose per cell in the presence of serial drug dilutions. Each dilution was tested in triplicate. Culture supernatants were harvested 4 or 7 days after infection and were assayed for p24 antigen production by a commercial assay (DuPont, NEN Products, Wilmington, Del.). Virus susceptibility to DXG was also determined in PBMCs as described previously (25).
Cloning and DNA sequencing of HIV-1 RT.HIV-1 RT was cloned from infected cell lysates as described previously (17, 18). The full-length RT-coding region was amplified by PCR and was ligated into a PCRII TA cloning vector (InVitrogen, San Diego, Calif.).Escherichia coli INVαF′ was subsequently transfected with the vector, and transformants were identified via EcoRI digestion. Plasmid DNA from transformants was purified (Qiagen Inc., Chatsworth, Calif.) and was sequenced as described previously (19, 21).
Production of mutant recombinant HIV-1.Generation of mutant recombinant HIV-1 was performed as described previously (18). Mutant RT was generated via oligonucleotide-directed mutagenesis and was ligated into the pXXHIV-1LAI cloning vector. Cloning was facilitated by the presence of two silent restriction sites at the 5′ and 3′ ends (XmaI andXbaI, respectively) of the pXXHIV-1LAI RT (17, 18). Infectious mutant recombinant HIV-1 was generated by electroporating MT-2 cells with mutated pXXHIV-1LAIclones. The presence of the desired mutations was verified by DNA sequencing.
Statistical analysis.The Wilcoxon rank sum test for independent sample means was used to test the significance of the observed differences in mean EC50s. One-tailed tests were used in comparing the EC50s for wild-type and recombinant viruses as the hypothesis tested was that the EC50s for recombinant viruses were greater than the EC50s for wild-type viruses.
Molecular modeling of DXG.Calculations were performed on a Silicon Graphics Indy workstation with SYBYL, version 6.1a, software (Tripos Associates, St. Louis, Mo.). ddI, DAPD, and DXG were constructed by using the routines available in SYBYL. The geometries of the molecules were optimized, and charges were generated by using PM3 hamiltonian in a MOPAC interface. A grid search was performed on each molecule with 300 increments about the rotatable bonds C-2→N-1→C-1′→O-4′ and C-3′→C-4′→C-5′→O-5′. The global minimum conformation was further relaxed. Energy minimizations were performed by using the Tripos force field with MOPAC charges by the conjugate gradient method. These molecules were superimposed for comparisons. The molecules were fit by use of the “fit atoms” option in SYBYL by using N-1, C-2, N-3, C-4, N-7, C-8, N-9, and C-1′ atoms in such a way that the root mean square deviation is minimized.
RESULTS
DXG-resistant HIV-1LAI isolates were selected in vitro in two separate experiments by serial passage of virus in MT-2 cells in escalating drug concentrations. In the first experiment, the EC50 of DXG increased only about twofold following the first seven passages (Fig. 2). By passage 10, however, a fivefold increase in the EC50 was apparent compared to the EC50 for control HIV-1LAIpassaged in the absence of drug. The level of DXG resistance reached a maximum of 7.3-fold by passage 13 and did not increase further with subsequent passages. In the second selection experiment, DXG-resistant virus emerged at a rate similar to that observed in the first selection, reaching its maximal resistance level of 12.2-fold by passage 14 (data not shown).
In vitro resistance to DXG.
Genotypic analysis.To determine the molecular basis for DXG resistance, the RT region (amino acids 1 to 560) of proviral DNA was amplified by PCR and was cloned from the infected cells of passage 13 (selection 1) and passage 14 (selection 2). As a control, proviral DNA was amplified and cloned from cells infected with virus passaged in the absence of drug. DNA sequencing identified the presence of a lysine (K)-to-arginine (R) mutation at position 65 (AAA to AGA) in 10 of 10 RT clones derived from passage 13 infected cells (selection 1). The K65R mutation was not detected in any of the control clones (0 of 10). Sequence analysis of RT clones derived from the second selection (passage 14 cells) revealed a leucine (L)-to-valine (V) mutation at position 74 (ATT to AGT) in 8 of 10 RT clones but in none of the control clones (0 of 10). No other mutations in the RT region were found in more than one clone from either selection experiment.
Susceptibilities of recombinant viruses.To confirm the suspected roles of the K65R and L74V mutations in DXG resistance, each mutation was introduced individually by site-specific mutagenesis into a wild-type proviral clone (HIV-1LAI). Additionally, each mutation was introduced into a genetic background for AZT resistance (67N, 70R, 215Y, 219Q) to assess the effect of each mutation on AZT susceptibility. Table 1 summarizes the susceptibilities of the K65R mutant recombinant viruses to DXG and AZT assayed in MT-2 cells (inhibition of p24 antigen production). The K65R mutation conferred 8.7-fold DXG resistance, which was similar to the level of DXG resistance observed for virus isolated in the first selection experiment (7.3-fold). A similar level of DXG resistance (8.0-fold) was observed when the K65R mutant was assayed in PBMCs (data not shown). The K65R mutation had a minimal effect on AZT susceptibility in a wild-type background, increasing the AZT EC50 by only 3.3-fold (Table 1). By contrast, the K65R mutation markedly increased susceptibility to AZT when the mutation was introduced into an AZT resistance background, reducing the level of AZT resistance from 47.7- to 1.3-fold (P < 0.05). The AZT resistance mutations had little effect on the level of DXG resistance conferred by the K65R mutation, reducing this from 8.7- to 5.0-fold (P < 0.05).
DXG-AZT mutational interactions
The L74V mutation conferred ∼4.2-fold DXG resistance when introduced into the HIV-1LAI clone. This was ∼3-fold lower DXG resistance than that observed for virus isolated in the second selection experiment (12.2-fold increase in the EC50). Tenfold resistance was observed when the L74V mutant virus was assayed in PBMCs (data not shown). Similar to K65R, the L74V mutation did not reduce susceptibility to AZT in a wild-type background (P = 0.19) (Table 1). The L74V mutation did, however, reduce the level of AZT resistance from 47.7- to 6.3-fold when the mutation was introduced into an AZT resistance background (P < 0.05), although this reduction was less than that observed with the K65R mutation (Table 1). The AZT resistance mutations had essentially no effect on the level of DXG resistance conferred by the L74V mutation (4.2- to 3.7-fold).
Cross-resistance of the K65R mutant.The K65R recombinant virus was tested for susceptibility to different antiretroviral agents in MT-2 cells to assess the level of cross-resistance conferred by this mutation. Table 2 shows that the K65R mutation significantly reduced the susceptibility to the other dioxolane derivatives tested (d-ACPD andd-DAPD), as well as to ddI, ddC, PMEA, d-FDOC, PFA, and 3TC. The K65R mutation had minimal effects on susceptibility to AZT (3.3-fold increase) (Table 1).
Cross-resistance of HIV-1LAI encoding K65R
Computer modeling.Figure 3shows the results of structural modeling of DXG, DAPD, and ddI. The pictures of the fit of the atoms show conformational differences among ddI, DAPD, and DXG. DAPD and DXG are closely related when the pseudorotation angle and other torsional angles are compared. However, ddI has a sugar conformation different from those of DAPD and DXG. The longer bond length of C-2′→O-3′ in DAPD and DXG than that of C-2′→C-3′ in ddI and the electronic effects of O-3′ play a key role in providing the unique conformations of DAPD and DXG.
Computer modeling of DXG, DAPD, and ddI.
DISCUSSION
HIV-1 variants resistant to DXG, a nucleoside analog RT inhibitor, were selected twice by serial passage of HIV-1LAI in increasing concentrations of compound. After 13 and 14 passages, variants with 7.3- and 12.2-fold increases in EC50s, respectively, were isolated. DNA sequencing and site-specific mutagenesis demonstrated that the K65R mutation was responsible for the DXG resistance observed in the first selection experiment. The K65R mutation also conferred cross-resistance to ddI, ddC, 3TC, PMEA, PFA,d-ACPD, d-FDOC, and d-DAPD but produced little change in AZT susceptibility. The K65R mutation was found to reverse the AZT resistance when the mutation was introduced into an AZT resistance background by ∼37-fold. These results are similar to those reported by Zhang et al. (30), who found that the K65R mutation resulted in resistance to ddC and cross-resistance to ddI but no change in susceptibility to AZT.
DNA sequence analysis of resistant virus from the second selection identified a leucine-to-valine substitution at codon 74. An L74V recombinant virus exhibited 4.2-fold resistance, which was ∼3-fold less than the level of resistance observed in the selected virus (12.2-fold). This difference was reproducible in three separate experiments and thus was not likely due to assay variation. Because no other mutations were consistently found in RT clones derived from this selection, mutations outside of the gene for RT, such as in the envelope gene or another HIV gene, could have contributed to the observed DXG resistance in the selected virus. The L74V mutation was also found to reduce the level of AZT resistance when the mutation was introduced into an AZT resistance background by ∼8-fold. A study by St. Clair et al. (28) identified the L74V mutation in ddI-resistant virus isolated from patients who had changed from AZT to ddI therapy and determined that L74V conferred cross-resistance to ddC and increased susceptibility to AZT. Our finding that DXG selects for the K65R and L74V mutations in vitro has been confirmed in recent work by Borroto-Esoda et al. [K. Borroto-Esoda, J. Mewshaw, D. Wakefield, B. Hooper, J. Jeffrey, P. Furman, and B. McCreedy, Antivir. Ther.4(Suppl. 1), abstr. 3, 1999] and Gu et al. (8). Borroto-Esoda et al. [Antivir. Ther. 4(Suppl. 1), abstr. 3, 1999] performed similar selection experiments with DAPD, a prodrug of DXG. Their initial selection with DAPD isolated the K65R mutant, whereas a second selection isolated the L74V mutant. Gu et al. (8) have shown that virus with the K65R, L74V, or M184V mutation exhibits two- to fivefold resistance to DXG or DAPD. The M184V mutation was not selected in our study, but low-level DXG cross-resistance from this mutation is not surprising because it is well known that it confers low-level cross-resistance to ddI and ddC. According to the recent “closed” RT structure published by Huang et al. (9), K65R, L74V, and M184I/V are each located in the region of RT that interacts with the incoming nucleotide. Therefore, mutations at each of these positions could potentially change the interaction of RT with DXG triphosphate.
Differences in in vitro selection methods have been reported to yield different resistance mutations (12). In this study, however, the same selection methods yielded two different resistance mutations (K65R and L74V). As mentioned above, HIV-1LAI was passaged 10 times in the absence of DXG prior to the start of the selection procedures, which likely increased the genetic heterogeneity of the starting virus population. Kinjerski and Buckheit (11) have reported that the genetic composition of the starting virus population may result in the selection of different genetic variants in vitro. In our study, the isolation of two distinct resistant mutants may have been a chance event as a result of the presence of different resistant quasispecies in the starting virus population or as a result of evolution under drug selective pressure.
That ddI, ddC, and DXG select similar resistance mutations suggests that the triphosphate forms of these analogs interact similarly with the active site of RT. Although the three compounds differ from each other in their nucleoside base compositions, they are all β-d enantiomers and contain similar although not identical sugar components. ddI and ddC are both 2′,3′-dideoxy nucleosides, whereas DXG is a 1,3-dioxolane nucleoside. Computer modeling suggests that the triphosphate forms of DXG and ddI (ddATP) have distinct structural configurations (Fig. 3), although these differences are apparently not sufficient to alter the resistance mutations that are selected. By contrast, AZT selects different resistance mutations than DXG or ddI, probably as a result of the large 3′ azido group of AZT. These observations indicate that the base component of a nucleoside has little influence on the resistance mutations that are selected, whereas the sugar component can be an important determinant (e.g., ddI and DXG versus AZT). These findings are supported by the “closed” RT structure published by Huang et al. (9), which showed that the large 3′ azido group of AZT is likely responsible for the selection of mutations on the “rear face” of RT. Conversely, the lack of 3′ substituents on compounds such as DXG, ddI, and ddC is likely responsible for the selection of mutations on the “front face” of RT.
While the suppressive effects of the L74V and M184V mutations have been documented, this property of the K65R mutation has not been previously recognized. Mutations selected by nonnucleoside RT inhibitors (L100I [3], Y181C [13, 14]) and the foscarnet resistance mutations (W88G, E89K, L92I, Q161L [29]), have also been shown to suppress AZT resistance when introduced into an AZT resistance background. The K65R and L74V mutations confer resistance to 2′,3′-dideoxy (7, 16, 27, 28, 30) and dioxolane nucleoside RT inhibitors, while they suppress the phenotypic effects of AZT resistance mutations. This is a significant finding given that DXG and AZT target the same site of HIV-1 RT but appear to have antagonistic resistance mutations. As mentioned above, the lack of AZT cross-resistance conferred by K65R and L74V is likely the result of the different interactions of the two compounds with RT. The suppression of AZT resistance by these mutations, however, is more difficult to explain. Gu et al. (7, 9) reported that the K65R mutation confers resistance to ddC and noted that positions 65 to 71 of RT play important roles in RT binding of deoxynucleoside triphosphate substrates. Additionally, Sluis-Cremer et al. [N. Sluis-Cremer, D. Arion, N. Kaushik, H. Lim, and M. A. Parniak, Antivir. Ther. 4(Suppl. 1), abstr. 29, 1999] recently reported that the K65R mutation results in a 10-fold decrease in the binding of ddCTP to HIV-1 RT. Given that K65R is an important site for deoxynucleoside triphosphate binding [Sluis-Cremer et al., Antivir. Ther. 4(Suppl. 1), abstr. 29, 1999], it is possible that mutations at this position in RT interfere with the process of pyrophosphorolysis, a proposed mechanism of AZT resistance (1, 20). There is in fact some evidence to suggest that a K65R mutant is less able than wild-type virus to undergo pyrophosphorolysis [Sluis-Cremer et al., Antivir. Ther.4(Suppl. 1), abstr. 29, 1999], thus potentially explaining how this mutation reverses AZT resistance. Additional studies on the interactions between the AZT resistance mutations and K65R and L74V would provide more insight into this observation. Of note is that the AZT resistance background used in this study did not include M41L. St Clair et al. (28) showed that while L74V could reverse the AZT resistance caused by T215Y, this effect was lost in the presence of the M41L mutation. Because we did not include the M41L mutation in our studies, it is unknown whether M41L would have a similar effect on the reversal of the AZT resistance caused by K65R.
The results of this study yield several important conclusions. First, there are common patterns of resistance to DXG and 2′,3′-dideoxynucleoside analogs, which suggest that some of these compounds should not be used in combination with each other. Second, that DXG, ddI, and ddC select for similar resistance mutations indicates that the nucleoside base component has little effect on the resistance mutations selected for by a given nucleoside analogue RT inhibitor. Third, the lack of cross-resistance to AZT and the reversal of AZT resistance by DXG resistance mutations provide strong rationales for the use of these compounds in combination with one another.
ACKNOWLEDGMENTS
This work was supported by grants from the Medical Research Service of the U.S. Department of Veterans Affairs (to R.F.S. and J.W.M.), the U.S. Department of Defense (to R.F.S. and J.W.M.), and the National Institutes of Health (grant RO1 AI-25899 [to C.K.C. and R.F.S.] and grant AI-41980 [to R.F.S.]).
FOOTNOTES
- Received 28 October 1999.
- Returned for modification 22 February 2000.
- Accepted 20 March 2000.
- Copyright © 2000 American Society for Microbiology
