Identification of biased amino acid substitution patterns in human immunodeficiency virus type 1 isolates from patients treated with protease inhibitors

Non-active site mutations in the HIV protease: Diminished drug binding affinity is achieved through modulating the hydrophobic sliding mechanism

The global HIV/AIDS epidemic still currently affects approximately 38 million individuals globally. The protease enzyme of the human immunodeficiency virus is a major drug target in antiviral therapy, however, under the influence of reverse transcriptase and in the context of drug pressure, the rapid PR mutation rate contributes significantly to clinical failure. The set of cooperative non-active site mutations, I13V/I62V/V77I, have been associated with reduced inhibitor susceptibility and are the focus of the current study. When compared to the wild-type protease the mutant protease exhibited decreased binding affinities towards ATV and DRV by 64- and 12-fold, respectively, and decreased the overall favourable Gibbs free energy for ATV, DRV, RTV and SQV. Moreover, these mutations decreased the thermal stability of the protease when in complex with ATV and DRV by approximately 6.4 and 4.2 °C, respectively. The crystal structure of the mutant protease revealed that the location of these mutations and their effect on the hydrophobic sliding mechanism may be crucial in their role in resistance.

Cantilever-centric mechanism of cooperative non-active site mutations in HIV protease: Implications for flap dynamics

The HIV-1 protease is an important drug target in antiretroviral therapy due to the crucial role it plays in viral maturation. A greater understanding of the dynamics of the protease as a result of drug-induced mutations has been successfully elucidated using computational models in the past. We performed induced-fit docking studies and molecular dynamics simulations on the wild-type South African HIV-1 subtype C protease and two non-active site mutation-containing protease variants; HP3 PR and HP4 PR. The HP3 PR contained the I13V, I62V, and V77I mutations while HP4 PR contained the same mutations with the addition of the L33F mutation. The simulations were initiated in a cubic cell universe containing explicit solvent, with the protease variants beginning in the fully closed conformation. The trajectory for each simulation totalled 50 ns. The results indicate that the mutations increase the dynamics of the flap, hinge, fulcrum and cantilever regions when compared to the wild-type protease while in complex with protease inhibitors. Specifically, these mutations result in the protease favouring the semi-open conformation when in complex with inhibitors. Moreover, the HP4 PR adopted curled flap tip conformers which coordinated several water molecules into the active site in a manner that may reduce inhibitor binding affinity. The mutations affected the thermodynamic landscape of inhibitor binding as there were fewer observable chemical contacts between the mutated variants and saquinavir, atazanavir and darunavir. These data help to elucidate the biophysical basis for the selection of cooperative non-active site mutations by the HI virus.

Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis

Background

The correlations of genotypic and phenotypic tests with treatment, clinical history and the significance of mutations in viruses of HIV-infected patients are used to establish resistance mutations to protease inhibitors (PIs). Emerging mutations in human immunodeficiency virus type 1 (HIV-1) protease confer resistance to PIs by inducing structural changes at the ligand interaction site. The aim of this study was to establish an in silico structural relationship between natural HIV-1 polymorphisms and unusual HIV-1 mutations that confer resistance to PIs.

Results

Protease sequences isolated from 151 Mexican HIV-1 patients that were naïve to, or subjected to antiretroviral therapy, were examined. We identified 41 unrelated resistance mutations with a prevalence greater than 1%. Among these mutations, nine exhibited positive selection, three were natural polymorphisms (L63S/V/H) in a codon associated with drug resistance, and six were unusual mutations (L5F, D29V, L63R/G, P79L and T91V). The D29V mutation, with a prevalence of 1.32% in the studied population, was only found in patients treated with antiretroviral drugs. Using in silico modelling, we observed that D29V formed unstable protease complexes when were docked with lopinavir, saquinavir, darunavir, tipranavir, indinavir and atazanavir.

Conclusions

The structural correlation of natural polymorphisms and unusual mutations with drug resistance is useful for the identification of HIV-1 variants with potential resistance to PIs. The D29V mutation likely confers a selection advantage in viruses; however, in silico, presence of this mutation results in unstable enzyme/PI complexes, that possibly induce resistance to PIs.

Application of a new informatics tool for contamination screening in the HIV sequencing laboratory

BACKGROUND

Current HIV-1 sequencing-based methods for detecting drug resistance-associated mutations are open and susceptible to contamination. Informatic identification of clinical sequences that are nearly identical to one another may indicate specimen-to-specimen contamination or another laboratory-associated issue.

OBJECTIVES

To design an informatic tool to rapidly identify potential contamination in the clinical laboratory using sequence analysis and to establish reference ranges for sequence variation in the HIV-1 protease and reverse transcriptase regions among a U.S. patient population.

STUDY DESIGN

We developed an open-source tool named HIV Contamination Detection (HIVCD). HIVCD was utilized to make pairwise comparisons of nearly 8000 partial HIV-1 pol gene sequences from patients across the United States and to calculate percent identities (PIDs) for each pair. ROC analysis and standard deviations of PID data were used to determine reference ranges for between-patient and within-patient comparisons and to guide selection of a threshold for identifying abnormally high PID between two unrelated sequences.

RESULTS

The PID reference range for between-patient comparisons ranged from 83.8 to 95.7% while within-patient comparisons ranged from 96 to 100%. Interestingly, 48% of between-patient sequence pairs with a PID>96.5 were geographically related. The selected threshold for abnormally high PIDs was 96 (AUC=0.993, sensitivity=0.980, specificity=0.999). During routine use, HIVCD identified a specimen mix-up and the source of contamination of a negative control.

CONCLUSIONS

In our experience, HIVCD is easily incorporated into laboratory workflow, useful for identifying potential laboratory errors, and contributes to quality testing. This type of analysis should be incorporated into routine laboratory practice.

Cooperative effects of drug-resistance mutations in the flap region of HIV-1 protease

Understanding the interdependence of multiple mutations in conferring drug resistance is crucial to the development of novel and robust inhibitors. As HIV-1 protease continues to adapt and evade inhibitors while still maintaining the ability to specifically recognize and efficiently cleave its substrates, the problem of drug resistance has become more complicated. Under the selective pressure of therapy, correlated mutations accumulate throughout the enzyme to compromise inhibitor binding, but characterizing their energetic interdependency is not straightforward. A particular drug resistant variant (L10I/G48V/I54V/V82A) displays extreme entropy-enthalpy compensation relative to wild-type enzyme but a similar variant (L10I/G48V/I54A/V82A) does not. Individual mutations of sites in the flaps (residues 48 and 54) of the enzyme reveal that the thermodynamic effects are not additive. Rather, the thermodynamic profile of the variants is interdependent on the cooperative effects exerted by a particular combination of mutations simultaneously present.

Multiple routes and milestones in the folding of HIV-1 protease monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease.

Protease polymorphisms in HIV-1 subtype CRF01_AE represent selection by antiretroviral therapy and host immune pressure

BACKGROUND

Most of our knowledge about how antiretrovirals and host immune responses influence the HIV-1 protease gene is derived from studies of subtype B virus. We investigated the effect of protease resistance-associated mutations (PRAMs) and population-based HLA haplotype frequencies on polymorphisms found in CRF01_AE pro.

METHODS

We used all CRF01_AE protease sequences retrieved from the LANL database and obtained regional HLA frequencies from the dbMHC database. Polymorphisms and major PRAMs in the sequences were identified using the Stanford Resistance Database, and we performed phylogenetic and selection analyses using HyPhy. HLA binding affinities were estimated using the Immune Epitope Database and Analysis.

RESULTS

Overall, 99% of CRF01_AE sequences had at least 1 polymorphism and 10% had at least 1 major PRAM. Three polymorphisms (L10 V, K20RMI and I62 V) were associated with the presence of a major PRAM (P < 0.05). Compared to the subtype B consensus, six additional polymorphisms (I13 V, E35D, M36I, R41K, H69K, L89M) were identified in the CRF01_AE consensus; all but L89M were located within epitopes recognized by HLA class I alleles. Of the predominant HLA haplotypes in the Asian regions of CRF01_AE origin, 80% were positively associated with the observed polymorphisms, and estimated HLA binding affinity was estimated to decrease 19-40 fold with the observed polymorphisms at positions 35, 36 and 41.

CONCLUSION

Polymorphisms in CRF01_AE protease gene were common, and polymorphisms at residues 10, 20 and 62 most likely represent selection by use of protease inhibitors, whereas R41K and H69K were more likely attributable to recognition of epitopes by the HLA haplotypes of the host population.

Therapeutic strategies underpinning the development of novel techniques for the treatment of HIV infection

The HIV replication cycle offers multiple targets for chemotherapeutic intervention, including the viral exterior envelope glycoprotein, gp120; viral co-receptors CXCR4 and CCR5; transmembrane glycoprotein, gp41; integrase; reverse transcriptase; protease and so on. Most currently used anti-HIV drugs are reverse transcriptase inhibitors or protease inhibitors. The expanding application of simulation to drug design combined with experimental techniques have developed a large amount of novel inhibitors that interact specifically with targets besides transcriptase and protease. This review presents details of the anti-HIV inhibitors discovered with computer-aided approaches and provides an overview of the recent five-year achievements in the treatment of HIV infection and the application of computational methods to current drug design.

New approaches to HIV protease inhibitor drug design II: testing the substrate envelope hypothesis to avoid drug resistance and discover robust inhibitors

PURPOSE OF REVIEW

Drug resistance results when the balance between the binding of inhibitors and the turnover of substrates is perturbed in favor of the substrates. Resistance is quite widespread to the HIV-1 protease inhibitors permitting the protease to process its 10 different substrates. This processing of the substrates permits the virus HIV-1 to mature and become infectious. The design of HIV-1 protease inhibitors that closely fit within the substrate-binding region is proposed to be a strategy to avoid drug resistance.

RECENT FINDINGS

Cocrystal structures of HIV-1 protease with its substrates define an overlapping substrate-binding region or substrate envelope. Novel HIV-1 protease inhibitors that were designed to fit within this substrate envelope were found to retain high binding affinity and have a flat binding profile against a panel of drug-resistant HIV-1 proteases.

SUMMARY

The avoidance of drug resistance needs to be considered in the initial design of inhibitors to quickly evolving targets such as HIV-1 protease. Using a detailed knowledge of substrate binding appears to be a promising strategy for achieving this goal to obtain robust HIV-1 protease inhibitors.

Resilience to resistance of HIV-1 protease inhibitors: profile of darunavir

The current effectiveness of HAART in the management of HIV infection is compromised by the emergence of extensively cross-resistant strains of HIV-1, requiring a significant need for new therapeutic agents. Due to its crucial role in viral maturation and therefore HIV-1 replication and infectivity, the HIV-1 protease continues to be a major development target for antiretroviral therapy. However, new protease inhibitors must have higher thresholds to the development of resistance and cross-resistance. Research has demonstrated that the binding characteristics between a protease inhibitor and the active site of the HIV-1 protease are key factors in the development of resistance. More specifically, the way in which a protease inhibitor fits within the substrate consensus volume, or "substrate envelope", appears to be critical. The currently available inhibitors are not only smaller than the native substrates, but also have a different shape. This difference in shape underlies observed patterns of resistance because primary drug-resistant mutations often arise at positions in the protease where the inhibitors protrude beyond the substrate envelope but are still in contact with the enzyme. Since all currently available protease inhibitors occupy a similar space (in spite of their structural differences) in the active site of the enzyme, the specific positions where the inhibitors protrude and contact the enzyme correspond to the locations where most mutations occur that give rise to multidrug-resistant HIV-1 strains. Detailed investigation of the structure, thermodynamics, and dynamics of the active site of the protease enzyme is enabling the identification of new protease inhibitors that more closely fit within the substrate envelope and therefore decrease the risk of drug resistance developing. The features of darunavir, the latest FDA-approved protease inhibitor, include its high binding affinity (Kd = 4.5 x 10-12 M) for the protease active site, the presence of hydrogen bonds with the backbone, and its ability to fit closely within the substrate envelope (or consensus volume). Darunavir is potent against both wild-type and protease inhibitor-resistant viruses in vitro, including a broad range of over 4,000 clinical isolates. Additionally, in vitro selection studies with wild-type HIV-1 strains have shown that resistance to darunavir develops much more slowly and is more difficult to generate than for existing protease inhibitors. Clinical studies have shown that darunavir administered with low-dose ritonavir (darunavir/ritonavir) provides highly potent viral suppression (including significant decreases in HIV viral load in patients with documented protease inhibitor resistance) together with favorable tolerability. In conclusion, as a result of its high binding affinity for and overall fit within the active site of HIV-1 protease, darunavir has a higher genetic barrier to the development of resistance and better clinical efficacy against multidrug-resistant HIV relative to current protease inhibitors. The observed efficacy, safety and tolerability of darunavir in highly treatment-experienced patients makes darunavir an important new therapeutic option for HIV-infected patients.

Insight into the folding inhibition of the HIV-1 protease by a small peptide

It has recently been shown that the highly protected segments 24-34 (S2) and 83-93 (S8) of each of the two 99-mers of human immunodeficiency virus type 1 protease play an essential role in the folding of the monomers, giving rise to the so-called (postcritical) folding nucleus ((FN) minimum condensation unit ensuring folding) when they dock. This scenario received further support from model calculations that demonstrated that the peptide p-S8, displaying an amino acid sequence identical to the corresponding (83-93) segment of the monomer, can be used to interfere with the formation of the FN and eventually to inhibit folding by docking the fragment 24-34. Experiments in vitro and in cells infected with ex vivo wild-type and multiresistant HIV isolates confirm that the inhibition power of p-S8 is robust. On the other hand, there is no direct evidence demonstrating the validity of the proposed mechanism of inhibition associated with p-S8. To shed light on this question and to provide the basis for the design of a molecule mimetic to p-S8, to be used as lead of an eventual drug against AIDS, we study, in this paper, with the help of all-atom simulations in explicit solvent and the novel method of metadynamics combined with parallel tempering: a), the free energy and the equilibrium structure of each of the peptides p-S2 and p-S8; b), the details of the docking mechanism of the two peptides and the free energy associated with this process. Whereas p-S8 is found to be well structured, p-S2 is rather flexible, wrapping itself around p-S8 to give rise to the FN, which is stabilized by three particular hydrogen bonds.

Insight into analysis of interactions of saquinavir with HIV-1 protease in comparison between the wild-type and G48V and G48V/L90M mutants based on QM and QM/MM calculations

Saquinavir (SQV) was the first HIV-1 PR inhibitor licensed for clinical use and widely used for acquired immunodeficiency syndrome (AIDS) therapy. Its effectiveness, however, has been hindered by the emergence of resistant mutations. The two most important HIV-1 PR mutants are G48V and G48V/L90M. Inhibition studies of SQV on these mutants demonstrated 13.5- and 419-fold reductions of susceptibility, respectively. In this study, an analysis of energetic binding affinity between saquinavir and the HIV-1 PR wild-type and these two mutants has been performed in detail based on density functional theory and the hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. We have found that the interaction of SQV with flap residue 48 of the mutants is significantly perturbed, as shown by the reduced stability of binding between SQV and residue 48 for the G48V and G48V/L90M mutants over the wild-type. This was associated with conformational changes of the inhibitor and the enzyme, leading to the loss of hydrogen bonding between the binding subsite P2 and the backbone carbonyl of residue 48. Moreover, the G48V/L90M mutations cause the repositioning of the residues close to residues 48 and 90, at important locations as a part of the flap and catalytic regions, respectively. The repositioning of these residues consequently perturbed the binding affinity of SQV in the pocket as indicated by the decreasing interaction energies. In addition to the loss of inhibitor/enzyme binding, it is interesting to observe that the mutation leads significantly to an increase of the stability of the enzyme.

Hydrophobic sliding: a possible mechanism for drug resistance in human immunodeficiency virus type 1 protease

Hydrophobic residues outside the active site of HIV-1 protease frequently mutate in patients undergoing protease inhibitor therapy; however, the mechanism by which these mutations confer drug resistance is not understood. From analysis of molecular dynamics simulations, 19 core hydrophobic residues appear to facilitate the conformational changes that occur in HIV-1 protease. The hydrophobic core residues slide by each other, exchanging one hydrophobic van der Waal contact for another, with little energy penalty, while maintaining many structurally important hydrogen bonds. Such hydrophobic sliding may represent a general mechanism by which proteins undergo conformational changes. Mutation of these residues in HIV-1 protease would alter the packing of the hydrophobic core, affecting the conformational flexibility of the protease. Therefore these residues impact the dynamic balance between processing substrates and binding inhibitors, and thus contribute to drug resistance.

Emergence of human immunodeficiency virus type 1 variants containing the Q151M complex in children receiving long-term antiretroviral chemotherapy

We examined 28 children with HIV-1 infection who were not responding to existing antiviral regimens and were enrolled into clinical trials conducted at the National Cancer Institute to receive salvage therapy. In 3 of the 28 patients (10.7%), the Q151M complex amino acid substitutions were identified. The three patients had received nucleoside reverse transcriptase inhibitor (NRTI) monotherapy and/or combination regimens with multiple NRTIs for 4.3-8.6 years prior to the study. Recombinant infectious clones generated by incorporating the RT-encoding region of HIV-1 isolated from patients' plasma samples were highly resistant to zidovudine, didanosine and stavudine, while they were moderately resistant to lamivudine and tenofovir disoproxil fumarate (TDF). TDF-containing regimens reduced HIV-1 viremia in two of the three children carrying the Q151M complex. These data suggest that the Q151M could be prevalent in pediatric patients with long-term NRTI monotherapy and/or dual NRTI regimens and that HAART regimens containing TDF may be meritorious in such patients.

Low-throughput model design of protein folding inhibitors

The stabilization energy of proteins in their native conformation is not distributed uniformly among all the amino acids, but is concentrated in few (short) fragments, fragments which play a key role in the folding process and in the stability of the protein. Peptides displaying the same sequence as these key fragments can compete with the formation of the most important native contacts, destabilizing the protein and thus inhibiting its biological activity. We present an essentially automatic method to individuate such peptidic inhibitors based on a low-throughput screening of the fragments which build the target protein. The efficiency and generality of the method is tested on proteins Src-SH3, G, CI2, and HIV-1-PR with the help of a simplified computational model. In each of the cases studied, we find few peptides displaying strong inhibitory properties, properties which are quite robust with respect to point mutations. The possibility of implementing the method through low-throughput experimental screening of the target protein is discussed.

Genotypic changes in human immunodeficiency virus type 1 protease associated with reduced susceptibility and virologic response to the protease inhibitor tipranavir

Tipranavir is a novel, nonpeptidic protease inhibitor of human immunodeficiency virus type 1 (HIV-1) with activity against clinical HIV-1 isolates from treatment-experienced patients. HIV-1 genotypic and phenotypic data from phase II and III clinical trials of tipranavir with protease inhibitor-experienced patients were analyzed to determine the association of protease mutations with reduced susceptibility and virologic response to tipranavir. Specific protease mutations were identified based on stepwise multiple-regression analyses of phase II study data sets. Validation included analyses of phase III study data sets to determine if the same mutations would be selected and to assess how these mutations contribute to multiple-regression models of tipranavir-related phenotype and of virologic response. A tipranavir mutation score was developed from these analyses, which consisted of a unique string of 16 protease positions and 21 mutations (10V, 13V, 20M/R/V, 33F, 35G, 36I, 43T, 46L, 47V, 54A/M/V, 58E, 69K, 74P, 82L/T, 83D, and 84V). HIV-1 isolates displaying an increasing number of these tipranavir resistance-associated mutations had a reduced phenotypic susceptibility and virologic response to tipranavir. Regression models for predicting virologic response in phase III trials revealed that each point in the tipranavir score was associated with a 0.16-log10 copies/ml-lower virologic response to tipranavir at week 24 of treatment. A lower number of points in the tipranavir score and a greater number of active drugs in the background regimen were predictive of virologic success. These analyses demonstrate that the tipranavir mutation score is a potentially valuable tool for predicting the virologic response to tipranavir in protease inhibitor-experienced patients.

Biological characterization of human immunodeficiency virus type 1 subtype C protease carrying indinavir drug-resistance mutations

Human immunodeficiency virus type 1 subtype C isolates belong to one of the most prevalent strains circulating worldwide and are responsible for the majority of new infections in the sub-Saharan region and other highly populated areas of the globe. In this work, the impact of drug-resistance mutations in the protease gene of subtype C viruses was analysed and compared with that of subtype B counterparts. A series of recombinant subtype C and B viruses was constructed carrying indinavir (IDV)-resistance mutations (M46V, I54V, V82A and L90M) and their susceptibility to six FDA-approved protease inhibitor compounds (amprenavir, indinavir, lopinavir, ritonavir, saquinavir and nelfinavir) was determined. A different impact of these mutations was found when nelfinavir and lopinavir were tested. The IDV drug-resistance mutations in the subtype C protease backbone were retained for a long period in culture without selective pressure when compared with those in subtype B counterparts in washout experiments.

Structural and dynamical properties of different protonated states of mutant HIV-1 protease complexed with the saquinavir inhibitor studied by molecular dynamics simulations

To understand the basis of drug resistance, particularly of the HIV-1 PR, three molecular dynamics (MD) simulations of HIV-1 PR mutant species, G48V, complexed with saquinavir (SQV) in explicit aqueous solution with three protonation states, diprotonation on Asp25 and Asp25' (Di-pro) and monoprotonation on each Asp residue (Mono-25 and Mono-25'). For all three states, H-bonds between saquinavir and HIV-1 PR were formed only in the two regions, flap and active site. It was found that conformation of P2 subsite of SQV in the Mono-25 state differs substantially from the other two states. The rotation about 177 degrees from the optimal structure of the wild type was observed, the hydrogen bond between P2 and the flap residue (Val48) was broken and indirect hydrogen bonds with the three residues (Asp29, Gly27, and Asp30) were found instead. In terms of complexation energies, interaction energy of -37.3 kcal/mol for the Mono-25 state is significantly lower than those of -30.7 and -10.7kcal/mol for the Mono-25' and Di-pro states, respectively. It was found also that protonation at the Asp25 leads to a better arrangement in the catalytic dyad, i.e., the Asp25-Asp25' interaction energy of -8.8 kcal/mol of the Mono-25 is significantly lower than that of -2.6kcal/mol for the Mono-25' state. The above data suggest us to conclude that interaction in the catalytic area should be used as criteria to enhance capability in drug designing and drug screening instead of using the total inhibitor/enzyme interaction.

Design of HIV-1-PR inhibitors that do not create resistance: blocking the folding of single monomers

The main problems found in designing drugs are those of optimizing the drug-target interaction and of avoiding the insurgence of resistance. We suggest a scheme for the design of inhibitors that can be used as leads for the development of a drug and that do not face either of these problems, and then apply it to the case of HIV-1-PR. It is based on the knowledge that the folding of single-domain proteins, such as each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES), stabilized by local contacts among hydrophobic, strongly interacting, and highly conserved amino acids that play a central role in the folding process. Because LES have evolved over many generations to recognize and strongly interact with each other so as to make the protein fold fast and avoid aggregation with other proteins, highly specific (and thus little toxic) as well as effective folding-inhibitor molecules suggest themselves: short peptides (or eventually their mimetic molecules) displaying the same amino acid sequence of that of LES (p-LES). Aside from being specific and efficient, these inhibitors are expected not to induce resistance; in fact, mutations in HIV-1-PR that successfully avoid the action of p-LES imply the destabilization of one or more LES and thus should lead to protein denaturation. Making use of Monte Carlo simulations, we first identify the LES of the HIV-1-PR and then show that the corresponding p-LES peptides act as effective inhibitors of the folding of the protease.

Conservation of functional domains and limited heterogeneity of HIV-1 reverse transcriptase gene following vertical transmission

BACKGROUND

The reverse transcriptase (RT) enzyme of human immunodeficiency virus type 1 (HIV-1) plays a crucial role in the life cycle of the virus by converting the single stranded RNA genome into double stranded DNA that integrates into the host chromosome. In addition, RT is also responsible for the generation of mutations throughout the viral genome, including in its own sequences and is thus responsible for the generation of quasi-species in HIV-1-infected individuals. We therefore characterized the molecular properties of RT, including the conservation of functional motifs, degree of genetic diversity, and evolutionary dynamics from five mother-infant pairs following vertical transmission.

RESULTS

The RT open reading frame was maintained with a frequency of 87.2% in five mother-infant pairs' sequences following vertical transmission. There was a low degree of viral heterogeneity and estimates of genetic diversity in mother-infant pairs' sequences. Both mothers and infants RT sequences were under positive selection pressure, as determined by the ratios of non-synonymous to synonymous substitutions. Phylogenetic analysis of 132 mother-infant RT sequences revealed distinct clusters for each mother-infant pair, suggesting that the epidemiologically linked mother-infant pairs were evolutionarily closer to each other as compared with epidemiologically unlinked mother-infant pairs. The functional domains of RT which are responsible for reverse transcription, DNA polymerization and RNase H activity were mostly conserved in the RT sequences analyzed in this study. Specifically, the active sites and domains required for primer binding, template binding, primer and template positioning and nucleotide recruitment were conserved in all mother-infant pairs' sequences.

CONCLUSION

The maintenance of an intact RT open reading frame, conservation of functional domains for RT activity, preservation of several amino acid motifs in epidemiologically linked mother-infant pairs, and a low degree of genetic variability following vertical transmission is consistent with an indispensable role of RT in HIV-1 replication in infected mother-infant pairs.

Combating susceptibility to drug resistance: lessons from HIV-1 protease

Drug resistance is a major obstacle in modern medicine. However, resistance is rarely considered in drug development and may inadvertently be facilitated, as many designed inhibitors contact residues that can mutate to confer resistance, without significantly impairing function. Contemporary drug design often ignores the detailed atomic basis for function and primarily focuses on disrupting the target's activity, which is necessary but not sufficient for developing a robust drug. In this study, we examine the impact of drug-resistant mutations in HIV-1 protease on substrate recognition and demonstrate that most primary active site mutations do not extensively contact substrates, but are critical to inhibitor binding. We propose a general, structure-based strategy to reduce the probability of drug resistance by designing inhibitors that interact only with those residues that are essential for function.

Use of the l1 norm for selection of sparse parameter sets that accurately predict drug response phenotype from viral genetic sequences

We describe the use of the l1 norm for selection of a sparse set of model parameters that are used in the prediction of viral drug response, based on genetic sequence data of the Human Immunodeficiency Virus (HIV) reverse-transcriptase enzyme. We discuss the use of the l1 norm in the Least Absolute Selection and Shrinkage Operator (LASSO) regression model and the Support Vector Machine model. When tested by cross-validation with laboratory measurements, these models predict viral phenotype, or resistance, in response to Reverse-Transcriptase Inhibitors (RTIs) more accurately than other known models. The l1 norm is the most selective convex function, which sets a large proportion of the parameters to zero and also assures that a single optimal solution will be found, given a particular model formulation and training data set. A statistical model that reliably predicts viral drug response is an important tool in the selection of Anti-Retroviral Therapy. These techniques have general application to modeling phenotype from complex genetic data.

Structural and thermodynamic basis for the binding of TMC114, a next-generation human immunodeficiency virus type 1 protease inhibitor

TMC114, a newly designed human immunodeficiency virus type 1 (HIV-1) protease inhibitor, is extremely potent against both wild-type (wt) and multidrug-resistant (MDR) viruses in vitro as well as in vivo. Although chemically similar to amprenavir (APV), the potency of TMC114 is substantially greater. To examine the basis for this potency, we solved crystal structures of TMC114 complexed with wt HIV-1 protease and TMC114 and APV complexed with an MDR (L63P, V82T, and I84V) protease variant. In addition, we determined the corresponding binding thermodynamics by isothermal titration calorimetry. TMC114 binds approximately 2 orders of magnitude more tightly to the wt enzyme (K(d) = 4.5 x 10(-12) M) than APV (K(d) = 3.9 x 10(-10) M). Our X-ray data (resolution ranging from 2.2 to 1.2 A) reveal strong interactions between the bis-tetrahydrofuranyl urethane moiety of TMC114 and main-chain atoms of D29 and D30. These interactions appear largely responsible for TMC114's very favorable binding enthalpy to the wt protease (-12.1 kcal/mol). However, TMC114 binding to the MDR HIV-1 protease is reduced by a factor of 13.3, whereas the APV binding constant is reduced only by a factor of 5.1. However, even with the reduction in binding affinity to the MDR HIV protease, TMC114 still binds with an affinity that is more than 1.5 orders of magnitude tighter than the first-generation inhibitors. Both APV and TMC114 fit predominantly within the substrate envelope, a property that may be associated with decreased susceptibility to drug-resistant mutations relative to that of first-generation inhibitors. Overall, TMC114's potency against MDR viruses is likely a combination of its extremely high affinity and close fit within the substrate envelope.

Altered replicative capacity of recombinant HIV type 1 containing mutations at codons 26 and 29 of the viral protease

HIV-1 protease activity is essential for viral replication. In spite of the high rates of HIV mutation, the active site of protease (residues 24-29) is a conserved site, where mutations have not been previously described. To determine the effect of mutations at positions T26 and A29 of the viral protease and its viability in recombinant HIV-1 strains, Sup-t1 cells were transfected by electroporation with PCR products from a protease containing the 26X or 29X mutation. These mutations were constructed by mutagenesis site directed with degenerative primers. Both changes modified the replicative capacity of the virus: viruses containing the 26X mutation have delayed replication as compared to the control virus HTLVIII B; the presence of the 29X mutation in the viral protease results in the absence of HIV-1 replication. These findings confirm that this region of the viral protease appears to be necessary for the viability of HIV-1, and it could be a good target for antiviral therapy.

Molecular dynamics simulations of 14 HIV protease mutants in complexes with indinavir

The molecular mechanisms of HIV drug resistance were studied using molecular dynamics simulations of HIV-1 protease complexes with the clinical inhibitor indinavir. One nanosecond molecular dynamics simulations were run for solvated complexes of indinavir with wild type protease, a control variant and 12 drug resistant mutants. The quality of the simulations was assessed by comparison with crystallographic and inhibition data. Molecular mechanisms that contribute to drug resistance include structural stability and affinity for inhibitor. The mutants showed a range of structural variation from 70 to 140% of the wild type protease. The protease affinity for indinavir was estimated by calculating the averaged molecular mechanics interaction energy. A correlation coefficient of 0.96 was obtained with observed inhibition constants for wild type and four mutants. Based on this good agreement, the trends in binding were predicted for the other mutants and discussed in relation to the clinical data for indinavir resistance. [figure: see text]. Poincare map representation for WT protease-indinavir complex. The side chain of Tyr 59 showing the positions of hydrogen atoms.

HIV type 1 pol gene diversity and archived nevirapine resistance mutation in pregnant women in Rwanda

This study aimed to find out whether genetic polymorphisms were present in positions potentially affecting susceptibility to antiretrovirals in non-B subtypes from HIV-1-infected patients in Rwanda. Viral pol gene diversity was investigated by direct sequencing in 43 treatment-naive women. In addition, 10 DNA sequences from uncultured peripheral blood mononuclear cells were analyzed 6 weeks after a single dose of nevirapine (prevention of mother-to-child transmission program). Phylogenetic analyses have shown 34 subtype A1, 6 subtype C, and 2 subtype D strains. In addition, an A/C recombinant between the protease (PR) (subtype A1) and the reverse transcriptase (RT) (subtype C) was identified. In the PR coding region, high numbers of polymorphisms were found, including substitutions in secondary PR resistance sites. PR 35D, 36I, and 37N were always present within subtype A as were PR 93L in subtype C strains. PR 10I/V, 20R, 33F, and 77V were found in subtype A whereas PR 36I was highly prevalent in subtype C strains. The A/C recombinant displayed substitutions related to resistance (PR 10, 33, 36 and RT 118). One nevirapine resistance mutation (RT 181Y/C) was found in proviral DNA after 6 weeks. In conclusion, subtypes A and C are predominant in this cohort in Rwanda. Substitutions similar to secondary protease inhibitor resistance mutations are common before treatment whereas major resistance mutation may be archived after a single dose of nevirapine. Accordingly, the hypothesis of a genetic background effect in non-B strains has to be further addressed in programs of introduction of antivirals in Africa.

Covariation of amino acid positions in HIV-1 protease

We have examined patterns of sequence variability for evidence of linked sequence changes in HIV-1 subtype B protease using translated sequences from protease inhibitor (PI) treated and untreated subjects downloaded from the Stanford HIV RT and Protease Sequence Database (http://hivdb.stanford.edu). The final data set size was 648 sequences from untreated subjects (notx) and 531 for PI-treated subjects (tx). Each subject was uniquely represented by a single sequence. Mutual information was calculated for all pairwise comparisons of positions with nonconsensus amino acids in at least 5% of sequences; significance of pairwise association was assessed using permutation tests. In addition pairs of positions were assessed for linkage by comparing the observed occurrences of amino acid combinations to expected values. The mutual information statistic indicated linkage between nine pairs of sites in the untreated data set (10:93, 12:19, 35:38, 37:41, 62:71, 63:64, 71:77, 71:93, 77:93). Strong statistical support for linkage in the treated data set was seen for 32 pairs, eight involving position 10:7 involving position 71, with the rest being 12:19, 15:77, 20:36, 30:88, 35:36, 35:37, 36:62, 36:77, 46:82, 46:84, 48:54, 48:82, 54:82, 63:64, 63:90, 73:90, 77:93, and 84:90. Most associations were positive, although negative associations were seen for five pairs of interactions. Structural proximity suggests that numerous pairs may interact within a local environment. These interactions include two distinct clusters around 36/77 and 71/93. While some of these interactions may reflect fortuitous linkage in heavily treated subjects with many resistance mutations, others will likely represent important cooperative interactions that are amenable to experimental validation.

Molecular cloning and analysis of full-length genome of HIV type 1 strains prevalent in countries of the former Soviet Union

The HIV-1 epidemic among injecting drug users (IDUs) in countries of the former Soviet Union (FSU) was caused mainly by two HIV-1 variants: subtype A and CRF03-AB. To date only three full-length HIV-I genomes from the FSU have been sequenced: one subtype A from Byelorussia and two CRF03-AB from Russia. We report the full-length genome cloning and analysis of two more HIV-1 strains from the FSU countries (98UA0116 of subtype A and 98BY10443 of CRF03-AB). Isolate 98UA0116 is the second cloned and sequenced full-length HIV-1 genome of subtype A lineage from the FSU, which may be a novel subsubtype within sub-type A. Isolate 98BY10443 is the third full-length HIV-1 genome of CRF03-AB in the world to be cloned and sequenced. Additionally, it is the first CRF03-AB strain discovered in Byelorussia. Cloned genomic sequences of the FSU HIV-1 isolates are being used for the development of a region-specific HIV-1 vaccine.

Elucidation of HIV-1 protease resistance by characterization of interaction kinetics between inhibitors and enzyme variants

The kinetics of the interaction between drug-resistant variants of HIV-1 protease (G48V, V82A, L90M, I84V/L90M, and G48V/V82A/I84V/L90M) and clinically used inhibitors (amprenavir, indinavir, nelfinavir, ritonavir, and saquinavir) were determined using biosensor technology. The enzyme variants were immobilized on a biosensor chip and the association and dissociation rate constants (k(on) and k(off)) and affinities (K(D)) for interactions with inhibitors were determined. A unique interaction kinetic profile was observed for each variant/inhibitor combination. Substitution of single amino acids in the protease primarily resulted in reduced affinity through increased k(off) for the inhibitors. For inhibitors characterized by fast association rates to wild-type protease (ritonavir, amprenavir, and indinavir), additional substitutions resulted in a further reduction of affinity by a combination of decreased k(on) and increased k(off). For inhibitors characterized by slow dissociation rates to wild-type enzyme (saquinavir and nelfinavir), the decrease of affinity conferred by additional mutations was attributed to increased k(off) values. Development of resistance thus appears to be associated with a change of the distinctive kinetic parameter contributing to high affinity. Further inhibitor design should focus on improving the "weak point" of the lead compound, that being either k(on) or k(off).

Development and evaluation of a phenotypic assay monitoring resistance formation to protease inhibitors in HIV-1-infected patients

A novel phenotypic assay, based on recombinant expression of the HIV-1-protease was developed and evaluated; it monitors the formation of resistance to protease inhibitors. The HIV-1 protease-encoding region from the blood sample of patients was amplified, ligated into the expression vector pBD2, and recombinantly expressed in Escherichia coli TG1 cells. The resulting recombinant enzyme was purified by a newly developed one-step acid extraction protocol. The protease activity was determined in presence of five selected HIV protease inhibitors and the 50% inhibitory concentration (IC(50)) to the respective protease inhibitors determined. The degree of resistance was expressed in terms of x-fold increase in IC(50) compared to the IC(50) value of an HIV-1 wild type protease preparation. The established test system showed a reproducible recombinant expression of each individual patients' HIV-1 protease population. Samples of nine clinically well characterised HIV-1-infected patients with varying degrees of resistance were analysed. There was a good correlation between clinical parameters and the results obtained by this phenotypic assay. For the majority of patients a blind genotypic analysis of the patients' protease domain revealed a fair correlation to the results of the phenotypic assay. In a minority of patients our phenotypic results diverged from the genotypic ones. This novel phenotypic assay can be carried out within 8-10 days, and offers a significant advantage in time to the current employed phenotypic tests.

Mutation patterns and structural correlates in human immunodeficiency virus type 1 protease following different protease inhibitor treatments

Although many human immunodeficiency virus type 1 (HIV-1)-infected persons are treated with multiple protease inhibitors in combination or in succession, mutation patterns of protease isolates from these persons have not been characterized. We collected and analyzed 2,244 subtype B HIV-1 isolates from 1,919 persons with different protease inhibitor experiences: 1,004 isolates from untreated persons, 637 isolates from persons who received one protease inhibitor, and 603 isolates from persons receiving two or more protease inhibitors. The median number of protease mutations per isolate increased from 4 in untreated persons to 12 in persons who had received four or more protease inhibitors. Mutations at 45 of the 99 amino acid positions in the protease-including 22 not previously associated with drug resistance-were significantly associated with protease inhibitor treatment. Mutations at 17 of the remaining 99 positions were polymorphic but not associated with drug treatment. Pairs and clusters of correlated (covarying) mutations were significantly more likely to occur in treated than in untreated persons: 115 versus 23 pairs and 30 versus 2 clusters, respectively. Of the 115 statistically significant pairs of covarying residues in the treated isolates, 59 were within 8 A of each other-many more than would be expected by chance. In summary, nearly one-half of HIV-1 protease positions are under selective drug pressure, including many residues not previously associated with drug resistance. Structural factors appear to be responsible for the high frequency of covariation among many of the protease residues. The presence of mutational clusters provides insight into the complex mutational patterns required for HIV-1 protease inhibitor resistance.

Selection of high-level resistance to human immunodeficiency virus type 1 protease inhibitors

Protease inhibitors represent some of the most potent agents available for therapeutic strategies designed to inhibit human immunodeficiency virus type 1 (HIV-1) replication. Under certain circumstances the virus develops resistance to the inhibitor, thereby negating the benefits of this therapy. We have carried out selections for high-level resistance to each of three protease inhibitors (indinavir, ritonavir, and saquinavir) in cell culture. Mutations accumulated over most of the course of the increasing selective pressure. There was significant overlap in the identity of the mutations selected with the different inhibitors, and this gave rise to high levels of cross-resistance. Virus particles from the resistant variants all showed defects in processing at the NC/p1 protease cleavage site in Gag. Selections with pairs of inhibitors yielded similar patterns of resistance mutations. A virus that could replicate at near-toxic levels of the three protease inhibitors combined was selected. The pro sequence of this virus was similar to that of the viruses that had been selected for high-level resistance to each of the drugs singly. Finally, a molecular clone carrying the eight most common resistance mutations seen in these selections was characterized. The sequence of this virus was relatively stable during selection for revertants in spite of displaying poor processing at the NC/p1 site and having significantly reduced fitness. These results reveal patterns of drug resistance that extend to near the limits of attainable selective pressure with these inhibitors and confirm the patterns of cross-resistance for these three inhibitors and the attenuation of virion protein processing and fitness that accompanies high-level resistance.

HIV-1 subtypes in Luxembourg, 1983-2000

OBJECTIVES

To study the prevalence of HIV-1 subtypes in Luxembourg between 1983 and 2000. To compare the drug susceptibility of non-B and B clade viruses and the prevalence of resistance-associated mutations and polymorphisms before antiretroviral treatment.

DESIGN

A retrospective study on plasma samples of HIV-infected patients registered at the National Service of Infectious Diseases, Luxembourg, between 1983 and 2000.

METHODS

Genotyping was performed by sequencing of the reverse transcriptase (RT) and protease coding region of the pol gene. Drug susceptibility was assessed in a recombinant virus assay.

RESULTS

A total of 20.1% of the HIV-positive patients were infected with non-B subtypes, and since 1990 the proportion of non-B viruses has increased ninefold. Eleven out of 14 F1 subtypes occurred in patients native to Luxembourg. Major resistance mutations related to protease inhibitors (PI), nucleoside reverse transcriptase inhibitors (NRTI) and non-nucleoside reverse transcriptase inhibitors (NNRTI) occurred in less than 3% of treatment-naive viruses; however, 87% of the viruses had at least one PI-associated mutation. Natural polymorphism of the protease and RT coding region was observed more frequently among non-B than B viruses. Significantly more B viruses displayed resistance to the tested PI, NRTI and NNRTI (P = 0.044).

CONCLUSION

The proportion of non-B viruses has increased dramatically since 1990. Non-B subtypes showed no decreased susceptibility to antiretroviral drugs, but displayed minor mutations and polymorphisms at higher frequency in their protease and RT coding region. In contrast, a significantly higher proportion of B viruses showed resistance to a range of antiretroviral drugs.

Catalytic efficiency and phenotype of HIV-1 proteases encoding single critical resistance substitutions

We have shown that a bacteriophage lambda genetic screen system may be useful in predicting the activity and phenotype of HIV-1 protease in the course of viral infection and antiretroviral therapy. This simple and rapid genetic screening system has been used here to characterize HIV-1 proteases encoding single primary resistance substitutions. Except for proteases with amino acid changes at positions 46 and 84, proteases containing single-resistance substitutions displayed a lower catalytic efficiency than the WT enzyme. Single mutants could be identified by their efficiency, demonstrating that modest differences in protease activity can be monitored with this simple assay. Overall, drug susceptibility could be reduced by introduction of single mutations. However, high-level protease inhibitor (PI) resistance was only achieved by multiple mutated proteases. The small but reproducible increase in resistance displayed by single mutants also demonstrated the ability of this genetic screen system for detecting minor reductions in drug susceptibility. These results show that the bacteriophage lambda genetic screen system used here is a useful tool in the analysis of specific contribution of mutations in the HIV protease-coding region or in specific cleavage sites that affect the process of PI resistance.

Changes in human immunodeficiency virus type 1 Gag at positions L449 and P453 are linked to I50V protease mutants in vivo and cause reduction of sensitivity to amprenavir and improved viral fitness in vitro

Human immunodeficiency virus type 1 (HIV-1) Gag protease cleavage sites (CS) undergo sequence changes during the development of resistance to several protease inhibitors (PIs). We have analyzed the association of sequence variation at the p7/p1 and p1/p6 CS in conjunction with amprenavir (APV)-specific protease mutations following PI combination therapy with APV. Querying a central resistance data repository resulted in the detection of significant associations (P < 0.001) between the presence of APV protease signature mutations and Gag L449F (p1/p6 LP1'F) and P453L (p1/p6 PP5'L) CS changes. In population-based sequence analyses the I50V mutant was invariably linked to either L449F or P453L. Clonal analysis revealed that both CS mutations were never present in the same genome. Sequential plasma samples from one patient revealed a transition from I50V M46L P453L viruses at early time points to I50V M46I L449F viruses in later samples. Various combinations of the protease and Gag mutations were introduced into the HXB2 laboratory strain of HIV-1. In both single- and multiple-cycle assay systems and in the context of I50V, the L449F and P453L changes consistently increased the 50% inhibitory concentration of APV, while the CS changes alone had no measurable effect on inhibitor sensitivity. The decreased in vitro fitness of the I50V mutant was only partially improved by addition of either CS change (I50V M46I L449F mutant replicative capacity approximately 16% of that of wild-type virus). Western blot analysis of Pr55 Gag precursor cleavage products from infected-cell cultures indicated accumulation of uncleaved Gag p1-p6 in all I50V viruses without coexisting CS changes. Purified I50V protease catalyzed cleavage of decapeptides incorporating the L449F or P453L change 10-fold and 22-fold more efficiently than cleavage of the wild-type substrate, respectively. HIV-1 protease CS changes are selected during PI therapy and can have effects on both viral fitness and phenotypic resistance to PIs.

Genotypic Correlates of Resistance to HIV-1 Protease Inhibitors on Longitudinal Data: The Role of Secondary Mutations

Direct sequencing of the pol gene was assessed retrospectively with protease inhibitor susceptibility in a longitudinal study. A total of 134 samples from 26 patients were analysed at regular intervals up to 2 years. Patients were included in virological failure despite indinavir, ritonavir or saquinavir based triple-drug therapy. Both the type and number of certain secondary protease mutations modulated the effect of primary mutations on phenotypic resistance. This was notably applicable to L10I/V, and to lesser extents to A71V/T. However, combinations of primary mutations, including I54V could predict resistance to the drug used and nelfinavir in more than 80%. In contrast, in vitro cross-resistance to amprenavir was rarely encountered. In addition, there was a relationship between a higher number of key mutations and poorer virological and clinical outcomes, respectively, from 6 and 3 months on. The key mutations were the protease mutations independently conferring phenotypic resistance and/or the reverse transcriptase mutations predicting treatment outcome. This relationship was independent from drug history, viral load and CD4 cell count measurements. In summary, even on a small sample size, sequence-based genotyping seems to be a good prognostic marker when performed longitudinally. In the context of primary resistance mutations, including additional secondary mutations, it may be useful in the prediction of phenotypic and clinical resistance. This should be assessed to optimize treatment monitoring before emergence of broadly cross-resistant virus.

Rapid and sensitive oligonucleotide ligation assay for detection of mutations in human immunodeficiency virus type 1 associated with high-level resistance to protease inhibitors

A sensitive, specific, and high-throughput oligonucleotide ligation assay (OLA) for the detection of genotypic human immunodeficiency virus type 1 (HIV-1) resistance to Food and Drug Administration-approved protease inhibitors was developed and evaluated. This ligation-based assay uses differentially modified oligonucleotides specific for wild-type or mutant sequences, allowing sensitive and simple detection of both genotypes in a single well of a microtiter plate. Oligonucleotides were designed to detect primary mutations associated with high-level resistance to amprenavir, nelfinavir, indinavir, ritonavir, saquinavir, and lopinavir, including amino acid substitutions D30N, I50V, V82A/S/T, I84V, N88D, and L90M. Plasma HIV-1 RNA from 54 infected patients was amplified by reverse transcription-PCR and sequenced by using dideoxynucleotide chain terminators for evaluation of mutations associated with drug resistance. These same amplicons were genotyped by the OLA at positions 30, 50, 82, 88, 84, and 90 for a total of 312 codons. The sensitivity of detection of drug-resistant genotypes was 96.7% (87 of 90 mutant codons) in the OLA compared to 92.2% (83 of 90) in consensus sequencing, presumably due to the increased sensitivity of the OLA. The OLA detected genetic subpopulations more often than sequencing, detecting 30 mixtures of mutant and wild-type sequences and two mixtures of drug-resistant sequences compared to 15 detected by DNA sequencing. Reproducible and semiquantitative detection of the mutant and the wild-type genomes by the OLA was observed by analysis of wild-type and mutant plasmid mixtures containing as little as 5% of either genotype in a background of the opposite genome. This rapid, simple, economical, and highly sensitive assay provides a practical alternative to dideoxy sequencing for genotypic evaluation of HIV-1 resistance to antiretrovirals.

Genotypic testing for human immunodeficiency virus type 1 drug resistance

There are 16 approved human immunodeficiency virus type 1 (HIV-1) drugs belonging to three mechanistic classes: protease inhibitors, nucleoside and nucleotide reverse transcriptase (RT) inhibitors, and nonnucleoside RT inhibitors. HIV-1 resistance to these drugs is caused by mutations in the protease and RT enzymes, the molecular targets of these drugs. Drug resistance mutations arise most often in treated individuals, resulting from selective drug pressure in the presence of incompletely suppressed virus replication. HIV-1 isolates with drug resistance mutations, however, may also be transmitted to newly infected individuals. Three expert panels have recommended that HIV-1 protease and RT susceptibility testing should be used to help select HIV drug therapy. Although genotypic testing is more complex than typical antimicrobial susceptibility tests, there is a rich literature supporting the prognostic value of HIV-1 protease and RT mutations. This review describes the genetic mechanisms of HIV-1 drug resistance and summarizes published data linking individual RT and protease mutations to in vitro and in vivo resistance to the currently available HIV drugs.

Analysis of HIV drug-resistant quasispecies in plasma, peripheral blood mononuclear cells and viral isolates from treatment-naive and HAART patients

The pattern of HIV-1 reverse transcriptase and protease mutations conferring resistance to antiretroviral drugs was studied in five treatment-naive patients and five HIV-infected patients receiving HAART [two reverse transcriptase inhibitors + one protease inhibitor] for > or = 1 year. Direct sequencing was performed on plasma HIV RNA, HIV DNA from peripheral blood mononuclear cells (PBMCs), and RNA from viral isolates. In addition, reverse transcriptase and protease PCR products from PBMCs HIV DNA, plasma HIV RNA, and viral isolate RNA were cloned in a plasmid to study the quasispecies distribution of drug-resistance associated mutations. Direct sequencing of HIV DNA from PBMCs and HIV RNA from plasma and viral isolates did not show the presence of drug resistance associated mutations in both reverse transcriptase and protease of HIV from all five treatment-naive patients. On the contrary, mutation analysis obtained by cloning plasma HIV RNA and PBMCs DNA showed the presence of drug-resistance related mutations at a low frequency in both HIV enzymes of four out of five treatment-naive patients. On the other hand, direct sequencing of plasma HIV RNA showed the presence of several reverse transcriptase and protease mutations in all five treated patients. Mutation analysis performed by cloning PBMCs HIV DNA, and HIV RNA from plasma and viral isolates, revealed additional reverse transcriptase and protease changes compared to direct sequencing of the relevant biological samples. All the additional changes were observed in a minority of clones. In conclusion, the data suggest that less frequent drug-resistant viral variants not detected by direct sequencing of PBMCs, plasma samples, or viral isolates are present in both treatment-naive and treatment-experienced HIV patients. These findings may have important implications in the understanding of the selection process of drug-resistant variants under drug pressure.

Viral evolution in response to the broad-based retroviral protease inhibitor TL-3

TL-3 is a protease inhibitor developed using the feline immunodeficiency virus protease as a model. It has been shown to efficiently inhibit replication of human, simian, and feline immunodeficiency viruses and therefore has broad-based activity. We now demonstrate that TL-3 efficiently inhibits the replication of 6 of 12 isolates with confirmed resistance mutations to known protease inhibitors. To dissect the spectrum of molecular changes in protease and viral properties associated with resistance to TL-3, a panel of chronological in vitro escape variants was generated. We have virologically and biochemically characterized mutants with one (V82A), three (M46I/F53L/V82A), or six (L24I/M46I/F53L/L63P/V77I/V82A) changes in the protease and structurally modeled the protease mutant containing six changes. Virus containing six changes was found to be 17-fold more resistant to TL-3 in cell culture than was wild-type virus but maintained similar in vitro replication kinetics compared to the wild-type virus. Analyses of enzyme activity of protease variants with one, three, and six changes indicated that these enzymes, compared to wild-type protease, retained 40, 47, and 61% activity, respectively. These results suggest that deficient protease enzymatic activity is sufficient for function, and the observed protease restoration might imply a selective advantage, at least in vitro, for increased protease activity.

Rapid genetic selection of inhibitor-resistant protease mutants: clinically relevant and novel mutants of the HIV protease

Site-specific proteolysis is an important regulatory mechanism in many basic cellular processes as well as playing critical roles in the life cycle of viruses and other pathogenic organisms. For the human immunodeficiency virus (HIV) the encoded protease is required for the replication of the virus and has been the target of novel antiviral therapeutics. However, the emergence of inhibitor-resistant viral strains has become an increasingly significant clinical problem. Using a bacteriophage-based genetic selection, a bank of inhibitor-resistant mutants in this protease has been isolated that includes mutations that correlate with resistant clinical isolates. The rapid selection of such mutations has implications for the prediction of relevant mutations and may be applicable to other viral systems.

International perspectives on antiretroviral resistance. Resistance to protease inhibitors

The availability of protease inhibitors (PIs) and their combination with nucleoside reverse transcriptase inhibitors marked the passage of antiretroviral therapy (ART) from potential for control to effective suppression and thus substantially reduced rates of morbidity and mortality related to HIV. Even so, what was first hoped to be an immutable HIV DNA treatment target has proved to be prone to resistance mutations, with substitutions identified at more than 20 amino acid sites, which reduces PI susceptibility and increases resistance to treatment. The mutation patterns associated with each PI have been defined, and have been observed to occur at one of two locations: at or near the active site, or in the substrate cleavage site. The natural history of PI resistance has been extensively studied, and the genetic and cellular pathways are described in detail in this article. In addition, cross-resistance among PIs is now recognized to be fairly extensive, although the degree of cross-resistance varies with the number of mutations and the variants selected by drug pressure. Thus, it is still possible to salvage a response with another PI after a first regimen with another PI has failed. The extensive basic science and clinical experience with PIs in the fight against HIV are reviewed in this article, which provides data on resistance-mutation profiles, cellular resistance mechanisms, viral fitness studies, and clinical outcome trials with various first-line and subsequent regimens that contain PIs. It is hoped that the information provided will guide physicians in best using PIs as part of a logical and successful ART strategy.

Variant human immunodeficiency virus type 1 proteases and response to combination therapy including a protease inhibitor

The objective of this observational study was to assess the genetic variability in the human immunodeficiency virus (HIV) protease gene from HIV type 1 (HIV-1)-positive (clade B), protease inhibitor-naïve patients and to evaluate its association with the subsequent effectiveness of a protease inhibitor-containing triple-drug regimen. The protease gene was sequenced from plasma-derived virus from 116 protease inhibitor-naïve patients. The virological response to a triple-drug regimen containing indinavir, ritonavir, or saquinavir was evaluated every 3 months for as long as 2 years (n = 40). A total of 36 different amino acid substitutions compared to the reference sequence (HIV-1 HXB2) were detected. No substitutions at the active site similar to the primary resistance mutations were found. The most frequent substitutions (prevalence, >10%) at baseline were located at codons 15, 13, 12, 62, 36, 64, 41, 35, 3, 93, 77, 63, and 37 (in ascending order of frequency). The mean number of polymorphisms was 4.2. A relatively poorer response to therapy was associated with a high number of baseline polymorphisms and, to a lesser extent, with the presence of I93L at baseline in comparison with the wild-type virus. A71V/T was slightly associated with a poorer response to first-line ritonavir-based therapy. In summary, within clade B viruses, protease gene natural polymorphisms are common. There is evidence suggesting that treatment response is associated with this genetic background, but most of the specific contributors could not be firmly identified. I93L, occurring in about 30% of untreated patients, may play a role, as A71V/T possibly does in ritonavir-treated patients.

Curling of flap tips in HIV-1 protease as a mechanism for substrate entry and tolerance of drug resistance

BACKGROUND

The human immunodeficiency virus type 1 (HIV-1) protease is an essential viral protein that is a major drug target in the fight against Acquired Immune Deficiency Syndrome (AIDS). Access to the active site of this homodimeric enzyme is gained when two large flaps, one from each monomer, open. The flap movements are therefore central to the function of the enzyme, yet determining how these flaps move at an atomic level has not been experimentally possible.

RESULTS

In the present study, we observe the flaps of HIV-1 protease completely opening during a 10 ns solvated molecular dynamics simulation starting from the unliganded crystal structure. This movement is on the time scale observed by Nuclear Magnetic Resonance (NMR) relaxation data. The highly flexible tips of the flaps, with the sequence Gly-Gly-Ile-Gly-Gly, are seen curling back into the protein and thereby burying many hydrophobic residues.

CONCLUSIONS

This curled-in conformational change has never been previously described. Previous models of this movement, with the flaps as rigid levers, are not consistent with the experimental data. The residues that participate in this hydrophobic cluster as a result of the conformational change are highly sensitive to mutation and often contribute to drug resistance when they do change. However, several of these residues are not part of the active site cavity, and their essential role in causing drug resistance could possibly be rationalized if this conformational change actually occurs. Trapping HIV-1 protease in this inactive conformation would provide a unique opportunity for future drug design.