Different Evolution of Genotypic Resistance Profiles to Emtricitabine Versus Lamivudine in Tenofovir-Containing Regimens - pdf attached
Download the PDF here |
JAIDS Journal of Acquired Immune Deficiency Syndromes:
1 November 2010
Svicher, Valentina PhD*; Alteri, Claudia PhD*; Artese, Anna PhD; Forbici, Federica PhD; Santoro, Maria Mercedes PhD*; Schols, Dominique PhD; Van Laethem, Kristel PhD; Alcaro, Stefano PhD; Costa, Giosue; Tommasi, Chiara MD; Zaccarelli, Mauro MD; Narciso, Pasquale MD; Antinori, Andrea MD; Ceccherini-Silberstein, Francesca PhD*; Balzarini, Jan MD, PhD; Perno, Carlo Federico MD, PhD*
From the *Department of Experimental Medicine, University of Rome "Tor Vergata", Rome, Italy; Department of Pharmacobiological Sciences, University of Catanzaro "Magna Graecia", Roccelletta di Borgia (CZ), Italy; I.N.M.I. "L. Spallanzani" Rome, Italy; and Department of Microbiology and Immunology, Rega Institute for Medical Research, Katholieke Universiteit Leuven, Belgium.
Supported by grants from the Italian National Institute of Health, the Ministry of University and Scientific Research, Current and Finalized Research of the Italian Ministry of Health. It was also supported from the Concerted Actions of the K.U.Leuven (GOA no. 05/19) (to JB & DS) and the "Fonds voor Wetenschappelijk Onderzoek Vlaanderen" (FWO no. G.0611.09), and from CHAIN, Collaborative HIV and Anti-HIV Drug Resistance Network, Integrated Project no.223131, funded by the European Commission Framework 7 Program.
Part of these data have been presented at the 16th Conference on Retroviruses and Opportunistic Infections, February 8-11, 2009, Montreal, Canada, as poster presentation; XVIII International HIV Drug Resistance Workshop: Basic Principles and Clinical Implications, June 9-12, 2009, Fort Myers, FL, as oral presentation.
Background: To investigate genotypic resistance profiles to emtricitabine + tenofovir (FTC + TDF) in-vivo and in-vitro, and compare them with lamivudine + tenofovir (3TC + TDF).
Methods: Three hundred fifty-two HIV-1 B-subtype pol sequences from 42 FTC + TDF-treated patients, 40 3TC + TDF-treated patients, and 270 patients treated with 3TC plus another nucleoside reverse transcriptase inhibitor (but not TDF). All patients never received FTC, 3TC, and TDF in their previous therapeutic regimen. 3TC/FTC ± TDF resistance was investigated using in vitro selection experiments and docking simulations.
Results: The M184V mutation is less prevalent in FTC + TDF-treated patients than in 3TC + TDF-treated, and 3TC-treated/TDF-naive patients (14.3% versus 40.0%, P = 0.01 and 55.6%, P < 0.001). Multivariable analysis shows that factors correlated with a lower probability of M184V emergence at failure were the use of FTC compared with 3TC [odds ratio (OR): 0.32 (95% confidence interval (CI): 0.10 to 0.99), P = 0.04], the use of boosted protease inhibitor, and the use of TDF [OR: 0.20 (95% CI: 0.11 to 0.37), P < 0.001, and OR: 0.47 (95%CI: 0.22 to 1.01), P = 0.05, respectively]. In vitro selection experiments and docking analysis show that other reverse transcriptase (RT) mutations, even localized in RT connection domain, can be selected by 3TC + TDF or FTC + TDF in M184V absence and can affect RT affinity for 3TC/FTC and/or TDF.
Conclusions: Our study shows lower rates of M184V development in FTC + TDF regimens versus 3TC + TDF and suggests a potential role of boosted protease inhibitors and TDF in delaying the M184V emergence. Novel RT mutational patterns, more complex than currently known, can contribute to 3TC, FTC, and TDF resistance.
Emergence of HIV-1 drug resistance is associated with increased mortality or faster progression of the HIV-1 infection.1-4 In this light, defining factors that can minimize the onset of drug resistance is crucial to maximize the opportunity to set up a long-term successful antiretroviral treatment. Recent studies have highlighted the protective role of boosted protease inhibitors (PIs) on the emergence of mutations associated with resistance not only to protease but also to reverse transcriptase (RT) inhibitors.5-7 In addition, another study has shown a lower emergence of the M184V mutations in patients receiving emtricitabine (FTC) + Tenofovir (TDF) + Efavirenz (EFV) than in patients receiving Lamivudine (3TC) + azidothymidine + EFV.8 These studies highlight the importance to identify those drug combinations that can prevent drug resistance emergence.
3TC and emtricitabine (FTC) are 2 structurally related nucleoside reverse transcriptase inhibitors (NRTIs), often proposed in combination with tenofovir (TDF) as NRTI backbone of antiretroviral regimens against HIV replication in both drug-naive and multiexperienced patients.9 Both 3TC and FTC usually select for the methionine-to-valine mutation in codon 184 (M184V) of RT. This mutation arises rapidly in failing therapy with 3TC and FTC and results in high-level (>100-fold) resistance to these drugs.10-12 This is often preceded by emergence of a methionine-to-isoleucine mutation in codon 184 (M184I), which is quickly replaced by M184V.13,14 Codon M184 is located in the highly conserved YMDD motif of RT, which is directly involved in binding the incoming nucleotide during reverse transcription. Although several studies have clarified the mechanism of action of the M184V/I mutations, very few information is so far available regarding (1) the antiviral potential of 3TC or FTC when combined with TDF; (2) the factors modulating the rate of M184I/V emergence in vivo; (3) the role of other mutations (beyond M184I/V) in resistance to FTC or 3TC when used not only as single drugs but also in combination with other antiretrovirals.
In the attempt to fulfill the above-mentioned gaps, the goal of this study is to define, by the analysis of clinical data supported by in-vitro experiments and modelling studies, characteristics and developments of genotypic resistance profiles to FTC + TDF versus 3TC + TDF.
This study shows that the use of FTC significantly correlates with a decreased probability of M184V emergence at HAART failure compared with the use of 3TC. This result, confirmed both in univariable and multivariable analysis, is fully consistent with a previous study showing a lower occurrence of M184I/V in patients who started FTC + TDF + EFV compared with those who received 3TC + ZDV + EFV.8 Our result has been also confirmed in an independent dataset showing that virological rebound in patients treated with 3TC/TDF plus a NNRTI or a PI (with or without ritonavir) is associated with a higher probability of developing 3 NRTI resistance mutations including the M184V and 1 NNRTI-associated mutation when compared with individuals receiving FTC + TDF.28
The lower prevalence of M184V can be explained by the higher potency of FTC than 3TC, as suggested by the present and previous in vitro studies10,29 and/or by the different intracellular pharmacokinetics of FTC triphosphate versus 3TC triphosphate.30 Indeed, FTC triphosphate has been shown to have a >2-fold higher intracellular half-life than 3TC triphosphate.31,32 This is due to the ability of FTC to inhibit the activity of the cellular efflux proteins, such as the multidrug resistance-associated proteins, that extrude the drugs out of the CD4+ cells.33 The higher intracellular half-life of FTC can increase its ability to inhibit the RT enzyme, thus reducing the probability of emergence of the M184V mutation. In addition to this, the lower prevalence of M184V may be also related to the coformulation of FTC + TDF that could favor the adherence of the patients to the therapeutic regimen.
Our multivariable model also shows that the use of boosted PIs leads to lower rates of M184V emergence at the failure of 2'-deoxycytidine analogue-containing HAART regimen. This is consistent with previous studies highlighting a lower emergence of mutations associated with resistance to both protease and RT inhibitors at the failure of a boosted PI-containing HAART.5-7 In particular, the negative correlation between the use of boosted PIs and the emergence of M184V shown in our study is fully consistent with Gupta et al7 showing a lower occurrence of M184V in patients who started boosted PI-based HAART compared with those who received a NNRTI regimen (21.0% versus 35.3%, P < 0.001).
In addition, our study suggests a potential role of tenofovir in delaying the emergence of M184V at 3TC failure compared with the use of thymidine analogues. Although this correlation is less strong than that observed for boosted PIs, and needs further confirmation in larger cohort of patients, it is consistent with previous studies showing the complete deselection of the M184V mutation during the pressure imposed by tenofovir in vitro selection experiments.34 Another study also showed the reversion of the M184V in rhesus macaques within 24 weeks of treatment with TDF, even in the presence of 3TC or FTC.35 The negative correlation between TDF and the M184V detection at 3TC/FTC failure may be explained by the fact that the M184V mutation has been shown to confer hypersusceptibility to TDF.36 Thus, the co-administration of a 2' deoxycytidine analogue with TDF may potentially delay the emergence of viral strains that failed by the M184V mutation. Further studies are necessary to verify this hypothesis.
Thus, overall findings highlight the importance to identify those drug combinations that can minimize drug resistance to HAART. Understanding this point is crucial to maximize the opportunity for successful and subsequent therapies after viral rebound not only in the western countries but also in resource-limited settings.
Furthermore, this is the first in vitro study showing that, although FTC was in vitro only slightly superior to 3TC alone in suppressing a high virus load as a single compound, the combinations of FTC + TDF were clearly superior to the combinations of 3TC + TDF (statistical analysis: P < 0.0005). This is consistent with the synergism of FTC and TDF reported by previous biochemical studies. In particular, it has been demonstrated that the levels of both TDF diphosphate and FTC triphosphate were significantly higher when these 2 drugs are combined compared with the levels seen with the individual drugs.37 Another recent study showed that the presence of FTC triphosphate can enhance the formation of the dead-end complexes by HIV-1 RT and TDF-terminated DNA, thus resulting in a reduced terminal NRTI excision.38
This study also highlights the appearance of novel mutations involved in resistance to the cytidine analogues 3TC and FTC and also to TDF. This is the case for T69I (occurring in 0.5% of drug-experienced and in 0.2% of drug-naive patients from HIV-1 Stanford Database), whose presence correlated with a significant increase in resistance to 3TC and FTC; a remarkable increase in resistance was observed also for TDF, though it remained not significant. In addition, our docking simulations show the ability of this mutation to affect the binding affinity of the above-mentioned drugs to the RT enzyme. This structural result is consistent with a previous study showing that T69I confers a selective advantage during chemotherapy with 3TC, presumably by narrowing the template-primer binding cleft and the 2' deoxynucleoside triphosphate-binding pocket.39
Of note, our study highlights that other mutations, localized in the connection domain of HIV-1 RT, can be selected under 3TC/FTC pressure. This is the case of P345S mutation selected under the pressure of 3TC + TDF. This mutation (occurring in 0.2% of drug-experienced and completely absent in drug-naive patients from HIV-1 Stanford Database) has been shown to strongly hamper the binding affinity of both NRTIs to RT. Thus, although T69I and P345S are rare mutations in vivo, overall findings support their role in resistance to NRTI.
In addition, other mutations, such as R356K and S379G, selected in the connection and RNAse-H domain of HIV-1 RT (in our in vitro selection experiments) are natural polymorphisms in non-B subtypes [present in of 49.4% and 4.9% of HIV-1 non-B subtypes-infected drug-naive patients, and in 68.0% and 6.0% of HIV-1 non-B subtypes-infected drug-experienced patients (data from HIV-1 Stanford Database)], thus suggesting that their presence could potentially modulate the virological response to 3TC + TDF or FTC + TDF, without a necessary M184V development. Further studies are necessary to verify this hypothesis.
In our in vitro selection experiments, we did not observe the selection of the K65R under the pressure imposed by tenofovir. This can be explained by the fact that we have used the HIV-1 B subtype laboratory strain HIV-1IIIB. Previous studies have shown a slower K65R emergence in HIV-1 B subtype due to the requirement of 2 nucleotide passages from the Lysine (K) to Arginine (R).40 It is conceivable that the in vitro K65R selection in HIV-1 B subtype requires higher concentration of tenofovir as observed in the study by Stone et al.34
In conclusion, our study shows lower rates of M184V development in FTC + TDF regimens versus 3TC + TDF, and suggests a potential role of boosted PIs and TDF in delaying the M184V emergence. Novel RT mutational patterns, more complex than currently known, can contribute to resistance to NRTI, including FTC/3TC. These results may have implications in selecting NRTIs both at the time of starting HAART and at therapeutic failure.
Three hundred fifty-two HIV-1 B-subtype-infected patients treated with FTC + TDF (n = 42), 3TC + TDF (n = 40), and 3TC plus other NRTI (but not TDF) (n = 270) met the inclusion criteria for this study (Table 1). All patients have started treatment after 2000 (year of introduction of TDF in clinical practice); none of them was receiving a mono or dual therapy at the time of resistance testing, and none of them received 3TC, FTC, and TDF in their previous therapeutic regimens. These 3 groups of patients are also comparable in terms of the median baseline viremia and CD4 cell count, of the median change of viremia and CD4 cell count from baseline to failure, of the median time spent under virological failure, and of the percentage of patients treated with PIs or NNRTIs (Table 1).
Prevalence of M184V at Virological Failure In Vivo
Genotypic analysis showed a significantly lower M184V prevalence in patients failing FTC + TDF than in patients failing 3TC + TDF, and in patients failing 3TC plus another NRTI (but not TDF) [6 of 42 (14.3%) versus 16 of 40 (40.0%) versus 150 of 270 (55.6%), P = 0.01 and P < 0.001, respectively] (Fig. 1). This result was confirmed when the analysis was restricted to the subset of patients failing first-line HAART [5 of 39 (12.8%) versus 10 of 29 (34.5%) versus 104 of 203 (51.2%), P = 0.04 and P < 0.001, respectively] (data not shown).
Multivariable analysis confirms that the use of FTC is an independent factor significantly correlated with a lower probability of M184V emergence at failure than 3TC [odds ratio (OR): 0.32 (95% confidence interval (CI): 0.10 to 0.99), P = 0.04] (Table 2). Other factors independently correlated with a lower probability of M184V emergence at 3TC/FTC failure were the use of boosted PIs [OR: 0.19 (95% CI: 0.10 to 0.37), P < 0.001, compared with NNRTI], and the use of TDF [OR: 0.47 (95% CI: 0.21 to 1.01), P = 0.05, compared with thymidine analogues] (Table 2). The use of TDF significantly correlated with a lower M184V emergence compared with abacavir only at univariable analysis [OR: 0.29 (95% CI: 0.12 to 0.69), P = 0.004].
The negative correlation of boosted PIs and TDF with M184V emergence is confirmed when the analysis was restricted only to patients receiving 3TC ± TDF [OR: 0.20 (95% CI: 0.10 to 0.38), P < 0.001, and OR: 0.47 (95% CI: 0.21 to 1.01), P = 0.05] (data not shown).
In Vitro Selection Experiments of Mutations Associated With Resistance to 3TC/FTC and TDF
Activity of FTC, 3TC, and TDF (as Single Compounds or in Combination) in Suppressing HIV-1 Replication
TDF and 3TC or FTC were administered to CEM cell cultures exposed to a high input of HIV-1(IIIB) (190,000 pg of p24). Under these experimental conditions, TDF, FTC, and 3TC inhibited the virus-induced CPE at an EC50 (50% effective concentration) of -10 µM, 0.075 µM, and 0.20 µM, respectively (Supplemental Digital Content 1, http://links.lww.com/QAI/A74). This is in line with previously published data.16,26
Serial subcultivations of all drug-treated virus-infected cell cultures were then performed by keeping the (combination of) drug concentrations unaltered in the new cell cultures. After 10 subcultivations (4-5 days for each subcultivation) each culture was examined for viral cytopathicity and p24 content. Interestingly, only 2 of 16 different 3TC + TDF combinations had resulted in no cytopathicity and no p24 (virus-free cultures) (TDF 10 µM in combination with 3TC 0.2 µM and 0.1 µM). By contrast, 7 of 16 different FTC + TDF combinations resulted virus-free (FTC 0.1 µM in combination with TDF 10 µM, 5 µM, and 2.5 µM, and FTC 0.2 µM in combination with TDF 10 µM, 5 µM, and 2.5 µM, and 0 µM) (Supplemental Digital Content 1, http://links.lww.com/QAI/A74) as confirmed by lack of p24 detection (Supplemental Digital Content 2, http://links.lww.com/QAI/A75). The microscopical reading fully correlated with the p24 measurements. In the p24-negative cultures, no viral DNA could be amplified. In addition, for 5 of 7 p24-negative cultures under treatment with FTC + TDF, the same concentrations of 3TC + TDF correlated with a viral production >450 pg/mL.
Also, it should be noticed that combination of several fixed drug concentrations proved clearly superior in suppressing virus replication when compared with the use of 2-fold higher concentrations of the single drugs. For example, single administration of 0.05 µM FTC and 5 µM TDF did not result in any viral suppression (after the tenth subcultivation) (p24: >1000 pg/mL), whereas the combination of these drug concentrations (0.05 µM FTC + 5 µM TDF) markedly suppressed virus replication (p24: 2.07 pg/mL). Double single-drug concentrations (ie, FTC at 0.1 µM or TDF at 10 µM) allowed virus replication at 818.98 and 540.54 pg/mL, respectively.
Thus, the check-board combinations revealed that although FTC was only slightly superior to 3TC alone in suppressing a high virus load as a single compound, the combinations of FTC + TDF were superior to the combinations of 3TC + TDF in suppressing virus replication and rendering the HIV-exposed cell cultures free of virus (significance calculated by a Student t test, with a 2-tailed distribution: P < 0.0005).
Analysis of the RT gene was performed on all cell cultures showing the presence of HIV replication (detectable CPE and p24). The different combinations of FTC + TDF never selected for the M184I or M184V mutation (Fig. 2). The M184I mutation was observed only under the selective pressure of FTC when used as single compound without TDF (Fig. 2; FTC 0.1 µM). In the case of 3TC + TDF, M184I was selected when 3TC was combined with the lowest concentration of TDF (1.25 µM, about 8-fold lower than the TDF EC50) (Fig. 2). No M184I/V emergence at all was observed by the other combinations of 3TC + TDF and FTC + TDF.
Other RT mutations were selected in-vitro under the selective pressure of 3TC + TDF or FTC + TDF in the absence of M184V/I. In particular, 3TC + TDF selected for the pair of mutations T69I + P345S and also for mutations P176A, H361Y, and S379G, whereas FTC + TDF selected for the pair of mutations P19S + E169K and for T472I and N494S (Fig. 2). Different mutations were selected by TDF when used without 3TC or FTC: S162N, F214L, and K385N (localized in the connection domain of HIV-1 RT) (Fig. 2). No K65R was found in these cell cultures.
We then investigated the role of mutations detected in in vitro selection experiments by analyzing sequences with known NRTI susceptibility (measured by Virco Antivirogram assay) from the Stanford HIV Drug Resistance Database. We found that the presence of T69I (in the absence of M184V and TAMs) significantly correlated with a 63.3-fold (P = 0.03) and 37.8-fold (P = 0.04) increase in 3TC and FTC resistance, respectively (Table 3). This mutation is also associated with a 14.4-fold increase in TDF resistance (though not statistically significant).
Based on the above-mentioned results, we further investigated the role of T69I in resistance to 3TC, FTC, and TDF by molecular modeling. We also used this approach to characterize the role of P345S in resistance to the drugs because this mutation was coselected (in the same viral genome) together with the T69I during the in vitro selection experiments.
After our previously reported approach,27 we applied a distance-filtering criterion to analyze the docking results, assuming a major "reagent-like" control in the NRTI mechanism of action. A distance >6Å was considered to be not compatible with the binding of the drug with RT.
After docking simulations, we found that the distance between 3TC and the oxydrilic group of the nascent DNA chain is increased (>6Å) in the presence of T69I, likely hampering the affinity between the RT and the nucleotide analogue (Fig. 3). Conversely, P345S, alone or in combination with T69I, did not abolish the nucleophilic attack between 3TC and the nascent DNA chain but was associated with an unfavorable binding affinity (ΔG: -2.74 Kcal/mol for wild-type, -1.85 Kcal/mol for P345S, -2.28 Kcal/mol for T69I + P345S).
Similarly, the presence of either T69I or P345S caused a larger distance between the upcoming TDF and the oxydrilic group of the nascent DNA chain (distance >6Å), with a further increase when these 2 mutations were combined (Fig. 4). Such a result was also confirmed by the energetic analysis of the best poses, that, in presence of both substitutions, were always incompatible with the bond formation and correlated with an unproductive energetic profile (ΔG: -0.88 Kcal/mol for wild-type, +0.15 Kcal/mol for T69I, +1.05 Kcal/mol for P345S, -0.25 Kcal/mol for T69I + P345S).
Finally, in the FTC/RT complexes, the copresence of T69I + P345S correlated with an increased distance between FTC and the oxydrilic group of DNA and with a reduced binding affinity with respect to the WT sequence (data not shown). When present as a single mutation, P345S correlated only with an unfavorable energetic profile, whereas T69I did not seem to negatively affect FTC binding (ΔG: -2.74 Kcal/mol for wild-type, -3.24 Kcal/mol for T69I, -2.02 Kcal/mol for P345S, -2.19 Kcal/mol for T69I + P345S).
Overall results support the role of T69I and P345S, alone or in combination, in resistance to the cytidine analogues 3TC and FTC and to TDF.
MATERIALS AND METHODS
Of 1370 HIV-1 B-subtype-infected patients failing a highly active antiretroviral therapy (HAART) regimen containing a 2'-deoxy-cytidine analogues (3TC or FTC), 352 patients treated with FTC + TDF (n = 42), 3TC + TDF (n = 40), and 3TC plus another NRTI (but not TDF) (n = 270) were selected for this study. The following inclusion criteria were used: starting of HAART after 2000; no exposure to 3TC, FTC, and TDF in previous antiretroviral strategies; absence of M184V in a prior genotypic resistance test. Among 352 patients, 271 patients, distributed as follows, failed their first-line HAART regimen: 39 FTC + TDF-treated patients, 29 3TC + TDF-treated patients, and 203 3TC + another NRTI (but not TDF) patients. Virological failure was defined as 2 consecutive determinations of HIV-1 RNA >50 copies per milliliter. All patients enrolled give written informed consent. Data for all patients were stored in a specifically designed anonymous database that included genotypic, demographic, immunologic, virologic, and therapeutic parameters.
HIV genotype analysis was performed on plasma samples by means of a commercially available kit (ViroSeq HIV-1 genotyping system; Abbott Laboratories, Abbott Park, IL). Briefly, RNA was extracted, retrotranscribed by murine leukemia virus RT, and amplified with Amplitaq-Gold polymerase enzyme by using 2 different sequence-specific primers for 40 cycles. Pol-amplified products (containing the entire protease and the first 335 amino acids of the RT open reading frame) were full-length sequenced in sense and antisense orientations by an automated sequencer (ABI 3100) by using 7 different overlapping sequence-specific primers.15 Sequences having a mixture of wild-type and mutant amino acid residues at single positions were considered to have the mutant(s) at that position.
1. Mutation prevalence: The prevalence of mutations associated with 3TC and FTC resistance was calculated in isolates from patients failing their first HAART regimen containing FTC + TDF (n = 42), 3TC + TDF (n = 40), and 3TC + another NRTI but not TDF (n = 270). χ2 tests of independence were used to verify whether the differences in frequency between the different groups of patients were statistically significant.
2. Factor affecting the presence of M184V at failure: A multivariable logistic regression was performed to define factors positively and negatively correlated with the presence of M184V at 3TC/FTC failure. The following variables were considered: subject demographic, year of treatment, use of cytidine analogues, drugs co-administered with 3TC/FTC, viremia and CD4 cell count at baseline, years under treatment, number of previous treatment failures, number of drugs received in the past, previous exposure to nonnucleoside reverse transcriptase inhibitor (NNRTI), previous exposure to NNRTI. The analysis was performed in the whole cohort of 352 HAART-treated patients, and in the subset of 270 patients receiving 3TC ± TDF.
In Vitro Selection Experiments
In a 48-well tray, 105 CEM cells were seeded in 1 ml RPMI-1640 culture medium-containing 10% fetal calf serum and 2 mM L-glutamine. TDF [(R)-PMPA] (at 10, 5, 2.5, 1.25 µM), FTC (at 0.2, 0.1, 0.05 and 0.025 µM) or 3TC (at 0.2, 0.1, 0.05 and 0.025 µM), and combinations of TDF with FTC or with 3TC were added to the wells. Then, 50 µL of HIV-1(IIIB) (190,000 pg p24) was added and the cultures further incubated for 4-5 days. At each subcultivation, 25-50 µL cell suspension was transferred to fresh cell cultures containing the same fixed concentrations of the test compounds. After the 10th subcultivation (2 months), microscopically visible cytopathicity (CPE) and p24 (pg/mL/well) measurements for all subcultures were determined. The cultures were thereafter further subcultivated for 4 additional passages in the absence of the compounds to verify potential virus breakthrough.16 The methodology of RT sequencing is reported in Snoeck et al17.
Association With 3TC, FTC, and TDF Susceptibility
We analyzed genotype-phenotype correlations from the HIV Stanford Drug Resistance Database18 to assess the association of mutations with 3TC, FTC, and TDF susceptibility. In particular, for each drug, we compared the median changes in resistance (n-fold) in relation to sequences with or without specific mutation. The change in resistance (n-fold) was measured by Virco Antivirogram assay. Mann-Whitney test was used to assess statistically significant difference.
Structural and Docking Analysis
The X-ray crystallographic coordinates of RT complexed to DNA deposited in the Protein Data Bank19 with code 1RTD20 were used for the structural analysis. For the 3 inhibitors, 3TC, TDF, and FTC, the heterocyclic base was appropriately modified before docking into the enzyme. Also the heterocyclic moiety of the n + 1th nucleotide in template overhang was modified according to the base complementary with respect to the incoming NRTI. Thus, in 3TC-TP and FTC-TP complexes, the adenine moiety, which was in the original X-ray structure, was changed into guanine, whereas in TDF-TP modeling, it was modified into thymine. Finally, the NRTI triphosphate was manually located in such an orientation paired with its complementary base in the template strand. From this starting model, the analyzed mutants were generated by single-residue replacement in both chains and energy minimized using the united atom AMBER* force field21 and the GB/SA water implicit solvation model22 as implemented in MacroModel ver 7.2 (Schrdinger).23
The fully optimized receptor coordinates were used for the docking simulations, carried out by AutoDock version 4.24 The map box was fixed equal to 70.190Å3, AMBER was used as force field, the receptor was kept rigid, and the drug flexible. This method applied a Lamarckian model of genetics, generating 50 allowed configurations per ligand. We carried out the maps for A, C, HD, N, NA, OA, P, SA, and F probes; the dielectric constant was set equal to 80.
To select the AutoDock generated poses able to receive the nucleophilic attack by the free 3'OH of the DNA template (NRTI-template distance), a distance cut-off equal to 6.0Å was adopted and the best pose, as ranked by the free-energy value, was identified.
All 3dimensional figures were performed using PyMOL graphics and modeling package ver 0.98.25