Epistasis as a Determinant of the HIV-1 Protease's Robustness to Mutation

The robustness of phenotypes to mutation is critical to protein evolution; robustness may be an adaptive trait if it promotes evolution. We hypothesised that native proteins subjected to natural selection in vivo should be more robust than proteins generated in vitro in the absence of natural selection. We compared the mutational robustness of two human immunodeficiency virus type 1 (HIV-1) proteases with comparable catalytic efficiencies, one isolated from an infected individual and the second generated in vitro via random mutagenesis. Single mutations in the protease (82 and 60 in the wild-type and mutant backgrounds, respectively) were randomly generated in vitro and the catalytic efficiency of each mutant was determined. No differences were observed between these two protease variants when lethal, neutral, and deleterious mutations were compared (P = 0.8025, chi-squared test). Similarly, average catalytic efficiency (−72.6% and −64.5%, respectively) did not significantly differ between protease mutant libraries (P = 0.3414, Mann Whitney test). Overall, the two parental proteins displayed similar mutational robustness. Importantly, strong and widespread epistatic interactions were observed when the effect of the same mutation was compared in both proteases, suggesting that epistasis can be a key determinant of the robustness displayed by the in vitro generated protease.     

Epistasis among Deleterious Mutations in the HIV-1 Protease

A central goal in molecular evolution is to understand how genetic interactions between protein mutations shape protein function and fitness. While intergenic epistasis has been extensively explored in eukaryotes, bacteria, and viruses, intragenic epistatic interactions have been insufficiently studied. Here, we employ a model system in which lambda phage fitness correlates with the enzymatic activity of human immunodeficiency virus type 1 (HIV-1) protease to systematically determine the epistatic interactions between intragenic pairs of deleterious protein substitutions. We generated 114 genotypes of the HIV-1 protease, each carrying pairs of nucleotide substitution mutations whose separated and combined deleterious effects on fitness were then determined. A high proportion (39%) of pairs displayed lethality. Several pairs exhibited significant interactions for fitness, including positive and negative epistasis. Significant negative epistatic interactions predominated (15%) over positive interactions (2%). However, the average ± SD epistatic effect, ē = 0.0025 ± 0.1334, was not significantly different from zero (p = 0.8368). Notably, epistatic interactions, regardless of epistatic direction, tend to be more frequent in the context of less deleterious mutations. In the present study, the high frequencies of lethality and negative epistasis indicate that the HIV-1 protease is highly sensitive to the effects of deleterious mutations. Therefore, proteins may not be as robust to mutational change as is usually expected.     

Small fitness effects and weak genetic interactions between deleterious mutations in heterozygous loci of the yeast Saccharomyces cerevisiae

Rare, random mutations were induced in budding yeast by ethyl methanesulfonate (EMS). Clones known to bear a single non-neutral mutation were used to obtain mutant heterozygotes and mutant homozygotes that were later compared with wild-type homozygotes. The average homozygous effect of mutation was an approximately 2% decrease in the growth rate. In heterozygotes, the harmful effect of these relatively mild mutations was reduced approximately fivefold. In a test of epistasis, two heterozygous mutant loci were paired at random. Fitness of the double mutants was best explained by multiplicative action of effects at single loci, with little evidence for epistasis and essentially excluding synergism. In other experiments, the same mutations in haploid and heterozygous diploid clones were compared. Regardless of the haploid phenotypes, mildly deleterious or lethal, fitness of the heterozygotes was decreased by less than half a per cent on average. In general, the results presented here suggest that most mutations tend to exhibit small and weakly interacting effects in heterozygous loci regardless of how harmful they are in haploids or homozygotes.     

Muller''s Ratchet, Epistasis and Mutation Effects

In this study, computer simulation is used to show that despite synergistic epistasis for fitness, Muller's ratchet can lead to lethal fitness loss in a population of asexuals through the accumulation of deleterious mutations. This result contradicts previous work that indicated that epistasis will halt the ratchet. The present results show that epistasis will not halt the ratchet provided that rather than a single deleterious mutation effect, there is a distribution of deleterious mutation effects with sufficient density near zero. In addition to epistasis and mutation distribution, the ability of Muller's ratchet to lead to the extinction of an asexual population under epistasis for fitness depends strongly on the expected number of offspring that survive to reproductive age. This strong dependence is not present in the nonepistatic model and suggests that interpreting the population growth parameter as fecundity is inadequate. Because a continuous distribution of mutation effects is used in this model, an emphasis is placed on the dynamics of the mutation effect distribution rather than on the dynamics of the number of least mutation loaded individuals. This perspective suggests that current models of gene interaction are too simple to apply directly to long-term prediction for populations undergoing the ratchet.     

Epistasis and the adaptability of an RNA virus

We have explored the patterns of fitness recovery in the vesicular stomatitis RNA virus. We show that, in our experimental setting, reversions to the wild-type genotype were rare and fitness recovery was at least partially driven by compensatory mutations. We compared compensatory adaptation for genotypes carrying (1) mutations with varying deleterious fitness effects, (2) one or two deleterious mutations, and (3) pairs of mutations showing differences in the strength and sign of epistasis. In all cases, we found that the rate of fitness recovery and the proportion of reversions were positively affected by population size. Additionally, we observed that mutations with large fitness effect were always compensated faster than mutations with small fitness effect. Similarly, compensatory evolution was faster for genotypes carrying a single deleterious mutation than for those carrying pairs of mutations. Finally, for genotypes carrying two deleterious mutations, we found evidence of a negative correlation between the epistastic effect and the rate of compensatory evolution.     

Fitness ranking of individual mutants drives patterns of epistatic interactions in HIV-1

Fitness interactions between mutations, referred to as epistasis, can strongly impact evolution. For RNA viruses and retroviruses with their high mutation rates, epistasis may be particularly important to overcome fitness losses due to the accumulation of deleterious mutations and thus could influence the frequency of mutants in a viral population. As human immunodeficiency virus type 1 (HIV-1) resistance to azidothymidine (AZT) requires selection of sequential mutations, it is a good system to study the impact of epistasis. Here we present a thorough analysis of a classical AZT-resistance pathway (the 41-215 cluster) of HIV-1 variants by fitness measurements in single round infection assays covering physiological drug concentrations ex vivo. The sign and value of epistasis varied and did not predict the epistatic effect on the mutant frequency. This complex behavior is explained by the fitness ranking of the variants that strongly depends on environmental factors, i.e., the presence and absence of drugs and the host cells used. Although some interactions compensate fitness losses, the observed small effect on the relative mutant frequencies suggests that epistasis might be inefficient as a buffering mechanism for fitness losses in vivo. While the use of epistasis-based hypotheses to make general assumptions on the evolutionary dynamics of viral populations is appealing, our data caution their interpretation without further knowledge on the characteristics of the viral mutant spectrum under different environmental conditions.     

Secondary mutations M36I and A71V in the human immunodeficiency virus type 1 protease can provide an advantage for the emergence of the primary mutation D30N

Development of resistance mutations in enzymatic targets of human immunodeficiency virus 1 (HIV-1) hampers the ability to provide adequate therapy. Of special interest is the effect mutations outside the active site of HIV-1 protease have on inhibitor binding and virus viability. We engineered protease mutants containing the active site mutation D30N alone and with the nonactive site polymorphisms M36I and/or A71V. We determined the K(i) values for the inhibitors nelfinavir, ritonavir, indinavir, KNI272, and AG1776 as well as the catalytic efficiency of the mutants. Single and double mutation combinations exhibited a decrease in catalytic efficiency, while the triple mutant displayed catalytic efficiency greater than that of the wild type. Variants containing M36I or A71V alone did not display a significant change in binding affinities to the inhibitors tested. The variant containing mutation D30N displayed a 2-6-fold increase in K(i) for all inhibitors tested, with nelfinavir showing the greatest increase. The double mutants containing a combination of mutations D30N, M36I, and A71V displayed -0.5-fold to +6-fold changes in the K(i) of all inhibitors tested, with ritonavir and nelfinavir most affected. Only the triple mutant showed a significant increase (>10-fold) in K(i) for inhibitor nelfinavir, ritonavir, or AG-1776 displaying 22-, 19-, or 15-fold increases, respectively. Our study shows that the M36I and A71V mutations provide a greater level of inhibitor cross-resistance combined with active site mutation D30N. M36I and A71V, when present as natural polymorphisms, could aid the virus in developing active site mutations to escape inhibitor binding while maintaining catalytic efficiency.     

Impact of Mutations in Highly Conserved Amino Acids of the HIV-1 Gag-p24 and Env-gp120 Proteins on Viral Replication in Different Genetic Backgrounds

It has been hypothesized that a single mutation at a highly conserved amino acid site (HCS) can be severely deleterious to HIV in most if not all isolate-specific genetic backgrounds. Consequently, potentially universal HIV-1 vaccines exclusively targeting highly conserved regions of the viral proteome have been proposed. To test this hypothesis, we examined the impact of 10 Gag-p24 and 9 Env-gp120 HCS single mutations on viral fitness. In the original founder sequence of the subject in whom these mutations were identified, all Gag-p24 HCS mutations significantly reduced viral replication fitness, including 7 that were lethal. Similar results were obtained at 9/10 sites when the same mutations were introduced into the founder sequences of two epidemiologically unlinked subjects. In contrast, none of the 9 Env-gp120 HCS mutations were lethal in the original founder sequence, and four had no fitness cost. Hence, HCS mutations in Gag-p24 are likely to be severely deleterious in different HIV-1 subtype B backgrounds; however, some HCS mutations in both Gag-p24 and Env-gp120 fragments can be well tolerated. Therefore, when designing HIV-1 immunogens that are intended to force the virus to nonviable escape pathways, the fitness constraints on the HIV segments included should be considered beyond their conservation level.     

Rapid evolutionary escape by large populations from local fitness peaks is likely in nature

Fitness interactions between loci in the genome, or epistasis, can result in mutations that are individually deleterious but jointly beneficial. Such epistasis gives rise to multiple peaks on the genotypic fitness landscape. The problem of evolutionary escape from such local peaks has been a central problem of evolutionary genetics for at least 75 years. Much attention has focused on models of small populations, in which the sequential fixation of valley genotypes carrying individually deleterious mutations operates most quickly owing to genetic drift. However, valley genotypes can also be subject to mutation while transiently segregating, giving rise to copies of the high fitness escape genotype carrying the jointly beneficial mutations. In the absence of genetic recombination, these mutations may then fix simultaneously. The time for this process declines sharply with increasing population size, and it eventually comes to dominate evolutionary behavior. Here we develop an analytic expression for N(crit), the critical population size that defines the boundary between these regimes, which shows that both are likely to operate in nature. Frequent recombination may disrupt high-fitness escape genotypes produced in populations larger than N(crit) before they reach fixation, defining a third regime whose rate again slows with increasing population size. We develop a novel expression for this critical recombination rate, which shows that in large populations the simultaneous fixation of mutations that are beneficial only jointly is unlikely to be disrupted by genetic recombination if their map distance is on the order of the size of single genes. Thus, counterintuitively, mass selection alone offers a biologically realistic resolution to the problem of evolutionary escape from local fitness peaks in natural populations.     

The role of deleterious mutations in allopatric speciation

Allopatric speciation is often assumed to occur as a consequence of adaptive divergence between two isolated populations. However, there are some scenarios in which reproductive isolation can be favored due to accumulated unconditionally deleterious mutations. If deleterious mutations have synergistic epistatic effects, it is shown here that the average fitness of recombinants between two parental lines with a given number of fixed mutations is lower than that of the parents in both the F1 and F2 generations. If individual mutations are only slightly deleterious, then they will tend to fixation at a high enough rate to cause lower hybrid fitness. If the fitness effects of mutation give rise to antagonistic epistasis, the hybrids tend to have a higher average fitness than the parental lines, suggesting a possible scenario for the origin of hybrid vigor. The other model of deleterious mutations investigated is the accumulation of knockout mutants in a duplicated gene family. While neutral in the parental lines, upon contact the F1 and later generations have a significant probability of carrying double knockouts. Under this scenario, selection may also favor reproductive isolation between the two lines. Even when the selection coefficients generated are too low to drive speciation, epistatic interactions between deleterious mutations offer a possible explanation for both outbreeding depression and hybrid vigor.     

FIXATION OF MUTATORS IN ASEXUAL POPULATIONS: THE ROLE OF GENETIC DRIFT AND EPISTASIS

We study the evolutionary dynamics of an asexual population of nonmutators and mutators on a class of epistatic fitness landscapes. We consider the situation in which all mutations are deleterious and mutators are produced from nonmutators continually at a constant rate. We find that in an infinitely large population, a minimum nonmutator‐to‐mutator conversion rate is required to fix the mutators but an arbitrarily small conversion rate results in the fixation of mutators in a finite population. We calculate analytical expressions for the mutator fraction at mutation‐selection balance and fixation time for mutators in a finite population when the difference between the mutation rate for mutator and nonmutator is smaller (regime I) and larger (regime II) than the selection coefficient. Our main result is that in regime I, the mutator fraction and the fixation time are independent of epistasis but in regime II, mutators are rarer and take longer to fix when the decrease in fitness with the number of deleterious mutations occurs at an accelerating rate (synergistic epistasis) than at a diminishing rate (antagonistic epistasis). Our analytical results are compared with numerics and their implications are discussed.     

Factors affecting the genetic load in Drosophila: synergistic epistasis and correlations among fitness components

Two factors that can affect genetic load, synergistic epistasis and sexual selection, were investigated in Drosophila melanogaster. A set of five chromosomal regions containing visible recessive mutations were put together in all combinations to create a full set of 32 homozygous lines fixed for different numbers of known mutations. Two measures of fitness were made for each line: productivity (a combined measure of fecundity and egg-to-adult survivorship) and competitive male mating success. Productivity, but not male mating success, showed a pattern of strong average synergistic epistasis, such that the log fitness declined nonlinearly with increasing numbers of mutations. Synergistic epistasis is known to reduce the mutation load. Both fitness components show some positive and some negative interactions between specific sets of mutations. Furthermore, alleles with deleterious effects on productivity tend to also diminish male mating success. Given that male mating success can affect relative fitness without changing the mean productivity of a population, these additional effects would lead to lower frequencies and lower fixation rates of deleterious alleles without higher costs to the mean fitness of the population.     

Pervasive Pairwise Intragenic Epistasis among Sequential Mutations in TEM-1 β-Lactamase

Interactions between mutations play a central role in shaping the fitness landscape, but a clear picture of intragenic epistasis has yet to emerge. To further reveal the prevalence and patterns of intragenic epistasis, we present a survey of epistatic interactions between sequential mutations in TEM-1 β-lactamase. We measured the fitness effect of ~ 12,000 pairs of consecutive amino acid substitutions and used our previous study of the fitness effects of single amino acid substitutions to calculate epistasis for over 8000 mutation pairs. Since sequential mutations are prone to physically interact, we postulated that our study would be surveying specific epistasis instead of nonspecific epistasis. We found widespread negative epistasis, especially in beta-strands, and a high frequency of negative sign epistasis among individually beneficial mutations. Negative epistasis (52%) occurred 7.6 times as frequently as positive epistasis (6.8%). Buried residues experienced more negative epistasis that surface-exposed residues. However, TEM-1 exhibited a couple of hotspots for positive epistasis, most notably L221/ R222 at which many combinations of mutations positively interacted. This study is the first to systematically examine pairwise epistasis throughout an entire protein performing its native function in its native host.     

Epistasis and its relationship to canalization in the RNA virus phi 6

Although deleterious mutations are believed to play a critical role in evolution, assessing their realized effect has been difficult. A key parameter governing the effect of deleterious mutations is the nature of epistasis, the interaction between the mutations. RNA viruses should provide one of the best systems for investigating the nature of epistasis because the high mutation rate allows a thorough investigation of mutational effects and interactions. Nonetheless, previous investigations of RNA viruses by S. Crotty and co-workers and by S. F. Elena have been unable to detect a significant effect of epistasis. Here we provide evidence that positive epistasis is characteristic of deleterious mutations in the RNA bacteriophage phi 6. We estimated the effects of deleterious mutations by performing mutation-accumulation experiments on five viral genotypes of decreasing fitness. We inferred positive epistasis because viral genotypes with low fitness were found to be less sensitive to deleterious mutations. We further examined environmental sensitivity in these genotypes and found that low-fitness genotypes were also less sensitive to environmental perturbations. Our results suggest that even random mutations impact the degree of canalization, the buffering of a phenotype against genetic and environmental perturbations. In addition, our results suggest that genetic and environmental canalization have the same developmental basis and finally that an understanding of the nature of epistasis may first require an understanding of the nature of canalization.     

Impaired fitness of human immunodeficiency virus type 1 variants with high-level resistance to protease inhibitors.

One hope to maintain the benefits of antiviral therapy against the human immunodeficiency virus type 1 (HIV-1), despite the development of resistance, is the possibility that resistant variants will show decreased viral fitness. To study this possibility, HIV-1 variants showing high-level resistance (up to 1,500-fold) to the substrate analog protease inhibitors BILA 1906 BS and BILA 2185 BS have been characterized. Active-site mutations V32I and I84V/A were consistently observed in the protease of highly resistant viruses, along with up to six other mutations. In vitro studies with recombinant mutant proteases demonstrated that these mutations resulted in up to 10(4)-fold increases in the Ki values toward BILA 1906 BS and BILA 2185 BS and a concomitant 2,200-fold decrease in catalytic efficiency of the enzymes toward a synthetic substrate. When introduced into viral molecular clones, the protease mutations impaired polyprotein processing, consistent with a decrease in enzyme activity in virions. Despite these observations, however, most mutations had little effect on viral replication except when the active-site mutations V32I and I84V/A were coexpressed in the protease. The latter combinations not only conferred a significant growth reduction of viral clones on peripheral blood mononuclear cells but also caused the complete disappearance of mutated clones when cocultured with wild-type virus on T-cell lines. Furthermore, the double nucleotide mutation I84A rapidly reverted to I84V upon drug removal, confirming its impact on viral fitness. Therefore, high-level resistance to protease inhibitors can be associated with impaired viral fitness, suggesting that antiviral therapies with such inhibitors may maintain some clinical benefits.     

Quantitative SNP-detection method for estimating HIV-1 replicative fitness: application to protease inhibitor-resistant viruses

We have improved the methods for the standard competitive growth assay of human immunodeficiency virus type 1 (HIV-1). The cloning step for the mixed viral population and subsequent genotype analysis for arbitrary numbers of clones were excluded from procedures. Instead, a single nucleotide polymorphism (SNP)-detection step was devised for the determination of viral populations. The quantitative SNP-detection method can rapidly estimate the proportion of wild-type and mutant populations with high reproducibility. Consequently, this method allows manipulation of many samples within a short period. Using this new competitive growth assay, replicative fitness of drug-resistant HIV-1 containing an M46I amino acid mutation in the protease was assessed in the presence or absence of indinavir. Without indinavir, replicative fitness of wild-type HIV-1 surpassed that of M46I-mutated HIV-1, and the fraction of mutated virus was reduced to about 10% at passage #9. In contrast, the fraction of M46I-mutated virus increased to >90% at passage #5 in the presence of 26.4 nM indinavir. Almost identical results were obtained for L90M-mutated HIV-1 with or without saquinavir. HIV-1 can survive under indinavir pressure by acquiring M46I mutation, as with acquisition of the L90M mutation under saquinavir pressure. However, these mutations damage viral replicative fitness under natural conditions without any drugs. Subtle differences between wild-type and mutant viruses are thus easily detected using the improved method.     

Compensatory mutations are repeatable and clustered within proteins

Compensatory mutations improve fitness in genotypes that contain deleterious mutations but have no beneficial effects otherwise. As such, compensatory mutations represent a very specific form of epistasis. We show that intragenic compensatory mutations occur non-randomly over gene sequence. Compensatory mutations are more likely to appear at some sites than others. Moreover, the sites of compensatory mutations are more likely than expected by chance to be near the site of the original deleterious mutation. Furthermore, compensatory mutations tend to occur more commonly in certain regions of the protein even when controlling for clustering around the site of the deleterious mutation. These results suggest that compensatory evolution at the protein level is partially predictable and may be convergent.     

Mutational and structural studies aimed at characterizing the monomer of HIV-1 protease and its precursor

An experimental protocol for folding the mature human immunodeficiency virus-1 (HIV-1) protease is presented that facilitates NMR studies at a low protein concentration of approximately 20 micoM. Under these conditions, NMR spectra show that the mature protease lacking its terminal beta-sheet residues 1-4 and 96-99 (PR(5-95)) exhibits a stable monomer fold spanning the region 10-90 that is similar to that of the single subunit of the wild-type dimer and the dimer bearing a D25N mutation (PR(D25N)). Urea-induced unfolding monitored both by changes in (1)H-(15)N heteronuclear single quantum correlation spectra and by protein fluorescence indicates that although PR(5-95) monomer displays a transition profile similar to that of the PR(D25N) dimer (50% unfolded (U(50)) = approximately 1.9 M), extending the protease with 4 residues (SFNF) of its N-terminally flanking sequence in the Gag-Pol precursor ((SFNF)PR(D25N)) decreases the stability of the fold (U(50) = approximately 1.5 M). Assigned backbone chemical shifts were used to elucidate differences in the stability of the PR(T26A) (U(50) = 2.5 M) and (SFNF)PR(D25N) monomers and compared with PR(D25N/T26A) monomer. Discernible differences in the backbone chemical shifts were observed for N-terminal protease residues 3-6 of (SFNF)PR(D25N) that may relate to the increase in the equilibrium dissociation constant (K(d)) and the very low catalytic activity of the protease prior to its autoprocessing at its N terminus from the Gag-Pol precursor.     

A major role for a set of non-active site mutations in the development of HIV-1 protease drug resistance

A major problem in the chemotherapy of HIV-1 infection is the appearance of drug resistance. In the case of HIV-1 protease inhibitors, resistance originates from mutations in the protease molecule that lower the affinity of inhibitors while still maintaining a viable enzymatic profile. Drug resistance mutations can be classified as active site or non-active site mutations depending on their location within the protease molecule. Active site mutations directly affect drug/target interactions, and their action can be readily understood in structural terms. Non-active site mutations influence binding from distal locations, and their mechanism of action is not immediately apparent. In this paper, we have characterized a mutant form of the HIV-1 protease, ANAM-11, identified in clinical isolates from HIV-1 infected patients treated with protease inhibitors. This mutant protease contains 11 mutations, 10 of which are located outside the active site (L10I/M36I/S37D/M46I/R57K/L63P/A71V/G73S/L90M/I93L) and 1 within the active site (I84V). ANAM-11 lowers the binding affinity of indinavir, nelfinavir, saquinavir, and ritonavir by factors of 4000, 3300, 5800, and 80000, respectively. Surprisingly, most of the loss in inhibitor affinity is due to the non-active site mutations as demonstrated by additional experiments performed with a protease containing only the 10 non-active site mutations (NAM-10) and another containing only the active site mutation (A-1). Kinetic analysis with two different substrates yielded comparable catalytic efficiencies for A-1, ANAM-11, NAM-10, and the wild-type protease. These studies demonstrate that non-active site mutations can be the primary source of resistance and that their role is not necessarily limited to compensate deleterious effects of active site mutations. Analysis of the structural stability of the proteases by differential scanning calorimetry reveals that ANAM-11 and NAM-10 are structurally more stable than the wild-type protease while A-1 is less stable. Together, the binding and structural thermodynamic results suggest that the non-active site mutants affect inhibitor binding by altering the geometry of the binding site cavity through the accumulation of mutations within the core of the protease molecule.