Exploring the Complexity of the HIV-1 Fitness Landscape

Although fitness landscapes are central to evolutionary theory, so far no biologically realistic examples for large-scale fitness landscapes have been described. Most currently available biological examples are restricted to very few loci or alleles and therefore do not capture the high dimensionality characteristic of real fitness landscapes. Here we analyze large-scale fitness landscapes that are based on predictive models for in vitro replicative fitness of HIV-1. We find that these landscapes are characterized by large correlation lengths, considerable neutrality, and high ruggedness and that these properties depend only weakly on whether fitness is measured in the absence or presence of different antiretrovirals. Accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. While the relative extent to which mutations affect fitness on their own (main effects) or in combination with other mutations (epistasis) is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged, less neutral, and more correlated than expected from the distribution of main effects and epistatic interactions alone. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place

Assessing Predicted HIV-1 Replicative Capacity in a Clinical Setting

HIV-1 replicative capacity (RC) provides a measure of within-host fitness and is determined in the context of phenotypic drug resistance testing. However it is unclear how these in-vitro measurements relate to in-vivo processes. Here we assess RCs in a clinical setting by combining a previously published machine-learning tool, which predicts RC values from partial pol sequences with genotypic and clinical data from the Swiss HIV Cohort Study. The machine-learning tool is based on a training set consisting of 65000 RC measurements paired with their corresponding partial pol sequences. We find that predicted RC values (pRCs) correlate significantly with the virus load measured in 2073 infected but drug naïve individuals. Furthermore, we find that, for 53 pairs of sequences, each pair sampled in the same infected individual, the pRC was significantly higher for the sequence sampled later in the infection and that the increase in pRC was also significantly correlated with the increase in plasma viral load and with the length of the time-interval between the sampling points. These findings indicate that selection within a patient favors the evolution of higher replicative capacities and that these in-vitro fitness measures are indicative of in-vivo HIV virus load

Estimating the Fitness Cost of Escape from HLA Presentation in HIV-1 Protease and Reverse Transcriptase

Human immunodeficiency virus (HIV-1) is, like most pathogens, under selective pressure to escape the immune system of its host. In particular, HIV-1 can avoid recognition by cytotoxic T lymphocytes (CTLs) by altering the binding affinity of viral peptides to human leukocyte antigen (HLA) molecules, the role of which is to present those peptides to the immune system. It is generally assumed that HLA escape mutations carry a replicative fitness cost, but these costs have not been quantified. In this study, we assess the replicative cost of mutations which are likely to escape presentation by HLA molecules in the region of HIV-1 protease and reverse transcriptase. Specifically, we combine computational approaches for prediction of in vitro replicative fitness and peptide binding affinity to HLA molecules. We find that mutations which impair binding to HLA-A molecules tend to have lower in vitro replicative fitness than mutations which do not impair binding to HLA-A molecules, suggesting that HLA-A escape mutations carry higher fitness costs than non-escape mutations. We argue that the association between fitness and HLA-A binding impairment is probably due to an intrinsic cost of escape from HLA-A molecules, and these costs are particularly strong for HLA-A alleles associated with efficient virus control. Counter-intuitively, we do not observe a significant effect in the case of HLA-B, but, as discussed, this does not argue against the relevance of HLA-B in virus control. Overall, this article points to the intriguing possibility that HLA-A molecules preferentially target more conserved regions of HIV-1, emphasizing the importance of HLA-A genes in the evolution of HIV-1 and RNA viruses in general

Amino Acid Substitutions at Position 190 of Human Immunodeficiency Virus Type 1 Reverse Transcriptase Increase Susceptibility to Delavirdine and Impair Virus Replication

Suboptimal treatment of human immunodeficiency virus type 1 (HIV-1) infection with nonnucleoside reverse transcriptase inhibitors (NNRTI) often results in the rapid selection of drug-resistant virus. Several amino acid substitutions at position 190 of reverse transcriptase (RT) have been associated with reduced susceptibility to the NNRTI, especially nevirapine (NVP) and efavirenz (EFV). In the present study, the effects of various 190 substitutions observed in viruses obtained from NNRTI-experienced patients were characterized with patient-derived HIV isolates and confirmed with a panel of isogenic viruses. Compared to wild-type HIV, which has a glycine at position 190 (G190), viruses with 190 substitutions (A, C, Q, S, V, E, or T, collectively referred to as G190X substitutions) were markedly less susceptible to NVP and EFV. In contrast, delavirdine (DLV) susceptibility of these G190X viruses increased from 3 to 300-fold (hypersusceptible) or was only slightly decreased. The replication capacity of viruses with certain 190 substitutions (C, Q, V, T, and E) was severely impaired and was correlated with reduced virion-associated RT activity and incomplete protease (PR) processing of the viral p55(gag) polyprotein. These defects were the result of inadequate p160(gagpol) incorporation into virions. Compensatory mutations within RT and PR improved replication capacity, p55(gag) processing, and RT activity, presumably through increased incorporation of p160(gagpol) into virions. We observe an inverse relationship between the degree of NVP and EFV resistance and the impairment of viral replication in viruses with substitutions at 190 in RT. These observations may have important implications for the future design and development of antiretroviral drugs that restrict the outgrowth of resistant variants with high replication capacity

Suppression of Dualtropic Human Immunodeficiency Virus Type 1 by the CXCR4 Antagonist AMD3100 Is Associated with Efficiency of CXCR4 Use and Baseline Virus Composition▿

In a phase I/II evaluation of the CXCR4 antagonist AMD3100, human immunodeficiency virus RNA levels were significantly reduced in a single study subject who harbored CXCR4 (X4)-tropic virus, but not in subjects who harbored either dual/mixed (DM)-tropic or CCR5 (R5)-tropic virus (C. W. Hendrix et al., J. Acquir. Immune Defic. Syndr. 37:1253-1262, 2004). In this study, we analyzed the envelope clones of DM-tropic virus in baseline and treated virus populations from 14 subjects. Ten subjects exhibited significant reductions in CXCR4-mediated infectivity after 10 days of AMD3100 therapy relative to baseline (X4 suppressor group), while four subjects had no reduction of CXCR4-mediated infectivity (X4 nonsuppressor group). The baseline viruses of the X4 suppressor group infected CXCR4-expressing cells less efficiently than those of the X4 nonsuppressor group. Clonal analysis indicated that the baseline viruses from the X4 suppressor group contained a higher proportion of R5-tropic variants mixed with CXCR4-using variants, while the X4 nonsuppressor group was enriched for CXCR4-using variants. AMD3100 suppressed X4-tropic variants in all subjects studied, but not all dualtropic variants. Furthermore, dualtropic variants that used CXCR4 efficiently were suppressed by AMD3100, while dualtropic variants that used CXCR4 poorly were not. This study demonstrated that AMD3100 has the ability to suppress both X4-tropic and certain dualtropic variants in vivo. The suppression of CXCR4-using variants by AMD3100 is dependent on both the tropism composition of the virus population and the efficiency of CXCR4 usage of individual variants

Coreceptor Tropism Can Be Influenced by Amino Acid Substitutions in the gp41 Transmembrane Subunit of Human Immunodeficiency Virus Type 1 Envelope Protein▿

Many studies have demonstrated that the third variable region (V3) of the human immunodeficiency virus type 1 (HIV-1) envelope protein (Env) is a major determinant of coreceptor tropism. Other regions in the surface gp120 subunit of Env can modulate coreceptor tropism in a manner that is not fully understood. In this study, we evaluated the effect of env determinants outside of V3 on coreceptor usage through the analysis of (i) patient-derived env clones that differ in coreceptor tropism, (ii) chimeric env sequences, and (iii) site-directed mutants. The introduction of distinct V3 sequences from CXCR4-using clones into an R5-tropic env backbone conferred the inefficient use of CXCR4 in some but not all cases. Conversely, in many cases, X4- and dual-tropic env backbones containing the V3 sequences of R5-tropic clones retained the ability to use CXCR4, suggesting that sequences outside of the V3 regions of these CXCR4-using clones were responsible for CXCR4 use. The determinants of CXCR4 use in a set of dual-tropic env sequences with V3 sequences identical to those of R5-tropic clones mapped to the gp41 transmembrane (TM) subunit. In one case, a single-amino-acid substitution in the fusion peptide of TM was able to confer CXCR4 use; however, TM substitutions associated with CXCR4 use varied among different env sequences. These results demonstrate that sequences in TM can modulate coreceptor specificity and that env sequences other than that of V3 may facilitate efficient CXCR4-mediated entry. We hypothesize that the latter plays an important role in the transition from CCR5 to CXCR4 coreceptor use

Emergence of CXCR4-Using Human Immunodeficiency Virus Type 1 (HIV-1) Variants in a Minority of HIV-1-Infected Patients following Treatment with the CCR5 Antagonist Maraviroc Is from a Pretreatment CXCR4-Using Virus Reservoir

Antagonists of the human immunodeficiency virus type 1 (HIV-1) coreceptor, CCR5, are being developed as the first anti-HIV agents acting on a host cell target. We monitored the coreceptor tropism of circulating virus, screened at baseline for coreceptor tropism, in 64 HIV-1-infected patients who received maraviroc (MVC, UK-427,857) as monotherapy for 10 days. Sixty-two patients harbored CCR5-tropic virus at baseline and had a posttreatment phenotype result. Circulating virus remained CCR5 tropic in 60/62 patients, 51 of whom experienced an HIV RNA reduction from baseline of >1 log(10) copies/ml, indicating that CXCR4-using variants were not rapidly selected despite CCR5-specific drug pressure. In two patients, viral load declined during treatment and CXCR4-using virus was detected at day 11. No pretreatment factor predicted the emergence of CXCR4-tropic virus during maraviroc therapy in these two patients. Phylogenetic analysis of envelope (Env) clones from pre- and posttreatment time points indicated that the CXCR4-using variants probably emerged by outgrowth of a pretreatment CXCR4-using reservoir, rather than via coreceptor switch of a CCR5-tropic clone under selection pressure from maraviroc. Phylogenetic analysis was also performed on Env clones from a third patient harboring CXCR4-using virus prior to treatment. This patient was enrolled due to a sample labeling error. Although this patient experienced no overall reduction in viral load in response to treatment, the CCR5-tropic components of the circulating virus did appear to be suppressed while receiving maraviroc as monotherapy. Importantly, in all three patients, circulating virus reverted to predominantly CCR5 tropic following cessation of maraviroc

Crystal structure of HIV-1 reverse transcriptase in complex with a polypurine tract RNA:DNA

We have determined the 3.0 Å resolution structure of wild-type HIV-1 reverse transcriptase in complex with an RNA:DNA oligonucleotide whose sequence includes a purine-rich segment from the HIV-1 genome called the polypurine tract (PPT). The PPT is resistant to ribonuclease H (RNase H) cleavage and is used as a primer for second DNA strand synthesis. The ‘RNase H primer grip’, consisting of amino acids that interact with the DNA primer strand, may contribute to RNase H catalysis and cleavage specificity. Cleavage specificity is also controlled by the width of the minor groove and the trajectory of the RNA:DNA, both of which are sequence dependent. An unusual ‘unzipping’ of 7 bp occurs in the adenine stretch of the PPT: an unpaired base on the template strand takes the base pairing out of register and then, following two offset base pairs, an unpaired base on the primer strand re-establishes the normal register. The structural aberration extends to the RNase H active site and may play a role in the resistance of PPT to RNase H cleavage