1
|
Handa Y, Okuwaki K, Kawashima Y, Hatada R, Mochizuki Y, Komeiji Y, Tanaka S, Furuishi T, Yonemochi E, Honma T, Fukuzawa K. Prediction of Binding Pose and Affinity of Nelfinavir, a SARS-CoV-2 Main Protease Repositioned Drug, by Combining Docking, Molecular Dynamics, and Fragment Molecular Orbital Calculations. J Phys Chem B 2024; 128:2249-2265. [PMID: 38437183 PMCID: PMC10946393 DOI: 10.1021/acs.jpcb.3c05564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 02/09/2024] [Accepted: 02/14/2024] [Indexed: 03/06/2024]
Abstract
A novel in silico drug design procedure is described targeting the Main protease (Mpro) of the SARS-CoV-2 virus. The procedure combines molecular docking, molecular dynamics (MD), and fragment molecular orbital (FMO) calculations. The binding structure and properties of Mpro were predicted for Nelfinavir (NFV), which had been identified as a candidate compound through drug repositioning, targeting Mpro. Several poses of the Mpro and NFV complexes were generated by docking, from which four docking poses were selected by scoring with FMO energy. Then, each pose was subjected to MD simulation, 100 snapshot structures were sampled from each of the generated MD trajectories, and the structures were evaluated by FMO calculations to rank the pose based on binding energy. Several residues were found to be important in ligand recognition, including Glu47, Asp48, Glu166, Asp187, and Gln189, all of which interacted strongly with NFV. Asn142 is presumably regarded to form hydrogen bonds or CH/π interaction with NFV; however, in the present calculation, their interactions were transient. Moreover, the tert-butyl group of NFV had no interaction with Mpro. Identifying such strong and weak interactions provides candidates for maintaining and substituting ligand functional groups and important suggestions for drug discovery using drug repositioning. Besides the interaction between NFV and the amino acid residues of Mpro, the desolvation effect of the binding pocket also affected the ranking order. A similar procedure of drug design was applied to Lopinavir, and the calculated interaction energy and experimental inhibitory activity value trends were consistent. Our approach provides a new guideline for structure-based drug design starting from a candidate compound whose complex crystal structure has not been obtained.
Collapse
Affiliation(s)
- Yuma Handa
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Koji Okuwaki
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Yusuke Kawashima
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Ryo Hatada
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Yuji Mochizuki
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Institute
of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Yuto Komeiji
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
- Department
of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Health
and Medical Research Institute, AIST, Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
- RIKEN
Center
for Biosystems Dynamics Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Shigenori Tanaka
- Graduate
School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Takayuki Furuishi
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Etsuo Yonemochi
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Teruki Honma
- RIKEN
Center
for Biosystems Dynamics Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Kaori Fukuzawa
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
- Department
of Biomolecular Engineering, Graduate School of Engineering, Tohoku University, 6-6-11 Aoba, Aramaki, Aoba-ku, Sendai 980-8579, Japan
| |
Collapse
|
2
|
Puławski W, Koliński A, Koliński M. Integrative modeling of diverse protein-peptide systems using CABS-dock. PLoS Comput Biol 2023; 19:e1011275. [PMID: 37405984 DOI: 10.1371/journal.pcbi.1011275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2023] [Accepted: 06/15/2023] [Indexed: 07/07/2023] Open
Abstract
The CABS model can be applied to a wide range of protein-protein and protein-peptide molecular modeling tasks, such as simulating folding pathways, predicting structures, docking, and analyzing the structural dynamics of molecular complexes. In this work, we use the CABS-dock tool in two diverse modeling tasks: 1) predicting the structures of amyloid protofilaments and 2) identifying cleavage sites in the peptide substrates of proteolytic enzymes. In the first case, simulations of the simultaneous docking of amyloidogenic peptides indicated that the CABS model can accurately predict the structures of amyloid protofilaments which have an in-register parallel architecture. Scoring based on a combination of symmetry criteria and estimated interaction energy values for bound monomers enables the identification of protofilament models that closely match their experimental structures for 5 out of 6 analyzed systems. For the second task, it has been shown that CABS-dock coarse-grained docking simulations can be used to identify the positions of cleavage sites in the peptide substrates of proteolytic enzymes. The cleavage site position was correctly identified for 12 out of 15 analyzed peptides. When combined with sequence-based methods, these docking simulations may lead to an efficient way of predicting cleavage sites in degraded proteins. The method also provides the atomic structures of enzyme-substrate complexes, which can give insights into enzyme-substrate interactions that are crucial for the design of new potent inhibitors.
Collapse
Affiliation(s)
- Wojciech Puławski
- Bioinformatics Laboratory, Mossakowski Medical Research Institute, Polish Academy of Sciences, Warsaw, Poland
| | | | - Michał Koliński
- Bioinformatics Laboratory, Mossakowski Medical Research Institute, Polish Academy of Sciences, Warsaw, Poland
| |
Collapse
|
3
|
Chan HT, Oliveira ASF, Schofield CJ, Mulholland AJ, Duarte F. Dynamical Nonequilibrium Molecular Dynamics Simulations Identify Allosteric Sites and Positions Associated with Drug Resistance in the SARS-CoV-2 Main Protease. JACS AU 2023; 3:1767-1774. [PMID: 37384148 PMCID: PMC10262681 DOI: 10.1021/jacsau.3c00185] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Revised: 05/19/2023] [Accepted: 05/19/2023] [Indexed: 06/30/2023]
Abstract
The SARS-CoV-2 main protease (Mpro) plays an essential role in the coronavirus lifecycle by catalyzing hydrolysis of the viral polyproteins at specific sites. Mpro is the target of drugs, such as nirmatrelvir, though resistant mutants have emerged that threaten drug efficacy. Despite its importance, questions remain on the mechanism of how Mpro binds its substrates. Here, we apply dynamical nonequilibrium molecular dynamics (D-NEMD) simulations to evaluate structural and dynamical responses of Mpro to the presence and absence of a substrate. The results highlight communication between the Mpro dimer subunits and identify networks, including some far from the active site, that link the active site with a known allosteric inhibition site, or which are associated with nirmatrelvir resistance. They imply that some mutations enable resistance by altering the allosteric behavior of Mpro. More generally, the results show the utility of the D-NEMD technique for identifying functionally relevant allosteric sites and networks including those relevant to resistance.
Collapse
Affiliation(s)
- H. T.
Henry Chan
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| | - A. Sofia F. Oliveira
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, UK
- School
of Biochemistry, University of Bristol, Bristol BS8 1TD, UK
| | - Christopher J. Schofield
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| | - Adrian J. Mulholland
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, UK
| | - Fernanda Duarte
- Chemistry
Research Laboratory, Department of Chemistry and the Ineos Oxford
Institute for Antimicrobial Research, University
of Oxford, 12 Mansfield Road, Oxford OX1 3TA, UK
| |
Collapse
|
4
|
Cheepurupalli L, Diaz A, Gopal AC, Rathore SS, Ramakrishnan V, Ramakrishnan J. In vitro and in silico screening of Klebsiella pneumoniae new Delhi metallo-β-lactamase-1 inhibitors from endophytic Streptomyces spp. J Biomol Struct Dyn 2022; 40:13593-13605. [PMID: 34657563 DOI: 10.1080/07391102.2021.1990132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The increase in drug resistance over the last two decades is a big threat in health care settings. More importantly, the dissemination of carbapenem-resistant Enterobacteriaceae is the major threat to public health with an increase in morbidity and mortality. β-lactamase is known to confer enteric bacteria with nearly complete resistance to all β-lactam antibiotics including the late-generation carbapenems. The commercially available β-lactamase inhibitors, clavulanic acid, sulbactam, and tazobactam are being met with an increasing number of resistant phenotypes and are ineffective against pathogens harbouring New Delhi metallo-β-lactamase (NDM-1). Inhibition of New Delhi metallo-β-lactamase-1 activity is one potential way to treat metallo β-lactamase (MBL) producing multi drug resistant (MDR) pathogen. The present study focused on screening of Klebsiella pneumoniae New Delhi metallo-β-lactamase-1 (BLIs) from endophytic Streptomyces spp. using in vitro and in silico methods. The study identified three potential inhibitors of New Delhi metallo-β-lactamase-1, namely dodecanoic acid, dl-alanyl-l-leucine and phenyl propanedioic acid. These molecules were found to bind to other MBLs namely, IMP-1 and VIM-2. To the best of our knowledge, this is the first kind of study reporting the binding mode of these molecules with New Delhi metallo-β-lactamase-1.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Lalitha Cheepurupalli
- Centre for Research in Infectious Diseases (CRID), School of Chemical and Biotechnology (SCBT), SASTRA Deemed University, Tirumalaisamudram, Thanjavur, Tamil Nadu, India
| | - Aathithya Diaz
- Computational Molecular Biophysics Laboratory (CMBL), School of Chemical and Biotechnology (SCBT), SASTRA Deemed University, Tirumalaisamudram, Thanjavur, Tamil Nadu, India
| | - Adithya Conjeevaram Gopal
- Centre for Research in Infectious Diseases (CRID), School of Chemical and Biotechnology (SCBT), SASTRA Deemed University, Tirumalaisamudram, Thanjavur, Tamil Nadu, India
| | - Sudarshan Singh Rathore
- Centre for Research in Infectious Diseases (CRID), School of Chemical and Biotechnology (SCBT), SASTRA Deemed University, Tirumalaisamudram, Thanjavur, Tamil Nadu, India
| | - Vigneshwar Ramakrishnan
- School of Chemical and Biotechnology (SCBT), SASTRA Deemed University, Tirumalaisamudram, Thanjavur, Tamil Nadu, India
| | - Jayapradha Ramakrishnan
- Centre for Research in Infectious Diseases (CRID), School of Chemical and Biotechnology (SCBT), SASTRA Deemed University, Tirumalaisamudram, Thanjavur, Tamil Nadu, India
| |
Collapse
|
5
|
Yu YX, Wang W, Sun HB, Zhang LL, Wang LF, Yin YY. Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2022; 33:805-831. [PMID: 36322686 DOI: 10.1080/1062936x.2022.2140708] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Accepted: 10/20/2022] [Indexed: 06/16/2023]
Abstract
Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L F Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Y Y Yin
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
6
|
Masoumzadeh E, Grozdanov PN, Jetly A, MacDonald CC, Latham MP. Electrostatic Interactions between CSTF2 and pre-mRNA Drive Cleavage and Polyadenylation. Biophys J 2022; 121:607-619. [PMID: 35090899 PMCID: PMC8873925 DOI: 10.1016/j.bpj.2022.01.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2021] [Revised: 12/20/2021] [Accepted: 01/07/2022] [Indexed: 11/25/2022] Open
Abstract
Nascent pre-mRNA 3'-end cleavage and polyadenylation (C/P) involves numerous proteins that recognize multiple RNA elements. Human CSTF2 binds to a downstream U- or G/U-rich sequence through its RNA recognition motif (RRM) regulating C/P. We previously reported the only known disease-related CSTF2 RRM mutant (CSTF2D50A) and showed that it changed the on-rate of RNA binding, leading to alternative polyadenylation in brains of mice carrying the same mutation. In this study, we further investigated the role of electrostatic interactions in the thermodynamics and kinetics of RNA binding for the CSTF2 RRM and the downstream consequences for regulation of C/P. By combining mutagenesis with NMR spectroscopy and biophysical assays, we confirmed that electrostatic attraction is the dominant factor in RRM binding to a naturally occurring U-rich RNA sequence. Moreover, we demonstrate that RNA binding is accompanied by an enthalpy-entropy compensation mechanism that is supported by changes in pico-to-nanosecond timescale RRM protein dynamics. We suggest that the dynamic binding of the RRM to U-rich RNA supports the diversity of sequences it encounters in the nucleus. Lastly, in vivo C/P assays demonstrate a competition between fast, high affinity RNA binding and efficient, correct C/P. These results highlight the importance of the surface charge of the RRM in RNA binding and the balance between nascent mRNA binding and C/P in vivo.
Collapse
|
7
|
Abstract
Viral proteases are diverse in structure, oligomeric state, catalytic mechanism, and substrate specificity. This chapter focuses on proteases from viruses that are relevant to human health: human immunodeficiency virus subtype 1 (HIV-1), hepatitis C (HCV), human T-cell leukemia virus type 1 (HTLV-1), flaviviruses, enteroviruses, and coronaviruses. The proteases of HIV-1 and HCV have been successfully targeted for therapeutics, with picomolar FDA-approved drugs currently used in the clinic. The proteases of HTLV-1 and the other virus families remain emerging therapeutic targets at different stages of the drug development process. This chapter provides an overview of the current knowledge on viral protease structure, mechanism, substrate recognition, and inhibition. Particular focus is placed on recent advances in understanding the molecular basis of diverse substrate recognition and resistance, which is essential toward designing novel protease inhibitors as antivirals.
Collapse
Affiliation(s)
- Jacqueto Zephyr
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, United States
| | - Nese Kurt Yilmaz
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, United States
| | - Celia A Schiffer
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, United States.
| |
Collapse
|
8
|
Foloppe N, Chen IJ. The reorganization energy of compounds upon binding to proteins, from dynamic and solvated bound and unbound states. Bioorg Med Chem 2021; 51:116464. [PMID: 34798378 DOI: 10.1016/j.bmc.2021.116464] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 09/26/2021] [Accepted: 10/01/2021] [Indexed: 11/30/2022]
Abstract
The intramolecular reorganization energy (ΔEReorg) of compounds upon binding to proteins is a component of the binding free energy, which has long received particular attention, for fundamental and practical reasons. Understanding ΔEReorg would benefit the science of molecular recognition and drug design. For instance, the tolerable strain energy of compounds upon binding has been elusive. Prior studies found some large ΔEReorg values (e.g. > 10 kcal/mol), received with skepticism since they imply excessive opposition to binding. Indeed, estimating ΔEReorg is technically difficult. Typically, ΔEReorg has been approached by taking two energy-minimized conformers representing the bound and unbound states, and subtracting their conformational energy. This is a drastic oversimplification, liable to conformational collapse of the unbound conformer. Instead, the present work applies extensive molecular dynamics (MD) and the modern OPLS3 force-field to simulate compounds bound and unbound states, in explicit solvent under physically relevant conditions. The thermalized unbound compounds populate multiple conformations, not reducible to one or a few energy-minimized conformers. The intramolecular energies in the bound and unbound states were averaged over pertinent conformational ensembles, and the reorganization enthalpy upon binding (ΔHReorg) deduced by subtraction. This was applied to 76 systems, including 43 approved drugs, carefully selected for i) the quality of the bioactive X-ray structures and ii) the diversity of the chemotypes, their properties and protein targets. It yielded comparatively low ΔHReorg values (median = 1.4 kcal/mol, mean = 3.0 kcal/mol). A new finding is the observation of negative ΔHReorg values. Indeed, reorganization energies do not have to oppose binding, e.g. when intramolecular interactions stabilize preferentially the bound state. Conversely, even with competing water molecules, intramolecular interactions can occur predominantly in the unbound compound, and be replaced by intermolecular counterparts upon protein binding. Such disruption of intramolecular interactions upon binding gives rise to occasional larger ΔHReorg values. Such counterintuitive larger ΔHReorg values may be rationalized as a redistribution of interactions upon binding, qualitatively compatible with binding.
Collapse
Affiliation(s)
- Nicolas Foloppe
- Vernalis (R&D) Ltd, Granta Park, Abington, Cambridge CB21 6GB, UK.
| | - I-Jen Chen
- Vernalis (R&D) Ltd, Granta Park, Abington, Cambridge CB21 6GB, UK.
| |
Collapse
|
9
|
Neela YI, Guruprasad L. Structures and energetics of darunavir and active site amino acids of native and mutant HIV–1 protease: a computational study. Struct Chem 2021. [DOI: 10.1007/s11224-021-01852-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
10
|
Yu YX, Liu WT, Li HY, Wang W, Sun HB, Zhang LL, Wu SL. Decoding molecular mechanism underlying binding of drugs to HIV-1 protease with molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2021; 32:889-915. [PMID: 34551634 DOI: 10.1080/1062936x.2021.1979647] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2021] [Accepted: 09/08/2021] [Indexed: 06/13/2023]
Abstract
HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W T Liu
- Shuifa Qilu Cultural Tourism Development Co., Ltd, Shuifa Ecological Industry Group, Jinan, China
| | - H Y Li
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - S L Wu
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
11
|
Burnaman SH, Kneller DW, Wang YF, Kovalevsky A, Weber IT. Revertant mutation V48G alters conformational dynamics of highly drug resistant HIV protease PRS17. J Mol Graph Model 2021; 108:108005. [PMID: 34419931 DOI: 10.1016/j.jmgm.2021.108005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2021] [Revised: 08/06/2021] [Accepted: 08/08/2021] [Indexed: 11/27/2022]
Abstract
Drug resistance is a serious problem for controlling the HIV/AIDS pandemic. Current antiviral drugs show several orders of magnitude worse inhibition of highly resistant clinical variant PRS17 of HIV-1 protease compared with wild-type protease. We have analyzed the effects of a common resistance mutation G48V in the flexible flaps of the protease by assessing the revertant PRS17V48G for changes in enzyme kinetics, inhibition, structure, and dynamics. Both PRS17 and the revertant showed about 10-fold poorer catalytic efficiency than wild-type enzyme (0.55 and 0.39 μM-1min-1 compared to 6.3 μM-1min-1). Clinical inhibitors, amprenavir and darunavir, showed 2-fold and 8-fold better inhibition, respectively, of the revertant than of PRS17, although the inhibition constants for PRS17V48G were still 25 to 1,200-fold worse than for wild-type protease. Crystal structures of inhibitor-free revertant and amprenavir complexes with revertant and PRS17 were solved at 1.3-1.5 Å resolution. The amprenavir complexes of PRS17V48G and PRS17 showed no significant differences in the interactions with inhibitor, although changes were observed in the conformation of Phe53 and the interactions of the flaps. The inhibitor-free structure of the revertant showed flaps in an open conformation, however, the flap tips do not have the unusual curled conformation seen in inhibitor-free PRS17. Molecular dynamics simulations were run for 1 μs on the two inhibitor-free mutants and wild-type protease. PRS17 exhibited higher conformational fluctuations than the revertant, while the wild-type protease adopted the closed conformation and showed the least variation. The second half of the simulations captured the transition of the flaps of PRS17 from a closed to a semi-open state, whereas the flaps of PRS17V48G tucked into the active site and the wild-type protease retained the closed conformation. These results suggest that mutation G48V contributes to drug resistance by altering the conformational dynamics of the flaps.
Collapse
Affiliation(s)
| | - Daniel W Kneller
- Neutron Scattering Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, TN, 37831, USA
| | - Yuan-Fang Wang
- Department of Biology, Georgia State University, Atlanta, GA, USA
| | - Andrey Kovalevsky
- Neutron Scattering Division, Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, TN, 37831, USA
| | - Irene T Weber
- Department of Biology, Georgia State University, Atlanta, GA, USA.
| |
Collapse
|
12
|
Shin YH, Park CM, Yoon CH. An Overview of Human Immunodeficiency Virus-1 Antiretroviral Drugs: General Principles and Current Status. Infect Chemother 2021; 53:29-45. [PMID: 34409780 PMCID: PMC8032919 DOI: 10.3947/ic.2020.0100] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Accepted: 02/22/2021] [Indexed: 12/11/2022] Open
Abstract
Treatment with highly active antiretroviral therapy (HAART) can prolong a patient's life-span by disrupting pivotal steps in the replication cycle of the human immunodeficiency virus-1 (HIV-1). However, drug resistance is emerging as a major problem worldwide due to the prolonged period of treatment undergone by HIV-1 patients. Since the approval of zidovudine in 1987, over thirty antiretroviral drugs have been categorized into the following six distinct classes based on their biological function and resistance profiles: (1) nucleoside analog reverse-transcriptase inhibitors; (2) non–nucleoside reverse transcriptase inhibitors; (3) integrase strand transferase inhibitors; (4) protease inhibitors; (5) fusion inhibitors; and (6) co-receptor antagonists. Additionally, several antiretroviral drugs have been developed recently, such as a long active drug, humanized antibody and pro-drug metabolized into an active form in the patient's body. Although plenty of antiretroviral drugs are beneficially used to treat patients with HIV-1, the ongoing efforts to develop antiretroviral drugs have overcome the drug resistances, adverse effects, and limited adherence of drugs observed in previous drugs to some extent. Furthermore, studies focused on agents targeting latent HIV-1 reservoirs should be strengthened, as that may lead to eradication of HIV-1.
Collapse
Affiliation(s)
- Young Hyun Shin
- Division of Chronic Viral Disease Research, Center for Emerging Virus Research, Korea National Institute of Health, Chungbuk, Korea
| | - Chul Min Park
- Center for Convergent Research of Emerging Virus Infection, Korea Research Institute of Chemical Technology, Daejeon, Korea
| | - Cheol Hee Yoon
- Division of Chronic Viral Disease Research, Center for Emerging Virus Research, Korea National Institute of Health, Chungbuk, Korea.
| |
Collapse
|
13
|
Ota R, So K, Tsuda M, Higuchi Y, Yamashita F. Prediction of HIV drug resistance based on the 3D protein structure: Proposal of molecular field mapping. PLoS One 2021; 16:e0255693. [PMID: 34347839 PMCID: PMC8336827 DOI: 10.1371/journal.pone.0255693] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2021] [Accepted: 07/21/2021] [Indexed: 11/19/2022] Open
Abstract
A method for predicting HIV drug resistance by using genotypes would greatly assist in selecting appropriate combinations of antiviral drugs. Models reported previously have had two major problems: lack of information on the 3D protein structure and processing of incomplete sequencing data in the modeling procedure. We propose obtaining the 3D structural information of viral proteins by using homology modeling and molecular field mapping, instead of just their primary amino acid sequences. The molecular field potential parameters reflect the physicochemical characteristics associated with the 3D structure of the proteins. We also introduce the Bayesian conditional mutual information theory to estimate the probabilities of occurrence of all possible protein candidates from an incomplete sequencing sample. This approach allows for the effective use of uncertain information for the modeling process. We applied these data analysis techniques to the HIV-1 protease inhibitor dataset and developed drug resistance prediction models with reasonable performance.
Collapse
Affiliation(s)
- Ryosaku Ota
- Department of Drug Delivery Research, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Kanako So
- Department of Applied Pharmaceutics and Pharmacokinetics, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Masahiro Tsuda
- Department of Applied Pharmaceutics and Pharmacokinetics, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Yuriko Higuchi
- Department of Drug Delivery Research, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
| | - Fumiyoshi Yamashita
- Department of Drug Delivery Research, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
- Department of Applied Pharmaceutics and Pharmacokinetics, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
- * E-mail:
| |
Collapse
|
14
|
Yu YX, Wang W, Sun HB, Zhang LL, Wu SL, Liu WT. Insights into effect of the Asp25/Asp25' protonation states on binding of inhibitors Amprenavir and MKP97 to HIV-1 protease using molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2021; 32:615-641. [PMID: 34157882 DOI: 10.1080/1062936x.2021.1939149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Accepted: 06/02/2021] [Indexed: 06/13/2023]
Abstract
The protonation states of two aspartic acids in the catalytic strands of HIV-1 protease (PR) remarkably affect bindings of inhibitors to PR. It is requisite for the design of potent inhibitors towards PR to investigate the influences of Asp25/Asp25' protonated states on dynamics behaviour of PR and binding mechanism of inhibitors to PR. In this work, molecular dynamics (MD) simulations, MM-GBSA method and principal component (PC) analysis were coupled to explore the effect of Asp25/Asp25' protonation states on conformational changes of PR and bindings of Amprenavir and MKP97 to PR. The results show that the Asp25/Asp25' protonation states exert different impacts on structural fluctuations, flexibility and motion modes of PR. Dynamics analysis verifies that Asp25/Asp25' protonated states highly affect conformational dynamics of two flaps in PR. The binding free energy calculations results suggest that the Asp25/Asp25' protonated states obviously strengthen bindings of inhibitors to PR compared to the non-protonation state. Calculations of residue-based free energy decomposition indicate that the Asp25/Asp25' protonation not only disturbs the interaction network of inhibitors with PR but also stabilizes bindings of inhibitors to PR by cancelling the electrostatic repulsive interaction. Therefore, special attentions should be paid to the Asp25/Asp25' protonation in the design of potent inhibitors towards PR.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - S L Wu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W T Liu
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
15
|
Peccati F, Jiménez-Osés G. Enthalpy-Entropy Compensation in Biomolecular Recognition: A Computational Perspective. ACS OMEGA 2021; 6:11122-11130. [PMID: 34056267 PMCID: PMC8153931 DOI: 10.1021/acsomega.1c00485] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2021] [Accepted: 04/07/2021] [Indexed: 05/08/2023]
Abstract
This mini-review provides an overview of the enthalpy-entropy compensation phenomenon in the simulation of biomacromolecular recognition, with particular emphasis on ligand binding. We approach this complex phenomenon from the point of view of practical computational chemistry. Without providing a detailed description of the plethora of existing methodologies already reviewed in depth elsewhere, we present a series of examples to illustrate different approaches to interpret and predict compensation phenomena at an atomistic level, which is far from trivial to predict using canonical, classic textbook assumptions.
Collapse
Affiliation(s)
- Francesca Peccati
- Center
for Cooperative Research in Biosciences (CIC bioGUNE), Basque Research
and Technology Alliance (BRTA), Bizkaia
Technology Park, Building
801A, 48160 Derio, Spain
| | - Gonzalo Jiménez-Osés
- Center
for Cooperative Research in Biosciences (CIC bioGUNE), Basque Research
and Technology Alliance (BRTA), Bizkaia
Technology Park, Building
801A, 48160 Derio, Spain
- Ikerbasque, Basque
Foundation for Science, 48013 Bilbao, Spain
| |
Collapse
|
16
|
Peng C, Wang J, Shi Y, Xu Z, Zhu W. Increasing the Sampling Efficiency of Protein Conformational Change by Combining a Modified Replica Exchange Molecular Dynamics and Normal Mode Analysis. J Chem Theory Comput 2020; 17:13-28. [PMID: 33351613 DOI: 10.1021/acs.jctc.0c00592] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Understanding conformational change at an atomic level is significant when determining a protein functional mechanism. Replica exchange molecular dynamics (REMD) is a widely used enhanced sampling method to explore protein conformational space. However, REMD with an explicit solvent model requires huge computational resources, immensely limiting its application. In this study, a variation of parallel tempering metadynamics (PTMetaD) with the omission of solvent-solvent interactions in exchange attempts and the use of low-frequency modes calculated by normal-mode analysis (NMA) as collective variables (CVs), namely ossPTMetaD, is proposed with the aim to accelerate MD simulations simultaneously in temperature and geometrical spaces. For testing the performance of ossPTMetaD, five protein systems with diverse biological functions and motion patterns were selected, including large-scale domain motion (AdK), flap movement (HIV-1 protease and BACE1), and DFG-motif flip in kinases (p38α and c-Abl). The simulation results showed that ossPTMetaD requires much fewer numbers of replicas than temperature REMD (T-REMD) with a reduction of ∼70% to achieve a similar exchange ratio. Although it does not obey the detailed balance condition, ossPTMetaD provides consistent results with T-REMD and experimental data. The high accessibility of the large conformational change of protein systems by ossPTMetaD, especially in simulating the very challenging DFG-motif flip of protein kinases, demonstrated its high efficiency and robustness in the characterization of the large-scale protein conformational change pathway and associated free energy profile.
Collapse
Affiliation(s)
- Cheng Peng
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| | - Jinan Wang
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China
| | - Yulong Shi
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| | - Zhijian Xu
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| | - Weiliang Zhu
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.,Open Studio for Druggability Research of Marine Lead Compounds, Qingdao National Laboratory for Marine Science and Technology, 1 Wenhai Road, Aoshanwei, Jimo, Qingdao 266237, China.,University of Chinese Academy of Sciences, No.19A Yuquan Road, Beijing, 100049, China
| |
Collapse
|
17
|
Khan SN, Persons JD, Guerrero M, Ilina TV, Oda M, Ishima R. A synergy of activity, stability, and inhibitor-interaction of HIV-1 protease mutants evolved under drug-pressure. Protein Sci 2020; 30:571-582. [PMID: 33314454 DOI: 10.1002/pro.4013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Revised: 12/04/2020] [Accepted: 12/06/2020] [Indexed: 01/01/2023]
Abstract
A clinically-relevant, drug-resistant mutant of HIV-1 protease (PR), termed Flap+(I54V) and containing L10I, G48V, I54V and V82A mutations, is known to produce significant changes in the entropy and enthalpy balance of drug-PR interactions, compared to wild-type PR. A similar mutant, Flap+(I54A) , which evolves from Flap+(I54V) and contains the single change at residue 54 relative to Flap+(I54V) , does not. Yet, how Flap+(I54A) behaves in solution is not known. To understand the molecular basis of V54A evolution, we compared nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy, isothermal titration calorimetry, and enzymatic assay data from four PR proteins: PR (pWT), Flap+(I54V) , Flap+(I54A) , and Flap+(I54) , a control mutant that contains only L10I, G48V and V82A mutations. Our data consistently show that selection to the smaller side chain at residue 54, not only decreases inhibitor affinity, but also restores the catalytic activity.
Collapse
Affiliation(s)
- Shahid N Khan
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, USA.,Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto, Japan
| | - John D Persons
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, USA
| | - Michel Guerrero
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, USA
| | - Tatiana V Ilina
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, USA
| | - Masayuki Oda
- Graduate School of Life and Environmental Sciences, Kyoto Prefectural University, Kyoto, Japan
| | - Rieko Ishima
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, Pennsylvania, USA
| |
Collapse
|
18
|
Affiliation(s)
- Mark Aldren M. Feliciano
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Brian Gold
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
19
|
Peng C, Wang J, Xu Z, Cai T, Zhu W. Accurate prediction of relative binding affinities of a series of HIV-1 protease inhibitors using semi-empirical quantum mechanical charge. J Comput Chem 2020; 41:1773-1780. [PMID: 32352193 DOI: 10.1002/jcc.26218] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2019] [Revised: 04/03/2020] [Accepted: 04/18/2020] [Indexed: 11/05/2022]
Abstract
A major challenge in computer-aided drug design is the accurate estimation of ligand binding affinity. Here, a new approach that combines the adaptive steered molecular dynamics (ASMD) and partial atomic charges calculated by semi-empirical quantum mechanics (SQMPC), namely ASMD-SQMPC, is suggested to predict the ligand binding affinities, with 24 HIV-1 protease inhibitors as testing examples. In the ASMD-SQMPC, the relative binding free energy (ΔG) is reflected by the average maximum potential of mean force (<PMF>max ) between bound and unbound states. The correlation coefficient (R2 ) between the <PMF>max and experimentally determined ΔG is 0.86, showing a significant improvement compared with the conventional ASMD (R2 = 0.52). Therefore, this study provides an efficient approach to predict the relative ΔG and reveals the significance of precise partial atomic charges in the theoretical simulations.
Collapse
Affiliation(s)
- Cheng Peng
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, China.,University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing, China
| | - Jinan Wang
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, China
| | - Zhijian Xu
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, China.,University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing, China
| | - Tingting Cai
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, China
| | - Weiliang Zhu
- CAS Key Laboratory of Receptor Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, China.,University of Chinese Academy of Sciences, No. 19A Yuquan Road, Beijing, China.,Open Studio for Druggability Research of Marine Natural Products, Pilot National Laboratory for Marine Science and Technology (Qingdao), 1 Wenhai Road, Aoshanwei, Jimo, Qingdao, China
| |
Collapse
|
20
|
Melo R, Lemos A, Preto AJ, Bueschbell B, Matos-Filipe P, Barreto C, Almeida JG, Silva RDM, Correia JDG, Moreira IS. An Overview of Antiretroviral Agents for Treating HIV Infection in Paediatric Population. Curr Med Chem 2020; 27:760-794. [PMID: 30182840 DOI: 10.2174/0929867325666180904123549] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2018] [Revised: 07/11/2018] [Accepted: 07/11/2018] [Indexed: 12/19/2022]
Abstract
Paediatric Acquired ImmunoDeficiency Syndrome (AIDS) is a life-threatening and infectious disease in which the Human Immunodeficiency Virus (HIV) is mainly transmitted through Mother-To- Child Transmission (MTCT) during pregnancy, labour and delivery, or breastfeeding. This review provides an overview of the distinct therapeutic alternatives to abolish the systemic viral replication in paediatric HIV-1 infection. Numerous classes of antiretroviral agents have emerged as therapeutic tools for downregulation of different steps in the HIV replication process. These classes encompass Non- Nucleoside Analogue Reverse Transcriptase Inhibitors (NNRTIs), Nucleoside/Nucleotide Analogue Reverse Transcriptase Inhibitors (NRTIs/NtRTIs), INtegrase Inhibitors (INIs), Protease Inhibitors (PIs), and Entry Inhibitors (EIs). Co-administration of certain antiretroviral drugs with Pharmacokinetic Enhancers (PEs) may boost the effectiveness of the primary therapeutic agent. The combination of multiple antiretroviral drug regimens (Highly Active AntiRetroviral Therapy - HAART) is currently the standard therapeutic approach for HIV infection. So far, the use of HAART offers the best opportunity for prolonged and maximal viral suppression, and preservation of the immune system upon HIV infection. Still, the frequent administration of high doses of multiple drugs, their inefficient ability to reach the viral reservoirs in adequate doses, the development of drug resistance, and the lack of patient compliance compromise the complete HIV elimination. The development of nanotechnology-based drug delivery systems may enable targeted delivery of antiretroviral agents to inaccessible viral reservoir sites at therapeutic concentrations. In addition, the application of Computer-Aided Drug Design (CADD) approaches has provided valuable tools for the development of anti-HIV drug candidates with favourable pharmacodynamics and pharmacokinetic properties.
Collapse
Affiliation(s)
- Rita Melo
- Centro de Ciencias e Tecnologias Nucleares, Instituto Superior Tecnico, Universidade de Lisboa, CTN, Estrada Nacional 10 (km 139,7), Bobadela LRS 2695-066, Portugal.,CNC - Center for Neuroscience and Cell Biology; Rua Larga, FMUC, Polo I, 1ºandar, Universidade de Coimbra, Coimbra 3004-517, Portugal
| | - Agostinho Lemos
- CNC - Center for Neuroscience and Cell Biology; Rua Larga, FMUC, Polo I, 1ºandar, Universidade de Coimbra, Coimbra 3004-517, Portugal.,GIGA Cyclotron Research Centre In Vivo Imaging, University of Liège, Liège 4000, Belgium
| | - António J Preto
- CNC - Center for Neuroscience and Cell Biology; Rua Larga, FMUC, Polo I, 1ºandar, Universidade de Coimbra, Coimbra 3004-517, Portugal
| | - Beatriz Bueschbell
- Pharmaceutical Chemistry I, PharmaCenter, Pharmaceutical Institute, University of Bonn, Bonn, Germany
| | - Pedro Matos-Filipe
- CNC - Center for Neuroscience and Cell Biology; Rua Larga, FMUC, Polo I, 1ºandar, Universidade de Coimbra, Coimbra 3004-517, Portugal
| | - Carlos Barreto
- CNC - Center for Neuroscience and Cell Biology; Rua Larga, FMUC, Polo I, 1ºandar, Universidade de Coimbra, Coimbra 3004-517, Portugal
| | - José G Almeida
- CNC - Center for Neuroscience and Cell Biology; Rua Larga, FMUC, Polo I, 1ºandar, Universidade de Coimbra, Coimbra 3004-517, Portugal
| | - Rúben D M Silva
- Centro de Ciencias e Tecnologias Nucleares, Instituto Superior Tecnico, Universidade de Lisboa, CTN, Estrada Nacional 10 (km 139,7), Bobadela LRS 2695-066, Portugal
| | - João D G Correia
- Centro de Ciencias e Tecnologias Nucleares, Instituto Superior Tecnico, Universidade de Lisboa, CTN, Estrada Nacional 10 (km 139,7), Bobadela LRS 2695-066, Portugal
| | - Irina S Moreira
- CNC - Center for Neuroscience and Cell Biology; Rua Larga, FMUC, Polo I, 1ºandar, Universidade de Coimbra, Coimbra 3004-517, Portugal.,Bijvoet Center for Biomolecular Research, Faculty of Science - Chemistry, Utrecht University, Utrecht 3584CH, Netherland
| |
Collapse
|
21
|
Triki D, Kermarrec M, Visseaux B, Descamps D, Flatters D, Camproux AC, Regad L. Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures. J Biomol Struct Dyn 2019; 38:5014-5026. [PMID: 31830870 DOI: 10.1080/07391102.2019.1704877] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Dhoha Triki
- Department of Chemistry, Bioinformatics Research Center, College of Science, North Carolina State University, Raleigh, North Carolina, USA.,Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| | - Maxime Kermarrec
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| | - Benoît Visseaux
- Université de Paris, IAME, UMR 1137, INSERM, AP-HP, Laboratoire de Virologie, Hôpital Bichat-Claude Bernard, Paris, France
| | - Diane Descamps
- Université de Paris, IAME, UMR 1137, INSERM, AP-HP, Laboratoire de Virologie, Hôpital Bichat-Claude Bernard, Paris, France
| | - Delphine Flatters
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| | | | - Leslie Regad
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| |
Collapse
|
22
|
Optimizing ligand conformations in flexible protein targets: a multi-objective strategy. Soft comput 2019. [DOI: 10.1007/s00500-019-04575-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
23
|
Ishima R, Kurt Yilmaz N, Schiffer CA. NMR and MD studies combined to elucidate inhibitor and water interactions of HIV-1 protease and their modulations with resistance mutations. JOURNAL OF BIOMOLECULAR NMR 2019; 73:365-374. [PMID: 31243634 PMCID: PMC6941145 DOI: 10.1007/s10858-019-00260-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2019] [Accepted: 06/19/2019] [Indexed: 06/09/2023]
Abstract
Over the last two decades, both the sensitivity of NMR and the time scale of molecular dynamics (MD) simulation have increased tremendously and have advanced the field of protein dynamics. HIV-1 protease has been extensively studied using these two methods, and has presented a framework for cross-evaluation of structural ensembles and internal dynamics by integrating the two methods. Here, we review studies from our laboratories over the last several years, to understand the mechanistic basis of protease drug-resistance mutations and inhibitor responses, using NMR and crystal structure-based parallel MD simulations. Our studies demonstrate that NMR relaxation experiments, together with crystal structures and MD simulations, significantly contributed to the current understanding of structural/dynamic changes due to HIV-1 protease drug resistance mutations.
Collapse
Affiliation(s)
- Rieko Ishima
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - Nese Kurt Yilmaz
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, USA
| | - Celia A Schiffer
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, USA.
| |
Collapse
|
24
|
Molecular alteration in drug susceptibility against subtype B and C-SA HIV-1 proteases: MD study. Struct Chem 2019. [DOI: 10.1007/s11224-019-01305-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
|
25
|
Characterizing early drug resistance-related events using geometric ensembles from HIV protease dynamics. Sci Rep 2018; 8:17938. [PMID: 30560871 PMCID: PMC6298995 DOI: 10.1038/s41598-018-36041-8] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2018] [Accepted: 11/14/2018] [Indexed: 02/07/2023] Open
Abstract
The use of antiretrovirals (ARVs) has drastically improved the life quality and expectancy of HIV patients since their introduction in health care. Several millions are still afflicted worldwide by HIV and ARV resistance is a constant concern for both healthcare practitioners and patients, as while treatment options are finite, the virus constantly adapts via complex mutation patterns to select for resistant strains under the pressure of drug treatment. The HIV protease is a crucial enzyme for viral maturation and has been a game changing drug target since the first application. Due to similarities in protease inhibitor designs, drug cross-resistance is not uncommon across ARVs of the same class. It is known that resistance against protease inhibitors is associated with a wider active site, but results from our large scale molecular dynamics simulations combined with statistical tests and network analysis further show, for the first time, that there are regions of local expansions and compactions associated with high levels of resistance conserved across eight different protease inhibitors visible in their complexed form within closed receptor conformations. The observed conserved expansion sites may provide an alternative drug-targeting site. Further, the method developed here is novel, supplementary to methods of variation analysis at sequence level, and should be applicable in analysing the structural consequences of mutations in other contexts using molecular ensembles.
Collapse
|
26
|
Chen J, Peng C, Wang J, Zhu W. Exploring molecular mechanism of allosteric inhibitor to relieve drug resistance of multiple mutations in HIV-1 protease by enhanced conformational sampling. Proteins 2018; 86:1294-1305. [PMID: 30260044 DOI: 10.1002/prot.25610] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2018] [Revised: 09/02/2018] [Accepted: 09/23/2018] [Indexed: 12/12/2022]
Abstract
Recently, allosteric regulations of HIV-1 protease (PR) are suggested as a promising approach to relieve drug resistance of mutations toward inhibitors targeting the active site of PR. Replica-exchange molecular dynamics (REMD) simulations and normal mode analysis (NMA) are integrated to enhance conformational sampling of PR. Molecular mechanics generalized Born surface area (MM-GBSA) method was applied to calculate binding free energies of three inhibitors APV, DRV, and NIT to the wild-type (WT) and multidrug resistance (MDR) PRs. The results suggest that binding free energies of APV and DRV are decreased in the MDR PR relative to the WT PR, suggesting drug resistance of mutations on these two inhibitors. However, the binding ability of the allosteric inhibitor NIT is not impaired in the MDR PR. In addition, internal dynamics analysis based on REMD simulations proves that mutations hardly produce obvious effect on the conformation of the MDR PR in comparison to the WT PR. Scanning of hydrophobic contacts and hydrogen bond contacts of inhibitors with residues of PRs on the concatenated trajectories of REMD demonstrates that mutations change the symmetric interaction networks of APV and DRV with PR, but do not generate obvious influence on the asymmetric interaction network of NIT with PR. In summary, allosteric inhibitor NIT can adapt the MDR PR better than those inhibitors toward the active site of PR, thus allosteric inhibitors of PR may be a possible channel to overcome drug resistance of PR.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China.,Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Cheng Peng
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China.,School of Chemistry, University of Chinese Academy of Sciences, Beijing, China
| | - Jinan Wang
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| | - Weiliang Zhu
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China
| |
Collapse
|
27
|
Yokoyama H, Sawada JI, Sato K, Ogo N, Kamei N, Ishikawa Y, Hara K, Asai A, Hashimoto H. Structural and Thermodynamic Basis of the Enhanced Interaction between Kinesin Spindle Protein Eg5 and STLC-type Inhibitors. ACS OMEGA 2018; 3:12284-12294. [PMID: 31459302 PMCID: PMC6644766 DOI: 10.1021/acsomega.8b00778] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/21/2018] [Accepted: 08/22/2018] [Indexed: 06/10/2023]
Abstract
For a better understanding of protein-inhibitor interactions, we report structural, thermodynamic, and biological analyses of the interactions between S-trityl-l-cysteine (STLC) derivatives and the motor domain of kinesin spindle protein Eg5. Binding of STLC-type inhibitors to Eg5 was enthalpically driven and entropically unfavorable. The introduction of a para-methoxy substituent in one phenyl ring of STLC enhances its inhibitory activity resulting from a larger enthalpy gain possibly due to the increased shape complementarity. The substituent fits to a recess in the binding pocket. To avoid steric hindrance, the substituted STLC is nudged toward the side opposite to the recess, which enhances the interaction of Eg5 with the remaining part of the inhibitor. Further introduction of an ethylene linkage between two phenyl rings enhances Eg5 inhibitory activity by reducing the loss of entropy in forming the complex. This study provides valuable examples of enhancing protein-inhibitor interactions without forming additional hydrogen bonds.
Collapse
Affiliation(s)
- Hideshi Yokoyama
- Faculty
of Pharmaceutical Sciences, Tokyo University
of Science, 2641 Yamazaki, Noda, Chiba 278-8510, Japan
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Jun-ichi Sawada
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Kohei Sato
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Naohisa Ogo
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Nanami Kamei
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Yoshinobu Ishikawa
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Kodai Hara
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Akira Asai
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| | - Hiroshi Hashimoto
- Department of Physical Biochemistry,
School of Pharmaceutical Sciences and Center for Drug
Discovery, Graduate School of Pharmaceutical Sciences, University of Shizuoka, 52-1 Yada, Suruga-ku, Shizuoka 422-8526, Japan
| |
Collapse
|
28
|
Maphumulo SI, Halder AK, Govender T, Maseko S, Maguire GEM, Honarparvar B, Kruger HG. Exploring the flap dynamics of the South African HIV subtype C protease in presence of FDA-approved inhibitors: MD study. Chem Biol Drug Des 2018; 92:1899-1913. [PMID: 30003668 DOI: 10.1111/cbdd.13364] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2017] [Accepted: 07/09/2018] [Indexed: 01/01/2023]
Abstract
HIV-1 protease (HIV PR) is considered as one of the most attractive targets for the treatment of HIV and the impact of flap dynamics of HIV PR on the binding affinities of protease inhibitors (PIs) is a crucial ongoing research field. Recently, our research group evaluated the binding affinities of different FDA approved PIs against the South African HIV-1 subtype C (C-SA) protease (PR). The CSA-HIV PR displayed weaker binding affinity for most of the clinical PIs compared to HIV-1 B subtype for West and Central Europe, the Americas. In the current work, the flap dynamics of four different systems of HIV-1 C-SA PR complexed to FDA approved second generation PIs and its impact on binding was explored over the molecular dynamic trajectories. It was observed that the interactions of the selected drugs with the binding site residues of the protease may not be the major contributor for affinity towards PIs. Various post-MD analyses were performed, also entropic contributions, solvation free energies and hydrophobic core formation interactions were studied to assess how the flap dynamics of C-SA PR which is affected by such factors. From these contributions, large van der Waals interactions and low solvation free energies were found to be major factors for the higher activity of ATV against C-SA HIV PR. Furthermore, a comparatively stable hydrophobic core may be responsible for higher stability of the PR flaps of the ATV complex. The outcome of this study provides significant guidance to how the flap dynamics of C-SA PR is affected by various factors as a result of the binding affinity of various protease inhibitors. It will also assist with the design of potent inhibitors against C-SA HIV PR that apart from binding in the active site of PR can interacts with the flaps to prevent opening of the flaps resulting in inactivation of the protease.
Collapse
Affiliation(s)
- Siyabonga I Maphumulo
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Amit K Halder
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Thavendran Govender
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Sibusiso Maseko
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Glenn E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa.,School of Chemistry and Physics, University of KwaZulu-Natal, Durban, South Africa
| | - Bahareh Honarparvar
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Hendrik G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| |
Collapse
|
29
|
Sun H, Duan L, Chen F, Liu H, Wang Z, Pan P, Zhu F, Zhang JZH, Hou T. Assessing the performance of MM/PBSA and MM/GBSA methods. 7. Entropy effects on the performance of end-point binding free energy calculation approaches. Phys Chem Chem Phys 2018; 20:14450-14460. [PMID: 29785435 DOI: 10.1039/c7cp07623a] [Citation(s) in RCA: 217] [Impact Index Per Article: 36.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Entropy effects play an important role in drug-target interactions, but the entropic contribution to ligand-binding affinity is often neglected by end-point binding free energy calculation methods, such as MM/GBSA and MM/PBSA, due to the expensive computational cost of normal mode analysis (NMA). Here, we systematically investigated entropy effects on the prediction power of MM/GBSA and MM/PBSA using >1500 protein-ligand systems and six representative AMBER force fields. Two computationally efficient methods, including NMA based on truncated structures and the interaction entropy approach, were used to estimate the entropic contributions to ligand-target binding free energies. In terms of the overall accuracy, we found that, for the minimized structures, in most cases the inclusion of the conformational entropies predicted by truncated NMA (enthalpynmode_min_9Å) compromises the overall accuracy of MM/GBSA and MM/PBSA compared with the enthalpies calculated based on the minimized structures (enthalpymin). However, for the MD trajectories, the binding free energies can be improved by the inclusion of the conformation entropies predicted by either truncated-NMA for a relatively high dielectric constant (εin = 4) or the interaction entropy method for εin = 1-4. In terms of reproducing the absolute binding free energies, the binding free energies estimated by including the truncated-NMA entropies based on the MD trajectories (ΔGnmode_md_9Å) give the lowest average absolute deviations against the experimental data among all the tested strategies for both MM/GBSA and MM/PBSA. Although the inclusion of the truncated NMA based on the MD trajectories (ΔGnmode_md_9Å) for a relatively high dielectric constant gave the overall best result and the lowest average absolute deviations against the experimental data (for the ff03 force field), it needs too much computational time. Alternatively, considering that the interaction entropy method does not incur any additional computational cost and can give comparable (at high dielectric constant, εin = 4) or even better (at low dielectric constant, εin = 1-2) results than the truncated-NMA entropy (ΔGnmode_md_9Å), the interaction entropy approach is recommended to estimate the entropic component for MM/GBSA and MM/PBSA based on MD trajectories, especially for a diverse dataset. Furthermore, we compared the predictions of MM/GBSA with six different AMBER force fields. The results show that the ff03 force field (ff03 for proteins and gaff with AM1-BCC charges for ligands) performs the best, but the predictions given by the tested force fields are comparable, implying that the MM/GBSA predictions are not very sensitive to force fields.
Collapse
Affiliation(s)
- Huiyong Sun
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | | | | | | | | | | | | | | | | |
Collapse
|
30
|
Fox JM, Zhao M, Fink MJ, Kang K, Whitesides GM. The Molecular Origin of Enthalpy/Entropy Compensation in Biomolecular Recognition. Annu Rev Biophys 2018; 47:223-250. [DOI: 10.1146/annurev-biophys-070816-033743] [Citation(s) in RCA: 91] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Biomolecular recognition can be stubborn; changes in the structures of associating molecules, or the environments in which they associate, often yield compensating changes in enthalpies and entropies of binding and no net change in affinities. This phenomenon—termed enthalpy/entropy (H/S) compensation—hinders efforts in biomolecular design, and its incidence—often a surprise to experimentalists—makes interactions between biomolecules difficult to predict. Although characterizing H/S compensation requires experimental care, it is unquestionably a real phenomenon that has, from an engineering perspective, useful physical origins. Studying H/S compensation can help illuminate the still-murky roles of water and dynamics in biomolecular recognition and self-assembly. This review summarizes known sources of H/ S compensation (real and perceived) and lays out a conceptual framework for understanding and dissecting—and, perhaps, avoiding or exploiting—this phenomenon in biophysical systems.
Collapse
Affiliation(s)
- Jerome M. Fox
- Department of Chemical and Biological Engineering, University of Colorado, Boulder, Colorado 80309, USA
| | - Mengxia Zhao
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA;, ,
| | - Michael J. Fink
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA;, ,
| | - Kyungtae Kang
- Department of Applied Chemistry, Kyung Hee University, Yongin, Gyeonggi 17104, Republic of Korea
| | - George M. Whitesides
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA;, ,
- Wyss Institute for Biologically Inspired Engineering, Harvard University, Cambridge, Massachusetts 02138, USA
- The Kavli Institute for Bionano Science and Technology, Harvard University, Cambridge, Massachusetts 02138, USA
| |
Collapse
|
31
|
Leidner F, Kurt Yilmaz N, Paulsen J, Muller YA, Schiffer CA. Hydration Structure and Dynamics of Inhibitor-Bound HIV-1 Protease. J Chem Theory Comput 2018; 14:2784-2796. [PMID: 29570286 DOI: 10.1021/acs.jctc.8b00097] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Water is essential in many biological processes, and the hydration structure plays a critical role in facilitating protein folding, dynamics, and ligand binding. A variety of biophysical spectroscopic techniques have been used to probe the water solvating proteins, often complemented with molecular dynamics (MD) simulations to resolve the spatial and dynamic features of the hydration shell, but comparing relative water structure is challenging. In this study 1 μs MD simulations were performed to identify and characterize hydration sites around HIV-1 protease bound to an inhibitor, darunavir (DRV). The water density, hydration site occupancy, extent and anisotropy of fluctuations, coordinated water molecules, and hydrogen bonds were characterized and compared to the properties of bulk water. The water density of the principal hydration shell was found to be higher than bulk, dependent on the topology and physiochemical identity of the biomolecular surface. The dynamics of water molecules occupying principal hydration sites was highly dependent on the number of water-water interactions and inversely correlated with hydrogen bonds to the protein-inhibitor complex. While many waters were conserved following the symmetry of homodimeric HIV protease, the asymmetry induced by DRV resulted in asymmetric lower-occupancy hydration sites at the concave surface of the active site. Key interactions between water molecules and the protease, that stabilize the protein in the inhibited form, were altered in a drug resistant variant of the protease indicating that modulation of solvent-solute interactions might play a key role in conveying drug resistance. Our analysis provides insights into the interplay between an enzyme inhibitor complex and the hydration shell and has implications in elucidating water structure in a variety of biological processes and applications including ligand binding, inhibitor design, and resistance.
Collapse
Affiliation(s)
- Florian Leidner
- Department of Biochemistry and Molecular Pharmacology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| | - Nese Kurt Yilmaz
- Department of Biochemistry and Molecular Pharmacology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| | - Janet Paulsen
- Department of Biochemistry and Molecular Pharmacology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| | - Yves A Muller
- Division of Biotechnology , Friedrich-Alexander-Universität Erlangen-Nürnberg , Erlangen 91052 , Germany
| | - Celia A Schiffer
- Department of Biochemistry and Molecular Pharmacology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| |
Collapse
|
32
|
Triki D, Billot T, Visseaux B, Descamps D, Flatters D, Camproux AC, Regad L. Exploration of the effect of sequence variations located inside the binding pocket of HIV-1 and HIV-2 proteases. Sci Rep 2018; 8:5789. [PMID: 29636521 PMCID: PMC5893546 DOI: 10.1038/s41598-018-24124-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2018] [Accepted: 03/26/2018] [Indexed: 01/01/2023] Open
Abstract
HIV-2 protease (PR2) is naturally resistant to most FDA (Food and Drug Administration)-approved HIV-1 protease inhibitors (PIs), a major antiretroviral class. In this study, we compared the PR1 and PR2 binding pockets extracted from structures complexed with 12 ligands. The comparison of PR1 and PR2 pocket properties showed that bound PR2 pockets were more hydrophobic with more oxygen atoms and fewer nitrogen atoms than PR1 pockets. The structural comparison of PR1 and PR2 pockets highlighted structural changes induced by their sequence variations and that were consistent with these property changes. Specifically, substitutions at residues 31, 46, and 82 induced structural changes in their main-chain atoms that could affect PI binding in PR2. In addition, the modelling of PR1 mutant structures containing V32I and L76M substitutions revealed a cooperative mechanism leading to structural deformation of flap-residue 45 that could modify PR2 flexibility. Our results suggest that substitutions in the PR1 and PR2 pockets can modify PI binding and flap flexibility, which could underlie PR2 resistance against PIs. These results provide new insights concerning the structural changes induced by PR1 and PR2 pocket variation changes, improving the understanding of the atomic mechanism of PR2 resistance to PIs.
Collapse
Affiliation(s)
- Dhoha Triki
- Molécules thérapeutiques in silico (MTi), INSERM UMR-, S973, Paris, France.,Université Paris Diderot, Sorbonne Paris Cité, Paris, France
| | - Telli Billot
- Molécules thérapeutiques in silico (MTi), INSERM UMR-, S973, Paris, France.,Université Paris Diderot, Sorbonne Paris Cité, Paris, France
| | - Benoit Visseaux
- IAME, UMR 1137, INSERM, Laboratoire de Virologie, Hôpital Bichât, AP-HP, Paris, France.,Université Paris Diderot, Sorbonne Paris Cité, Paris, France
| | - Diane Descamps
- IAME, UMR 1137, INSERM, Laboratoire de Virologie, Hôpital Bichât, AP-HP, Paris, France.,Université Paris Diderot, Sorbonne Paris Cité, Paris, France
| | - Delphine Flatters
- Molécules thérapeutiques in silico (MTi), INSERM UMR-, S973, Paris, France.,Université Paris Diderot, Sorbonne Paris Cité, Paris, France
| | - Anne-Claude Camproux
- Molécules thérapeutiques in silico (MTi), INSERM UMR-, S973, Paris, France.,Université Paris Diderot, Sorbonne Paris Cité, Paris, France
| | - Leslie Regad
- Molécules thérapeutiques in silico (MTi), INSERM UMR-, S973, Paris, France. .,Université Paris Diderot, Sorbonne Paris Cité, Paris, France.
| |
Collapse
|
33
|
Khan SN, Persons JD, Paulsen JL, Guerrero M, Schiffer CA, Kurt-Yilmaz N, Ishima R. Probing Structural Changes among Analogous Inhibitor-Bound Forms of HIV-1 Protease and a Drug-Resistant Mutant in Solution by Nuclear Magnetic Resonance. Biochemistry 2018; 57:1652-1662. [PMID: 29457713 PMCID: PMC5850901 DOI: 10.1021/acs.biochem.7b01238] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the era of state-of-the-art inhibitor design and high-resolution structural studies, detection of significant but small protein structural differences in the inhibitor-bound forms is critical to further developing the inhibitor. Here, we probed differences in HIV-1 protease (PR) conformation among darunavir and four analogous inhibitor-bound forms and compared them with a drug-resistant mutant using nuclear magnetic resonance chemical shifts. Changes in amide chemical shifts of wild-type (WT) PR among these inhibitor-bound forms, ΔCSP, were subtle but detectable and extended >10 Å from the inhibitor-binding site, asymmetrically between the two subunits of PR. Molecular dynamics simulations revealed differential local hydrogen bonding as the molecular basis of this remote asymmetric change. Inhibitor-bound forms of the drug-resistant mutant also showed a similar long-range ΔCSP pattern. Differences in ΔCSP values of the WT and the mutant (ΔΔCSPs) were observed at the inhibitor-binding site and in the surrounding region. Comparing chemical shift changes among highly analogous inhibitors and ΔΔCSPs effectively eliminated local environmental effects stemming from different chemical groups and enabled exploitation of these sensitive parameters to detect subtle protein conformational changes and to elucidate asymmetric and remote conformational effects upon inhibitor interaction.
Collapse
Affiliation(s)
- Shahid N Khan
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - John D Persons
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - Janet L. Paulsen
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, USA
| | - Michel Guerrero
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| | - Celia A. Schiffer
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, USA
| | - Nese Kurt-Yilmaz
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA, USA
| | - Rieko Ishima
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, USA
| |
Collapse
|
34
|
Sanusi ZK, Govender T, Maguire GEM, Maseko SB, Lin J, Kruger HG, Honarparvar B. An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant. J Comput Aided Mol Des 2018; 32:459-471. [PMID: 29397520 DOI: 10.1007/s10822-018-0099-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2017] [Accepted: 01/16/2018] [Indexed: 01/12/2023]
Abstract
The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.
Collapse
Affiliation(s)
- Zainab K Sanusi
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Thavendran Govender
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Glenn E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.,School of Chemistry and Physics, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Sibusiso B Maseko
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Johnson Lin
- School of Life Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Hendrik G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| | - Bahareh Honarparvar
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| |
Collapse
|
35
|
Sanusi ZK, Govender T, Maguire GEM, Maseko SB, Lin J, Kruger HG, Honarparvar B. Investigation of the binding free energies of FDA approved drugs against subtype B and C-SA HIV PR: ONIOM approach. J Mol Graph Model 2017; 76:77-85. [PMID: 28711760 DOI: 10.1016/j.jmgm.2017.06.026] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2017] [Revised: 06/08/2017] [Accepted: 06/09/2017] [Indexed: 01/15/2023]
Abstract
Human immune virus subtype C is the most widely spread HIV subtype in Sub-Sahara Africa and South Africa. A profound structural insight on finding potential lead compounds is therefore necessary for drug discovery. The focus of this study is to rationalize the nine Food and Drugs Administration (FDA) HIV antiviral drugs complexed to subtype B and C-SA PR using ONIOM approach. To achieve this, an integrated two-layered ONIOM model was used to optimize the geometrics of the FDA approved HIV-1 PR inhibitors for subtype B. In our hybrid ONIOM model, the HIV-1 PR inhibitors as well as the ASP 25/25' catalytic active residues were treated at high level quantum mechanics (QM) theory using B3LYP/6-31G(d), and the remaining HIV PR residues were considered using the AMBER force field. The experimental binding energies of the PR inhibitors were compared to the ONIOM calculated results. The theoretical binding free energies (?Gbind) for subtype B follow a similar trend to the experimental results, with one exemption. The computational model was less suitable for C-SA PR. Analysis of the results provided valuable information about the shortcomings of this approach. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide much improved binding energies for complex enzyme drug interactions.
Collapse
Affiliation(s)
- Z K Sanusi
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - T Govender
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - G E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa; School of Chemistry and Physics, University of KwaZulu-Natal, 4001 Durban, South Africa
| | - S B Maseko
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - J Lin
- School of Life Sciences, University of KwaZulu-Natal, Durban 4001, South Africa
| | - H G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa.
| | - B Honarparvar
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban 4001, South Africa.
| |
Collapse
|
36
|
Gopal SM, Klumpers F, Herrmann C, Schäfer LV. Solvent effects on ligand binding to a serine protease. Phys Chem Chem Phys 2017; 19:10753-10766. [DOI: 10.1039/c6cp07899k] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
ITC experiments and MD simulations reveal the mechanism behind enthalpy/entropy compensation upon trypsin-benzamidine binding at different solvation conditions.
Collapse
Affiliation(s)
- Srinivasa M. Gopal
- Center for Theoretical Chemistry
- Faculty of Chemistry and Biochemistry
- Ruhr-University Bochum
- D-44780 Bochum
- Germany
| | - Fabian Klumpers
- Physical Chemistry I
- Faculty of Chemistry and Biochemistry
- Ruhr-University Bochum
- D-44780 Bochum
- Germany
| | - Christian Herrmann
- Physical Chemistry I
- Faculty of Chemistry and Biochemistry
- Ruhr-University Bochum
- D-44780 Bochum
- Germany
| | - Lars V. Schäfer
- Center for Theoretical Chemistry
- Faculty of Chemistry and Biochemistry
- Ruhr-University Bochum
- D-44780 Bochum
- Germany
| |
Collapse
|
37
|
Yu Y, Wang J, Chen Z, Wang G, Shao Q, Shi J, Zhu W. Structural insights into HIV-1 protease flap opening processes and key intermediates. RSC Adv 2017. [DOI: 10.1039/c7ra09691g] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
The study provided an integrated view of the transition pathway of the flap opening of HIV-1 protease using MD simulation.
Collapse
Affiliation(s)
- Yuqi Yu
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| | - Jinan Wang
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| | - Zhaoqiang Chen
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| | - Guimin Wang
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| | - Qiang Shao
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| | - Jiye Shi
- UCB Biopharma SPRL
- Chemin du Foriest
- Belgium
| | - Weiliang Zhu
- Drug Discovery and Design Center
- CAS Key Laboratory of Receptor Research
- Shanghai Institute of Materia Medica
- Chinese Academy of Sciences
- Shanghai
| |
Collapse
|
38
|
Liu Z, Huang X, Hu L, Pham L, Poole KM, Tang Y, Mahon BP, Tang W, Li K, Goldfarb NE, Dunn BM, McKenna R, Fanucci GE. Effects of Hinge-region Natural Polymorphisms on Human Immunodeficiency Virus-Type 1 Protease Structure, Dynamics, and Drug Pressure Evolution. J Biol Chem 2016; 291:22741-22756. [PMID: 27576689 DOI: 10.1074/jbc.m116.747568] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2016] [Revised: 08/30/2016] [Indexed: 11/06/2022] Open
Abstract
Multidrug resistance to current Food and Drug Administration-approved HIV-1 protease (PR) inhibitors drives the need to understand the fundamental mechanisms of how drug pressure-selected mutations, which are oftentimes natural polymorphisms, elicit their effect on enzyme function and resistance. Here, the impacts of the hinge-region natural polymorphism at residue 35, glutamate to aspartate (E35D), alone and in conjunction with residue 57, arginine to lysine (R57K), are characterized with the goal of understanding how altered salt bridge interactions between the hinge and flap regions are associated with changes in structure, motional dynamics, conformational sampling, kinetic parameters, and inhibitor affinity. The combined results reveal that the single E35D substitution leads to diminished salt bridge interactions between residues 35 and 57 and gives rise to the stabilization of open-like conformational states with overall increased backbone dynamics. In HIV-1 PR constructs where sites 35 and 57 are both mutated (e.g. E35D and R57K), x-ray structures reveal an altered network of interactions that replace the salt bridge thus stabilizing the structural integrity between the flap and hinge regions. Despite the altered conformational sampling and dynamics when the salt bridge is disrupted, enzyme kinetic parameters and inhibition constants are similar to those obtained for subtype B PR. Results demonstrate that these hinge-region natural polymorphisms, which may arise as drug pressure secondary mutations, alter protein dynamics and the conformational landscape, which are important thermodynamic parameters to consider for development of inhibitors that target for non-subtype B PR.
Collapse
Affiliation(s)
- Zhanglong Liu
- From the Department of Chemistry, University of Florida, Gainesville, Florida 32611 and
| | - Xi Huang
- From the Department of Chemistry, University of Florida, Gainesville, Florida 32611 and
| | - Lingna Hu
- From the Department of Chemistry, University of Florida, Gainesville, Florida 32611 and
| | - Linh Pham
- From the Department of Chemistry, University of Florida, Gainesville, Florida 32611 and
| | - Katye M Poole
- the Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, Florida 32610
| | - Yan Tang
- the Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, Florida 32610
| | - Brian P Mahon
- the Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, Florida 32610
| | - Wenxing Tang
- the Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, Florida 32610
| | - Kunhua Li
- From the Department of Chemistry, University of Florida, Gainesville, Florida 32611 and
| | - Nathan E Goldfarb
- the Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, Florida 32610
| | - Ben M Dunn
- the Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, Florida 32610
| | - Robert McKenna
- the Department of Biochemistry and Molecular Biology, University of Florida, Gainesville, Florida 32610
| | - Gail E Fanucci
- From the Department of Chemistry, University of Florida, Gainesville, Florida 32611 and
| |
Collapse
|
39
|
Ahmad M, Helms V, Lengauer T, Kalinina OV. How Molecular Conformational Changes Affect Changes in Free Energy. J Chem Theory Comput 2016; 11:2945-57. [PMID: 26575732 DOI: 10.1021/acs.jctc.5b00235] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A simple quantitative relationship between the molecular conformational changes and the corresponding changes in the free energy is presented. The change in free energy is the sum of that part of the enthalpic change that is due to the externally applied work (perturbation) and of that part of the entropic change, termed dissipative entropy, that is related to the conformational changes. The dissipative entropy is equivalent to the relative entropy, a concept from information theory, between the distributions of the conformations in the initial and the final states. The remaining change in entropy (nondissipative) cancels exactly with the remaining enthalpic change. The calculation of the dissipative entropy is demonstrated to pose the main difficulty in free energy computation. The straightforward decomposition of the dissipative entropy into contributions from different parts of the system promises to improve the understanding of the role of conformational changes in biochemical reactions.
Collapse
Affiliation(s)
- Mazen Ahmad
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics , Campus E1 4, 66123 Saarbrücken, Germany
| | - Volkhard Helms
- Center for Bioinformatics, Saarland University , 66123 Saarbrücken, Germany
| | - Thomas Lengauer
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics , Campus E1 4, 66123 Saarbrücken, Germany
| | - Olga V Kalinina
- Department for Computational Biology and Applied Algorithmics, Max Planck Institute for Informatics , Campus E1 4, 66123 Saarbrücken, Germany
| |
Collapse
|
40
|
Park JH, Sayer JM, Aniana A, Yu X, Weber IT, Harrison RW, Louis JM. Binding of Clinical Inhibitors to a Model Precursor of a Rationally Selected Multidrug Resistant HIV-1 Protease Is Significantly Weaker Than That to the Released Mature Enzyme. Biochemistry 2016; 55:2390-400. [PMID: 27039930 DOI: 10.1021/acs.biochem.6b00012] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have systematically validated the activity and inhibition of a HIV-1 protease (PR) variant bearing 17 mutations (PR(S17)), selected to represent high resistance by machine learning on genotype-phenotype data. Three of five mutations in PR(S17) correlating with major drug resistance, M46L, G48V, and V82S, and five of 11 natural variations differ from the mutations in two clinically derived extreme mutants, PR20 and PR22 bearing 19 and 22 mutations, respectively. PR(S17), which forms a stable dimer (<10 nM), is ∼10- and 2-fold less efficient in processing the Gag polyprotein than the wild type and PR20, respectively, but maintains the same cleavage order. Isolation of a model precursor of PR(S17) flanked by the 56-amino acid transframe region (TFP-p6pol) at its N-terminus, which is impossible upon expression of an analogous PR20 precursor, allowed systematic comparison of inhibition of TFP-p6pol-PR(S17) and mature PR(S17). Resistance of PR(S17) to eight protease inhibitors (PIs) relative to PR (Ki) increases by 1.5-5 orders of magnitude from 0.01 to 8.4 μM. Amprenavir, darunavir, atazanavir, and lopinavir, the most effective of the eight PIs, inhibit precursor autoprocessing at the p6pol/PR site with IC50 values ranging from ∼7.5 to 60 μM. Thus, this process, crucial for stable dimer formation, shows inhibition ∼200-800-fold weaker than that of the mature PR(S17). TFP/p6pol cleavage, which occurs faster, is inhibited even more weakly by all PIs except darunavir (IC50 = 15 μM); amprenavir shows a 2-fold increase in IC50 (∼15 μM), and atazanavir and lopinavir show increased IC50 values of >42 and >70 μM, respectively.
Collapse
Affiliation(s)
- Joon H Park
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Department of Health and Human Services , Bethesda, Maryland 20892, United States
| | - Jane M Sayer
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Department of Health and Human Services , Bethesda, Maryland 20892, United States
| | - Annie Aniana
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Department of Health and Human Services , Bethesda, Maryland 20892, United States
| | | | | | | | - John M Louis
- Laboratory of Chemical Physics, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Department of Health and Human Services , Bethesda, Maryland 20892, United States
| |
Collapse
|
41
|
Abstract
How to design a ligand to bind multiple targets, rather than to a single target, is the focus of this review. Rational polypharmacology draws on knowledge that is both broad ranging and hierarchical. Computer-aided multitarget ligand design methods are described according to their nested knowledge level. Ligand-only and then receptor-ligand strategies are first described; followed by the metabolic network viewpoint. Subsequently strategies that view infectious diseases as multigenomic targets are discussed, and finally the disease level interpretation of medicinal therapy is considered. As yet there is no consensus on how best to proceed in designing a multitarget ligand. The current methodologies are bought together in an attempt to give a practical overview of how polypharmacology design might be best initiated.
Collapse
|
42
|
Anighoro A, Bajorath J. Three-Dimensional Similarity in Molecular Docking: Prioritizing Ligand Poses on the Basis of Experimental Binding Modes. J Chem Inf Model 2016; 56:580-7. [DOI: 10.1021/acs.jcim.5b00745] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Andrew Anighoro
- Department of Life Science
Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal
Chemistry, Rheinische Friedrich-Wilhelms-Universität, Dahlmannstrasse 2, D-53113 Bonn, Germany
| | - Jürgen Bajorath
- Department of Life Science
Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal
Chemistry, Rheinische Friedrich-Wilhelms-Universität, Dahlmannstrasse 2, D-53113 Bonn, Germany
| |
Collapse
|
43
|
Lockhat HA, Silva JRA, Alves CN, Govender T, Lameira J, Maguire GEM, Sayed Y, Kruger HG. Binding Free Energy Calculations of Nine FDA-approved Protease Inhibitors Against HIV-1 Subtype C I36T↑T Containing 100 Amino Acids Per Monomer. Chem Biol Drug Des 2016; 87:487-98. [PMID: 26613568 DOI: 10.1111/cbdd.12690] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2015] [Revised: 09/28/2015] [Accepted: 10/22/2015] [Indexed: 12/19/2022]
Abstract
In this work, have investigated the binding affinities of nine FDA-approved protease inhibitor drugs against a new HIV-1 subtype C mutated protease, I36T↑T. Without an X-ray crystal structure, homology modelling was used to generate a three-dimensional model of the protease. This and the inhibitor models were employed to generate the inhibitor/I36T↑T complexes, with the relative positions of the inhibitors being superimposed and aligned using the X-ray crystal structures of the inhibitors/HIV-1 subtype B complexes as a reference. Molecular dynamics simulations were carried out on the complexes to calculate the average binding free energies for each inhibitor using the molecular mechanics generalized Born surface area (MM-GBSA) method. When compared to the binding free energies of the HIV-1 subtype B and subtype C proteases (calculated previously by our group using the same method), it was clear that the I36T↑T proteases mutations and insertion had a significant negative effect on the binding energies of the non-pepditic inhibitors nelfinavir, darunavir and tipranavir. On the other hand, ritonavir, amprenavir and indinavir show improved calculated binding energies in comparison with the corresponding data for wild-type C-SA protease. The computational model used in this study can be used to investigate new mutations of the HIV protease and help in establishing effective HIV drug regimes and may also aid in future protease drug design.
Collapse
Affiliation(s)
- Husain A Lockhat
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - José R A Silva
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, Belém, PA, 66075-110, Brazil
| | - Cláudio N Alves
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, Belém, PA, 66075-110, Brazil
| | - Thavendran Govender
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Jerônimo Lameira
- Laboratório de Planejamento e Desenvolvimento de Fármacos, Instituto de Ciências Exatas e Naturais, Universidade Federal do Pará, CP 11101, Belém, PA, 66075-110, Brazil
| | - Glenn E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.,School of Chemistry and Physics, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Yasien Sayed
- Protein Structure-Function Research Unit, School of Molecular and Cell Biology, University of the Witwatersrand, Wits, 2050, South Africa
| | - Hendrik G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| |
Collapse
|
44
|
Chen J. Drug resistance mechanisms of three mutations V32I, I47V and V82I in HIV-1 protease toward inhibitors probed by molecular dynamics simulations and binding free energy predictions. RSC Adv 2016. [DOI: 10.1039/c6ra09201b] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
Molecular dynamics simulation and binding free energy calculations were used to probe drug resistance of HIV-1 protease mutations toward inhibitors.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- China
| |
Collapse
|
45
|
Four Amino Acid Changes in HIV-2 Protease Confer Class-Wide Sensitivity to Protease Inhibitors. J Virol 2015; 90:1062-9. [PMID: 26559830 DOI: 10.1128/jvi.01772-15] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2015] [Accepted: 11/02/2015] [Indexed: 12/30/2022] Open
Abstract
UNLABELLED Protease is essential for retroviral replication, and protease inhibitors (PI) are important for treating HIV infection. HIV-2 exhibits intrinsic resistance to most FDA-approved HIV-1 PI, retaining clinically useful susceptibility only to lopinavir, darunavir, and saquinavir. The mechanisms for this resistance are unclear; although HIV-1 and HIV-2 proteases share just 38 to 49% sequence identity, all critical structural features of proteases are conserved. Structural studies have implicated four amino acids in the ligand-binding pocket (positions 32, 47, 76, and 82). We constructed HIV-2ROD9 molecular clones encoding the corresponding wild-type HIV-1 amino acids (I32V, V47I, M76L, and I82V) either individually or together (clone PRΔ4) and compared the phenotypic sensitivities (50% effective concentration [EC50]) of mutant and wild-type viruses to nine FDA-approved PI. Single amino acid replacements I32V, V47I, and M76L increased the susceptibility of HIV-2 to multiple PI, but no single change conferred class-wide sensitivity. In contrast, clone PRΔ4 showed PI susceptibility equivalent to or greater than that of HIV-1 for all PI. We also compared crystallographic structures of wild-type HIV-1 and HIV-2 proteases complexed with amprenavir and darunavir to models of the PRΔ4 enzyme. These models suggest that the amprenavir sensitivity of PRΔ4 is attributable to stabilizing enzyme-inhibitor interactions in the P2 and P2' pockets of the protease dimer. Together, our results show that the combination of four amino acid changes in HIV-2 protease confer a pattern of PI susceptibility comparable to that of HIV-1, providing a structural rationale for intrinsic HIV-2 PI resistance and resolving long-standing questions regarding the determinants of differential PI susceptibility in HIV-1 and HIV-2. IMPORTANCE Proteases are essential for retroviral replication, and HIV-1 and HIV-2 proteases share a great deal of structural similarity. However, only three of nine FDA-approved HIV-1 protease inhibitors (PI) are active against HIV-2. The underlying reasons for intrinsic PI resistance in HIV-2 are not known. We examined the contributions of four amino acids in the ligand-binding pocket of the enzyme that differ between HIV-1 and HIV-2 by constructing HIV-2 clones encoding the corresponding HIV-1 amino acids and testing the PI susceptibilities of the resulting viruses. We found that the HIV-2 clone containing all four changes (PRΔ4) was as susceptible as HIV-1 to all nine PI. We also modeled the PRΔ4 enzyme structure and compared it to existing crystallographic structures of HIV-1 and HIV-2 proteases complexed with amprenavir and darunavir. Our findings demonstrate that four positions in the ligand-binding cleft of protease are the primary cause of HIV-2 PI resistance.
Collapse
|
46
|
Shen Y, Radhakrishnan ML, Tidor B. Molecular mechanisms and design principles for promiscuous inhibitors to avoid drug resistance: lessons learned from HIV-1 protease inhibition. Proteins 2015; 83:351-72. [PMID: 25410041 PMCID: PMC4829108 DOI: 10.1002/prot.24730] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2014] [Revised: 10/14/2014] [Accepted: 11/06/2014] [Indexed: 11/16/2022]
Abstract
Molecular recognition is central to biology and ranges from highly selective to broadly promiscuous. The ability to modulate specificity at will is particularly important for drug development, and discovery of mechanisms contributing to binding specificity is crucial for our basic understanding of biology and for applications in health care. In this study, we used computational molecular design to create a large dataset of diverse small molecules with a range of binding specificities. We then performed structural, energetic, and statistical analysis on the dataset to study molecular mechanisms of achieving specificity goals. The work was done in the context of HIV‐1 protease inhibition and the molecular designs targeted a panel of wild‐type and drug‐resistant mutant HIV‐1 protease structures. The analysis focused on mechanisms for promiscuous binding to bind robustly even to resistance mutants. Broadly binding inhibitors tended to be smaller in size, more flexible in chemical structure, and more hydrophobic in nature compared to highly selective ones. Furthermore, structural and energetic analyses illustrated mechanisms by which flexible inhibitors achieved binding; we found ligand conformational adaptation near mutation sites and structural plasticity in targets through torsional flips of asymmetric functional groups to form alternative, compensatory packing interactions or hydrogen bonds. As no inhibitor bound to all variants, we designed small cocktails of inhibitors to do so and discovered that they often jointly covered the target set through mechanistic complementarity. Furthermore, using structural plasticity observed in experiments, and potentially in simulations, is suggested to be a viable means of designing adaptive inhibitors that are promiscuous binders. Proteins 2015; 83:351–372. © 2014 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Yang Shen
- Department of Biological EngineeringMassachusetts Institute of TechnologyCambridgeMassachusetts02139
- Department of Electrical Engineering and Computer ScienceMassachusetts Institute of TechnologyCambridgeMassachusetts02139
- Computer Science and Artificial Intelligence LaboratoryMassachusetts Institute of TechnologyCambridgeMassachusetts02139
- Present address:
Center for Bioinformatics and Genomic Systems EngineeringDepartment of Electrical and Computer EngineeringTexas A&M UniversityCollege StationTexas77843
| | | | - Bruce Tidor
- Department of Biological EngineeringMassachusetts Institute of TechnologyCambridgeMassachusetts02139
- Department of Electrical Engineering and Computer ScienceMassachusetts Institute of TechnologyCambridgeMassachusetts02139
- Computer Science and Artificial Intelligence LaboratoryMassachusetts Institute of TechnologyCambridgeMassachusetts02139
| |
Collapse
|
47
|
Chen J, Wang X, Zhu T, Zhang Q, Zhang JZH. A Comparative Insight into Amprenavir Resistance of Mutations V32I, G48V, I50V, I54V, and I84V in HIV-1 Protease Based on Thermodynamic Integration and MM-PBSA Methods. J Chem Inf Model 2015; 55:1903-13. [DOI: 10.1021/acs.jcim.5b00173] [Citation(s) in RCA: 98] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Jianzhong Chen
- School
of Science, Shandong Jiaotong University, Jinan 250357 China
| | - Xingyu Wang
- NYU−ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Tong Zhu
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| | - Qinggang Zhang
- Collage
of Physics and Electronic Science, Shandong Normal University, Jinan 250014, China
| | - John Z. H. Zhang
- NYU−ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| |
Collapse
|
48
|
Ekins S, Litterman NK, Lipinski CA, Bunin BA. Thermodynamic Proxies to Compensate for Biases in Drug Discovery Methods. Pharm Res 2015; 33:194-205. [DOI: 10.1007/s11095-015-1779-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2015] [Accepted: 08/13/2015] [Indexed: 11/24/2022]
|
49
|
Laco GS. HIV-1 protease substrate-groove: Role in substrate recognition and inhibitor resistance. Biochimie 2015; 118:90-103. [PMID: 26300060 DOI: 10.1016/j.biochi.2015.08.009] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2015] [Accepted: 08/18/2015] [Indexed: 11/17/2022]
Abstract
A key target in the treatment of HIV-1/AIDS has been the viral protease. Here we first studied in silico the evolution of protease resistance. Primary active site resistance mutations were found to weaken interactions between protease and both inhibitor and substrate P4-P4' residues. We next studied the effects of secondary resistance mutations, often distant from the active site, on protease binding to inhibitors and substrates. Those secondary mutations contributed to the rise of multi-drug resistance while also enhancing viral replicative capacity. Here many secondary resistance mutations were found in the HIV-1 protease substrate-grooves, one on each face of the symmetrical protease dimer. The protease active site binds substrate P4-P4' residues, while the substrate-groove allows the protease to bind residues P12-P5/P5'-P12', for a total of twenty-four residues. The substrate-groove secondary resistance mutations were found to compensate for the loss of interactions between the inhibitor resistant protease active site and substrate P4-P4' residues, due to primary resistance mutations, by increasing interactions with substrate P12-P5/P5'-P12' residues. In vitro experiments demonstrated that a multi-drug resistant protease with substrate-groove resistance mutations was slower than wild-type protease in cleaving a peptide substrate, which did not allow for substrate-groove interactions, while it had similar activity as wild-type protease when using a Gag polyprotein in which cleavage-site P12-P5/P5'-P12' residues could be bound by the protease substrate-grooves. When the Gag MA/CA cleavage site P12-P5/P5'-P12' residues were mutated the multi-drug resistant protease cleaved the mutant Gag significantly slower, indicating the importance of the protease S-grooves in binding to substrate.
Collapse
Affiliation(s)
- Gary S Laco
- Laboratory of Computational and Molecular Biochemistry, The Roskamp Institute, Sarasota, FL, USA.
| |
Collapse
|
50
|
Effects of drug-resistant mutations on the dynamic properties of HIV-1 protease and inhibition by Amprenavir and Darunavir. Sci Rep 2015; 5:10517. [PMID: 26012849 PMCID: PMC4444956 DOI: 10.1038/srep10517] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2015] [Accepted: 04/16/2015] [Indexed: 12/20/2022] Open
Abstract
Molecular dynamics simulations are performed to investigate the dynamic properties of wild-type HIV-1 protease and its two multi-drug-resistant variants (Flap + (L10I/G48V/I54V/V82A) and Act (V82T/I84V)) as well as their binding with APV and DRV inhibitors. The hydrophobic interactions between flap and 80 s (80’s) loop residues (mainly I50-I84’ and I50’-I84) play an important role in maintaining the closed conformation of HIV-1 protease. The double mutation in Act variant weakens the hydrophobic interactions, leading to the transition from closed to semi-open conformation of apo Act. APV or DRV binds with HIV-1 protease via both hydrophobic and hydrogen bonding interactions. The hydrophobic interactions from the inhibitor is aimed to the residues of I50 (I50’), I84 (I84’), and V82 (V82’) which create hydrophobic core clusters to further stabilize the closed conformation of flaps, and the hydrogen bonding interactions are mainly focused with the active site of HIV-1 protease. The combined change in the two kinds of protease-inhibitor interactions is correlated with the observed resistance mutations. The present study sheds light on the microscopic mechanism underlying the mutation effects on the dynamics of HIV-1 protease and the inhibition by APV and DRV, providing useful information to the design of more potent and effective HIV-1 protease inhibitors.
Collapse
|