Evolution of human immunodeficiency virus type 1 protease genotypes and phenotypes in vivo under selective pressure of the protease inhibitor ritonavir

Inferring viral population structures using heteroduplex mobility and DNA sequence analyses

Heteroduplex mobility (HMA) and tracking assays (HTA) are used to assess genetic relationships between DNA molecules. While distinguishing relationships between clonal or nearly clonal molecules is relatively straightforward, inferring population structures is more complex. To address this issue, HIV-1 quasispecies with varying levels of diversity were studied using both HTA and DNA sequencing. Viral diversity estimates and the temporal features of virus evolution were found to be generally concordant between HTA and DNA sequencing. In addition, the distribution of pairwise differences and the rates of virus divergence were similar between the two methods. These findings support the use of HTA to characterize variant populations of DNA and strengthen previous inferences concerning the evolution of HIV-1 over the course of infection.

Selection of Drug-Resistant Feline Immunodeficiency Virus (FIV) Encoding FIV/HIV Chimeric Protease in the Presence of HIV-Specific Protease Inhibitors

An infectious chimeric feline immunodeficiency virus (FIV)/HIV strain carrying six HIV-like protease (PR) mutations (I37V/N55M/V59I/I98S/Q99V/P100N) was subjected to selection in culture against the PR inhibitor lopinavir (LPV), darunavir (DRV), or TL-3. LPV selection resulted in the sequential emergence of V99A (strain S-1X), I59V (strain S-2X), and I108V (strain S-3X) mutations, followed by V37I (strain S-4X). Mutant PRs were analyzed in vitro , and an isogenic virus producing each mutant PR was analyzed in culture for LPV sensitivity, yielding results consistent with the original selection. The 50% inhibitory concentrations (IC 50 s) for S-1X, S-2X, S-3X, and S-4X were 95, 643, 627, and 1,543 nM, respectively. The primary resistance mutations, V99 82 A, I59 50 V, and V37 32 I, are consistent with the resistance pattern developed by HIV-1 under similar selection conditions. While resistance to LPV emerged readily, similar PR mutations causing resistance to either DRV or TL-3 failed to emerge after passage for more than a year. However, a G37D mutation in the nucleocapsid (NC) was observed in both selections and an isogenic G37D mutant replicated in the presence of 100 nM DRV or TL-3, whereas parental chimeric FIV could not. An additional mutation, L92V, near the PR active site in the folded structure recently emerged during TL-3 selection. The L92V mutant PR exhibited an IC 50 of 50 nM, compared to 35 nM for 6s-98S PR, and processed the NC-p2 junction more efficiently, consistent with increased viral fitness. These findings emphasize the role of mutations outside the active site of PR in increasing viral resistance to active-site inhibitors and suggest additional targets for inhibitor development.

Human Immunodeficiency Virus Gag and protease: partners in resistance

Human Immunodeficiency Virus (HIV) maturation plays an essential role in the viral life cycle by enabling the generation of mature infectious virus particles through proteolytic processing of the viral Gag and GagPol precursor proteins. An impaired polyprotein processing results in the production of non-infectious virus particles. Consequently, particle maturation is an excellent drug target as exemplified by inhibitors specifically targeting the viral protease (protease inhibitors; PIs) and the experimental class of maturation inhibitors that target the precursor Gag and GagPol polyproteins. Considering the different target sites of the two drug classes, direct cross-resistance may seem unlikely. However, coevolution of protease and its substrate Gag during PI exposure has been observed both in vivo and in vitro. This review addresses in detail all mutations in Gag that are selected under PI pressure. We evaluate how polymorphisms and mutations in Gag affect PI therapy, an aspect of PI resistance that is currently not included in standard genotypic PI resistance testing. In addition, we consider the consequences of Gag mutations for the development and positioning of future maturation inhibitors.

Single genome sequencing of HIV-1 gag and protease resistance mutations at virologic failure during the OK04 trial of simplified versus standard maintenance therapy

BACKGROUND

Ritonavir-boosted lopinavir (LPV/RTV) alone has been evaluated as simplified maintenance therapy for HIV-1 infection, but there are concerns about greater potential for emergence of protease inhibitor (PI) resistance. The OK04 trial evaluated maintenance therapy with LPV/RTV alone versus standard therapy (ST) with two NRTIs plus LPV/RTV in 205 patients, of whom 15 had virological rebound by week 48 (11 versus 4 patients, respectively). We developed a single genome sequencing (SGS) assay of HIV-1 gag and protease to assess the emergence of low frequency drug-resistant variants during virological rebound.

METHODS

Plasma samples from 15 subjects at virological rebound were analysed by SGS of HIV-1 gag and protease genes. A total of 45 SGS sequences were planned per sample, providing 90% power to detect variants comprising >5% of the virus population.

RESULTS

Overall, 521 single sequences obtained from 13 patients (range 4-48 sequences/patient) revealed similar frequencies of major protease resistance mutations in samples from the LPV/RTV alone (3/11) and ST (3/4) arms (P=0.10), with a median number of minor protease resistance mutations of 3.0 versus 3.5, respectively (P=0.23). Median number of gag PI resistance mutations were similar between the LPV/RTV alone and ST arms at cleavage sites (3.0 versus 2.5; P=0.83), non-cleavage sites (21 versus 16.5; P=0.71) and the transframe protein-p6 pol region cleavage sites (4.0 versus 3.0; P=0.6).

CONCLUSIONS

Although more subjects with simplified maintenance therapy with LPV/RTV alone had virological rebound compared to the ST arm, this was not associated with more frequent emergence of variants encoding PI resistance mutations in gag or protease detected by SGS.

Accurate sampling and deep sequencing of the HIV-1 protease gene using a Primer ID

Viruses can create complex genetic populations within a host, and deep sequencing technologies allow extensive sampling of these populations. Limitations of these technologies, however, potentially bias this sampling, particularly when a PCR step precedes the sequencing protocol. Typically, an unknown number of templates are used in initiating the PCR amplification, and this can lead to unrecognized sequence resampling creating apparent homogeneity; also, PCR-mediated recombination can disrupt linkage, and differential amplification can skew allele frequency. Finally, misincorporation of nucleotides during PCR and errors during the sequencing protocol can inflate diversity. We have solved these problems by including a random sequence tag in the initial primer such that each template receives a unique Primer ID. After sequencing, repeated identification of a Primer ID reveals sequence resampling. These resampled sequences are then used to create an accurate consensus sequence for each template, correcting for recombination, allelic skewing, and misincorporation/sequencing errors. The resulting population of consensus sequences directly represents the initial sampled templates. We applied this approach to the HIV-1 protease (pro) gene to view the distribution of sequence variation of a complex viral population within a host. We identified major and minor polymorphisms at coding and noncoding positions. In addition, we observed dynamic genetic changes within the population during intermittent drug exposure, including the emergence of multiple resistant alleles. These results provide an unprecedented view of a complex viral population in the absence of PCR resampling.

Generation of infectious feline immunodeficiency virus (FIV) encoding FIV/human immunodeficiency virus chimeric protease

Feline immunodeficiency virus (FIV) and human immunodeficiency virus type 1 (HIV-1) proteases (PRs) share only 23% amino acid identity and exhibit distinct specificities yet have very similar 3-dimensional structures. Chimeric PRs in which HIV residues were substituted in structurally equivalent positions in FIV PR were prepared in order to study the molecular basis of PR specificity. Previous in vitro analyses showed that such substitutions dramatically altered the inhibitor specificity of mutant PRs but changed the rate and specificity of Gag cleavage so that chimeric FIVs were not infectious. Chimeric PRs encoding combinations of the I37V, N55M, M56I, V59I, L97T, I98P, Q99V, and P100N mutations were cloned into FIV Gag-Pol, and those constructs that best approximated the temporal cleavage pattern generated by wild-type FIV PR, while maintaining HIV-like inhibitor specificity, were selected. Two mutations, M56I and L97T, were intolerant to change and caused inefficient cleavage at NC-p2. However, a mutant PR with six substitutions (I37V, N55M, V59I, I98P, Q99V, and P100N) was selected and placed in the context of full-length FIV-34TF10. This virus, termed YCL6, had low-level infectivity ex vivo, and after passage, progeny that exhibited a higher growth rate emerged. The residue at the position of one of the six mutations, I98P, further mutated on passage to either P98H or P98S. Both PRs were sensitive to the HIV-1 PR inhibitors lopinavir (LPV) and darunavir (DRV), as well as to the broad-based inhibitor TL-3, with 50% inhibitory concentrations (IC(50)) of 30 to 40 nM, consistent with ex vivo results obtained using mutant FIVs. The chimeras offer an infectivity system with which to screen compounds for potential as broad-based PR inhibitors, define structural parameters that dictate specificity, and investigate pathways for drug resistance development.

HIV protease resistance and viral fitness

PURPOSE OF REVIEW

This review focuses on the evolution of protease inhibitor resistance and replication capacity in the presence and absence of protease inhibitor pressure.

RECENT FINDINGS

Classically, HIV escapes through mutations in the protease itself causing a decrease in affinity to the inhibitor, leading to resistance. These changes also affect the binding of the enzyme to the natural substrate, and as a consequence cause a decrease in replication capacity of the virus. Continuous replication of these viruses may result in the acquisition of compensatory changes, which will fixate the drug-resistant variant in the viral population. Furthermore, novel treatment strategies have been developed to combat the development of classic protease inhibitor resistance. Using these strategies, the development of resistance in the viral protease is blocked because single or double mutations do not confer significant resistance. Alternative protease inhibitor resistance pathways are described, which enable the virus to escape these novel strategies.

SUMMARY

Suboptimal protease inhibitor pressure clearly results in the selection of mutations conferring resistance and in the acquisition of mutations compensating the initial reduction in viral replicative capacity. The major implications of the selection of these compensatory changes on evolution in the absence of protease inhibitor pressure are discussed.

A phylogenetic and Markov model approach for the reconstruction of mutational pathways of drug resistance

MOTIVATION

Modern HIV-1, hepatitis B virus and hepatitis C virus antiviral therapies have been successful at keeping viruses suppressed for prolonged periods of time, but therapy failures attributable to the emergence of drug resistant mutations continue to be a distressing reminder that no therapy can fully eradicate these viruses from their host organisms. To better understand the emergence of drug resistance, we combined phylogenetic and statistical models of viral evolution in a 2-phase computational approach that reconstructs mutational pathways of drug resistance.

RESULTS

The first phase of the algorithm involved the modeling of the evolution of the virus within the human host environment. The inclusion of longitudinal clonal sequence data was a key aspect of the model due to the progressive fashion in which multiple mutations become linked in the same genome creating drug resistant genotypes. The second phase involved the development of a Markov model to calculate the transition probabilities between the different genotypes. The proposed method was applied to data from an HIV-1 Efavirenz clinical trial study. The obtained model revealed the direction of evolution over time with greater detail than previous models. Our results show that the mutational pathways facilitate the identification of fast versus slow evolutionary pathways to drug resistance.

AVAILABILITY

Source code for the algorithm is publicly available at http://biorg.cis.fiu.edu/vPhyloMM/

Macrophage delivery of nanoformulated antiretroviral drug to the brain in a murine model of neuroAIDS

Antiretroviral therapy (ART) shows variable blood-brain barrier penetration. This may affect the development of neurological complications of HIV infection. In attempts to attenuate viral growth for the nervous system, cell-based nanoformulations were developed with the focus on improving drug pharmacokinetics. We reasoned that ART carriage could be facilitated within blood-borne macrophages traveling across the blood-brain barrier. To test this idea, an HIV-1 encephalitis (HIVE) rodent model was used where HIV-1-infected human monocyte-derived macrophages were stereotactically injected into the subcortex of severe combined immunodeficient mice. ART was prepared using indinavir (IDV) nanoparticles (NP, nanoART) loaded into murine bone marrow macrophages (BMM, IDV-NP-BMM) after ex vivo cultivation. IDV-NP-BMM was administered i.v. to mice resulting in continuous IDV release for 14 days. Rhodamine-labeled IDV-NP was readily observed in areas of HIVE and specifically in brain subregions with active astrogliosis, microgliosis, and neuronal loss. IDV-NP-BMM treatment led to robust IDV levels and reduced HIV-1 replication in HIVE brain regions. We conclude that nanoART targeting to diseased brain through macrophage carriage is possible and can be considered in developmental therapeutics for HIV-associated neurological disease.

Drug-associated changes in amino acid residues in Gag p2, p7(NC), and p6(Gag)/p6(Pol) in human immunodeficiency virus type 1 (HIV-1) display a dominant effect on replicative fitness and drug response

Regions of HIV-1 gag between p2 and p6(Gag)/p6(Pol), in addition to protease (PR), develop genetic diversity in HIV-1 infected individuals who fail to suppress virus replication by combination protease inhibitor (PI) therapy. To elucidate functional consequences for viral replication and PI susceptibility by changes in Gag that evolve in vivo during PI therapy, a panel of recombinant viruses was constructed. Residues in Gag p2/p7(NC) cleavage site and p7(NC), combined with residues in the flap of PR, defined novel fitness determinants that restored replicative capacity to the posttherapy virus. Multiple determinants in Gag have a dominant effect on PR phenotype and increase susceptibility to inhibitors of drug-resistant or drug-sensitive PR genes. Gag determinants of drug sensitivity and replication alter the fitness landscape of the virus, and viral replicative capacity can be independent of drug sensitivity. The functional linkage between Gag and PR provides targets for novel therapeutics to inhibit drug-resistant viruses.

Plasma HIV RNA decline and emergence of drug resistance mutations among patients with multiple virologic failures receiving resistance testing-guided HAART

Early recognition of virologic failure in patients with extensive drug resistance receiving salvage-HAART is essential to avoid exposure to subinhibitory regimens. We studied plasma viral load (PVL) decline and rates of drug-resistance mutation (DRM) accumulation in such patients. A prospective, 48 week study of 38 heavily pretreated patients receiving genotypic resistance testing (GRT)-guided HAART was conducted. The rate of PVL decline was studied by weekly PVL determinations. To assess DRM accumulation, serial GRTs were performed in all nonresponders (never reaching PVL <50 or two PVLs >50 copies/ml after suppression). Over 48 weeks, 10 patients (26%) were nonresponders. Receiving less then two fully active drugs and having an elevated number of PI and NRTI mutations at baseline were strongly associated with virologic failure. There was no evidence of a difference in the change from baseline PVL to week 1 and 2 between responders and nonresponders. By contrast, PVL reductions from week 2 to week 3 and thereafter were significantly greater for responders (p < 0.01). Among nonresponders, the incidence rates per patient-month (95% CI) of emergent DRM were 0.67 (0.13-1.20), 0.40 (0.00-0.74), and 0.37 (0.00-0.75) at weeks 4, 8, and 24, respectively. Having limited baseline resistance, receiving at least two fully active drugs, and showing constant PVL reductions from week 2 to week 3 and thereafter were predictive of virologic response. In contrast, early changes in PVL levels were not. Virologic failure was associated with detection of emergent DRMs. Virologic rebound in patients on salvage-HAART should be addressed aggressively.

Therapeutic vaccination reduces HIV sequence variability

With HIV persisting lifelong in infected persons, therapeutic vaccination is a novel alternative concept to control virus replication. Even though CD8 and CD4 cell responses to such immunizations have been demonstrated, their effects on virus replication are still unclear. In view of this fact, we studied the impact of a therapeutic vaccination with HIV nef delivered by a recombinant modified vaccinia Ankara vector on viral diversity. We investigated HIV sequences derived from chronically infected persons before and after therapeutic vaccination. Before immunization the mean +/- se pairwise variability of patient-derived Nef protein sequences was 0.1527 +/- 0.0041. After vaccination the respective value was 0.1249 +/- 0.0042, resulting in a significant (P<0.0001) difference between the two time points. The genes vif and 5'gag tested in parallel and nef sequences in control persons yielded a constant amino acid sequence variation. The data presented suggest that Nef immunization induced a selective pressure, limiting HIV sequence variability. To our knowledge this is the first report directly linking therapeutic HIV vaccination to decreasing diversity in patient-derived virus isolates.

Fitness landscape of human immunodeficiency virus type 1 protease quasispecies

Here we show, at a high resolution (1%), the human immunodeficiency virus type 1 (HIV-1) protease gene quasispecies landscape from three infected naïve individuals. A huge range of genetic configurations was found (67%, 71%, and 80% of the nucleotide clones from the three individuals, respectively, were different), and these configurations created a dense net that linked different parts of the viral population. Similarly, a vast diversity of different protease activities was also found. Importantly, 65% of the analyzed enzymes had detectable protease activity, and 11% of the minority individual variants showed similar or better fitness than the master (most abundant) enzyme, suggesting that the viral complexity in this genomic region does not exclusively depend on the enzyme's catalytic efficiency. Several high-fitness minority variants had only one substitution compared to the master sequence, supporting the possibility that the rugged HIV-1 protease quasispecies fitness landscape may be formed by a continuous network that can be traversed by single mutational steps without passing through defective or less-adapted proteins.

HIV-1 matrix protein: a mysterious regulator of the viral life cycle

Significant progress has been achieved in the last few years concerning the human immunodeficiency virus (HIV-1) life cycle, mostly in the fields of cellular receptors for the virus, virus assembly and budding of virus particles from the cell surface. Meanwhile, some aspects, such as postentry events, virus maturation and the regulatory role of individual viral proteins remain poorly defined. This review summarizes some recent findings concerning the role of Gag Pr55 and its proteolytic processing in the HIV-1 life cycle with particular emphasis on the functions of matrix protein p17 (MA), the protein which plays a key role in regulation of the early and late steps of viral morphogenesis. Based on our recent observations, the possibility is discussed that two subsets of MA exist, one cleaved from the Gag precursor in the host cell (cMA), and the other cleaved in the virions (vMA). It is suggested that two MA fractions possess diverse functions and are involved in different stages of virus morphogenesis as key regulators of the viral life cycle.

Altered gag polyprotein cleavage specificity of feline immunodeficiency virus/human immunodeficiency virus mutant proteases as demonstrated in a cell-based expression system

We have used feline immunodeficiency virus (FIV) protease (PR) as a mutational system to study the molecular basis of substrate-inhibitor specificity for lentivirus PRs, with a focus on human immunodeficiency virus type 1 (HIV-1) PR. Our previous mutagenesis studies demonstrated that discrete substitutions in the active site of FIV PR with structurally equivalent residues of HIV-1 PR dramatically altered the specificity of the mutant PRs in in vitro analyses. Here, we have expanded these studies to analyze the specificity changes in each mutant FIV PR expressed in the context of the natural Gag-Pol polyprotein ex vivo. Expression mutants were prepared in which 4 to 12 HIV-1-equivalent substitutions were made in FIV PR, and cleavage of each Gag-Pol polyprotein was then assessed in pseudovirions from transduced cells. The findings demonstrated that, as with in vitro analyses, inhibitor specificities of the mutants showed increased HIV-1 PR character when analyzed against the natural substrate. In addition, all of the mutant PRs still processed the FIV polyprotein but the apparent order of processing was altered relative to that observed with wild-type FIV PR. Given the importance of the order in which Gag-Pol is processed, these findings likely explain the failure to produce infectious FIVs bearing these mutations.

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

Treatment of human immunodeficiency virus type 1 with protease inhibitors (PIs) is associated with the emergence of resistance-associated mutations. Treatment-characterized datasets have been used to identify novel treatment-associated protease mutations. In this study, we utilized two large reference laboratory databases (>115,000 viral sequences) to identify non-established resistance-associated protease mutations. We found 20 non-established protease mutations occurring in 82% of viruses with a PI resistance score of 4-7, 62% of viruses with a resistance score of 1-3, and 35% of viruses with no predicted PI resistance. We correlated mutational prevalence to treatment duration in a treatment-characterized dataset of 2161 patients undergoing non-suppressive PI therapy. In the non-suppressed dataset, 24 mutations became more prevalent and three mutations became less prevalent after more than 48 months of non-suppressive PI-therapy. Longer durations of non-suppressive treatment correlated with higher PI resistance scores. Mutations at eight non-established positions that were more common in viruses with the longest duration of non-suppressive therapy were also more common in viruses with the highest PI resistance score. Covariation analysis of 3036 protease amino acid substitutions identified 75 positive and nine negative correlations between resistance associated positions. Our findings support the utility of reference laboratory datasets for surveillance of mutation prevalence and covariation.

Effect of a protease inhibitor-induced genetic bottleneck on human immunodeficiency virus type 1 env gene populations

The initiation of drug therapy or the addition of a new drug to preexisting therapy can have a significant impact on human immunodeficiency virus type 1 (HIV-1) populations within a person. Drug therapy directed at reverse transcriptase and protease can result in dramatic decreases in virus load, causing a contraction in the virus population that represents a potential genetic bottleneck as a subset of virus with genomes carrying resistance mutations repopulate the host. While this bottleneck exerts an effect directly on the region that is being targeted by the drugs, it also affects other regions of the viral genome. We have applied heteroduplex tracking assays (HTA) specific to variable regions 1 and 2 (V1/V2) and variable region 3 (V3) of the HIV-1 env gene to analyze the effect of a genetic bottleneck created by the selection of resistance to ritonavir, a protease inhibitor. Subjects were classified into groups on the basis of the extent of the initial drop in virus load and the duration of virus load reduction. Subjects with a strong initial drop in virus load exhibited a loss of heterogeneity in the env region at virus load rebound; in contrast, subjects with a weak initial drop in virus load exhibited little to no loss of heterogeneity at virus load rebound in either region of env examined. The duration of virus load reduction also affected env populations. Subjects that had prolonged reductions exhibited slower population diversification and the appearance of new V1/V2 species after rebound. The longer reduction of virus load in these subjects may have allowed for improved immune system function, which in turn could have selected for new escape mutants. Subjects with rapid rebound quickly reequilibrated the entry env variants back into the resistant population. When the pro gene developed further resistance mutations subsequent to virus load rebound, no changes were observed in V1/V2 or V3 populations, suggesting that the high virus loads allowed the env populations to reequilibrate rapidly. The rapid equilibration of env variants during pro gene sequence transitions at high virus load suggests that recombination is active in defining the HIV-1 sequence population. Conversely, part of the success of suppressive antiviral therapy may be to limit the potential for evolution through recombination, which requires dually infected cells.