Role of Gag mutations in PI resistance in the Swiss HIV cohort study: bystanders or contributors?

HIV-1 envelope facilitates the development of protease inhibitor resistance through acquiring mutations associated with viral entry and immune escape

Introduction

There is increasing evidence supporting a role for HIV-1 envelope in the development of Protease Inhibitor drug resistance, and a recent report from our group suggested that Env mutations co-evolve with Gag-Protease mutations in the pathway to Lopinavir resistance. In this study, we investigated the effect of co-evolving Env mutations on virus function and structure.

Methods

Co-receptor usage and n-linked glycosylation were investigated using Geno2Pheno as well as tools available at the Los Alamos sequence database. Molecular dynamics simulations were performed using Amber 18 and analyzed using Cpptraj, and molecular interactions were calculated using the Ring server.

Results

The results showed that under Protease Inhibitor drug selection pressure, the envelope gene modulates viral entry by protecting the virus from antibody recognition through the increased length and number of N-glycosylation sites observed in V1/V2 and to some extent V5. Furthermore, gp120 mutations appear to modulate viral entry through a switch to the CXCR4 coreceptor, induced by higher charge in the V3 region and specific mutations at the coreceptor binding sites. In gp41, S534A formed a hydrogen bond with L602 found in the disulfide loop region between the Heptad Repeat 1 and Heptad Repeat 2 domains and could negatively affect the association of gp120-gp41 during viral entry. Lastly, P724Q/S formed both intermolecular and intramolecular interactions with residues within the Kennedy loop, a known epitope.

Discussion

In conclusion, the results suggest that mutations in envelope during Protease Inhibitor treatment failure are related to immune escape and that S534A mutants could preferentially use the cell-to-cell route of infection.

HIV-1 Evolutionary Dynamics under Nonsuppressive Antiretroviral Therapy

Prolonged virologic failure on 2nd-line protease inhibitor (PI)-based antiretroviral therapy (ART) without emergence of major protease mutations is well recognized and provides an opportunity to study within-host evolution in long-term viremic individuals. Using next-generation sequencing and in silico haplotype reconstruction, we analyzed whole-genome sequences from longitudinal plasma samples of eight chronically infected HIV-1-positive individuals failing 2nd-line regimens from the French National Agency for AIDS and Viral Hepatitis Research (ANRS) 12249 Treatment as Prevention (TasP) trial. On nonsuppressive ART, there were large fluctuations in synonymous and nonsynonymous variant frequencies despite stable viremia. Reconstructed haplotypes provided evidence for selective sweeps during periods of partial adherence, and viral haplotype competition, during periods of low drug exposure. Drug resistance mutations in reverse transcriptase (RT) were used as markers of viral haplotypes in the reservoir, and their distribution over time indicated recombination. We independently observed linkage disequilibrium decay, indicative of recombination. These data highlight dramatic changes in virus population structure that occur during stable viremia under nonsuppressive ART. IMPORTANCE HIV-1 infections are most commonly initiated with a single founder virus and are characterized by extensive inter- and intraparticipant genetic diversity. However, existing literature on HIV-1 intrahost population dynamics is largely limited to untreated infections, predominantly in subtype B-infected individuals. The manuscript characterizes viral population dynamics in long-term viremic treatment-experienced individuals, which has not been previously characterized. These data are particularly relevant for understanding HIV dynamics but can also be applied to other RNA viruses. With this unique data set we propose that the virus is highly unstable, and we have found compelling evidence of HIV-1 within-host viral diversification, recombination, and haplotype competition during nonsuppressive ART.

Potential Associations of Mutations within the HIV-1 Env and Gag Genes Conferring Protease Inhibitor (PI) Drug Resistance

An increasing number of patients in Africa are experiencing virological failure on a second-line antiretroviral protease inhibitor (PI)-containing regimen, even without resistance-associated mutations in the protease region, suggesting a potential role of other genes in PI resistance. Here, we investigated the prevalence of mutations associated with Lopinavir/Ritonavir (LPV/r) failure in the Envelope gene and the possible coevolution with mutations within the Gag-protease (gag-PR) region. Env and Gag-PR sequences generated from 24 HIV-1 subtype C infected patients failing an LPV/r inclusive treatment regimen and 344 subtype C drug-naïve isolates downloaded from the Los Alamos Database were analyzed. Fisher’s exact test was used to determine the differences in mutation frequency. Bayesian network probability was applied to determine the relationship between mutations occurring within the env and gag-PR regions and LPV/r treatment. Thirty-five mutations in the env region had significantly higher frequencies in LPV/r-treated patients. A combination of Env and Gag-PR mutations was associated with a potential pathway to LPV/r resistance. While Env mutations were not directly associated with LPV/r resistance, they may exert pressure through the Gag and minor PR mutation pathways. Further investigations using site-directed mutagenesis are needed to determine the impact of Env mutations alone and in combination with Gag-PR mutations on viral fitness and LPV/r efficacy.

Systematic screening of viral and human genetic variation identifies antiretroviral resistance and immune escape link

Background:

Considering the remaining threat of drug-resistantmutations (DRMs) to antiretroviral treatment (ART) efficacy, we investigated how the selective pressure of human leukocyte antigen (HLA)-restricted cytotoxic T lymphocytes drives certain DRMs’ emergence and retention.

Methods:

We systematically screened DRM:HLA class I allele combinations in 3997 ART-naïve Swiss HIV Cohort Study (SHCS) patients. For each pair, a logistic regression model preliminarily tested for an association with the DRM as the outcome. The three HLA:DRM pairs remaining after multiple testing adjustment were analyzed in three ways: cross-sectional logistic regression models to determine any HLA/infection time interaction, survival analyses to examine if HLA type correlated with developing specific DRMs, and via NetMHCpan to find epitope binding evidence of immune escape.

Results:

Only one pair, RT-E138:HLA-B18, exhibited a significant interaction between infection duration and HLA. The survival analyses predicted two pairs with an increased hazard of developing DRMs: RT-E138:HLA-B18 and RT-V179:HLA-B35. RT-E138:HLA-B18 exhibited the greatest significance in both analyses (interaction term odds ratio [OR] 1.169 [95% confidence interval (CI) 1.075–1.273]; p-value<0.001; survival hazard ratio 12.211 [95% CI 3.523–42.318]; p-value<0.001). The same two pairs were also predicted by netMHCpan to have epitopic binding.

Conclusions:

We identified DRM:HLA pairs where HLA presence is associated with the presence or emergence of the DRM, indicating that the selective pressure for these mutations alternates direction depending on the presence of these HLA alleles.

Funding:

Funded by the Swiss National Science Foundation within the framework of the SHCS, and the University of Zurich, University Research Priority Program: Evolution in Action: From Genomes Ecosystems, in Switzerland.

Reviewing HIV-1 Gag Mutations in Protease Inhibitors Resistance: Insights for Possible Novel Gag Inhibitor Designs

HIV protease inhibitors against the viral protease are often hampered by drug resistance mutations in protease and in the viral substrate Gag. To overcome this drug resistance and inhibit viral maturation, targeting Gag alongside protease rather than targeting protease alone may be more efficient. In order to successfully inhibit Gag, understanding of its drug resistance mutations and the elicited structural changes on protease binding needs to be investigated. While mutations on Gag have already been mapped to protease inhibitor resistance, there remain many mutations, particularly the non-cleavage mutations, that are not characterized. Through structural studies to unravel how Gag mutations contributes to protease drug resistance synergistically, it is thus possible to glean insights to design novel Gag inhibitors. In this review, we discuss the structural role of both novel and previously reported Gag mutations in PI resistance, and how new Gag inhibitors can be designed.

Novel Protease Inhibitors Containing C-5-Modified bis-Tetrahydrofuranylurethane and Aminobenzothiazole as P2 and P2' Ligands That Exert Potent Antiviral Activity against Highly Multidrug-Resistant HIV-1 with a High Genetic Barrier against the Emergence of Drug Resistance

Combination antiretroviral therapy has achieved dramatic reductions in the mortality and morbidity in people with HIV-1 infection. Darunavir (DRV) represents a most efficacious and well-tolerated protease inhibitor (PI) with a high genetic barrier to the emergence of drug-resistant HIV-1. However, highly DRV-resistant variants have been reported in patients receiving long-term DRV-containing regimens. Here, we report three novel HIV-1 PIs (GRL-057-14, GRL-058-14, and GRL-059-14), all of which contain a P2-amino-substituted-bis-tetrahydrofuranylurethane (bis-THF) and a P2'-cyclopropyl-amino-benzothiazole (Cp-Abt). These PIs not only potently inhibit the replication of wild-type HIV-1 (50% effective concentration [EC₅₀], 0.22 nM to 10.4 nM) but also inhibit multi-PI-resistant HIV-1 variants, including highly DRV-resistant HIV_DRV ^R _P51 (EC₅₀, 1.6 nM to 30.7 nM). The emergence of HIV-1 variants resistant to the three compounds was much delayed in selection experiments compared to resistance to DRV, using a mixture of 11 highly multi-PI-resistant HIV-1 isolates as a starting HIV-1 population. GRL-057-14 showed the most potent anti-HIV-1 activity and greatest thermal stability with wild-type protease, and potently inhibited HIV-1 protease's proteolytic activity (K_i value, 0.10 nM) among the three PIs. Structural models indicate that the C-5-isopropylamino-bis-THF moiety of GRL-057-14 forms additional polar interactions with the active site of HIV-1 protease. Moreover, GRL-057-14's P1-bis-fluoro-methylbenzene forms strong hydrogen bonding and effective van der Waals interactions. The present data suggest that the combination of C-5-aminoalkyl-bis-THF, P1-bis-fluoro-methylbenzene, and P2'-Cp-Abt confers highly potent activity against wild-type and multi-PI-resistant HIV strains and warrant further development of the three PIs, in particular, that of GRL-057-14, as potential therapeutic for HIV-1 infection and AIDS.

The N-Terminus of the HIV-1 p6 Gag Protein Regulates Susceptibility to Degradation by IDE

As part of the Pr55^Gag polyprotein, p6 fulfills an essential role in the late steps of the replication cycle. However, almost nothing is known about the functions of the mature HIV-1 p6 protein. Recently, we showed that p6 is a bona fide substrate of the insulin-degrading enzyme (IDE), a ubiquitously expressed zinc metalloprotease. This phenomenon appears to be specific for HIV-1, since p6 homologs of HIV-2, SIV and EIAV were IDE-insensitive. Furthermore, abrogation of the IDE-mediated degradation of p6 reduces the replication capacity of HIV-1 in an Env-dependent manner. However, it remained unclear to which extent the IDE mediated degradation is phylogenetically conserved among HIV-1. Here, we describe two HIV-1 isolates with IDE resistant p6 proteins. Sequence comparison allowed deducing one single amino acid regulating IDE sensitivity of p6. Exchanging the N-terminal leucine residue of p6 derived from the IDE sensitive isolate HIV-1_NL4-3 with proline enhances its stability, while replacing Pro-1 of p6 from the IDE insensitive isolate SG3 with leucine restores susceptibility towards IDE. Phylogenetic analyses of this natural polymorphism revealed that the N-terminal leucine is characteristic for p6 derived from HIV-1 group M except for subtype A, which predominantly expresses p6 with an N-terminal proline. Consequently, p6 peptides derived from subtype A are not degraded by IDE. Thus, IDE mediated degradation of p6 is specific for HIV-1 group M isolates and not occasionally distributed among HIV-1.

Gp41 and Gag amino acids linked to HIV-1 protease inhibitor-based second-line failure in HIV-1 subtype A from Western Kenya

Introduction

Failure of protease‐inhibitor (PI)‐based second‐line antiretroviral therapy (ART) with medication adherence but no protease drug resistance mutations (DRMs) is not well understood. This study investigated the involvement of gp41 and gag as alternative mechanisms, not captured by conventional resistance testing and particularly relevant in resource‐limited settings where third‐line ART is limited.

Methods

We evaluated gp41 and gag for unique amino acids in seven subtype A infected Kenyans failing second‐line therapy with no PI resistance yet detectable lopinavir (query dataset), compared to seven similar‐setting patients with PI resistance or undetectable lopinavir and 69 publically available subtype A Kenyan whole‐genomes sequences.

Results

Three gp41 (607T, 641L, 721I) and four gag (124S, 143V, 339P, 357S) amino acids were significantly more frequent in the query dataset compared to the other datasets, with significantly high co‐occurrence.

Conclusion

The genotypic analysis of a unique group of HIV‐1 subtype A infected patients, identified seven amino acids that could potentially contribute to a multi‐gene mechanism of PI‐based ART failure in the absence of PI DR mutations.

Evolution of gag and gp41 in Patients Receiving Ritonavir-Boosted Protease Inhibitors

Several groups have proposed that genotypic determinants in gag and the gp41 cytoplasmic domain (gp41-CD) reduce protease inhibitor (PI) susceptibility without PI-resistance mutations in protease. However, no gag and gp41-CD mutations definitively responsible for reduced PI susceptibility have been identified in individuals with virological failure (VF) while receiving a boosted PI (PI/r)-containing regimen. To identify gag and gp41 mutations under selective PI pressure, we sequenced gag and/or gp41 in 61 individuals with VF on a PI/r (n = 40) or NNRTI (n = 20) containing regimen. We quantified nonsynonymous and synonymous changes in both genes and identified sites exhibiting signal for directional or diversifying selection. We also used published gag and gp41 polymorphism data to highlight mutations displaying a high selection index, defined as changing from a conserved to an uncommon amino acid. Many amino acid mutations developed in gag and in gp41-CD in both the PI- and NNRTI-treated groups. However, in neither gene, were there discernable differences between the two groups in overall numbers of mutations, mutations displaying evidence of diversifying or directional selection, or mutations with a high selection index. If gag and/or gp41 encode PI-resistance mutations, they may not be confined to consistent mutations at a few sites.