1
|
Wang J, Yang W, Zhao L, Wei B, Chen J. Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning. Molecules 2024; 29:1857. [PMID: 38675678 PMCID: PMC11054187 DOI: 10.3390/molecules29081857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2024] [Revised: 04/16/2024] [Accepted: 04/16/2024] [Indexed: 04/28/2024] Open
Abstract
Bromodomain 4 and 9 (BRD4 and BRD9) have been regarded as important targets of drug designs in regard to the treatment of multiple diseases. In our current study, molecular dynamics (MD) simulations, deep learning (DL) and binding free energy calculations are integrated to probe the binding modes of three inhibitors (H1B, JQ1 and TVU) to BRD4 and BRD9. The MD trajectory-based DL successfully identify significant functional function domains, such as BC-loop and ZA-loop. The information from the post-processing analysis of MD simulations indicates that inhibitor binding highly influences the structural flexibility and dynamic behavior of BRD4 and BRD9. The results of the MM-GBSA calculations not only suggest that the binding ability of H1B, JQ1 and TVU to BRD9 are stronger than to BRD4, but they also verify that van der Walls interactions are the primary forces responsible for inhibitor binding. The hot spots of BRD4 and BRD9 revealed by residue-based free energy estimation provide target sites of drug design in regard to BRD4 and BRD9. This work is anticipated to provide useful theoretical aids for the development of selective inhibitors over BRD family members.
Collapse
Affiliation(s)
- Jian Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Wanchun Yang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Lu Zhao
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
2
|
Chen J, Wang W, Sun H, He W. Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters. Mini Rev Med Chem 2024; 24:1323-1333. [PMID: 38265367 DOI: 10.2174/0113895575252165231122095555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Revised: 09/05/2023] [Accepted: 10/16/2023] [Indexed: 01/25/2024]
Abstract
Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| | - Weikai He
- School of Science, Shandong Jiaotong University, Jinan-250357, China
| |
Collapse
|
3
|
Sun J, Liang S, Liu X, Zhang S, Li M, Zhang Q, Chen J. Insights into the selectivity of a brain-penetrant CDK4/6 vs CDK1/2 inhibitor for glioblastoma used in multiple replica molecular dynamics simulations. J Biomol Struct Dyn 2023:1-20. [PMID: 38112295 DOI: 10.1080/07391102.2023.2294175] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2023] [Accepted: 11/23/2023] [Indexed: 12/21/2023]
Abstract
Cyclin dependent kinases (CDKs) play an important role in cell cycle regulation and their dysfunction is associated with many cancers. That is why CDKs have been attractive targets for the treatment of cancer. Glioblastoma is a cancer caused by the aberrant expression of CDK4/6, so exploring the mechanism of the selection of CDK4/6 toward inhibitors relative to the other family members CDK1/2 is essential. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), molecular mechanics Poisson-Boltzmann/Generalized Born surface area (MM-PB/GBSA) and other methods were integrated to decipher the selectively binding mechanism of the inhibitor N1J to CDK4/6 and CDK1/2. Molecular electrostatic potential (MESP) analysis provides an explanation for the N1J selectivity. Residue-based free energy decomposition reveals that most of the hot residues are located at the same location of CDKs proteins, but the different types of residues in different proteins cause changes in binding energy, which is considered as a potential developmental direction to improve the selectivity of inhibitors to CDK4/6. These results provide insights into the source of inhibitor and CDK4/6 selectivity for the future development of more selective inhibitors.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Jiahao Sun
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shanshan Liang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
4
|
Alharthy RD, Fatima G, Yousaf N, Iqbal MS, Sattar S, Alanzi AR, Ali I, Muddassar M. Binding selectivity analysis of AURKs inhibitors through molecular dynamics simulation studies. PLoS One 2023; 18:e0295741. [PMID: 38113210 PMCID: PMC10729953 DOI: 10.1371/journal.pone.0295741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2023] [Accepted: 11/29/2023] [Indexed: 12/21/2023] Open
Abstract
Aurora kinases (AURKs) have been identified as promising biological targets for the treatment of cancer. In this study, molecular dynamics simulations were employed to investigate the binding selectivity of three inhibitors (HPM, MPY, and VX6) towards AURKA and AURKB by predicting their binding free energies. The results show that the inhibitors HPM, MPY, and VX6 have more favorable interactions with AURKB as compared to AURKA. The binding energy decomposition analysis revealed that four common residue pairs (L139, L83), (V147, V91), (L210, L154), and (L263, L207) showed significant binding energies with HPM, MPY, and VX6, hence responsible for the binding selectivity of AURKA and AURKB to the inhibitors. The MD trajectory analysis also revealed that the inhibitors affect the dynamic flexibility of protein structure, which is also responsible for the partial selectivity of HPM, MPY, and VX6 towards AURKA and AURKB. As expected, this study provides useful insights for the design of potential inhibitors with high selectivity for AURKA and AURKB.
Collapse
Affiliation(s)
- Rima D. Alharthy
- Department of Chemistry, Science and Arts College, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Ghulam Fatima
- Department of Biosciences, COMSATS University Islamabad, Islamabad, Pakistan
| | - Numan Yousaf
- Department of Biosciences, COMSATS University Islamabad, Islamabad, Pakistan
| | | | - Sadia Sattar
- Department of Biosciences, COMSATS University Islamabad, Islamabad, Pakistan
| | - Abdullah R. Alanzi
- Department of Pharmacognosy, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Ijaz Ali
- Centre for Applied Mathematics and Bioinformatics (CAMB), Gulf University for Science and Technology, Hawally, Kuwait
| | - Muhammad Muddassar
- Department of Biosciences, COMSATS University Islamabad, Islamabad, Pakistan
| |
Collapse
|
5
|
Zaib S, Rana N, Ali HS, Ur Rehman M, Awwad NS, Ibrahium HA, Khan I. Identification of potential inhibitors targeting yellow fever virus helicase through ligand and structure-based computational studies. J Biomol Struct Dyn 2023:1-18. [PMID: 38109183 DOI: 10.1080/07391102.2023.2294839] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 12/05/2023] [Indexed: 12/19/2023]
Abstract
Yellow fever is a flavivirus having plus-sensed RNA which encodes a single polyprotein. Host proteases cut this polyprotein into seven nonstructural proteins including a vital NS3 protein. The present study aims to identify the most effective inhibitor against the helicase (NS3) using different advanced ligand and structure-based computational studies. A set of 300 ligands was selected against helicase by chemical structural similarity model, which are similar to S-adenosyl-l-cysteine using infiniSee. This tool screens billions of compounds through a similarity search from in-built chemical spaces (CHEMriya, Galaxi, KnowledgeSpace and REALSpace). The pharmacophore was designed from ligands in the library that showed same features. According to the sequence of ligands, six compounds (29, 87, 99, 116, 148, and 208) were taken for pharmacophore designing against helicase protein. Subsequently, compounds from the library which showed the best pharmacophore shared-features were docked using FlexX functionality of SeeSAR and their optibrium properties were analyzed. Afterward, their ADME was improved by replacing the unfavorable fragments, which resulted in the generation of new compounds. The selected best compounds (301, 302, 303 and 304) were docked using SeeSAR and their pharmacokinetics and toxicological properties were evaluated using SwissADME. The optimal inhibitor for yellow fever helicase was 2-amino-N-(4-(dimethylamino)thiazol-2-yl)-4-methyloxazole-5-carboxamide (302), which exhibits promising potential for drug development.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sumera Zaib
- Department of Basic and Applied Chemistry, Faculty of Science and Technology, University of Central Punjab, Lahore, Pakistan
| | - Nehal Rana
- Department of Basic and Applied Chemistry, Faculty of Science and Technology, University of Central Punjab, Lahore, Pakistan
| | - Hafiz Saqib Ali
- Chemistry Research Laboratory, Department of Chemistry and the INEOS Oxford Institute for Antimicrobial Research, University of Oxford, Oxford, UK
| | - Mujeeb Ur Rehman
- Department of Basic and Applied Chemistry, Faculty of Science and Technology, University of Central Punjab, Lahore, Pakistan
| | - Nasser S Awwad
- Department of Chemistry, King Khalid University, Abha, Saudi Arabia
| | - Hala A Ibrahium
- Biology Department, Faculty of Science, King Khalid University, Abha, Saudi Arabia
| | - Imtiaz Khan
- Manchester Institute of Biotechnology, The University of Manchester, Manchester, UK
| |
Collapse
|
6
|
Sun J, Liu X, Zhang S, Li M, Zhang Q, Chen J. Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study. Phys Chem Chem Phys 2023; 25:28479-28496. [PMID: 37846774 DOI: 10.1039/d3cp03392a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.
Collapse
Affiliation(s)
- Jiahao Sun
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| |
Collapse
|
7
|
Wang Y, Yang F, Yan D, Zeng Y, Wei B, Chen J, He W. Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies. Molecules 2023; 28:4773. [PMID: 37375328 DOI: 10.3390/molecules28124773] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Revised: 06/12/2023] [Accepted: 06/12/2023] [Indexed: 06/29/2023] Open
Abstract
β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.
Collapse
Affiliation(s)
- Yiwen Wang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| | - Fen Yang
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Dongliang Yan
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Yalin Zeng
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
| | - Benzheng Wei
- Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
| | - Jianzhong Chen
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Weikai He
- School of Information Science and Electrical Engineering, Shandong Jiaotong University, Jinan 250357, China
- School of Aeronautics, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
8
|
Nguyen TH, Tam NM, Tuan MV, Zhan P, Vu VV, Quang DT, Ngo ST. Searching for potential inhibitors of SARS-COV-2 main protease using supervised learning and perturbation calculations. Chem Phys 2023; 564:111709. [PMID: 36188488 PMCID: PMC9511900 DOI: 10.1016/j.chemphys.2022.111709] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Revised: 08/11/2022] [Accepted: 09/21/2022] [Indexed: 11/28/2022]
Abstract
Inhibiting the biological activity of SARS-CoV-2 Mpro can prevent viral replication. In this context, a hybrid approach using knowledge- and physics-based methods was proposed to characterize potential inhibitors for SARS-CoV-2 Mpro. Initially, supervised machine learning (ML) models were trained to predict a ligand-binding affinity of ca. 2 million compounds with the correlation on a test set of R=0.748±0.044. Atomistic simulations were then used to refine the outcome of the ML model. Using LIE/FEP calculations, nine compounds from the top 100 ML inhibitors were suggested to bind well to the protease with the domination of van der Waals interactions. Furthermore, the binding affinity of these compounds is also higher than that of nirmatrelvir, which was recently approved by the US FDA to treat COVID-19. In addition, the ligands altered the catalytic triad Cys145 - His41 - Asp187, possibly disturbing the biological activity of SARS-CoV-2.
Collapse
|
9
|
Wang L, Lu D, Wang Y, Xu X, Zhong P, Yang Z. Binding selectivity-dependent molecular mechanism of inhibitors towards CDK2 and CDK6 investigated by multiple short molecular dynamics and free energy landscapes. J Enzyme Inhib Med Chem 2023; 38:84-99. [DOI: 10.1080/14756366.2022.2135511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Affiliation(s)
- Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Dan Lu
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| | - Yan Wang
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Xiaoyan Xu
- School of Science, Shandong Jiaotong University, Jinan, PR China
| | - Peihua Zhong
- College of Computer Information and Engineering, Jiangxi Agriculture University, Nanchang, PR China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University, Nanchang, PR China
| |
Collapse
|
10
|
Pyasi S, Jonniya NA, Sk MF, Nayak D, Kar P. Finding potential inhibitors against RNA-dependent RNA polymerase (RdRp) of bovine ephemeral fever virus (BEFV): an in- silico study. J Biomol Struct Dyn 2022; 40:10403-10421. [PMID: 34238122 DOI: 10.1080/07391102.2021.1946714] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The bovine ephemeral fever virus (BEFV) is an enzootic agent that affects millions of bovines and causes major economic losses. Though the virus is seasonally reported with a very high morbidity rate (80-100%) from African, Australian, and Asiatic continents, it remains a neglected pathogen in many of its endemic areas, with no proper therapeutic drugs or vaccines presently available for treatment. The RNA-dependent RNA polymerase (RdRp) catalyzes the viral RNA synthesis and is an appropriate candidate for antiviral drug developments. We utilized integrated computational tools to build the 3D model of BEFV-RdRp and then predicted its probable active binding sites. The virtual screening and optimization against these active sites, using several small-molecule inhibitors from a different category of Life Chemical database and FDA-approved drugs from the ZINC database, was performed. We found nine molecules that have docking scores varying between -6.84 to -10.43 kcal/mol. Furthermore, these complexes were analyzed for their conformational dynamics and thermodynamic stability using molecular dynamics simulations in conjunction with the molecular mechanics generalized Born surface area (MM-GBSA) scheme. The binding free energy calculations depict that the electrostatic interactions play a dominant role in the RdRp-inhibitor binding. The hot spot residues, such as Arg565, Asp631, Glu633, Asp740, and Glu707, were found to control the RdRp-inhibitor interaction. The ADMET analysis strongly suggests favorable pharmacokinetics of these compounds that may prove useful for treating the BEFV ailment. Overall, we anticipate that these findings would help explore and develop a wide range of anti-BEFV therapy.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Shruti Pyasi
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Debasis Nayak
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| |
Collapse
|
11
|
Ferreiro D, Khalil R, Gallego MJ, Osorio NS, Arenas M. The evolution of the HIV-1 protease folding stability. Virus Evol 2022; 8:veac115. [PMID: 36601299 PMCID: PMC9802575 DOI: 10.1093/ve/veac115] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 10/10/2022] [Accepted: 12/03/2022] [Indexed: 12/11/2022] Open
Abstract
The evolution of structural proteins is generally constrained by the folding stability. However, little is known about the particular capacity of viral proteins to accommodate mutations that can potentially affect the protein stability and, in general, the evolution of the protein stability over time. As an illustrative model case, here, we investigated the evolution of the stability of the human immunodeficiency virus (HIV-1) protease (PR), which is a common HIV-1 drug target, under diverse evolutionary scenarios that include (1) intra-host virus evolution in a cohort of seventy-five patients sampled over time, (2) intra-host virus evolution sampled before and after specific PR-based treatments, and (3) inter-host evolution considering extant and ancestral (reconstructed) PR sequences from diverse HIV-1 subtypes. We also investigated the specific influence of currently known HIV-1 PR resistance mutations on the PR folding stability. We found that the HIV-1 PR stability fluctuated over time within a constant and wide range in any studied evolutionary scenario, accommodating multiple mutations that partially affected the stability while maintaining activity. We did not identify relationships between change of PR stability and diverse clinical parameters such as viral load, CD4+ T-cell counts, and a surrogate of time from infection. Counterintuitively, we predicted that nearly half of the studied HIV-1 PR resistance mutations do not significantly decrease stability, which, together with compensatory mutations, would allow the protein to adapt without requiring dramatic stability changes. We conclude that the HIV-1 PR presents a wide structural plasticity to acquire molecular adaptations without affecting the overall evolution of stability.
Collapse
Affiliation(s)
- David Ferreiro
- CINBIO, Universidade de Vigo, Vigo 36310, Spain,Departamento de Bioquímica, Genética e Inmunología, Universidade de Vigo, Vigo 36310, Spain
| | - Ruqaiya Khalil
- CINBIO, Universidade de Vigo, Vigo 36310, Spain,Departamento de Bioquímica, Genética e Inmunología, Universidade de Vigo, Vigo 36310, Spain
| | - María J Gallego
- CINBIO, Universidade de Vigo, Vigo 36310, Spain,Departamento de Bioquímica, Genética e Inmunología, Universidade de Vigo, Vigo 36310, Spain
| | - Nuno S Osorio
- Life and Health Sciences Research Institute, School of Medicine, University of Minho, Braga 4710-057, Portugal,ICVS/3Bs—PT Government Associate Laboratory, Guimarães 4806-909, Portugal
| | | |
Collapse
|
12
|
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
|
13
|
Nguyen TH, Tran PT, Pham NQA, Hoang VH, Hiep DM, Ngo ST. Identifying Possible AChE Inhibitors from Drug-like Molecules via Machine Learning and Experimental Studies. ACS OMEGA 2022; 7:20673-20682. [PMID: 35755364 PMCID: PMC9219098 DOI: 10.1021/acsomega.2c00908] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Accepted: 05/27/2022] [Indexed: 05/30/2023]
Abstract
Acetylcholinesterase (AChE) is one of the most important drug targets for Alzheimer's disease (AD) treatment. In this work, a machine learning model was trained to rapidly and accurately screen large chemical databases for the potential inhibitors of AChE. The obtained results were then validated via in vitro enzyme assay. Moreover, atomistic simulations including molecular docking and molecular dynamics simulations were then used to understand molecular insights into the binding process of ligands to AChE. In particular, two compounds including benzyl trifluoromethyl ketone and trifluoromethylstyryl ketone were indicated as highly potent inhibitors of AChE because they established IC50 values of 0.51 and 0.33 μM, respectively. The obtained IC50 of two compounds is significantly lower than that of galantamine (2.10 μM). The predicted log(BB) suggests that the compounds may be able to traverse the blood-brain barrier. A good agreement between computational and experimental studies was observed, indicating that the hybrid approach can enhance AD therapy.
Collapse
Affiliation(s)
- Trung Hai Nguyen
- Laboratory
of Theoretical and Computational Biophysics, Advanced Institute of
Materials Science, Ton Duc Thang
University, Ho Chi Minh City, Vietnam
- Faculty
of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Phuong-Thao Tran
- Hanoi
University of Pharmacy, 13-15 Le Thanh Tong, Hanoi 008404, Vietnam
| | - Ngoc Quynh Anh Pham
- Faculty
of Chemical Engineering, Ho Chi Minh City
University of Technology (HCMUT), Ho Chi Minh City 700000, Vietnam
| | - Van-Hai Hoang
- Faculty
of Pharmacy, Phenikka University, Hanoi 008404, Vietnam
- Phenikka
Institute for Advanced Study, Phenikka University, Hanoi 008404, Vietnam
| | - Dinh Minh Hiep
- Department
of Agriculture and Rural Development, Ho Chi Minh City 700000, Vietnam
| | - Son Tung Ngo
- Laboratory
of Theoretical and Computational Biophysics, Advanced Institute of
Materials Science, Ton Duc Thang
University, Ho Chi Minh City, Vietnam
- Faculty
of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| |
Collapse
|
14
|
Palli V, Leonis G, Zoupanou N, Georgiou N, Chountoulesi M, Naziris N, Tzeli D, Demetzos C, Valsami G, Marousis KD, Spyroulias GA, Mavromoustakos T. Losartan Interactions with 2-Hydroxypropyl-β-CD. Molecules 2022; 27:molecules27082421. [PMID: 35458617 PMCID: PMC9027256 DOI: 10.3390/molecules27082421] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 04/03/2022] [Accepted: 04/05/2022] [Indexed: 02/05/2023] Open
Abstract
Losartan potassium salt (LSR) is a well-known antihypertensive drug with proven beneficial effects on human health. Its formulation with the non-toxic 2-hydroxypropyl-β-cyclodextrin (2-HP-β-CD) could improve its pharmacological profile. Thus, its molecular interactions are studied using a combination of Differential Scanning Calorimetry (DSC), Nuclear Magnetic Resonance (NMR) and Molecular Dynamics (MD). First, its complexation is shown through Differential Scanning Calorimetry as lyophilization provided distinct thermal properties in comparison to the mixture. The complexation is further proved by utilizing the chemical shift changes in the complexation and T1 values. Furthermore, the reversible favorable complexation was shown by MD calculations. Such physical chemical properties provide evidence that this formulation must be further explored through biological experiments.
Collapse
Affiliation(s)
- Vasiliki Palli
- Department of Chemistry, School of Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (V.P.); (G.L.); (N.Z.); (N.G.); (D.T.)
| | - Georgios Leonis
- Department of Chemistry, School of Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (V.P.); (G.L.); (N.Z.); (N.G.); (D.T.)
| | - Nikoletta Zoupanou
- Department of Chemistry, School of Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (V.P.); (G.L.); (N.Z.); (N.G.); (D.T.)
| | - Nikitas Georgiou
- Department of Chemistry, School of Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (V.P.); (G.L.); (N.Z.); (N.G.); (D.T.)
| | - Maria Chountoulesi
- Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (M.C.); (N.N.); (C.D.)
| | - Nikolaos Naziris
- Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (M.C.); (N.N.); (C.D.)
- Department of General Biophysics, Faculty of Biology and Environmental Protection, University of Lodz, Pomorska 141/143, 90-236 Lodz, Poland
| | - Demeter Tzeli
- Department of Chemistry, School of Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (V.P.); (G.L.); (N.Z.); (N.G.); (D.T.)
- Theoretical and Physical Chemistry Institute, National Hellenic Research Foundation, 48 Vassileos Constantinou Ave., 11635 Athens, Greece
| | - Costas Demetzos
- Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (M.C.); (N.N.); (C.D.)
| | - Georgia Valsami
- Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (M.C.); (N.N.); (C.D.)
- Correspondence: (G.V.); (T.M.)
| | | | | | - Thomas Mavromoustakos
- Department of Chemistry, School of Sciences, National and Kapodistrian University of Athens, 15784 Athens, Greece; (V.P.); (G.L.); (N.Z.); (N.G.); (D.T.)
- Correspondence: (G.V.); (T.M.)
| |
Collapse
|
15
|
Wang C, Chen L, Wang R, Tang W, Zhao B. Effects of the G48M mutant on the dynamics properties and binding mechanism of PR with SQV and ATV. MOLECULAR SIMULATION 2022. [DOI: 10.1080/08927022.2022.2055013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Chao Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Lin Chen
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Ruige Wang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Wanxia Tang
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
| | - Bing Zhao
- College of Chemistry and Chemical Engineering, Qiqihar University, Qiqihar, People’s Republic of China
- Heilongjiang Provincial Key Laboratory of Surface Active Agent and Auxiliary, Qiqihar, People’s Republic of China
| |
Collapse
|
16
|
Feng X, Li F, Ding M, Zhang R, Shi T, Lu Y, Jiang W. Molecular dynamic simulation: Study on the recognition mechanism of linear β-(1 → 3)-D-glucan by Dectin-1. Carbohydr Polym 2022; 286:119276. [DOI: 10.1016/j.carbpol.2022.119276] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Revised: 02/13/2022] [Accepted: 02/18/2022] [Indexed: 12/26/2022]
|
17
|
Liang S, Liu X, Zhang S, Li M, Zhang Q, Chen J. Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations. Phys Chem Chem Phys 2022; 24:1743-1759. [PMID: 34985081 DOI: 10.1039/d1cp04361g] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (Mpro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, Mpro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to Mpro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of Mpro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to Mpro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-Mpro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against Mpro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of Mpro.
Collapse
Affiliation(s)
- Shanshan Liang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Meng Li
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan, 250358, China.
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| |
Collapse
|
18
|
Wang L, Wang Y, Zhao J, Yu Y, Kang N, Yang Z. Theoretical exploration of the binding selectivity of inhibitors to BRD7 and BRD9 with multiple short molecular dynamics simulations. RSC Adv 2022; 12:16663-16676. [PMID: 35754900 PMCID: PMC9169554 DOI: 10.1039/d2ra02637f] [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: 04/25/2022] [Accepted: 05/29/2022] [Indexed: 12/18/2022] Open
Abstract
Bromodomain-containing proteins 7 and 9 (BRD7 and BRD9) have been considered as potential targets of clinical drug design toward treatment of human cancers and other diseases. Multiple short molecular dynamics simulations and binding free energy predictions were carried out to decipher the binding selectivity of three inhibitors 4L2, 5U6, and 6KT toward BRD7 and BRD9. The results show that 4L2 has more favorable binding ability to BRD7 over BRD9 compared to 5U6 and 6KT, while 5U6 and 6KT possess more favorable associations with BRD9 than BRD7. Furthermore, estimations of residue-based free energy decompositions further identify that four common residue pairs, including (F155, F44), (V160, V49), (Y168, Y57) and (Y217, Y106) in (BRD7, BRD9) generate obvious binding differences with 4L2, 5U6, and 6KT, which mostly drives the binding selectivity of 4L2, 5U6, and 6KT to BRD7 and BRD9. Dynamic information arising from trajectory analysis also suggests that inhibitor bindings affect structural flexibility and motion modes, which is responsible for the partial selectivity of 4L2, 5U6, and 6KT toward BRD7 and BRD9. As per our expectation, this study theoretically provides useful hints for design of dual inhibitors with high selectivity on BRD7 and BRD9. Bromodomains (BRDs) are structurally conserved epigenetic reader modules observed in numerous chromatin- and transcription-associated proteins that have a capability to identify acetylated lysine residues.![]()
Collapse
Affiliation(s)
- Lifei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Yan Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Juan Zhao
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Yingxia Yu
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Nianqian Kang
- Department of Physics, Jiangxi Agricultural University, Nanchang 330045, China
| | - Zhiyong Yang
- Department of Physics, Jiangxi Agricultural University, Nanchang 330045, China
| |
Collapse
|
19
|
Guo W, Sun S, Sanchez JE, Lopez-Hernandez AE, Ale TA, Chen J, Afrin T, Qiu W, Xie Y, Li L. Using a comprehensive approach to investigate the interaction between Kinesin-5/Eg5 and the microtubule. Comput Struct Biotechnol J 2022; 20:4305-4314. [PMID: 36051882 PMCID: PMC9396395 DOI: 10.1016/j.csbj.2022.08.020] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Revised: 08/01/2022] [Accepted: 08/08/2022] [Indexed: 01/02/2023] Open
Abstract
Kinesins are microtubule-based motor proteins that play important roles ranging from intracellular transport to cell division. Human Kinesin-5 (Eg5) is essential for mitotic spindle assembly during cell division. By combining molecular dynamics (MD) simulations with other multi-scale computational approaches, we systematically studied the interaction between Eg5 and the microtubule. We find the electrostatic feature on the motor domains of Eg5 provides attractive interactions to the microtubule. Additionally, the folding and binding energy analysis reveals that the Eg5 motor domain performs its functions best when in a weak acidic environment. Molecular dynamics analyses of hydrogen bonds and salt bridges demonstrate that, on the binding interfaces of Eg5 and the tubulin heterodimer, salt bridges play the most significant role in holding the complex. The salt bridge residues on the binding interface of Eg5 are mostly positive, while salt bridge residues on the binding interface of tubulin heterodimer are mostly negative. Such salt bridge residue distribution is consistent with electrostatic potential calculations. In contrast, the interface between α and β-tubulins is dominated by hydrogen bonds rather than salt bridges. Compared to the Eg5/α-tubulin interface, the Eg5/β-tubulin interface has a greater number of salt bridges and higher occupancy for salt bridges. This asymmetric salt bridge distribution may play a significant role in Eg5′s directionality. The residues involved in hydrogen bonds and salt bridges are identified in this work and may be helpful for anticancer drug design.
Collapse
Affiliation(s)
- Wenhan Guo
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Shengjie Sun
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Jason E. Sanchez
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | | | - Tolulope A. Ale
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
| | - Jiawei Chen
- Department of Physics, University of Texas at El Paso, El Paso, TX, USA
| | - Tanjina Afrin
- Department of Physics, University of Texas at El Paso, El Paso, TX, USA
- Department of Physics, Oregon State University, Corvallis, OR, USA
| | - Weihong Qiu
- Department of Physics, Oregon State University, Corvallis, OR, USA
| | - Yixin Xie
- Department of Information Technology, Kennesaw State University, Kennesaw, GA, USA
| | - Lin Li
- Computational Science Program, University of Texas at El Paso, El Paso, TX, USA
- Department of Physics, University of Texas at El Paso, El Paso, TX, USA
- Corresponding author.
| |
Collapse
|
20
|
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
|
21
|
Sk MF, Haridev S, Roy R, Kar P. Investigating potency of TMC-126 against wild-type and mutant variants of HIV-1 protease: a molecular dynamics and free energy study. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2021; 32:941-962. [PMID: 34787532 DOI: 10.1080/1062936x.2021.1999318] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 10/25/2021] [Indexed: 06/13/2023]
Abstract
A detailed computational study was performed to investigate the conformational changes of flap region and the mechanism underlying the binding of the inhibitor TMC-126 to HIV-1 protease (PR1) and its mutant variants through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation. Further, we have studied the effectiveness of the inhibitor against HIV-2 protease (PR2). The MM-PBSA calculation suggests that TMC-126 loses its potency against mutant variants and PR2 compared to wild-type PR1 mainly due to the loss in intermolecular electrostatic interactions. The potency of the inhibitor decreases in the order: wild type PR1 > M46L > MDR20 > I50V > PR2 > V32I > A28S. Our study reveals that the flap of PR1 adopts a semi-open conformation due to the mutation I50V or MDR20. The dissimilar nature of the movement of the flap tip of both monomers is evident from the dynamic cross-correlation map. The protein structural network analysis displays that mutation causes structural rearrangements and changes the communication path between residues. Overall, we believe our study may help explore and accelerate the development of novel HIV-1/HIV-2 protease inhibitors with better potency.
Collapse
Affiliation(s)
- M F Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| | - S Haridev
- Department of Physics, Indian Institute of Technology Indore, Indore, India
| | - R Roy
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| | - P Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| |
Collapse
|
22
|
Jakhmola S, Jonniya NA, Sk MF, Rani A, Kar P, Jha HC. Identification of Potential Inhibitors against Epstein-Barr Virus Nuclear Antigen 1 (EBNA1): An Insight from Docking and Molecular Dynamic Simulations. ACS Chem Neurosci 2021; 12:3060-3072. [PMID: 34340305 DOI: 10.1021/acschemneuro.1c00350] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Epstein-Barr virus (EBV), a known tumorigenic virus, is associated with various neuropathies, including multiple sclerosis (MS). However, there is no anti-EBV FDA-approved drug available in the market. Our study targeted EBV protein EBV nuclear antigen 1 (EBNA1), crucial in virus replication and expressed in all the stages of viral latencies. This dimeric protein binds to an 18 bp palindromic DNA sequence and initiates the process of viral replication. We chose phytochemicals and FDA-approved MS drugs based on literature survey followed by their evaluation efficacies as anti-EBNA1 molecules. Molecular docking revealed FDA drugs ozanimod, siponimod, teriflunomide, and phytochemicals; emodin; protoapigenone; and EGCG bound to EBNA1 with high affinities. ADMET and Lipinski's rule analysis of the phytochemicals predicted favorable druggability. We supported our assessments of pocket druggability with molecular dynamics simulations and binding affinity predictions by the molecular mechanics generalized Born surface area (MM/GBSA) method. Our results establish a stable binding for siponimod and ozanimod with EBNA1 mainly via van der Waals interactions. We identified hot spot residues like I481', K477', L582', and K586' in the binding of ligands. In particular, K477' at the amino terminal of EBNA1 is known to establish interaction with two bases at the major groove of the DNA. Siponimod bound to EBNA1 engaging K477', thus plausibly making it unavailable for DNA interaction. Computational alanine scanning further supported the significant roles of K477', I481', and K586' in the binding of ligands with EBNA1. Conclusively, the compounds showed promising results to be used against EBNA1.
Collapse
Affiliation(s)
- Shweta Jakhmola
- Infection Bioengineering Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh 453552, India
| | - Nisha Amarnath Jonniya
- Computational Biophysics Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh 453552, India
| | - Md Fulbabu Sk
- Computational Biophysics Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh 453552, India
| | - Annu Rani
- Infection Bioengineering Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh 453552, India
| | - Parimal Kar
- Computational Biophysics Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh 453552, India
| | - Hem Chandra Jha
- Infection Bioengineering Group, Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh 453552, India
| |
Collapse
|
23
|
Dostál J, Brynda J, Vaňková L, Zia SR, Pichová I, Heidingsfeld O, Lepšík M. Structural determinants for subnanomolar inhibition of the secreted aspartic protease Sapp1p from Candida parapsilosis. J Enzyme Inhib Med Chem 2021; 36:914-921. [PMID: 33843395 PMCID: PMC8043539 DOI: 10.1080/14756366.2021.1906664] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Pathogenic Candida albicans yeasts frequently cause infections in hospitals. Antifungal drugs lose effectiveness due to other Candida species and resistance. New medications are thus required. Secreted aspartic protease of C. parapsilosis (Sapp1p) is a promising target. We have thus solved the crystal structures of Sapp1p complexed to four peptidomimetic inhibitors. Three potent inhibitors (Ki: 0.1, 0.4, 6.6 nM) resembled pepstatin A (Ki: 0.3 nM), a general aspartic protease inhibitor, in terms of their interactions with Sapp1p. However, the weaker inhibitor (Ki: 14.6 nM) formed fewer nonpolar contacts with Sapp1p, similarly to the smaller HIV protease inhibitor ritonavir (Ki: 1.9 µM), which, moreover, formed fewer H-bonds. The analyses have revealed the structural determinants of the subnanomolar inhibition of C. parapsilosis aspartic protease. Because of the high similarity between Saps from different Candida species, these results can further be used for the design of potent and specific Sap inhibitor-based antimycotic drugs.
Collapse
Affiliation(s)
- Jiří Dostál
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Jiří Brynda
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Lucie Vaňková
- Laboratory of Ligand Engineering, Institute of Biotechnology, Czech Academy of Sciences, v.v.i., BIOCEV Research Center, Vestec, Czech Republic
| | - Syeda Rehana Zia
- Department of Chemistry, University of Karachi, Karachi, Pakistan
| | - Iva Pichová
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Olga Heidingsfeld
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic.,Department of Biochemistry, Faculty of Science, Charles University in Prague, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| |
Collapse
|
24
|
Chen J, Li N, Wang X, Chen J, Zhang JZH, Zhu T. Molecular mechanism related to the binding of fluorophores to Mango-II revealed by multiple-replica molecular dynamics simulations. Phys Chem Chem Phys 2021; 23:10636-10649. [PMID: 33904542 DOI: 10.1039/d0cp06438f] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Recently, RNA aptamers activating small-molecule fluorophores have been successfully applied to tag and track RNAs in vivo. It is of significance to investigate the molecular mechanism of the fluorophore-RNA aptamer bindings at the atomic level to seek a possible pathway to enhance the fluorescence efficiency of fluorophores. In this work, multiple replica molecular dynamics (MRMD) simulations, essential dynamics (ED) analysis, and hierarchical clustering analysis were coupled to probe the effect of A22U mutation on the binding of two fluorophores, TO1-Biotin (TO1) and TO3-Biotin (TO3), to the Mango-II RNA aptamer (Mango-II). ED analysis reveals that A22U induces alterations in the binding pocket and sites of TO1 and TO3 to the Mango-II, which in turn tunes the fluorophore-RNA interface and changes the interactions of TO1 and TO3 with separate nucleotides of Mango-II. Dynamics analyses also uncover that A22U exerts the opposite impact on the molecular surface areas of the Mango-II and sugar puckers of nucleotides 22 and 23 in Mango-II complexed with TO1 and TO3. Moreover, the calculations of binding free energies suggest that A22U strengthens the binding ability of TO1 to the mutated Mango-II but weakens TO3 to the mutated Mango-II when compared with WT. These findings imply that point mutation in nucleotides possibly tune the fluorescence of fluorophores binding to RNA aptamers, providing a possible scheme to enhance the fluorescence of fluorophores.
Collapse
Affiliation(s)
- Junxiao Chen
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and School of Chemistry and Pharmaceutical Engineering, Qilu University of Technology, Jinan, 250353, People's Republic of China
| | - Na Li
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China.
| | - Xingyu Wang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China
| | - Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, People's Republic of China.
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China
| | - Tong Zhu
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, People's Republic of China. and NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, People's Republic of China and Shandong Key Laboratory of Biophysics, Institute of Biophysics, Dezhou University, Dezhou 253023, China
| |
Collapse
|
25
|
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
|
26
|
Chen J, Zhang S, Wang W, Sun H, Zhang Q, Liu X. Binding of Inhibitors to BACE1 Affected by pH-Dependent Protonation: An Exploration from Multiple Replica Gaussian Accelerated Molecular Dynamics and MM-GBSA Calculations. ACS Chem Neurosci 2021; 12:2591-2607. [PMID: 34185514 DOI: 10.1021/acschemneuro.0c00813] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
To date, inhibiting the activity of β-amyloid cleaving enzyme 1 (BACE1) has been considered an efficient approach for treating Alzheimer's disease (AD). In the current work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics general Born surface area (MM-GBSA) method were combined to investigate the effect of pH-dependent protonation on the binding of the inhibitors CS9, C6U, and 6WE to BACE1. Dynamic analyses based on the MR-GaMD trajectory show that pH-dependent protonation strongly affects the structural flexibility, correlated motions, and dynamic behavior of inhibitor-bound BACE1. According to the constructed free energy profiles, in the protonated state at low pH, inhibitor-bound BACE1 tends to populate at more conformations than in high pH. The binding free energies calculated by MM-GBSA suggest that inhibitors possess stronger binding abilities under the protonation conditions at high pH than under the protonation conditions at low pH. Moreover, pH-dependent protonation exerts a significant effect on the hydrogen bonding interactions of CS9, C6U, and 6WE to BACE1, which correspondingly alters the binding abilities of the three inhibitors to BACE1. Furthermore, in different protonated environments, three inhibitors share common interaction clusters and similar binding sites in BACE1, which are reliably used as efficient targets for the design of potent inhibitors of BACE1.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Shaolong Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Qinggang Zhang
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China
| | - Xinguo Liu
- School of Physics and Electronics, Shandong Normal University, Jinan 250358, China
| |
Collapse
|
27
|
Amarnath Jonniya N, Sk MF, Kar P. Elucidating specificity of an allosteric inhibitor WNK476 among With‐No‐Lysine kinase isoforms using molecular dynamic simulations. Chem Biol Drug Des 2021; 98:405-420. [DOI: 10.1111/cbdd.13863] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2021] [Revised: 04/19/2021] [Accepted: 05/01/2021] [Indexed: 02/06/2023]
Affiliation(s)
- Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering Indian Institute of Technology Indore Indore India
| | - Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering Indian Institute of Technology Indore Indore India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering Indian Institute of Technology Indore Indore India
| |
Collapse
|
28
|
Wang L, Wang Y, Yang Z, Xu S, Li H. Binding Selectivity of Inhibitors toward Bromodomains BAZ2A and BAZ2B Uncovered by Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations. ACS OMEGA 2021; 6:12036-12049. [PMID: 34056358 PMCID: PMC8154142 DOI: 10.1021/acsomega.1c00687] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Accepted: 04/19/2021] [Indexed: 06/12/2023]
Abstract
Two Bromodomain-Containing proteins BAZ2A and BAZ2B are responsible for remodeling chromatin and regulating noncoding RNAs. As for our current studies, integration of multiple short molecular dynamics simulations (MSMDSs) with molecular mechanics generalized Born surface area (MM-GBSA) method is adopted for insights into binding selectivity of three small molecules D8Q, D9T and UO1 to BAZ2A against BAZ2B. The calculations of MM-GBSA unveil that selectivity of inhibitors toward BAZ2A and BAZ2B highly depends on the enthalpy changes and the details uncover that D8Q has better selectivity toward BAZ2A than BAZ2B, D9T more favorably bind to BAZ2B than BAZ2A, and UO1 does not show obvious selectivity toward these two proteins. The analysis of interaction network between residues and inhibitors indicates that seven residues are mainly responsible for the selectivity of D8Q, six residues for D9T and four residues provide significant contributions to associations of UO1 with two proteins. Moreover the analysis of interaction network not only reveals warm spots of inhibitor bindings to BAZ2A and BAZ2B but also unveils that common residue pairs, including (W1816, W1887), (P1817, P1888), (F1818, F1889), (V1822, V1893), (N1823, N1894),(L1826, L1897), (V1827, V1898), (F1872, F1943), (N1873, N1944) and (V1879, I1950) belonging to (BAZ2A, BAZ2B), induce mainly binding differences of inhibitors to BAZ2A and BAZ2B. Hence, insights from our current studies offer useful dynamics information relating with conformational alterations and structure-affinity relationship at atomistic levels for novel therapeutic strategies toward BAZ2A and BAZ2B.
Collapse
Affiliation(s)
- Lifei Wang
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Yan Wang
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Zhiyong Yang
- Department
of Physics, Jiangxi Agricultural University, 1101 Zhimin Road, Economic and Technological
Development Zone, Nanchang, Jiangxi Province 330045, China
| | - Shuobo Xu
- School
of Information Science and Electrical Engineering, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| | - Hongyun Li
- School
of Science, Shandong Jiaotong University, 5001 Haitang Road, Changqing District, Jinan, Shandong Province 250357, China
| |
Collapse
|
29
|
Chen J, Wang W, Sun H, Pang L, Bao H. Binding mechanism of inhibitors to p38α MAP kinase deciphered by using multiple replica Gaussian accelerated molecular dynamics and calculations of binding free energies. Comput Biol Med 2021; 134:104485. [PMID: 33993013 DOI: 10.1016/j.compbiomed.2021.104485] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2021] [Revised: 05/07/2021] [Accepted: 05/07/2021] [Indexed: 12/12/2022]
Abstract
The p38α MAP Kinase has been an important target of drug design for treatment of inflammatory diseases and cancers. This work applies multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations and the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding mechanism of inhibitors L51, R24 and 1AU to p38α. Dynamics analyses show that inhibitor bindings exert significant effect on conformational changes of the active helix α2 and the conserved DFG loop. The rank of binding free energies calculated with MM-GBSA not only agrees well with that determined by the experimental IC50 values but also suggests that mutual compensation between the enthalpy and entropy changes can improve binding of inhibitors to p38α. The analyses of free energy landscapes indicate that the L51, R24 and 1AU bound p38α display a DFG-out conformation. The residue-based free energy decomposition method is used to evaluate contributions of separate residues to the inhibitor-p38α binding and the results imply that residues V30, V38, L74, L75, I84, T106, H107, L108, M109, L167, F169 and D168 can be utilized as efficient targets of potent inhibitors toward p38α.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan, 250357, China.
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan, 250357, China
| | - Huayin Bao
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan, 250355, China.
| |
Collapse
|
30
|
Wu S, Wang L, Zhang L, Xu X, Zhao J. Molecular dynamics simulations data of six compounds F3J-BRD4/CBP, EX1-BRD4/CBP, and E2T-BRD4/CBP. Data Brief 2021; 36:107009. [PMID: 33898668 PMCID: PMC8058528 DOI: 10.1016/j.dib.2021.107009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 03/13/2021] [Accepted: 03/24/2021] [Indexed: 11/18/2022] Open
Abstract
The data here described are related to the research article entitled “Molecular dynamics insights into binding selectivity of inhibitors toward BRD4 and CBP” [1]. Bromodomain-containing protein 4 (BRD4) and CREB binding protein (CBP) play important roles in tumorigenesis and development. We performed 200-ns molecular dynamics (MD) simulations on three pairs of inhibitor-BRD4 and inhibitor-CBP complexes to clarify binding selectivity of inhibitors toward BRD4 and CBP. Principal component (PC) analysis was used to probe changes in internal dynamics and conformations of BRD4 and CBP due to inhibitor bindings. Analysis of residue-based free energy decomposition was employed to explore the roles of separate residues in binding selectivity of inhibitors to BRD4 versus CBP.
Collapse
|
31
|
Jonniya NA, Sk MF, Kar P. Characterizing an allosteric inhibitor-induced inactive state in with-no-lysine kinase 1 using Gaussian accelerated molecular dynamics simulations. Phys Chem Chem Phys 2021; 23:7343-7358. [DOI: 10.1039/d0cp05733a] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
The binding of an allosteric inhibitor in WNK1 leads to the inactive state.
Collapse
Affiliation(s)
- Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, MP
- India
| | - Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, MP
- India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, MP
- India
| |
Collapse
|
32
|
Leonis G, Ntountaniotis D, Christodoulou E, Mavromoustakos T. Molecular Dynamics Protocols for the Study of Cyclodextrin Drug Delivery Systems. Methods Mol Biol 2021; 2207:109-125. [PMID: 33113131 DOI: 10.1007/978-1-0716-0920-0_9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Hypertension treatment is a current therapeutic priority as there is a constantly increasing part of the population that suffers from this risk factor, which may lead to cardiovascular and encephalic episodes and eventually to death. A number of marketed medicines consist of active ingredients that may be relatively potent; however, there is plenty of room to enhance their pharmacological profile and therapeutic index by improving specific physicochemical properties. In this work, we focus on a class of blood pressure regulators, called sartans, and we present the computational scheme for the pharmacological improvement of irbesartan (IRB) as a representative example. IRB has been shown to exert increased pharmacological action compared with other sartans, but it appears to be highly lipophilic and violates Lipinski rule (MLogP >4.15). To circumvent this drawback, proper hydrophilic molecules, such as cyclodextrins, can be used as drug carriers. This chapter describes the combinatory use of computational methods, namely molecular docking, quantum mechanics, molecular dynamics, and free energy calculations, to study the interactions and the energetic contributions that govern the IRB:cyclodextrin association. We provide a detailed computational protocol, which aims to assist the improvement of the pharmacological properties of sartans. This protocol can also be applied to any other drug molecule with diminished hydrophilic character.
Collapse
Affiliation(s)
- Georgios Leonis
- Department of Chemistry, National and Kapodistrian University of Athens, Zografou, Greece.
| | | | - Eirini Christodoulou
- Department of Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| | - Thomas Mavromoustakos
- Department of Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| |
Collapse
|
33
|
Falls Z, Fine J, Chopra G, Samudrala R. Accurate Prediction of Inhibitor Binding to HIV-1 Protease Using CANDOCK. Front Chem 2021; 9:775513. [PMID: 35111726 PMCID: PMC8801943 DOI: 10.3389/fchem.2021.775513] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Accepted: 11/25/2021] [Indexed: 12/27/2022] Open
Abstract
The human immunodeficiency virus 1 (HIV-1) protease is an important target for treating HIV infection. Our goal was to benchmark a novel molecular docking protocol and determine its effectiveness as a therapeutic repurposing tool by predicting inhibitor potency to this target. To accomplish this, we predicted the relative binding scores of various inhibitors of the protease using CANDOCK, a hierarchical fragment-based docking protocol with a knowledge-based scoring function. We first used a set of 30 HIV-1 protease complexes as an initial benchmark to optimize the parameters for CANDOCK. We then compared the results from CANDOCK to two other popular molecular docking protocols Autodock Vina and Smina. Our results showed that CANDOCK is superior to both of these protocols in terms of correlating predicted binding scores to experimental binding affinities with a Pearson coefficient of 0.62 compared to 0.48 and 0.49 for Vina and Smina, respectively. We further leveraged the Database of Useful Decoys: Enhanced (DUD-E) HIV protease set to ascertain the effectiveness of each protocol in discriminating active versus decoy ligands for proteases. CANDOCK again displayed better efficacy over the other commonly used molecular docking protocols with area under the receiver operating characteristic curve (AUROC) of 0.94 compared to 0.71 and 0.74 for Vina and Smina. These findings support the utility of CANDOCK to help discover novel therapeutics that effectively inhibit HIV-1 and possibly other retroviral proteases.
Collapse
Affiliation(s)
- Zackary Falls
- Department of Biomedical Informatics, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo, State University of New York, Buffalo, NY, United States
| | - Jonathan Fine
- Department of Chemistry, Purdue University, West Lafayette, IN, United States
| | - Gaurav Chopra
- Department of Chemistry, Purdue University, West Lafayette, IN, United States.,Purdue Institute for Drug Discovery, West Lafayette, IN, United States.,Purdue Center for Cancer Research, West Lafayette, IN, United States.,Purdue Institute for Inflammation, Immunology and Infectious Disease, West Lafayette, IN, United States.,Purdue Institute for Integrative Neuroscience, West Lafayette, IN, United States
| | - Ram Samudrala
- Department of Biomedical Informatics, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo, State University of New York, Buffalo, NY, United States
| |
Collapse
|
34
|
Ngo ST, Quynh Anh Pham N, Thi Le L, Pham DH, Vu VV. Computational Determination of Potential Inhibitors of SARS-CoV-2 Main Protease. J Chem Inf Model 2020; 60:5771-5780. [PMID: 32530282 PMCID: PMC7323056 DOI: 10.1021/acs.jcim.0c00491] [Citation(s) in RCA: 94] [Impact Index Per Article: 23.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Indexed: 12/13/2022]
Abstract
The novel coronavirus (SARS-CoV-2) has infected several million people and caused thousands of deaths worldwide since December 2019. As the disease is spreading rapidly all over the world, it is urgent to find effective drugs to treat the virus. The main protease (Mpro) of SARS-CoV-2 is one of the potential drug targets. Therefore, in this context, we used rigorous computational methods, including molecular docking, fast pulling of ligand (FPL), and free energy perturbation (FEP), to investigate potential inhibitors of SARS-CoV-2 Mpro. We first tested our approach with three reported inhibitors of SARS-CoV-2 Mpro, and our computational results are in good agreement with the respective experimental data. Subsequently, we applied our approach on a database of ∼4600 natural compounds, as well as 8 available HIV-1 protease (PR) inhibitors and an aza-peptide epoxide. Molecular docking resulted in a short list of 35 natural compounds, which was subsequently refined using the FPL scheme. FPL simulations resulted in five potential inhibitors, including three natural compounds and two available HIV-1 PR inhibitors. Finally, FEP, the most accurate and precise method, was used to determine the absolute binding free energy of these five compounds. FEP results indicate that two natural compounds, cannabisin A and isoacteoside, and an HIV-1 PR inhibitor, darunavir, exhibit a large binding free energy to SARS-CoV-2 Mpro, which is larger than that of 13b, the most reliable SARS-CoV-2 Mpro inhibitor recently reported. The binding free energy largely arises from van der Waals interaction. We also found that Glu166 forms H-bonds to all of the inhibitors. Replacing Glu166 by an alanine residue leads to ∼2.0 kcal/mol decreases in the affinity of darunavir to SARS-CoV-2 Mpro. Our results could contribute to the development of potential drugs inhibiting SARS-CoV-2.
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and
Computational Biophysics, Ton Duc Thang
University, Ho Chi Minh City 700000,
Vietnam
- Faculty of Applied Sciences,
Ton Duc Thang University, Ho Chi Minh
City 700000, Vietnam
| | - Ngoc Quynh Anh Pham
- Faculty of Chemical Engineering,
Ho Chi Minh City University of Technology
(HCMUT), Ho Chi Minh City 700000,
Vietnam
| | - Ly Thi Le
- School of Biotechnology,
International University, Ho Chi Minh
Ciy 700000, Vietnam
| | - Duc-Hung Pham
- Division of Immunobiology,
Cincinnati Children’s Hospital Medical
Center, Cincinnati, Ohio 45229, United
States
| | - Van V. Vu
- NTT Hi-Tech Institute, Nguyen
Tat Thanh University, Ho Chi Minh City 700000,
Vietnam
| |
Collapse
|
35
|
Prediction of Drug Potencies of BACE1 Inhibitors: A Molecular Dynamics Simulation and MM_GB(PB)SA Scoring. COMPUTATION 2020. [DOI: 10.3390/computation8040106] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Alzheimer’s disease (AD) is a progressive neurodegenerative brain disorder. One of the important therapeutic approaches of AD is the inhibition of β-site APP cleaving enzyme-1 (BACE1). This enzyme plays a central role in the synthesis of the pathogenic β-amyloid peptides (Aβ) in Alzheimer’s disease. A group of potent BACE1 inhibitors with known X-ray structures (PDB ID 5i3X, 5i3Y, 5iE1, 5i3V, 5i3W, 4LC7, 3TPP) were studied by molecular dynamics simulation and binding energy calculation employing MM_GB(PB)SA. The calculated binding energies gave Kd values of 0.139 µM, 1.39 nM, 4.39 mM, 24.3 nM, 1.39 mM, 29.13 mM, and 193.07 nM, respectively. These inhibitors showed potent inhibitory activities in enzymatic and cell assays. The Kd values are compared with experimental values and the structures are discussed in view of the energy contributions to binding. Drug likeness of these inhibitors is also discussed. Accommodation of ligands in the catalytic site of BACE1 is discussed depending on the type of fragment involved in each structure. Molecular dynamics (MD) simulations and energy studies were used to explore the recognition of the selected BACE1 inhibitors by Asp32, Asp228, and the hydrophobic flap. The results show that selective BACE1 inhibition may be due to the formation of strong electrostatic interactions with Asp32 and Asp228 and a large number of hydrogen bonds, in addition to π–π and van der Waals interactions with the amino acid residues located inside the catalytic cavity. Interactions with the ligands show a similar binding mode with BACE1. These results help to rationalize the design of selective BACE1 inhibitors.
Collapse
|
36
|
Chen J, Wang W, Sun H, Pang L, Yin B. Mutation-mediated influences on binding of anaplastic lymphoma kinase to crizotinib decoded by multiple replica Gaussian accelerated molecular dynamics. J Comput Aided Mol Des 2020; 34:1289-1305. [DOI: 10.1007/s10822-020-00355-5] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2020] [Accepted: 10/14/2020] [Indexed: 12/19/2022]
|
37
|
Jonniya NA, Sk MF, Kar P. A comparative study of structural and conformational properties of WNK kinase isoforms bound to an inhibitor: insights from molecular dynamic simulations. J Biomol Struct Dyn 2020; 40:1400-1415. [DOI: 10.1080/07391102.2020.1827035] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Nisha Amarnath Jonniya
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Md Fulbabu Sk
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| | - Parimal Kar
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, Madhya Pradesh, India
| |
Collapse
|
38
|
Discovery of a stable tripeptide targeting the N-domain of CRF1 receptor. Amino Acids 2020; 52:1337-1351. [PMID: 32996057 DOI: 10.1007/s00726-020-02895-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2019] [Accepted: 09/19/2020] [Indexed: 10/23/2022]
Abstract
The corticotropin-releasing factor (CRF) and its CRF1 receptor (CRF1R) play a central role in the maintenance of homeostasis. Malfunctioning of the CRF/CRF1R unit is associated with several disorders, such as anxiety and depression. Non-peptide CRF1R-selective antagonists have been shown to exert anxiolytic and antidepressant effects on experimental animals. However, none of them is in clinical use today because of several side effects, thus demonstrating the need for the development of other more suitable CRF1R antagonists. In an effort to develop novel CRF1R antagonists we designed, synthesized and chemically characterized two tripeptide analogues of CRF, namely (R)-LMI and (S)-LMI, having their Leu either in R (or D) or in S (or L) configuration, respectively. Their design was based on the crystal structure of the N-extracellular domain (N-domain) of CRF1R/CRF complex, using a relevant array of computational methods. Experimental evaluation of the stability of synthetic peptides in human plasma has revealed that (R)-LMI is proteolytically more stable than (S)-LMI. Based on this finding, (R)-LMI was selected for pharmacological characterization. We have found that (R)-LMI is a CRF antagonist, inhibiting (1) the CRF-stimulated accumulation of cAMP in HEK 293 cells expressing the CRF1R, (2) the production of interleukins by adipocytes and (3) the proliferation rate of RAW 264.7 cells. (R)-LMI likely blocked agonist actions by interacting with the N-domain of CRF1R as suggested by data using a constitutively active chimera of CRF1R. We propose that (R)-LMI can be used as an optimal lead compound in the rational design of novel CRF antagonists.
Collapse
|
39
|
Sk MF, Jonniya NA, Kar P. Exploring the energetic basis of binding of currently used drugs against HIV-1 subtype CRF01_AE protease via molecular dynamics simulations. J Biomol Struct Dyn 2020; 39:5892-5909. [DOI: 10.1080/07391102.2020.1794965] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Md Fulbabu Sk
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| | - Nisha Amarnath Jonniya
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| | - Parimal Kar
- Discipline of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| |
Collapse
|
40
|
Leonis G, Christodoulou E, Ntountaniotis D, Chatziathanasiadou MV, Mavromoustakos T, Naziris N, Chountoulesi M, Demetzos C, Valsami G, Damalas DE, Tzakos AG, Thomaidis NS, Karageorgos V, Liapakis G. Antihypertensive activity and molecular interactions of irbesartan in complex with 2-hydroxypropyl-β-cyclodextrin. Chem Biol Drug Des 2020; 96:668-683. [PMID: 32691965 DOI: 10.1111/cbdd.13664] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2019] [Revised: 11/28/2019] [Accepted: 12/07/2019] [Indexed: 12/19/2022]
Abstract
Irbesartan (IRB) exerts beneficial effects either alone or in combination with other drugs on numerous diseases, such as cancer, diabetes, and hypertension. However, due to its high lipophilicity, IRB does not possess the optimum pharmacological efficiency. To circumvent this problem, a drug delivery system with 2-hydroxypropyl-β-cyclodextrin (2-HP-β-CD) was explored. The 1:1 complex between IRB and 2-HP-β-CD was identified through ESI QTF HRMS. Dissolution studies showed a higher dissolution rate of the lyophilized IRB-2-HP-β-CD complex than the tablet containing IRB at pH = 1.2. DSC results revealed the differences of the thermal properties between the complex and various mixtures consisting of the two components, namely IRB and 2-HP-β-CD. Interestingly, depending on the way the mixture preparation was conducted, different association between the two components was observed. Molecular dynamics (MD) simulations predicted the favorable formation of the above complex and identified the dominant interactions between IRB and 2-HP-β-CD. In vitro pharmacological results verified that the inclusion complex not only preserves the binding affinity of IRB for AT1R receptor, but also it slightly increases it. As the complex formulation lacks the problems of the tablet, our approach is a promising new way to improve the efficiency of IRB.
Collapse
Affiliation(s)
- Georgios Leonis
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| | - Eirini Christodoulou
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| | - Dimitrios Ntountaniotis
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| | - Maria V Chatziathanasiadou
- Department of Chemistry, Section of Organic Chemistry and Biochemistry, University of Ioannina, Ioannina, Greece
| | - Thomas Mavromoustakos
- Laboratory of Organic Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| | - Nikolaos Naziris
- Section of Pharmaceutical Technology, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Athens, Greece
| | - Maria Chountoulesi
- Section of Pharmaceutical Technology, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Athens, Greece
| | - Costas Demetzos
- Section of Pharmaceutical Technology, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Athens, Greece
| | - Georgia Valsami
- Section of Pharmaceutical Technology, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Athens, Greece
| | - Dimitrios E Damalas
- Department of Chemistry, Laboratory of Analytical Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| | - Andreas G Tzakos
- Department of Chemistry, Section of Organic Chemistry and Biochemistry, University of Ioannina, Ioannina, Greece
| | - Nikolaos S Thomaidis
- Department of Chemistry, Laboratory of Analytical Chemistry, National and Kapodistrian University of Athens, Zografou, Greece
| | - Vlasios Karageorgos
- Department of Basic Sciences, School of Medicine, University of Crete, Heraklion, Crete, Greece
| | - Georgios Liapakis
- Department of Basic Sciences, School of Medicine, University of Crete, Heraklion, Crete, Greece
| |
Collapse
|
41
|
Wu SL, Wang LF, Sun HB, Wang W, Yu YX. Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2020; 31:547-570. [PMID: 32657160 DOI: 10.1080/1062936x.2020.1777584] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Accepted: 05/31/2020] [Indexed: 06/11/2023]
Abstract
It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.
Collapse
Affiliation(s)
- S L Wu
- School of Science, Shandong Jiaotong University , Jinan, China
| | - L F Wang
- School of Science, Shandong Jiaotong University , Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University , Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University , Jinan, China
| | - Y X Yu
- School of Science, Shandong Jiaotong University , Jinan, China
| |
Collapse
|
42
|
Rahman SU, Ali HS, Jafari B, Zaib S, Hameed A, Al-Kahraman YMSA, Langer P, Iqbal J. Structure-based virtual screening of dipeptidyl peptidase 4 inhibitors and their in vitro analysis. Comput Biol Chem 2020; 91:107326. [PMID: 32739275 DOI: 10.1016/j.compbiolchem.2020.107326] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2019] [Revised: 06/23/2020] [Accepted: 06/25/2020] [Indexed: 10/24/2022]
Abstract
Type 2 diabetes mellitus (T2DM) is one of the most widely prevalent metabolic disorders with no cure to date thus remains the most challenging task in the current drug discovery. Therefore, the only strategy to control diabetes prevalence is to develop novel efficacious therapeutics. Dipeptidyl Peptidase 4 (DPP-4) inhibitors are currently used as anti-diabetic drugs for the inhibition of incretins. This study aims to construct the chemical feature based on pharmacophore models for dipeptidyl peptidase IV. The structure-based pharmacophore modeling has been employed to evaluate new inhibitors of DPP-4. A four-featured pharmacophore model was developed from crystal structure of DPP-4 enzyme with 4-(2-aminoethyl) benzenesulfonyl fluoride in its active site via pharmacophore constructing tool of Molecular Operating Environment (MOE) consisting F1 Hyd (hydrophobic region), F2 Hyd|Cat|Don (hydrophobic cationic and donor region), F3 Acc (acceptor region) and F4 Hyd (hydrophobic region). The generated pharmacophore model was used for virtual screening of in-house compound library (the available compounds which were used for initial screening to get the few compounds for the current studies). The resultant selected compounds, after virtual screening were further validated using in vitro assay. Furthermore, structure-activity relationship was carried out for the compounds possessing significant inhibition potential after docking studies. The binding free energy of analogs was evaluated via molecular mechanics generalized Born surface area (MM-GBSA) and Poisson-Boltzmann surface area (MM-PBSA) methods using AMBER 16 as a molecular dynamics (MD) simulation package. Based on potential findings, we report that selected candidates are more likely to be used as DPP-4 inhibitors or as starting leads for the development of novel and potent DPP-4 inhibitors.
Collapse
Affiliation(s)
- Shafiq Ur Rahman
- Centre for Advanced Drug Research, COMSATS University Islamabad, Abbottabad Campus, Abbottabad, 22060, Pakistan
| | - Hafiz Saqib Ali
- Department of Chemistry & Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, United Kingdom
| | - Behzad Jafari
- Institut für Chemie, Universität Rostock, Albert-Einstein-Str. 3a, 18059 Rostock, Germany
| | - Sumera Zaib
- Centre for Advanced Drug Research, COMSATS University Islamabad, Abbottabad Campus, Abbottabad, 22060, Pakistan
| | - Abdul Hameed
- Centre for Advanced Drug Research, COMSATS University Islamabad, Abbottabad Campus, Abbottabad, 22060, Pakistan
| | - Yasser M S A Al-Kahraman
- Centre for Advanced Drug Research, COMSATS University Islamabad, Abbottabad Campus, Abbottabad, 22060, Pakistan
| | - Peter Langer
- Institut für Chemie, Universität Rostock, Albert-Einstein-Str. 3a, 18059 Rostock, Germany; Leibniz Institut für Katalyse an der Universität Rostock e.V. (LIKAT), Albert-Einstein-Str. 29a, 18059 Rostock, Germany
| | - Jamshed Iqbal
- Centre for Advanced Drug Research, COMSATS University Islamabad, Abbottabad Campus, Abbottabad, 22060, Pakistan.
| |
Collapse
|
43
|
Chen J, Yin B, Wang W, Sun H. Effects of Disulfide Bonds on Binding of Inhibitors to β-Amyloid Cleaving Enzyme 1 Decoded by Multiple Replica Accelerated Molecular Dynamics Simulations. ACS Chem Neurosci 2020; 11:1811-1826. [PMID: 32459964 DOI: 10.1021/acschemneuro.0c00234] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The β-amyloid cleaving enzyme 1 (BACE1) has been thought to be an efficient target for treatment of Alzheimer's disease (AD). Deep insight into inhibitor-BACE1 binding mechanism is of significance for design of potent drugs toward BACE1. In this work, multiple replica accelerated molecular dynamics (MR-aMD) simulations, principal component (PC) analysis, and free energy landscapes were integrated to decode the effect of disulfide bonds (SSBs) in BACE1 on bindings of three inhibitors 3KO, 3KT, and 779 to BACE1. The results from cross-correlation analysis suggest that the breaking of SSBs exerts significant influence on structural flexibility and internal dynamics of inhibitor-bound BACE1. PC analysis and free energy landscapes reveal that the breaking of SSBs not only evidently induces the conformational rearrangement of BACE1 but also highly changes binding poses of three inhibitors in BACE1 and leads to more disordered binding of three inhibitors to BACE1, which is further supported by the increase in binding entropy of inhibitors to BACE1 due to the breaking of SSBs. Residue-based free energy decomposition method was utilized to evaluate contributions of separate residues to inhibitor-BACE1 binding. The results suggest that although the breaking of SSBs in BACE1 does not destroy the interaction network of inhibitors with BACE1 it changes interaction strength of some residues with inhibitors. Meanwhile, the information from residue-based free energy decomposition indicates that residues L91, S96, V130, Y132, Q134, W137, F169, I171, and I179 can be used as efficient targets of drug design toward BACE1.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Baohua Yin
- School of Pharmacy, Shandong University of Traditional Chinese Medicine, Jinan 250355, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Haibo Sun
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| |
Collapse
|
44
|
Khan MT, Ali S, Zeb MT, Kaushik AC, Malik SI, Wei DQ. Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance. Front Mol Biosci 2020; 7:52. [PMID: 32328498 PMCID: PMC7160322 DOI: 10.3389/fmolb.2020.00052] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 03/16/2020] [Indexed: 12/16/2022] Open
Abstract
A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Sajid Ali
- Department of Microbiology, Quaid-i-Azam University Islamabad, Islamabad, Pakistan
| | | | - Aman Chandra Kaushik
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
| | - Shaukat Iqbal Malik
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Dong-Qing Wei
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
- Peng Cheng Laboratory, Shenzhen, China
| |
Collapse
|
45
|
Whitfield TW, Ragland DA, Zeldovich KB, Schiffer CA. Characterizing Protein-Ligand Binding Using Atomistic Simulation and Machine Learning: Application to Drug Resistance in HIV-1 Protease. J Chem Theory Comput 2020; 16:1284-1299. [PMID: 31877249 DOI: 10.1021/acs.jctc.9b00781] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Over the past several decades, atomistic simulations of biomolecules, whether carried out using molecular dynamics or Monte Carlo techniques, have provided detailed insights into their function. Comparing the results of such simulations for a few closely related systems has guided our understanding of the mechanisms by which changes such as ligand binding or mutation can alter the function. The general problem of detecting and interpreting such mechanisms from simulations of many related systems, however, remains a challenge. This problem is addressed here by applying supervised and unsupervised machine learning techniques to a variety of thermodynamic observables extracted from molecular dynamics simulations of different systems. As an important test case, these methods are applied to understand the evasion by human immunodeficiency virus type-1 (HIV-1) protease of darunavir, a potent inhibitor to which resistance can develop via the simultaneous mutation of multiple amino acids. Complex mutational patterns have been observed among resistant strains, presenting a challenge to developing a mechanistic picture of resistance in the protease. In order to dissect these patterns and gain mechanistic insight into the role of specific mutations, molecular dynamics simulations were carried out on a collection of HIV-1 protease variants, chosen to include highly resistant strains and susceptible controls, in complex with darunavir. Using a machine learning approach that takes advantage of the hierarchical nature in the relationships among the sequence, structure, and function, an integrative analysis of these trajectories reveals key details of the resistance mechanism, including changes in the protein structure, hydrogen bonding, and protein-ligand contacts.
Collapse
Affiliation(s)
- Troy W Whitfield
- Department of Medicine , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States.,Program in Bioinformatics and Integrative Biology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| | - Debra A Ragland
- Department of Biochemistry and Molecular Pharmacology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| | - Konstantin B Zeldovich
- Program in Bioinformatics and Integrative Biology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| | - Celia A Schiffer
- Department of Biochemistry and Molecular Pharmacology , University of Massachusetts Medical School , Worcester , Massachusetts 01605 , United States
| |
Collapse
|
46
|
Chen J, Wang X, Pang L, Zhang JZH, Zhu T. Effect of mutations on binding of ligands to guanine riboswitch probed by free energy perturbation and molecular dynamics simulations. Nucleic Acids Res 2020; 47:6618-6631. [PMID: 31173143 PMCID: PMC6649850 DOI: 10.1093/nar/gkz499] [Citation(s) in RCA: 118] [Impact Index Per Article: 29.5] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2019] [Revised: 05/22/2019] [Accepted: 05/29/2019] [Indexed: 12/14/2022] Open
Abstract
Riboswitches can regulate gene expression by direct and specific interactions with ligands and have recently attracted interest as potential drug targets for antibacterial. In this work, molecular dynamics (MD) simulations, free energy perturbation (FEP) and molecular mechanics generalized Born surface area (MM-GBSA) methods were integrated to probe the effect of mutations on the binding of ligands to guanine riboswitch (GR). The results not only show that binding free energies predicted by FEP and MM-GBSA obtain an excellent correlation, but also indicate that mutations involved in the current study can strengthen the binding affinity of ligands GR. Residue-based free energy decomposition was applied to compute ligand-nucleotide interactions and the results suggest that mutations highly affect interactions of ligands with key nucleotides U22, U51 and C74. Dynamics analyses based on MD trajectories indicate that mutations not only regulate the structural flexibility but also change the internal motion modes of GR, especially for the structures J12, J23 and J31, which implies that the aptamer domain activity of GR is extremely plastic and thus readily tunable by nucleotide mutations. This study is expected to provide useful molecular basis and dynamics information for the understanding of the function of GR and possibility as potential drug targets for antibacterial.
Collapse
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
| | - Laixue Pang
- School of Science, Shandong Jiaotong University, Jinan 250357 China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Tong Zhu
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| |
Collapse
|
47
|
Chen J, Liu X, Zhang S, Chen J, Sun H, Zhang L, Zhang Q. Molecular mechanism with regard to the binding selectivity of inhibitors toward FABP5 and FABP7 explored by multiple short molecular dynamics simulations and free energy analyses. Phys Chem Chem Phys 2020; 22:2262-2275. [DOI: 10.1039/c9cp05704h] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Recently, fatty acid binding proteins 5 and 7 (FABP5 and FABP7) have been regarded as the prospective targets for clinically treating multiple diseases related to FABPs.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Xinguo Liu
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| | - Shaolong Zhang
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| | - Junxiao Chen
- School of Chemistry and Pharmaceutical Engineering
- Qilu University of Technology
- Jinan
- People's Republic of China
| | - Haibo Sun
- School of Science
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Lin Zhang
- School of Construction Machinery
- Shandong Jiaotong University
- Jinan 250357
- People's Republic of China
| | - Qinggang Zhang
- School of Physics and Electronics
- Shandong Normal University
- Jinan
- People's Republic of China
| |
Collapse
|
48
|
Ngo ST, Hong ND, Quynh Anh LH, Hiep DM, Tung NT. Effective estimation of the inhibitor affinity of HIV-1 protease via a modified LIE approach. RSC Adv 2020; 10:7732-7739. [PMID: 35492181 PMCID: PMC9049864 DOI: 10.1039/c9ra09583g] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2019] [Accepted: 02/06/2020] [Indexed: 01/07/2023] Open
Abstract
The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses. Computer-aided drug design (CADD) would enhance the discovery of new therapies, through which the estimation of ligand-binding affinity is critical to predict the most efficient inhibitor. A time-consuming binding free energy method would reduce the usefulness of CADD. The modified linear interaction energy (LIE) approach emerges as an appropriate protocol that performs this task. In particular, the polar interaction free energy, which is obtained via numerically resolving the linear Poisson–Boltzmann equation, plays as an important role in driving the binding mechanism of the HIV-1 PR + inhibitor complex. The electrostatic interaction energy contributes to the attraction between two molecules, but the vdW interaction acts as a repulsive factor between the ligand and the HIV-1 PR. Moreover, the ligands were found to adopt a very strong hydrophobic interaction with the HIV-1 PR. Furthermore, the results obtained corroborate the high accuracy and precision of computational studies with a large correlation coefficient value R = 0.83 and a small RMSE δRMSE = 1.25 kcal mol−1. This method is less time-consuming than the other end-point methods, such as the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and free energy perturbation (FEP) approaches. Overall, the modified LIE approach would provide ligand-binding affinity with HIV-1 PR accurately, precisely, and rapidly, resulting in a more efficient design of new inhibitors. The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses.![]()
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics
- Ton Duc Thang University
- Ho Chi Minh City
- Vietnam
- Faculty of Applied Sciences
| | - Nam Dao Hong
- University of Medicine and Pharmacy at Ho Chi Minh City
- Ho Chi Minh City
- Vietnam
| | - Le Huu Quynh Anh
- Department of Climate Change and Renewable Energy
- Ho Chi Minh City University of Natural Resources and Environment
- Ho Chi Minh City
- Vietnam
| | | | - Nguyen Thanh Tung
- Institute of Materials Science & Graduate University of Science and Technology
- Vietnam Academy of Science and Technology
- Hanoi
- Vietnam
| |
Collapse
|
49
|
Nguyen NT, Nguyen TH, Pham TNH, Huy NT, Bay MV, Pham MQ, Nam PC, Vu VV, Ngo ST. Autodock Vina Adopts More Accurate Binding Poses but Autodock4 Forms Better Binding Affinity. J Chem Inf Model 2019; 60:204-211. [DOI: 10.1021/acs.jcim.9b00778] [Citation(s) in RCA: 131] [Impact Index Per Article: 26.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Nguyen Thanh Nguyen
- Department of Theoretical Physics, Ho Chi Minh City University of Science, Ho Chi Minh City 700000, Vietnam
| | - Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
| | - T. Ngoc Han Pham
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
| | - Nguyen Truong Huy
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
| | - Mai Van Bay
- Department of Chemical Engineering, The University of Da Nang, University of Science and Technology, Da Nang City 550000, Vietnam
| | - Minh Quan Pham
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi 100000, Vietnam
| | - Pham Cam Nam
- Department of Chemical Engineering, The University of Da Nang, University of Science and Technology, Da Nang City 550000, Vietnam
| | - Van V. Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City 700000, Vietnam
| | - Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
| |
Collapse
|
50
|
Theoretical analyses on enantiospecificity of L-2-haloacid dehalogenase (DehL) from Rhizobium sp. RC1 towards 2-chloropropionic acid. J Mol Graph Model 2019; 92:131-139. [DOI: 10.1016/j.jmgm.2019.07.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2019] [Revised: 07/02/2019] [Accepted: 07/20/2019] [Indexed: 11/24/2022]
|