Genetic Changes in HIV-1 Gag-Protease Associated with Protease Inhibitor-Based Therapy Failure in Pediatric Patients

Understanding the co-evolutionary molecular mechanisms of resistance in the HIV-1 Gag and protease

Due to high human immunodeficiency virus type 1 (HIV-1) subtype C infections coupled with increasing antiretroviral treatment failure, the elucidation of complex drug resistance mutational patterns arising through protein co-evolution is required. Despite the inclusion of potent protease inhibitors Lopinavir (LPV) and Darunavir (DRV) in second- and third-line therapies, many patients still fail treatment due to the accumulation of mutations in protease (PR) and recently, Gag. To understand the co-evolutionary molecular mechanisms of resistance in the HIV-1 PR and Gag, we performed 100 ns molecular dynamic simulations on multidrug resistant PR's when bound to LPV, DRV or a mutated A431V NC|p1 Gag cleavage site (CS). Here we showed that distinct changes in PR's active site, flap and elbow regions due to several PR resistance mutations (L10F, M46I, I54V, L76V, V82A) were found to alter LPV and DRV drug binding. However, binding was significantly exacerbated when the mutant PRs were bound to the NC|p1 Gag CS. Although A431V was shown to coordinate several residues in PR, the L76V PR mutation was found to have a significant role in substrate recognition. Consequently, a greater binding affinity was observed when the mutated substrate was bound to an L76V-inclusive PR mutant (G_bind: -62.46 ± 5.75 kcal/mol) than without (Gbind: -50.34 ± 6.28 kcal/mol). These data showed that the co-selection of resistance mutations in the enzyme and substrate can simultaneously constrict regions in PR's active site whilst flexing the flaps to allow flexible movement of the substrate and multiple, complex mechanisms of resistance to occur. Communicated by Ramaswamy H. Sarma.

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.

Protease and gag diversity and drug resistance mutations among treatment-naive Mexican people living with HIV

Introduction

In Mexico, HIV genotyping is performed in people living with HIV (PLWH) failing their first-line antiretroviral (ARV) regimen; it is not routinely done for all treatment-naive PLWH before ARV initiation. The first nationally representative survey published in 2016 reported that the prevalence of pretreatment drug mutations in treatment-naive Mexican PLWH was 15.5% to any antiretroviral drug and 10.6% to non-nucleoside reverse transcriptase inhibitors (NNRTIs) using conventional Sanger sequencing. Most reports in Mexico focus on HIV pol gene and nucleoside and non-nucleoside reverse transcriptase inhibitor (NRTI and NNRTI) drug resistance mutations (DRMs) prevalence, using Sanger sequencing, next-generation sequencing (NGS) or both. To our knowledge, NGS has not be used to detect pretreatment drug resistance mutations (DRMs) in the HIV protease (PR) gene and its substrate the Gag polyprotein.

Methods

Treatment-naive adult Mexican PLWH were recruited between 2016 and 2019. HIV Gag and protease sequences were obtained by NGS and DRMs were identified using the WHO surveillance drug resistance mutation (SDRM) list.

Results

One hundred PLWH attending a public national reference hospital were included. The median age was 28 years-old, and most were male. The median HIV viral load was 4.99 [4.39–5.40] log copies/mL and median CD4 cell count was 150 [68.0–355.78] cells/mm³. As expected, most sequences clustered with HIV-1 subtype B (97.9%). Major PI resistance mutations were detected: 8 (8.3%) of 96 patients at a detection threshold of 1% and 3 (3.1%) at a detection threshold of 20%. A total of 1184 mutations in Gag were detected, of which 51 have been associated with resistance to PI, most of them were detected at a threshold of 20%. Follow-up clinical data was available for 79 PLWH at 6 months post-ART initiation, seven PLWH failed their first ART regimen; however no major PI mutations were identified in these individuals at baseline.

Conclusions

The frequency of DRM in the HIV protease was 7.3% at a detection threshold of 1% and 3.1% at a detection threshold of 20%. NGS-based HIV drug resistance genotyping provide improved detection of DRMs. Viral load was used to monitor ARV response and treatment failure was 8.9%.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07446-8.

HIV-1 drug resistance in adults and adolescents on protease inhibitor-based antiretroviral treatment in KwaZulu-Natal Province, South Africa

BACKGROUND

In low- and middle-income countries, increasing levels of HIV drug resistance (HIVDR) on second-line protease inhibitor (PI)-based regimens are a cause for concern, given limited drug options for third-line antiretroviral therapy (ART).

OBJECTIVES

We conducted a retrospective analysis of routine HIV-1 genotyping laboratory data from KwaZulu-Natal, in South Africa, to describe the frequency and patterns of HIVDR mutations and their consequent impact on standardized third-line regimens.

METHODS

This was a cross-sectional analysis of all HIV-1 genotypic resistance tests conducted by the National Health Laboratory Service in KwaZulu-Natal, South Africa (Jan 2015 - Dec 2016), for adults and adolescents (age ≥10 years) on second-line PI-based ART with virological failure. We assigned a third-line regimen to each record, based on a national treatment algorithm and calculated the genotypic susceptibility score (GSS) for that regimen.

RESULTS

Of 348 samples analyzed, 287 (83%) had at least one drug resistance mutation (DRM) and 114 (33%) had at least one major PI DRM. Major PI resistance was associated with longer duration on second-line ART (aOR per 6-months, 1.11, 95% CI 1.04-1.19) and older age (aOR 1.03, 95% CI 1.01-1.05). Of 112 patients requiring third-line ART, 12 (11%) had a GSS of <2 for the algorithm-assigned third-line regimen.

CONCLUSIONS

One in three people failing second-line ART had significant PI DRMs. A subgroup of these individuals had extensive HIVDR, where the predicted activity of third-line ART was suboptimal, highlighting the need for continuous evaluation of outcomes on third-line regimens and close monitoring for emergent HIV-1 integrase-inhibitor resistance.

Baseline PI susceptibility by HIV-1 Gag-protease phenotyping and subsequent virological suppression with PI-based second-line ART in Nigeria

Objectives

Previous work showed that gag-protease-derived phenotypic susceptibility to PIs differed between HIV-1 subtype CRF02_AG/subtype G-infected patients who went on to successfully suppress viral replication versus those who experienced virological failure of lopinavir/ritonavir monotherapy as first-line treatment in a clinical trial. We analysed the relationship between PI susceptibility and outcome of second-line ART in Nigeria, where subtypes CRF02_AG/G dominate the epidemic.

Methods

Individuals who experienced second-line failure with ritonavir-boosted PI-based ART were matched (by subtype, sex, age, viral load, duration of treatment and baseline CD4 count) to those who achieved virological response (‘successes’). Successes were defined by viral load <400 copies of HIV-1 RNA/mL by week 48. Full-length Gag-protease was amplified from patient samples for in vitro phenotypic susceptibility testing, with PI susceptibility expressed as IC₅₀ fold change (FC) relative to a subtype B reference strain.

Results

The median (IQR) lopinavir IC₅₀ FC was 4.04 (2.49–7.89) for virological failures and 4.13 (3.14–8.17) for virological successes (P = 0.94). One patient had an FC >10 for lopinavir at baseline and experienced subsequent virological failure with ritonavir-boosted lopinavir as the PI. There was no statistically significant difference in single-round replication efficiency between the two groups (P = 0.93). There was a moderate correlation between single-round replication efficiency and FC for lopinavir (correlation coefficient 0.32).

Conclusions

We found no impact of baseline HIV-1 Gag-protease-derived phenotypic susceptibility on outcomes of PI-based second-line ART in Nigeria.

Drug susceptibility and replication capacity of a rare HIV-1 subtype C protease hinge region variant

BACKGROUND

Protease inhibitors form the main component of second-line antiretroviral treatment in South Africa. Despite their efficacy, mutations arising within the HIV-1 gag and protease coding regions contribute to the development of resistance against this class of drug. In this paper we investigate a South African HIV-1 subtype C Gag-protease that contains a hinge region mutation and insertion (N37T↑V).

METHODS

In vitro single-cycle drug susceptibility and viral replication capacity assays were performed on W1201i, a wild-type reference isolate (MJ4) and a chimeric construct (MJ4GagN37T↑VPR). Additionally, enzyme assays were performed on the N37T↑V protease and a wild-type reference protease.

RESULTS

W1201i showed a small (threefold), but significant (P<0.0001) reduction in drug susceptibility to darunavir compared with MJ4. Substitution of W1201i-Gag with MJ4-Gag resulted in an additional small (twofold), but significant (P<0.01) reduction in susceptibility to lopinavir and atazanavir. The W1201i pseudovirus had a significantly (P<0.01) reduced replication capacity (16.4%) compared with the MJ4. However, this was dramatically increased to 164% (P<0.05) when W1201i-Gag was substituted with MJ4-Gag. Furthermore, the N37T↑V protease displayed reduced catalytic processing compared with the SK154 protease.

CONCLUSIONS

Collectively, these data suggest that the N37T↑V mutation and insertion increases viral infectivity and decreases drug susceptibility. These variations are classified as secondary mutations, and indirectly impact inhibitor binding, enzyme fitness and enzyme stability. Additionally, polymorphisms arising in Gag can modify the impact of protease with regards to viral replication and susceptibility to protease inhibitors.

Gag P2/NC and pol genetic diversity, polymorphism, and drug resistance mutations in HIV-1 CRF02_AG- and non-CRF02_AG-infected patients in Yaoundé, Cameroon

In HIV-1 subtype-B, specific mutations in Gag cleavage sites (CS) are associated with treatment failure, with limited knowledge among non-B subtypes. We analyzed non-B HIV-1 gag and pol (protease/reverse-transcriptase) sequences from Cameroonians for drug resistance mutations (DRMs) in the gag P2/NC CS, and pol major DRMs. Phylogeny of the 141 sequences revealed a high genetic diversity (12 subtypes): 67.37% CRF02_AG versus 32.6% non-CRF02_AG. Overall, 7.3% transmitted and 34.3% acquired DRMs were found, including M184V, thymidine analogue mutations (T215F, D67N, K70R, K219Q), NNRTIs (L100I, Y181C, K103N, V108I, Y188L), and PIs (V82L). Twelve subjects [10 with HIV-1 CRF02_AG, 8 treatment-naïve and 4 on 3TC-AZT-NVP] showed 3 to 4 mutations in the Gag P2/NC CS: S373Q/T/A, A374T/S/G/N, T375S/A/N/G, I376V, G381S, and R380K. Subjects with or without Gag P2/NC CS mutations showed no significant difference in viral loads. Treatment-naïve subjects harboring NRTI-DRMs had significantly lower CD4 cells than those with NRTI-DRMs on ART (p = 0.042). Interestingly, two subjects had major DRMs to NRTIs, NNRTIs, and 4 mutations in the Gag P2/NC CS. In this prevailing CRF02_AG population with little exposure to PIs (~3%), mutations in the Gag P2/NC CS could increase the risk of treatment failure if there is increased use of PIs-based therapy.

Virological Outcomes of Second-line Protease Inhibitor-Based Treatment for Human Immunodeficiency Virus Type 1 in a High-Prevalence Rural South African Setting: A Competing-Risks Prospective Cohort Analysis

Background

Second-line antiretroviral therapy (ART) based on ritonavir-boosted protease inhibitors (bPIs) represents the only available option after first-line failure for the majority of individuals living with human immunodeficiency virus (HIV) worldwide. Maximizing their effectiveness is imperative.

Methods

This cohort study was nested within the French National Agency for AIDS and Viral Hepatitis Research (ANRS) 12249 Treatment as Prevention (TasP) cluster-randomized trial in rural KwaZulu-Natal, South Africa. We prospectively investigated risk factors for virological failure (VF) of bPI-based ART in the combined study arms. VF was defined by a plasma viral load >1000 copies/mL ≥6 months after initiating bPI-based ART. Cumulative incidence of VF was estimated and competing risk regression was used to derive the subdistribution hazard ratio (SHR) of the associations between VF and patient clinical and demographic factors, taking into account death and loss to follow-up.

Results

One hundred one participants contributed 178.7 person-years of follow-up. Sixty-five percent were female; the median age was 37.4 years. Second-line ART regimens were based on ritonavir-boosted lopinavir, combined with zidovudine or tenofovir plus lamivudine or emtricitabine. The incidence of VF on second-line ART was 12.9 per 100 person-years (n = 23), and prevalence of VF at censoring was 17.8%. Thirteen of these 23 (56.5%) virologic failures resuppressed after a median of 8.0 months (interquartile range, 2.8-16.8 months) in this setting where viral load monitoring was available. Tuberculosis treatment was associated with VF (SHR, 11.50 [95% confidence interval, 3.92-33.74]; P < .001).

Conclusions

Second-line VF was frequent in this setting. Resuppression occurred in more than half of failures, highlighting the value of viral load monitoring of second-line ART. Tuberculosis was associated with VF; therefore, novel approaches to optimize the effectiveness of PI-based ART in high-tuberculosis-burden settings are needed.

Clinical Trials Registration

NCT01509508.

Wide variation in susceptibility of transmitted/founder HIV-1 subtype C Isolates to protease inhibitors and association with in vitro replication efficiency

The gag gene is highly polymorphic across HIV-1 subtypes and contributes to susceptibility to protease inhibitors (PI), a critical class of antiretrovirals that will be used in up to 2 million individuals as second-line therapy in sub Saharan Africa by 2020. Given subtype C represents around half of all HIV-1 infections globally, we examined PI susceptibility in subtype C viruses from treatment-naïve individuals. PI susceptibility was measured in a single round infection assay of full-length, replication competent MJ4/gag chimeric viruses, encoding the gag gene and 142 nucleotides of pro derived from viruses in 20 patients in the Zambia-Emory HIV Research Project acute infection cohort. Ten-fold variation in susceptibility to PIs atazanavir and lopinavir was observed across 20 viruses, with EC50s ranging 0.71-6.95 nM for atazanvir and 0.64-8.54 nM for lopinavir. Ten amino acid residues in Gag correlated with lopinavir EC50 (p < 0.01), of which 380 K and 389I showed modest impacts on in vitro drug susceptibility. Finally a significant relationship between drug susceptibility and replication capacity was observed for atazanavir and lopinavir but not darunavir. Our findings demonstrate large variation in susceptibility of PI-naïve subtype C viruses that appears to correlate with replication efficiency and could impact clinical outcomes.

Contribution of Gag and Protease to HIV-1 Phenotypic Drug Resistance in Pediatric Patients Failing Protease Inhibitor-Based Therapy

Protease inhibitors (PIs) are used as a first-line regimen in HIV-1-infected children. Here we investigated the phenotypic consequences of amino acid changes in Gag and protease on lopinavir (LPV) and ritonavir (RTV) susceptibility among pediatric patients failing PI therapy. The Gag-protease from isolates from 20 HIV-1 subtype C-infected pediatric patients failing an LPV and/or RTV-based regimen was phenotyped using a nonreplicativein vitroassay. Changes in sensitivity to LPV and RTV relative to that of the matched baseline (pretherapy) sample were calculated. Gag and protease amino acid substitutions associated with PI failure were created in a reference clone by site-directed mutagenesis and assessed. Predicted phenotypes were determined using the Stanford drug resistance algorithm. Phenotypic resistance or reduced susceptibility to RTV and/or LPV was observed in isolates from 10 (50%) patients, all of whom had been treated with RTV. In most cases, this was associated with protease resistance mutations, but substitutions at Gag cleavage and noncleavage sites were also detected. Gag amino acid substitutions were also found in isolates from three patients with reduced drug susceptibilities who had wild-type protease. Site-directed mutagenesis confirmed that some amino acid changes in Gag contributed to PI resistance but only in the presence of major protease resistance-associated substitutions. The isolates from all patients who received LPV exclusively were phenotypically susceptible. Baseline isolates from the 20 patients showed a large (47-fold) range in the 50% effective concentration of LPV, which accounted for most of the discordance seen between the experimentally determined and the predicted phenotypes. Overall, the inclusion of thegaggene and the use of matched baseline samples provided a more comprehensive assessment of the effect of PI-induced amino acid changes on PI resistance. The lack of phenotypic resistance to LPV supports the continued use of this drug in pediatric patients.

Evidence for Reduced Drug Susceptibility without Emergence of Major Protease Mutations following Protease Inhibitor Monotherapy Failure in the SARA Trial

Background

Major protease mutations are rarely observed following failure with protease inhibitors (PI), and other viral determinants of failure to PI are poorly understood. We therefore characterized Gag-Protease phenotypic susceptibility in subtype A and D viruses circulating in East Africa following viral rebound on PIs.

Methods

Samples from baseline and treatment failure in patients enrolled in the second line LPV/r trial SARA underwent phenotypic susceptibility testing. Data were expressed as fold-change in susceptibility relative to a LPV-susceptible reference strain.

Results

We cloned 48 Gag-Protease containing sequences from seven individuals and performed drug resistance phenotyping from pre-PI and treatment failure timepoints in seven patients. For the six patients where major protease inhibitor resistance mutations did not emerge, mean fold-change EC₅₀ to LPV was 4.07 fold (95% CI, 2.08–6.07) at the pre-PI timepoint. Following viral failure the mean fold-change in EC₅₀ to LPV was 4.25 fold (95% CI, 1.39–7.11, p = 0.91). All viruses remained susceptible to DRV. In our assay system, the major PI resistance mutation I84V, which emerged in one individual, conferred a 10.5-fold reduction in LPV susceptibility. One of the six patients exhibited a significant reduction in susceptibility between pre-PI and failure timepoints (from 4.7 fold to 9.6 fold) in the absence of known major mutations in protease, but associated with changes in Gag: V7I, G49D, R69Q, A120D, Q127K, N375S and I462S. Phylogenetic analysis provided evidence of the emergence of genetically distinct viruses at the time of treatment failure, indicating ongoing viral evolution in Gag-protease under PI pressure.

Conclusions

Here we observe in one patient the development of significantly reduced susceptibility conferred by changes in Gag which may have contributed to treatment failure on a protease inhibitor containing regimen. Further phenotype-genotype studies are required to elucidate genetic determinants of protease inhibitor failure in those who fail without traditional resistance mutations whilst PI use is being scaled up globally.

Mutations Proximal to Sites of Autoproteolysis and the α-Helix That Co-evolve under Drug Pressure Modulate the Autoprocessing and Vitality of HIV-1 Protease

N-Terminal self-cleavage (autoprocessing) of the HIV-1 protease precursor is crucial for liberating the active dimer. Under drug pressure, evolving mutations are predicted to modulate autoprocessing, and the reduced catalytic activity of the mature protease (PR) is likely compensated by enhanced conformational/dimer stability and reduced susceptibility to self-degradation (autoproteolysis). One such highly evolved, multidrug resistant protease, PR20, bears 19 mutations contiguous to sites of autoproteolysis in retroviral proteases, namely clusters 1-3 comprising residues 30-37, 60-67, and 88-95, respectively, accounting for 11 of the 19 mutations. By systematically replacing corresponding clusters in PR with those of PR20, and vice versa, we assess their influence on the properties mentioned above and observe no strict correlation. A 10-35-fold decrease in the cleavage efficiency of peptide substrates by PR20, relative to PR, is reflected by an only ∼4-fold decrease in the rate of Gag processing with no change in cleavage order. Importantly, optimal N-terminal autoprocessing requires all 19 PR20 mutations as evaluated in vitro using the model precursor TFR-PR20 in which PR is flanked by the transframe region. Substituting PR20 cluster 3 into TFR-PR (TFR-PR(PR20-3)) requires the presence of PR20 cluster 1 and/or 2 for autoprocessing. In accordance, substituting PR clusters 1 and 2 into TFR-PR20 affects the rate of autoprocessing more drastically (>300-fold) compared to that of TFR-PR(PR20-3) because of the cumulative effect of eight noncluster mutations present in TFR-PR20(PR-12). Overall, these studies imply that drug resistance involves a complex synchronized selection of mutations modulating all of the properties mentioned above governing PR regulation and function.