Gag-protease coevolution analyses define novel structural surfaces in the HIV-1 matrix and capsid involved in resistance to Protease Inhibitors

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.

Selection among site-dependent structurally constrained substitution models of protein evolution by approximate Bayesian computation

MOTIVATION

The selection among substitution models of molecular evolution is fundamental for obtaining accurate phylogenetic inferences. At the protein level, evolutionary analyses are traditionally based on empirical substitution models but these models make unrealistic assumptions and are being surpassed by structurally constrained substitution (SCS) models. The SCS models often consider site-dependent evolution, a process that provides realism but complicates their implementation into likelihood functions that are commonly used for substitution model selection.

RESULTS

We present a method to perform selection among site-dependent SCS models, also among empirical and site-dependent SCS models, based on the approximate Bayesian computation (ABC) approach and its implementation into the computational framework ProteinModelerABC. The framework implements ABC with and without regression adjustments and includes diverse empirical and site-dependent SCS models of protein evolution. Using extensive simulated data, we found that it provides selection among SCS and empirical models with acceptable accuracy. As illustrative examples, we applied the framework to analyze a variety of protein families observing that SCS models fit them better than the corresponding best-fitting empirical substitution models.

AVAILABILITY AND IMPLEMENTATION

ProteinModelerABC is freely available from https://github.com/DavidFerreiro/ProteinModelerABC, can run in parallel and includes a graphical user interface. The framework is distributed with detailed documentation and ready-to-use examples.

Antiretroviral Imprints and Genomic Plasticity of HIV-1 pol in Non-clade B: Implications for Treatment

Combinational antiretroviral therapy (cART) is the most effective tool to prevent and control HIV-1 infection without an effective vaccine. However, HIV-1 drug resistance mutations (DRMs) and naturally occurring polymorphisms (NOPs) can abrogate cART efficacy. Here, we aimed to characterize the HIV-1 pol mutation landscape in Cameroon, where highly diverse HIV clades circulate, and identify novel treatment-associated mutations that can potentially affect cART efficacy. More than 8,000 functional Cameroonian HIV-1 pol sequences from 1987 to 2020 were studied for DRMs and NOPs. Site-specific amino acid frequencies and quaternary structural features were determined and compared between periods before (≤2003) and after (2004-2020) regional implementation of cART. cART usage in Cameroon induced deep mutation imprints in reverse transcriptase (RT) and to a lower extent in protease (PR) and integrase (IN), according to their relative usage. In the predominant circulating recombinant form (CRF) 02_AG (CRF02_AG), 27 canonical DRMs and 29 NOPs significantly increased or decreased in RT during cART scale-up, whereas in IN, no DRM and only seven NOPs significantly changed. The profound genomic imprints and higher prevalence of DRMs in RT compared to PR and IN mirror the dominant use of reverse transcriptase inhibitors (RTIs) in sub-Saharan Africa and the predominantly integrase strand transfer inhibitor (InSTI)-naïve study population. Our results support the potential of InSTIs for antiretroviral treatment in Cameroon; however, close surveillance of IN mutations will be required to identify emerging resistance patterns, as observed in RT and PR. Population-wide genomic analyses help reveal the presence of selective pressures and viral adaptation processes to guide strategies to bypass resistance and reinstate effective treatment.

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.

Viral and Cellular factors leading to the Loss of CD4 Homeostasis in HIV-1 Viremic Nonprogressors

Human immunodeficiency virus type 1 (HIV-1) viremic nonprogressors (VNPs) represent a very rare HIV-1 extreme phenotype. VNPs are characterized by persistent high plasma viremia and maintenance of CD4⁺ T-cell counts in the absence of treatment. However, the causes of nonpathogenic HIV-1 infection in VNPs remain elusive. Here, we identified for the first time two VNPs who experienced the loss of CD4⁺ homeostasis (LoH) after more than 13 years. We characterized in deep detail viral and host factors associated with the LoH and compared with standard VNPs and healthy controls. The viral factors determined included HIV-1 coreceptor usage and replicative capacity. Changes in CD4⁺ and CD8⁺ T-cell activation, maturational phenotype, and expression of CCR5 and CXCR6 in CD4⁺ T-cells were also evaluated as host-related factors. Consistently, we determined a switch in HIV-1 coreceptor use to CXCR4 concomitant with an increase in replicative capacity at the LoH for the two VNPs. Moreover, we delineated an increase in the frequency of HLA-DR+CD38⁺ CD4⁺ and CD8⁺ T cells and traced the augment of naive T-cells upon polyclonal activation with LoH. Remarkably, very low and stable levels of CCR5 and CXCR6 expression in CD4⁺ T-cells were measured over time. Overall, our results demonstrated HIV-1 evolution toward highly pathogenic CXCR4 strains in the context of very limited and stable expression of CCR5 and CXCR6 in CD4⁺ T cells as potential drivers of LoH in VNPs. These data bring novel insights into the correlates of nonpathogenic HIV-1 infection.

IMPORTANCE The mechanism behind nonpathogenic human immunodeficiency virus type 1 (HIV-1) infection remains poorly understood, mainly because of the very low frequency of viremic nonprogressors (VNPs). Here, we report two cases of VNPs who experienced the loss of CD4⁺ T-cell homeostasis (LoH) after more than 13 years of HIV-1 infection. The deep characterization of viral and host factors supports the contribution of viral and host factors to the LoH in VNPs. Thus, HIV-1 evolution toward highly replicative CXCR4 strains together with changes in T-cell activation and maturational phenotypes were found. Moreover, we measured very low and stable levels of CCR5 and CXCR6 in CD4⁺ T-cells over time. These findings support viral evolution toward X4 strains limited by coreceptor expression to control HIV-1 pathogenesis and demonstrate the potential of host-dependent factors, yet to be fully elucidated in VNPs, to control HIV-1 pathogenesis.

Elasticity-Associated Functionality and Inhibition of the HIV Protease

HIV protease plays a critical role in the life cycle of the virus through the generation of mature and infectious virions. Detailed knowledge of the structure of the enzyme and its substrate has led to the development of protease inhibitors. However, the development of resistance to all currently available protease inhibitors has contributed greatly to the decreased success of antiretroviral therapy. When therapy failure occurs, multiple mutations are found within the protease sequence starting with primary mutations, which directly impact inhibitor binding, which can also negatively impact viral fitness and replicative capacity by decreasing the binding affinity of the natural substrates to the protease. As such, secondary mutations which are located outside of the active site region accumulate to compensate for the recurrently deleterious effects of primary mutations. However, the resistance mechanism of these secondary mutations is not well understood, but what is known is that these secondary mutations contribute to resistance in one of two ways, either through increasing the energetic penalty associated with bringing the protease into the closed conformation, or, through decreasing the stability of the protein/drug complex in a manner that increases the dissociation rate of the drug, leading to diminished inhibition. As a result, the elasticity of the enzyme-substrate complex has been implicated in the successful recognition and catalysis of the substrates which may be inferred to suggest that the elasticity of the enzyme/drug complex plays a role in resistance. A realistic representation of the dynamic nature of the protease may provide a more powerful tool in structure-based drug design algorithms.

The impact of Gag non-cleavage site mutations on HIV-1 viral fitness from integrative modelling and simulations

The high mutation rate in retroviruses is one of the leading causes of drug resistance. In human immunodeficiency virus type-1 (HIV-1), synergistic mutations in its protease and the protease substrate - the Group-specific antigen (Gag) polyprotein - work together to confer drug resistance against protease inhibitors and compensate the mutations affecting viral fitness. Some Gag mutations can restore Gag-protease binding, yet most Gag-protease correlated mutations occur outside of the Gag cleavage site. To investigate the molecular basis for this, we now report multiscale modelling approaches to investigate various sequentially cleaved Gag products in the context of clinically relevant mutations that occur outside of the cleavage sites, including simulations of the largest Gag proteolytic product in its viral membrane-bound state. We found that some mutations, such as G123E and H219Q, involve direct interaction with cleavage site residues to influence their local environment, while certain mutations in the matrix domain lead to the enrichment of lipids important for Gag targeting and assembly. Collectively, our results reveal why non-cleavage site mutations have far-reaching implications outside of Gag proteolysis, with important consequences for drugging Gag maturation intermediates and tackling protease inhibitor resistance.

Nucleocapsid Protein Precursors NCp9 and NCp15 Suppress ATP-Mediated Rescue of AZT-Terminated Primers by HIV-1 Reverse Transcriptase

In HIV-1, development of resistance to AZT (3'-azido-3'-deoxythymidine) is mediated by the acquisition of thymidine analogue resistance mutations (TAMs) (i.e., M41L, D67N, K70R, L210W, T215F/Y, and K219E/Q) in the viral reverse transcriptase (RT). Clinically relevant combinations of TAMs, such as M41L/T215Y or D67N/K70R/T215F/K219Q, enhance the ATP-mediated excision of AZT monophosphate (AZTMP) from the 3' end of the primer, allowing DNA synthesis to continue. Additionally, during HIV-1 maturation, the Gag polyprotein is cleaved to release a mature nucleocapsid protein (NCp7) and two intermediate precursors (NCp9 and NCp15). NC proteins interact with the viral genome and facilitate the reverse transcription process. Using wild-type and TAM-containing RTs, we showed that both NCp9 and NCp15 inhibited ATP-mediated rescue of AZTMP-terminated primers annealed to RNA templates but not DNA templates, while NCp7 had no effect on rescue activity. RNase H inactivation by introducing the active-site mutation E478Q led to the loss of the inhibitory effect shown by NCp9. NCp15 had a stimulatory effect on the RT's RNase H activity not observed with NCp7 and NCp9. However, analysis of RNase H cleavage patterns revealed that in the presence of NCp9, RNA/DNA complexes containing duplexes of 12 bp had reduced stability in comparison with those obtained in the absence of NC or with NCp7 or NCp15. These effects are expected to have a strong influence on the inhibitory action of NCp9 and NCp15 by affecting the efficiency of RNA-dependent DNA polymerization after unblocking DNA primers terminated with AZTMP and other nucleotide analogues.

HIV-1 Gag mutations alone are sufficient to reduce darunavir susceptibility during virological failure to boosted PI therapy

BACKGROUND

Virological failure (VF) to boosted PIs with a high genetic barrier is not usually linked to the development of resistance-associated mutations in the protease gene.

METHODS

From a cohort of 520 HIV-infected subjects treated with lopinavir/ritonavir or darunavir/ritonavir monotherapy, we retrospectively identified nine patients with VF. We sequenced the HIV-1 Gag-protease region and generated clonal virus from plasma samples. We characterized phenotypically clonal variants in terms of replicative capacity and susceptibility to PIs. Also, we used VESPA to identify signature mutations and 3D molecular modelling information to detect conformational changes in the Gag region.

RESULTS

All subjects analysed harboured Gag-associated polymorphisms in the absence of resistance mutations in the protease gene. Most Gag changes occurred outside Gag cleavage sites. VESPA analyses identified K95R and R286K (P < 0.01) as signature mutations in Gag present at VF. In one out of four patients with clonal analysis available, we identified clonal variants with high replicative capacity and 8- to 13-fold reduction in darunavir susceptibility. These clonal variants harboured K95R, R286K and additional mutations in Gag. Low susceptibility to darunavir was dependent on the Gag sequence context. All other clonal variants analysed preserved drug susceptibility and virus replicative capacity.

CONCLUSIONS

Gag mutations may reduce darunavir susceptibility in the absence of protease mutations while preserving viral fitness. This effect is Gag-sequence context dependent and may occur during boosted PI failure.

Gag-protease coevolution shapes the outcome of lopinavir-inclusive treatment regimens in chronically infected HIV-1 subtype C patients

MOTIVATION

Commonly, protease inhibitor failure is characterized by the development of multiple protease resistance mutations (PRMs). While the impact of PRMs on therapy failure are understood, the introduction of Gag mutations with protease remains largely unclear.

RESULTS

Here, we utilized phylogenetic analyses and Bayesian network learning as tools to understand Gag-protease coevolution and elucidate the pathways leading to Lopinavir failure in HIV-1 subtype C infected patients. Our analyses indicate that while PRMs coevolve in response to drug selection pressure within protease, the Gag mutations added to the existing network while specifically interacting with known Lopinavir failure PRMs. Additionally, the selection of mutations at specific positions in Gag-protease suggests that these coevolving mutational changes occurs to maintain structural integrity during Gag cleavage.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Shall I trust the report? Variable performance of Sanger sequencing revealed by deep sequencing on HIV drug resistance mutation detection

BACKGROUND

The clinical utilisation of deep sequencing in HIV treatment has been hindered due to its unknown correlation with standard Sanger genotyping and the undetermined value of minority drug resistance mutation (DRM) detection.

OBJECTIVES

To compare deep sequencing performance to standard Sanger genotyping with clinical samples, in an effort to delineate the correlation between the results from the two methods and to find the optimal deep sequencing threshold for clinical utilisation.

METHODS

We conducted a retrospective study using stored plasma collected from August 2014 to March 2018 for HIV genotyping with the commercial Sanger genotyping kit. Samples with available Sanger genotyping reports were further deep sequenced. Drug resistance was interpreted according to the Stanford HIV drug resistance database algorithm.

RESULTS

At 15-25% minority detection thresholds, 9-15% cases had underestimated DRMs by Sanger sequencing. The concordance between the Sanger and deep sequencing reports was 68-82% in protease-reverse transcriptase region and 88-97% in integrase region at 5-25% thresholds. The undetected drug resistant minority variants by Sanger sequencing contributed to the lower negative predictive value of Sanger genotyping in cases harbouring DRMs.

CONCLUSIONS

Use of deep sequencing improved detection of antiretroviral resistance mutations especially in cases with virological failure or previous treatment interruption. Deep sequencing with 10-15% detection thresholds may be considered a suitable substitute for Sanger sequencing on antiretroviral DRM detection.

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.

Selection analyses of paired HIV-1 gag and gp41 sequences obtained before and after antiretroviral therapy

Most HIV-1-infected individuals with virological failure on a pharmacologically-boosted protease inhibitor (PI) regimen do not develop PI-resistance protease mutations. One proposed explanation is that HIV-1 gag or gp41 cytoplasmic domain mutations might also reduce PI susceptibility. In a recent study of paired gag and gp41 sequences from individuals with virological failure on a PI regimen, we did not identify PI-selected mutations and concluded that if such mutations existed, larger numbers of paired sequences from multiple studies would be needed for their identification. In this study, we generated site-specific amino acid profiles using gag and gp41 published sequences from 5,338 and 4,242 ART-naïve individuals, respectively, to assist researchers identify unusual mutations arising during therapy and to provide scripts for performing established and novel maximal likelihood estimates of dN/dS substitution rates in paired sequences. The pipelines used to generate the curated sequences, amino acid profiles, and dN/dS analyses will facilitate the application of consistent methods to paired gag and gp41 sequence datasets and expedite the identification of potential sites under PI-selection pressure.

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.