Volume 16 Supplement 19
Structural studies on molecular mechanisms of Nelfinavir resistance caused by non-active site mutation V77I in HIV-1 protease
© Gupta et al. 2015
Published: 16 December 2015
The human immunodeficiency virus (HIV-1) is a retrovirus causing acquired immunodeficiency syndrome (AIDS), which has become a serious problem across the world and has no cure reported to date. Human immunodeficiency virus (HIV-1) protease is an attractive target for antiviral treatment and a number of therapeutically useful inhibitors have been designed against it. The emergence of drug resistant mutants of HIV-1 poses a serious problem for conventional therapies that have been used so far. Until now, thirteen protease inhibitors (PIs), major mutation sites and many secondary mutations have been listed in the HIV Drug Resistance Database. In this study, we have studied the effect of the V77I mutation in HIV-PR along with the co-occurring mutations L33F and K20T through multi-nanosecond molecular dynamics simulations. V77I is known to cause Nelfinavir (NFV) resistance in the subtype B population of HIV-1 protease. We have for the first time reported the effect of this clinically relevant mutation on the binding of Nelfinavir and the conformational flexibility of the protease.
Two HIV-PR mutants have been considered in this study - the Double Mutant Protease (DBM) V77I-L33F and Triple Mutant Protease (TPM) V77I-K20T-L33F. The molecular dynamics simulation studies were carried out and the RMSD trajectories of the unliganded wild type and mutated protease were found to be stable. The binding affinity of NFV with wild type HIV-PR was very high with a Glide XP docking score of -9.3 Kcal/mol. NFV showed decreased affinity towards DBM with a docking score of -8.0 Kcal/mol, whereas its affinity increased towards TPM (Glide XP score: -10.3). Prime/MM-GBSA binding free energy of the wild type, DBM and TPM HIV-PR docked structures were calculated as -38.9, -11.1 and -42.6 Kcal/mol respectively. The binding site cavity volumes of wild type, DBM and TPM protease were 1186.1, 1375.5 and 1042.5 Å3 respectively.
In this study, we have studied the structural roles of the two HIV-PR mutations by conducting molecular dynamics simulation studies of the wild type and mutant HIV-1 PRs. The present study proposes that DBM protease showed greater flexibility and the flap separation was greater with respect to the wild type protease. The cavity size of the MD-stabilized DBM was also found to be increased, which may be responsible for the decreased interaction of Nelfinavir with the cavity residues, thus explaining the decreased binding affinity. On the other hand, the binding affinity of NFV for TPM was found to be enhanced, accounted for by the decrease in cavity size of the mutant which facilitated strong interactions with the flap residues. The flap separation of TPM was less than the wild type protease and the decreased cavity size may be responsible for its lower resistance, and hence, may be the reason for its lower clinical relevance.
KeywordsHIV Nelfinavir mutation mechanism resistance docking molecular dynamics simulations
Human immunodeficiency virus type 1 (HIV-1), a member of the retrovirus family, causes acquired immunodeficiency syndrome (AIDS) which progressively destroys the body's natural immune system leading to deadly infections and cancers . According to the World Health Organization (WHO), HIV has affected 35 million people so far and continues to be a major public health burden globally, with 69% of the population affected alone in sub-Saharan Africa [2, 3]. The current therapeutic options available include highly active antiretroviral therapy (HAART) combination, a set of antiretroviral drugs which inhibit the replication of virus in the body and reduce the burden of the disease, but drugs or vaccines which can eradicate the viruses from the human body still remain a question to be addressed . Additionally, the existing antiretroviral drugs are very expensive and are associated with the risk of non-AIDS disorders that include cardiovascular, liver, kidney and neurological disease . Owing to the huge global burden of the disease, the lack of effective drugs and vaccines which can kill the pathogen, the high costs and the serious side effects related to the existing drugs and the accumulating drug resistance has made the search for anti-HIV drugs a foremost research priority.
Functional inhibition of HIV-1 PR leads to incomplete viral replication and therefore makes it an attractive target for anti-HIV drugs . To date, seven protease inhibitors (PIs) have been approved by the U.S. Food and Drug Administration . However, emergence of several resistant viral species due to genetic mutations in active and non-active sites of HIV-1 PR has made the current therapies inefficient and inappropriate for use. The active site mutations also referred to as primary mutations cause direct resistance, while the secondary mutations often accompany primary mutations (accessory mutations) or they show synergistic resistance in the presence of other secondary mutations . These mutations are a consequence of the selective pressure rising due to antiviral agents. Another driving force of the resistance causing mutations is the recently reported immunological pressure and the mutations are described as 'Escape mutations' .
V77I is one of the important non active site secondary mutations in HIV-PR, causing resistance against Nelfinavir (NFV). It is a highly polymorphic mutation present near the cheek sheet of protease, with a marked presence in subtype B virus. This minor mutation is accompanied by other primary and secondary mutations . L33F is a major mutation present in the active site of the protease , which results in reduced susceptibility towards NFV in the presence of other mutations. It is a non-polymorphic mutation which provides resistance against all PIs except Indinavir and Saquinavir. L33I is a less commonly occurring mutation with similar effects to L33F, while L33V mutation has not been related to any kind of drug resistance in PI therapy. L33F co-occurs with V77I in large number of HIV-1 subtype B infected patient samples, as reported in Stanford's HIV Drug Resistance Database . The mutation at the 20th residue is another non-polymorphic site present in the cheek turn and is involved in rendering resistance against all PIs except Saquinavir and Tipranavir. K20T is the most prominent mutation occurring at the 20th residue  and has been found to co-occur with V77I in subtype B population. It is interesting to note that the 77th, 33rd and 20th amino acids form a set of residues interacting with the 36th residue of the protease, which itself is present on non-active site and its mutation causes resistance to NFV in non-subtype B viruses . L23I, D30N, E35G, M46I/L/V, G48V, I54L, G73S/T/C/A, T74S, V82A/F/S/T, I84V, N88D/S and L90M are other mutations correlated to NFV resistance.
Number of subtype-b clinical isolates reported in HIV drug resistance database
Materials and methods
Preparation of wild type and mutated HIV-1 protease structures
The crystal structure of protease-NFV complex (PDB: 1OHR)  was downloaded from the Protein Data Bank  and was pre-processed using ViewerLite, a visualizing tool from Accelrys (Accelrys, Inc., San Diego, CA, USA). The ligand NFV and the water molecules were removed from the structure and the protein was further prepared and optimized using Schrödinger's protein preparation wizard [28, 29]. The residues V77 and L33 were mutated to isoleucine and phenyalanine to obtain double mutant DBM (Figure 1A) through the protein preparation wizard. Similarly triple mutant TPM (Figure 1C) was obtained by mutating residues V77, L33 and K20 to isoleucine, phenylalanine and threonine respectively. The preparation of structures involved addition of hydrogen bonds, creation of disulfide bonds, removal of bad contacts, capping of protein terminals, optimization of bond lengths, conversion of selanomethionine to methionine and cleaning the geometry of overlapping residues. Side chain prediction and refinement of selected residues was carried out using PRIME module of Schrodinger (Prime, version 2.1, Schrödinger, LLC, New York, NY, 2009). The study of flap movements is crucial to understand and compare the molecular dynamics of wild type and mutants HIV-1 PR. In the ligand-bound crystal structure, the flap residues were involved in interactions with the inhibitor or natural substrate [30, 31] while the flaps of unliganded protease were found to be highly flexible . To analyze the flap motions of our mutants with respect to the wild type protease, we selected the unliganded structure  of HIV-1 protease [PDB: 1HVP]. The double (DBM_M) and triple (TPM_M) mutants were prepared similarly to the procedure described above.
The drug Nelfinavir (NFV) (CID 64143), was processed before docking using LigPrep's ligand preparation protocol (Ligprep v2.5; Schrödinger, Inc.: Portland, OR, 2011). The three-dimensional coordinates (tautomeric, stereochemical, ionizing variants) were generated along with their energy minimization and flexible filtering.
Molecular dynamics simulations
MD simulations of the docked and unliganded complexes (both wild type and mutant) were accomplished using Desmond Molecular Dynamics system, with Optimized Potentials for Liquid Simulations (OPLS) all-atom force field 2005 [34, 35] as described in some of our previous studies [36–39]. The prepared protein molecules were solvated in the presence of explicit solvent on a fully hydrated model with TIP4P water model in a triclinic periodic boundary box (distance between box wall and protein complex was kept at 10 Å to avoid the direct interaction with its own periodic image) to generate required systems for MD simulations. The energy of prepared systems for MD simulations was minimized to 5000 steps maximum using the steepest descent method until a gradient threshold (25 kcal/mol/Å) was reached, followed by L-BFGS (Low-memory Broyden-Fletcher- Goldfarb Shanno quasi-Newtonian minimizer) until a convergence threshold of 1 kcal/mol/Å was met. The default parameters in Desmond were applied for systems equilibration. The so equilibrated systems were then used for simulations at 300 K temperature and a constant pressure of 1atm, with a time step of 2fs. The long range electrostatic interactions were handled using Smooth Particle Mesh Ewald Method. Cutoff method was selected to define the short range electrostatic interactions. A cutoff of 9 Å radius (default), was used.
The prepared conformations of NFV were docked to the stabilized mutants DBM and TPM using Glide docking software [40, 41]. We performed semi-flexible docking where the ligand was made flexible keeping the receptor macromolecule rigid. The flexibility of ligand molecule gives it the freedom to search from six degrees of rotational and translational freedom and an arbitrary number of torsional degrees of freedom. A random perturbation to each was applied at each time step, and the interaction energy was evaluated for the new location and conformation.
A scoring grid was prepared on the active site of the homodimer i.e. at the interface of both the subunits, using receptor grid generation platform of Schrödinger [40, 41]. Keeping all the parameters default, a grid of size 20 × 20 × 20 Å with an inner box size of 10 × 10 × 10 Å was generated.
Calculation of binding energies using MM/GBSA
Where ER + EL is the sum of energies of unbound ligand and receptor, and ER:L is the energy of the docked complex. ΔGSA is the difference of surface area energy of the protein-ligand complex and the sum of surface area energies of protein and ligand individually. ΔGSOLV is the difference in the GBSA solvation energy of the complex and summation of individual salvation energies of protein and ligand. Energies of the complex were calculated using the OPLS-All Atom force field  and GB/SA continuum solvent model.
Hydrogen bond and hydrophobic interaction analysis
The hydrophobic interactions and H-bonds of the docked complexes were analysed using the Ligplot program . The parameters defining the H-bonds between ligand and the protein complexes were as follows: acceptor-donor atoms distances less than 3.3 Å, hydrogen acceptor atom distances less than 2.7 Å and an acceptor-donor angle of 90° or more. Ligand-bound protease structures obtained from Glide and the MD-stabilized representative structures from Desmond were selected for carrying out interaction studies. A representative structure was prepared by averaging the coordinates of various frames extracted from the most stable region of the trajectory, which persisted until the end of the simulation run.
SYFPEITHI epitope prediction analysis
The human immune system proteins human leukocyte antigens (HLAs) bind to the intracellular epitopes arising from digestion of viral proteins. HLAs are responsible for presenting these epitopes on the cell surface and hence triggering an immune response against the virus. Escape mutations hinder strong binding of HLA to the epitopes and thereby assist in bypassing the immune response. Epitope prediction was done using Syfpeithi Database of MHC ligands and peptide motifs . The algorithm relies on the scoring of binding motifs. From the first amino acid of the protein, its sequence is divided into octamers, nonamers and decamers. The score of each oligomer is then calculated according to the summation of scores of individual amino acids. The amino acids are scored based on their observed frequencies. Most frequently occurring residues in the anchor positions are given a value 10, followed by a value of 8 given to the residues occurring in significant number of ligands. Likewise, residues regarded unfavorable for binding have a coefficient of -1 to -3. We have used MHC class HLA-A3 for our analysis.
Calculation of volume and surface area of HIV-1 PR cavity
We used CASTp server  to estimate the cavity volume and surface area of wild type and mutant proteases. CASTp works on the principles of Alpha Shape Theory  for detection and measurement of pockets in a protein which are inaccessible to the solvent outside. A probe of radius 1.4 Å was used for cavity measurement.
Results and discussion
Molecular dynamics simulations and energy stabilization of unliganded wild type, DBM and TPM protease
The representative structures from the MD simulation trajectories of DBM and TPM were then selected for studying their interactions with NFV. The mutants were docked with all the stable conformations of NFV using Glide and were compared with the NFV-docked crystal structure of wild type HIV-1 PR.
Docking of NFV with wild type protease
Arg 8(B), Leu 23(B), Gly 27(B), Asp 25 (A)(B), Ala 28 (A)(B), Asp 29(A)(B), Asp 30 (A)(B), Val 32(A), Ile 47(A), Gly 48(B), Ile 50 (B), Val 82(B), Ile 84(B)
Leu 23(B), Gly 27(B), Ala 28 (A), Asp 29(A), Val 32(B), Ile 47(B), Gly 48(B), Gly 49(B), Ile 50 (A)(B), Val 82(A)(B), Pro 81(A)(B), Ile 84(A)(B)
After MD Simulation
Arg 8(A)(B), Ala 28(A)(B), Asp 30(B), Val 32(B), Ile 47(A), Gly 48(B), Gly 29(B), Ile 50 (A)(B), Pro 81(A), Val 82(A), Ile 84(A)
Arg 8(A), Leu 23(A)(B), Gly 27(A), Ala 28 (A)(B), Asp 29(B), Asp 30(B), Gly 48(B), Ile 49(B), Ile 50(A), Ile 54(B), Thr 80(B), Phe 81(B), Ile 84(A)(B)
After MD Simulation
Gly 27(B), Asp 25(B), Ala 28 (B), Asp 29(B), Asp 30(B), Val 32(B), Ile 47(B), Gly 49(A)(B), Ile 50 (B), Thr 80(A), Val 81(A)(B), Ile 84(A)(B)
Leu 23(A), Gly 27(A)(B), Asp 25(B), Ala 28 (A)(B), Val 32(B), Gly 48(B), Val 82(B), Pro 81(A), Ile 84(A)
Docking of NFV with DBM and energy stabilization of NFV-DBM docked complex
Docking of NFV with TPM and energy stabilization of NFV-TPM docked complex
The results of docking and simulations studies of triple mutant TPM were notably different from other resistant proteases and therefore justified its lower clinical presence with respect to DBM. The docking scores of NFV with stable mutant TPM with NFV was found to be -10.31 Kcal/mol. The Prime/MM-GBSA free binding energy was stabilized by 3.68 Kcal/mol to -42 Kcal/mol with respect to the wild type protease. This binding involved five hydrogen bonds with residues Arg 8(A), Asp 25(A), Ile 50(A) and Gly 48(B) along with a number of hydrophobic interactions (Figure 4E, 5E). However, the number of interactions got significantly reduced after MD simulations. The NFV-TPM complex was found to be stable throughout the simulations trajectory with RMSD of only 0.183. The RMSF values of individual residues were also found lowered with respect to the undocked TPM. NFV was observed forming only two hydrogen bonds which were made with Asp 25(A). From this it becomes evident that the otherwise resistant K20T mutation has the potential to reduce the resistance of TPM in comparison to DBM (Figure 4F, 5F).
Possible intra-molecular interactions by V77I, L33F, K20T and the neighboring residues, and their effect on cavity size
Binding cavity size and area
VOLUME (Å 3 )
AREA (Å 2 )
VOLUME (Å 3 )
AREA (Å 2 )
VOLUME (Å 3 )
AREA (Å 2 )
V77I mutation in combination with L33F (DBM) displayed increase in the size of binding cavity (representative structure) with volume and area as 1375.5 Å3 and 732.10 Å2 respectively. This increase in cavity size is the probable reason behind decreased binding affinity of NFV to DBM (due to decrease in contact surface area of ligand and active site residues). L33F mutation caused positional change of neighboring residue Glu34, which resulted in formation of an extra hydrogen bond between Glu34 and Lys20 in DBM.
The binding pocket volume and area of the triple mutant representative structure got reduced with respect to the wild type protease, thereby increasing the contact surface area between ligand and active site residues. The pocket volume and surface area of TPM was found to be 1042.5 Å3 and 634.3 Å2 respectively. Thr20 formed an additional hydrogen bond with Gly21 in B chain, whereas hydrogen bonds between 33rd residue and Gly34 in B chain, and Leu76 in A chain were lost as a result of mutations.
Comparison of flap movements of double and triple mutant with the wild type protease
Putative selective mechanism of resistant mutations
Along with the selective pressure arising with anti-viral treatment, immunological pressure has also been reported as a sound theory behind the emergence of resistant mutations in HIV-1 PR. This is also credible in case of V77I mutation as it is located as an anchoring residue in the epitope recognized by HLA-A3 , (nonamer-LIGPTPVNI). The score representing probability of binding and presenting of the peptide by HLA-A3 was seen to be decreased to a small extent, suggesting that the emergence of mutation V77I is preferentially due to selective pressures imposed by anti-viral therapy, and less likely due to immunological pressure.
The present study explains the molecular mechanisms through which the V77I mutation in HIV-1 protease cause resistance towards NFV. Since this is a non-active site accessory mutation and clinically occurs with other resistant mutants, we have considered two types of mutant proteases double (DBM) V77I-L33F and triple mutant (TPM) V77I-L33F-K20T in the study. NFV showed a lower binding affinity towards DBM, and this mutant was found to be more stable than the wild type. The flap opening conformation of DBM_M suggested wider separation of flaps and higher flexibility, thus showing that the effect of mutation on the equilibrium of closed and semiopen conformations of protease which could be one reason behind the resistance showed by DBM. Further, the increased cavity size of DBM explained decreased binding affinity of NFV for the mutant protease being accounted by reduced contact surface area. TPM showed increased affinity towards drug and this explains the reason behind its less clinical prevalence. The decreased pocket size and stable flaps suggested that the combination of three mutations made HIV-1 PR non-resistant towards NFV and hence would not have been selected by nature. However the clinical presence of these three mutations together suggests that the mutant protease in nature may have been made resistant due to the presence of other mutations.
AG is thankful to Jawaharlal Nehru University for usage of all computational facilities. AG would like thank University Grants Commission, India for the Faculty Recharge position.
AG would like to acknowledge financial support from Department of Science and Technology, Government of India towards publication of this article.
This article has been published as part of BMC Bioinformatics Volume 16 Supplement 19, 2015: Brazilian Symposium on Bioinformatics 2014. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/16/S19
- Mitsuya H, Yarchoan R, Broder S: Molecular targets for AIDS therapy. Science. 1990, 249 (4976): 1533-1544.View ArticlePubMedGoogle Scholar
- Krausslich HG, Wimmer E: Viral proteinases. Annu Rev Biochem. 1988, 701-54. 57Google Scholar
- Kohl NE, Emini EA, Schleif WA, Davis LJ, Heimbach JC, Dixon RA, et al: Active human immunodeficiency virus protease is required for viral infectivity. Proc Natl Acad Sci USA. 1988, 85 (13): 4686-4690.View ArticlePubMedPubMed CentralGoogle Scholar
- Vondrasek J, Wlodawer A: HIVdb: a database of the structures of human immunodeficiency virus protease. Proteins. 2002, 49 (4): 429-431.View ArticlePubMedGoogle Scholar
- Navia MA, Fitzgerald PM, McKeever BM, Leu CT, Heimbach JC, Herber WK, et al: Three-dimensional structure of aspartyl protease from human immunodeficiency virus HIV-1. Nature. 1989, 337 (6208): 615-620.View ArticlePubMedGoogle Scholar
- Wlodawer A, Miller M, Jaskolski M, Sathyanarayana BK, Baldwin E, Weber IT, et al: Conserved folding in retroviral proteases: crystal structure of a synthetic HIV-1 protease. Science. 1989, 245 (4918): 616-621.View ArticlePubMedGoogle Scholar
- Louis JM, Dyda F, Nashed NT, Kimmel AR, Davies DR: Hydrophilic peptides derived from the transframe region of Gag-Pol inhibit the HIV-1 protease. Biochemistry. 1998, 37 (8): 2105-2110.View ArticlePubMedGoogle Scholar
- Prabu-Jeyabalan M, Nalivaika E, Schiffer CA: How does a symmetric dimer recognize an asymmetric substrate? a substrate complex of HIV-1 protease. J Mol Biol. 2000, 301 (5): 1207-1220.View ArticlePubMedGoogle Scholar
- Turner BG, Summers MF: Structural biology of HIV. J Mol Biol. 1999, 285 (1): 1-32.View ArticlePubMedGoogle Scholar
- Wlodawer A, Vondrasek J: Inhibitors of HIV-1 protease: a major success of structure-assisted drug design. Annu Rev Biophys Biomol Struct. 1998, 249-284. 27Google Scholar
- Arvieux C, Tribut O: Amprenavir or fosamprenavir plus ritonavir in HIV infection: pharmacology, efficacy and tolerability profile. Drugs. 2005, 65 (5): 633-659.View ArticlePubMedGoogle Scholar
- Leslie A, Kavanagh D, Honeyborne I, Pfafferott K, Edwards C, Pillay T, et al: Transmission and accumulation of CTL escape variants drive negative associations between HIV polymorphisms and HLA. J Exp Med. 2005, 201 (6): 891-902.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 (1): 298-303.View ArticlePubMedPubMed CentralGoogle Scholar
- Jallow S, Alabi A, Sarge-Njie R, Peterson K, Whittle H, Corrah T, et al: Virological response to highly active antiretroviral therapy in patients infected with human immunodeficiency virus type 2 (HIV-2) and in patients dually infected with HIV-1 and HIV-2 in the Gambia and emergence of drug-resistant variants. J Clin Microbiol. 2009, 47 (7): 2200-2208.View ArticlePubMedPubMed CentralGoogle Scholar
- Van Marck H, Dierynck I, Kraus G, Hallenberger S, Pattery T, Muyldermans G, et al: The impact of individual human immunodeficiency virus type 1 protease mutations on drug susceptibility is highly influenced by complex interactions with the background protease sequence. J Virol. 2009, 83 (18): 9512-9520.View ArticlePubMedPubMed CentralGoogle Scholar
- Ode H, Matsuyama S, Hata M, Neya S, Kakizawa J, Sugiura W, Hoshino T: Computational characterization of structural role of the non-active site mutation M36I of human immunodeficiency virus type 1 protease. J Mol Biol. 2007, 370 (3): 598-607.View ArticlePubMedGoogle Scholar
- Soni S, Tyagi C, Grover A, Goswami SK: Molecular modeling and molecular dynamics simulations based structural analysis of the SG2NA protein variants. BMC Res Notes. 2014, 7 (1): 446-View ArticlePubMedPubMed CentralGoogle Scholar
- Goyal M, Grover S, Dhanjal JK, Goyal S, Tyagi C, Chacko S, Grover A: Novel Natural Structure Corrector of ApoE4 for Checking Alzheimer's Disease: Benefits from High Throughput Screening and Molecular Dynamics Simulations. BioMed Res Int. 2013, Article ID 620793 doi:10.1155/2013/620793Google Scholar
- Ode H, Ota M, Neya S, Hata M, Sugiura W, Hoshino T: Resistant mechanism against Nelfinavir of human immunodeficiency virus type 1 proteases. J Phys Chem B. 2005, 109 (1): 565-574.View ArticlePubMedGoogle Scholar
- Skalova T, Dohnalek J, Duskova J, Petrokova H, Hradilek M, Soucek M, et al: HIV-1 protease mutations and inhibitor modifications monitored on a series of complexes. Structural basis for the effect of the A71V mutation on the active site. J Med Chem. 2006, 49 (19): 5777-5784.View ArticlePubMedGoogle Scholar
- Piana S, Carloni P, Rothlisberger U: Drug resistance in HIV-1 protease: Flexibility-assisted mechanism of compensatory mutations. Protein Science. 2002, 11 (10): 2393-2402.View ArticlePubMedPubMed CentralGoogle Scholar
- Ode H, Neya S, Hata M, Sugiura W, Hoshino T: Computational simulations of HIV-1 proteases--multi-drug resistance due to nonactive site mutation L90M. J Am Chem Soc. 2006, 128 (24): 7887-7895.View ArticlePubMedGoogle Scholar
- Meiselbach H, Horn AH, Harrer T, Sticht H: Insights into amprenavir resistance in E35D HIV-1 protease mutation from molecular dynamics and binding free-energy calculations. J Mol Model. 2007, 13 (2): 297-304.View ArticlePubMedGoogle Scholar
- Batista PR, Wilter A, Durham EH, Pascutti PG: Molecular dynamics simulations applied to the study of subtypes of HIV-1 protease common to Brazil, Africa, and Asia. Cell Biochem Biophys. 2006, 44 (3): 395-404.View ArticlePubMedGoogle Scholar
- Ode H, Matsuyama S, Hata M, Hoshino T, Kakizawa J, Sugiura W: Mechanism of drug resistance due to N88S in CRF01_AE HIV-1 protease, analyzed by molecular dynamics simulations. J Med Chem. 2007, 50 (8): 1768-1777.View ArticlePubMedGoogle Scholar
- Kaldor SW, Kalish VJ, Davies JF, Shetty BV, Fritz JE, Appelt K, et al: Viracept (Nelfinavir mesylate, AG1343): a potent, orally bioavailable inhibitor of HIV-1 protease. J Med Chem. 1997, 40 (24): 3979-3985.View ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al: The Protein Data Bank. Nucleic Acids Res. 2000, 28 (1): 235-242.View ArticlePubMedPubMed CentralGoogle Scholar
- Schrodinger L: Schrodinger Software Suite. 2011, New York: Schrödinger, LLCGoogle Scholar
- Sastry GM, Adzhigirey M, Day T, Annabhimoju R, Sherman W: Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments. J Comput Aided Mol Des. 2013, 27 (3): 221-234.View ArticlePubMedGoogle Scholar
- Miller M, Schneider J, Sathyanarayana BK, Toth MV, Marshall GR, Clawson L, et al: Structure of complex of synthetic HIV-1 protease with a substrate-based inhibitor at 2.3 A resolution. Science. 1989, 246 (4934): 1149-1152.View ArticlePubMedGoogle Scholar
- Hornak V, Okur A, Rizzo RC, Simmerling C: HIV-1 protease flaps spontaneously open and reclose in molecular dynamics simulations. Proc Natl Acad Sci USA. 2006, 103 (4): 915-920.View ArticlePubMedPubMed CentralGoogle Scholar
- Ishima R, Freedberg DI, Wang YX, Louis JM, Torchia DA: Flap opening and dimer-interface flexibility in the free and inhibitor-bound HIV protease, and their implications for function. Structure. 1999, 7 (9): 1047-1055.View ArticlePubMedGoogle Scholar
- Weber IT, Miller M, Jaskolski M, Leis J, Skalka AM, Wlodawer A: Molecular modeling of the HIV-1 protease and its substrate binding site. Science. 1989, 243 (4893): 928-931.View ArticlePubMedGoogle Scholar
- Kaminski GA, Friesner RA, Tirado-Rives J, Jorgensen WL: Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides†. J Phys Chem B. 2001, 105 (28): 6474-6487.View ArticleGoogle Scholar
- Jorgensen WL, Maxwell DS, Tirado-Rives J: Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J Am Chem Soc. 1996, 118 (45): 11225-11236.View ArticleGoogle Scholar
- Goyal S, Grover S, Dhanjal JK, Goyal M, Tyagi C, Chacko S, Grover A: Mechanistic insights into mode of actions of novel oligopeptidase B inhibitors for combating leishmaniasis. J Mol Model. 2014, 20 (3): 1-9.View ArticleGoogle Scholar
- Goyal M, Dhanjal JK, Goyal S, Tyagi C, Hamid R, Grover A: Development of dual inhibitors against Alzheimer's disease using fragment-based QSAR and molecular docking. Biomed Research International. 2014, Article ID 979606Google Scholar
- Tyagi C, Grover S, Dhanjal JK, Goyal S, Goyal M, Grover A: Mechanistic insights into mode of action of novel natural cathepsin L inhibitors. BMC Genomics. 2013, 14 (Suppl 8): S10.-View ArticlePubMedPubMed CentralGoogle Scholar
- Dhanjal JK, Goyal S, Sharma S, Hamid R, Grover A: Mechanistic insights into mode of action of potent natural antagonists of BACE-1 for checking Alzheimer's plaque pathology. Biochem Biophys Res Commun. 2014, 443 (3): 1054-1059.View ArticlePubMedGoogle Scholar
- Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, et al: Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem. 2004, 47 (7): 1739-1749.View ArticlePubMedGoogle Scholar
- Halgren TA, Murphy RB, Friesner RA, Beard HS, Frye LL, Pollard WT, Banks JL: Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J Med Chem. 2004, 47 (7): 1750-1759.View ArticlePubMedGoogle Scholar
- Lyne PD, Lamb ML, Saeh JC: Accurate prediction of the relative potencies of members of a series of kinase inhibitors using molecular docking and MM-GBSA scoring. J Med Chem. 2006, 49 (16): 4805-4808.View ArticlePubMedGoogle Scholar
- Wallace AC, Laskowski RA, Thornton JM: LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Engineering. 1995, 8 (2): 127-134.View ArticlePubMedGoogle Scholar
- Rammensee H, Bachmann J, Emmerich NP, Bachor OA, Stevanovic S: SYFPEITHI: database for MHC ligands and peptide motifs. Immunogenetics. 1999, 50 (3-4): 213-219.View ArticlePubMedGoogle Scholar
- Dundas J, Ouyang Z, Tseng J, Binkowski A, Turpaz Y, Liang J: CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues. Nucleic Acids Res. 2006, 34 (Web Server): W116-W118.View ArticlePubMedPubMed CentralGoogle Scholar
- Liang J, Edelsbrunner H, Woodward C: Anatomy of protein pockets and cavities: measurement of binding site geometry and implications for ligand design. Protein Science. 1998, 7 (9): 1884-1897.View ArticlePubMedPubMed CentralGoogle Scholar
- Toth G, Borics A: Flap opening mechanism of HIV-1 protease. J Mol Graphics Modell. 2006, 24 (6): 465-474.View ArticleGoogle Scholar
- John M, Moore CB, James IR, Mallal SA: Interactive selective pressures of HLA-restricted immune responses and antiretroviral drugs on HIV-1. Antiviral Therapy. 2005, 10 (4): 551-555.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.