Skip to main content
  • Research article
  • Open access
  • Published:

Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis



The correlations of genotypic and phenotypic tests with treatment, clinical history and the significance of mutations in viruses of HIV-infected patients are used to establish resistance mutations to protease inhibitors (PIs). Emerging mutations in human immunodeficiency virus type 1 (HIV-1) protease confer resistance to PIs by inducing structural changes at the ligand interaction site. The aim of this study was to establish an in silico structural relationship between natural HIV-1 polymorphisms and unusual HIV-1 mutations that confer resistance to PIs.


Protease sequences isolated from 151 Mexican HIV-1 patients that were naïve to, or subjected to antiretroviral therapy, were examined. We identified 41 unrelated resistance mutations with a prevalence greater than 1%. Among these mutations, nine exhibited positive selection, three were natural polymorphisms (L63S/V/H) in a codon associated with drug resistance, and six were unusual mutations (L5F, D29V, L63R/G, P79L and T91V). The D29V mutation, with a prevalence of 1.32% in the studied population, was only found in patients treated with antiretroviral drugs. Using in silico modelling, we observed that D29V formed unstable protease complexes when were docked with lopinavir, saquinavir, darunavir, tipranavir, indinavir and atazanavir.


The structural correlation of natural polymorphisms and unusual mutations with drug resistance is useful for the identification of HIV-1 variants with potential resistance to PIs. The D29V mutation likely confers a selection advantage in viruses; however, in silico, presence of this mutation results in unstable enzyme/PI complexes, that possibly induce resistance to PIs.


Diversity of viral populations is a result of sophisticated recombination, replication and/or selection events that induce drug-resistant human immunodeficiency virus type 1 (HIV-1) variants. The lack of reverse transcription corrections, transitional printing and transversion mutations, along with viral recombination, has resulted in the emergence of HIV-1 variants with high resistance to pharmacological stressors [1, 2]. These variants form populations that evade antiretroviral agents, due to emerging phenotypic changes within and around the active enzyme site [3]. These mutations, which give rise to drug resistance, result in reduced efficacy of highly active antiretroviral therapy (HAART) [4]. Correlations between genotypic and phenotypic tests with treatment, clinical history, and significance of mutations identified in HIV-1 of infected patients are used to determine the presence of mutations that confer resistance to protease inhibitors (PIs) [1].

Disruption at interaction sites causes an alteration in affinity between proteins and their inhibitors, and has been recognized as a property of drug resistant HIV-1 proteins [5, 6]. Protein folding simulation models can create Local Elementary Structures (LES). These secondary structures are stabilized by amino acids that interact with the polypeptide chain [7]. Using the Gromacs software (version 3.0), LES were found to form in protease (PR) regions 23–33, 74–78, and 83–92, and also docked in a folding nucleus [8]. Other studies have shown that mutations further from the active site can alter the flexibility of HIV-1 PR, inducing structural changes that affect the efficacy of most PIs currently used [9]. Theoretical studies, either alone or in combination with experimental methods, have pointed to an increase in the flexibility of mutant enzymes at various sites, including the active site, as a resistance mechanism that causes a decrease in the affinity of PIs [10]. Part of the cause of such flexibility could be the unusual mutations that generally emerge only after "major" and "minor" resistance mutations have been introduced [11]. Other mutations that can affect the interaction between PR and PIs are natural polymorphisms and unusual mutations in positions that confer drug resistance. Although the main mutations associated with drug resistance have been characterized [12, 13], little is known about the influence of natural polymorphisms and unusual mutations with respect to the development of drug resistance. The aim of this study was to describe an in silico experiment that showed structural correlations between natural HIV-1 polymorphisms and unusual HIV-1mutations in the PR region of HIV-1 pol with potential PIs resistance.


Sequence data

We analysed 151 HIV-1 sequences from Mexican patients who had been tested for resistance to antiretroviral drugs between 2005 and 2011 in the Laboratory of Immunodeficiencies and Human Retroviruses, Western Biomedical Research Center, Mexican Institute of Social Security. Sequences were obtained from 22 naïve, and 129 treated patients that were not responsive to drugs. Sequences were registered in the GenBank database [14], with the following accession numbers: [EU045452–EU045489; GU382757–GU382851; GU437199–GU437200; and KC416212–KC416227]. All sequences were analysed for the presence or absence of highly mutated sequences using HYPERMUT software (version 2.0) [15]. For a reference sequence, we used the subtype B consensus sequence, which was derived from an alignment of subtype B sequences maintained at the Los Alamos HIV Sequence Database (LANL), and available from the HIV Drug Resistance Database (HIVDB), Stanford University [16].

Phylogenetic analysis

Nucleotide homology analysis for HIV-1 sequences was conducted using the NCBI Genotyping Tool program [17]. Subtype determinations were further confirmed by phylogenetic analysis performed with the Molecular Evolution Genetics Analysis (MEGA) software package (version 5.0) [18], which includes the recommended reference sequence sets, available from the Los Alamos HIV Sequence Database [19]. Prior to all phylogenetic analyses, HIV-1 pol sequences were aligned using Clustal X (European Bioinformatics Institute, EMBL) [20]. Sequences with 100% homology were excluded from the analysis. The nucleotide distance matrix was generated using the Kimura two-parameter Neighbour-joining method [21]. The statistical robustness of the generated trees was verified by bootstrap analysis of 1000 replicates.

Detection of multidrug resistance phenotypes in HIV-1 protease

The genetic changes associated with drug resistance in viral sequences were established according to HIVdb algorithm version 6.0.9 ( [22]. The interpretation of drug resistance was performed at various levels of susceptibility for the following USA Food and Drug Administration (FDA)-approved PIs: atazanavir (ATV); darunavir (DRV); amprenavir (APV); indinavir (IDV); lopinavir (LPV); saquinavir (SQV); tipranavir (TPV); nelfinavir (NFV);and ritonavir (RTV). The resistance mutations were classified as major or minor according to HIVdb criteria, or as natural polymorphisms or unusual mutations if they were not associated with resistance [16]. The prevalence (p) for each mutation in the protease region of pol was quantitatively determined as the frequency of the mutation (M) among total sequences evaluated for each position (N), p = M/N, using Microsoft Excel 2010. The genetic variation was calculated as the total number of mutations at a nucleotide position divided by the number of evaluated sequences. The Phenotypic Variation (PV) was defined as the percentage (%) of amino acid substitutions for each position relative to the consensus sequence. For each region, the PV was classified as follows: conserved, <1%; semi-conserved, 1 to <5%; variable, 5 to <10%; and highly variable, ≥10%. Values found below the 15th percentile and above the 75th percentile were not considered. Phenotypic mutations with a prevalence of ≥1.0% among 151 amino acid sequences were compared for each PI against the IAS–USA drug resistance mutations list [12]. The structural conservation of PR was defined in a complementary way to that of PV.

Analysis of selective pressure

The selective pressure and reconstruction of the ancestral state for each PR codon was determined using a maximum likelihood (ML) substitution model and the HyPhy algorithm, included in MEGA5 package [23, 24]. The synonymous site divergence (dS) and nonsynonymous site divergence (dN) per branch was estimated using the Muse–Gaut codon model [25]. The values of the ML model were estimated from the topology of the phylogenetic tree. The probability of rejecting the hypothesis of neutral evolution was significant with p ≤ 0.05. The standardized values of dN–dS were obtained by the total number of substitutions in the tree (calculated as substitutions expected by site). To distinguish between drug pressure and immune system pressure, results were compared using the HIV positive selection mutation database [26].

Molecular modelling

Once the natural polymorphisms and unusual mutation codons with positive selection (dN/dS > 1) and prevalence >1% were obtained, homology modelling was used to predict changes in the PR structure. Homology modelling of natural polymorphisms and unusual mutations followed these steps: (i) template selection; (ii) structural alignment; (iii) model construction; and (iv) refinement [27]. To select the template, HIV-1 protease X-ray crystal structure FASTA sequences available from the Protein Data Bank (PDB) [28] and the HIV-1 subtype B consensus sequence available from the HIVdb were aligned using ClustalX [20]. The PR sequence exhibiting greatest identity with HIV-1 subtype B consensus (wild-type template) was chosen as the template for modelling mutant proteases (PRs). The resistant PRs used for reference were modelled with each major PI resistance mutation. Every mutant protein was modelled using a Swiss-Model workspace, which showed the identity (%). The expected alignment value with the template sequence (E) and the Qualitative Model Energy Analysis (QMEAN4), which estimates the absolute quality model, ranged from 0–1 [28, 29].

Estimation of the free energy of binding

Using the Autodock/Vina application on a LINUX platform, which had the PyMOL (version 1.4.1) molecular graphics system installed, we estimated the free energy of binding of the complex between mutant PR structures and Pls [30]. Rectangular boxes were used to define the binding sites and these were adjusted by providing specific coordinates of active PR sites before each docking.

Receptor and ligand representations in the Protein Data Bank, Partial Charge & Atom Type formats (pdbqt) containing atomic charges, atom type definitions and topological information, were produced using Autodock/Vina [30]. To determine if the differences caused by natural polymorphisms and unusual mutations had any effect on the free energy of binding of Pls, the free energy values obtained for the resistant protease/ligand complexes were compared. Natural polymorphisms or unusual mutations with lower or equal affinity to PIs compared with reference proteins containing drug resistance mutations indicate positive resistance. Higher affinity was considered to favour susceptibility of the HIV-1 variant to PIs. The coupled proteases included the wild-type PR [PDB: 1GNO], PRs with major drug-resistance mutations and PRs with natural polymorphisms or unusual mutations at codons having positive selection.

Measurement of distances between protease residues and PIs

To evaluate the natural polymorphism and unusual mutation atoms that affect the affinity of PIs, we measured the distances (Å) between the amino acid residue Cα-atoms implicated in drug resistance, and the closest heteroatoms of the PIs. Complexes that showed signs of free energy of binding were analysed, suggesting increased drug resistance because of the presence of natural polymorphisms and unusual mutations. Distances were compared with those obtained for the same pair of atoms in the wild-type and resistant PR structures available from the PDB [28]. All interatomic distances were measured with PyMOL (version 1.4.1) [31].

Results and discussion

Genetic relationships of HIV-1 variants isolated from Mexican patients

Phylogenetic analysis of the 151 HIV-1 protease fragment nucleotide sequences was conducted using a Neighbour-joining tree. Phylogenetic relationships were grouped into the internal nodes of the tree, using subtype B reference sequences [GenBank: U63632 and U21135]. The HIV-1 variants isolated from Mexican patients, and confirmed by analysis with the NCBI Genotyping Tool, were subtype B. This result is consistent with other molecular epidemiology studies of Mexican HIV-1 patients, with or without antiretroviral intervention, where subtype B prevails [32, 33].

Drug resistant phenotypes and genotypes of HIV-1 protease

Non-synonymous genetic changes largely contribute to phenotypic changes [34]. Because of degeneracy in the genetic code, transcription and translation errors during the viral replication cycle, along with functional, structural, pharmacological and immunological selection pressures, there is no absolute mathematical relationship between genetic and phenotypic variations [35]. Variations in the primary structure of the 151 PR sequences are presented in Figure 1. Of the 36 codons associated with major or minor resistance, 19 showed PV <5% (L11, L24, D30, V32, L33, E34, K43, I47, G48, I50, F53, D60, G73, T74, L76, N83, I85, N88 and L89). Among these codons, G48, I50 and F53 were present in conserved regions, and L24, D30, L76 and N88 in semi-conserved regions. Six codons (G16, K20, I54, Q58, H69 and I84) displayed PV of 5–10%. Of these, only codons I54 and Q58 were located in conserved and semi-conserved regions, respectively. The remaining 11 codons associated with drug resistance (L10, M36, M46, I62, L63, I64, A71, V77, V82, L90 and I93) had variation >10%, with L10 and I93 present in conserved regions, and M46, V77 and L90 in semi-conserved regions. There was also a variable PV of 5–10%, in the codons neighbouring the drug resistance positions (T12, I15, L19, G68, and K70), with codons E35, N37, R41, R57, and I72 highly variable. Figure 2 illustrates the mutations with prevalence ≥1% found in the protease region of HIV-1 that were present in the 151 PR sequences examined. According to the IAS-USA, the mutations associated with drug resistance, with a p >10%, were L10I, M36I, I62V, L63P, I64V, A71V/T, V77I, L90M, and I93L[12, 36, 37].

Figure 1
figure 1

Genetic and phenotypic representation of primary structure variation within the HIV-1 protease consensus. Codons 1–3 are shown in grey and were not included in our analysis. Conserved regions are shown in yellow, semi-conserved regions in ochre, variable regions in orange and highly variable regions in red. The major (dark green) and minor (light green) resistance mutations are indicated for each codon. aGenetic variation = total number of mutations at the nucleotide position/number of sequences evaluated. bPhenotypic variation = total number of mutations at the amino acid position/number of sequences evaluated.

Figure 2
figure 2

Prevalence of mutations within HIV-1 PR pol . Red bars represent mutations associated with drug resistance, and the green bars represent natural polymorphisms and unusual mutations not associated with drug resistance.

Structural studies of PRs have reported a slight widening of the active site due to mutations associated with drug resistance for the majority of PIs [9, 10, 38]. However for other inhibitors, such as IDV which is characterized by three aromatic rings, structural changes are caused by mutations at the active site and adjacent positions [39].

Prevalence of natural polymorphisms and unusual mutations in PRs without established drug resistance

Table 1 shows the natural polymorphisms or unusual mutations with a p >1% that were found in the PR sequences of HIV-1 isolated from the Mexican patients. These are weakly associated with PI resistance, but are not included in the IAS–USA guides or the HIVdb as accessory or minor mutations [16, 40, 41]. Of the 14 mutations, only L63A and H69Y were found in drug resistance positions, and T12A/I, I15V, E35D, N37D/E, R57K, K70E and I72V were contiguous to positions associated with resistance. Overall, these mutations have little effect on drug susceptibility; however, a phenotypic change in any of them could have relevance to the affinity to one or more PIs [6, 42]. These mutations, in combination with resistance mutations, might have an effect on the dynamics of the evolution of cross-resistance [43].

Table 1 Polymorphisms or unusual mutations ( p> 1%) weakly associated with PI resistance in HIV-1 protease from treated and naïve individuals according to the HIVdb

The I13V (17.33%), E35D (18.21%), R41K (19%) and R57K (17.88%) mutations had a p ≥ 10% and were located in polymorphic positions observed in non-B subtypes [35, 55, 56]. In the HIVdb, W6R and I72R are unusual mutations with a frequency <0.05% that only emerge after multiple major and minor resistance mutations [57]. Table 2 shows 41 mutations with a p >1% that have not been associated with resistance, 25 are natural polymorphisms and the remaining 16 were unusual mutations. According to phenotypic conservation analysis, the L5F and Q7E mutations were within the conserved regions, while D29V, P39S, K43R, Q61E, E65D, C67F, P79L, T91V and Q92G/K were within semi-conserved regions. The T12P/S, K14R, G17D/E, Q18H, L19I/V/T, G68E, H69Q and K70R/T/I mutations were within the variable regions, and N37S/T/C/H/I, L63S/V/R/G/H, I72V/T/E/M and I93F were in highly variable regions.

Table 2 Natural polymorphisms and unusual mutations of HIV-1 protease (p> 1%) without evidence of resistance to PIs

Among the codons presented in the Table 2, the mutations in positions K43, L63, H69 and I93 were located in sites associated with minor resistance, but the distance between its localization and the enzyme’s active site reduces the possibility of the structure contributing to drug resistance. All the described mutations could be due to random transcriptional errors, or positive selection from drug and/or immunological stressors [37, 58]. Generally, natural polymorphisms occur in remote regions away from the active site, and form domains that define the shape of the homodimer. However, unusual mutations are found in positions associated with drug resistance and possibly generate allosteric changes in the binding site that favour enzymatic function, or decrease the affinity with certain PIs [59]. Therefore, the study of such structural changes produced by these emerging mutations may help in determining the new effects of PIs with different affinities.

Figure 3 shows PR tertiary structure positions that are: not associated with PI resistance; weakly associated with PI resistance; associated with PI resistance. We have also presented the locations of natural polymorphisms and unusual mutations (Figure 3). The codons T12, N37, L63, H69, K70 and I72 include mutations weakly associated with PI resistance (T12A/I, N37D/E, L63A, H69Y, K70E, and I72R), and mutations lacking evidence of PI resistance (T12P/S, N37S/T/C/H/I, L63S/V/R/G/H, H69Q, K70R/T/I, and I72V/T/E/M).

Figure 3
figure 3

Codons with natural polymorphisms and unusual mutations in the HIV-1 PR tertiary structure. Codons in the PR that were not associated with PI resistance (cyan), weakly associated with PI resistance (yellow), and associated with PI resistance (red).

The D29V and P79L mutations are located near the active site of the protease, and therefore possibly contribute to the generation of PI resistance. It is of interest to evaluate these unusual mutations in silico, and establish their association with resistance to PIs.

Phenotypic conservation of HIV-1 protease

Figure 4 shows the conserved, semi-conserved, variable and highly variable regions of PRs according to PV. Mutations were clustered into 15 regions, for amino acids 4–99 of the protease. For average PV calculation, when the asymmetry in the distribution was greater than 1.4 between the 15th and 75th percentiles, the residues were not considered. We found three conserved, three variable, three highly variable and six semi-conserved regions for each chain. The positions excluded from the PV calculated for each region were W6, L10, I13, K14, G17, Q18, E35, G40, R41, M46, I54, V56, R57, I63, V77, N83, L90, Q92, I93 and K97. The PV in these codons had very different values from those presented by the codons in their respective regions. According to our model of protease conservation, the LES formed by fragments 23–33 and 74–78 were in semi-conserved regions (E21–L34 and G73–P81, except for V77). The LES formed by the 83–92 fragment involved two codons with variable PV, I84 (6.29%) and L90 (12.33%), and two codons with semi-conserved PV, T91 (3.33%) and Q92 (4.05%) [8, 60]. Codon 90 contained a drug resistance mutation (L90M) common for most PIs, with the exception of DRV and TPV, while T91 and Q92 contained the T91V, Q92G, and Q92K mutations, which are classified in the literature as unusual mutations. The prevalence of the L90M, T91V, Q92G, and Q92K mutations was 12.0, 3.33, 2.03 and 2.03%, respectively. Although the effectiveness and specificity of PR proteolytic activity is determined by its active site (amino acids 25–29), these characteristics are influenced by mutations in neighbouring structures, which mainly affect intramolecular interactions with the active site [5, 38, 42, 61]. Contiguous regions and the active site have a semi-conserved state, with a PV of 1.2%. It has been shown that active sites with poor capacity to carry out structural changes help adjust the specificity of natural substrates without losing proteolytic effectiveness [45]. A study that identified the minimal conserved structure of HIV-1 PR, in the presence or absence of drug stress, showed that most of the PV is a product of pharmacological stress [62]. In contrast, the peripheral structural regions have a relatively high PV (for variable and highly variable regions) courtesy of negative selection, and to a lesser extent through resistance of HIV-1 to immune stress [63, 64].

Figure 4
figure 4

Phenotypic conservation of HIV-1 PR isolated from Mexican patients. A consensus sequence was obtained from 151 individual sequences. Regions are shown from red to yellow in proportion to their phenotypic variation (PV,%). Positions with natural polymorphisms and unusual mutations at drug resistance codons are shown in the protease model. Mutations were clustered into 15 regions between codons 4–99 of the protease.

Selective pressure in the pr fragment of HIV-1 pol

Antiretroviral treatment can exert strong selective pressures within pol, which transcribes PR, reverse transcriptase and integrase [62, 65]. We have presented the selection pressure results for 10 codons with natural polymorphisms and unusual mutations (Table 3). According to these results, codons 5, 29, 63, 79, 91 and 93 represent positive pressure (dN–dS > 1) through the ML substitution model using the HyPhy algorithm. When these results are compared with the data available in the UCLA HIV Positive Selection Mutation Database, only codons 63 and 93 were consistent with positive selection under immunological and/or pharmacological stressors [26]. The difference in selection pressure for codons 5, 29, 79 and 91 could be due to variability in the antiretroviral regimen sequences administered to Mexican patients. In addition to positive selection, the aligned sites often evolve at different rates because of other biological factors that include site-specific mutation rates and functional constraints of amino acid substitutions [66].

Table 3 Selection pressure for codons with unusual mutations and natural polymorphisms

The codons that were not associated with resistance due to pharmacological stress, and had PV ≥2% were D29 (2.65%) and P79 (2.48%). These were located near the active site of the enzyme; T91 (3.33%) was also found to be necessary for the establishment of the PR dimer. Codons associated with resistance due to pharmacological stress and PV ≥2% were I47 (2.65%), V82 (10.6%) and I84 (6.29%). Only one of these sequences belonged to a naïve individual, with a mutation at V82; the remaining sequences were from treated individuals (three were treated with reverse transcriptase inhibitors only, the remainder were given reverse transcriptase inhibitors and at least one PI).

Structural prediction of mutant HIV-1 PRs

The molecular structure of all mutant HIV-1 PRs was predicted by comparative homology modelling using the wild-type HIV-1 PR as a template [PDB: 1GNO]. This structure had higher sequence identity compared with the HIV-1 subtype B consensus PR sequence available from the HIVdb. Additional file 1: Table S1 shows the % identity, the expected value of the alignment with the template sequence (E), and the score for the absolute quality of the models. We modelled the proteins with unusual mutations (L5F, D29V, L63G, L63R, P79L and T91V), natural polymorphisms (L63H andL63S), and drug-resistant mutant PRs with single mutations or patterns of mutations (D30N, V32I, M36I, M46I, I47V, G48V, I50V, I50L, I54M, Q58E, T74P, L76V, V82A, V82L, N83D, N88S, I84V, and L90M).

The model’s accuracy was increased because of the identity between the mutant and template sequences; therefore, we concluded that the model was suitable for all structures. The low E values obtained from the modelled proteins indicate template identification, and adequate target template alignment [27]. The reliability of the predicted structures with natural polymorphisms and unusual mutations in drug resistance positions ranged 0.87–0.89, while positions for major mutation proteins ranged 0.83–0.91. The lower QMEAN4 values correlated with mutants containing patterns of resistance, as a result of the reduced identity of these proteins with respect to the template structure [67]. The QMEAN4 values were acceptable for all the modelled structures. Figure 5 shows the overlapping structures of the wild-type PR [PDB: 1GNO] and the D29V mutant, with high similarity between both structures, as well as a difference in the location of the mutation site (position 29).

Figure 5
figure 5

Wild-type and D29V mutant protease structures. The structure of the wild-type HIV-1 protease (WT) was obtained by X-ray crystallography. [PDB: 1GNO] (blue), and the mutant protease (green) can be clearly seen, with the red structures corresponding to the oxygen atoms of D29.

D29 plays a crucial role for the folding of retroviral PRs [38, 68]. Using crystallography, it has been shown that D29 forms hydrogen bonds with R87, which partially constitutes the highly conserved triad, G86 − R87 − D88 [62]. The loss of this specific interaction between α − helix 1 (residues 87–91) and D29 destabilizes the dimer interface [69]. The PR structures of related viruses such as HIV − 2, equine infectious anaemia virus (EIAV), feline immunodeficiency virus (FIV), rous sarcoma virus (RSV) and simian immunodeficiency virus (SIV), also demonstrate a proximity between side chains R87 and D29 [70].

The in silico modelling of mutant proteins generated structures very similar to those obtained by X-ray crystallography. The structures with natural polymorphisms, unusual mutations in drug resistance positions, and drug resistance mutations obtained by comparative homology modelling, were appropriate for molecular docking with their respective inhibitors.

Natural polymorphisms and unusual mutations in PRs and their effects on the free energy of binding by PIs

We have presented the free energy of binding (kcal/mol), as well as the average value of the five lowest energy conformations for the complexes formed by PRs and the main PIs (Table 4). The wild-type PR had the lowest free energy of binding for all PIs, except for IDV, compared with mutant PRs containing major and multiple drug resistance mutations. The magnitudes corresponding to the minor values of the free energy of binding to the reference protein were: wild-type protease-1GNO < major drug-resistant mutant proteases < multiple drug-resistant mutant proteases. The PIs had the greatest degree of affinity for PR 1GNO, consistent with the wild-type PR, whereas reduced affinity for mutant PRs was proportional to the number of mutations [46].

Table 4 Free energy of binding for protease-PI complexes obtained in silico

Among the PIs, IDV demonstrated a higher affinity for mutant proteins than PR [PDB: 1GNO]. Additionally, a study that correlated the in vivo genetic resistance of HIV-1 to IDV indicated that the development of resistance occurs through the combined effects of several mutations, which do not confer a measurable degree of resistance when they occur on their own [39]. For the other PIs, significant viral resistance has been shown to be a result of the appearance of one or two substitutions in drug-resistance positions [40, 71].

The difference between affinities of complexes formed by wild-type and drug-resistant PRs indicates some contribution of phenotypic changes towards PI resistance [72, 73]. The complexes with the largest differences involved ATV and DRV, both with a difference of −1.2 kcal/mol. This indicates high susceptibility of both compounds to drug resistance mutations. Lower differences were observed (between −0.7 and −0.3 kcal/mol) for other complexes, indicating these drug resistance mutations have a minor or supplementary effect [72, 73].

We obtained a positive value when we calculated the difference of the free energy of binding between the wild-type-IDV complex and the drug-resistant mutant-IDV complex. This is consistent with a high genetic barrier to resistance for IDV, which has lower susceptibility to drug-resistance mutations compared with other PIs [39]. When comparing the free energy of binding between the complexes with drug resistance mutations versus natural polymorphisms and unusual mutation complexes, resistance to ATV, LPV, NFV, and TPV was always observed. The PRs with L5F, D29V, L63G, L63H, L63R, L63S, and P79L mutations had lower or equal free energy of binding to ATV, LPV, NFV and TPV, than those with wild-type or drug-resistant PRs.

The complex formed by the D29V mutant showed considerable differences between the distance of the V29 and D30 Cα -atoms and the heteroatoms closest to the PIs (Table 4). This is probably because of the absence of the Cβ carboxyl group in the valine compared with the wild-type D29. The electrostatic interactions exercised by the D29 carboxyl oxygens provide stronger affinity to PIs in the active site, resulting in greater affinity compared with the V29 mutant [5, 74, 75]. The absence of V29 carboxyl oxygens decreases the level of interactions, thus decreasing affinity. Such differences can be observed when measuring the distance between the functional groups of the wild-type, resistant and D29V mutant PRs docked to DRV and TPV (Figure 6). For each natural polymorphism and unusual mutation, Table 5 shows the degree of resistance to PIs based on free energy of binding differences when compared with reference PRs.

Figure 6
figure 6

HIV-1 protease structures. Wild-type (WT), resistant and D29V mutant proteases coupled to darunavir (top) and tipranavir (bottom). The numbers in the left upper corner are the PDB ID numbers used to model darunavir and tipranavir with the protease and to measure the distance (Å) between functional groups (purple).

Table 5 Distances (Å) between the amino acid of the protease and the PI heteroatoms

Of the emerging mutations, D29V appears to favour resistance in silico in seven of nine PIs. Designing more effective DRV analogues requires an interaction between D29 and the bis-tetrahydrofuran ring, as this contributes to complex stability [5, 42]. All complexes that formed among the PRs with natural polymorphisms, unusual mutations and drug resistance mutations to TPV and LPV had similar free energies of binding. TPV mainly forms hydrogen bonds with residues D25, D29, D30, G48 and I50, while LPV interacts with G27, D29 and D30. A study that elucidated the mechanism by which PIs minimize the harmful effects of resistance mutations, showed that TPV, ATV, LPV, APV, IDV and DRV conserve their antiretroviral activity in the presence of drug resistance mutations. This phenomenon is due to the compensation of the loss of enthalpy (ΔH) with an entropy gain (−TΔS), except in the case of TPV [75]. Our results are consistent with another report that showed isolated strains with a high level of phenotypic resistance to LPV were susceptible to other PIs [76]. This corresponds with PR resistance to TPV and LPV that contain emerging mutations whose free energies of binding were greater than those obtained with wild-type PR.

We found a high prevalence (89%) of L63PGHRS mutations in HIV-1 variants isolated from Mexican individuals, probably because of the prevalence of HIV-1 subtype B [32, 33]. In the present study, among the functional groups found at position 63 (L63G, L63S, L63H and L63R), only glycine had hydrophobic characteristics, while serine was hydrophilic, and histidine and arginine were alkaline. These four mutations conferred resistance to NFV, ATV, TPV and LPV, most probably through an allosteric effect, given that the substitutions were not located close to the residues where the PIs bind [10, 74]. Few mutations at position 63 have been examined for their resistance effects to PIs. The L63P mutation has a compensatory effect that increases catalytic activity from 110% to 530%; when L63P is associated with M46I, it forms a combination that is resistant to APV, IDV, LPV or NFV [9, 77]. Residue 63 provides hydrophobic contacts between the slit of the loop formed by amino acids 38–42 and a β-sheet (residues 59–63) [74]. Although the study of mutations in this position has been limited to L63P to assess the effect of mutations that provide non-hydrophobic characteristics, alternative mechanisms could be shown by which HIV-1 PR compensates for pharmacological stressors.

Clinical characteristics of patients with unusual mutations at resistance sites and/or natural polymorphisms

Of the participating individuals, 48 of 151 (31.8%) showed resistance to at least one PI. Of these, 34 (70.8%) showed a high level of resistance, four (8.3%) showed intermediate levels of resistance, and 10 (20.8%) showed low level resistance.

Of the 151 sequences, 24 (15.9%) had one or more unusual mutations at resistance sites and/or natural polymorphisms. They were isolated from 21 (87.5%) male and three (12.5%) female patients; 70.8% of whom lived in the central-east of Mexico, and 29.2% in the north-west. Of these 24 patients, 23 (95.8%) received antiretroviral therapy, and one (4.2%) was naïve to treatment. The nucleoside reverse transcriptase inhibitors (NRTIs) and non-nucleoside reverse transcriptase inhibitors (NNRTIs) used, in order of frequency, were AZT, 3TC, ddI, d4T, ddC, NVP, EFV and ABC. The main PIs used were IDV, RTV and SQV. The average viral load in this group of patients was 228,225 virus copies/mL and a mean CD4+ lymphocyte count of 223 cells/μL. According to the case definition of HIV infection and AIDS by the Centers for Disease Control and Prevention (Atlanta, USA) [78], patients were classified as asymptomatic (n = 2, 8.5%), symptomatic (n = 6, 25%), AIDS (n = 12, 50%), and of unknown clinical category (n = 4, 16.5%).


The use of bioinformatics to identify potential mutations that confer resistance to antiretroviral drugs allows researchers to develop realistic three-dimensional models that illustrate the atomic interactions between an enzyme and its substrate. In silico, the structural correlation of natural polymorphisms and unusual mutations of drug resistance codons, allows the identification of HIV-1 variants resistant to PIs. The D29V mutation increases the probability of resistance to PIs as it generates unstable complexes at the HIV-1 protease active site. The prevalence of this mutation in different populations should be further studied, and parallel crystallographic studies are required to confirm our in silico findings.

Among mutant PRs-PIs complexes evaluated, TPV and LPV had free energies of binding greater than those obtained with wild-type PRs.

Furthermore, the presence of a high rate of L63P, I93L, V77I and I62V polymorphisms among the Mexican population is similar to that observed in patients that underwent antiretroviral treatments in other American and western European countries. These data reinforced the knowledge regarding the molecular epidemiology of the HIV-1 subtype B in Mexico through the presence of HIV polymorphisms.


The Contents of this publication are the authors responsibility and do not necessarily represent the official views of the Instituto Mexicano del Seguro Social.



USA Food and Drug Administration


Highly Active Antiretroviral Therapy


Human Immunodeficiency Virus


Local Elementary Structures


Molecular Evolution Genetics Analysis


Protein Date Bank


Protease Inhibitors


Phenotypic Variation


Qualitative Model Energy Analysis




Los Alamos HIV sequence Database


HIV Drug Resistance Database






Maximun Likelihood


synonimous site divergence


nonsynonimous site divergence.


  1. Shafer RW, Kantor R, Gonzales MJ: The genetic basis of HIV-1 resistance to reverse transcriptase and protease inhibitors. AIDS Rev. 2000, 2: 211-228.

    PubMed Central  PubMed  Google Scholar 

  2. Rambaut A, Posada D, Crandall KA, Holmes EC: The causes and consequences of HIV evolution. Nat Rev Genet. 2004, 5: 52-61. 10.1038/nrg1246.

    Article  PubMed  CAS  Google Scholar 

  3. Ribeiro RM, Bonhoeffer S: Production of resistant HIV mutants during antiretroviral therapy. Proc Natl Acad Sci USA. 2000, 97: 7681-7686. 10.1073/pnas.97.14.7681.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  4. Adamson CS, Freed EO: Recent progress in antiretrovirals - lessons from resistance. Drug Discov Today. 2008, 13: 424-432. 10.1016/j.drudis.2008.02.003.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  5. Ghosh AK, Sridhar PR, Leshchenko S, Hussain AK, Li J, Kovalevsky AY, Walters DE, Wedekind JE, Grum-Tokars V, Das D, Koh Y, Maeda K, Gatanaga H, Weber IT, Mitsuya H: Structure-based design of novel HIV-1 protease inhibitors to combat drug resistance. J Med Chem. 2006, 49: 5252-5261. 10.1021/jm060561m.

    Article  PubMed  CAS  Google Scholar 

  6. Ohtaka H, Schon A, Freire E: Multidrug resistance to HIV-1 protease inhibition requires cooperative coupling between distal mutations. Biochemistry. 2003, 42: 13659-13666. 10.1021/bi0350405.

    Article  PubMed  CAS  Google Scholar 

  7. Broglia RA, Provasi D, Vasile F, Ottolina G, Longhi R, Tiana G: A folding inhibitor of the HIV-1 protease. Proteins. 2006, 62: 928-933.

    Article  PubMed  CAS  Google Scholar 

  8. Broglia RA, Tiana G, Sutto L, Provasi D, Simona F: Design of HIV-1-PR inhibitors that do not create resistance: blocking the folding of single monomers. Protein Sci. 2005, 14: 2668-2681. 10.1110/ps.051670905.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  9. Piana S, Carloni P, Rothlisberger U: Drug resistance in HIV-1 protease: Flexibility-assisted mechanism of compensatory mutations. Protein Sci. 2002, 11: 2393-2402.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  10. Hornak V, Simmerling C: Targeting structural flexibility in HIV-1 protease inhibitor binding. Drug Discov Today. 2007, 12: 132-138. 10.1016/j.drudis.2006.12.011.

    Article  PubMed  CAS  Google Scholar 

  11. Shafer RW, Rhee SY, Pillay D, Miller V, Sandstrom P, Schapiro JM, Kuritzkes DR, Bennett D: HIV-1 protease and reverse transcriptase mutations for drug resistance surveillance. AIDS. 2007, 21: 215-223. 10.1097/QAD.0b013e328011e691.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  12. Johnson VA, Calvez V, Günthard HF, Paredes R, Pillay D, Shafer R, Wensing AM, Richman DD: Update of the drug resistance mutations in HIV-1: March 2013. Top Antivir Med. 2013, 21: 4-12.

    Google Scholar 

  13. Johnson JA, Li JF, Wei X, Lipscomb J, Irlbeck D, Craig C, Smith A, Bennett DE, Monsour M, Sandstrom P, Lanier ER, Heneine W: Minority HIV-1 drug resistance mutations are present in antiretroviral treatment-naive populations and associate with reduced treatment efficacy. PLoS Med. 2008, 5: e158-10.1371/journal.pmed.0050158.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Rapp BA, Wheeler DL: GenBank. Nucleic Acids Res. 2002, 30: 17-20. 10.1093/nar/30.1.17.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  15. Rose PP, Korber BT: Detecting hypermutations in viral sequences with an emphasis on G – > A hypermutation. Bioinformatics. 2000, 16: 400-401. 10.1093/bioinformatics/16.4.400.

    Article  PubMed  CAS  Google Scholar 

  16. Rhee SY, Gonzales MJ, Kantor R, Betts BJ, Ravela J, Shafer RW: Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. 2003, 31: 298-303. 10.1093/nar/gkg100.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  17. Rozanov M, Plikat U, Chappey C, Kochergin A, Tatusova T: A web-based genotyping resource for viral sequences. Nucleic Acids Res. 2004, 32: W654-W659. 10.1093/nar/gkh419.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  18. 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. 10.1093/molbev/msr121.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  19. Kuiken C, Korber B, Shafer RW: HIV sequence databases. AIDS Rev. 2003, 5: 52-61.

    PubMed Central  PubMed  Google Scholar 

  20. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.

    Article  PubMed  CAS  Google Scholar 

  21. Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.

    Article  PubMed  CAS  Google Scholar 

  22. Rhee SY, Kantor R, Katzenstein DA, Camacho R, Morris L, Sirivichayakul S, Jorgensen L, Brigido LF, Schapiro JM, Shafer RW, for the International Non Subtype: HIV-1 pol mutation frequency by subtype and treatment experience: extension of the HIVseq program to seven non-B subtypes. AIDS. 2006, 20: 643-651. 10.1097/01.aids.0000216363.36786.2b.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Pond SL, Frost SD, Muse SV: HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005, 21: 676-679. 10.1093/bioinformatics/bti079.

    Article  PubMed  CAS  Google Scholar 

  24. Yang Z, Goldman N, Friday A: Comparison of models for nucleotide substitution used in maximum-likelihood phylogenetic estimation. Mol Biol Evol. 1994, 11: 316-324.

    PubMed  CAS  Google Scholar 

  25. Muse SV, Gaut BS: A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome. Mol Biol Evol. 1994, 11: 715-724.

    PubMed  CAS  Google Scholar 

  26. Pan C, Kim J, Chen L, Wang Q, Lee C: The HIV positive selection mutation database. Nucleic Acids Res. 2007, 35: D371-D375. 10.1093/nar/gkl855.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  27. Bordoli L, Kiefer F, Arnold K, Benkert P, Battey J, Schwede T: Protein structure homology modeling using SWISS-MODEL workspace. Nat Protocols. 2008, 4: 1-13. 10.1038/nprot.2008.197.

    Article  Google Scholar 

  28. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The protein data bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  29. Benkert P, Tosatto SCE, Schomburg D: QMEAN: a comprehensive scoring function for model quality assessment. Proteins. 2008, 71: 261-277. 10.1002/prot.21715.

    Article  PubMed  CAS  Google Scholar 

  30. Trott O, Olson AJ: 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.

    PubMed Central  PubMed  CAS  Google Scholar 

  31. Delano WL: The PyMOL Molecular Graphics System. 2002, []

    Google Scholar 

  32. Vazquez-Valls E, Escoto-Delgadillo M, Lopez-Marquez FC, Castillero-Manzano M, Echegaray-Guerrero E, Bitzer-Quintero OK, Kobayashi-Gutierrez A, Torres-Mendoza BM: Molecular epidemiology of HIV type 1 in Mexico: emergence of BG and BF intersubtype recombinants. AIDS Res Hum Retroviruses. 2010, 26: 777-781. 10.1089/aid.2009.0195.

    Article  PubMed  CAS  Google Scholar 

  33. Avila-Ríos S, García-Morales C, Garrido-Rodríguez D, Ormsby CE, Hernández-Juan R, Andrade-Villanueva J, González-Hernández LA, Torres-Escobar I, Navarro-Álvarez S, Reyes-Terán G, Mexican HIV Molecular Epidemiology Project Group: National prevalence and trends of HIV transmitted drug resistance in Mexico. PLoS One. 2011, 6: e27812-10.1371/journal.pone.0027812.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Hartman JL, Garvik B, Hartwell L: Principles for the buffering of genetic variation. Science. 2001, 9 (291): 1001-4.

    Article  Google Scholar 

  35. Bürger R, Willensdorfer M, Nowak MA: Why are phenotypic mutation rates much higher than genotypic mutation rates?. Genetics. 2006, 172: 197-206.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Rhee SY, Liu TF, Holmes SP, Shafer RW: HIV-1 subtype b protease and reverse transcriptase amino acid covariation. PLoS Comput Biol. 2007, 3: e87-10.1371/journal.pcbi.0030087.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Chen L, Lee C: Distinguishing HIV-1 drug resistance, accessory, and viral fitness mutations using conditional selection pressure analysis of treated versus untreated patient samples. Biol Direct. 2006, 1: 14-10.1186/1745-6150-1-14.

    Article  PubMed Central  PubMed  Google Scholar 

  38. Bonomi M, Gervasio FL, Tiana G, Provasi D, Broglia RA, Parrinello M: Insight into the folding inhibition of the HIV-1 protease by a small peptide. Biophys J. 2007, 93: 2813-2821. 10.1529/biophysj.107.106369.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  39. Condra JH, Holder DJ, Schleif WA, Blahy OM, Danovich RM, Gabryelski LJ, Graham DJ, Laird D, Quintero JC, Rhodes A, Robbins HL, Roth E, Shivaprakash M, Yang T, Chodakewitz JA, Deutsch PJ, Leavitt RY, Massari FE, Mellors JW, Squires KE, Steigbigel RT, Teppler H, Emini EA: Genetic correlates of in vivo viral resistance to indinavir, a human immunodeficiency virus type 1 protease inhibitor. J Virol. 1996, 70: 8270-8276.

    PubMed Central  PubMed  CAS  Google Scholar 

  40. Zhang J, Hou T, Wang W, Liu JS: Detecting and understanding combinatorial mutation patterns responsible for HIV drug resistance. Proc Natl Acad Sci USA. 2010, 107: 1321-1326. 10.1073/pnas.0907304107.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  41. Shafer RW, Rhee SY, Bennett DE: Consensus drug resistance mutations for epidemiological surveillance: basic principles and potential controversies. Antivir Ther. 2008, 13 (Suppl 2): 59-68.

    PubMed Central  PubMed  Google Scholar 

  42. Ghosh AK, Chapsal BD, Weber IT, Mitsuya H: Design of HIV protease inhibitors targeting protein backbone: an effective strategy for combating drug resistance. Acc Chem Res. 2007, 41: 78-86.

    Article  PubMed  Google Scholar 

  43. Haq O, Levy R, Morozov A, Andrec M: Pairwise and higher-order correlations among drug-resistance mutations in HIV-1 subtype B protease. BMC Bioinformatics. 2009, 10: S10-

    Article  PubMed Central  PubMed  Google Scholar 

  44. Melnick L, Yang SS, Rossi R, Zepp C, Heefner D: An Escherichia coli expression assay and screen for human immunodeficiency virus protease variants with decreased susceptibility to indinavir. Antimicrob Agents Chemother. 1998, 42: 3256-65.

    PubMed Central  PubMed  CAS  Google Scholar 

  45. Svicher V, Ceccherini-Silberstein F, Erba F, Santoro M, Gori C, Bellocchi MC, Giannella S, Trotta MP, Monforte A, Antinori A, Perno CF: Novel human immunodeficiency virus type 1 protease mutations potentially involved in resistance to protease inhibitors. Antimicrob Agents Chemother. 2005, 49: 2015-2025. 10.1128/AAC.49.5.2015-2025.2005.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  46. Jacobsen H, Yasargil K, Winslow DL, Craig JC, Krohn A, Duncan IB, Mous J: Characterization of human immunodeficiency virus type 1 mutants with decreased sensitivity to proteinase inhibitor Ro 31–8959. Virology. 1995, 206: 527-534. 10.1016/S0042-6822(95)80069-7.

    Article  PubMed  CAS  Google Scholar 

  47. Markowitz M, Conant M, Hurley A, Schluger R, Duran M, Peterkin J, Chapman S, Patick A, Hendricks A, Yuen GJ, Hoskins W, Clendeninn N, Ho DD: A preliminary evaluation of nelfinavir mesylate, an inhibitor of human immunodeficiency virus (HIV)-1 protease, to treat HIV infection. J Infect Dis. 1998, 177: 1533-40. 10.1086/515312.

    Article  PubMed  CAS  Google Scholar 

  48. Marcelin AG, Flandre P, De Mendoza C, Roquebert B, Peytavin G, Valer L, Wirden M, Abbas S, Katlama C, Soriano V, Calvez V: Clinical validation of saquinavir/ritonavir genotypic resistance score in protease-inhibitor-experienced patients. Antivir Ther. 2007, 12: 247-52.

    PubMed  CAS  Google Scholar 

  49. Marcelin AG, Lamotte C, Delaugerre C, Ktorza N, Ait Mohand H, Cacace R, Bonmarchand M, Wirden M, Simon A, Bossi P, Bricaire F, Costagliola D, Katlama C, Peytavin G, Calvez V, Genophar Study Group: Genotypic inhibitory quotient as predictor of virological response to ritonavir-amprenavir in human immunodeficiency virus type 1 protease inhibitor-experienced patients. Antimicrob Agents Chemother. 2003, 47: 594-600. 10.1128/AAC.47.2.594-600.2003.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  50. Molla A, Brun S, Mo H, Real K, Podding J, Bernstein B, Hertogs K, Larder B, Lie Y, Hellmann N, Vasavanonda S, Chernyavskiy T, Freimuth W, Japour A, Sun E, Kempf D: Genotypic and phenotypic analysis of viral isolates from subjects with detectable viral load on therapy with ABT-378/ritonavir. Antiviral Ther. 2000, 5 (Suppl 3): 30-

    Google Scholar 

  51. Harrigan PR, Hertogs K, Verbiest W, Pauwels R, Larder B, Kemp S, Bloor S, Yip B, Hogg R, Alexander C, Montaner JS: Baseline HIV drug resistance profile predicts response to ritonavir-saquinavir protease inhibitor therapy in a community setting. AIDS. 1999, 13: 1863-1871. 10.1097/00002030-199910010-00008.

    Article  PubMed  CAS  Google Scholar 

  52. Shafer RW, Hsu P, Patick AK, Craig C, Brendel V: Identification of biased amino acid substitution patterns in human immunodeficiency virus type 1 isolates from patients treated with protease inhibitors. J Virol. 1999, 73: 6197-6202.

    PubMed Central  PubMed  CAS  Google Scholar 

  53. Carrillo A, Stewart KD, Sham HL, Norbeck DW, Kohlbrenner WE, Leonard JM, Kempf DJ, Molla A: In vitro selection and characterization of human immunodeficiency virus type 1 variants with increased resistance to ABT-378, a novel protease inhibitor. J Virol. 1998, 72: 7532-7541.

    PubMed Central  PubMed  CAS  Google Scholar 

  54. De Meyer S, Azijn H, Van Ginderen M, De Baere I, Pauwels R, De Béthune MP: In vitro selection experiments demonstrate an increased genetic barrier to resistance development to TMC114 as compared with currently licensed protease inhibitors. Antiviral Ther. 2002, 7: S5-

    Google Scholar 

  55. Coman RM, Robbins AH, Fernandez MA, Gilliland CT, Sochet AA, Goodenow MM, McKenna R, Dunn BM: The contribution of naturally occurring polymorphisms in altering the biochemical and structural characteristics of HIV-1 subtype C protease. Biochemistry. 2008, 47: 731-743. 10.1021/bi7018332.

    Article  PubMed  CAS  Google Scholar 

  56. Kantor R, Katzenstein DA, Efron B, Carvalho AP, Wynhoven B, Cane P, Clarke J, Sirivichayakul S, Soares MA, Snoeck J, Pillay C, Rudich H, Rodrigues R, Holguin A, Ariyoshi K, Bouzas MB, Cahn P, Sugiura W, Soriano V, Brigido LF, Grossman Z, Morris L, Vandamme AM, Tanuri A, Phanuphak P, Weber JN, Pillay D, Harrigan PR, Camacho R, Schapiro JM, et al: Impact of HIV-1 subtype and antiretroviral therapy on protease and reverse transcriptase genotype: results of a global collaboration. PLoS Med. 2005, 2: e112-10.1371/journal.pmed.0020112.

    Article  PubMed Central  PubMed  Google Scholar 

  57. The HIVdb Release Notes - HIV Drug Resistance Database. []

  58. Ali A, Bandaranayake RM, Cai Y, King NM, Kolli M, Mittal S, Murzycki JF, Nalam MN, Nalivaika EA, Ozen A, Prabu-Jeyabalan MM, Thayer K, Schiffer CA: Molecular basis for drug resistance in HIV-1 protease. Viruses. 2010, 2: 2509-2535. 10.3390/v2112509.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  59. Kozisek M, Bray J, Rezacova P, Saskova K, Brynda J, Pokorna J, Mammano F, Rulisek L, Konvalinka J: Molecular analysis of the HIV-1 resistance development: enzymatic activities, crystal structures, and thermodynamics of nelfinavir-resistant HIV protease mutants. J Mol Biol. 2007, 374: 1005-1016. 10.1016/j.jmb.2007.09.083.

    Article  PubMed  CAS  Google Scholar 

  60. Broglia RA: Learning to design resistance proof drugs from folding. Eur Phys J D At Mol Opt Plasma Phys. 2009, 51: 137-151.

    CAS  Google Scholar 

  61. Parera M, Fernández G, Clotet B, Mart¡nez MA: HIV-1 protease catalytic efficiency effects caused by random single amino acid substitutions. Mol Biol Evol. 2007, 24: 382-387.

    Article  PubMed  CAS  Google Scholar 

  62. Ceccherini-Silberstein F, Erba F, Gago F, Bertoli A, Forbici F, Bellocchi MC, Gori C, D'Arrigo R, Marcon L, Balotta C, Antinori A, Monforte A, Perno CF: Identification of the minimal conserved structure of HIV-1 protease in the presence and absence of drug pressure. AIDS. 2004, 18: F11-F19.

    Article  PubMed  CAS  Google Scholar 

  63. Boeri E, Gianotti N, Canducci F, Hasson H, Giudici B, Castagna A, Lazzarin A, Clementi M: Evolutionary characteristics of HIV type 1 variants resistant to protease inhibitors in the absence of drug-selective pressure. AIDS Res Hum Retroviruses. 2003, 19: 1151-1153. 10.1089/088922203771881257.

    Article  PubMed  CAS  Google Scholar 

  64. Karlsson AC, Deeks SG, Barbour JD, Heiken BD, Younger SR, Hoh R, Lane M, Sällberg M, Ortiz GM, Demarest JF, Liegler T, Grant RM, Martin JN, Nixon D: Dual pressure from antiretroviral therapy and cell-mediated immune response on the human immunodeficiency virus type 1 protease gene. J Virol. 2003, 77: 6743-6752. 10.1128/JVI.77.12.6743-6752.2003.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  65. Prabu-Jeyabalan M, Nalivaika EA, King NM, Schiffer CA: Viability of a drug-resistant human immunodeficiency virus type 1 protease variant: structural insights for better antiviral therapy. J Virol. 2003, 77: 1306-1315. 10.1128/JVI.77.2.1306-1315.2003.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  66. De Oliveira T, Salemi M, Gordon M, Vandamme AM, Van Rensburg EJ, Engelbrecht S, Coovadia HM, Cassol S: Mapping sites of positive selection and amino acid diversification in the HIV genome. Genetics. 2004, 167: 1047-1058. 10.1534/genetics.103.018135.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  67. Arnold K, Bordoli L, Kopp J+, Schwede T: The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006, 22: 195-201. 10.1093/bioinformatics/bti770.

    Article  PubMed  CAS  Google Scholar 

  68. Bonomi M, Barducci A, Gervasio FL, Parrinello M: Multiple routes and milestones in the folding of HIV-1 protease monomer. PLoS ONE. 2010, 5: e13208-10.1371/journal.pone.0013208.

    Article  PubMed Central  PubMed  Google Scholar 

  69. Levy Y, Caflisch A, Onuchic JN, Wolynes PG: The folding and dimerization of HIV-1 protease: evidence for a stable monomer from simulations. J Mol Biol. 2004, 340: 67-79. 10.1016/j.jmb.2004.04.028.

    Article  PubMed  CAS  Google Scholar 

  70. Sharp PM, Hahn BH: The evolution of HIV-1 and the origin of AIDS. Philos Trans R Soc Lond B Biol Sci. 2010, 365: 2487-2494. 10.1098/rstb.2010.0031.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  71. Wu TD, Schiffer CA, Gonzales MJ, Taylor J, Kantor R, Chou S, Israelski D, Zolopa AR, Fessel WJ, Shafer RW: Mutation patterns and structural correlates in human immunodeficiency virus type 1 protease following different protease inhibitor treatments. J Virol. 2003, 77: 4836-4847. 10.1128/JVI.77.8.4836-4847.2003.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  72. Bromberg Y, Rost B: Correlating protein function and stability through the analysis of single amino acid substitutions. BMC Bioinformatics. 2009, 10: S8-

    Article  PubMed Central  PubMed  Google Scholar 

  73. Gao C, Park MS, Stern HA: Accounting for ligand conformational restriction in calculations of protein-ligand binding affinities. Biophys J. 2010, 98: 901-910. 10.1016/j.bpj.2009.11.018.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  74. Heaslet H, Rosenfeld R, Giffin M, Lin YC, Tam K, Torbett BE, Elder JH, McRee DE, Stout CD: Conformational flexibility in the flap domains of ligand-free HIV protease. Acta Crystallogr D Biol Crystallogr. 2007, 63: 866-875. 10.1107/S0907444907029125.

    Article  PubMed  CAS  Google Scholar 

  75. Muzammil S, Armstrong AA, Kang LW, Jakalian A, Bonneau PR, Schmelmer V, Amzel LM, Freire E: Unique thermodynamic response of tipranavir to human immunodeficiency virus type 1 protease drug resistance mutations. J Virol. 2007, 81: 5144-5154. 10.1128/JVI.02706-06.

    Article  PubMed Central  PubMed  CAS  Google Scholar 

  76. Friend J, Parkin N, Liegler T, Martin JN, Deeks SG: Isolated lopinavir resistance after virological rebound of a ritonavir/lopinavir-based regimen. AIDS. 2004, 18: 1965-1966. 10.1097/00002030-200409240-00016.

    Article  PubMed  Google Scholar 

  77. Weikl TR, Hemmateenejad B: How conformational changes can affect catalysis, inhibition and drug resistance of enzymes with induced-fit binding mechanism such as the HIV-1 protease. Biochim Biophys Acta. 2013, 1834: 867-873. 10.1016/j.bbapap.2013.01.027.

    Article  PubMed  CAS  Google Scholar 

  78. Centers for Disease Control and Prevention: 1993 revised classification system for HIV infection and expanded surveillance case definition for AIDS among adolescents and adults. MMWR Recomm Rep. 1992, 41 (RR-17): 1-19.

    Google Scholar 

Download references


The study was funded in part by the Instituto Mexicano del Seguro Social, Grant Number FIS/IMSS/G09/752. The authors wish to thank Ms. Luz Segovia-Santos for her editorial services and to Greg Davies for his English revision services.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Eduardo Vázquez-Valls.

Additional information

Competing interests

The authors declare that they have no competing interest.

Authors’ contributions

CMM, MED, MVT and MFS performed sequences and bioinformatics analyses. EVV designed the study. AVO, FGG, MCM and LGFR provided clinical samples and collected the data. BTM performed the data analysis, CMM and EVV wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1: Identity value, expected value and QMEAN analyses of mutant proteases models tested to estimate the quality of the predicted structure. (DOCX 20 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mata-Munguía, C., Escoto-Delgadillo, M., Torres-Mendoza, B. et al. Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis. BMC Bioinformatics 15, 72 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: