Next Article in Journal
Human Papillomavirus and Tonsillar and Base of Tongue Cancer
Previous Article in Journal
The Cryptophlebia Leucotreta Granulovirus—10 Years of Commercial Field Use
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 7
Issue 3
10.3390/v7031313
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Selection Pressure in CD8+ T-cell Epitopes in the pol Gene of HIV-1 Infected Individuals in Colombia. A Bioinformatic Approach

by
Liliana Acevedo-Sáenz
1,
Rodrigo Ochoa
2,
Maria Teresa Rugeles
1,
Patricia Olaya-García
3,
Paula Andrea Velilla-Hernández
1 and
Francisco J. Diaz
1,*
1
Grupo Inmunovirología, Facultad de Medicina, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, 050010, Colombia
2
Programa de Estudio y Control de Enfermedades Tropicales—PECET, Universidad de Antioquia, Medellín, 050010, Colombia
3
Centro de Análisis Molecular, Bogotá D.C., 111156, Colombia
*
Author to whom correspondence should be addressed.
Viruses 2015, 7(3), 1313-1331; https://doi.org/10.3390/v7031313
Submission received: 17 September 2014 / Revised: 23 February 2015 / Accepted: 24 February 2015 / Published: 20 March 2015
(This article belongs to the Section Animal Viruses)
Browse Figure
Versions Notes

Abstract

:
One of the main characteristics of the human immunodeficiency virus is its genetic variability and rapid adaptation to changing environmental conditions. This variability, resulting from the lack of proofreading activity of the viral reverse transcriptase, generates mutations that could be fixed either by random genetic drift or by positive selection. Among the forces driving positive selection are antiretroviral therapy and CD8+ T-cells, the most important immune mechanism involved in viral control. Here, we describe mutations induced by these selective forces acting on the pol gene of HIV in a group of infected individuals. We used Maximum Likelihood analyses of the ratio of non-synonymous to synonymous mutations per site (dN/dS) to study the extent of positive selection in the protease and the reverse transcriptase, using 614 viral sequences from Colombian patients. We also performed computational approaches, docking and algorithmic analyses, to assess whether the positively selected mutations affected binding to the HLA molecules. We found 19 positively-selected codons in drug resistance-associated sites and 22 located within CD8+ T-cell epitopes. A high percentage of mutations in these epitopes has not been previously reported. According to the docking analyses only one of those mutations affected HLA binding. However, algorithmic methods predicted a decrease in the affinity for the HLA molecule in seven mutated peptides. The bioinformatics strategies described here are useful to identify putative positively selected mutations associated with immune escape but should be complemented with an experimental approach to define the impact of these mutations on the functional profile of the CD8+ T-cells.

1. Introduction

Human immunodeficiency virus 1 (HIV-1) is characterized by high levels of diversity and rapid evolution. For HIV-1, four genetic groups have been described: group M causes most of the infections worldwide [1] and is divided into nine subtypes: A–D, F–H, J and K [2]; groups O, N and P which are restricted to the West and Central Africa and have a lower epidemiological impact [3,4,5]. Intersubtype recombinants originated in individuals co-infected with two or more subtypes are also common [6]. Two studies carried out in the cities of Medellín and Bogotá (Colombia) with HIV sequences of 2001 and 2002 classified the circulating strains within subtype B [7,8].
The genetic diversity of HIV is partially due to the high rates of replication and mutation [9]. The lack of proofreading activity of the viral reverse transcriptase (RT) is the major mechanism of mutation [10,11]. The mutation rate is defined as the number of substitutions per nucleotide site per replicative cycle (s/s/c); in HIV, this rate varies between 2.4 × 10−5 and 3.4 × 10−5 (s/s/c) [12,13,14]. From this, it follows that between one third and one fourth of the new viral particles carry at least one mutation.
These mutations contribute to HIV diversification and have been associated with evasion of the immune response [15], defective virus particles [16,17], changes of cellular tropism and resistance to antiretroviral drugs [18]. Some mutations are positively selected by the antiretroviral therapy or the immune response, driving HIV evolution [19,20]. Specifically, CD8+ T-cell responses to HIV-1, restricted through human leukocyte antigen (HLA) molecules, are responsible for the positive selection of escape mutations observed in HIV [21,22]. These mutations are often located in sites of inter- and intra-subtype variability [23]. The rapid generation of escape mutants by HIV-1 strains is an enormous challenge to be overcome by the immune system, and also in the design of a globally effective HIV vaccine [24].
Moore and collaborators determined the association between specific HLA alleles in HIV-1-infected individuals and the presence or absence of some amino acid replacements in the virus sequences at particular positions [22]; this study confirmed selection of escape mutations restricted by HLA class I molecules, and underscores the high proportion of viral codons whose evolution is driven by immune escape. In fact, between 10% and 40% of the variation in codons of HIV-1 subtype B strains is driven by immune pressure associated with class I HLA alleles [25]. The clinical implication of these mutations associated with the immune escape is the increase in the rate of progression of the infection; this phenomenon has most frequently been observed in individuals expressing HLA-B*27 and HLA-B*51 in which amino acid substitutions in immunodominant epitopes led to high levels of HIV replication and rapid progression to AIDS [26,27,28].
The aim of our study was to explore the presence of mutations presumably associated with immune escape in the protease (PRO) and RT proteins of HIV-1 of chronically infected individuals in Colombia. For this purpose, we designed an in silico strategy consisting in detecting positively selected mutations in T CD8+ epitopes, followed by the prediction of changes in the affinity between the mutated peptide and the corresponding HLA molecule. We identified some selected peptides not previously reported as escape mutants within epitopes restricted to class I HLA molecules with high prevalence in the Colombian population and predicted their consequences in the process of antigenic presentation.

2. Materials and Methods

2.1. Data Sources and Sequence Alignments

We analyzed 614 sequences from the pol gene (genome positions 2262–2549 and 2661–3290) coding for the PRO and part of the RT. These sequences were provided by the Centro de Análisis Molecular, a specialized laboratory located in Bogotá, the capital, which performed antiviral sensitivity tests in samples from different regions of Colombia. These samples were collected between 2000 and 2007 from individuals receiving antiretroviral therapy. Sequencing was performed by the TRUGENE® kit (Siemens Healthcare Diagnostics, Tarrytown, NY, USA) [29]. A code was assigned to each sequence in order to be used anonymously; clinical data and result of additional laboratory studies were not available. Sequences were aligned by the Clustal W method implemented in the BioEdit package (Ibis Bioscience, Carlsbad, CA, USA), using the reference sequence HXB2 of HIV-1 (GenBank accession number K03455.1). Analysis of antiviral resistance data, subtyping and phylodynamics of this database have already been published [30,31]. Most (82%) sequences exhibit resistance to at least one antiviral drug [30]. Almost all of these pol gene sequences were subtype B; the exception was one subtype F. No evidence of inter-subtype recombination was found [31].

2.2. Tests of Positive Selection

We first used the Z-test, as implemented in the MEGA 5.2.2 software [32], to explore evidence of natural selection in the dataset. The statistics (dN–dS), where dS and dN are the numbers of synonymous substitutions per synonymous site and nonsynonymous substitutions per nonsynonymous site, respectively, were estimated using the Nei-Gojobori method [33]. The variance of the difference was computed using the bootstrap method (10,000 replicates).
We also searched for natural selection by HyPhy, a Maximum Likelihood (ML) based method [34]. For each codon, estimates of dS and dN were produced using the joint ML reconstructions of ancestral states under the Muse-Gaut model [35] of codon substitution and the General Time Reversible model [36] of nucleotide substitution. Positive selection was inferred for codons with dN/dS values significantly greater than 1 and the probability of rejecting the null hypothesis of neutral evolution (p-value) was calculated [37,38]; p values < 0.05 were considered significant.
In both the Z-test and HyPhy analyses, all positions with less than 90% site coverage were eliminated. That is, fewer than 10% alignment gaps, missing data, and ambiguous bases were allowed at any position. There was a total of 304 codons in the final dataset.
To confirm the results obtained with the MEGA implementation of HyPhy, the Single Likelihood Ancestor Counting (SLAC) method implemented in the Datamonkey webserver [37] was used [39]. Since Datamonkey does not process datasets with more than 500 sequences, for this analysis the sequences were arbitrarily divided in two groups. The nucleotide and codons substitution models were the same used in the analysis performed in MEGA 5.2.2 software.

2.3. Identification of Peptides

We identified HLA-class I restricted epitopes using the HIV Molecular Immunology Database of Los Alamos [40]. Each epitope was located in the reference sequence HXB2 of HIV-1 (GenBank accession number: K03455.1) and compared to the homologous sequences in our dataset. For each epitope sequence, we identified the mutations, measured their frequency and classified them as susceptible forms (SF), immune escape (IE) or drug resistance associated (DRAS) mutations. Susceptible forms, defined as epitopes with mutations that are still able to induce a CD8+ T-cell response, and immune escape, defined as mutated epitopes with documented or inferred escape capabilities, were derived from the Los Alamos Database. Substitutions were considered as drug resistance mutations if they were included in the Stanford HIV drug resistance database [41,42] or in the Update of the Drug Resistance Mutations in HIV-1 of the International AIDS Society [43]; other substitutions observed but not included in any of the previous definitions were dubbed as “non-classified” (NC) mutations.

2.4. Collection and Preparation of the Three-Dimensional Structure of HLA Molecules

From the Protein Data Bank (PDB) database [44] we downloaded the three-dimensional structure of certain HLA molecules found with high frequency in Colombia. We only considered proteins with a resolution ≤ 2.0 Å. The proteins were edited using the AutoDock Tools program [45]; external ligands and crystallized water were eliminated. AutoDock tools were used to add hydrogen to polar atoms of the molecules, and Gasteiger changes to HLA proteins, facilitating its potential interactions with the peptides. The PDB codes were 1E27, 1DUZ, 3UTQ, 3MRK, 1EFX, 1M6O, 3BO8, 3RL1, 1XR9, 3C9N and 4HX1.

2.5. Structural Prediction of Peptides

Since we required the 3D structures of all peptides included in the study for the in silico docking simulation assay, we used the homology-modeling I-TASSER [46,47] and ESyPred3D [48] servers. With this strategy, we predicted the 3D structure of the wild type (WT) and mutated peptides in order to assess how the mutations affect their folding and the affinity of their binding to the MHC molecules.

2.6. HLA-Peptide Binding Predictions

We predicted the affinities between the HLA protein and the WT or mutant peptides using the receptor-ligand docking simulation AutoDockVina software [49] and three different algorithms of affinity estimation: NetMHCpan 2.8 [50,51], NetMHC [52] and stabilized matrix method (SMM) [53,54]. For AutoDockVina, we evaluated the affinity constant based on the indirect score of the free energy of binding. The conformational search space in the AutoDockVina simulation was limited to the vicinity of the major histocompatibility complex (MHC) pocket. For each peptide, the docking simulations were run 20 times and the best scores of each run were selected and averaged; the interactions between residues from both, ligand and receptor were verified in detail. Statistical comparison of the free energy of binding of the WT and mutated epitopes were performed by the Kruskal-Wallis and Mann-Whitney U tests.
For the NetMHCpan, NetMHC and SMM algorithms, we considered as good, intermediate and poor bindings when the Inhibitory Concentration 50 (IC50) was < 50 nM, 50–500 nM or > 500 nM, respectively. Two-fold or larger differences between the affinities of the WT and mutated peptides with the corresponding HLA molecule were considered as significant [55]. Correlation analyses among the different methods of affinity assessment were performed by the Spearman test. All statistics were performed using GraphPrism 5.0 package; p-values lower than 0.05 were considered significant.

3. Results

3.1. Z-Test and Maximum Likelihood Analysis of Positive Selection

We analyzed 614 pol gene sequences corresponding to the PRO (codons 4–99) and the RT (codons 38–247) obtained from HIV-infected individuals living in 16 Colombian cities, mainly Bogotá (44.4%), Cali (12.6%) and Medellín (10.6%). There were 280 variable sites out of a total of 304 codons.
Initially, we searched for evidence of selection using the Z-test in the overall average of all sequences. The null hypothesis of strict neutrality (dN = dS) was rejected in favor of the alternative hypothesis (dN ≠ dS) with a p = 0.000. The value for the statistic dN–dS was −10.5, indicating a higher rate of synonymous than non-synonymous substitutions. When tested for positive selection the null hypothesis of neutrality in favor of the alternative (dN > dS) could not be rejected (p = 1.000). Conversely, when tested for purifying selection the null hypothesis of neutrality was rejected (p = 0.000) in favor of the alternative (dN < dS), indicating that negative selection predominates along the studied segments of pol.
However, it is well known that a purifying selection in most parts of the sequence could hide positive selection occurring at specific codons [33,36]. To obtain an assessment of codon-specific selective pressures, we performed an analysis of dN/dS by the Maximum Likelihood method using the HyPhy implementation in MEGA. The test of positive selection was significant (p < 0.05) in 35 codons (18 in PRO and 17 in RT) (Table 1). For both proteins, sites under positive selection were spread throughout the gene.
Table
Table 1. Mutations identified as positively selected by the Maximum Likelihood method in protease (PRO) and reverse transcriptase (RT) of 614 human immunodeficiency virus (HIV) infected individuals.
Table 1. Mutations identified as positively selected by the Maximum Likelihood method in protease (PRO) and reverse transcriptase (RT) of 614 human immunodeficiency virus (HIV) infected individuals.
Codons Position (HXB2)SubstitutiondN/dSap ValueLocation b
Protease
12T → P/S4.40.0000 *Epitope
13I → V5.00.0002Epitope
19L → I/Q/V1.70.0015Epitope
35E → D1.90.0038Epitope
37S→N/E/D5.70.0000 *Epitope
41R→K2.90.0000 *Epitope
54I → L/V/M/T/A/S3.00.0000 *Epitope/DRAS
62I → V2.90.0002DRAS
64I → L/M/V6.80.0000 *DRAS
71A → V/I/T/L2.30.0000 *DRAS
72I → V/T/L/R2.10.0115Epitope
73G → C/S/T/A2.10.0146DRAS
74T → P15.60.0000 *DRAS
77V → I5.30.0000 *Epitope/DRAS
82V → A/F/S/T4.20.0000 *Epitope/DRAS
85I → V360.0460DRAS
90L → M4.40.0000 *DRAS
93I → L/M8.20.0000 *DRAS
Reverse transcriptase
39T → A/K/S/L4.40.0000 *Epitope
48S → T2.60.0203Epitope
69T → S/N/D/A/G2.20.0008DRAS
74L → V/I3.00.0000 *Epitope/DRAS
75V → I2.40.0024DRAS
98A → S2.00.0106Epitope
102K → R/Q/E/N/H20.10.0000 *DRAS
103K → N/S2.00.0003DRAS
118V → I1.70.0179Epitope
135I → T/V/L/R/M/K3.30.0000 *Epitope
162S → C/A/Y/D/N/H2.10.0003Epitope
184M → V2.50.0000 *Epitope/DRAS
188Y → L4.80.0002DRAS
200T → A/I/E17.40.0000 *Epitope
202I → V606.10.0083Epitope
211R → K/Q/G/T1.50.0051Epitope
215T → I3.40.0000 *Epitope
a Ratio of nonsynonymous (dN) to synonymous (dS) substitutions per site. b DRAS: Drug resistance associated sites according to the Stanford HIV drug resistance database [41,42] or the last Update of the Drug Resistance Mutations in HIV-1 of the International AIDS Society [43]. Epitope: peptide recognized by CD8+ T-cells according to the HIV Molecular Immunology Database of Los Alamos [40]. * p-values below the Bonferroni corrected (0.00016) significance level.
Eleven out of 18 selected codons in the PRO and 7 out of 17 in the RT have been previously associated with drug resistance; one of them (position 102 in the RT) was associated with drug resistance in a previous report, although it is not included in the last update of the Drug Resistance Mutations in HIV-1 or in the Stanford HIV drug resistance database [43,56,57]. On the other hand, 10 of the selected codons in the PRO and 12 in the RT were located in recognized CD8+ T-cell epitopes. Three codons in the PRO and two in the RT were located in drug resistance-associated sites that are also located within CD8+ T-cell epitopes. Similar results were obtained when the selection analysis was performed in the divided dataset using the SLAC method in Datamonkey; this method identified as positively selected the same positions previously identified by the HyPhy method implemented in MEGA 5.2.2, except for codon 202 of the RT.

3.2. Identification of CD8+ T-cell Epitopes with Amino Acid Substitutions

Based on the previous analysis, we selected the mutations located within CD8+ T-cell epitopes with significant (p < 0.05) evidence of positive selection and with a frequency greater than 5% among the 614 sequences. There were 20 such mutations, 11 in the PRO and 9 in the RT (Table 2).
Table
Table 2. Identification of epitopes affected by the positively selected mutations in in PRO and RT.
Table 2. Identification of epitopes affected by the positively selected mutations in in PRO and RT.
MutationsFrequency (%)Epitope Affected (HLA Alleles) aAssociation b
Protease
I13V20.4QRPLVTIKI (A*01:01)NC
QRPLVTIKIG (B51)NC
VTIKIGGQLK (A*11:01)SF
TIKIGGQLK (A3 supertype)NC
L19I9.0VTIKIGGQLK (A*11:01, A*03:01)SF
TIKIGGQLK (A3 supertype)NC
E35D28.7DTVLEEMSL (A*68:02)NC
EEMSLPGRW (B*44:02, B*44:03, B18, B40)IE
S37N56.0DTVLEEMSL (A*68:02)SF
EEMSLPGRW (B*44:02, B*44:03, B18, B40)SF
S37D14.1DTVLEEMSL (A*68:02)NC
EEMSLPGRW (B*44:02, B*44:03, B18, B40)SF
R41K42.8EEMSLPGRW (B*44:02, B*44:03, B18, B40)NC
LPGRWKPKMI (Cw3)NC
I54Vc19.3KMIGGIGGFI (B62)IE
I72V9.3IEICGHKAIG (B18, B40, B44)NC
GHKAIGTVL (B15)NC
I72T5.5IEICGHKAIG (B18, B40, B44)NC
GHKAIGTVL (B15)NC
V77Ic27.8LVGPTPVNI (A2)NC
V82A c13.8LVGPTPVNI (A2)IE
Reverse Transcriptase
T39A7.1ALVEICTEM (A*02, A*02:01, A2)NC
A98S8.7GIPHPAGLK (A*03:01)NC
V118I19.1VLDVGDAYFSV (A*02:01)NC
DAYFSVPL (A24, B*51:01)NC
I135T30.7KYTAFTIPSI (A2)NC
TAFTIPSI (B*51)IE
I35V7.7KYTAFTIPSI (A2)NC
TAFTIPSI (B*51)IE
S162C10.1SPAIFQSSM (B7, B35)SF
AIFQSSMTK (A*03:01)SF
T200A19.0DLEIGQHRTK (A3)NC
I202V6.8KIEELRQHL (A2)NC
KIEELRQHLL (B58)NC
IEELRQHLL (B*40:01, B60, B61)IE
R211K49.3EELRQHLLRW (B44)NC
a The residues in the epitope affected by the substitutions are in bold. b IE: Immune escape, i.e., mutations with documented or inferred escape; SF: susceptibility forms, defined as epitope sequences with mutations able to induce a CTL response; NC: none classified, i.e., mutations not included in the previous definitions. c Mutations associated with both, immune escape and drug resistance.
Some of these mutations were located in more than one overlapping epitope; for instance, the I13V mutation was found in four peptides, although only VTIKIGGQLK had been previously recognized as an epitope binding to different HLA molecules and influencing the immune response [58]. As a consequence 31 epitopes were affected by mutations at a high frequency. In addition, there was more than one mutation in the some epitopes, for example in DTVLEDMNL and EDMNLPGRW.
Only five mutant epitopes were previously recognized as immune escape (IE) and seven as susceptible forms (SF). Most (19) of the mutated epitopes have not been classified as IE or SF. All epitopes identified in this study as positively selected are presented by HLA alleles that have been reported in our population, except the epitope LPGRWKPKMI which is presented by HLA-Cw3; frequencies of alleles at the HLA-C locus have not been reported for Colombia.

3.3. Docking Simulation and Algorithmic Estimation of the Affinity of Peptides Binding to HLA Molecules

We implemented an in silico screening system to assess the affinity between the HLA proteins and the corresponding WT or mutated peptides, using a combination of different computational tools. This protocol was applied to peptides selected in Table 2. First, each peptide with the respective HLA molecule was subjected to docking analysis in AutoDockVina and the average of the data was plotted (Figure 1). All the evaluated peptides, either WT or mutated, were docked to the respective HLA molecule (Figure 1). Average free energy of binding ranged from −4.0 to −7.7 Kcal/mol for the PRO and RT peptides. No significant differences in values of free energy between the WT and the mutated forms of the epitopes where found, except for the DTVLEEMNL variant, which has a decrease of 0.44 Kcal/mol (p < 0.05) from the wild type DTVLEEMSL form (Figure 1A,B).
Viruses 07 01313 g001 1024
Figure 1. Free energy of binding to HLA molecules of peptides using the in silico assay. The mean of the lowest free energy of binding in 20 runs was plotted for each epitope derived from the PRO (A); and RT (B).
Figure 1. Free energy of binding to HLA molecules of peptides using the in silico assay. The mean of the lowest free energy of binding in 20 runs was plotted for each epitope derived from the PRO (A); and RT (B).
Viruses 07 01313 g001
Next, we estimated the binding affinity using three algorithmic tools: SMM, NetMHC and NetMHCpan. Many epitopes had an IC50 > 500 nM, indicating a low binding affinity; only one epitope in the PRO and six in the RT had good predicted binding affinities (IC50 < 50 nM) by at least one of the algorithms (Table 3). Four mutated epitopes in the PRO and three in the RT exhibited a significant decrease in the predicted affinity with respect to the WT form, as assessed by a two-fold or greater increase in the IC50 in at least two algorithms. Four of these epitopes, two in the PRO (TIKIGGQIK, DTVLEEMDL) and two in the RT (VLDVGDAYFSI, KYTAFTIPST) had not been reported previously as IE mutations (Table 3).
The IC50 estimated by the three algorithmic methods were highly correlated, with r ≥ 0.6991 and p < 0.0001. Conversely, the correlations of these values with the docking’s ΔGs were all negatively but not significantly correlated (Supplementary Table S1).
Table
Table 3. Epitope-HLA binding affinities by three computational algorithms (SMM, NetMHC and NetMHCpan).
Table 3. Epitope-HLA binding affinities by three computational algorithms (SMM, NetMHC and NetMHCpan).
Amino Acid SequenceAllelesSMMNetMHCNetMHCpan
Affinity (nM)Affinity (nM)Affinity (nM)
LPPVVAKEI aB*51172102797
NLVPMVATV aA*02662921
Protease
QRPLVTIKIA*01:011943342183736562
QRPLVTVKI2314992163937194
QRPLVTIKIGB511663603075143536
QRPLVTVKIG1682863070843446
TIKIGGQLKA3432582537
TVKIGGQLK456720856
TIKIGGQIK1041 b2126 b1654 b
DTVLEEMSLA*68:0287426861709
DTVLEEMNL101237952078
DTVLEEMDL2278 b12370 b6676 b
EEMSLPGRWB*44:02302514
EEMNLPGRW322821
EEMSLPGKW302214
EDMSLPGRW431 b565 b578 b
IEICGHKAIGB44209369696674
IEICGHKAVG212257925601
IEICGHKATG210350605162
GHKAIGTVLB15130911497215778
GHKAVGTVL98401365714015
GHKATGTVL82221294714673
LVGPTPVNIA2302738294005
LIGPTPVNI194521951596
LVGPTPANI355530143912
KMIGGIGGFIB62514415769
KMIGGIGGFV1493 b937 b1744 b
Reverse transcriptase
ALVEICTEMA21167041
ALVEICAEM1015032
GIPHPAGLKA*03:01290108316
GIPHPSGLK26699334
VLDVGDAYFSVA*02:0143142849
VLDVGDAYFSI8788 b53433 b
DAYFSVPLB*51:01762865274202
DAYFSIPL1515078054148
KYTAFTIPSIA2147051999216
KYTAFTIPST4925 b15627 b27884 b
KYTAFTIPSV38115095743
Reverse transcriptase
TAFTIPSIB5199623991153
TAFTIPST112816898 b17009 b
TAFTIPSV137240522887 b
DLEIGQHRTKA3798818815564
DLEIGQHRAK839901115773
KIEELRQHLA25689100728180
KVEELRQHL88531364414234
KIEELRQHLA25689100728180
KVEELRQHL 88531364414234
EELRQHLLRWB447810429
EELRQHLLKW788429
a Known epitopes of HIV-1 GAG and cytomegalovirus pp65, respectively, included as positive controls of the analyses performance. b Two-fold or larger differences between the affinities of the WT and the mutated peptide.

4. Discussion

In the present study, we explored immune escape mutations in the PRO and the RT, using HIV-1 sequences from 614 Colombian patients who were genotyped for antiviral sensitivity purposes. We found evidence of positive selection in many codons, associated or not with drug resistance (Table 1). All 17 selected substitutions not previously associated to drug resistance were located in CD8+ T-cell epitopes. Of these, 14 have not been previously identified as immune escape in the HIV Molecular Immunology Database of Los Alamos [40], opening the possibility of being considered as newly identified escape mutations, and prompting us to consider alternative explanations.
One possibility is that some of these mutations are indeed associated to drug resistance, but they have not been included in the drug resistance databases [43] due to the lack of strong evidence of this association. For instance, the I13V mutation in the PRO, that had been reported as a minor mutation associated with resistance to tipranavir and ritonavir, was removed from the 2010 update [59]. A related explanation is that these mutations could have compensatory effects that had not been previously detected. Since most HIV drug resistance mutations are naturally deleterious, such compensatory mutations are frequently selected. Similarly, compensatory mutations that restore fitness losses have previously been reported in CD8+ T-cell escape mutants [26,60].
An alternative explanation is that the selection test could be detecting significant elevations in the dN/dS ratio just for stochastic reasons; in a 304 codon-long sequence with a significance threshold of 0.05 one could expect 15 codons with a dN/dS < 0.05 just by chance. However, even with the strict Bonferroni correction for multiple comparisons, most (11/18 in the PRO and 7/17 in the RT) codons remain with a dN/dS ratio with a p-value below the corrected (0.00016) significance level in the test for positive selection, as shown in Table 1.
A more arguable explanation for this finding is the potential “hitch-hiking” effect of selective sweeps induced by the selective pressure exerted by antiviral drugs. In the presence of antiviral therapy, a drug resistance mutation will increase its frequency in the viral population up to fixation; other mutations arising in the same molecule, could hitch-hike along with it to fixation even if they are slightly deleterious [61]. Previous reports indicate that such hitch-hiking mutations persist for about 18 months in HIV-infected patients receiving antiviral therapy [62]. This effect is attenuated by the high rate of recombination, which is even higher than the mutation rate in terms of events per genome per HIV-1 cycle [63]; this could limit the effect of selective sweeps to closely linked sites. However, the extent of this effect is difficult to assess and we cannot rule out the possibility that this phenomenon has played a role in the mutations observed in our dataset.
Even so, the consistent location of non-synonymous substitutions in CD8+ T-cell epitopes, and not in other sites that could also be affected by hitch-hiking, still makes the immune escape hypothesis a plausible explanation for the high frequency of such mutations. Although not all mutations were located at the second position or at the C-terminal amino acid of the peptide chain that are the critical positions for the stability of the HLA-epitope interaction [64], the immune escape could also be associated with lack of T-cell receptor (TCR) recognition. Previous reports indicate that some mutations in the epitope can alter TCR recognition [65,66], resulting in a reduced functional capacity of CD8+ T-cells [67,68]. Functional studies are required to explore whether or not these mutations alter the strength of the interaction between the TCR and the HLA-peptide complex.
CD8+ T-cells play a crucial role in the control of HIV replication [69], although complete viral suppression is not achieved for a number of different reasons, including the appearance of mutations in immunodominant epitopes [70]. Several studies have demonstrated that HIV-specific CD8+ T-cell response selects escape mutants during the acute and chronic phases of HIV and simian immunodeficiency virus (SIV) infections [68,71,72], inducing loss of immune control and disease progression [73]. The selection of escape mutations depends on the HLA molecules. Indeed, the escape mutations could be predicted based on the host MHC alleles [74]. Most escape mutations revert when the virus is transmitted to another host with a different HLA make up. However, for HLA alleles that are common in a population, the corresponding escape mutations may become predominant in the strains circulating in that particular region, contributing to viral evolution and inducing an imprinting effect [25,75].
In this study, most CD8+ T-cell epitopes affected by positively selected mutations are restricted to the HLA alleles found in high frequency in the Colombian population (Supplementary Table S2). For instance, several mutated peptides are presented by the HLA-A*02 and -A*03 molecules, which in Colombia have frequencies of 22.2% and 7.9%, respectively [76]. In the case of HLA-B alleles, we observed mutations in epitopes restricted to the molecules -B*18 (7.1%), -B*35 (17.8%), -B*44 (9.1%) and -B*51 (6.4%) that are found in high frequencies (shown in parenthesis) [76]. However, the correlation between the frequency of the HLA alleles and the frequency of the mutation is not perfect. Only one HLA-A*24-restricted and one HLA-A*68-restricted epitopes showed substitutions in our study; these alleles had frequencies of 19.1% and 5.6%, respectively in Colombia; this could be due to a higher functional constraint against mutations in these peptides. Moreover, some epitopes restricted to rare HLA alleles in our population, such as HLA-B*15 (0.6%), were mutated with a significant (5.5%) frequency; these mutants could have appeared in a viral lineage before its introduction to Colombia and its frequency could have persisted or even increased by random genetic drift.
To predict the possible functional consequences of the observed amino acid substitutions in CD8+ T-cell epitopes, and further substantiate their immune escape potential, we performed an in silico docking simulation assay between the HLA molecules and the WT and mutated peptides. In addition, we predicted changes in the affinities of binding using three computational algorithms. We observed discrepancies between the inferred affinities obtained by docking simulation and by the computational algorithms; no correlation was observed in the values obtained by these two strategies. While the docking showed almost no changes in the free energy of binding in the mutant peptides, the algorithms predicted a significant decrease in affinities in seven peptides, four of which (TIKIGGQIK, DTVLEEMDL, VLDVGDAYFSI and KYTAFTIPST) had not been previously reported as immune escape mutants. Moreover, the results of the three algorithmic methods were strongly correlated (Supplementary Table S1).
The discrepancies between strategies could be explained by their different approaches. In general, the algorithms generate a prediction from a set of experimental data [77] derived from previous reports in-house experiments or specialized databases containing epitope-related information; estimates are expressed as a function of the concentration of a binding inhibitor [78]. The docking strategy is based on computational structural analysis; it predicts the affinity based on the spatial disposition of the atoms and chemical bonding, treating the protein as a rigid structure and estimating the free energy of binding; the accuracy of this kind of prediction could be affected by the size of the ligand, the dimension of the motion space or the rotational degrees of freedom of the ligand [79].
In this study, the algorithmic methods detected some (three out of five) of the previously reported IE mutations and allowed to infer four other mutations not previously reported as such. The docking strategy did not identify any of these IE mutations; instead, it detected a modest change of 0.44 kcal/mol between the WT and mutated forms in an epitope catalogued as SF. This suggests that algorithmic methods are more sensitive as tools for predicting potential IE mutants; however, we observed that some epitopes that experimentally had high binding to HLA molecules and high immunogenic capacity, exhibited low predicted affinities by the computational algorithms. For example, LVGPTPVNI was previously reported with a high HLA-A2 binding affinity and cytotoxic response of CD8+ T cells [80]; however, all the predicted affinities had IC50s > 500 nM (SMM: 3027.68 nM, NetMHC: 3829.00 nM and NetMHCpan: 4005.25 nM). Taking this into account, the four mutations found here with both, strong evidence of positive selection and significant decrease in the affinity of the epitopes, should be regarded only as candidate IE mutations; in vitro assays are required to determine whether they lead to destabilization of the HLA-peptide-TCR complex or to changes in the functional profiles in CD8+ T-cells.
The main limitations of this study were the absence of HLA typing, CD4+ T-cell counts and clinical data of the patients from whom the HIV-1 sequences were obtained, which would have allowed to make associations between mutations and disease progression. Nevertheless we can assume that all patients had viral loads ≥ 1000 RNA copies/µL, the minimal viral load that permits sequencing with the TRUGENE® kit. On the other hand, the strength of the bioinformatic approach to predict immune escape mutants that we introduce here is the combination of several analytical tools of different nature, and its application to a large dataset.
In summary, we introduced a strategy aimed at screening immune escape mutations that combine methods for positive selection, protein-protein docking and computational algorithms for affinity prediction. This allowed us to postulate four possible new mutations associated with immune escape in HIV-1-infected patients. Functional studies should be performed to confirm the effect of these mutations on the response of CD8+ T-cells; prospective clinical studies could test the biological relevance of these findings.

Supplementary Files

Supplementary File 1

Acknowledgments

This work was financed by Colciencias (grant 1115-569-33380) and the program “Estrategia de Sostenibilidad del Grupo Inmunovirología 2014–2015 de la Universidad de Antioquia”. L.A.-S. is recipient of a doctoral scholarship from Colciencias. We are indebted to two anonymous reviewers for many useful suggestions and to Anne-Lise Haenni for her grammar corrections.

Author Contributions

F.J.D., R.O., M.T.R., P.O.-G., P.A.V.-H. and L.A.-S. conceived and designed the experiment; F.J.D., R.O and L.A.-S. performed the analyses; F.J.D., R.O. and L.A.-S. analyzed the date; F.J.D. and L.A.-S. wrote the manuscript; F.J.D., M.T.R., R.O., P.A.V.-H. and L.A.-S. proof-read the manuscript. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References and Notes

  1. Hemelaar, J.; Gouws, E.; Ghys, P.D.; Osmanov, S. Global and regional distribution of HIV-1 genetic subtypes and recombinants in 2004. Aids 2006, 20, W13–W23. [Google Scholar] [CrossRef] [PubMed]
  2. Robertson, D.L.; Anderson, J.P.; Bradac, J.A.; Carr, J.K.; Foley, B.; Funkhouser, R.K.; Gao, F.; Hahn, B.H.; Kalish, M.L.; Kuiken, C.; et al. Hiv-1 nomenclature proposal. Science 2000, 288, 55–56. [Google Scholar]
  3. Ayouba, A.; Mauclere, P.; Martin, P.M.; Cunin, P.; Mfoupouendoun, J.; Njinku, B.; Souquieres, S.; Simon, F. HIV-1 group o infection in cameroon, 1986 to 1998. Emerg. Infect. Dis. 2001, 7, 466–467. [Google Scholar] [CrossRef] [PubMed]
  4. Ayouba, A.; Souquieres, S.; Njinku, B.; Martin, P.M.; Muller-Trutwin, M.C.; Roques, P.; Barre-Sinoussi, F.; Mauclere, P.; Simon, F.; Nerrienet, E. HIV-1 group n among HIV-1-seropositive individuals in cameroon. Aids 2000, 14, 2623–2625. [Google Scholar] [CrossRef] [PubMed]
  5. Plantier, J.C.; Leoz, M.; Dickerson, J.E.; De Oliveira, F.; Cordonnier, F.; Lemee, V.; Damond, F.; Robertson, D.L.; Simon, F. A new human immunodeficiency virus derived from gorillas. Nat. Med. 2009, 15, 871–872. [Google Scholar] [CrossRef] [PubMed]
  6. Robertson, D.L.; Sharp, P.M.; McCutchan, F.E.; Hahn, B.H. Recombination in hiv-1. Nature 1995, 374, 124–126. [Google Scholar] [CrossRef] [PubMed]
  7. Sanchez, G.I.; Bautista, C.T.; Eyzaguirre, L.; Carrion, G.; Arias, S.; Sateren, W.B.; Negrete, M.; Montano, S.M.; Sanchez, J.L.; Carr, J.K. Molecular epidemiology of human immunodeficiency virus-infected individuals in Medellin, Colombia. Am. J. Trop. Med. Hyg. 2006, 74, 674–677. [Google Scholar] [PubMed]
  8. Montano, S.M.; Sanchez, J.L.; Laguna-Torres, A.; Cuchi, P.; Avila, M.M.; Weissenbacher, M.; Serra, M.; Vinoles, J.; Russi, J.C.; Aguayo, N.; et al. Prevalences, genotypes, and risk factors for HIV transmission in south america. J. Acquir. Immune Defic. Syndr. 2005, 40, 57–64. [Google Scholar] [CrossRef]
  9. Ho, D.D.; Neumann, A.U.; Perelson, A.S.; Chen, W.; Leonard, J.M.; Markowitz, M. Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection. Nature 1995, 373, 123–126. [Google Scholar] [CrossRef] [PubMed]
  10. O'Neil, P.K.; Sun, G.; Yu, H.; Ron, Y.; Dougherty, J.P.; Preston, B.D. Mutational analysis of HIV-1 long terminal repeats to explore the relative contribution of reverse transcriptase and RNA polymerase II to viral mutagenesis. J. Biol. Chem. 2002, 277, 38053–38061. [Google Scholar] [CrossRef] [PubMed]
  11. Ji, J.P.; Loeb, L.A. Fidelity of HIV-1 reverse transcriptase copying RNA in vitro. Biochemistry 1992, 31, 954–958. [Google Scholar] [CrossRef] [PubMed]
  12. Abram, M.E.; Ferris, A.L.; Shao, W.; Alvord, W.G.; Hughes, S.H. Nature, position, and frequency of mutations made in a single cycle of HIV-1 replication. J. Virol. 2010, 84, 9864–9878. [Google Scholar] [CrossRef] [PubMed]
  13. Mansky, L.M.; Temin, H.M. Lower in vivo mutation rate of human immunodeficiency virus type 1 than that predicted from the fidelity of purified reverse transcriptase. J. Virol. 1995, 69, 5087–5094. [Google Scholar] [PubMed]
  14. Simon, V.; Ho, D.D. HIV-1 dynamics in vivo: Implications for therapy. Nat. Rev. Microbiol. 2003, 1, 181–190. [Google Scholar] [CrossRef] [PubMed]
  15. Leslie, A.J.; Pfafferott, K.J.; Chetty, P.; Draenert, R.; Addo, M.M.; Feeney, M.; Tang, Y.; Holmes, E.C.; Allen, T.; Prado, J.G.; et al. HIV evolution: CTL escape mutation and reversion after transmission. Nat. Med. 2004, 10, 282–289. [Google Scholar]
  16. Marozsan, A.J.; Fraundorf, E.; Abraha, A.; Baird, H.; Moore, D.; Troyer, R.; Nankja, I.; Arts, E.J. Relationships between infectious titer, capsid protein levels, and reverse transcriptase activities of diverse human immunodeficiency virus type 1 isolates. J. Virol. 2004, 78, 11130–11141. [Google Scholar] [CrossRef] [PubMed]
  17. Dimitrov, D.S.; Willey, R.L.; Sato, H.; Chang, L.J.; Blumenthal, R.; Martin, M.A. Quantitation of human immunodeficiency virus type 1 infection kinetics. J. Virol. 1993, 67, 2182–2190. [Google Scholar] [PubMed]
  18. Hemelaar, J. Implications of HIV diversity for the HIV-1 pandemic. J. Infect. 2013, 66, 391–400. [Google Scholar] [CrossRef] [PubMed]
  19. De, S.L.E.; Holmes, E.C.; Zanotto, P.M. Distinct patterns of natural selection in the reverse transcriptase gene of HIV-1 in the presence and absence of antiretroviral therapy. Virology 2004, 325, 181–191. [Google Scholar] [CrossRef] [PubMed]
  20. McMichael, A.J.; Borrow, P.; Tomaras, G.D.; Goonetilleke, N.; Haynes, B.F. The immune response during acute HIV-1 infection: Clues for vaccine development. Nat. Rev. Immunol. 2010, 10, 11–23. [Google Scholar] [CrossRef] [PubMed]
  21. Kawashima, Y.; Pfafferott, K.; Frater, J.; Matthews, P.; Payne, R.; Addo, M.; Gatanaga, H.; Fujiwara, M.; Hachiya, A.; Koizumi, H.; et al. Adaptation of HIV-1 to human leukocyte antigen class I. Nature 2009, 458, 641–645. [Google Scholar]
  22. Moore, C.B.; John, M.; James, I.R.; Christiansen, F.T.; Witt, C.S.; Mallal, S.A. Evidence of HIV-1 adaptation to HLA-restricted immune responses at a population level. Science 2002, 296, 1439–1443. [Google Scholar] [CrossRef] [PubMed]
  23. Matthews, P.C.; Leslie, A.J.; Katzourakis, A.; Crawford, H.; Payne, R.; Prendergast, A.; Power, K.; Kelleher, A.D.; Klenerman, P.; Carlson, J.; et al. HLA footprints on human immunodeficiency virus type 1 are associated with interclade polymorphisms and intraclade phylogenetic clustering. J. Virol. 2009, 83, 4605–4615. [Google Scholar]
  24. Hemelaar, J. The origin and diversity of the HIV-1 pandemic. Trends Mol. Med. 2012, 18, 182–192. [Google Scholar] [CrossRef] [PubMed]
  25. Brumme, Z.L.; Brumme, C.J.; Heckerman, D.; Korber, B.T.; Daniels, M.; Carlson, J.; Kadie, C.; Bhattacharya, T.; Chui, C.; Szinger, J.; et al. Evidence of differential HLA class I-mediated viral evolution in functional and accessory/regulatory genes of HIV-1. PLoS Pathog. 2007, 3, e94. [Google Scholar]
  26. Kelleher, A.D.; Long, C.; Holmes, E.C.; Allen, R.L.; Wilson, J.; Conlon, C.; Workman, C.; Shaunak, S.; Olson, K.; Goulder, P.; et al. Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses. J. Exp. Med. 2001, 193, 375–386. [Google Scholar]
  27. Feeney, M.E.; Tang, Y.; Roosevelt, K.A.; Leslie, A.J.; McIntosh, K.; Karthas, N.; Walker, B.D.; Goulder, P.J. Immune escape precedes breakthrough human immunodeficiency virus type 1 viremia and broadening of the cytotoxic T-lymphocyte response in an HLA-B27-positive long-term-nonprogressing child. J. Virol. 2004, 78, 8927–8930. [Google Scholar] [CrossRef] [PubMed]
  28. Kawashima, Y.; Kuse, N.; Gatanaga, H.; Naruto, T.; Fujiwara, M.; Dohki, S.; Akahoshi, T.; Maenaka, K.; Goulder, P.; Oka, S.; et al. Long-term control of HIV-1 in hemophiliacs carrying slow-progressing allele HLA-B*5101. J. Virol. 2010, 84, 7151–7160. [Google Scholar]
  29. Grant, R.M.; Kuritzkes, D.R.; Johnson, V.A.; Mellors, J.W.; Sullivan, J.L.; Swanstrom, R.; D’Aquila, R.T.; van Gorder, M.; Holodniy, M.; Lloyd, R.M., Jr.; et al. Accuracy of the trugene HIV-1 genotyping kit. J. Clin. Microbiol. 2003, 41, 1586–1593. [Google Scholar]
  30. Gómez, S.; Olaya-García, P.; Díaz, F. Resistencia a los medicamentos antirretrovirales en pacientes que reciben tratamiento para vih-sida en colombia. Infectio 2010, 14, 248–256. [Google Scholar] [CrossRef]
  31. Pineda-Peña, A.; Rodrigues Faria, N.; Diaz, F.; Olaya-García, P.; Moller Frederiksen, C.; Guangdi, L.; Gomez-Lopez, A.; Lemey, P.; Vandamme, A. The colombian pidemic is dominated by HIV-1subtype b over time: A molecular epidemiolog and phylodynamics study. In Proceedings of the 20th International HIV Dynamics and Evolution, Utrecht, The Netherlands, 8–11 May 2013.
  32. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. Mega5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed]
  33. Nei, M.; Gojobori, T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol. Biol. Evol. 1986, 3, 418–426. [Google Scholar] [PubMed]
  34. Pond, S.L.; Frost, S.D.; Muse, S.V. Hyphy: Hypothesis testing using phylogenies. Bioinformatics 2005, 21, 676–679. [Google Scholar] [CrossRef] [PubMed]
  35. Muse, S.V.; Gaut, B.S. A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome. Mol. Biol. Evol. 1994, 11, 715–724. [Google Scholar] [PubMed]
  36. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University Press, Inc: New York, NY, USA, 2000; pp. 147–163. [Google Scholar]
  37. Kosakovsky Pond, S.L.; Frost, S.D. Not so different after all: A comparison of methods for detecting amino acid sites under selection. Mol. Biol. Evol. 2005, 22, 1208–1222. [Google Scholar] [CrossRef]
  38. Suzuki, Y.; Gojobori, T. A method for detecting positive selection at single amino acid sites. Mol. Biol. Evol. 1999, 16, 1315–1328. [Google Scholar] [CrossRef] [PubMed]
  39. Delport, W.; Poon, A.F.; Frost, S.D.; Kosakovsky Pond, S.L. Datamonkey 2010: A suite of phylogenetic analysis tools for evolutionary biology. Bioinformatics 2010, 26, 2455–2457. [Google Scholar] [CrossRef] [PubMed]
  40. HIV molecular immunology database. Available online: http://www.hiv.lanl.gov/ (accessed on 12 February 2014).
  41. Liu, T.F.; Shafer, R.W. Web Resources for HIV type 1 genotypic-resistance Test Interpretation. Clin. Infect. Dis. 2006, 42, 1608–1618. [Google Scholar] [CrossRef] [PubMed]
  42. HIV drug resistance database – Stanford University. Available online: http://hivdb.stanford.edu/ (accessed on 16 May 2014).
  43. Johnson, V.A.; Calvez, V.; Gunthard, H.F.; Paredes, R.; Pillay, D.; Shafer, R.W.; Wensing, A.M.; Richman, D.D. Update of the drug resistance mutations in HIV-1: March 2013. Top. Antivir. Med. 2013, 21, 6–14. [Google Scholar] [PubMed]
  44. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The protein data bank. Nucleic. Acids. Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed]
  45. Morris, G.M.; Green, L.G.; Radic, Z.; Taylor, P.; Sharpless, K.B.; Olson, A.J.; Grynszpan, F. Automated docking with protein flexibility in the design of femtomolar “click chemistry” inhibitors of acetylcholinesterase. J. Chem. Inf. Model. 2013, 53, 898–906. [Google Scholar] [CrossRef] [PubMed]
  46. Zhang, Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9, 40. [Google Scholar] [CrossRef]
  47. Roy, A.; Kucukural, A.; Zhang, Y. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar] [CrossRef] [PubMed]
  48. Lambert, C.; Leonard, N.; de Bolle, X.; Depiereux, E. ESyPred3D: Prediction of proteins 3D structures. Bioinformatics 2002, 18, 1250–1256. [Google Scholar] [CrossRef] [PubMed]
  49. Trott, O.; Olson, A.J. Autodock vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [PubMed]
  50. Hoof, I.; Peters, B.; Sidney, J.; Pedersen, L.E.; Sette, A.; Lund, O.; Buus, S.; Nielsen, M. Netmhcpan, a method for mhc class I binding prediction beyond humans. Immunogenetics 2009, 61, 1–13. [Google Scholar] [CrossRef] [PubMed]
  51. Nielsen, M.; Lundegaard, C.; Blicher, T.; Lamberth, K.; Harndahl, M.; Justesen, S.; Roder, G.; Peters, B.; Sette, A.; Lund, O.; et al. Netmhcpan, a method for quantitative predictions of peptide binding to any HLA-A and -B locus protein of known sequence. PLoS One 2007, 2, e796. [Google Scholar]
  52. Buus, S.; Lauemoller, S.L.; Worning, P.; Kesmir, C.; Frimurer, T.; Corbet, S.; Fomsgaard, A.; Hilden, J.; Holm, A.; Brunak, S. Sensitive quantitative predictions of peptide-MHC binding by a “query by committee” artificial neural network approach. Tissue Antigens 2003, 62, 378–384. [Google Scholar] [CrossRef] [PubMed]
  53. Peters, B.; Sette, A. Generating quantitative models describing the sequence specificity of biological processes with the stabilized matrix method. BMC Bioinform. 2005, 6, 132. [Google Scholar] [CrossRef] [Green Version]
  54. Immune epitope database and analysis resource. Available online: http://www.iedb.org/ (accessed on 12 January 2015).
  55. Sette, A.; Vitiello, A.; Reherman, B.; Fowler, P.; Nayersina, R.; Kast, W.M.; Melief, C.J.; Oseroff, C.; Yuan, L.; Ruppert, J.; et al. The relationship between class I binding affinity and immunogenicity of potential cytotoxic T cell epitopes. J. Immunol. 1994, 153, 5586–5592. [Google Scholar]
  56. Hachiya, A.; Gatanaga, H.; Kodama, E.; Ikeuchi, M.; Matsuoka, M.; Harada, S.; Mitsuya, H.; Kimura, S.; Oka, S. Novel patterns of nevirapine resistance-associated mutations of human immunodeficiency virus type 1 in treatment-naive patients. Virology 2004, 327, 215–224. [Google Scholar] [CrossRef] [PubMed]
  57. Alteri, C.; Santoro, M.M.; Abbate, I.; Rozera, G.; Bruselles, A.; Bartolini, B.; Gori, C.; Forbici, F.; Orchi, N.; Tozzi, V.; et al. “Sentinel” mutations in standard population sequencing can predict the presence of HIV-1 reverse transcriptase major mutations detectable only by ultra-deep pyrosequencing. J. Antimicrob. Chemother. 2011, 66, 2615–2623. [Google Scholar]
  58. Sidney, J.; Southwood, S.; Sette, A. HLA Binding Peptides and Their Uses. U.S. Patent 20,060,079,453, 13 April 2006. [Google Scholar]
  59. Manosuthi, W.; Thongyen, S.; Nilkamhang, S.; Manosuthi, S.; Sungkanuparph, S. HIV-1 drug resistance-associated mutations among antiretroviral-naive thai patients with chronic HIV-1 infection. J. Med. Virol. 2013, 85, 194–199. [Google Scholar] [CrossRef] [PubMed]
  60. Brockman, M.A.; Schneidewind, A.; Lahaie, M.; Schmidt, A.; Miura, T.; Desouza, I.; Ryvkin, F.; Derdeyn, C.A.; Allen, S.; Hunter, E.; et al. Escape and compensation from early HLA-B57-mediated cytotoxic T-lymphocyte pressure on human immunodeficiency virus type 1 gag alter capsid interactions with cyclophilin a. J. Virol. 2007, 81, 12608–12618. [Google Scholar]
  61. Smith, J.M.; Haigh, J. The hitch-hiking effect of a favourable gene. Genet. Res. 2007, 89, 391–403. [Google Scholar] [CrossRef] [PubMed]
  62. Pennings, P.S.; Kryazhimskiy, S.; Wakeley, J. Loss and recovery of genetic diversity in adapting populations of HIV. PLoS Genet. 2014, 10, e1004000. [Google Scholar] [CrossRef] [PubMed]
  63. Neher, R.A.; Leitner, T. Recombination rate and selection strength in HIV intra-patient evolution. PLoS Comput. Biol. 2010, 6, e1000660. [Google Scholar] [CrossRef] [PubMed]
  64. Madden, D.R. The three-dimensional structure of peptide-MHC complexes. Annu. Rev. Immunol. 1995, 13, 587–622. [Google Scholar] [CrossRef] [PubMed]
  65. Ladell, K.; Hashimoto, M.; Iglesias, M.C.; Wilmann, P.G.; McLaren, J.E.; Gras, S.; Chikata, T.; Kuse, N.; Fastenackels, S.; Gostick, E.; et al. A molecular basis for the control of preimmune escape variants by HIV-specific CD8+ T cells. Immunity 2013, 38, 425–436. [Google Scholar]
  66. Reid, S.W.; McAdam, S.; Smith, K.J.; Klenerman, P.; O'Callaghan, C.A.; Harlos, K.; Jakobsen, B.K.; McMichael, A.J.; Bell, J.I.; Stuart, D.I.; et al. Antagonist HIV-1 Gag peptides induce structural changes in HLA B8. J. Exp. Med. 1996, 184, 2279–2286. [Google Scholar]
  67. Davenport, M.P. Antagonists or altruists: Do viral mutants modulate T-cell responses? Immunol. Today 1995, 16, 432–436. [Google Scholar] [CrossRef] [PubMed]
  68. Allen, T.M.; O’Connor, D.H.; Jing, P.; Dzuris, J.L.; Mothe, B.R.; Vogel, T.U.; Dunphy, E.; Liebl, M.E.; Emerson, C.; Wilson, N.; et al. Tat-specific cytotoxic T lymphocytes select for siv escape variants during resolution of primary viraemia. Nature 2000, 407, 386–390. [Google Scholar]
  69. Koup, R.A.; Safrit, J.T.; Cao, Y.; Andrews, C.A.; McLeod, G.; Borkowsky, W.; Farthing, C.; Ho, D.D. Temporal association of cellular immune responses with the initial control of viremia in primary human immunodeficiency virus type 1 syndrome. J. Virol. 1994, 68, 4650–4655. [Google Scholar] [PubMed]
  70. Borrow, P.; Lewicki, H.; Wei, X.; Horwitz, M.S.; Peffer, N.; Meyers, H.; Nelson, J.A.; Gairin, J.E.; Hahn, B.H.; Oldstone, M.B.; et al. Antiviral pressure exerted by HIV-1-specific cytotoxic T lymphocytes (CTLs) during primary infection demonstrated by rapid selection of CTL escape virus. Nat. Med. 1997, 3, 205–211. [Google Scholar]
  71. Brumme, Z.L.; John, M.; Carlson, J.M.; Brumme, C.J.; Chan, D.; Brockman, M.A.; Swenson, L.C.; Tao, I.; Szeto, S.; Rosato, P.; et al. HLA-associated immune escape pathways in HIV-1 subtype B Gag, Pol and Nef proteins. PLoS One 2009, 4, e6687. [Google Scholar]
  72. Prado, J.G.; Honeyborne, I.; Brierley, I.; Puertas, M.C.; Martinez-Picado, J.; Goulder, P.J. Functional consequences of human immunodeficiency virus escape from an HLA-B*13-restricted CD8+ T-cell epitope in p1 Gag protein. J. Virol. 2009, 83, 1018–1025. [Google Scholar] [CrossRef] [PubMed]
  73. Walker, B.D.; Korber, B.T. Immune control of HIV: The obstacles of HLA and viral diversity. Nat. Immunol. 2001, 2, 473–475. [Google Scholar] [CrossRef]
  74. Carlson, J.M.; Brumme, Z.L. HIV evolution in response to HLA-restricted CTL selection pressures: A population-based perspective. Microbes Infect./Inst. Pasteur 2008, 10, 455–461. [Google Scholar] [CrossRef]
  75. Goulder, P.J.; Brander, C.; Tang, Y.; Tremblay, C.; Colbert, R.A.; Addo, M.M.; Rosenberg, E.S.; Nguyen, T.; Allen, R.; Trocha, A.; et al. Evolution and transmission of stable CTL escape mutations in HIV infection. Nature 2001, 412, 334–338. [Google Scholar]
  76. Rodriguez, L.M.; Giraldo, M.C.; Garcia, N.; Velasquez, L.; Paris, S.C.; Alvarez, C.M.; Garcia, L.F. Human leucocyte antigen gene (HLA-A, HLA-B, HLA-DRB1) frequencies in deceased organ donors. Biomedica 2007, 27, 537–547. [Google Scholar] [CrossRef] [PubMed]
  77. Kim, Y.; Ponomarenko, J.; Zhu, Z.; Tamang, D.; Wang, P.; Greenbaum, J.; Lundegaard, C.; Sette, A.; Lund, O.; Bourne, P.E.; et al. Immune epitope database analysis resource. Nucleic Acids Res. 2012, 40, W525–W530. [Google Scholar]
  78. Peters, B.; Bui, H.H.; Frankild, S.; Nielson, M.; Lundegaard, C.; Kostem, E.; Basch, D.; Lamberth, K.; Harndahl, M.; Fleri, W.; et al. A community resource benchmarking predictions of peptide binding to MHC-I molecules. PLoS Comput. Biol. 2006, 2, e65. [Google Scholar]
  79. Dhanik, A.; McMurray, J.S.; Kavraki, L.E. Dinc: A new autodock-based protocol for docking large ligands. BMC Struct. Biol. 2013, 13, S11. [Google Scholar] [CrossRef] [PubMed]
  80. Konya, J.; Stuber, G.; Bjorndal, A.; Fenyo, E.M.; Dillner, J. Primary induction of human cytotoxic lymphocytes against a synthetic peptide of the human immunodeficiency virus type 1 protease. J. Gen. Virol. 1997, 78, 2217–2224. [Google Scholar] [PubMed]

Share and Cite

MDPI and ACS Style

Acevedo-Sáenz, L.; Ochoa, R.; Rugeles, M.T.; Olaya-García, P.; Velilla-Hernández, P.A.; Diaz, F.J. Selection Pressure in CD8+ T-cell Epitopes in the pol Gene of HIV-1 Infected Individuals in Colombia. A Bioinformatic Approach. Viruses 2015, 7, 1313-1331. https://doi.org/10.3390/v7031313

AMA Style

Acevedo-Sáenz L, Ochoa R, Rugeles MT, Olaya-García P, Velilla-Hernández PA, Diaz FJ. Selection Pressure in CD8+ T-cell Epitopes in the pol Gene of HIV-1 Infected Individuals in Colombia. A Bioinformatic Approach. Viruses. 2015; 7(3):1313-1331. https://doi.org/10.3390/v7031313

Chicago/Turabian Style

Acevedo-Sáenz, Liliana, Rodrigo Ochoa, Maria Teresa Rugeles, Patricia Olaya-García, Paula Andrea Velilla-Hernández, and Francisco J. Diaz. 2015. "Selection Pressure in CD8+ T-cell Epitopes in the pol Gene of HIV-1 Infected Individuals in Colombia. A Bioinformatic Approach" Viruses 7, no. 3: 1313-1331. https://doi.org/10.3390/v7031313

APA Style

Acevedo-Sáenz, L., Ochoa, R., Rugeles, M. T., Olaya-García, P., Velilla-Hernández, P. A., & Diaz, F. J. (2015). Selection Pressure in CD8+ T-cell Epitopes in the pol Gene of HIV-1 Infected Individuals in Colombia. A Bioinformatic Approach. Viruses, 7(3), 1313-1331. https://doi.org/10.3390/v7031313

Article Metrics

Back to TopTop