1
|
Wang C, Chen L, Wang R, Tang W, Zhao B. Effects of the G48M mutant on the dynamics properties and binding mechanism of PR with SQV and ATV. MOLECULAR SIMULATION 2022. [DOI: 10.1080/08927022.2022.2055013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Chao Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Wanxia Tang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
- Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar, People’s Republic of China
| |
Collapse
|
2
|
Bastys T, Gapsys V, Walter H, Heger E, Doncheva NT, Kaiser R, de Groot BL, Kalinina OV. Non-active site mutants of HIV-1 protease influence resistance and sensitisation towards protease inhibitors. Retrovirology 2020; 17:13. [PMID: 32430025 PMCID: PMC7236880 DOI: 10.1186/s12977-020-00520-6] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2020] [Accepted: 05/04/2020] [Indexed: 02/07/2023] Open
Abstract
Background HIV-1 can develop resistance to antiretroviral drugs, mainly through mutations within the target regions of the drugs. In HIV-1 protease, a majority of resistance-associated mutations that develop in response to therapy with protease inhibitors are found in the protease’s active site that serves also as a binding pocket for the protease inhibitors, thus directly impacting the protease-inhibitor interactions. Some resistance-associated mutations, however, are found in more distant regions, and the exact mechanisms how these mutations affect protease-inhibitor interactions are unclear. Furthermore, some of these mutations, e.g. N88S and L76V, do not only induce resistance to the currently administered drugs, but contrarily induce sensitivity towards other drugs. In this study, mutations N88S and L76V, along with three other resistance-associated mutations, M46I, I50L, and I84V, are analysed by means of molecular dynamics simulations to investigate their role in complexes of the protease with different inhibitors and in different background sequence contexts. Results Using these simulations for alchemical calculations to estimate the effects of mutations M46I, I50L, I84V, N88S, and L76V on binding free energies shows they are in general in line with the mutations’ effect on \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$IC_{50}$$\end{document}IC50 values. For the primary mutation L76V, however, the presence of a background mutation M46I in our analysis influences whether the unfavourable effect of L76V on inhibitor binding is sufficient to outweigh the accompanying reduction in catalytic activity of the protease. Finally, we show that L76V and N88S changes the hydrogen bond stability of these residues with residues D30/K45 and D30/T31/T74, respectively. Conclusions We demonstrate that estimating the effect of both binding pocket and distant mutations on inhibitor binding free energy using alchemical calculations can reproduce their effect on the experimentally measured \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$IC_{50}$$\end{document}IC50 values. We show that distant site mutations L76V and N88S affect the hydrogen bond network in the protease’s active site, which offers an explanation for the indirect effect of these mutations on inhibitor binding. This work thus provides valuable insights on interplay between primary and background mutations and mechanisms how they affect inhibitor binding.
Collapse
Affiliation(s)
- Tomas Bastys
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics, 66123, Saarbrücken, Germany.,Saarbrücken Graduate School of Computer Science, University of Saarland, 66123, Saarbrücken, Germany
| | - Vytautas Gapsys
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, 37077, Göttingen, Germany
| | - Hauke Walter
- Medizinisches Labor Stendal, 39576, Stendal, Germany
| | - Eva Heger
- Institute of Virology, University of Cologne, 50935, Cologne, Germany
| | - Nadezhda T Doncheva
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics, 66123, Saarbrücken, Germany.,Faculty of Health and Medical Sciences, University of Copenhagen, 2200, Copenhagen, Denmark
| | - Rolf Kaiser
- Institute of Virology, University of Cologne, 50935, Cologne, Germany
| | - Bert L de Groot
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, 37077, Göttingen, Germany
| | - Olga V Kalinina
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics, 66123, Saarbrücken, Germany. .,Helmholtz Institute for Pharmaceutical Research Saarland (HIPS), Helmholtz Centre for Infection Research (HZI), 66123, Saarbrücken, Germany. .,Faculty of Medicine, Saarland University, 66421, Homburg, Germany.
| |
Collapse
|
3
|
Sk MF, Roy R, Kar P. Exploring the potency of currently used drugs against HIV-1 protease of subtype D variant by using multiscale simulations. J Biomol Struct Dyn 2020; 39:988-1003. [PMID: 32000612 DOI: 10.1080/07391102.2020.1724196] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Acquired immune deficiency syndrome (AIDS) is caused by the human immunodeficiency virus (HIV), type 1 and 2. Further, the diversity in HIV-1 has given rise to many serotypes and recombinant strains. The currently used protease inhibitors have been developed for subtype B, although non-B subtype strains account for ∼ 90% of the global HIV infections. Subtype D is spreading rapidly and infecting a large population in North Africa and the Middle East. In the current study, molecular dynamics simulations in conjunction with the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) scheme was used to investigate the potency of four drugs, namely atazanavir (ATV), darunavir (DRV), lopinavir (LPV) and tipranavir (TPV) against the subtype D variant. Our calculations predicted that the potency of the inhibitors decreased in the order TPV > ATV > DRV > LPV. TPV was found to be the most potent against subtype D due to an increase in van der Waals and electrostatic interactions and reduction in the desolvation energy compared to other inhibitors. This result is further supported by the hydrogen bond interactions between inhibitors and protease. Furthermore, our analyses suggested that the binding of TPV induced a more closed conformation of the flap compared to apo or other complexes. It was observed that TPV/PRD has a lower cavity volume relative to the other three complexes leading to a tighter binding. The open conformation of the flap was observed for LPV/PRD. We expect that this study might be useful for designing more potent inhibitors against HIV-1 subtype D. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Md Fulbabu Sk
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Simrol Campus, Indore, Madhya Pradesh, India
| | - Rajarshi Roy
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Simrol Campus, Indore, Madhya Pradesh, India
| | - Parimal Kar
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Simrol Campus, Indore, Madhya Pradesh, India
| |
Collapse
|
4
|
Henes M, Kosovrasti K, Lockbaum GJ, Leidner F, Nachum GS, Nalivaika EA, Bolon DN, Yilmaz NK, Schiffer CA, Whitfield TW. Molecular Determinants of Epistasis in HIV-1 Protease: Elucidating the Interdependence of L89V and L90M Mutations in Resistance. Biochemistry 2019; 58:3711-3726. [PMID: 31386353 PMCID: PMC6941756 DOI: 10.1021/acs.biochem.9b00446] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Protease inhibitors have the highest potency among antiviral therapies against HIV-1 infections, yet the virus can evolve resistance. Darunavir (DRV), currently the most potent Food and Drug Administration-approved protease inhibitor, retains potency against single-site mutations. However, complex combinations of mutations can confer resistance to DRV. While the interdependence between mutations within HIV-1 protease is key for inhibitor potency, the molecular mechanisms that underlie this control remain largely unknown. In this study, we investigated the interdependence between the L89V and L90M mutations and their effects on DRV binding. These two mutations have been reported to be positively correlated with one another in HIV-1 patient-derived protease isolates, with the presence of one mutation making the probability of the occurrence of the second mutation more likely. The focus of our investigation is a patient-derived isolate, with 24 mutations that we call "KY"; this variant includes the L89V and L90M mutations. Three additional KY variants with back-mutations, KY(V89L), KY(M90L), and the KY(V89L/M90L) double mutation, were used to experimentally assess the individual and combined effects of these mutations on DRV inhibition and substrate processing. The enzymatic assays revealed that the KY(V89L) variant, with methionine at residue 90, is highly resistant, but its catalytic function is compromised. When a leucine to valine mutation at residue 89 is present simultaneously with the L90M mutation, a rescue of catalytic efficiency is observed. Molecular dynamics simulations of these DRV-bound protease variants reveal how the L90M mutation induces structural changes throughout the enzyme that undermine the binding interactions.
Collapse
Affiliation(s)
- Mina Henes
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Klajdi Kosovrasti
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Gordon J. Lockbaum
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Florian Leidner
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Gily S. Nachum
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Ellen A. Nalivaika
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Daniel N.A. Bolon
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Nese Kurt Yilmaz
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
| | - Celia A. Schiffer
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA,Corresponding Author Celia A. Schiffer: Phone: +1 508 856 8008; , Troy W. Whitfield: Phone: +1 508 856 4401;
| | - Troy W. Whitfield
- Department of Medicine, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA,Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA,Corresponding Author Celia A. Schiffer: Phone: +1 508 856 8008; , Troy W. Whitfield: Phone: +1 508 856 4401;
| |
Collapse
|
5
|
C.S. V, Tamizhselvi R, Munusami P. Exploring the drug resistance mechanism of active site, non-active site mutations and their cooperative effects in CRF01_AE HIV-1 protease: molecular dynamics simulations and free energy calculations. J Biomol Struct Dyn 2019; 37:2608-2626. [DOI: 10.1080/07391102.2018.1492459] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Vasavi C.S.
- School of Biosciences and Technology, VIT University, Vellore, India
| | | | - Punnagai Munusami
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad, India
| |
Collapse
|
6
|
Bastys T, Gapsys V, Doncheva NT, Kaiser R, de Groot BL, Kalinina OV. Consistent Prediction of Mutation Effect on Drug Binding in HIV-1 Protease Using Alchemical Calculations. J Chem Theory Comput 2018; 14:3397-3408. [PMID: 29847122 DOI: 10.1021/acs.jctc.7b01109] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Despite a large number of antiretroviral drugs targeting HIV-1 protease for inhibition, mutations in this protein during the course of patient treatment can render them inefficient. This emerging resistance inspired numerous computational studies of the HIV-1 protease aimed at predicting the effect of mutations on drug binding in terms of free binding energy Δ G, as well as in mechanistic terms. In this study, we analyze ten different protease-inhibitor complexes carrying major resistance-associated mutations (RAMs) G48V, I50V, and L90M using molecular dynamics simulations. We demonstrate that alchemical free energy calculations can consistently predict the effect of mutations on drug binding. By explicitly probing different protonation states of the catalytic aspartic dyad, we reveal the importance of the correct choice of protonation state for the accuracy of the result. We also provide insight into how different mutations affect drug binding in their specific ways, with the unifying theme of how all of them affect the crucial drug binding regions of the protease.
Collapse
Affiliation(s)
- Tomas Bastys
- Department for Computational Biology and Applied Algorithmics , Max Planck Institute for Informatics , D-66123 Saarbrücken , Germany.,Saarbrücken Graduate School of Computer Science , University of Saarland , D-66123 Saarbrücken , Germany
| | - Vytautas Gapsys
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics , Max Planck Institute for Biophysical Chemistry , D-37077 Göttingen , Germany
| | - Nadezhda T Doncheva
- Department for Computational Biology and Applied Algorithmics , Max Planck Institute for Informatics , D-66123 Saarbrücken , Germany.,Faculty of Health and Medical Sciences , University of Copenhagen , 2200 Copenhagen , Denmark
| | - Rolf Kaiser
- Institute for Virology , University Clinic of Cologne , D-50935 Köln , Germany
| | - Bert L de Groot
- Computational Biomolecular Dynamics Group, Department of Theoretical and Computational Biophysics , Max Planck Institute for Biophysical Chemistry , D-37077 Göttingen , Germany
| | - Olga V Kalinina
- Department for Computational Biology and Applied Algorithmics , Max Planck Institute for Informatics , D-66123 Saarbrücken , Germany
| |
Collapse
|
7
|
Vasavi C, Tamizhselvi R, Munusami P. Drug Resistance Mechanism of L10F, L10F/N88S and L90M mutations in CRF01_AE HIV-1 protease: Molecular dynamics simulations and binding free energy calculations. J Mol Graph Model 2017. [DOI: 10.1016/j.jmgm.2017.06.007] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
8
|
Nishiyama K. Exploration of peptides that fit into the thermally vibrating active site of cathepsin K protease by alternating artificial intelligence and molecular simulation. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.06.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
9
|
Flynn WF, Haldane A, Torbett BE, Levy RM. Inference of Epistatic Effects Leading to Entrenchment and Drug Resistance in HIV-1 Protease. Mol Biol Evol 2017; 34:1291-1306. [PMID: 28369521 PMCID: PMC5435099 DOI: 10.1093/molbev/msx095] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Understanding the complex mutation patterns that give rise to drug resistant viral strains provides a foundation for developing more effective treatment strategies for HIV/AIDS. Multiple sequence alignments of drug-experienced HIV-1 protease sequences contain networks of many pair correlations which can be used to build a (Potts) Hamiltonian model of these mutation patterns. Using this Hamiltonian model, we translate HIV-1 protease sequence covariation data into quantitative predictions for the probability of observing specific mutation patterns which are in agreement with the observed sequence statistics. We find that the statistical energies of the Potts model are correlated with the fitness of individual proteins containing therapy-associated mutations as estimated by in vitro measurements of protein stability and viral infectivity. We show that the penalty for acquiring primary resistance mutations depends on the epistatic interactions with the sequence background. Primary mutations which lead to drug resistance can become highly advantageous (or entrenched) by the complex mutation patterns which arise in response to drug therapy despite being destabilizing in the wildtype background. Anticipating epistatic effects is important for the design of future protease inhibitor therapies.
Collapse
Affiliation(s)
- William F. Flynn
- Department of Physics and Astronomy, Rutgers University, New Brunswick, NJ
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA
| | - Allan Haldane
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA
- Department of Chemistry, Temple University, Philadelphia, PA
| | - Bruce E. Torbett
- Department of Molecular and Experimental Medicine, The Scripps Research Institute, La Jolla, CA
| | - Ronald M. Levy
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, PA
- Department of Chemistry, Temple University, Philadelphia, PA
| |
Collapse
|
10
|
Duan LL, Zhu T, Li YC, Zhang QG, Zhang JZH. Effect of polarization on HIV-1protease and fluoro-substituted inhibitors binding energies by large scale molecular dynamics simulations. Sci Rep 2017; 7:42223. [PMID: 28155907 PMCID: PMC5290483 DOI: 10.1038/srep42223] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Accepted: 01/06/2017] [Indexed: 01/01/2023] Open
Abstract
Molecular dynamics simulations in explicit water are carried out to study the binding of six inhibitors to HIV-1 protease (PR) for up to 700 ns using the standard AMBER force field and polarized protein-specific charge (PPC). PPC is derived from quantum mechanical calculation for protein in solution and therefore it includes electronic polarization effect. Our results show that in all six systems, the bridging water W301 drifts away from the binding pocket in AMBER simulation. However, it is very stable in all six complexes systems using PPC. Especially, intra-protease, protease-inhibitor hydrogen bonds are dynamic stabilized in MD simulation. The computed binding free energies of six complexes have a significantly linear correlation with those experiment values and the correlation coefficient is found to be 0.91 in PPC simulation. However, the result from AMBER simulation shows a weaker correlation with the correlation coefficient of −0.51 due to the lack of polarization effect. Detailed binding interactions of W301, inhibitors with PR are further analyzed and discussed. The present study provides important information to quantitative understanding the interaction mechanism of PR-inhibitor and PR-W301 and these data also emphasizes the importance of both the electronic polarization and the bridging water molecule in predicting precisely binding affinities.
Collapse
Affiliation(s)
- Li L Duan
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - T Zhu
- Department of Chemistry, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Yu C Li
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - Qing G Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250014, China
| | - John Z H Zhang
- Department of Chemistry, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
11
|
Decision tree for the binding of dipeptides to the thermally fluctuating surface of cathepsin K. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2016.01.035] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
12
|
Gupta A, Jamal S, Goyal S, Jain R, Wahi D, Grover A. Structural studies on molecular mechanisms of Nelfinavir resistance caused by non-active site mutation V77I in HIV-1 protease. BMC Bioinformatics 2015; 16 Suppl 19:S10. [PMID: 26695135 PMCID: PMC4686784 DOI: 10.1186/1471-2105-16-s19-s10] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Background 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. Results 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. Conclusion 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.
Collapse
|
13
|
Takaoka T, Mori K, Okimoto N, Neya S, Hoshino T. Prediction of the Structure of Complexes Comprised of Proteins and Glycosaminoglycans Using Docking Simulation and Cluster Analysis. J Chem Theory Comput 2015; 3:2347-56. [PMID: 26636224 DOI: 10.1021/ct700029q] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A typical docking simulation provides information on the structure of ligand-receptor complexes and their binding affinity in terms of a docking energy. We have developed a potent method combining a docking simulation with cluster analysis to extract adequate docking structures from the many possible output structures of the simulation. First, we tried to predict the structure of basic fibroblast growth factor (bFGF) bound to heparin, using the docking simulation program AutoDock 3.0. Two X-ray crystal structures had already been obtained for bFGF. One was a complex of the protein and heparin, a kind of glycosaminoglycan, and the other, only the protein itself, hereafter called a simplex. We docked a heparin molecule onto the protein simplex and generated many trial structures for the bFGF-heparin complex. The structures of those docked complexes were optimized through energy minimization by AMBER8. Although neither the docking energy calculated by AMBER8 nor that calculated by AutoDock 3.0 could be used satisfactorily by themselves to select a proper heparin-binding complex from the output structures, the majority of the structures generated by AutoDock 3.0 were fairly close to each other in atom geometry, and the averaged geometry over these structures was also close to that of the crystal. Hence, we utilized only the atom geometry for evaluation and carried out cluster analysis with the collection of geometries. This procedure enabled selection of a structure considerably close to the crystal's. We applied this approach to two other heparin-binding proteins: antithrombin and annexin V. Two crystal structures, a complex and a simplex, had been elucidated for these proteins as well as for bFGF. Our trials gave an exact prediction of the heparin-binding structures of these proteins, showing the approach in this study is effective in studying the docking of ligands that have a variety of docking conformations due to the presence of multiple rotatable bonds and charged chemical groups.
Collapse
Affiliation(s)
- Tsubasa Takaoka
- Graduate School of Pharmaceutical Sciences, Chiba University, Chiba 263-8522, Japan, Bioinformatics Group, GSC, RIKEN, Yokohama, Kanagawa 230-0046, Japan, and PRESTO, Japan Science and Technology Agency, Kawaguchi, Saitama 332-0012, Japan
| | - Kenichi Mori
- Graduate School of Pharmaceutical Sciences, Chiba University, Chiba 263-8522, Japan, Bioinformatics Group, GSC, RIKEN, Yokohama, Kanagawa 230-0046, Japan, and PRESTO, Japan Science and Technology Agency, Kawaguchi, Saitama 332-0012, Japan
| | - Noriaki Okimoto
- Graduate School of Pharmaceutical Sciences, Chiba University, Chiba 263-8522, Japan, Bioinformatics Group, GSC, RIKEN, Yokohama, Kanagawa 230-0046, Japan, and PRESTO, Japan Science and Technology Agency, Kawaguchi, Saitama 332-0012, Japan
| | - Saburo Neya
- Graduate School of Pharmaceutical Sciences, Chiba University, Chiba 263-8522, Japan, Bioinformatics Group, GSC, RIKEN, Yokohama, Kanagawa 230-0046, Japan, and PRESTO, Japan Science and Technology Agency, Kawaguchi, Saitama 332-0012, Japan
| | - Tyuji Hoshino
- Graduate School of Pharmaceutical Sciences, Chiba University, Chiba 263-8522, Japan, Bioinformatics Group, GSC, RIKEN, Yokohama, Kanagawa 230-0046, Japan, and PRESTO, Japan Science and Technology Agency, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
14
|
Mahdi M, Szojka Z, Mótyán JA, Tőzsér J. Inhibition Profiling of Retroviral Protease Inhibitors Using an HIV-2 Modular System. Viruses 2015; 7:6152-62. [PMID: 26633459 PMCID: PMC4690855 DOI: 10.3390/v7122931] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2015] [Revised: 11/11/2015] [Accepted: 11/13/2015] [Indexed: 12/13/2022] Open
Abstract
Retroviral protease inhibitors (PIs) are fundamental pillars in the treatment of HIV infection and acquired immunodeficiency syndrome (AIDS). Currently used PIs are designed against HIV-1, and their effect on HIV-2 is understudied. Using a modular HIV-2 protease cassette system, inhibition profiling assays were carried out for protease inhibitors both in enzymatic and cell culture assays. Moreover, the treatment-associated resistance mutations (I54M, L90M) were introduced into the modular system, and comparative inhibition assays were performed to determine their effect on the susceptibility of the protease. Our results indicate that darunavir, saquinavir, indinavir and lopinavir were very effective HIV-2 protease inhibitors, while tipranavir, nelfinavir and amprenavir showed a decreased efficacy. I54M, L90M double mutation resulted in a significant reduction in the susceptibility to most of the inhibitors with the exception of tipranavir. To our knowledge, this modular system constitutes a novel approach in the field of HIV-2 protease characterization and susceptibility testing.
Collapse
Affiliation(s)
- Mohamed Mahdi
- Laboratory of Retroviral Biochemistry, Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Debrecen, H-4010 Debrecen, Hungary.
| | - Zsófia Szojka
- Laboratory of Retroviral Biochemistry, Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Debrecen, H-4010 Debrecen, Hungary.
| | - János András Mótyán
- Laboratory of Retroviral Biochemistry, Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Debrecen, H-4010 Debrecen, Hungary.
| | - József Tőzsér
- Laboratory of Retroviral Biochemistry, Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Debrecen, H-4010 Debrecen, Hungary.
| |
Collapse
|
15
|
Qi F, Fudo S, Neya S, Hoshino T. A Dominant Factor for Structural Classification of Protein Crystals. J Chem Inf Model 2015; 55:1673-85. [DOI: 10.1021/acs.jcim.5b00052] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fei Qi
- Graduate School of Pharmaceutical
Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba 260-8675, Japan
| | - Satoshi Fudo
- Graduate School of Pharmaceutical
Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba 260-8675, Japan
| | - Saburo Neya
- Graduate School of Pharmaceutical
Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba 260-8675, Japan
| | - Tyuji Hoshino
- Graduate School of Pharmaceutical
Sciences, Chiba University, Inohana 1-8-1, Chuo-ku, Chiba 260-8675, Japan
| |
Collapse
|
16
|
Dependence of glycine peptide behavior on thermal fluctuations on the surface of cathepsin K. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.05.003] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
17
|
Meher BR, Wang Y. Exploring the drug resistance of V32I and M46L mutant HIV-1 protease to inhibitor TMC114: flap dynamics and binding mechanism. J Mol Graph Model 2014; 56:60-73. [PMID: 25562662 DOI: 10.1016/j.jmgm.2014.11.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2014] [Revised: 10/14/2014] [Accepted: 11/07/2014] [Indexed: 12/22/2022]
Abstract
Inhibitors of HIV-1 protease (HIV-1-pr) generally only bind to the active site of the protease. However, for some mutants such as V32I and M46L the TMC114 can bind not only to the active cavity but also to the groove of the flexible flaps. Although the second binding site suggests the higher efficiency of the drug against HIV-1-pr, the drug resistance in HIV-1-pr due to mutations cannot be ignored, which prompts us to investigate the molecular mechanisms of drug resistance and behavior of double bound TMC114 (2T) to HIV-1-pr. The conformational dynamics of HIV-1-pr and the binding of TMC114 to the WT, V32I and M46L mutants were investigated with all-atom molecular dynamic (MD) simulation. The 20 ns MD simulation shows many fascinating effects of the inhibitor binding to the WT and mutant proteases. MM-PBSA calculations explain the binding free energies unfavorable for the M46L and V32I mutants as compared to the WT. For the single binding (1T) the less binding affinity can be attributed to the entropic loss for both V32I-1T and M46L-1T. Although the second binding of TMC114 with flap does increase binding energy for the mutants (V32I-2T and M46L-2T), the considerable entropy loss results in the lower binding Gibbs free energies. Thus, binding of TMC114 in the flap region does not help much in the total gain in binding affinity of the system, which was verified from this study and thereby validating experiments.
Collapse
Affiliation(s)
- Biswa Ranjan Meher
- Computational Chemistry Laboratory, Department of Natural Sciences, Albany State University, Albany, GA 31705, USA
| | - Yixuan Wang
- Computational Chemistry Laboratory, Department of Natural Sciences, Albany State University, Albany, GA 31705, USA.
| |
Collapse
|
18
|
Ragland DA, Nalivaika EA, Nalam MNL, Prachanronarong KL, Cao H, Bandaranayake RM, Cai Y, Kurt-Yilmaz N, Schiffer CA. Drug resistance conferred by mutations outside the active site through alterations in the dynamic and structural ensemble of HIV-1 protease. J Am Chem Soc 2014; 136:11956-63. [PMID: 25091085 PMCID: PMC4151706 DOI: 10.1021/ja504096m] [Citation(s) in RCA: 67] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
![]()
HIV-1 protease inhibitors are part
of the highly active antiretroviral
therapy effectively used in the treatment of HIV infection and AIDS.
Darunavir (DRV) is the most potent of these inhibitors, soliciting
drug resistance only when a complex combination of mutations occur
both inside and outside the protease active site. With few exceptions,
the role of mutations outside the active site in conferring resistance
remains largely elusive. Through a series of DRV–protease complex
crystal structures, inhibition assays, and molecular dynamics simulations,
we find that single and double site mutations outside the active site
often associated with DRV resistance alter the structure and dynamic
ensemble of HIV-1 protease active site. These alterations correlate
with the observed inhibitor binding affinities for the mutants, and
suggest a network hypothesis on how the effect of
distal mutations are propagated to pivotal residues at the active
site and may contribute to conferring drug resistance.
Collapse
Affiliation(s)
- Debra A Ragland
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School , Worcester, Massachusetts 01605, United States
| | | | | | | | | | | | | | | | | |
Collapse
|
19
|
Tzoupis H, Leonis G, Avramopoulos A, Mavromoustakos T, Papadopoulos MG. Systematic molecular dynamics, MM-PBSA, and ab initio approaches to the saquinavir resistance mechanism in HIV-1 PR due to 11 double and multiple mutations. J Phys Chem B 2014; 118:9538-52. [PMID: 25036111 DOI: 10.1021/jp502687q] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Mutations in the human immunodeficiency virus (HIV) enable virus replication even when appropriate antiretroviral therapy is followed, thus leading to the emergence of drug resistance. In a previous work, we systematically examined seven single mutations that are associated with saquinavir (SQV) resistance in HIV-1 protease (Tzoupis, H.; Leonis, G.; Mavromoustakos, T.; Papadopoulos, M. G. J. Chem. Theory Comput. 2013, 9, 1754-1764). Herein, we extend our analysis, which includes seven double (G48V-V82A, L10I-G48V, G48V-L90M, I84V-L90M, L10I-V82A, L10I-L63P, A71V-G73S) and four multiple (L10I-L63P-A71V, L10I-G48V-V82A, G73S-I84V-L90M, L10I-L63P-A71V-G73S-I84V-L90M) SQV-HIV-1 PR mutant complexes, in an attempt to generalize our findings and formulate the main elements of the SQV resistance mechanism in the protease. On the basis of molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and ab initio computational approaches, we identified specific features that constitute the HIV-1 PR mechanism of resistance at the molecular level: the low flexibility of SQV in the binding cavity and the preservation of hydrogen bonding (HB) and van der Waals interactions between SQV and several active-site (Gly27/27', Asp29/29'/30/30', especially Asp25/25') and flap (Ile50/50', Gly48/48') residues of the protease contribute significantly to efficient binding. The total enthalpy loss in all mutants is mostly due to the loss in enthalpy of the active-site region. Furthermore, it was observed that mutation accumulation may induce stabilization to SQV and to the flaps through enhanced HB interactions that lead to improved inhibition (e.g., accumulation of mutations in complexes containing L10I, G48V, L63P, I84V, or L90M single mutations). It was also concluded that permanent flap closure is obtained independently of mutations and SQV binding is mostly driven by van der Waals, nonpolar, and exchange-energy contributions. Importantly, it was indicated that the optimal positioning of SQV and the structure of the binding cavity are tightly coupled, since small changes in geometry may affect the binding energy greatly. The results of our theoretical approaches are in agreement with experimental evidence and provide a reliable description of SQV resistance in HIV-1 PR.
Collapse
Affiliation(s)
- Haralambos Tzoupis
- Institute of Biology, Pharmaceutical Chemistry and Biotechnology, National Hellenic Research Foundation , 48 Vas. Constantinou Avenue, Athens 11635, Greece
| | | | | | | | | |
Collapse
|
20
|
A study of the interaction between HIV-1 protease and C 2-symmetric inhibitors by computational methods. J Mol Model 2014; 20:2369. [DOI: 10.1007/s00894-014-2369-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2014] [Accepted: 06/29/2014] [Indexed: 12/14/2022]
|
21
|
Leonis G, Steinbrecher T, Papadopoulos MG. A Contribution to the Drug Resistance Mechanism of Darunavir, Amprenavir, Indinavir, and Saquinavir Complexes with HIV-1 Protease Due to Flap Mutation I50V: A Systematic MM–PBSA and Thermodynamic Integration Study. J Chem Inf Model 2013; 53:2141-53. [DOI: 10.1021/ci4002102] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Georgios Leonis
- Institute of Biology, Medicinal
Chemistry and Biotechnology, National Hellenic Research Foundation, 48 Vas. Constantinou Avenue, Athens 11635,
Greece
| | - Thomas Steinbrecher
- Institute of Physical
Chemistry, Department of Theoretical Chemical Biology, KIT, Kaiserstrasse 12, 76131 Karlsruhe,
Germany
| | - Manthos G. Papadopoulos
- Institute of Biology, Medicinal
Chemistry and Biotechnology, National Hellenic Research Foundation, 48 Vas. Constantinou Avenue, Athens 11635,
Greece
| |
Collapse
|
22
|
Meher BR, Kumar MVS, Bandyopadhyay P. Interchain hydrophobic clustering promotes rigidity in HIV-1 protease flap dynamics: new insights from molecular dynamics. J Biomol Struct Dyn 2013; 32:899-915. [PMID: 23782135 DOI: 10.1080/07391102.2013.795873] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Abstract
The dynamics of HIV-1 protease (HIV-pr), a drug target for HIV infection, has been studied extensively by both computational and experimental methods. The flap dynamics of HIV-pr is considered to be more important for better ligand binding and enzymatic actions. Moreover, it has been demonstrated that the drug-induced mutations can change the flap dynamics of HIV-pr affecting the binding affinity of the ligands. Therefore, detailed understanding of flap dynamics is essential for designing better inhibitors. Previous computational investigations observed significant variation in the flap opening in nanosecond time scale indicating that the dynamics is highly sensitive to the simulation protocols. To understand the sensitivity of the flap dynamics on the force field and simulation protocol, molecular dynamics simulations of HIV-pr have been performed with two different AMBER force fields, ff99 and ff02. Two different trajectories (20 ns each) were obtained using the ff99 and ff02 force field. The results showed polarizable force field (ff02) make the flap tighter than the nonpolarizable force field (ff99). Some polar interactions and hydrogen bonds involving flap residues were found to be stronger with ff02 force field. The formation of interchain hydrophobic cluster (between flap tip of one chain and active site wall of another chain) was found to be dominant in the semi-open structures obtained from the simulations irrespective of the force field. It is proposed that an inhibitor, which will promote this interchain hydrophobic clustering, may make the flaps more rigid, and presumably the effect of mutation would be small on ligand binding.
Collapse
Affiliation(s)
- Biswa Ranjan Meher
- a Computational Biology Research Laboratory, Department of Biotechnology , Indian Institute of Technology , Guwahati , Assam , 781 039 , India
| | | | | |
Collapse
|
23
|
Tzoupis H, Leonis G, Mavromoustakos T, Papadopoulos MG. A Comparative Molecular Dynamics, MM–PBSA and Thermodynamic Integration Study of Saquinavir Complexes with Wild-Type HIV-1 PR and L10I, G48V, L63P, A71V, G73S, V82A and I84V Single Mutants. J Chem Theory Comput 2013; 9:1754-64. [DOI: 10.1021/ct301063k] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Haralambos Tzoupis
- Institute of Biology, Medicinal
Chemistry and Biotechnology, National Hellenic Research Foundation,
48 Vas. Constantinou Ave., Athens 11635, Greece
- Department of Chemistry, National
and Kapodistrian University of Athens, Panepistimioupolis Zographou
15771, Greece
| | - Georgios Leonis
- Institute of Biology, Medicinal
Chemistry and Biotechnology, National Hellenic Research Foundation,
48 Vas. Constantinou Ave., Athens 11635, Greece
| | - Thomas Mavromoustakos
- Department of Chemistry, National
and Kapodistrian University of Athens, Panepistimioupolis Zographou
15771, Greece
| | - Manthos G. Papadopoulos
- Institute of Biology, Medicinal
Chemistry and Biotechnology, National Hellenic Research Foundation,
48 Vas. Constantinou Ave., Athens 11635, Greece
| |
Collapse
|
24
|
|
25
|
How conformational changes can affect catalysis, inhibition and drug resistance of enzymes with induced-fit binding mechanism such as the HIV-1 protease. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:867-73. [PMID: 23376188 DOI: 10.1016/j.bbapap.2013.01.027] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2012] [Revised: 01/21/2013] [Accepted: 01/24/2013] [Indexed: 11/21/2022]
Abstract
A central question is how the conformational changes of proteins affect their function and the inhibition of this function by drug molecules. Many enzymes change from an open to a closed conformation upon binding of substrate or inhibitor molecules. These conformational changes have been suggested to follow an induced-fit mechanism in which the molecules first bind in the open conformation in those cases where binding in the closed conformation appears to be sterically obstructed such as for the HIV-1 protease. In this article, we present a general model for the catalysis and inhibition of enzymes with induced-fit binding mechanism. We derive general expressions that specify how the overall catalytic rate of the enzymes depends on the rates for binding, for the conformational changes, and for the chemical reaction. Based on these expressions, we analyze the effect of mutations that mainly shift the conformational equilibrium on catalysis and inhibition. If the overall catalytic rate is limited by product unbinding, we find that mutations that destabilize the closed conformation relative to the open conformation increase the catalytic rate in the presence of inhibitors by a factor exp(ΔΔGC/RT) where ΔΔGC is the mutation-induced shift of the free-energy difference between the conformations. This increase in the catalytic rate due to changes in the conformational equilibrium is independent of the inhibitor molecule and, thus, may help to understand how non-active-site mutations can contribute to the multi-drug-resistance that has been observed for the HIV-1 protease. A comparison to experimental data for the non-active-site mutation L90M of the HIV-1 protease indicates that the mutation slightly destabilizes the closed conformation of the enzyme. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.
Collapse
|
26
|
Nishiyama K. Local fluctuation control of papain by changing a highly fluctuating residue. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2012.10.058] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
27
|
Meher BR, Patel S. Structural and dynamical aspects of HIV-1 protease and its role in drug resistance. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2013; 92:299-324. [PMID: 23954105 DOI: 10.1016/b978-0-12-411636-8.00008-0] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Acquired immunodeficiency syndrome (AIDS) caused by the retrovirus human immunodeficiency virus (HIV) has become a major epidemic afflicting mankind. The Joint United Nations Program on HIV/AIDS (UNAIDS) projection shows the existence of millions of AIDS patients at the end of 2012. All the Food and Drug Administration (FDA)-approved drugs are getting ineffective due to resistance offered by the mutation-prone HIV. Hence, there is an urgent need for developing new drugs with greater potential. HIV life cycle is controlled by the activities of its essential proteins like glycoproteins (gp41 and gp120), HIV reverse transcriptase (HIV-RT), HIV integrase (HIV-IN), and HIV-1 protease (HIV-pr). This chapter focuses on the protein HIV-pr, which is important for the cleavage of Gag and Gag-Pol polyproteins to form mature, structural, and functional virions. The conformation and dynamics of the protein HIV-pr play a pivotal role in ligand binding and the catalytic process, which is affected by the rapid point mutations and various physiological parameters. The effect of the mutations and the varied simulation protocols on conformational dynamics and drug resistance of HIV-pr is discussed.
Collapse
Affiliation(s)
- Biswa Ranjan Meher
- Department of Chemistry, Gottwald Center for the Sciences, University of Richmond, Richmond, Virginia, USA.
| | | |
Collapse
|
28
|
Fujisaki S, Takashita E, Yokoyama M, Taniwaki T, Xu H, Kishida N, Sato H, Tashiro M, Imai M, Odagiri T. A single E105K mutation far from the active site of influenza B virus neuraminidase contributes to reduced susceptibility to multiple neuraminidase-inhibitor drugs. Biochem Biophys Res Commun 2012; 429:51-6. [PMID: 23131559 DOI: 10.1016/j.bbrc.2012.10.095] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2012] [Accepted: 10/24/2012] [Indexed: 11/28/2022]
Abstract
Drugs inhibiting the enzymatic activity of influenza virus neuraminidase (NA) are the cornerstone of therapy for influenza virus infection. The emergence of drug-resistant variants may limit the benefits of antiviral therapy. Here we report the recovery of an influenza B virus with reduced susceptibilities to NA inhibitors from a human patient with no history of antiviral drug treatment. The virus, designated B/Kochi/61/2011, was isolated by inoculating Madin-Darby canine kidney (MDCK) cells with respiratory specimens from the patient. NA inhibition assays demonstrated that the B/Kochi/61/2011 isolate showed a remarkable reduction in susceptibility to peramivir. The isolate also exhibited low to moderately reduced sensitivity to oseltamivir, laninamivir, and zanamivir. A sequence analysis of viruses propagated in MDCK cells revealed that the isolate contained a mutation (E105K) not previously associated with reduced susceptibility to NA inhibitors. However, pyrosequencing analysis showed that the NA E105K mutation was below a detectable level in the original clinical specimens, suggesting that the mutant virus may be preferably selected during propagation in MDCK cells. Analysis of the three-dimensional model of E105 and K105 NAs with peramivir suggested that the E105K mutation at the monomer-monomer interface of the NA tetramer may destabilize the tetrameric form of NA, leading to decreased susceptibility to NA inhibitors. These results have implications for understanding the mechanism of resistance against NA-inhibitor drugs.
Collapse
Affiliation(s)
- Seiichiro Fujisaki
- Laboratory of Influenza Virus Surveillance, Influenza Virus Research Center, National Institute of Infectious Diseases, Tokyo 208-0011, Japan
| | | | | | | | | | | | | | | | | | | |
Collapse
|
29
|
MEHER BISWARANJAN, SATISH KUMAR MATTAPARTHIVENKATA, SHARMA SMRITI, BANDYOPADHYAY PRADIPTA. CONFORMATIONAL DYNAMICS OF HIV-1 PROTEASE: A COMPARATIVE MOLECULAR DYNAMICS SIMULATION STUDY WITH MULTIPLE AMBER FORCE FIELDS. J Bioinform Comput Biol 2012; 10:1250018. [DOI: 10.1142/s0219720012500187] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Flap dynamics of HIV-1 protease (HIV-pr) controls the entry of inhibitors and substrates to the active site. Dynamical models from previous simulations are not all consistent with each other and not all are supported by the NMR results. In the present work, the effect of force field on the dynamics of HIV-pr is investigated by MD simulations using three AMBER force fields ff99, ff99SB, and ff03. The generalized order parameters for amide backbone are calculated from the three force fields and compared with the NMR S2 values. We found that the ff99SB and ff03 force field calculated order parameters agree reasonably well with the NMR S2 values, whereas ff99 calculated values deviate most from the NMR order parameters. Stereochemical geometry of protein models from each force field also agrees well with the remarks from NMR S2 values. However, between ff99SB and ff03, there are several differences, most notably in the loop regions. It is found that these loops are, in general, more flexible in the ff03 force field. This results in a larger active site cavity in the simulation with the ff03 force field. The effect of this difference in computer-aided drug design against flexible receptors is discussed.
Collapse
Affiliation(s)
- BISWA RANJAN MEHER
- Computational Biology Research Laboratory, Department of Biotechnology, Indian Institute of Technology, Guwahati, Assam 781039, India
- Computational Chemistry Laboratory, Department of Natural Sciences, Albany State University, Albany, Georgia 31705, USA
| | - MATTAPARTHI VENKATA SATISH KUMAR
- Computational Biology Research Laboratory, Department of Biotechnology, Indian Institute of Technology, Guwahati, Assam 781039, India
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - SMRITI SHARMA
- Centre for Computational Biology and Bioinformatics, School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - PRADIPTA BANDYOPADHYAY
- Centre for Computational Biology and Bioinformatics, School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| |
Collapse
|
30
|
Oehme DP, Brownlee RTC, Wilson DJD. Effect of atomic charge, solvation, entropy, and ligand protonation state on MM-PB(GB)SA binding energies of HIV protease. J Comput Chem 2012; 33:2566-80. [DOI: 10.1002/jcc.23095] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2012] [Revised: 06/27/2012] [Accepted: 07/25/2012] [Indexed: 11/05/2022]
|
31
|
Ode H, Nakashima M, Kitamura S, Sugiura W, Sato H. Molecular dynamics simulation in virus research. Front Microbiol 2012; 3:258. [PMID: 22833741 PMCID: PMC3400276 DOI: 10.3389/fmicb.2012.00258] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2012] [Accepted: 07/02/2012] [Indexed: 01/24/2023] Open
Abstract
Virus replication in the host proceeds by chains of interactions between viral and host proteins. The interactions are deeply influenced by host immune molecules and anti-viral compounds, as well as by mutations in viral proteins. To understand how these interactions proceed mechanically and how they are influenced by mutations, one needs to know the structures and dynamics of the proteins. Molecular dynamics (MD) simulation is a powerful computational method for delineating motions of proteins at an atomic-scale via theoretical and empirical principles in physical chemistry. Recent advances in the hardware and software for biomolecular simulation have rapidly improved the precision and performance of this technique. Consequently, MD simulation is quickly extending the range of applications in biology, helping to reveal unique features of protein structures that would be hard to obtain by experimental methods alone. In this review, we summarize the recent advances in MD simulations in the study of virus–host interactions and evolution, and present future perspectives on this technique.
Collapse
Affiliation(s)
- Hirotaka Ode
- Clinical Research Center, National Hospital Organization Nagoya Medical Center Nagoya, Aichi, Japan
| | | | | | | | | |
Collapse
|
32
|
Acevedo O, Ambrose Z, Flaherty PT, Aamer H, Jain P, Sambasivarao SV. Identification of HIV inhibitors guided by free energy perturbation calculations. Curr Pharm Des 2012; 18:1199-216. [PMID: 22316150 DOI: 10.2174/138161212799436421] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2011] [Accepted: 12/06/2011] [Indexed: 01/14/2023]
Abstract
Free energy perturbation (FEP) theory coupled to molecular dynamics (MD) or Monte Carlo (MC) statistical mechanics offers a theoretically precise method for determining the free energy differences of related biological inhibitors. Traditionally requiring extensive computational resources and expertise, it is only recently that its impact is being felt in drug discovery. A review of computer-aided anti-HIV efforts employing FEP calculations is provided here that describes early and recent successes in the design of human immunodeficiency virus type 1 (HIV-1) protease and non-nucleoside reverse transcriptase inhibitors. In addition, our ongoing work developing and optimizing leads for small molecule inhibitors of cyclophilin A (CypA) is highlighted as an update on the current capabilities of the field. CypA has been shown to aid HIV-1 replication by catalyzing the cis/trans isomerization of a conserved Gly-Pro motif in the Nterminal domain of HIV-1 capsid (CA) protein. In the absence of a functional CypA, e.g., by the addition of an inhibitor such as cyclosporine A (CsA), HIV-1 has reduced infectivity. Our simulations of acylurea-based and 1-indanylketone-based CypA inhibitors have determined that their nanomolar and micromolar binding affinities, respectively, are tied to their ability to stabilize Arg55 and Asn102. A structurally novel 1-(2,6-dichlorobenzamido) indole core was proposed to maximize these interactions. FEP-guided optimization, experimental synthesis, and biological testing of lead compounds for toxicity and inhibition of wild-type HIV-1 and CA mutants have demonstrated a dose-dependent inhibition of HIV-1 infection in two cell lines. While the inhibition is modest compared to CsA, the results are encouraging.
Collapse
Affiliation(s)
- Orlando Acevedo
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA.
| | | | | | | | | | | |
Collapse
|
33
|
Binding stability of peptides derived from 1ALA residue and 7GLY residues to sites near active center of fluctuating papain. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.03.102] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
34
|
Meher BR, Wang Y. Interaction of I50V mutant and I50L/A71V double mutant HIV-protease with inhibitor TMC114 (darunavir): molecular dynamics simulation and binding free energy studies. J Phys Chem B 2012; 116:1884-900. [PMID: 22239286 DOI: 10.1021/jp2074804] [Citation(s) in RCA: 61] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In the present work, the binding of inhibitor TMC114 (darunavir) to wild-type (WT), single (I50V) as well as double (I50L/A71V) mutant HIV-proteases (HIV-pr) was investigated with all-atom molecular dynamics (MD) simulations as well as molecular mechanic-Poisson-Boltzmann surface area (MM-PBSA) calculation. For both the apo and complexed HIV-pr, many intriguing effects due to double mutant, I50L/A71V, are observed. For example, the flap-flap distance and the distance from the active site to the flap residues in the apo I50L/A71V-HIV-pr are smaller than those of WT- and I50V-HIV-pr, probably making the active site smaller in volume and closer movement of flaps. For the complexed HIV-pr with TMC114, the double mutant I50L/A71V shows a less curling of the flap tips and less flexibility than WT and the single mutant I50V. As for the other previous studies, the present results also show that the single mutant I50V decreases the binding affinity of I50V-HIV-pr to TMC, resulting in a drug resistance; whereas the double mutant I50L/A71V increases the binding affinity, and as a result of the stronger binding, the I50L/A71V may be well adapted by the TMC114. The energy decomposition analysis suggests that the increase of the binding for the double mutant I50L/A71V-HIV-pr can be mainly attributed to the increase in electrostatic energy by -5.52 kacl/mol and van der Waals by -0.42 kcal/mol, which are canceled out in part by the increase of polar solvation energy of 1.99 kcal/mol. The I50L/A71V mutant directly increases the binding affinity by approximately -0.88 (Ile50 to Leu50) and -0.90 (Ile50' to Leu50') kcal/mol, accounting 45% for the total gain of the binding affinity. Besides the direct effects from the residues Leu50 and Leu50', the residue Gly49' increases the binding affinity of I50L/A71V-HIV-pr to the inhibitor by -0.74 kcal/mol, to which the electrostatic interaction of Leu50's backbone contributes by -1.23 kcal/mol. Another two residues Ile84 and Ile47' also increase the binding affinity by -0.22 and -0.29 kcal/mol, respectively, which can be mainly attributed to van der Waals terms (ΔT(vdw) = -0.21 and -0.39 kcal/mol).
Collapse
Affiliation(s)
- Biswa Ranjan Meher
- Computational Chemistry Laboratory, Department of Natural Sciences, Albany State University, Albany, Georgia 31705, USA
| | | |
Collapse
|
35
|
Yang Y, Shen Y, Li S, Jin N, Liu H, Yao X. Molecular dynamics and free energy studies on Aurora kinase A and its mutant bound with MLN8054: insight into molecular mechanism of subtype selectivity. MOLECULAR BIOSYSTEMS 2012; 8:3049-60. [DOI: 10.1039/c2mb25217a] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|
36
|
The evolution of cefotaximase activity in the TEM β-lactamase. J Mol Biol 2011; 415:205-20. [PMID: 22075446 DOI: 10.1016/j.jmb.2011.10.041] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2011] [Revised: 10/18/2011] [Accepted: 10/25/2011] [Indexed: 11/21/2022]
Abstract
The development of a molecular-level understanding of drug resistance through β-lactamase is critical not only in designing newer-generation antibacterial agents but also in providing insight into the evolutionary mechanisms of enzymes in general. In the present study, we have evaluated the effect of four drug resistance mutations (A42G, E104K, G238S, and M182T) on the cefotaximase activity of the TEM-1 β-lactamase. Using computational methods, including docking and molecular mechanics calculations, we have been able to correctly identify the relative order of catalytic activities associated with these four single point mutants. Further analyses suggest that the changes in catalytic efficiency for mutant enzymes are correlated to structural changes within the binding site. Based on the energetic and structural analyses of the wild-type and mutant enzymes, structural rearrangement is suggested as a mechanism of evolution of drug resistance through TEM β-lactamase. The present study not only provides molecular-level insight into the effect of four drug resistance mutations on the structure and function of the TEM β-lactamase but also establishes a foundation for a future molecular-level analysis of complete evolutionary trajectory for this class of enzymes.
Collapse
|
37
|
[Bioinformatics studies on drug resistance against anti-HIV-1 drugs]. Uirusu 2011; 61:35-47. [PMID: 21972554 DOI: 10.2222/jsv.61.35] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
More than 20 drugs have been available for anti-HIV-1 treatment in Japan. Combination therapy with these drugs dramatically decreases in morbidity and mortality of AIDS. However, due to high mutation rate of HIV-1, treatment with ineffective drugs toward patients infected with HIV-1 causes accumulation of mutations in the virus, and emergence of drug resistant viruses. Thus, to achieve appropriate application of the drugs toward the respective patients living with HIV-1, methods for predicting the level of drug-resistance using viral sequence information has been developed on the basis of bioinformatics. Furthermore, ultra-deep sequencing by next-generation sequencer whose data analysis is also based on bioinformatics, or in silico structural modeling have been achieved to understand drug resistant mechanisms. In this review, I overview the bioinformatics studies about drug resistance against anti-HIV-1 drugs.
Collapse
|
38
|
Lee J, Goodey NM. Catalytic contributions from remote regions of enzyme structure. Chem Rev 2011; 111:7595-624. [PMID: 21923192 DOI: 10.1021/cr100042n] [Citation(s) in RCA: 68] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Jeeyeon Lee
- Department of Chemistry, 413 Wartik Laboratory, The Pennsylvania State University, University Park, Pennsylvania 16802, USA.
| | | |
Collapse
|
39
|
Abstract
Accurate computational methods for predicting electrostatic energies are of major importance for our understanding of protein energetics in general for computer-aided drug design as well as for the design of novel biocatalysts and protein therapeutics. Electrostatic energies are of particular importance in such applications as virtual screening, drug design and protein-protein docking due to the high charge density of protein ligands and small-molecule drugs, and the frequent protonation state changes observed when drugs bind to their protein targets. Therefore, the development of a reliable and fast algorithm for the evaluation of electrostatic free energies, as an important contributor to the overall protein energy function, has been the focus for many scientists over the past three decades. In this review we describe the current state-of-the-art in modeling electrostatic effects in proteins and protein-ligand complexes. We focus mainly on the merits and drawbacks of the continuum methodology, and speculate on future directions in refining algorithms for calculating electrostatic energies in proteins using experimental data.
Collapse
|
40
|
Nishiyama K. Specific bindings of glycine peptides of distinctly different chain length on dynamic papain surfaces. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.05.022] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
41
|
Oehme DP, Wilson DJD, Brownlee RTC. Effect of structural stress on the flexibility and adaptability of HIV-1 protease. J Chem Inf Model 2011; 51:1064-73. [PMID: 21500830 DOI: 10.1021/ci2000677] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Resistance remains a major issue with regards to HIV-1 protease, despite the availability of numerous HIV-1 protease inhibitors and copious amounts of structural and binding data. In an effort to improve our understanding of how HIV-1 protease is able to "outsmart" new drugs, we have investigated the flexibility of HIV-1 protease and in particular how it adapts to different structural stresses. Our analysis has highlighted the effects of space group on the variability between structures of HIV-1 protease and suggests that consideration of multiple structures and appropriate consideration of different conformations of the Ile50 residue is necessary in any structural analysis. Calculation of the root-mean-square deviation on a per-residue basis has been used to identify 'natural variation', while mutational and ligand analyses have been carried out to identify the effect on structure as a result of specific stresses. It was observed that mutations readily cause changes to occur at sites both close to and distant from a mutation site, with changes more likely to occur at residues that are sites of other major mutations. It is also revealed that HIV-1 protease adaption is dependent on the type and the structure of any bound ligand. Identification of the specific changes that occur due to these stresses will aid in the understanding of resistance and also aid in the design of new drugs.
Collapse
Affiliation(s)
- Daniel P Oehme
- Department of Chemistry and La Trobe Institute for Molecular Sciences (LIMS), La Trobe University , Bundoora, Victoria, 3086, Australia
| | | | | |
Collapse
|
42
|
|
43
|
A comparative study of HIV-1 and HTLV-I protease structure and dynamics reveals a conserved residue interaction network. J Mol Model 2011; 17:2693-705. [PMID: 21279524 DOI: 10.1007/s00894-011-0971-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2010] [Accepted: 01/11/2011] [Indexed: 12/14/2022]
Abstract
The two retroviruses human T-lymphotropic virus type I (HTLV-I) and human immunodeficiency virus type 1 (HIV-1) are the causative agents of severe and fatal diseases including adult T-cell leukemia and the acquired immune deficiency syndrome (AIDS). Both viruses code for a protease that is essential for replication and therefore represents a key target for drugs interfering with viral infection. The retroviral proteases from HIV-1 and HTLV-I share 31% sequence identity and high structural similarities. Yet, their substrate specificities and inhibition profiles differ substantially. In this study, we performed all-atom molecular dynamics (MD) simulations for both enzymes in their ligand-free states and in complex with model substrates in order to compare their dynamic behaviors and enhance our understanding of the correlation between sequence, structure, and dynamics in this protein family. We found extensive similarities in both local and overall protein dynamics, as well as in the energetics of their interactions with model substrates. Interestingly, those residues that are important for strong ligand binding are frequently not conserved in sequence, thereby offering an explanation for the differences in binding specificity. Moreover, we identified an interaction network of contacts between conserved residues that interconnects secondary structure elements and serves as a scaffold for the protein fold. This interaction network is conformationally stable over time and may provide an explanation for the highly similar dynamic behavior of the two retroviral proteases, even in the light of their rather low overall sequence identity.
Collapse
|
44
|
Fanfrlík J, Bronowska AK, Rezác J, Prenosil O, Konvalinka J, Hobza P. A reliable docking/scoring scheme based on the semiempirical quantum mechanical PM6-DH2 method accurately covering dispersion and H-bonding: HIV-1 protease with 22 ligands. J Phys Chem B 2011; 114:12666-78. [PMID: 20839830 DOI: 10.1021/jp1032965] [Citation(s) in RCA: 96] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
In this study, we introduce a fast and reliable rescoring scheme for docked complexes based on a semiempirical quantum mechanical PM6-DH2 method. The method utilizes a PM6-based Hamiltonian with corrections for dispersion energy and hydrogen bonds. The total score is constructed as the sum of the PM6-DH2 interaction enthalpy, the empirical force field (AMBER) interaction entropy, and the sum of the deformation (PM6-DH2, SMD) and the desolvation (SMD) energies of the ligand. The main advantage of the procedure is the fact that we do not add any empirical parameter for either an individual component of the total score or an individual protein-ligand complex. This rescoring method is applied to a very challenging system, namely, the HIV-1 protease with a set of ligands. As opposed to the conventional DOCK procedure, the PM6-DH2 rescoring based on all of the terms distinguishes between binders and nonbinders and provides a reliable correlation of the theoretical and experimental binding free energies. Such a dramatic improvement, resulting from the PM6-DH2 rescoring of all the complexes, provides a valuable yet inexpensive tool for rational drug discovery and de novo ligand design.
Collapse
Affiliation(s)
- Jindrich Fanfrlík
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic and Center for Biomolecules and Complex Molecular Systems, 166 10 Prague 6, Czech Republic
| | | | | | | | | | | |
Collapse
|
45
|
Substrate specificity of papain dynamic structures for peptides consisting of 8–10 GLY residues. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2010.11.057] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
46
|
Hao GF, Yang GF, Zhan CG. Computational mutation scanning and drug resistance mechanisms of HIV-1 protease inhibitors. J Phys Chem B 2010; 114:9663-76. [PMID: 20604558 DOI: 10.1021/jp102546s] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The drug resistance of various clinically available HIV-1 protease inhibitors has been studied using a new computational protocol, that is, computational mutation scanning (CMS), leading to valuable insights into the resistance mechanisms and structure-resistance correction of the HIV-1 protease inhibitors associated with a variety of active site and nonactive site mutations. By using the CMS method, the calculated mutation-caused shifts of the binding free energies linearly correlate very well with those derived from the corresponding experimental data, suggesting that the CMS protocol may be used as a generalized approach to predict drug resistance associated with amino acid mutations. Because it is essentially important for understanding the structure-resistance correlation and for structure-based drug design to develop an effective computational protocol for drug resistance prediction, the reasonable and computationally efficient CMS protocol for drug resistance prediction should be valuable for future structure-based design and discovery of antiresistance drugs in various therapeutic areas.
Collapse
Affiliation(s)
- Ge-Fei Hao
- Key Laboratory of Pesticide & Chemical Biology of Ministry of Education, College of Chemistry, Central China Normal University, Wuhan 430079, PR China
| | | | | |
Collapse
|
47
|
Sano E, Li W, Yuki H, Liu X, Furihata T, Kobayashi K, Chiba K, Neya S, Hoshino T. Mechanism of the decrease in catalytic activity of human cytochrome P450 2C9 polymorphic variants investigated by computational analysis. J Comput Chem 2010; 31:2746-58. [DOI: 10.1002/jcc.21568] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
48
|
Hu G, Wang D, Liu X, Zhang Q. A computational analysis of the binding model of MDM2 with inhibitors. J Comput Aided Mol Des 2010; 24:687-97. [PMID: 20490618 PMCID: PMC2907675 DOI: 10.1007/s10822-010-9366-0] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2010] [Accepted: 05/06/2010] [Indexed: 12/13/2022]
Abstract
It is a new and promising strategy for anticancer drug design to block the MDM2-p53 interaction using a non-peptide small-molecule inhibitor. We carry out molecular dynamics simulations to study the binding of a set of six non-peptide small-molecule inhibitors with the MDM2. The relative binding free energies calculated using molecular mechanics Poisson-Boltzmann surface area method produce a good correlation with experimentally determined results. The study shows that the van der Waals energies are the largest component of the binding free energy for each complex, which indicates that the affinities of these inhibitors for MDM2 are dominated by shape complementarity. The A-ligands and the B-ligands are the same except for the conformation of 2,2-dimethylbutane group. The quantum mechanics and the binding free energies calculation also show the B-ligands are the more possible conformation of ligands. Detailed binding free energies between inhibitors and individual protein residues are calculated to provide insights into the inhibitor-protein binding model through interpretation of the structural and energetic results from the simulations. The study shows that G1, G2 and G3 group mimic the Phe19, Trp23 and Leu26 residues in p53 and their interactions with MDM2, but the binding model of G4 group differs from the original design strategy to mimic Leu22 residue in p53.
Collapse
Affiliation(s)
- Guodong Hu
- College of Physics and Electronics, Shandong Normal University, 250014 Jinan, China
- Department of Physics, University of Texas at San Antonio, San Antonio, TX 78249, USA
| | - Dunyou Wang
- College of Physics and Electronics, Shandong Normal University, 250014 Jinan, China
| | - Xinguo Liu
- College of Physics and Electronics, Shandong Normal University, 250014 Jinan, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, 250014 Jinan, China,
| |
Collapse
|
49
|
Li D, Ji B, Hwang K, Huang Y. Crucial roles of the subnanosecond local dynamics of the flap tips in the global conformational changes of HIV-1 protease. J Phys Chem B 2010; 114:3060-9. [PMID: 20143801 DOI: 10.1021/jp1005549] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
To understand the underlying mechanisms of the open and closed conformational change of HIV-1 protease (HIV-1 PR) at multiple time scales, we performed serial fully unrestrained, extremely long time molecular dynamics simulations with an explicit solvent model. Spontaneous semiopen to closed conformational transition and inhibitor-collision-induced opening of the flaps were simulated in a real time scale. We found that the rapid, local subnanosecond fluctuations of the flap tips might be the mechanisms triggering the global open and close conformational transitions at the 100-ns time scale. The subnanosecond fluctuation is induced by the Phi-Psi rotations of the residues at the flap tips, mainly Psi of Gly49 and Phi of Ile50, disturbing the interactions between the two tips and then their stability. We further showed that the water molecule W301 is helpful for the stability of the PR-inhibitor complex by acting as a collision buffer for the dynamic interaction between flap tips and the inhibitor. These results might help gain a better insight into the dynamics of HIV-1 PR, especially the local dynamics of the flap tips, which may provide important guidelines for design of novel potent inhibitors.
Collapse
Affiliation(s)
- Dechang Li
- Department of Engineering Mechanics, School of Aerospace, Tsinghua University, Beijing 100084, China
| | | | | | | |
Collapse
|
50
|
Ode H, Yokoyama M, Kanda T, Sato H. Identification of folding preferences of cleavage junctions of HIV-1 precursor proteins for regulation of cleavability. J Mol Model 2010; 17:391-9. [PMID: 20480379 DOI: 10.1007/s00894-010-0739-z] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2010] [Accepted: 04/30/2010] [Indexed: 11/30/2022]
Abstract
Human immunodeficiency virus type 1 protease (HIV-1 PR) cleaves two viral precursor proteins, Gag and Gag-Pol, at multiple sites. Although the processing proceeds in the rank order to assure effective viral replication, the molecular mechanisms by which the order is regulated are not fully understood. In this study, we used bioinformatics approaches to examine whether the folding preferences of the cleavage junctions influence their cleavabilities by HIV-1 PR. The folding of the eight-amino-acid peptides corresponding to the seven cleavage junctions of the HIV-1(HXB2) Gag and Gag-Pol precursors were simulated in the PR-free and PR-bound states with molecular dynamics and homology modeling methods, and the relationships between the folding parameters and the reported kinetic parameters of the HIV-1(HXB2) peptides were analyzed. We found that a folding preference for forming a dihedral angle of Cβ (P1)-Cα (P1)- Cα (P1')-Cβ (P1') in the range of 150 to 180 degrees in the PR-free state was positively correlated with the 1/K(m) (R = 0.95, P = 0.0008) and that the dihedral angle of the O (P2)-C (P2)- C (P1)- O (P1) of the main chains in the PR-bound state was negatively correlated with k(cat) (R = 0.94, P = 0.001). We further found that these two folding properties influenced the overall cleavability of the precursor protein when the sizes of the side chains at the P1 site were similar. These data suggest that the dihedral angles at the specific positions around the cleavage junctions before and after binding to PR are both critical for regulating the cleavability of precursor proteins by HIV-1 PR.
Collapse
Affiliation(s)
- Hirotaka Ode
- Pathogen Genomics Center, National Institute of Infectious Diseases, Tokyo, Japan.
| | | | | | | |
Collapse
|