Dataset Information

Prediction of mutational tolerance in HIV-1 protease and reverse transcriptase using flexible backbone protein design.

ABSTRACT: Predicting which mutations proteins tolerate while maintaining their structure and function has important applications for modeling fundamental properties of proteins and their evolution; it also drives progress in protein design. Here we develop a computational model to predict the tolerated sequence space of HIV-1 protease reachable by single mutations. We assess the model by comparison to the observed variability in more than 50,000 HIV-1 protease sequences, one of the most comprehensive datasets on tolerated sequence space. We then extend the model to a second protein, reverse transcriptase. The model integrates multiple structural and functional constraints acting on a protein and uses ensembles of protein conformations. We find the model correctly captures a considerable fraction of protease and reverse-transcriptase mutational tolerance and shows comparable accuracy using either experimentally determined or computationally generated structural ensembles. Predictions of tolerated sequence space afforded by the model provide insights into stability-function tradeoffs in the emergence of resistance mutations and into strengths and limitations of the computational model.

SUBMITTER: Humphris-Narayanan E

PROVIDER: S-EPMC3426558 | biostudies-literature | 2012

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Prediction of mutational tolerance in HIV-1 protease and reverse transcriptase using flexible backbone protein design.

Humphris-Narayanan Elisabeth E Akiva Eyal E Varela Rocco R Ó Conchúir Shane S Kortemme Tanja T

PLoS computational biology 20120823 8

Predicting which mutations proteins tolerate while maintaining their structure and function has important applications for modeling fundamental properties of proteins and their evolution; it also drives progress in protein design. Here we develop a computational model to predict the tolerated sequence space of HIV-1 protease reachable by single mutations. We assess the model by comparison to the observed variability in more than 50,000 HIV-1 protease sequences, one of the most comprehensive data ...[more]

PMID: 22927804

Similar Datasets

Project description:Amino acid Lys(65) is part of the highly flexible beta3-beta4 loop in the fingers domain of the 66 kDa subunit of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT). Recent crystal data show that the epsilon-amino group of Lys(65) interacts with the gamma-phosphate of the bound deoxynucleoside triphosphate ('dNTP') substrate [Huang, Chopra, Verdine and Harrison (1998) Science 282, 1669-1675]. In order to biochemically define the function of RT Lys(65), we have used site-specific mutagenesis to generate RT with a variety of substitutions at this position, including K65E, K65Q, K65A and K65R. Kinetic analyses demonstrate that if Lys(65) in RT is substituted with an amino acid other than arginine the enzyme exhibits dramatic decreases in the binding affinity (K(m)) for all dNTP substrates, in RT catalytic efficiency (k(cat)/K(m)) and in the mutant enzyme's ability to carry out pyrophosphorolysis, the reverse reaction of DNA synthesis. The pH optimum for the DNA polymerase activity of K65E RT was 6.5, compared to 7.5 for the wild-type enzyme, and 8.0 for the K65R, K65A and K65Q mutants. Molecular modelling studies show that mutations of Lys(65) do not affect the geometry of the loop's alpha-carbon backbone, but rather lead to changes in positioning of the side chains of residues Lys(70) and Arg(72). In particular, Glu in K65E can form a salt bridge with Arg(72), leading to the diminution of the latter residue's interaction with the alpha-phosphate of the dNTP residue. This alteration in dNTP-binding may explain the large pH-dependent changes in both dNTP-binding and catalytic efficiency noted with the enzyme. Furthermore, the K65A, K65Q and K65E mutant enzymes are 100-fold less sensitive to all dideoxynucleoside triphosphate ('ddNTP') inhibitors, whereas the K65R mutation results in a selective 10-fold decrease in binding of ddCTP and ddATP only. This implies that mutations at position 65 in HIV-1 RT influence the nucleotide-binding specificity of the enzyme.

Project description:UnlabelledHIV-1 protease (PR), reverse transcriptase (RT), and integrase (IN) variability presents a challenge to laboratories performing genotypic resistance testing. This challenge will grow with increased sequencing of samples enriched for proviral DNA such as dried blood spots and increased use of next-generation sequencing (NGS) to detect low-abundance HIV-1 variants. We analyzed PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to characterize variation at each amino acid position, identify mutations indicating APOBEC-mediated G-to-A editing, and identify mutations resulting from selective drug pressure. Forty-seven percent of PR, 37% of RT, and 34% of IN positions had one or more amino acid variants with a prevalence of ≥1%. Seventy percent of PR, 60% of RT, and 60% of IN positions had one or more variants with a prevalence of ≥0.1%. Overall 201 PR, 636 RT, and 346 IN variants had a prevalence of ≥0.1%. The median intersubtype prevalence ratios were 2.9-, 2.1-, and 1.9-fold for these PR, RT, and IN variants, respectively. Only 5.0% of PR, 3.7% of RT, and 2.0% of IN variants had a median intersubtype prevalence ratio of ≥10-fold. Variants at lower prevalences were more likely to differ biochemically and to be part of an electrophoretic mixture compared to high-prevalence variants. There were 209 mutations indicative of APOBEC-mediated G-to-A editing and 326 mutations nonpolymorphic treatment selected. Identification of viruses with a high number of APOBEC-associated mutations will facilitate the quality control of dried blood spot sequencing. Identifying sequences with a high proportion of rare mutations will facilitate the quality control of NGS.ImportanceMost antiretroviral drugs target three HIV-1 proteins: PR, RT, and IN. These proteins are highly variable: many different amino acids can be present at the same position in viruses from different individuals. Some of the amino acid variants cause drug resistance and occur mainly in individuals receiving antiretroviral drugs. Some variants result from a human cellular defense mechanism called APOBEC-mediated hypermutation. Many variants result from naturally occurring mutation. Some variants may represent technical artifacts. We studied PR and RT sequences from >100,000 individuals and IN sequences from >10,000 individuals to quantify variation at each amino acid position in these three HIV-1 proteins. We performed analyses to determine which amino acid variants resulted from antiretroviral drug selection pressure, APOBEC-mediated editing, and naturally occurring variation. Our results provide information essential to clinical, research, and public health laboratories performing genotypic resistance testing by sequencing HIV-1 PR, RT, and IN.

Project description:BackgroundN348I in HIV-1 reverse transcriptase (RT) confers resistance to zidovudine (AZT) and nevirapine. Biochemical studies demonstrated that N348I indirectly increases AZT resistance by decreasing the frequency of secondary ribonuclease H (RNase H) cleavages that reduce the RNA/DNA duplex length of the template/primer (T/P) and diminish the efficiency of AZT-monophosphate (MP) excision. By contrast, there is some discrepancy in the literature in regard to the mechanisms associated with nevirapine resistance: one study suggested that it is due to decreased inhibitor binding while others suggest that it may be related to the decreased RNase H cleavage phenotype. From a structural perspective, N348 in both subunits of RT resides distal to the enzyme's active sites, to the T/P binding tract and to the nevirapine-binding pocket. As such, the structural mechanisms associated with the resistance phenotypes are not known.ResultsUsing a novel modelled structure of RT in complex with an RNA/DNA T/P, we identified a putative interaction between the β14-β15 loop in the p51 subunit of RT and the RNA template. Substitution of the asparagine at codon 348 in the p51 subunit with either isoleucine or leucine abrogated the observed protein-RNA interaction, thus, providing a possible explanation for the decreased RNase H phenotype. By contrast, alanine or glutamine substitutions exerted no effect. To validate this model, we introduced the N348I, N348L, N348A and N348Q mutations into RT and purified enzymes that contained subunit-specific mutations. N348I and N348L significantly decreased the frequency of secondary RNase H cleavages and increased the enzyme's ability to excise AZT-MP. As predicted by the modelling, this phenotype was due to the mutation in the p51 subunit of RT. By contrast, the N348A and N348Q RTs exhibited RNase H cleavage profiles and AZT-MP excision activities similar to the wild-type enzyme. All N348 mutant RTs exhibited decreased nevirapine susceptibility, although the N348I and N348L mutations conferred higher fold resistance values compared to N348A and N348Q. Nevirapine resistance was also largely due to the mutation present in the p51 subunit of RT.ConclusionsThis study demonstrates that N348I-mediated AZT and nevirapine resistance is due to the mutation in the p51 subunit of RT.

Project description:Despite the high degree of HIV-1 protease and reverse transcriptase (RT) mutation in the setting of antiretroviral therapy, the spectrum of possible virus variants appears to be limited by patterns of amino acid covariation. We analyzed patterns of amino acid covariation in protease and RT sequences from more than 7,000 persons infected with HIV-1 subtype B viruses obtained from the Stanford HIV Drug Resistance Database (http://hivdb.stanford.edu). In addition, we examined the relationship between conditional probabilities associated with a pair of mutations and the order in which those mutations developed in viruses for which longitudinal sequence data were available. Patterns of RT covariation were dominated by the distinct clustering of Type I and Type II thymidine analog mutations and the Q151M-associated mutations. Patterns of protease covariation were dominated by the clustering of nelfinavir-associated mutations (D30N and N88D), two main groups of protease inhibitor (PI)-resistance mutations associated either with V82A or L90M, and a tight cluster of mutations associated with decreased susceptibility to amprenavir and the most recently approved PI darunavir. Different patterns of covariation were frequently observed for different mutations at the same position including the RT mutations T69D versus T69N, L74V versus L74I, V75I versus V75M, T215F versus T215Y, and K219Q/E versus K219N/R, and the protease mutations M46I versus M46L, I54V versus I54M/L, and N88D versus N88S. Sequence data from persons with correlated mutations in whom earlier sequences were available confirmed that the conditional probabilities associated with correlated mutation pairs could be used to predict the order in which the mutations were likely to have developed. Whereas accessory nucleoside RT inhibitor-resistance mutations nearly always follow primary nucleoside RT inhibitor-resistance mutations, accessory PI-resistance mutations often preceded primary PI-resistance mutations.

Dataset Information

Prediction of mutational tolerance in HIV-1 protease and reverse transcriptase using flexible backbone protein design.

Publications

Prediction of mutational tolerance in HIV-1 protease and reverse transcriptase using flexible backbone protein design.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets