The HIV-1 late domain-2 S40A polymorphism in antiretroviral (or ART)-exposed individuals influences protease inhibitor susceptibility
© The Author(s) 2016
Received: 5 February 2016
Accepted: 21 August 2016
Published: 6 September 2016
The p6 region of the HIV-1 structural precursor polyprotein, Gag, contains two motifs, P7TAP11 and L35YPLXSL41, designated as late (L) domain-1 and -2, respectively. These motifs bind the ESCRT-I factor Tsg101 and the ESCRT adaptor Alix, respectively, and are critical for efficient budding of virus particles from the plasma membrane. L domain-2 is thought to be functionally redundant to PTAP. To identify possible other functions of L domain-2, we examined this motif in dominant viruses that emerged in a group of 14 women who had detectable levels of HIV-1 in both plasma and genital tract despite a history of current or previous antiretroviral therapy.
Remarkably, variants possessing mutations or rare polymorphisms in the highly conserved L domain-2 were identified in seven of these women. A mutation in a conserved residue (S40A) that does not reduce Gag interaction with Alix and therefore did not reduce budding efficiency was further investigated. This mutation causes a simultaneous change in the Pol reading frame but exhibits little deficiency in Gag processing and virion maturation. Whether introduced into the HIV-1 NL4-3 strain genome or a model protease (PR) precursor, S40A reduced production of mature PR. This same mutation also led to high level detection of two extended forms of PR that were fairly stable compared to the WT in the presence of IDV at various concentrations; one of the extended forms was effective in trans processing even at micromolar IDV.
Our results indicate that L domain-2, considered redundant in vitro, can undergo mutations in vivo that significantly alter PR function. These may contribute fitness benefits in both the absence and presence of PR inhibitor.
As HIV-1 evolves in infected individuals, genetic changes can alter the way the virus replicates, leading to lower replication rates, persistence and latency or, alternatively, higher replication rates and progressive disease. Mutations that confer resistance to protease (PR) inhibitors (PIs) are often accompanied by compensatory mutations in the gag gene at the PR cleavage sites in the encoded Gag substrate [1–4]. Mutations within the p6 region of Gag have also been linked to resistance [5–11], however, with the exception of a recent report , the evidence to date has been mainly correlative. The p6 region harbors determinants of viral budding; productive PR cleavages that ultimately lead to release of infectious particles are spatially and temporally linked to viral budding. Budding is directed by two motifs in p6 designated as the late (L) domains: The P7(T/S)AP motif (L domain-1) binds the ESCRT-I component Tsg101 [(ESCRT, endosomal sorting complex required for transport); [13–15]; reviewed in . The L35YPL(A/T/V/other)SL motif (L domain-2) binds the ESCRT adaptor protein Alix ([17–19]; reviewed in ). L domain-2 is thought to have a secondary redundant function when L domain-1 is intact because (1), all natural variants of the virus conserve PTAP while HIV subtype C lacks Y36, a major determinant of Alix binding and (2), mutations in PTAP impair virus release to a significant extent while the impact of LYPLXSL mutation is cell type-dependent.
We examined HIV-1 variants from 14 HIV-1 infected women with detectable HIV in both the genital tract and the plasma compartment despite a history of current or past antiretroviral therapy (ART). Eleven women reported receiving ART at the time point studied and three had been treated previously. We found, surprisingly, that half of the study participants harbored viruses that had mutations in L domain-2. Two substitutions that occurred at the same residue, S40F and S40A, were further investigated by engineering the individual substitutions into an HIV genome (pNL4-3). The substitution of F (S40F), where TCC = S to TTC = F, does not alter the sequence of the overlapping pol frame (TTC = F in the WT to TTT = F) but had an impact on several aspects of gag function [21–25]. The alanine (A) substitution (S40A), where TCC = S to GCC = A, in contrast, had a less deleterious effect on Gag processing, polyubiquitination or particle maturation morphology but altered the sequence at the p6*/PR cleavage site in the overlapping pol frame. In the wild type, the p6* reading frame ends with GTVSFNF/(N may be S, due to a natural polymorphism in subtype B); PR begins with\PQITLWQ. The S40A change in the gag frame is accompanied by an F to C change at the N-terminus of PR in the overlapping pol frame resulting in GTVSFN (or S) C/PQITLWQ. In the context of the NL4-3 strain of HIV-1, S40A altered sensitivity to the protease (PR) inhibitor Indinavir (IDV). In the context of a model precursor, the mutation inhibited mature PR liberation but produced extended PR forms with autoproteolytic activity that resisted IDV treatment. Thus, our investigation suggests that L domain-2, considered redundant in vitro, can evolve during disease progression to significantly alter PR function and contribute to viral fitness in both genital tract and plasma in vivo.
Polymorphisms* in L domain-2 of HIV-1 variants from infected individuals
Clinical history of variants encoding L domain-2 mutations
Viral RNA loadb
CD4 + countc
ABC, AZT, 3TC
[*As the parental sequence of the variants amplified from the clinical specimens is unknown, we defined here a polymorphism as a sequence variation identified in the Los Alamos National Laboratory (LANL) database at a frequency of at least 1 %. If the frequency is <1 %, the variation is regarded as a mutation.]
The S40A mutation in Gag alters the amino acid sequence at the p6*-PR cleavage site in the overlapping pol frame
Mutation of S40 to Ala (S40A) has little apparent effect on virus budding, maturation, or infectivity
The S40A mutation in NL4-3 alters PI sensitivity
IC50, IDV, relative change in IC50 compared to NL4-3 susceptibility
NL4-3 multi-drug resistant strain
The F56C mutation interferes with mature PR release from a model precursor
To test the supposition that the F56C mutation influences with PR autoprocessing in the presence of the PR inhibitor, we examined autoprocessing products in response to IDV treatment. A typical autoprocessing profile of the WT L-MBP precursor is shown in Fig. 6a (left). The upper Flag panel shows detection of L-MBP-Flag-p6* resulting from proximal site processing and L-MBP-Flag resulting from distal site processing. These processing reactions were not significantly suppressed by IDV at concentrations < 15 μM showing that precursor autoprocessing was significantly less sensitive than the mature PR to inhibition by PR inhibitors. The released PR-HA, which is the other product of proximal site processing, was readily detectable; its detection gradually decreased with increasing concentration of IDV mirroring the L-MBP-Flag-p6* profile (Fig. 6, panel B vs. C, solid lines). There were two p6*-PR-HA products reactive to anti-HA antibody (but not anti-Flag antibody, data not shown), referred to as p6*-PR-HA a and p6*-PR-HA b , respectively (Fig. 6a, lower HA panel). These two products were not detectable in the absence of PI but became detectable at high concentrations of IDV. Possibly, they were rapidly degraded in the absence of PR inhibitor but accumulated following treatment. Consequently, they displayed a bell-shaped detection profile reflecting a dynamic balance between production from distal processing with rapid degradation at low IDV, then preservation with increasing IDV, followed by loss of detection due to suppression of proximal processing at high IDV (Fig. 6, panels D and E, solid lines).
The F56C mutation showed several changes compared to the WT control. First, the proximal processing event was reduced as indicated by the fact that the corresponding products, L-MBP-Flag-p6* and PR-HA, were both reduced compared to the WT control (Fig. 6, panels B and C). Meanwhile, the distal processing event was not significantly affected or was slightly increased at low IDV concentrations, as indicated by L-MBP-Flag production (Fig. 6d). Interestingly, one of the distal processing products, p6*-PR-HA a , exhibited a profile distinct from that released from the WT control (Fig. 6e). This p6*-PR-HA a product appeared to resist self-degradation and remained at high detection levels with IDV up to 3 μM, unlike the one released from the WT precursor that was self-degraded in the absence of IDV and was only detectable at high concentrations of IDV. Thus, the F56C mutation altered the detection profile of the p6*-PR-HA a product. These data on their own do not rule out the possibility that mature PR degrades p6*-PR-HA a in trans as their detection profiles were inversely correlated. The other product, p6*-PR-HA b , displayed a response profile similar to the WT control in that it was not detectable at low IDV concentrations and became detectable with increasing IDV concentrations (Fig. 6f). Additionally, the p6*-PR-HA b released from the F56C mutant precursor appeared to maintain a high level of detection over a wider IDV range than that released from the WT control (Fig. 6f). These results support the supposition that the F56C mutation altered precursor autoprocessing leading to release of products with distinct enzymatic activities in the presence or absence of IDV. Collectively, the F56C mutation exerted multifaceted effects including (1) decreased production of free mature PR due to suboptimal substrate sequence, (2) release of an extended PR form (p6*-PR-HA a ) that is fairly stable over a wide range of IDV (Fig. 6e) and (3) increased detection of p6*-PR-HA b .
Trans cleavage by F56C-related PR products
The precursor PR has intrinsic but limited proteolytic activity [38, 39]. The initial production of mature PR must be catalyzed by the Gag-Pol precursor itself serving as both the substrate and enzyme. Therefore, we designed a trans cleavage assay to specifically assess how the precursor enzymes proteolytically process an engineered precursor substrate. The substrate (Myc-M2-PRD25N-V5) is a mini-precursor derivative carrying D25N mutation making it catalytically inactive, however, it has the WT proximal PR cleavage site and thus can be processed by a trans acting enzyme to release a PR(D25N)-V5 product (Fig. 7a). The amount of the released PR-V5 (band intensity) was quantified to represent trans processing efficiency (Fig. 7b).
We compared four L-MBP fusion precursors as the source of trans-acting enzymes in cells co-transfected with substrate- and enzyme- encoding plasmids. In addition to the full length fusion precursor, at least three species of HA-tagged enzymes were detected in cell lysates with each showing a different IDV response profile (Fig. 7d–g, upper panels). Trans cleavage of the substrate was observed with all four L-MBP enzymes, as indicated by detection of the 11 kDa PR-V5 product (Fig. 7d–g, lower panels). Surprisingly, the WT enzymes did not produce the highest level of the PR-V5 product, despite the high level of mature PR-HA detected in the lysate (Fig. 7d). This observation suggests that the mature PR-HA was not effective at trans processing the V5-tagged substrate under the test condition. Instead, PR-V5 production quantitatively coincides with p6*-PR-HAb detection signal (Additional file 1: Figure S1), suggesting that p6*-PRb might be responsible for trans processing the substrate. We were unable to exclusively define p6*-PRb as the sole enzyme for substrate trans processing due to simultaneous existing/detection of other HA antibody reactive enzymes. Possibilities exist that p6*-PRb enzyme functions synergistically with PR and/or p6*-PRa. Another interesting observation was that the two F56C-containing constructs released more PR-V5 products than the WT and F56 V controls and demonstrated consistently high levels of trans processing products even at micromolar concentrations of IDV (Fig. 7b). In contrast, the WT p6*-PR-HAb enzyme appeared to be less efficient at trans cleavage despite of its high level of detection (Fig. 7c). The F56V control was unable to release any mature PR and only produced p6*-PR-HAa and p6*-PR-HAb enzymes that demonstrated low level of trans processing (Fig. 7g). Thus, the data indicate that the p6*-PR-HAb enzyme was enzymatically active at trans processing the engineered substrate. Furthermore, the p6*-PR-HAb enzyme in the context of F56C content was more effective at trans processing than those made by the WT or F56 V controls while being resistant to IDV suppression. These findings are consistent with a previous report demonstrating that the p6*-PR fragment can process some Gag cleavage sites and many Pol cleavage sites in the absence of mature PR . Our results support the notion that p6*-PR enzymes that remain active could prevent ART from suppressing virus replication by functionally complementing reduced production and activity of mature PR, permitting effective Gag processing and formation of mature infectious viral particles.
The retention and conservation of L domain-2, which appears to be redundant for the budding function, implies that it provides the virus with an evolutionary advantage. It should be noted that treatment with ART or HAART has been shown to reduce the HIV load more effectively and render the viral load undetectable more often in the genital tract than in the plasma [41, 42]. The women whom we studied were specifically chosen because they had detectable HIV in the CVL and plasma even though they had a history of ART. Possibly, the L domain-2 mutations contributed to the viral loads in both compartments, especially the CVL, as suggested by the findings for the polymorphism at the S40 site.
In participant #36, variants carrying the S40A mutation were the predominant variants detected in both vaginal fluid and plasma at the time of sampling. It is possible that the mutation was selected under pressure of ART after infection with virus that had S at Gagp6 position 40 or, alternatively, existed in the region and as a virus with the S40A already in it. As the study participants were HA/ART-experienced, we hypothesized that the selection and maintenance of the mutation possibly contributed to escape from drug-mediated inhibition. Sequencing pre- and post- ART samples archived from this participant might provide evidence that the mutation was selected during ART, thereby supporting a relationship. Interestingly, of all the possible nineteen other amino acid changes, we found only Ala and Phe (S40A and S40F, respectively) as substitutions for the Ser40 residue. These residues do not have chemical or physical properties in common and one (Ala) altered the amino acid sequence in the overlapping pol frame while the other did not. In addition, while the substitution of F (TT/C) for S (TC/C) requires a pyrimidine to pyrimidine substitution (i.e., a transition), the change from S to A (GC/C) that was detected more frequently requires a pyrimidine to purine substitution (i.e., a transversion) which occurs much less frequently than transitions. We therefore hypothesized that the Ser40 to Ala change was not random but, rather, was strongly selected.
Our studies indicate that the Gag p6 S40A/Pol p6*-PR F56C mutation reduces mature PR production but enhances production, stability, and activity of extended forms of PR. Interestingly, these extended PRs were somehow resistant to (at least one) PR inhibitor, IDV and one of them (p6*-PRb) can effectively process an engineered substrate in trans in a model system. The role of p6*-PRa remains undefined at the moment but it might provide some fitness advantages of processing other cleavages sites in the presence of IDV and reduced amounts of mature PR. It should be noted that these extended PRs were also observed with the WT control but at low levels compared to the F56C construct (Fig. 6a, the lower panels). This is consistent with previous reports showing the distal site located in the N-terminal region of p6* as another autoprocessing target [32, 33, 36]. Both the proximal and distal sites appear to be processed equally in the context of WT fusion precursor. With the F56C mutation, the proximal site processing was suppressed and the distal site processing was slightly enhanced (Fig. 6a, the upper panels), which could contribute to the enhanced production of extended PRs. The fact that IDV was the PR inhibitor on record as prescribed during the period when virus suppression failed is highly suggestive that the observations are linked to treatment outcome, although caveats related to adherence remain. The reduced production of mature PR resulting from the mutation nevertheless permitted formation of infectious particles. This is consistent with results obtained in our previous study and herein where NL4-3 bearing the S40A mutation produced viral particles showing normal Gag processing and mature virion morphology . We speculate that the extended forms of PR could be attributed to this apparent paradox. There are at least 10 cleavage sites in the Gag and Gag/Pol polyproteins that are recognized and processed by the mature PR. Proper processing of these sites is vital to producing infectious progeny. The extended forms of PR are capable of processing some but not all of the cleavage sites . The F56C mutation somehow stabilizes the extended forms (both a and b forms) over a wide range of IDV concentrations (Fig. 6b). By being resistant to IDV inhibition, these enzymes are able to accomplish the “heavy lifting” of processing many cleavage sites in the presence of PIs, leaving much lower levels of cleavage sites that have to be processed by the mature PR. Such a tag-teaming strategy might be used by the virus to achieve ART resistance. We noticed that only F56C-containing enzymes demonstrated enhanced trans processing efficiency while the F56 V control showed a much reduced trans processing efficiency (Fig. 7b). Therefore, there might be something unique about the F56C mutation that not only permits it to stabilize the extended PR but also to enhance their proteolysis. Consistent with this speculation, a highly resistant patient-derived strain also contains the same F56C mutation in the proximal cleavage site . It of course will require more molecular and structural investigation to determine how the F56C mutation can regulate activities of the extended PRs while conferring drug resistance at the same time.
It should be borne in mind that the infectivity assay and the assays using model precursors do not test PR inhibitor susceptibility in the same manner. In the infectivity assay, the virus is producing both extended and mature PR forms. As the drug inhibits both but was optimized to target the mature PR, the outcome may mainly reflect impact on this form. In contrast, the assay model allows detection of the model precursors and uses conditions optimized to detect the extended forms. The observed reduced production of mature PR due to the F56C mutation compared to the level made by the parental NL4-3 virus provides an explanation for the mutant’s increased IDV resistance: Extended PR forms that are not the authentic target of the drug were generated instead of the mature PR. It should be noted that although our study characterizes the virologic consequences of the S40A polymorphism detected in the HIV-1 variants that were obtained directly from the WIHS participants, the limitations imposed by the number of subjects examined and details of treatment prevent us from drawing direct conclusions about the effect of the mutations on susceptibility to the participant’s treatment. Nevertheless, the S40A mutation impacted drug susceptibility. This change in L domain-2 could thus influence treatment success or failure through its influence on viral PR.
A final point is that, in contrast to the S40A mutation where much of the viral infectivity was maintained, the S40F mutation reduced the infectivity as shown here and previously [21, 22] yet the mutation was maintained in the infected host quasispecies. This suggests that both S40A, which impacts PR, and S40F, which does not, provide a fitness advantage. A possible explanation for an advantage of the S40F mutation might be the formation of filopodia-like structures, a phenotype that was exhibited by that mutant . This might support cell-to-cell transmission of the virus which might permit escape from ART surveillance. We examined this possibility and found that, indeed, the infectivity defect of the S40F mutant was rescued by cell-to-cell transmission (Additional file 2: Figure S2). This rescue did not prevent effective inhibition of virus transmission by IDV.
Changes in HIV’s L domain-2 can significantly alter PR function and affect fitness. The alteration of PR activity can contribute to productive replication in both the absence and presence of PR inhibitor.
The Women’s Interagency HIV Study (WIHS) is a multicenter, longitudinal observational cohort study of HIV-1 infection of women . To compare HIV-1 drug resistance in the plasma and genital tract of women treated with antiretroviral agents, we previously analyzed 308 HIV-1 sequences from 14 women in the Bronx and Brooklyn WIHS (described in ; GenBank accession numbers DQ372126-DQ372434). Three hundred eight unique sequences (280 single variants and 28 population sequences) encompassed the protease–reverse transcriptase (PR-RT) region of the pol gene. We examined the overlapping gag frame of these 308 sequences for L domain-2 mutations and identified such mutations in 7 of the 14 subjects. To examine the prevalence of L domain-2 mutations in a larger group of subjects in the WIHS during approximately the same time period, we examined population-based sequences of the overlapping gag frame from an additional 90 women who exhibited plasma viremia after ≥18 consecutive months of ART (Weiser, unpublished data). Details of sample collection and analysis; viral RNA isolation, reverse-transcriptase polymerase chain reaction (RT-PCR) and sequencing; and genotypic drug resistance are described in Kemal et al. .
Plasmids and reagents
Plasmids encoding pNL4-3∆Env , pNL4-3-S40A or pNL4-3-S40F , and l-maltose binding protein (MBP)-FLAG-p6*-PR-HA  have been previously described. Additional NL4-3 mutants were constructed by site-directed mutagenesis (Agilent). NL4-3 multiple drug resistance (MDR) encodes the following mutations that together contribute to high PR and reverse transcriptase (RT) drug resistance (PR: 3I, 10I, 33F, 41K, 46I, 54L, 62V, 63P, 66L, 73A, 77I, 84V, 89M, 90M and 93L; RT: 41L,43E, 44A,68G,74V, 118F, 165I, 181C, 190S, 203Q, 207D, 210W, 211K, 214F, 215Y, 219N, 250E). Fusion precursor mutations tested in this study were introduced to the L-MBP-Flag-p6*-PR-HA by site-directed mutagenesis as previously described . The trans processing substrate, Myc-M2-PRD25N-V5, was engineered the same way. In this construct, the Myc epitope (EQKLISEEDL) is led by a tripeptide (MAS) followed by M2 (a truncated version of p6* lacking the N-terminal 20 amino acids as reported previously ) fused to PR-V5 . All the constructs were verified by sequencing analysis.
Fusion precursors and the autoprocessing products were detected using antibodies raised against the HA tag (H9658, Sigma), the V5 tag (600-401-378, Rockland Laboratory) or the FLAG (F1804, Sigma). GAPDH levels were detected by using antibody from Millipore (MAB374) following manufacturer’s recommendation. The diluted primary antibodies were pre-absorbed against PVDF membranes coated with cell lysates made of mock transfected cells to remove background noise. Indinavir (Cat# 8145) was obtained from the AIDS Research and Reference Reagent Program, Division of AIDS, NIAID, NIH.
ELISA and MAGI assays
The concentration of p24 in the filtered medium was determined by ELISA (enzyme-linked immunosorbent assay; Immunodiagnostics, Inc.) and equivalent amounts of p24 were used to infect MAGI (multinuclear activation of a β-galactosidase indicator) cells to determine infectivity .
293T cells were co-transfected with pNL4-3∆Env and pVSV G. Where indicated, PR inhibitor was added. Viral supernatants were collected at 48 h post-transfection, clarified by centrifugation, and stored at −80 °C. Serially diluted viral supernatants were used to infect 1 × 104 TZM-bl cells in 96-well plates. Infections were done in triplicate with viral supernatants for each construct tested. Infectivity of the virus particles in TZM-bl cells was assessed at 48 h post-infection using a Galacto-Star System for detecting β-galactosidase activity (Applied Biosystems), as described previously [46, 47].
Cells grown on ACLAR film were fixed in 4 % paraformaldehyde/0.1 % EM grade glutaraldehyde in PBS, soaked in 2 % osmium tetroxide, dehydrated in a graded series of ethyl alcohol solutions and embedded in Durpan resin. Eighty nm ultrathin sections were counterstained with uranyl acetate and lead citrate and viewed with a FEI Tecanal BioTwinG2 electron microscope.
Cell culture, transfection, western blots, auto-processing and trans processing assays
HEK293T cells were maintained in DMEM (Dulbecco’s Modified Eagle’s Medium; Invitrogen, Carlsbad, CA) culture media containing penicillin and streptomycin and 10 % fetal bovine serum as previously described [32, 36]. Transfection of HEK293T cells was previously described [32, 36]. In brief, cells were seeded onto 6-well, 12-well or 24-well plates the day before transfection. Cell confluency at time of transfection was 40–60 %. The total amount of transfected DNA per well was 1 μg (6-well plate), 0.5 μg (12-well plate), or 0.25 μg (24 well plate). Plasmids pcDNA and peGFP (mixed at 19:1) were used as a mock transfection control. For trans processing analysis, a DNA ratio of enzyme to substrate was consistently maintained at 1:1 unless otherwise stated. For a 12-well plate, autoclaved H2O was added to bring the DNA + H2O volume to 65.7 μL. To this, 9.3 μL 2 M CaCl2 and 75 μL 2× HBS (50 mM HEPES, 280 mM NaCl, 10 mM KCl, 12 mM Dextrose, 1.5 mM Na2HPO4, pH 7.05) were added drop-wise for a total volume of 150 μL. The mixture was then added to plated cells and incubated for 7–11 h before the media was changed. For experiments involving drug treatment, drug was added at this point to the desired concentration. At 24–30 h post transfection, cells were washed with PBS once and lysed with 80 μL Lysis Buffer/PI solution (Lysis Buffer A: 25 mM Tris–HCl pH 8.0, 150 mM NaCl, 1 % sodium deoxycholate, 1 % Triton X-100 plus 1× PI cocktail). Cell debris was removed and lysate was transferred to a microcentrifuge tube containing 15 μL 6× SDS Loading Buffer (60 % glycerol, 0.6 M DTT powder, 6 % SDS, 0.006 % Bromophenol blue, 0.35 M Tris–HCl, water). Samples were boiled for 3–5 min, resolved by 13 % SDS-PAGE, and transferred onto a PVDF membrane (Millipore 0.45 µ Immobilon P). Fluorescently labeled secondary antibodies IR800 goat anti-mouse (Li-COR cat#926-32210), IR700 goat anti-mouse (Li-COR cat#926-68020), and IR800 goat anti-rabbit (Li-COR cat#926-32211) were used to visualize Western blot images by an Odyssey infrared dual laser scanning unit (LI-COR Biotechnology, Lincoln, Nebraska). ImageStudio® quantification program (Li-COR Biotechnology) was used to determine signal intensity for a given band by measuring the total value of the band area minus local background noise times band area. Trans processing efficiency was reflected by PR-V5 intensity.
SW, VS, ND, BRK, SM, KSK, BS, BF, HL and KA were involved in data acquisition; SW, VS, BKC, BW, HB, CC and CAC participated in data analysis and interpretation; SW, CC, and CAC designed the studies and wrote the manuscript; all authors read and approved the final manuscript
We thank Drs. James Konopka, Stony Brook University, and Wei Gao, Albert Einstein College of Medicine, for assistance with retrieval of databank sequences. We thank Dr. Lorna Ehrlich (Stony Brook University) for helpful discussions.
The authors declare that they have no competing interests.
This study was supported by NIH NIAID R01 AI068463, NIGMS R01 GM111028 and ARRA supplemental funding awards to CAC; NIH NIAID R21A1080351 and R03AI108392 to CC; NIH U3U 01AI35004 award to KA and NIH NIAID R01AI064001 to VS; NIH R01GM113885 to BKC. The content is solely the responsibility of the authors.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Adamson CS, Sakalian M, Salzwedel K, Freed EO. Polymorphisms in Gag spacer peptide 1 confer varying levels of resistance to the HIV-1 maturation inhibitor bevirimat. Retrovirology. 2010;7:36.View ArticlePubMedPubMed CentralGoogle Scholar
- Lefebvre E, Schiffer CA. Resilience to resistance of HIV-1 protease inhibitors: profile of darunavir. AIDS Rev. 2008;10:131–42.PubMedPubMed CentralGoogle Scholar
- Wensing AM, van Maarseveen NM, Nijhuis M. Fifteen years of HIV Protease Inhibitors: raising the barrier to resistance. Antiviral Res. 2010;85:59–74.View ArticlePubMedGoogle Scholar
- Fun A, Wensing AM, Verheyen J, Nijhuis M. Human immunodeficiency virus Gag and protease: partners in resistance. Retrovirology. 2012;9:63.View ArticlePubMedPubMed CentralGoogle Scholar
- Brumme ZL, Chan KJ, Dong WW, Wynhoven B, Mo T, Hogg RS, Montaner JS, O’Shaughnessy MV, Harrigan PR. Prevalence and clinical implications of insertions in the HIV-1 p6Gag N-terminal region in drug-naive individuals initiating antiretroviral therapy. Antivir Ther. 2003;8:91–6.PubMedGoogle Scholar
- Gallego O, de Mendoza C, Corral A, Soriano V. Changes in the human immunodeficiency virus p7-p1-p6 gag gene in drug-naive and pretreated patients. J Clin Microbiol. 2003;41:1245–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Ibe S, Shibata N, Utsumi M, Kaneda T. Selection of human immunodeficiency virus type 1 variants with an insertion mutation in the p6(gag) and p6(pol) genes under highly active antiretroviral therapy. Microbiol Immunol. 2003;47:71–9.View ArticlePubMedGoogle Scholar
- Lastere S, Dalban C, Collin G, Descamps D, Girard PM, Clavel F, Costagliola D, Brun-Vezinet F. Impact of insertions in the HIV-1 p6 PTAPP region on the virological response to amprenavir. Antivir Ther. 2004;9:221–7.PubMedGoogle Scholar
- Marlowe N, Flys T, Hackett J Jr, Schumaker M, Jackson JB, Eshleman SH. Analysis of insertions and deletions in the gag p6 region of diverse HIV type 1 strains. AIDS Res Hum Retroviruses. 2004;20:1119–25.View ArticlePubMedGoogle Scholar
- Martins AN, Arruda MB, Pires AF, Tanuri A, Brindeiro RM. Accumulation of P(T/S)AP late domain duplications in HIV type 1 subtypes B, C, and F derived from individuals failing ARV therapy and ARV drug-naive patients. AIDS Res Hum Retroviruses. 2011;27:687–92.View ArticlePubMedGoogle Scholar
- Peters S, Munoz M, Yerly S, Sanchez-Merino V, Lopez-Galindez C, Perrin L, Larder B, Cmarko D, Fakan S, Meylan P, Telenti A. Resistance to nucleoside analog reverse transcriptase inhibitors mediated by human immunodeficiency virus type 1 p6 protein. J Virol. 2001;75:9644–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Martins AN, Waheed AA, Ablan SD, Huang W, Newton A, Petropoulos CJ, Brindeiro RM, Freed EO. Elucidation of the molecular mechanism driving duplication of the HIV-1 PTAP Late domain. J Virol. 2015;90:768–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Garrus JE, von Schwedler UK, Pornillos OW, Morham SG, Zavitz KH, Wang HE, Wettstein DA, Stray KM, Cote M, Rich RL, Myszka DG, Sundquist WI. Tsg101 and the vacuolar protein sorting pathway are essential for HIV-1 budding. Cell Dev Biol. 2001;107:55–65.Google Scholar
- VerPlank L, Bouamr F, LaGrassa TJ, Agresta B, Kikonyogo A, Leis J, Carter CA. Tsg101, a homologue of ubiquitin-conjugating (E2) enzymes, binds the L domain in HIV type 1 Pr55Gag. Proc Natl Acad Sci USA. 2001;98:7724–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin-Serrano J, Neil SJ. Host factors involved in retroviral budding and release. Nat Rev Microbiol. 2011;9:519–31.View ArticlePubMedGoogle Scholar
- Hurley JH. ESCRTs are everywhere. EMBO J. 2015;34:2398–407.View ArticlePubMedGoogle Scholar
- Martin-Serrano J, Yarovoy A, Perez-Caballero D, Bieniasz PD. Divergent retroviral late-budding domains recruit vacuolar protein sorting factors by using alternative adaptor proteins. Proc Natl Acad Sci USA. 2003;100:12414–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Strack B, Calistri A, Craig S, Popova E, Gottlinger HG. AIP1/ALIX is a binding partner for HIV-1 p6 and EIAV p9 functioning in virus budding. Cell. 2003;114:689–99.View ArticlePubMedGoogle Scholar
- von Schwedler UK, Stuchell M, Muller B, Ward DM, Chung HY, Morita E, Wang HE, Davis T, He GP, Cimbora DM, et al. The protein network of HIV budding. Cell. 2003;114:701–13.View ArticleGoogle Scholar
- Bissig C, Gruenberg J. ALIX and the multivesicular endosome: ALIX in Wonderland. Trends Cell Biol. 2014;24:19–25.View ArticlePubMedGoogle Scholar
- Votteler J, Neumann L, Hahn S, Hahn F, Rauch P, Schmidt K, Studtrucker N, Solbak SM, Fossen T, Henklein P, et al. Highly conserved serine residue 40 in HIV-1 p6 regulates capsid processing and virus core assembly. Retrovirology. 2011;8:11.View ArticlePubMedPubMed CentralGoogle Scholar
- Watanabe SM, Chen MH, Khan M, Ehrlich L, Kemal KS, Weiser B, Shi B, Chen C, Powell M, Anastos K, et al. The S40 residue in HIV-1 Gag p6 impacts local and distal budding determinants, revealing additional late domain activities. Retrovirology. 2013;10:143.View ArticlePubMedPubMed CentralGoogle Scholar
- Hahn F, Setz C, Friedrich M, Rauch P, Solbak SM, Froystein NA, Henklein P, Votteler J, Fossen T, Schubert U. Mutation of the highly conserved Ser-40 of the HIV-1 p6 gag protein to Phe causes the formation of a hydrophobic patch, enhances membrane association, and polyubiquitination of Gag. Viruses. 2014;6:3738–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Kudoh A, Takahama S, Sawasaki T, Ode H, Yokoyama M, Okayama A, Ishikawa A, Miyakawa K, Matsunaga S, Kimura H, et al. The phosphorylation of HIV-1 Gag by atypical protein kinase C facilitates viral infectivity by promoting Vpr incorporation into virions. Retrovirology. 2014;11:9.View ArticlePubMedPubMed CentralGoogle Scholar
- Radestock B, Burk R, Muller B, Krausslich HG. Re-visiting the functional Relevance of the highly conserved Serine 40 Residue within HIV-1 p6(Gag). Retrovirology. 2014;11:114.View ArticlePubMedPubMed CentralGoogle Scholar
- Kemal KS, Burger H, Mayers D, Anastos K, Foley B, Kitchen C, Huggins P, Schroeder T, Picchio G, Back S, et al. HIV-1 drug resistance in variants from the female genital tract and plasma. J Infect Dis. 2007;195:535–45.View ArticlePubMedGoogle Scholar
- Vincent O, Rainbow L, Tilburn J, Arst HN Jr, Penalva MA. YPXL/I is a protein interaction motif recognized by aspergillus PalA and its human homologue, AIP1/Alix. Mol Cell Biol. 2003;23:1647–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Krausslich HG, Wimmer E. Viral proteinases. Annu Rev Biochem. 1988;57:701–54.View ArticlePubMedGoogle Scholar
- Konvalinka J, Krausslich HG, Muller B. Retroviral proteases and their roles in virion maturation. Virology. 2015;479–480:403–17.View ArticlePubMedGoogle Scholar
- Teira R, Camara MM, Escobar A, Munoz P, Lopez de Munain J, Santamaria JM. A randomized clinical trial to compare the effectiveness of indinavir, ritonavir and saquinavir. Med Clin (Barc). 2003;121:532–4.View ArticleGoogle Scholar
- Huang L, Li Y, Chen C. Flexible catalytic site conformations implicated in modulation of HIV-1 protease autoprocessing reactions. Retrovirology. 2011;8:79.View ArticlePubMedPubMed CentralGoogle Scholar
- Counts CJ, Ho PS, Donlin MJ, Tavis JE, Chen C. A functional interplay between human immunodeficiency virus type 1 protease residues 77 and 93 involved in differential regulation of precursor autoprocessing and mature protease activity. PLoS ONE. 2015;10:e0123561.View ArticlePubMedPubMed CentralGoogle Scholar
- Shafer RW, Schapiro JM. HIV-1 drug resistance mutations: an updated framework for the second decade of HAART. AIDS Rev. 2008;10:67–84.PubMedPubMed CentralGoogle Scholar
- Fisher R, van Zyl GU, Travers SA, Kosakovsky Pond SL, Engelbrech S, Murrell B, Scheffler K, Smith D. Deep sequencing reveals minor protease resistance mutations in patients failing a protease inhibitor regimen. J Virol. 2012;86:6231–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang L, Chen C. Autoprocessing of human immunodeficiency virus type 1 protease miniprecursor fusions in mammalian cells. AIDS Res Ther. 2010;7:27.View ArticlePubMedPubMed CentralGoogle Scholar
- Mildner AM, Rothrock DJ, Leone JW, Bannow CA, Lull JM, Reardon IM, Sarcich JL, Howe WJ, Tomich CS, Smith CW, et al. The HIV-1 protease as enzyme and substrate: mutagenesis of autolysis sites and generation of a stable mutant with retained kinetic properties. Biochemistry. 1994;33:9405–13.View ArticlePubMedGoogle Scholar
- Louis JM, Nashed NT, Parris KD, Kimmel AR, Jerina DM. Kinetics and mechanism of autoprocessing of human immunodeficiency virus type 1 protease from an analog of the Gag-Pol polyprotein. Proc Natl Acad Sci USA. 1994;91:7970–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Pettit SC, Clemente JC, Jeung JA, Dunn BM, Kaplan AH. Ordered processing of the human immunodeficiency virus type 1 GagPol precursor is influenced by the context of the embedded viral protease. J Virol. 2005;79:10601–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Ludwig C, Leiherer A, Wagner R. Importance of protease cleavage sites within and flanking human immunodeficiency virus type 1 transframe protein p6* for spatiotemporal regulation of protease activation. J Virol. 2008;82:4573–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Graham SM, Jalalian-Lechak Z, Shafi J, Chohan V, Deya RW, Jaoko W, Mandaliya KN, Peshu NM, Overbaugh J, McClelland RS. Antiretroviral treatment interruptions predict female genital shedding of genotypically resistant HIV-1 RNA. J Acquir Immune Defic Syndr. 2012;60:511–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Cu-Uvin S, Caliendo AM, Reinert S, Chang A, Juliano-Remollino C, Flanigan TP, Mayer KH, Carpenter CC. Effect of highly active antiretroviral therapy on cervicovaginal HIV-1 RNA. AIDS. 2000;14:415–21.View ArticlePubMedGoogle Scholar
- Kozisek M, Henke S, Saskova KG, Jacobs GB, Schuch A, Buchholz B, Muller V, Krausslich HG, Rezacova P, Konvalinka J, Bodem J. Mutations in HIV-1 gag and pol compensate for the loss of viral fitness caused by a highly mutated protease. Antimicrob Agents Chemother. 2012;56:4320–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Anastos K, Gange SJ, Lau B, Weiser B, Detels R, Giorgi JV, Margolick JB, Cohen M, Phair J, Melnick S, et al. Association of race and gender with HIV-1 RNA levels and immunologic progression. J Acquir Immune Defic Syndr. 2000;24:218–26.View ArticlePubMedGoogle Scholar
- He J, Landau N. Use of a novel human immunodeficiency virus type 1 reporter virus expressing human placental alkaline phosphatase to detect an alternative viral receptor. J Virol. 1995;69:4587–92.PubMedPubMed CentralGoogle Scholar
- Harari A, Ooms M, Mulder LC, Simon V. Polymorphisms and splice variants influence the antiretroviral activity of human APOBEC3H. J Virol. 2009;83:295–303.View ArticlePubMedGoogle Scholar
- Simon V, Zennou V, Murray D, Huang Y, Ho DD, Bieniasz PD. Natural variation in Vif: differential impact on APOBEC3G/3F and a potential role in HIV-1 diversification. PLoS Pathog. 2005;1:e6.View ArticlePubMedPubMed CentralGoogle Scholar
- Durham ND, Chen BK. Measuring T cell-to-T cell HIV-1 transfer, viral fusion, and infection using flow cytometry. Methods Mol Biol. 2016;1354:21–38.View ArticlePubMedGoogle Scholar