Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, March 2002, p. 708-715, Vol. 46, No. 3
0066-4804/02/$04.00+0 DOI: 10.1128/AAC.46.3.708-715.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Tsuyoshi Oishi,1 Aiko Okano,1 Teiichirou Shiino,1 Masashi Tatsumi,4 Masakazu Matsuda,1 Hanae Abumi,1 Noboru Takata,5 Satoshi Shirahata,6 Kaneo Yamada,7 Hiroshi Yoshikura,1 and Yoshiyuki Nagai1
AIDS Research Center, National Institute of Infectious Diseases, Musashimurayama, Tokyo 2080011,1 Department of Veterinary Science, National Institute of Infectious Diseases, Shinjyuku, Tokyo 1628640,4 Division of Blood Transfusion Services, Hiroshima University Medical Hospital, Hiroshima 7348551,5 Department of Pediatrics, University of Occupational and Environmental Health, Fukuoka 8128582,6 Institute of Medical Science, St. Marianna University School of Medicine, Kawasaki 2168511, Japan,7 Virco NV, Mechelen, Belgium ,2 Virco Ltd., Cambridge CB4 0GA, United Kingdom3
Received 27 June 2001/ Returned for modification 5 October 2001/ Accepted 3 December 2001
|
|
|---|
|
|
|---|
Today, six protease inhibitors (PIs) are available clinically (10, 16, 18, 25, 28, 30), and all induce drug resistance mutations (2, 17, 23, 24, 27). Interestingly many of these PI resistance mutations are located within the subsites (6, 8), indicating that the acquisition of these mutations might affect protease activity. Indeed, several PI resistance mutations have been reported to demonstrate impaired enzyme activity (14) and reduced viral fitness (15, 19, 34). This reduced activity could be due to some conformational change in the protease, causing a reduced affinity to the natural substrates, p55gag and p160pol precursors (32), or instability of the protease homodimer (33). As flexible as the protease is, accumulation of mutations in the protease continues after the acquisition of so-called primary mutations responsible for drug resistance. The most fit virus will be selected gradually along with the acquisition of additional mutations, which may complement reduced protease activity, and will become the predominant population. A well-known example of such complementary mutations is L63P in protease, which recovers the viral fitness in a background of multiple combinations of other resistance mutations (22). The interactions of accumulated mutations are still not well understood, but there should be more patterns of complementary interaction among the mutations.
It is also plausible that there could be combinations of mutations that are incompatible and might enhance the level of protease activity impairment. However, this type of mutational combination would be difficult to find, as virus with such mutations would be a minor or underrepresented population in vivo. We have been interested in identifying such combinations, as they would provide important information not only on the structure-function relationships of the protease but possibly on the strategic use of drugs. Here we report a specific pair of mutations that greatly impair protease activity, with almost complete loss of viral infectivity and replication capacity. This combination comprises two major drug resistance mutations, a substitution of asparagine (N) for aspartic acid (D) at codon 30 (D30N) and a substitution of methionine (M) for leucine (L) at codon 90 (L90M). D30N is known to be a primary nelfinavir resistance mutation, which appears to be very specific to this inhibitor (24), whereas L90M is a primary mutation responsible for resistance to both nelfinavir and saquinavir (7, 27) and also appears to be associated with resistance to other PIs (31). We demonstrate an extremely low incidence of these two mutations in combination in the clinical setting, suggesting a highly exclusive relationship between the two mutations in the evolutionary process within each patient. Furthermore, we demonstrate their antagonism in vitro, including in protease activity itself as well as viral infectivity and replication capacity.
|
|
|---|
The other study group was derived from patient samples in a large database (n = 16,996). We used the data to define the frequencies of D30N, L90M, N88D, and V82A and their combinations. The samples were obtained from patients who had received antiretroviral therapy. These were submitted in the United States during 1998 and 1999 for routine assessment of drug susceptibility. Due to the nature of sample collection, we were unable to obtain therapeutic and clinical histories of the registered patients. The patients with viral loads of >1,000 HIV-1 RNA copies/ml were selected, and the frequencies of the mutations D30N, L90M, N88D, and V82A as well as their combinations were determined.
Extraction, amplification, and sequencing of viral RNA. Viral RNA was extracted from 200 µl of patient plasma by using a commercially available kit (Roche Diagnostics, Mannheim, Germany). Extracted RNA was reverse transcribed using avian myeloblastosis virus reverse transcriptase (RT) (Takara, Otsu, Japan) and the specific primer DRPRO-2 (5"-ATTTTCAGGCCCATTTTTTGA). Subsequently, the protease region was amplified by nested PCR using Taq (Takara) and Pyrobest (Takara) DNA polymerases, and the specific outer and inner primers were used for outer and inner PCR, respectively. The outer primers used were DRPRO-1 (5"-CCAACAGCCCCACCAGA) and DRPRO-2 in the reverse transcription reaction, and the inner primer pairs used were DRPRO-3 (5"-AGCAGGAGACGATAGACAAGG) and DRPRO-4 (5"-CTGGCTTTAATTTTACTGGTA). Outer and inner PCRs were performed with the following program: a 2-min hold at 95°C and then 25 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 30 s. Nucleotide sequences of amplicons were analyzed by cycle sequencing methods using Big-dye terminator (PE Biosystems, Foster City, Calif.) and an ABI-377 autosequencer (PE Biosystems).
Clonal analysis for acquisition of PI resistance mutations. Clonal analyses of the protease sequences were performed for eight cases, for which plasma samples could be obtained at multiple points, chosen from the PI-experienced cases described above. These patients also demonstrated virological failure with a nelfinavir-containing regimen. The protease genes were cloned, and the nucleotide sequences were determined. For two patients (JPR-1 and JPR-3), the sequence results were aligned by the Clustal-W program of the Phylip package and then analyzed by both the neighbor-joining and the maximum-likelihood methods.
Generation of recombinant virus clones with PI resistance mutations. A recombinant virus was constructed using HXB2 and the NL4-3 clone as a backbone virus DNA. Initially, an NL4-3 DNA ApaI (restriction site at nucleotide 2011 of NL4-3)-KpnI (restriction site at nucleotide 3831) fragment, which contained the complete protease gene and the 5" half of the RT gene, was inserted into pGEM7zf(-) (Promega, Madison, Wis.). The construct was designated pGEM7/4-3ApaI-KpnI. Using pGEM7/4-3ApaI-KpnI as a template, two restriction sites, a NotI site at nucleotide 2275 and a SmaI site at nucleotide 2592, were introduced into the pol gene as synonymous mutations. Mutagenesis reactions were performed by employing the ExSite PCR-based site-directed mutagenesis kit (Stratagene, La Jolla, Calif.). The introduced NotI and SmaI sites corresponded to the 10th codon of protease and the 15th codon of RT, respectively. The following primer sets were used to introduce the restriction sites: NotI, 5"-GATCACTCTTTGGCAGCGGCCGCTGCTCGTCACAATAAAGATA and 5"-TATCTTTATTGTGACGAGCGGCCGCTGCCAAAGAGTGATC; SmaI,5"-AAAATTAAAGCCCGGGATGGATGGCCCAA and 5"-TTGGGCCATCCATCCCGGGCTTTAATTTT. The entire sequence was verified after the mutagenesis reaction, and a clone with introduced NotI and SmaI sites was designated pGEM7/4-3Not-SmaI. An ApaI (restriction site at nucleotide 2011)-to-AgeI (restriction site at nucleotide 3486) fragment, which contained the NotI-to-SmaI fragment, was excised from pGEM7/4-3Not-SmaI and cloned back into the wild-type molecular clone HXB2. This HXB2/NL4-3 recombinant virus clone with NotI and SmaI sites was designated HXB2cv and used as the host vector for the subsequent HIV-1 recombinant viruses with mutant proteases. The six protease mutants D25N, D30N, L90M, N88D, D30N plus L90M (D30N/L90M), and D30N plus N88D plus L90M (D30N/N88D/L90M) were constructed on the basis of the pGEM7/4-3ApaI-KpnI clone. The primer pairs used to construct each protease mutation were as follows: D25N, 5"-TAAAGGAAGCTATATTAAATACAGGAGCAGATG and 5"-CATCTGTCCTGTATTTAATAGAGCTTCCTTTA; D30N, 5"-GGAGCAGATAATACAGTA and 5"-TACTGTATTATCTGCTCC; L90M, 5"-AGAAATCTGATGACTCAG and 5"-CTGAGTCATCAGATTTCT; and N88D, 5"-ATAATTGGAAGAGATCTGTTGACTCAGATT and 5"-AATCTGAGTCAACAGATCTCTTCCAATTAT. After each mutagenesis reaction, the entire sequence was verified and cloned back to HXB2cv as a NotI-SmaI fragment. The full-length virus DNA was transfected into Cos7 cells by using Fugene6 (Roche Molecular Biochemicals), and culture supernatants were harvested 72 h after transfection and used as virus stocks.
Analysis of p55gag and p160pol processing of mutant viruses. The recombinant DNA clones were transfected into 0.5 x 107 Cos7 cells by electroporation using a Gene Pulser apparatus (Bio-Rad, Hercules, Calif.) at 250 V and 250 µF. The transfected cells were transferred to Falcon T-150 (Becton Dickinson, Lincoln Park, N.J.) culture flasks containing prewarmed Dulbecco modified Eagle medium with 10% fetal calf serum, and the whole medium was changed 24 h after transfection. Seventy-two hours after transfection, the culture supernatant was harvested, and virus lysates were prepared according to a previously described method (20) with a minor modification. In brief, cell debris was removed from the collected culture supernatant by centrifugation at 3,000 rpm (Beckman Allegra 6KR centrifuge) for 15 min at 4°C. The cleared medium was filtered through a 0.2-µm-pore-size filter (Millipore, Bedford, Mass.), overlaid onto a 20% sucrose cushion, and centrifuged at 25,000 rpm (Beckman SW28 rotor) for 1.5 h at 4°C by Optima L-90K (Beckman Coulter, Fullerton, Calif.). The virus pellet was dissolved in an appropriate amount of radioimmunoprecipitation assay lysis buffer (0.05 M Tris-Cl [pH 7.2], 0.15 M NaCl, 0.1% sodium dodecyl sulfate, 1.0% TritonX-100, 1.0% sodium deoxycholate) and subjected to sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The proteins in the gels were transferred to Immobilon P membranes (Millipore) by passive transfer. The membranes were incubated with an HIV-1-seropositive serum for 1 h at 37°C, followed by incubation with a biotinylated secondary antibody (Amersham Pharmacia Biotec, Uppsala, Sweden) and then with avidin-horseradish peroxidase (avidin-HRP) (Amersham Pharmacia Biotec). The HIV-1-positive control human serum used in this study was kindly supplied by T. Lee, Harvard School of Public Health. 3,3"-Diaminobenzidine (Sigma, St Louis, Mo.) was used as a substrate for HRP. For gp120 staining, sheep anti-gp120 C-terminal polyclonal antibody (International Enzymes, Inc., Fallbrook, Calif.) was used, followed by incubation with biotinylated anti-sheep immunoglobulin G secondary antibody (Amersham Pharmacia Biotec) and then avidin-HRP. 3,3"-Diaminobenzidine was used as a substrate for HRP.
Infectivity assay and replication kinetics of recombinant viruses with PI resistance mutations. The infectivity of each protease mutant virus was assayed using the reporter cell line MAGIC5.B1-4. Details of this cell line are described elsewhere (21). Briefly, the cell line expresses CD4, CXCR4, and CCR5 on the cell surface and also contains the HIV-1 long terminal repeat-driven ß-galactosidase gene. Thus, HIV-1-infected cells can be visualized as blue cells when stained with X-Gal (5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside). The stocks of PI-resistant viruses were produced from transfected Cos7 cells as described above. The amount of virus in each stock was determined by RT activity assay as described below. MAGIC5.B1-4 cells were plated in six-well culture plates at a density of 105 cells per well the day before infection and inoculated with each virus adjusted to contain same amount of RT activity. Each infection was done in triplicate wells. Forty-eight hours after infection, the cells were fixed with 1% formaldehyde-0.5% glutaraldehyde in phosphate-buffered saline and stained with an X-Gal staining kit (Roche Diagnostics). The number of blue cells in each well was counted under a light microscope, and mean values from three wells were calculated.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Baseline L90M incidence and acquisition of the D30N nelfinavir resistance mutation
|
|
View this table: [in a new window] |
TABLE 2. Clonal analyses of D30N acquisition pattern in eight patients undergoing protease inhibitor treatment
|
![]() View larger version (40K): [in a new window] |
FIG. 1. (a and c) Clinical courses of JPR-1 (a) and JPR-3 (c). Treatment protocols and durations are demonstrated with arrows. AZT, zidovudine; ddI, dideoxyinosine; IDV, idinavir; NFV, nelfinavir; SQV, saquinavir; RTV, ritonavir; ddC, dideoxycytosine. (b and d) Development of clones with various mutations, and their phylogenetic relationships, of JPR-1(b) and JPR-3(d). A to H (JPR-1) and a to f (JPR-3) indicate the sampling time points and are referred to the respective clone names with additional serial numbers. Codon 30 and 90 amino acid patterns of the clones are indicated.
|
It was quite late in the clinical process, after nelfinavir therapy had been switched to saquinavir-ritonavir double-protease therapy, that two D30N/L90M clones (H-4 and H-5) emerged, suggesting a requirement of strong selection pressure for their evolution. Both of the clones appeared to have descended from the D30N lineage. However, their evolutionary process was complex, because they contained several additional substitutions. Interestingly, these substitutions included L10V, K14R, M64I, L63P, and N88D in common, suggesting similar, albeit complex, pathways for the evolution of these two mutants.
A similar phylogenetic analysis of case JPR-3 is shown in Fig. 1d, and the clinical and therapeutic courses are shown in Fig. 1c. The patient is a 24-year old male hemophiliac patient infected in 1980s who also experienced a long treatment history with RT inhibitors and PIs. Again, viral load suppression was not complete during the observation period. It also appeared that the D30N and L90M clones were derived from 30D/90L wild type clones and evolved independently of each other (clusters III and IV). It was also late in nelfinavir therapy (sampling point d and thereafter) that the D30N and L90M double mutants were detected. The evolutionary tree was more complex than that of JPR-1, and hence it was difficult to define the lineage from which the D30N/L90M double mutants emerged. However, at least 8 out of 10 D30N/L90M clones (d-2, d-3, e-1, e-2, e-8, e-9, f-1, and f-2) appeared to be associated with a D30N lineage (Fig. 1d). This is compatible with the view that D30N clones could acquire an additional L90M mutation much more easily than L90M clones could acquire a D30N mutation. Another feature shared with case JPR-1 was that all D30N/L90M mutants contained additional mutations, including N88D and L63P. Overall, the longitudinal clonal analyses of viral populations in the two clinical cases suggested that D30N and L90M mutations tended to evolve independently of each other. The D30N/L90M double mutants could evolve only under strong selection pressure, possibly via acquisition of particular substitutions such as N88D and L63P.
Frequency of PI resistance mutations in a large database. The incidence of D30N, L90M, D30N/L90M, and other PI resistance mutations was estimated using genotype data derived from a large database (n = 16,996). As shown in Table 3, the frequency of the D30N/L90M double mutation was only 0.8%, whereas that of another double mutation, V82A/L90M, was 5.9%. The frequencies of D30N, L90M, N88D, and V82A single point mutations were found to be 6.3, 29.1, 5.2, and 12.3%, respectively. Given that two residues in either pair, D30N and L90M or V82 and L90M, can change independently of each other, one can predict a value of nearly 1.8% (6.3% x 12.3%) for the D30N/L90M double mutation, which is about half of that (3.6%) for V82A/L90M. The actual estimate (0.8%) for D30N/L90M, however, was less than this prediction. Thus, the results derived from the database also suggest mutual incompatibility between the D30N and L90M mutations. In addition, as these data were derived from population-based sequencing, the D30N/L90M double mutants observed may not all have been together on the same viral genome.
|
View this table: [in a new window] |
TABLE 3. Frequencies of PI resistance mutations estimated from a large database (n = 16,996)
|
As shown in Fig. 2, the products, CA (p24) and MA (p17), generated from p55gag precursor were clearly seen for the HXB2cv wild type (lane 2), whereas they were totally absent in nelfinavir-treated virions (lane 3) and in the D25N virions (lane 4). There was residual precursor p55gag in HXB2cv, but its intensity was much higher in nelfinavir-treated virions and D25N virions (lane 3). The mature Pol products, p66 RT and p51 RT, were seen in HXB2cv, whereas they were not detectable in nelfinavir-treated virions and D25N virions. With all three single point mutations, D30N, N88D, and L90M, both p24 and p17 were seen as clearly as they were in HXB2cv (Fig. 2, lanes 5, 6, and 7). At the same time, however, p55gag precursor was found in remarkable amounts in two (D30N and N88D) of the three mutants, whereas it was at a marginal level in the third mutant (L90M), as in HXB2cv. The 39-kDa band could represent an intermediate cleavage product of p55gag and was seen clearly in the first two mutants but only in a trace amount in the third and not at all in HXB2cv. The p160pol products, p66 RT and p51 RT, were detected in these three single point mutants as clearly as in the parental HXB2cv (Fig. 2). These results indicated that the protease could tolerate almost completely the L90M as well as the D30N or N88D mutation.
![]() View larger version (50K): [in a new window] |
FIG. 2. Analyses of Gag-Pol precursor protein processing of the wild type (HXB2cv) and various mutant viruses by Western blotting. Conditions and viruses used are indicated above each lane. The top part of the membrane (a) was probed with anti-gp120 C-terminal polyclonal antibody, and the bottom part (b) was probed with anti-HIV-1-positive human serum. Positions of molecular markers are indicated on the left, and positions of mature proteins, partially processed intermediate products, and unprocessed precursors are indicated on the right. nfv, nelfinavir.
|
Infectivities of recombinant viruses with protease mutations. We next examined the infectivities of protease mutant viruses by using the MAGIC5.B.1-4 reporter cell assay. The virus inocula were adjusted in each case to contain the same level of RT activity. Each infection was performed in triplicate, and mean blue cell numbers were determined. The wild-type and D30N mutant viruses gave 350 ± 28 and 308 ± 15 positive cells, and L90M displayed 463 ± 13 positive cells, the highest among those tested (Table 4). The virus with D30N/L90M mutations gave only 35 ± 9 positive cells, and this was the lowest among those tested. This is consistent with the lowest cleavage capacity of the mutant protease (Fig. 2). The addition of N88D to D30/L90M (Fig. 2) to give the corresponding D30N/L90M/N88D virus resulted in significantly higher infectivity (191 ± 7 positive cells) than for D30N/L90M, which was also consistent with the protease cleavage activities. Thus, the greater the overall protease impairment, the lower the viral infectivity.
|
View this table: [in a new window] |
TABLE 4. Infectivities of protease mutant viruses evaluated with MAGIC5B.1-4 reporter cells
|
|
|
|---|
Combinations of V32I plus I84V/A (3, 5) and G48V plus L90M (9) also have been reported to demonstrate low protease activity and reduced viral fitness; however, the levels of the impairment seem to be less than those we observed with the D30N/L90M pair. Thus, to our knowledge, D30N/L90M is the most incompatible mutation pair described. The double mutation caused almost complete elimination of infectivity and replication capacity and can be regarded as being a virtually mutually exclusive combination. On the other hand, the same series of assays for processing, infectivity, and replication suggested that either D30N or L90M did not compromise protease activity to any significant degree. These results explain why either mutation is selected so frequently but the combination of the two is so infrequent (both in the large database and in the smaller study group). We previously reported that the frequency of D30N acquisition was higher in patients who received nelfinavir as an initial PI than in those who received nelfinavir as an alternative (29). Here we demonstrated that the incidence of D30N acquisition was strikingly low when the L90M mutation already existed as a background mutation.
The longitudinal phylogenetic analysis of virus clones in a nelfinavir treatment failure case, JPR-1, strongly suggested that clones with D30N and those with L90M emerged directly from a 30D/90L wild-type clone and formed two independent clusters during almost the entire observation period. No clear evidence which suggested conversion of L90M to D30N or vice versa was obtained. These data strongly suggested mutual exclusion of D30N and L90M in the process of evolution of PI resistance mutants within a body. However, two clones with D30N/L90M (H-4 and H-5) eventually clearly appeared from the D30N lineage but not from the L90M lineage. Both of the clones possessed additional mutations, many of which were commonly shared. These included L10V, K14R, M64I, L63P, and N88D. Indeed, one such additional change, N88D, was found to restore, at least in part, the impaired processing activity and infectivity. In another nelfinavir treatment failure case, JPR-3, D30N/L90M clones also appeared late in the clinical course. In this case, their ascendants could not be defined well. However, most of them were associated with the D30N lineage. In addition, the L63P and N88D mutations were again commonly found in all D30N/L90M clones from JPR-3. Interestingly, correlation of the N88D mutation with D30N/L90M was also observed in the large database study (data not shown). These additional changes might have compensated sufficiently to allow the enzyme to tolerate the D30N/L90M mutational combination.
In the present study, we did not study drug susceptibility or replication of the viruses under drug selective pressure in vitro. However, in case JPR-1, the predominance of L90M clones (D-3, D-4, E-2 to -7, F-1, F-3, and F-4) during 3TC-d4T-nelfinavir treatment indicates an advantage of L90M for viral replication in vivo. Clones with D30N (G-1 to -6, H-2, and H-3) appeared after the treatment was switched to dideoxyinosine-d4T-saquinavir-ritonavir. It seems that D30N virus has enhanced fitness and higher drug resistance than clones with L90M.
The D30N/L90M protease processed p160pol precursor into p51 and p66 RT proteins, although a little less efficiently than HXB2cv, but it was hardly able to produce the final products of p55gag, MA and CA,. Thus, D30N/L90M protease behaved differently from the completely inactivated D25N protease and the HXB2cv protease whose function was inhibited by 5 µM nelfinavir. These results suggested that a certain level of catalytic activity was retained by D30N/L90M protease, but the substrate-enzyme interaction, particularly interaction with p55gag, might be significantly disturbed or altered.
The exact mechanism of the incompatibility between D30N and L90M is not understood. The three-dimensional structure of the HIV-1 protease dimer shows that the 30D residue is located in the S2 subsite (24), while the 90L residue is rather distal to the catalytic site. Thus, these two loci are not juxtaposed to each other. Therefore, direct steric hindrance or interaction between these residues is unlikely to occur. Interestingly, 88N is located between 30D and 90L, and the N88D mutation was compensatory to some extent. Thus, some surface area including 88N and 90L and the S2 subsite may be involved in or important for anchoring the p55gag substrate on the enzyme. A subtle change of the nonpolar side chain due to L90M and loss of negative charge due to D30N may somehow disturb the enzyme-substrate interaction. Acquisition of a negative charge due to N88D may compensate for the loss of a negative charge at codon 30 due to the D-to-N change.
Here we have identified an antagonistic relationship between the D30N and L90M mutations in the HIV-1 protease in terms of protease function and viral infectivity. We have also suggested that since the protease is highly flexible in structure, it can accommodate additional mutations to overcome incompatible mutations. Theoretically, evolutionary forces drive an impaired protease towards the recovery of its activity (1). Accumulation and selection of such mutations in vivo appeared to have taken a long time period in our clinical cases. The development of mutations that reduce viral fitness might cause a slowing of disease progression and possibly enable the partial restoration of host immune systems. As our analyses were performed mainly by using artificially constructed cloned viruses with limited mutation combinations, further investigations using clinically isolated viruses to identify interactions of drug resistance mutations are required. Accumulation of clinical data may provide guidance for a more rational multiple-drug protocol to help control HIV-1 replication and HIV-1 disease progression.
We thank the clinicians of the Research Committee on Prevention of Developing Illnesses and Therapy for HIV-Infected Patients of Japan for their collaboration and support.
Present address: Visible Genetics, Cambridge CB4 0GA, United Kingdom. ![]()
|
|
|---|
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»