Increasing prevalence of HIV-1 protease inhibitor-associated mutations correlates with long-term non-suppressive protease inhibitor treatment

Deep sequencing of protease inhibitor resistant HIV patient isolates reveals patterns of correlated mutations in Gag and protease

While the role of drug resistance mutations in HIV protease has been studied comprehensively, mutations in its substrate, Gag, have not been extensively cataloged. Using deep sequencing, we analyzed a unique collection of longitudinal viral samples from 93 patients who have been treated with therapies containing protease inhibitors (PIs). Due to the high sequence coverage within each sample, the frequencies of mutations at individual positions were calculated with high precision. We used this information to characterize the variability in the Gag polyprotein and its effects on PI-therapy outcomes. To examine covariation of mutations between two different sites using deep sequencing data, we developed an approach to estimate the tight bounds on the two-site bivariate probabilities in each viral sample, and the mutual information between pairs of positions based on all the bounds. Utilizing the new methodology we found that mutations in the matrix and p6 proteins contribute to continued therapy failure and have a major role in the network of strongly correlated mutations in the Gag polyprotein, as well as between Gag and protease. Although covariation is not direct evidence of structural propensities, we found the strongest correlations between residues on capsid and matrix of the same Gag protein were often due to structural proximity. This suggests that some of the strongest inter-protein Gag correlations are the result of structural proximity. Moreover, the strong covariation between residues in matrix and capsid at the N-terminus with p1 and p6 at the C-terminus is consistent with residue-residue contacts between these proteins at some point in the viral life cycle.

Understanding the structure of HIV proteins and the function of drug-resistant mutations of these proteins is critical for the development of effective HIV treatments. Selected gag mutations have been shown to provide compensatory functions for protease resistance mutations and may directly contribute to the development of drug resistance. To determine associations between protease inhibitor mutations and gag, we utilized deep sequencing of HIV gag and protease from a collection of viral isolates from patients treated with highly active retroviral protease inhibitors. Deep sequencing allows for accurate measurement of mutation frequencies at each position, allowing estimation, using a novel method we developed, of the covariation between any two residues on gag. Using this information, we characterize the variation within gag and protease and identify the most strongly correlated pairs of inter- and intra-protein residues. Our results suggest that matrix and p1/p6 mutations form the core of a network of strongly correlated gag mutations and contribute to recurrent treatment failure. Extracting gag residue covariation information from the deep sequencing of patient viral samples may provide insight into structural aspects of the Gag polyprotein as well new areas for small molecule targeting to disrupt Gag function.

A multifaceted analysis of HIV-1 protease multidrug resistance phenotypes

Background

Great strides have been made in the effective treatment of HIV-1 with the development of second-generation protease inhibitors (PIs) that are effective against historically multi-PI-resistant HIV-1 variants. Nevertheless, mutation patterns that confer decreasing susceptibility to available PIs continue to arise within the population. Understanding the phenotypic and genotypic patterns responsible for multi-PI resistance is necessary for developing PIs that are active against clinically-relevant PI-resistant HIV-1 variants.

Results

In this work, we use globally optimal integer programming-based clustering techniques to elucidate multi-PI phenotypic resistance patterns using a data set of 398 HIV-1 protease sequences that have each been phenotyped for susceptibility toward the nine clinically-approved HIV-1 PIs. We validate the information content of the clusters by evaluating their ability to predict the level of decreased susceptibility to each of the available PIs using a cross validation procedure. We demonstrate the finding that as a result of phenotypic cross resistance, the considered clinical HIV-1 protease isolates are confined to ~6% or less of the clinically-relevant phenotypic space. Clustering and feature selection methods are used to find representative sequences and mutations for major resistance phenotypes to elucidate their genotypic signatures. We show that phenotypic similarity does not imply genotypic similarity, that different PI-resistance mutation patterns can give rise to HIV-1 isolates with similar phenotypic profiles.

Conclusion

Rather than characterizing HIV-1 susceptibility toward each PI individually, our study offers a unique perspective on the phenomenon of PI class resistance by uncovering major multidrug-resistant phenotypic patterns and their often diverse genotypic determinants, providing a methodology that can be applied to understand clinically-relevant phenotypic patterns to aid in the design of novel inhibitors that target other rapidly evolving molecular targets as well.

Five-year follow up of genotypic resistance patterns in HIV-1 subtype C infected patients in Botswana after failure of thymidine analogue-based regimens

Objective

Our objective was to establish genotypic resistance profiles among the 4% of Batswana patients who experienced virologic failure while being followed within Botswana's National Antiretroviral Treatment Program between 2002 and 2007.

Methods

At the beginning of the national program in 2002, almost all patients received stavudine (d4T), together with didanosine (ddI), as part of their first nucleoside reverse transcriptase inhibitor (NRTI)-based regimen (Group 1). In contrast, the standard of care for all patients subsequently enrolled (2002-2007) included zidovudine/lamivudine (ZDV/3TC) (Group 2). Genotypes were analyzed in 26 patients from Group 1 and 37 patients from Group 2. Associations between mutations were determined using Pearson's correlation coefficient and Jaccard's coefficient of similarity.

Results

Seventy-eight percent of genotyped patients possessed mutations associated with protease inhibitor (PI) resistance while 87% and 90%, respectively, exhibited mutations associated with NRTIs and non-nucleoside reverse transcriptase inhibitors (NNRTIs). The most frequent PI mutations involving resistance to NFV were L90M (25.2%) and D30N (16.2%), but mutations at positions K45Q and D30N were often observed in tandem (P = 60.5, J = 50; p = 0.002; Group 2) alongside Q61E in 42.8% of patients who received ZDV/3TC. Both major patterns of thymidine analogue mutations, TAM 1 (48%) and TAM 2 (59%), were represented in patients from Group 1 and 2, although M184V was higher among individuals who had initially received ddI (61% versus 40.5%). In contrast, L74V was more frequent among individuals from Group 2 (16.2% versus 7.7%). Differences in regard to NNRTI mutations were also observed between Group 1 and Group 2 patients.

Conclusion

Despite a low rate of therapeutic failure (4%) among these patients, those who failed possessed high numbers of resistance mutations as well as novel resistance mutations and/or polymorphisms at sites within reverse transcriptase and protease.

Nonpolymorphic human immunodeficiency virus type 1 protease and reverse transcriptase treatment-selected mutations

The spectrum of human immunodeficiency virus type 1 (HIV-1) protease and reverse transcriptase (RT) mutations selected by antiretroviral (ARV) drugs requires ongoing reassessment as ARV treatment patterns evolve and increasing numbers of protease and RT sequences of different viral subtypes are published. Accordingly, we compared the prevalences of protease and RT mutations in HIV-1 group M sequences from individuals with and without a history of previous treatment with protease inhibitors (PIs) or RT inhibitors (RTIs). Mutations in protease sequences from 26,888 individuals and in RT sequences from 25,695 individuals were classified according to whether they were nonpolymorphic in untreated individuals and whether their prevalence increased fivefold with ARV therapy. This analysis showed that 88 PI-selected and 122 RTI-selected nonpolymorphic mutations had a prevalence that was fivefold higher in individuals receiving ARVs than in ARV-naïve individuals. This was an increase of 47% and 77%, respectively, compared with the 60 PI- and 69 RTI-selected mutations identified in a similar analysis that we published in 2005 using subtype B sequences obtained from one-fourth as many individuals. In conclusion, many nonpolymorphic mutations in protease and RT are under ARV selection pressure. The spectrum of treatment-selected mutations is changing as data for more individuals are collected, treatment exposures change, and the number of available sequences from non-subtype B viruses increases.

Molecular characterization of clinical isolates of human immunodeficiency virus resistant to the protease inhibitor darunavir

Darunavir is the most recently approved human immunodeficiency virus (HIV) protease (PR) inhibitor (PI) and is active against many HIV type 1 PR variants resistant to earlier-generation PIs. Darunavir shows a high genetic barrier to resistance development, and virus strains with lower sensitivity to darunavir have a higher number of PI resistance-associated mutations than viruses resistant to other PIs. In this work, we have enzymologically and structurally characterized a number of highly mutated clinically derived PRs with high levels of phenotypic resistance to darunavir. With 18 to 21 amino acid residue changes, the PR variants studied in this work are the most highly mutated HIV PR species ever studied by means of enzyme kinetics and X-ray crystallography. The recombinant proteins showed major defects in substrate binding, while the substrate turnover was less affected. Remarkably, the overall catalytic efficiency of the recombinant PRs (5% that of the wild-type enzyme) is still sufficient to support polyprotein processing and particle maturation in the corresponding viruses. The X-ray structures of drug-resistant PRs complexed with darunavir suggest that the impaired inhibitor binding could be explained by change in the PR-inhibitor hydrogen bond pattern in the P2' binding pocket due to a substantial shift of the aminophenyl moiety of the inhibitor. Recombinant virus phenotypic characterization, enzyme kinetics, and X-ray structural analysis thus help to explain darunavir resistance development in HIV-positive patients.

Alternative methods to analyse the impact of HIV mutations on virological response to antiviral therapy

Background

Principal component analysis (PCA) and partial least square (PLS) regression may be useful to summarize the HIV genotypic information. Without pre-selection each mutation presented in at least one patient is considered with a different weight. We compared these two strategies with the construction of a usual genotypic score.

Methods

We used data from the ANRS-CO3 Aquitaine Cohort Zephir sub-study. We used a subset of 87 patients with a complete baseline genotype and plasma HIV-1 RNA available at baseline and at week 12. PCA and PLS components were determined with all mutations that had prevalences >0. For the genotypic score, mutations were selected in two steps: 1) p-value < 0.01 in univariable analysis and prevalences between 10% and 90% and 2) backwards selection procedure based on the Cochran-Armitage Test. The predictive performances were compared by means of the cross-validated area under the receiver operating curve (AUC).

Results

Virological failure was observed in 46 (53%) patients at week 12. Principal components and PLS components showed a good performance for the prediction of virological response in HIV infected patients. The cross-validated AUCs for the PCA, PLS and genotypic score were 0.880, 0.868 and 0.863, respectively. The strength of the effect of each mutation could be considered through PCA and PLS components. In contrast, each selected mutation contributes with the same weight for the calculation of the genotypic score. Furthermore, PCA and PLS regression helped to describe mutation clusters (e.g. 10, 46, 90).

Conclusion

In this dataset, PCA and PLS showed a good performance but their predictive ability was not clinically superior to that of the genotypic score.

Enzymatic and structural analysis of the I47A mutation contributing to the reduced susceptibility to HIV protease inhibitor lopinavir

Lopinavir (LPV) is a second-generation HIV protease inhibitor (PI) designed to overcome resistance development in patients undergoing long-term antiviral therapy. The mutation of isoleucine at position 47 of the HIV protease (PR) to alanine is associated with a high level of resistance to LPV. In this study, we show that recombinant PR containing a single I47A substitution has the inhibition constant (K(i) ) value for lopinavir by two orders of magnitude higher than for the wild-type PR. The addition of the I47A substitution to the background of a multiply mutated PR species from an AIDS patient showed a three-order-of-magnitude increase in K(i) in vitro relative to the patient PR without the I47A mutation. The crystal structure of I47A PR in complex with LPV showed the loss of van der Waals interactions in the S2/S2' subsites. This is caused by the loss of three side-chain methyl groups due to the I47A substitution and by structural changes in the A47 main chain that lead to structural changes in the flap antiparallel beta-strand. Furthermore, we analyzed possible interaction of the I47A mutation with secondary mutations V32I and I54V. We show that both mutations in combination with I47A synergistically increase the relative resistance to LPV in vitro. The crystal structure of the I47A/I54V PR double mutant in complex with LPV shows that the I54V mutation leads to a compaction of the flap, and molecular modeling suggests that the introduction of the I54V mutation indirectly affects the strain of the bound inhibitor in the PR binding cleft.

Ninety-nine is not enough: molecular characterization of inhibitor-resistant human immunodeficiency virus type 1 protease mutants with insertions in the flap region

While the selection of amino acid insertions in human immunodeficiency virus (HIV) reverse transcriptase (RT) is a known mechanism of resistance against RT inhibitors, very few reports on the selection of insertions in the protease (PR) coding region have been published. It is still unclear whether these insertions impact protease inhibitor (PI) resistance and/or viral replication capacity. We show that the prevalence of insertions, especially between amino acids 30 to 41 of HIV type 1 (HIV-1) PR, has increased in recent years. We identified amino acid insertions at positions 33 and 35 of the PR of HIV-1-infected patients who had undergone prolonged treatment with PIs, and we characterized the contribution of these insertions to viral resistance. We prepared the corresponding mutated, recombinant PR variants with or without insertions at positions 33 and 35 and characterized them in terms of enzyme kinetics and crystal structures. We also engineered the corresponding recombinant viruses and analyzed the PR susceptibility and replication capacity by recombinant virus assay. Both in vitro methods confirmed that the amino acid insertions at positions 33 and 35 contribute to the viral resistance to most of the tested PIs. The structural analysis revealed local structural rearrangements in the flap region and in the substrate binding pockets. The enlargement of the PR substrate binding site together with impaired flap dynamics could account for the weaker inhibitor binding by the insertion mutants. Amino acid insertions in the vicinity of the binding cleft therefore represent a novel mechanism of HIV resistance development.

Identification of a novel resistance (E40F) and compensatory (K43E) substitution in HIV-1 reverse transcriptase

Background

HIV-1 nucleoside reverse transcriptase inhibitors (NRTIs) have been used in the clinic for over twenty years. Interestingly, the complete resistance pattern to this class has not been fully elucidated. Novel mutations in RT appearing during treatment failure are still being identified. To unravel the role of two of these newly identified changes, E40F and K43E, we investigated their effect on viral drug susceptibility and replicative capacity.

Results

A large database (Quest Diagnostics database) was analysed to determine the associations of the E40F and K43E changes with known resistance mutations. Both amino acid changes are strongly associated with the well known NRTI-resistance mutations M41L, L210W and T215Y. In addition, a strong positive association between these changes themselves was observed. A panel of recombinant viruses was generated by site-directed mutagenesis and phenotypically analysed. To determine the effect on replication capacity, competition and in vitro evolution experiments were performed. Introduction of E40F results in an increase in Zidovudine resistance ranging from nine to fourteen fold depending on the RT background and at the same time confers a decrease in viral replication capacity. The K43E change does not decrease the susceptibility to Zidovudine but increases viral replication capacity, when combined with E40F, demonstrating a compensatory role for this codon change.

Conclusion

In conclusion, we have identified a novel resistance (E40F) and compensatory (K43E) change in HIV-1 RT. Further research is indicated to analyse the clinical importance of these changes.

HIV-1 subtype B protease and reverse transcriptase amino acid covariation

Despite the high degree of HIV-1 protease and reverse transcriptase (RT) mutation in the setting of antiretroviral therapy, the spectrum of possible virus variants appears to be limited by patterns of amino acid covariation. We analyzed patterns of amino acid covariation in protease and RT sequences from more than 7,000 persons infected with HIV-1 subtype B viruses obtained from the Stanford HIV Drug Resistance Database (http://hivdb.stanford.edu). In addition, we examined the relationship between conditional probabilities associated with a pair of mutations and the order in which those mutations developed in viruses for which longitudinal sequence data were available. Patterns of RT covariation were dominated by the distinct clustering of Type I and Type II thymidine analog mutations and the Q151M-associated mutations. Patterns of protease covariation were dominated by the clustering of nelfinavir-associated mutations (D30N and N88D), two main groups of protease inhibitor (PI)-resistance mutations associated either with V82A or L90M, and a tight cluster of mutations associated with decreased susceptibility to amprenavir and the most recently approved PI darunavir. Different patterns of covariation were frequently observed for different mutations at the same position including the RT mutations T69D versus T69N, L74V versus L74I, V75I versus V75M, T215F versus T215Y, and K219Q/E versus K219N/R, and the protease mutations M46I versus M46L, I54V versus I54M/L, and N88D versus N88S. Sequence data from persons with correlated mutations in whom earlier sequences were available confirmed that the conditional probabilities associated with correlated mutation pairs could be used to predict the order in which the mutations were likely to have developed. Whereas accessory nucleoside RT inhibitor-resistance mutations nearly always follow primary nucleoside RT inhibitor-resistance mutations, accessory PI-resistance mutations often preceded primary PI-resistance mutations.