MBE Advance Access originally published online on November 7, 2006
Molecular Biology and Evolution 2007 24(2):382-387; doi:10.1093/molbev/msl168
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
HIV-1 Protease Catalytic Efficiency Effects Caused by Random Single Amino Acid Substitutions
Fundació irsiCaixa, Universitat Autònoma de Barcelona, Barcelona, Spain
E-mail: mmartinez{at}irsicaixa.es.
 |
Abstract |
|---|
 TOP Abstract IntroductionMaterials and MethodsResultsDiscussionAcknowledgementsReferences |
|---|
Protein evolution has occurred by successive fixation of individual mutations. The probability of fixation depends on the fitness of the mutation, and the arising variant can be deleterious, neutral, or beneficial. Despite its relevance, only few studies have estimated the distribution of fitness effects caused by random single mutations on protein function. The human immunodeficiency virus type 1 (HIV-1) protease was chosen as a model protein to quantify protein's tolerability to random single mutations. After determining the enzymatic activity of 107 single random mutants, we found that 86% of single mutations were deleterious for the enzyme catalytic efficiency and 54% lethal. Only 2% of the mutations significantly increased the catalytic efficiency of the enzyme. These data demonstrate the vulnerability of HIV-1 protease to single random mutations. When a second random mutagenesis library was constructed from an HIV-1 protease carrying a highly deleterious single mutation (D30N), a higher proportion of mutations with neutral or beneficial effect were found, 26% and 9%, respectively. Importantly, antagonist epistasis was observed between deleterious mutations. In particular, the mutation N88D, lethal for the wild-type protease, restored the wild-type catalytic efficiency when combined with the highly deleterious mutation D30N. The low tolerability to single random substitutions shown here for the wild-type HIV-1 protease contrasts with its in vivo ability to generate an adaptive variation. Thus, the antagonist epistasis between deleterious or lethal mutations may be responsible for increasing the protein mutational robustness and evolvability.
Key Words: protein evolution • robustness • epistasis
|
Introduction |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAcknowledgementsReferences |
|---|
Robustness defined as a protein's tolerance to substitutions is of fundamental importance to understand natural or artificial protein evolution. Recently, it has been shown that mutational robustness can enhance proteins' capacity to evolve (Aharoni et al. 2005
; Bloom, Labthavikul, et al. 2006). Protein thermodynamic stability seemed to be responsible for conferring protein robustness (Bloom et al. 2005; Wilke et al. 2005; Bloom, Drummond, et al. 2006; Bloom, Labthavikul, et al. 2006). It has also been claimed that selection might favor proteins that are robust to genetic change (Taverna and Goldstein 2002; Wilke and Adami 2003; Wagner 2005; Pal et al. 2006), which may lead to an increase in selectively neutral variants.
Proteins are to a certain extent tolerant to single mutations. Examples are bacteriophage T4 lysozyme and Escherichia coli lac repressor in which 84% and 65%, respectively, of single mutants were functional (Rennell et al. 1991; Markiewicz et al. 1994). We and others have also demonstrated the protein robustness to multiple substitutions (Martinez et al. 1996; Guo et al. 2004; Bloom et al. 2005). Although the former studies have provided important information about protein tolerance to amino acid changes, they may not reflect the protein tolerance to single substitutions because when several mutations are present simultaneously in the same genotype, they interact to determine the overall fitness, a phenomenon also termed epistasis (Phillips et al. 2000). Because the evolution of natural proteins is thought to have occurred by successive fixation of individual mutations (Smith 1970), it may be relevant to explore protein tolerance to single substitutions in order to understand the evolution of natural proteins.
To characterize protein tolerance to single random substitutions, a randomly mutated library of the human immunodeficiency virus type 1 (HIV-1) protease was constructed. We then determined the catalytic efficiency of 107 mutants carrying unique single amino acid substitutions. It has been suggested that genotypes of reduced fitness generally experience a higher fraction of mutations with beneficial effect than genotypes of high fitness (Wilke et al. 2003). Thus, a second randomly mutated library from a HIV-1 protease carrying a highly deleterious single mutation (D30N) was generated. Similarly, we analyzed the catalytic efficiency of single random variants. We decided to study the HIV-1 protease for 2 reasons. First, numerous studies have described the HIV-1 protease variability and polymorphisms found in naive or protease inhibitor–treated infected individuals (Wu et al. 2003; Ceccherini-Silberstein et al. 2004), which suggested the high evolvability of this protein. Second, HIV-1 proteins are subjected to a higher mutational burden than cellular proteins due to the error-prone nature of HIV-1 replication. It has been hypothesized that genomes experiencing high mutational burden may face selective pressure to evolve proteins that are more tolerant to change (Guo et al. 2004).
|
Materials and Methods |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAcknowledgementsReferences |
|---|
Construction of Random HIV-1 Protease Mutation Libraries in Lambda Phage
Two separate error-prone polymerase chain reaction (PCR) mutagenesis protocols that generate complementary mutational spectra were performed. In the first protocol, PCR was carried out in 10 mM Tris–HCl (pH 8.3), 50 mM KCl, 2.5 mM MgCl2, 0.5 mM MnCl2, 1 µM of each oligonucleotide, 5 units Taq polymerase (Promega, Madison, WI), and biased deoxynucleoside triphosphate concentrations (30 µM deoxycytidine triphosphate/1 mM deoxythymidine triphosphate/30 µM deoxyadenosine triphosphate/1 mM deoxy-guanosine triphosphate, Promega), as previously described (Vartanian et al. 1996
). Cycling parameters were as follows: 50x (95°C, 30 s; 55°C, 30 s; 72°C, 10 min). Input DNA was 1 ng of a pBluescript SK plasmid containing a wild-type HIV-1 protease DNA (HXB2 strain) or a mutated (D30N) HIV-1 protease DNA per 100 µl of reaction. The second PCR mutagenesis protocol was performed with Mutazyme in the GeneMorph kit (Stratagene, La Jolla, CA). Mutational spectra are described in table 1. In both protocols the following PCR oligonucleotides were employed: HIVPRLSFNF (sense, 5'-GGG GAATTCTTCCTTTAACTTCCCTCAG-3', HXB2 residues 2240–2258; underline indicates an EcoRI restriction site) and Xho8R (antisense, 5'-GGGAGGGGCTCGAGTCAAAGGCCATCCATTCCTGGC-3', HXB2 residues 2588–2604; underline indicates an XhoI restriction site; stop codons are indicated in bold face). The resulting PCR products were digested with EcoRI and XhoI, isolated, and ligated to lambda DNA (Uni-ZAP XR Vector Kit, Stratagene). The ligations were packaged (Uni-ZAP XR Giagapack Cloning Kit, Stratagene), tittered, and amplified according to standard procedures. The composition of the library was determined by nucleotide sequencing of the encoded HIV-1 protease gene included in individual phage colonies. Phage DNA from individual colonies was PCR amplified and sequenced with the flanking oligonucleotides T3 (5'-AATTAACCCTCACTAAAGGG-3') and T7 (5'-TCGAGGTCGACGGTATC-3') using the ABI PRISM dRhodamine terminator cycle sequencing kit (Applied Biosystems, Foster City, CA). Sequence alignment and editing was performed with the Sequencer version 4.1 (GeneCodes, Ann Arbor, MI) software program.
|
Determination of Protease Enzymatic Activities
The catalytic efficiencies of the different HIV-1 proteases were calculated using a phage lambda–based genetic screening, as previously described (Sices and Kristie 1998
; Martinez et al. 2000; Cabana et al. 2002; Martinez and Clotet 2003; Parera et al. 2004). Briefly, E. coli JM109 cells containing the plasmid p2X-cI.HIV were transformed with plasmid pcI.HIV-cro. The resulting cells were grown in the presence of 0.2% maltose, harvested by centrifugation, and suspended to 2.0 optical density at 600 nm (OD600)/ml in 10 mM MgSO4. Cells (200 µl) were infected with 5 x 107 pfu of phages containing the different HIV-1 proteases. After 15 min at 37°C, the cells were washed with 1 ml of 10 mM MgSO4, harvested by centrifugation, and suspended in 1 ml of Luria broth (LB) medium containing 12.5 µg of tetracycline, 0.2% maltose, 10 mM MgSO4, and 0.1 mM isopropyl-beta-D-1-thiogalactopyranoside (IPTG). The cell cultures were incubated at 37°C for 3 h and harvested by centrifugation. An additional cycle of selective growth was done by suspending the infected cells with a fresh aliquot (200 µl) of JM109 pcI.HIV-cro cells. After 2 selective growth cycles, the titer of the resulting phage was determined by coplating the cultures with 200 µl of E. coli XL-1 Blue cells (OD600 = 2.0/ml in 10 mM MgSO4) on LB plates using 3 ml top agar containing 12.5 µg of tetracycline/ml, 0.2% maltose, and 0.1 mM IPTG. After incubation at 37°C for 6 h, the plaques were counted for growth scores. The catalytic efficiency of every mutant was calculated as the mean ± standard error of at least 3 independent replicates.
Statistical Analysis
The classification of mutations as deleterious, neutral, or beneficial was based on a paired t-test included in the GraphPad Prism version 4.00 for Windows, San Diego, CA (http://www.graphpad.com/prism/Prism.htm).
|
Results |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAcknowledgementsReferences |
|---|
HIV-1 Protease Vulnerability to Single Random Amino Acid Substitutions
The HIV-1 protease is an aspartic protease consisting of 2 identical 99-amino acid monomers. The viral protease is the enzyme required for processing the gag and gag-pol polyproteins to yield mature, infectious virions. We used PCR mutagenesis to generate a randomly mutated library of the HIV-1 protease (HXB2 strain) (table 1) averaging 2.2 amino acid substitutions per mutant. Sequencing of 973 individual protease clones identified 107 mutants carrying unique single amino acid substitutions (fig. 1). The 107 single substitutions were distributed through the protein and affected 68 different residues (fig. 1). To determine the enzymatic activity of the different identified single-mutant proteases, a phage lambda–based genetic screen was used. This genetic screen is based on the phage lambda regulatory circuit in which the viral repressor cI is specifically cleaved to initiate the lysogenic to lytic switch (Sices and Kristie 1998
). Introduction of an HIV-1 protease in a wild-type phage will cleave a mutant cI repressor containing a specific HIV-1 protease cleavage site, allowing the phage to go into the lytic replication cycle. As we have previously demonstrated, the cI repressor cleavage is directly proportional to the protease catalytic efficiency (Martinez et al. 2000; Cabana et al. 2002; Martinez and Clotet 2003; Parera et al. 2004). Likewise, we have also shown that the enzyme activities determined with this genetic screen were affecting to a similar extent the ex vivo replication capacity of the virus (Cabana et al. 2002; Fernandez G, Clotet B, Martinez MA, submitted). The introduction of different protease variants in the HIV-1 infectious clone pNL4.3luc+ R–E– A confirmed that HIV-1 growth was proportional to the enzyme activity observed with the genetic screen used here (Fernandez G, Clotet B, Martinez MA, submitted) (data not shown). This HIV-1 protease genetic screen model has to be seen as a complement to the classical biochemical approach for monitoring protease proteolytic activity. The enzymatic activities of the different single-variant proteases analyzed in this study were related to the activity of the wild-type HXB2 protease (100%) (fig. 1). Fifty-eight (54.2%) mutations were lethal (t-test, P < 0.05), 34 (31.7%) were deleterious (t-test, P < 0.05), 13 (12.1%) had no significant effect (t-test, P < 0.05), and 2 (1.9%) were significantly beneficial (t-test, P < 0.05) (fig. 1). The average effect on the catalytic efficiency of the 34 deleterious mutations was –71.6%. For the 13 neutral and 2 beneficial mutations, the average effect was 1.0% and 52.5%, respectively. As expected, several lethal or highly deleterious mutants were located in conserved residues critical to the structure and function of the protease (Loeb et al. 1989; Rao et al. 1991; Ishima et al. 2001) as the active site (amino acids 22–34), the flap (amino acids 47–52), and amino acids 84–94 which encompass a single 
-helix that interacts with the substrate. Nevertheless, the former preserved coding regions were not the only residues in which lethal or deleterious mutations were found (fig. 1). Residues located outside the above 3 critical regions in which we found lethal or highly deleterious mutants (V11, E21, L38, G40, V56, Y59, G68, V75, L76, T80, P81, G94, C95, and T96) were absent or rarely mutated in drug-naive infected individuals (Wu et al. 2003; Ceccherini-Silberstein et al. 2004), confirming the possible in vivo lethal or deleterious effect caused by the former mutations. In contrast, neutral or beneficial mutations (L10, S37, L63, I64, and A71) were located in very well known polymorphic positions in the virus isolated from infected individuals (Wu et al. 2003; Ceccherini-Silberstein et al. 2004). Interestingly, 2 mutations at positions L10 and A71 could be strongly beneficial. It may be argued that our assay does not fully reflect the requirements for functionality of HIV-1 protease in vivo. Nevertheless, site-directed mutagenesis of residues L10 and A71 has demonstrated that viruses carrying these substitutions could be as infectious as the wild-type virus (Mammano et al. 2000). A 3-dimensional image representing the crystallographic HIV-1 protease structure showed that most of the neutral or beneficial mutations were located at the peripheral areas of the enzyme and almost entirely positioned in surface loops far from the active-site and substrate-binding regions (not shown). In short, these results demonstrate the low tolerability of the HIV-1 protease to single random mutations.
|
Tolerability to Single Random Amino Acid Substitutions of a Low-Fitness HIV-1 Protease
To explore the tolerance to single random amino acid substitutions of a low-fitness HIV-1 protease, we decided to generate a second randomly mutated protease library using, as starting template, a mutated HXB2 protease carrying the single substitution D30N. It has been demonstrated that proteases carrying the substitution D30N, a primary mutation associated to the protease inhibitor nelfinavir resistance, had a very low catalytic efficiency (Martinez-Picado et al. 1999
; Cabana et al. 2002). Similarly, the D30N substitution reduced the viral replication capacity (Martinez-Picado et al. 1999). In this study, we estimated that the mutant D30N HXB2 protease had a 2.6% ± 0.4% of the wild-type HXB2 proteolytic activity. Sequencing of 382 individual clones of this mutant library showed an average of 1.7 amino acid substitutions per clone (table 1).
Fifty-nine unique single mutants were identified, of which 23 were identical to single variants found in the first mutant library. Assessment of the catalytic efficiency of these 23 mutants was performed (fig. 2), and the results were compared with those obtained by the homologous 23 mutants obtained from the first library. In this second mutant library, the enzymatic activities of the different single mutants were compared with the activity of the mutant D30N HXB2 protease (100%) (fig. 2). Now, 12 (52.2%) mutations were lethal (t-test, P < 0.05), 3 (13.0%) were deleterious (t-test, P < 0.05), 6 (26.1%) had no significant effect (t-test, P<0.05), and 2 (8.7%) were significantly beneficial (t-test, P < 0.05) (fig. 2). The average effect on the catalytic efficiency of the 3 deleterious mutations was –73%; these percentages were 4.0% and 950% for the 6 neutral and the 2 beneficial mutations, respectively. When the same analysis was performed with the homologous 23 variants obtained from the wild-type HXB2 protease, a lower proportion of mutations with neutral or beneficial effect were observed. These values were 57%, 26%, 17%, and 0% for lethal, deleterious, neutral, and beneficial mutations, respectively. These results suggested that the same mutations, when expressed together with the highly deleterious mutation D30N, may have a different effect on the protease catalytic efficiency. Of note, the mutation N88D that was lethal for the wild-type HXB2 protease (fig. 1) improved the catalytic efficiency of the enzyme to numbers obtained with the wild-type HXB2 protease (fig. 2). Mutations I64V and I85V also improved the catalytic efficiency of the mutant D30N HXB2 protease, but their effect was modest. Similar to the N88D mutation, the I85V mutation was highly deleterious for the wild-type HXB2 protease (fig. 1). Both mutations, N88D and I85V, were rarely found in isolates from drug-naive infected individuals. Nevertheless, they do not confer drug resistance but can be detected in isolates from protease inhibitor–treated individuals (Wu et al. 2003). This result showed that the former mutations may act as compensatory or stabilizing substitutions in the presence of accompanying deleterious mutations, suggesting that antagonistic epistasis may buffer the effects of deleterious mutations.
|
|
Discussion |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAcknowledgementsReferences |
|---|
There are many examples of proteins highly robust to mutations. They include several enzymes that can tolerate many amino acid changes (Martinez et al. 1996
; Guo et al. 2004). Protein robustness seems to be a selectable trait because neutral mutations can be a key to future evolutionary innovations (Wagner 2005). A recent evolution experiment demonstrated the evolutionary advantages of neutral mutations by showing that human or bacterial enzymes can acquire new functions without loosing their original functions (Aharoni et al. 2005). This apparent selective advantage of robustness contrasts with the low tolerability to single substitutions demonstrated in this study for the HIV-1 protease, even with the high mutational burden faced by this protein. Here, we show that up to 54% of random single amino acid substitutions were lethal. Moreover, the average catalytic efficiency effect for the 32% of deleterious but nonlethal mutations was –72%. This high vulnerability to single mutations contrasts with the high HIV-1 protease genetic variability found within infected individuals (Wu et al. 2003; Ceccherini-Silberstein et al. 2004). Especially intriguing is the extremely rapid evolvability displayed by this protein after protease inhibitor treatment (Condra et al. 1995). HIV-1 protease can rapidly acquire mutations that lead to drug resistance but that barely affect its catalytic efficiency (Martinez-Picado et al. 1999; Nijhuis et al. 1999). Moreover, up to 22 (22%) different residues have been associated with drug resistance (Johnson et al. 2005). Although the high mutation rates and large population size of HIV-1 are obviously favoring the rapid evolvability of the HIV-1 protease (Domingo et al. 1997), protein evolution also depends on a reduced lethality of mutations (Aharoni et al. 2005; Wagner 2005; Bloom, Labthavikul, et al. 2006).
After showing the low tolerability to single random substitutions of a high-fitness protein, we decided to explore the tolerance to single random amino acid substitutions of a low-fitness HIV-1 protease. The chosen low-fitness protein was the HXB2 protease carrying the single mutation D30N. This substitution is primarily associated to resistance to the HIV-1 protease inhibitor nelfinavir (Johnson et al. 2005). It is very well known that this mutation impairs the enzyme's catalytic efficiency and the virus replication capacity (Martinez-Picado et al. 1999; Cabana et al. 2002). Interestingly, we found a higher proportion of mutations with a neutral or beneficial effect with the low-fitness protease. When epistatic relationships were searched between the highly deleterious mutation D30N and the different single random mutations, strong antagonistic epistasis was found with some mutations. In particular, the substitutions N88D and D30N turned out to be strongly epistatic. Both mutations reduced the catalytic efficiency of the wild-type protease (N88D was lethal and D30N strongly deleterious), but together, their effect was better. Because the mutation D30N confers some resistance to the drug and emerges first during in vivo protease inhibitor resistance development (Martinez-Picado et al. 1999; Cabana et al. 2002), the appearance of the mutation N88D would restore the catalytic defect originated by the D30N substitution. Indeed, the D30N and N88D mutations are specifically associated with resistance to nelfinavir in HIV-1–infected individuals (Patick et al. 1998).
In a number of recent studies carried out with RNA viruses, including HIV-1, a tendency toward antagonistic epistasis has been observed (Bonhoeffer et al. 2004; Burch and Chao 2004; Sanjuan et al. 2004). The existence of antagonistic epistasis is not restricted to RNA viruses, and it has also been found in bacteria, eukaryotic proteins, and digital organisms (Maisnier-Patin et al. 2002; Lenski et al. 2003; Bridgham et al. 2006). Especially interesting is the recent work performed by Bridham et al. (2006). They demonstrate the role of antagonistic epistasis in the evolution of eukaryotic proteins such as the hormone receptors. In this study the antagonist epistatic interaction of 2 mutations introduced into an ancestral sequence recapitulates the evolution of present-day receptor specificity. Similar to our findings, an occasional highly deleterious mutation was rescued by a partner mutation that conferred a beneficial trait.
It has been hypothesized that antagonistic epistasis is characteristic of hypersensitive genomes (Elena et al. 2006; Sanjuan and Elena 2006); our findings may explain how antagonist epistasis can increase the mutational robustness and evolvability of a protein highly vulnerable to single random mutations.
|
Acknowledgements |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAcknowledgementsReferences |
|---|
This work was supported by Spanish Minsiterio de Educación y Ciencia project BFU2006-01066/BMC and Fondo de Investigación Sanitaria project PI050022.
|
Footnotes |
|---|
Edward Holmes, Associate Editor
|
References |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAcknowledgementsReferences |
|---|
Aharoni A, Gaidukov L, Khersonsky O, Gould SM, Roodveldt C, Tawfik DS. (2005) The ‘evolvability’ of promiscuous protein functions. Nat Genet 37:73–76.[Web of Science][Medline]
Bloom JD, Drummond DA, Arnold FH, Wilke CO. (2006) Structural determinants of the rate of protein evolution in yeast. Mol Biol Evol 23:1751–1761.
Bloom JD, Labthavikul ST, Otey CR, Arnold FH. (2006) Protein stability promotes evolvability. Proc Natl Acad Sci USA 103:5869–5874.
Bloom JD, Silberg JJ, Wilke CO, Drummond DA, Adami C, Arnold FH. (2005) Thermodynamic prediction of protein neutrality. Proc Natl Acad Sci USA 102:606–611.
Bonhoeffer S, Chappey C, Parkin NT, Whitcomb JM, Petropoulos CJ. (2004) Evidence for positive epistasis in HIV-1. Science 306:1547–1550.
Bridgham JT, Carroll SM, Thornton JW. (2006) Evolution of hormone-receptor complexity by molecular exploitation. Science 312:97–101.
Burch CL and Chao L. (2004) Epistasis and its relationship to canalization in the RNA virus phi 6. Genetics 167:559–567.
Cabana M, Fernandez G, Parera M, Clotet B, Martinez MA. (2002) Catalytic efficiency and phenotype of HIV-1 proteases encoding single critical resistance substitutions. Virology 300:71–78.[CrossRef][Web of Science][Medline]
Ceccherini-Silberstein F, Erba F, Gago F, et al. (13 co-authors). (2004) Identification of the minimal conserved structure of HIV-1 protease in the presence and absence of drug pressure. Aids 18:F11–F19.[Web of Science][Medline]
Condra JH, Schleif WA, Blahy OM, et al. (16 co-authors). (1995) In vivo emergence of HIV-1 variants resistant to multiple protease inhibitors. Nature 374:569–571.[CrossRef][Medline]
Domingo E, Menendez-Arias L, Quinones-Mateu ME, Holguin A, Gutierrez-Rivas M, Martinez MA, Quer J, Novella IS, Holland JJ. (1997) Viral quasispecies and the problem of vaccine-escape and drug-resistant mutants. Prog Drug Res 48:99–128.[Medline]
Elena SF, Carrasco P, Daros JA, Sanjuan R. (2006) Mechanisms of genetic robustness in RNA viruses. EMBO Rep 7:168–173.[CrossRef][Web of Science][Medline]
Guo HH, Choe J, Loeb LA. (2004) Protein tolerance to random amino acid change. Proc Natl Acad Sci USA 101:9205–9210.
Ishima R, Ghirlando R, Tozser J, Gronenborn AM, Torchia DA, Louis JM. (2001) Folded monomer of HIV-1 protease. J Biol Chem 276:49110–49116.
Johnson VA, Brun-Vezinet F, Clotet B, Conway B, Kuritzkes DR, Pillay D, Schapiro JM, Telenti A, Richman DD. (2005) Update of the drug resistance mutations in HIV-1: Fall 2005. Top HIV Med 13:125–131.[Medline]
Lenski RE, Ofria C, Pennock RT, Adami C. (2003) The evolutionary origin of complex features. Nature 423:139–144.[CrossRef][Medline]
Loeb DD, Swanstrom R, Everitt L, Manchester M, Stamper SE, Hutchison CA 3rd. (1989) Complete mutagenesis of the HIV-1 protease. Nature 340:397–400.[CrossRef][Medline]
Maisnier-Patin S, Berg OG, Liljas L, Andersson DI. (2002) Compensatory adaptation to the deleterious effect of antibiotic resistance in Salmonella typhimurium. Mol Microbiol 46:355–366.[CrossRef][Web of Science][Medline]
Mammano F, Trouplin V, Zennou V, Clavel F. (2000) Retracing the evolutionary pathways of human immunodeficiency virus type 1 resistance to protease inhibitors: virus fitness in the absence and in the presence of drug. J Virol 74:8524–8531.
Markiewicz P, Kleina LG, Cruz C, Ehret S, Miller JH. (1994) Genetic studies of the lac repressor. XIV. Analysis of 4000 altered Escherichia coli lac repressors reveals essential and non-essential residues, as well as "spacers" which do not require a specific sequence. J Mol Biol 240:421–433.[CrossRef][Web of Science][Medline]
Martinez MA, Cabana M, Parera M, Gutierrez A, Este JA, Clotet B. (2000) A bacteriophage lambda-based genetic screen for characterization of the activity and phenotype of the human immunodeficiency virus type 1 protease. Antimicrob Agents Chemother 44:1132–1139.
Martinez MA and Clotet B. (2003) Genetic screen for monitoring hepatitis C virus NS3 serine protease activity. Antimicrob Agents Chemother 47:1760–1765.
Martinez MA, Pezo V, Marliere P, Wain-Hobson S. (1996) Exploring the functional robustness of an enzyme by in vitro evolution. EMBO J 15:1203–1210.[Web of Science][Medline]
Martinez-Picado J, Savara AV, Sutton L, D'Aquila RT. (1999) Replicative fitness of protease inhibitor-resistant mutants of human immunodeficiency virus type 1. J Virol 73:3744–3752.
Nijhuis M, Schuurman R, de Jong D, Erickson J, Gustchina E, Albert J, Schipper P, Gulnik S, Boucher CA. (1999) Increased fitness of drug resistant HIV-1 protease as a result of acquisition of compensatory mutations during suboptimal therapy. Aids 13:2349–2359.[CrossRef][Web of Science][Medline]
Pal C, Papp B, Lercher MJ. (2006) An integrated view of protein evolution. Nat Rev Genet 7:337–348.[CrossRef][Web of Science][Medline]
Parera M, Clotet B, Martinez MA. (2004) Genetic screen for monitoring severe acute respiratory syndrome coronavirus 3C-like protease. J Virol 78:14057–14061.
Patick AK, Duran M, Cao Y, Shugarts D, Keller MR, Mazabel E, Knowles M, Chapman S, Kuritzkes DR, Markowitz M. (1998) Genotypic and phenotypic characterization of human immunodeficiency virus type 1 variants isolated from patients treated with the protease inhibitor nelfinavir. Antimicrob Agents Chemother 42:2637–2644.
Phillips P, Otto S, Whitlock M. (2000) Beyond the average: the evolutionary importance of gene interactions and variability of epistatic effects. In Wolf J, Brodie E, Wade M (Eds.). Epistasis and the evolutionary process(Oxford University Press, Oxford) pp. 20–38.
Rao JK, Erickson JW, Wlodawer A. (1991) Structural and evolutionary relationships between retroviral and eucaryotic aspartic proteinases. Biochemistry 30:4663–4671.[CrossRef][Medline]
Rennell D, Bouvier SE, Hardy LW, Poteete AR. (1991) Systematic mutation of bacteriophage T4 lysozyme. J Mol Biol 222:67–88.[CrossRef][Web of Science][Medline]
Sanjuan R and Elena SF. (2006) Epistasis correlates to genomic complexity. Proc Natl Acad Sci USA 103:14402–14405.
Sanjuan R, Moya A, Elena SF. (2004) The contribution of epistasis to the architecture of fitness in an RNA virus. Proc Natl Acad Sci USA 101:15376–15379.
Sices HJ and Kristie TM. (1998) A genetic screen for the isolation and characterization of site-specific proteases. Proc Natl Acad Sci USA 95:2828–2833.
Smith JM. (1970) Natural selection and the concept of a protein space. Nature 225:563–564.[CrossRef][Medline]
Taverna DM and Goldstein RA. (2002) Why are proteins so robust to site mutations? J Mol Biol 315:479–484.[CrossRef][Web of Science][Medline]
Vartanian JP, Henry M, Wain-Hobson S. (1996) Hypermutagenic PCR involving all four transitions and a sizeable proportion of transversions. Nucleic Acids Res 24:2627–2631.
Wagner A. (2005) Robustness, evolvability, and neutrality. Febs Lett 579:1772–1778.[CrossRef][Web of Science][Medline]
Wilke CO and Adami C. (2003) Evolution of mutational robustness. Mutat Res 522:3–11.[Web of Science][Medline]
Wilke CO, Bloom JD, Drummond DA, Raval A. (2005) Predicting the tolerance of proteins to random amino acid substitution. Biophys J 89:3714–3720.
Wilke CO, Lenski RE, Adami C. (2003) Compensatory mutations cause excess of antagonistic epistasis in RNA secondary structure folding. BMC Evol Biol 3:3.[CrossRef][Medline]
Wu TD, Schiffer CA, Gonzales MJ, Taylor J, Kantor R, Chou S, Israelski D, Zolopa AR, Fessel WJ, Shafer RW. (2003) Mutation patterns and structural correlates in human immunodeficiency virus type 1 protease following different protease inhibitor treatments. J Virol 77:4836–4847.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  G. Fernandez, B. Clotet, and M. A. Martinez Fitness Landscape of Human Immunodeficiency Virus Type 1 Protease Quasispecies J. Virol., March 1, 2007; 81(5): 2485 - 2496. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||