1
|
Guo F, Yang H, Li S, Jiang Y, Bai X, Hu C, Li W, Han W. Using Gaussian accelerated molecular dynamics combined with Markov state models to explore the mechanism of action of new oral inhibitors on Complex I. Comput Biol Med 2024; 177:108598. [PMID: 38776729 DOI: 10.1016/j.compbiomed.2024.108598] [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: 03/18/2024] [Revised: 04/15/2024] [Accepted: 05/11/2024] [Indexed: 05/25/2024]
Abstract
In this study, our focus was on investigating H-1,2,3-triazole derivative HP661 as a novel and highly efficient oral OXPHOS inhibitor, with its molecular-level inhibitory mechanism not yet fully understood. We selected the ND1, NDUFS2, and NDUFS7 subunits of Mitochondrial Complex I as the receptor proteins and established three systems for comparative analysis: protein-IACS-010759, protein-lead compound 10, and protein-HP661. Through extensive analysis involving 500 ns Gaussian molecular dynamics simulations, we gained insights into these systems. Additionally, we constructed a Markov State Models to examine changes in secondary structures during the motion processes. The research findings suggest that the inhibitor HP661 enhances the extensibility and hydrophilicity of the receptor protein. Furthermore, HP661 induces the unwinding of the α-helical structure in the region of residues 726-730. Notably, key roles were identified for Met37, Phe53, and Pro212 in the binding of various inhibitors. In conclusion, we delved into the potential molecular mechanisms of triazole derivative HP661 in inhibiting Complex I. These research outcomes provide crucial information for a deeper understanding of the mechanisms underlying OXPHOS inhibition, offering valuable theoretical support for drug development and disease treatment design.
Collapse
Affiliation(s)
- Fangfang Guo
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China
| | - Hengzheng Yang
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China
| | - Shihong Li
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China
| | - Yongxin Jiang
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China
| | - Xue Bai
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China
| | - Chengxiang Hu
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China
| | - Wannan Li
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China.
| | - Weiwei Han
- Edmond H. Fischer Signal Transduction Laboratory and Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun, 130012, China.
| |
Collapse
|
2
|
Liu X, Sun S, Liu J, Dang Q, Gao Y, Fang L, Min W. Isolation, Virtual Screening, and Evaluation of Hazelnut-Derived Immunoactive Peptides for the Inhibition of SARS-CoV-2 Main Protease. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2024; 72:11561-11576. [PMID: 38739709 DOI: 10.1021/acs.jafc.4c01942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
The aim of this study is to validate the activity of hazelnut (Corylus avellana L.)-derived immunoactive peptides inhibiting the main protease (Mpro) of SARS-CoV-2 and further unveil their interaction mechanism using in vitro assays, molecular dynamics (MD) simulations, and binding free energy calculations. In general, the enzymatic hydrolysis components, especially molecular weight < 3 kDa, possess good immune activity as measured by the proliferation ability of mouse splenic lymphocytes and phagocytic activity of mouse peritoneal macrophages. Over 866 unique peptide sequences were isolated, purified, and then identified by nanohigh-performance liquid chromatography/tandem mass spectrometry (NANO-HPLC-MS/MS) from hazelnut protein hydrolysates, but Trp-Trp-Asn-Leu-Asn (WWNLN) and Trp-Ala-Val-Leu-Lys (WAVLK) in particular are found to increase the cell viability and phagocytic capacity of RAW264.7 macrophages as well as promote the secretion of the cytokines nitric oxide (NO), tumor necrosis factor-α (TNF-α), and interleukin-1β (IL-1β). Fluorescence resonance energy transfer assay elucidated that WWNLN and WAVLK exhibit excellent inhibitory potency against Mpro, with IC50 values of 6.695 and 16.750 μM, respectively. Classical all-atom MD simulations show that hydrogen bonds play a pivotal role in stabilizing the complex conformation and protein-peptide interaction. Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) calculation indicates that WWNLN has a lower binding free energy with Mpro than WAVLK. Furthermore, adsorption, distribution, metabolism, excretion, and toxicity (ADMET) predictions illustrate favorable drug-likeness and pharmacokinetic properties of WWNLN compared to WAVLK. This study provides a new understanding of the immunomodulatory activity of hazelnut hydrolysates and sheds light on peptide inhibitors targeting Mpro.
Collapse
Affiliation(s)
- Xiaoting Liu
- College of Food Science and Engineering, National Engineering Laboratory of Wheat and Corn Deep Processing, Jilin Agricultural University, Changchun 130118, Jilin, P. R. China
| | - Shuo Sun
- College of Food Science and Engineering, National Engineering Laboratory of Wheat and Corn Deep Processing, Jilin Agricultural University, Changchun 130118, Jilin, P. R. China
| | - Jiale Liu
- College of Food Science and Engineering, National Engineering Laboratory of Wheat and Corn Deep Processing, Jilin Agricultural University, Changchun 130118, Jilin, P. R. China
| | - Qiao Dang
- College of Food Science and Engineering, National Engineering Laboratory of Wheat and Corn Deep Processing, Jilin Agricultural University, Changchun 130118, Jilin, P. R. China
| | - Yawen Gao
- College of Food Science and Engineering, National Engineering Laboratory of Wheat and Corn Deep Processing, Jilin Agricultural University, Changchun 130118, Jilin, P. R. China
| | - Li Fang
- College of Food Science and Engineering, National Engineering Laboratory of Wheat and Corn Deep Processing, Jilin Agricultural University, Changchun 130118, Jilin, P. R. China
| | - Weihong Min
- College of Food Science and Engineering, National Engineering Laboratory of Wheat and Corn Deep Processing, Jilin Agricultural University, Changchun 130118, Jilin, P. R. China
- College of Food and Health, Zhejiang A&F University, Hangzhou 311300, P. R. China
| |
Collapse
|
3
|
Zhao W, Yang H, Cui H, Li W, Xing S, Han W. Elucidating the structural basis of vitamin B 12 derivatives as novel potent inhibitors of PTP1B: Insights from inhibitory mechanisms using Gaussian accelerated molecular dynamics (GaMD) and in vitro study. Int J Biol Macromol 2024; 268:131902. [PMID: 38692532 DOI: 10.1016/j.ijbiomac.2024.131902] [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: 11/01/2023] [Revised: 04/19/2024] [Accepted: 04/24/2024] [Indexed: 05/03/2024]
Abstract
Vitamin B12 is a group of biologically active cobalamin compounds. In this study, we investigated the inhibitory effects of methylcobalamin (MeCbl) and hydroxocobalamin acetate (OHCbl Acetate) on protein tyrosine phosphatase 1B (PTP1B). MeCbl and OHCbl Acetate exhibited an IC50 of approximately 58.390 ± 2.811 μM and 8.998 ± 0.587 μM, respectively. The Ki values of MeCbl and OHCbl Acetate were 25.01 μM and 4.04 μM respectively. To elucidate the inhibition mechanism, we conducted a 500 ns Gaussian accelerated molecular dynamics (GaMD) simulation. Utilizing PCA and tICA, we constructed Markov state models (MSM) to examine secondary structure changes during motion. Our findings revealed that the α-helix at residues 37-42 remained the most stable in the PTP1B-OHCbl Acetate system. Furthermore, upon binding of OHCbl Acetate or MeCbl, the WPD loop of PTP1B moved inward to the active pocket, forming a closed conformation and potentially obstructs substrate entry. Protein-ligand interaction analysis and MM-PBSA showed that OHCbl Acetate exhibited lower binding free energy and engaged in more residue interactions with PTP1B. In summary, our study confirmed the substantial inhibitory activity of OHCbl Acetate against PTP1B, with its inhibitory potency notably surpassing that of MeCbl. We demonstrated potential molecular mechanisms of OHCbl Acetate inhibiting PTP1B.
Collapse
Affiliation(s)
- Wencheng Zhao
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Hengzheng Yang
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Huizi Cui
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Wannan Li
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, 2699 Qianjin Street, Changchun 130012, China.
| | - Shu Xing
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, 2699 Qianjin Street, Changchun 130012, China.
| | - Weiwei Han
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, 2699 Qianjin Street, Changchun 130012, China.
| |
Collapse
|
4
|
Rabaan AA, Alshahrani FS, Garout M, Alissa M, Mashraqi MM, Alshehri AA, Alsaleh AA, Alwarthan S, Sabour AA, Alfaraj AH, AlShehail BM, Alotaibi N, Abduljabbar WA, Aljeldah M, Alestad JH. Repositioning of anti-infective compounds against monkeypox virus core cysteine proteinase: a molecular dynamics study. Mol Divers 2024:10.1007/s11030-023-10802-8. [PMID: 38652365 DOI: 10.1007/s11030-023-10802-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2023] [Accepted: 12/26/2023] [Indexed: 04/25/2024]
Abstract
Monkeypox virus (MPXV) core cysteine proteinase (CCP) is one of the major drug targets used to examine the inhibitory action of chemical moieties. In this study, an in silico technique was applied to screen 1395 anti-infective compounds to find out the potential molecules against the MPXV-CCP. The top five hits were selected after screening and processed for exhaustive docking based on the docked score of ≤ -9.5 kcal/mol. Later, the top three hits based on the exhaustive-docking score and interaction profile were selected to perform MD simulations. The overall RMSD suggested that two compounds, SC75741 and ammonium glycyrrhizinate, showed a highly stable complex with a standard deviation of 0.18 and 0.23 nm, respectively. Later, the MM/GBSA binding free energies of complexes showed significant binding strength with ΔGTOTAL from -21.59 to -15 kcal/mol. This report reported the potential inhibitory activity of SC75741 and ammonium glycyrrhizinate against MPXV-CCP by competitively inhibiting the binding of the native substrate.
Collapse
Affiliation(s)
- Ali A Rabaan
- Molecular Diagnostic Laboratory, Johns Hopkins Aramco Healthcare, 31311, Dhahran, Saudi Arabia.
- College of Medicine, Alfaisal University, 11533, Riyadh, Saudi Arabia.
- Department of Public Health and Nutrition, The University of Haripur, Haripur, 22610, Pakistan.
| | - Fatimah S Alshahrani
- Department of Internal Medicine, College of Medicine, King Saud University, 11362, Riyadh, Saudi Arabia
- Division of Infectious Diseases, Department of Internal Medicine, College of Medicine, King Saud University and King Saud University Medical City, 11451, Riyadh, Saudi Arabia
| | - Mohammed Garout
- Department of Community Medicine and Health Care for Pilgrims, Faculty of Medicine, Umm Al-Qura University, 21955, Makkah, Saudi Arabia
| | - Mohammed Alissa
- Department of Medical Laboratory, College of Applied Medical Sciences, Prince Sattam bin Abdulaziz University, 11942, Al-Kharj, Saudi Arabia
| | - Mutaib M Mashraqi
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Najran University, 61441, Najra, Saudi Arabia
| | - Ahmad A Alshehri
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Najran University, 61441, Najra, Saudi Arabia
| | - Abdulmonem A Alsaleh
- Clinical Laboratory Science Department, Mohammed Al-Mana College for Medical Sciences, 34222, Dammam, Saudi Arabia
| | - Sara Alwarthan
- Department of Internal Medicine, College of Medicine, Imam Abdulrahman Bin Faisal University, 34212, Dammam, Saudi Arabia
| | - Amal A Sabour
- Department of Botany and Microbiology, College of Science, King Saud University, 11451, Riyadh, Saudi Arabia
| | - Amal H Alfaraj
- Pediatric Department, Abqaiq General Hospital, First Eastern Health Cluster, 33261, Abqaiq, Saudi Arabia
| | - Bashayer M AlShehail
- Pharmacy Practice Department, College of Clinical Pharmacy, Imam Abdulrahman Bin Faisal University, 31441, Dammam, Saudi Arabia
| | - Nouf Alotaibi
- Clinical pharmacy Department, College of Pharmacy, Umm Al-Qura University, 21955, Makkah, Saudi Arabia
| | - Wesam A Abduljabbar
- Department of Medical laboratory sciences, Fakeeh College for Medical Science, 21134, Jeddah, Saudi Arabia
| | - Mohammed Aljeldah
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, University of Hafr Al Batin, 39831, Hafr Al Batin, Saudi Arabia
| | - Jeehan H Alestad
- Immunology and Infectious Microbiology Department, University of Glasgow, Glasgow, G1 1XQ, UK.
- Microbiology Department, Collage of Medicine, 46300, Jabriya, Kuwait.
| |
Collapse
|
5
|
Wu Z, Wang C, Li C, Xu N, Cao X, Chen S, Shi Y, He Y, Zhang P, Ji J. Integrated Computational Pipeline for the High-Throughput Discovery of Cell Adhesion Peptides. J Phys Chem Lett 2024; 15:3748-3756. [PMID: 38551401 DOI: 10.1021/acs.jpclett.4c00393] [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: 04/12/2024]
Abstract
Cell adhesion peptides (CAPs) often play a critical role in tissue engineering research. However, the discovery of novel CAPs for diverse applications remains a challenging and time-intensive process. This study presents an efficient computational pipeline integrating sequence embeddings, binding predictors, and molecular dynamics simulations to expedite the discovery of new CAPs. A Pro2vec model, trained on vast CAP data sets, was built to identify RGD-similar tripeptide candidates. These candidates were further evaluated for their binding affinity with integrin receptors using the Mutabind2 machine learning model. Additionally, molecular dynamics simulations were applied to model receptor-peptide interactions and calculate their binding free energies, providing a quantitative assessment of the binding strength for further screening. The resulting peptide demonstrated performance comparable to that of RGD in endothelial cell adhesion and spreading experimental assays, validating the efficacy of the integrated computational pipeline.
Collapse
Affiliation(s)
- Zhiyu Wu
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310058, China
- Institute of Zhejiang University-Quzhou, Quzhou 324000, China
| | - Cong Wang
- MOE Key Laboratory of Macromolecular Synthesis and Functionalization, Department of Polymer Science and Engineering, Zhejiang University, Hangzhou 310058, China
| | - Chen Li
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310058, China
- Institute of Zhejiang University-Quzhou, Quzhou 324000, China
| | - Nan Xu
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310058, China
- Institute of Zhejiang University-Quzhou, Quzhou 324000, China
| | - Xiaoyong Cao
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310058, China
- Institute of Zhejiang University-Quzhou, Quzhou 324000, China
| | - Shengfu Chen
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310058, China
| | - Yao Shi
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310058, China
- Key Laboratory of Biomass Chemical Engineering of Ministry of Education, Zhejiang University, Hangzhou 310058, China
| | - Yi He
- College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310058, China
- Institute of Zhejiang University-Quzhou, Quzhou 324000, China
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Peng Zhang
- MOE Key Laboratory of Macromolecular Synthesis and Functionalization, Department of Polymer Science and Engineering, Zhejiang University, Hangzhou 310058, China
- State Key Laboratory of Transvascular Implantation Devices, Qidi Road 456, Hangzhou 310058, China
| | - Jian Ji
- MOE Key Laboratory of Macromolecular Synthesis and Functionalization, Department of Polymer Science and Engineering, Zhejiang University, Hangzhou 310058, China
- State Key Laboratory of Transvascular Implantation Devices, Qidi Road 456, Hangzhou 310058, China
| |
Collapse
|
6
|
Demian MD, Amasiorah VI, Johnson TO, Ebenyi LN. Phytochemical identification and in silico elucidation of interactions of bioactive compounds from Citrullus lanatus with androgen receptor towards prostate cancer treatment. In Silico Pharmacol 2024; 12:27. [PMID: 38596366 PMCID: PMC10999405 DOI: 10.1007/s40203-024-00193-5] [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: 05/22/2023] [Accepted: 02/07/2024] [Indexed: 04/11/2024] Open
Abstract
Androgen receptor (AR) is known to play a crucial role in the development and progression of prostate cancer, and compounds that inhibit its activity are regarded as promising for the development of drugs to treat the disease. This study aimed to investigate the AR-inhibiting potential of Citrullus lanatus fruit compounds for prostate cancer drug development. Following HPLC identification, the binding energies, molecular interactions, and pharmacological potentials of the compounds against AR were elucidated using in silico techniques such as, molecular docking, induced-fit docking, molecular dynamics simulation, and ADMET prediction. Some of the compounds found to be present in Citrullus lanatus fruit included flavonoids such as proanthocyanin, naringin, flavan 3 ol, flavonones, naringenin, epicatechin, citrulline, and catechin. Naringenin exhibited the highest docking score in the molecular docking analysis, followed by resveratrol, ribalinidine, and epicatechin. These compounds share a common AR binding site with the standard ligand, dihydrotestosterone (DHT). Some of the compounds showed favorable ADMET profiles, while others showed at least one toxicity potential. The induced-fit docking of naringenin with AR yielded a higher docking score than the initial score obtained from standard docking while preserving stable molecular contacts with the interacting amino acids. Consistent hydrogen bond interactions of naringenin with PHE 764, ASN 705, and THR 877 of AR, including a persistent pi-pi stacking contact with PHE 764, were observed from the molecular dynamic simulation. The Citrullus lanatus compounds, particularly naringenin, may therefore be considered for further research towards the development of drugs for prostate cancer therapy.
Collapse
Affiliation(s)
| | | | - Titilayo Omolara Johnson
- Department of Biochemistry, Faculty of Basic Medical Science, College of Health Sciences, University of Jos, Jos, Nigeria
| | - Lilian N. Ebenyi
- Department of Biotechnology, Ebonyi State University, Abakaliki, Ebonyi State Nigeria
| |
Collapse
|
7
|
Wang M, Liu K. Molecular dynamics simulations to explore the binding mode between the amyloid-β protein precursor (APP) and adaptor protein Mint2. Sci Rep 2024; 14:7975. [PMID: 38575686 PMCID: PMC10995209 DOI: 10.1038/s41598-024-58584-9] [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: 12/04/2023] [Accepted: 04/01/2024] [Indexed: 04/06/2024] Open
Abstract
Alzheimer's disease (AD) presents a significant challenge in neurodegenerative disease management, with limited therapeutic options available for its prevention and treatment. At the heart of AD pathogenesis is the amyloid-β (Aβ) protein precursor (APP), with the interaction between APP and the adaptor protein Mint2 being crucial. Despite previous explorations into the APP-Mint2 interaction, the dynamic regulatory mechanisms by which Mint2 modulates APP binding remain poorly understood. This study undertakes molecular dynamics simulations across four distinct systems-free Mint2, Mint2 bound to APP, a mutant form of Mint2, and the mutant form bound to APP-over an extensive 400 ns timeframe. Our findings reveal that the mutant Mint2 experiences significant secondary structural transformations, notably the formation of an α-helix in residues S55-K65 upon APP binding, within the 400 ns simulation period. Additionally, we observed a reduction in the active pocket size of the mutant Mint2 compared to its wild-type counterpart, enhancing its APP binding affinity. These insights hold promise for guiding the development of novel inhibitors targeting the Mints family, potentially paving the way for new therapeutic strategies in AD prevention and treatment.
Collapse
Affiliation(s)
- Min Wang
- International Research Centre for Nano Handling and Manufacturing of China, Changchun University of Science and Technology, Changchun, 130022, China.
- Ministry of Education Key Laboratory for Cross-Scale Micro and Nano Manufacturing, Changchun University of Science and Technology, Changchun, 130022, China.
| | - Kaifeng Liu
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Changchun, 130012, China
| |
Collapse
|
8
|
Jiang D, Du H, Zhao H, Deng Y, Wu Z, Wang J, Zeng Y, Zhang H, Wang X, Wang E, Hou T, Hsieh CY. Assessing the performance of MM/PBSA and MM/GBSA methods. 10. Prediction reliability of binding affinities and binding poses for RNA-ligand complexes. Phys Chem Chem Phys 2024; 26:10323-10335. [PMID: 38501198 DOI: 10.1039/d3cp04366e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Ribonucleic acid (RNA)-ligand interactions play a pivotal role in a wide spectrum of biological processes, ranging from protein biosynthesis to cellular reproduction. This recognition has prompted the broader acceptance of RNA as a viable candidate for drug targets. Delving into the atomic-scale understanding of RNA-ligand interactions holds paramount importance in unraveling intricate molecular mechanisms and further contributing to RNA-based drug discovery. Computational approaches, particularly molecular docking, offer an efficient way of predicting the interactions between RNA and small molecules. However, the accuracy and reliability of these predictions heavily depend on the performance of scoring functions (SFs). In contrast to the majority of SFs used in RNA-ligand docking, the end-point binding free energy calculation methods, such as molecular mechanics/generalized Born surface area (MM/GBSA) and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA), stand as theoretically more rigorous approaches. Yet, the evaluation of their effectiveness in predicting both binding affinities and binding poses within RNA-ligand systems remains unexplored. This study first reported the performance of MM/PBSA and MM/GBSA with diverse solvation models, interior dielectric constants (εin) and force fields in the context of binding affinity prediction for 29 RNA-ligand complexes. MM/GBSA is based on short (5 ns) molecular dynamics (MD) simulations in an explicit solvent with the YIL force field; the GBGBn2 model with higher interior dielectric constant (εin = 12, 16 or 20) yields the best correlation (Rp = -0.513), which outperforms the best correlation (Rp = -0.317, rDock) offered by various docking programs. Then, the efficacy of MM/GBSA in identifying the near-native binding poses from the decoys was assessed based on 56 RNA-ligand complexes. However, it is evident that MM/GBSA has limitations in accurately predicting binding poses for RNA-ligand systems, particularly compared with notably proficient docking programs like rDock and PLANTS. The best top-1 success rate achieved by MM/GBSA rescoring is 39.3%, which falls below the best results given by docking programs (50%, PLNATS). This study represents the first evaluation of MM/PBSA and MM/GBSA for RNA-ligand systems and is expected to provide valuable insights into their successful application to RNA targets.
Collapse
Affiliation(s)
- Dejun Jiang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
- Hangzhou Carbonsilicon AI Technology Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Hongyan Du
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Huifeng Zhao
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
- Hangzhou Carbonsilicon AI Technology Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Yafeng Deng
- Hangzhou Carbonsilicon AI Technology Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Zhenxing Wu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Jike Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Yundian Zeng
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Haotian Zhang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Xiaorui Wang
- China State Key Laboratory of Quality Research in Chinese Medicines, Macau University of Science and Technology, Macau 999078, China
| | - Ercheng Wang
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, China.
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Chang-Yu Hsieh
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| |
Collapse
|
9
|
Liao J, Wu M, Gao J, Chen C. Calculation of solvation force in molecular dynamics simulation by deep-learning method. Biophys J 2024:S0006-3495(24)00163-2. [PMID: 38444159 DOI: 10.1016/j.bpj.2024.02.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Revised: 01/31/2024] [Accepted: 02/29/2024] [Indexed: 03/07/2024] Open
Abstract
Electrostatic calculations are generally used in studying the thermodynamics and kinetics of biomolecules in solvent. Generally, this is performed by solving the Poisson-Boltzmann equation on a large grid system, a process known to be time consuming. In this study, we developed a deep neural network to predict the decomposed solvation free energies and forces of all atoms in a molecule. To train the network, the internal coordinates of the molecule were used as the input data, and the solvation free energies along with transformed atomic forces from the Poisson-Boltzmann equation were used as labels. Both the training and prediction tasks were accelerated on GPU. Formal tests demonstrated that our method can provide reasonable predictions for small molecules when the network is well-trained with its simulation data. This method is suitable for processing lots of snapshots of molecules in a long trajectory. Moreover, we applied this method in the molecular dynamics simulation with enhanced sampling. The calculated free energy landscape closely resembled that obtained from explicit solvent simulations.
Collapse
Affiliation(s)
- Jun Liao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Mincong Wu
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Junyong Gao
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China.
| |
Collapse
|
10
|
Barbhuiya T, Beard S, Shah ET, Mason S, Bolderson E, O’Byrne K, Guddat LW, Richard DJ, Adams MN, Gandhi NS. Targeting the hSSB1-INTS3 Interface: A Computational Screening Driven Approach to Identify Potential Modulators. ACS OMEGA 2024; 9:8362-8373. [PMID: 38405517 PMCID: PMC10882649 DOI: 10.1021/acsomega.3c09267] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Revised: 01/20/2024] [Accepted: 01/25/2024] [Indexed: 02/27/2024]
Abstract
Human single-stranded DNA binding protein 1 (hSSB1) forms a heterotrimeric complex, known as a sensor of single-stranded DNA binding protein 1 (SOSS1), in conjunction with integrator complex subunit 3 (INTS3) and C9ORF80. This sensory protein plays an important role in homologous recombination repair of double-strand breaks in DNA to efficiently recruit other repair proteins at the damaged sites. Previous studies have identified elevated hSSB1-mediated DNA repair activities in various cancers, highlighting its potential as an anticancer target. While prior efforts have focused on inhibiting hSSB1 by targeting its DNA binding domain, this study seeks to explore the inhibition of the hSSB1 function by disrupting its interaction with the key partner protein INTS3 in the SOSS1 complex. The investigative strategy entails a molecular docking-based screening of a specific compound library against the three-dimensional structure of INTS3 at the hSSB1 binding interface. Subsequent assessments involve in vitro analyses of protein-protein interaction (PPI) disruption and cellular effects through co-immunoprecipitation and immunofluorescence assays, respectively. Moreover, the study includes an evaluation of the structural stability of ligands at the INTS3 hot-spot site using molecular dynamics simulations. The results indicate a potential in vitro disruption of the INTS3-hSSB1 interaction by three of the tested compounds obtained from the virtual screening with one impacting the recruitment of hSSB1 and INTS3 to chromatin following DNA damage. To our knowledge, our results identify the first set of drug-like compounds that functionally target INTS3-hSSB1 interaction, and this provides the basis for further biophysical investigations that should help to speed up PPI inhibitor discovery.
Collapse
Affiliation(s)
- Tabassum
Khair Barbhuiya
- Centre
for Genomics and Personalised Health, and School of Chemistry and
Physics, Faculty of Science, Queensland
University of Technology, 2 George Street, Brisbane, QLD 4000, Australia
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
| | - Sam Beard
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
- Centre
for Genomics and Personalised Health, and School of Biomedical Sciences,
Faculty of Health, Queensland University
of Technology, Kelvin Grove, QLD 4059, Australia
| | - Esha T. Shah
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
- Centre
for Genomics and Personalised Health, and School of Biomedical Sciences,
Faculty of Health, Queensland University
of Technology, Kelvin Grove, QLD 4059, Australia
| | - Steven Mason
- School
of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, QLD 4072, Australia
| | - Emma Bolderson
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
- Centre
for Genomics and Personalised Health, and School of Biomedical Sciences,
Faculty of Health, Queensland University
of Technology, Kelvin Grove, QLD 4059, Australia
| | - Ken O’Byrne
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
- Centre
for Genomics and Personalised Health, and School of Biomedical Sciences,
Faculty of Health, Queensland University
of Technology, Kelvin Grove, QLD 4059, Australia
| | - Luke W. Guddat
- School
of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, QLD 4072, Australia
| | - Derek J. Richard
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
- Centre
for Genomics and Personalised Health, and School of Biomedical Sciences,
Faculty of Health, Queensland University
of Technology, Kelvin Grove, QLD 4059, Australia
| | - Mark N. Adams
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
- Centre
for Genomics and Personalised Health, and School of Biomedical Sciences,
Faculty of Health, Queensland University
of Technology, Kelvin Grove, QLD 4059, Australia
| | - Neha S. Gandhi
- Centre
for Genomics and Personalised Health, and School of Chemistry and
Physics, Faculty of Science, Queensland
University of Technology, 2 George Street, Brisbane, QLD 4000, Australia
- Cancer
and Ageing Research Program, Woolloongabba, QLD 4102, Australia
| |
Collapse
|
11
|
Zhang X, Tong J, Wang T, Wang T, Xu L, Wang Z, Hou T, Pan P. Dissecting the role of ALK double mutations in drug resistance to lorlatinib with in-depth theoretical modeling and analysis. Comput Biol Med 2024; 169:107815. [PMID: 38128254 DOI: 10.1016/j.compbiomed.2023.107815] [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: 09/14/2023] [Revised: 11/03/2023] [Accepted: 12/03/2023] [Indexed: 12/23/2023]
Abstract
Anaplastic lymphoma kinase (ALK) is implicated in the genesis of multiple malignant tumors. Lorlatinib stands out as the most advanced and effective inhibitor currently used in the clinic for the treatment of ALK-positive non-small cell lung cancer. However, resistance to lorlatinib has inevitably manifested over time, with double/triple mutations of G1202, L1196, L1198, C1156 and I1171 frequently observed in clinical practice, and tumors regrow within a short time after treatment with lorlatinib. Therefore, elucidating the mechanism of resistance to lorlatinib is paramount in paving the way for innovative therapeutic strategies and the development of next-generation drugs. In this study, we leveraged multiple computational methodologies to delve into the resistance mechanisms of three specific double mutations of ALKG1202R/L1196M, ALKG1202R/L1198F and ALKI1171N/L1198F to lorlatinib. We analyzed these mechanisms through qualitative (PCA, DCCM) and quantitative (MM/GBSA, US) kinetic analyses. The qualitative analysis shows that these mutations exert minimal perturbations on the conformational dynamics of the structural domains of ALK. The energetic and structural assessments show that the van der Waals interactions, formed by the conserved residue Leu1256 within the ATP-binding site and the residues Glu1197 and Met1199 in the hinge domain with lorlatinib, play integral roles in the occurrence of drug resistance. Furthermore, the US simulation results elucidate that the pathways through which lorlatinib dissociates vary across mutant systems, and the distinct environments during the dissociation process culminate in diverse resistance mechanisms. Collectively, these insights provide important clues for the design of novel inhibitors to combat resistance.
Collapse
Affiliation(s)
- Xing Zhang
- College of Chemistry and Chemical Engineering, Shaanxi University of Science and Technology, Xi'an, 710021, Shaanxi, China; College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| | - Jianbo Tong
- College of Chemistry and Chemical Engineering, Shaanxi University of Science and Technology, Xi'an, 710021, Shaanxi, China.
| | - Tianhao Wang
- College of Chemistry and Chemical Engineering, Shaanxi University of Science and Technology, Xi'an, 710021, Shaanxi, China; College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| | - Tianyue Wang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| | - Lei Xu
- Institute of Bioinformatics and Medical Engineering, School of Electrical and Information Engineering, Jiangsu University of Technology, Changzhou, 213001, China
| | - Zhe Wang
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China
| | - Peichen Pan
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, 310058, Zhejiang, China.
| |
Collapse
|
12
|
Fan M, Zhou Z, Zhu W, Li M, Tu Y, Yu Z, Li J, Zhang M, Liang K. Reinforced dentin remineralization via a novel dual-affinity peptide. Dent Mater 2024; 40:254-266. [PMID: 37989605 DOI: 10.1016/j.dental.2023.11.009] [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: 07/25/2023] [Revised: 10/30/2023] [Accepted: 11/04/2023] [Indexed: 11/23/2023]
Abstract
OBJECTIVES In light of the constantly flowing saliva, anti-caries remineralization agents are inclined to be taken away. Owing to their limited residence time, the remineralization effect is not as desirable as expected. Hence, our study aimed to synthesize a novel peptide (DGP) with high affinity to both collagen fibrils and hydroxyapatite, and investigated its dentin remineralization efficacy in vitro and anti-caries capability in vivo. METHODS DGP was synthesized through Fmoc solid-phase reaction. The binding ability and interaction mechanism of DGP to demineralized dentin were investigated. Dentin specimens were demineralized, then treated with DGP and deionized water respectively. The specimens were incubated in artificial saliva and in-vitro remineralization effectiveness was analyzed after 14 days. The rat caries model was established to further scrutinize the in-vivo efficacy of caries prevention. RESULTS DGP possesses an enhanced adhesion force of 12.29 ± 1.12 nN to demineralized dentin. The favorable adsorption capacity is ascribed to the stable hydrogen bonds between S2P-101 and ASP-100 of DGP and GLY33 and PRO-16 of collagen fibers. Abundant mineral deposits and remarkable tubule occlusion were observed in the DGP group. DGP-treated dentin obtained notable microhardness recovery and higher mineral content after a 14-day remineralization regimen. DGP also demonstrated potent caries prevention in vivo, with substantially fewer carious lesions and significantly lower Keyes scoring. SIGNIFICANCE DGP proves to possess a high affinity to demineralized dentin regardless of saliva flowing, thus enhancing remineralization potency significantly in vitro and in vivo, potential for dental caries prevention and combatting initial dentin caries clinically.
Collapse
Affiliation(s)
- Menglin Fan
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China
| | - Zilin Zhou
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China
| | - Wanchun Zhu
- Chongqing Medical University Stomatology College, Chongqing 401147, China
| | - Meng Li
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China
| | - Yuanyuan Tu
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China
| | - Zhaohan Yu
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China
| | - Jiyao Li
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China
| | - Min Zhang
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China.
| | - Kunneng Liang
- State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, Department of Cariology and Endodontics, West China Hospital of Stomatology, Sichuan University, Chengdu 610041, China.
| |
Collapse
|
13
|
Rabaan AA, AlShehail BM, Halwani MA, Alshengeti A, Najim MA, Garout M, Bajunaid HA, Alshamrani SA, Al Fares MA, Alissa M, Alwashmi ASS. Investigation of Zika virus methyl transferase inhibitors using steered molecular dynamics. J Biomol Struct Dyn 2024; 42:1711-1724. [PMID: 37325855 DOI: 10.1080/07391102.2023.2224882] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Accepted: 04/05/2023] [Indexed: 06/17/2023]
Abstract
Zika virus (ZIKV) spread is considered a major public health threat by the World Health Organization (WHO). There are no vaccines or drugs available to control the infection of the Zika virus, therefore a highly effective medicinal molecule is urgently required. In this study, a computationally intensive investigation was performed to identify a potent natural compound that could inhibit the ZIKV NS5 methyltransferase. This research approach is based on target-based drug identification principles where the native inhibitor SAH (S-adenosylhomocysteine) of ZIKV NS5 methyltransferase was selected as a reference. High-throughput virtual screening and tanimoto similarity coefficient were applied to the natural compound library for ranking the potential candidates. The top five compounds were selected for interaction analysis, MD simulation, total binding free energy through MM/GBSA, and steered MD simulation. Among these compounds, Adenosine 5'-monophosphate monohydrate, Tubercidin, and 5-Iodotubercidin showed stable binding to the protein compared to the native compound, SAH. These three compounds also showed less fluctuations in RMSF in contrast to native compound. Additionally, the same interacting residues observed in SAH also made strong interactions with these three compounds. Adenosine 5'-monophosphate monohydrate and 5-Iodotubercidin had greater total binding free energies than the reference ligand. Moreover, the dissociation resistance of all three compounds was equivalent to that of the reference ligand. This study suggested binding properties of three-hit compounds that could be used to develop drugs against Zika virus infections.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Ali A Rabaan
- Molecular Diagnostic Laboratory, Johns Hopkins Aramco Healthcare, Dhahran, Saudi Arabia
- College of Medicine, Alfaisal University, Riyadh, Saudi Arabia
- Department of Public Health and Nutrition, The University of Haripur, Haripur, Pakistan
| | - Bashayer M AlShehail
- Pharmacy Practice Department, College of Clinical Pharmacy, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
| | - Muhammad A Halwani
- Department of Medical Microbiology, Faculty of Medicine, Al Baha University, Saudi Arabia
| | - Amer Alshengeti
- Department of Pediatrics, College of Medicine, Taibah University, Al-Madinah, Saudi Arabia
- Department of Infection Prevention and Control, Prince Mohammad Bin Abdulaziz Hospital, National Guard Health Affairs, Al-Madinah, Saudi Arabia
| | - Mustafa A Najim
- Department of Medical Laboratories Technology, College of Applied Medical Sciences, Taibah University, Madinah, Saudi Arabia
| | - Mohammed Garout
- Department of Community Medicine and Health Care for Pilgrims, Faculty of Medicine, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Huda A Bajunaid
- Makkah Specialized Laboratory, Fakeeh Care group, Hadda, Saudi Arabia
| | - Saleh A Alshamrani
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Najran University, Najran, Saudi Arabia
| | - Mona A Al Fares
- Department of Internal Medicine, King Abdulaziz University Hospital, Jeddah, Saudi Arabia
| | - Mohammed Alissa
- Department of Medical Laboratory Sciences, College of Applied Medical Sciences, Prince Sattam bin Abdulaziz University, Al-Kharj, Saudi Arabia
| | - Ameen S S Alwashmi
- Department of Medical Laboratories, College of Applied Medical Sciences, Qassim University, Buraydah, Saudi Arabia
| |
Collapse
|
14
|
Quesnel A, Martin LD, Tarzi C, Lenis VP, Coles N, Islam M, Angione C, Outeiro TF, Khundakar AA, Filippou PS. Uncovering potential diagnostic and pathophysiological roles of α-synuclein and DJ-1 in melanoma. Cancer Med 2024; 13. [PMID: 38189631 PMCID: PMC10807602 DOI: 10.1002/cam4.6900] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 11/20/2023] [Accepted: 12/19/2023] [Indexed: 01/09/2024] Open
Abstract
BACKGROUND Melanoma, the most lethal skin cancer type, occurs more frequently in Parkinson's disease (PD), and PD is more frequent in melanoma patients, suggesting disease mechanisms overlap. α-synuclein, a protein that accumulates in PD brain, and the oncogene DJ-1, which is associated with PD autosomal recessive forms, are both elevated in melanoma cells. Whether this indicates melanoma progression or constitutes a protective response remains unclear. We hereby investigated the molecular mechanisms through which α-synuclein and DJ-1 interact, suggesting novel biomarkers and targets in melanoma. METHODS The Cancer Genome Atlas (TCGA) expression profiles derived from UCSC Xena were used to obtain α-synuclein and DJ-1 expression and correlated with survival in skin cutaneous melanoma (SKCM). Immunohistochemistry determined the expression in metastatic melanoma lymph nodes. Protein-protein interactions (PPIs) and molecular docking assessed protein binding and affinity with chemotherapeutic drugs. Further validation was performed using in vitro cellular models and ELISA immunoassays. RESULTS α-synuclein and DJ-1 were upregulated in primary and metastatic SKCM. Aggregated α-synuclein was selectively detected in metastatic melanoma lymph nodes. α-synuclein overexpression in SK-MEL-28 cells induced the expression of DJ-1, supporting PPI and a positive correlation in melanoma patients. Molecular docking revealed a stable protein complex, with differential binding to chemotherapy drugs such as temozolomide, dacarbazine, and doxorubicin. Parallel reduction of both proteins in temozolomide-treated SK-MEL-28 spheroids suggests drug binding may affect protein interaction and/or stability. CONCLUSION α-synuclein, together with DJ-1, may play a role in melanoma progression and chemosensitivity, constituting novel targets for therapeutic intervention, and possible biomarkers for melanoma.
Collapse
Affiliation(s)
- Agathe Quesnel
- School of Health & Life Sciences, Teesside University, Middlesbrough, UK
- National Horizons Centre, Teesside University, Darlington, UK
| | - Leya Danielle Martin
- School of Health & Life Sciences, Teesside University, Middlesbrough, UK
- National Horizons Centre, Teesside University, Darlington, UK
| | - Chaimaa Tarzi
- School of Computing, Engineering & Digital Technologies, Teesside University, Middlesbrough, UK
- Centre for Digital Innovation, Teesside University, Middlesbrough, UK
| | - Vasileios P Lenis
- School of Health & Life Sciences, Teesside University, Middlesbrough, UK
- National Horizons Centre, Teesside University, Darlington, UK
| | - Nathan Coles
- School of Health & Life Sciences, Teesside University, Middlesbrough, UK
- National Horizons Centre, Teesside University, Darlington, UK
| | - Meez Islam
- School of Health & Life Sciences, Teesside University, Middlesbrough, UK
- National Horizons Centre, Teesside University, Darlington, UK
| | - Claudio Angione
- National Horizons Centre, Teesside University, Darlington, UK
- School of Computing, Engineering & Digital Technologies, Teesside University, Middlesbrough, UK
- Centre for Digital Innovation, Teesside University, Middlesbrough, UK
| | - Tiago F Outeiro
- Translational and Clinical Research Institute, Faculty of Medical Sciences, Newcastle University, Newcastle upon Tyne, UK
- Department of Experimental Neurodegeneration, Center for Biostructural Imaging of Neurodegeneration, University Medical Center, Göttingen, Germany
- Max Planck Institute for Multidisciplinary Sciences, Göttingen, Germany
- Deutsches Zentrum für Neurodegenerative Erkrankungen (DZNE), Göttingen, Germany
| | - Ahmad A Khundakar
- School of Health & Life Sciences, Teesside University, Middlesbrough, UK
- National Horizons Centre, Teesside University, Darlington, UK
- Translational and Clinical Research Institute, Faculty of Medical Sciences, Newcastle University, Newcastle upon Tyne, UK
| | - Panagiota S Filippou
- School of Health & Life Sciences, Teesside University, Middlesbrough, UK
- National Horizons Centre, Teesside University, Darlington, UK
| |
Collapse
|
15
|
Valdivia A, Luque FJ, Llabrés S. Binding of Cholesterol to the N-Terminal Domain of the NPC1L1 Transporter: Analysis of the Epimerization-Related Binding Selectivity and Loop Mutations. J Chem Inf Model 2024; 64:189-204. [PMID: 38152929 PMCID: PMC10777396 DOI: 10.1021/acs.jcim.3c01319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 12/14/2023] [Accepted: 12/15/2023] [Indexed: 12/29/2023]
Abstract
Cholesterol is a fat-like substance with a pivotal physiological relevance in humans, and its homeostasis is tightly regulated by various cellular processes, including the import in the small intestine and the reabsorption in the biliary ducts by the Niemann-Pick C1 Like 1 (NPC1L1) importer. NPC1L1 can mediate the absorption of a variety of sterols but strikingly exhibits a large sensitivity to cholesterol epimerization. This study examines the molecular basis of the epimerization-related selective binding of cholesterol by combining extended unbiased molecular dynamics simulations of the apo and holo species of the N-terminal domain of wild-type NPC1L1, in conjunction with relative binding free energy, umbrella sampling, and well-tempered metadynamics calculations. The analysis of the results discloses the existence of two distinct binding modes for cholesterol and epi-cholesterol. The former binds deeper in the cavity, forming key hydrogen-bond interactions with Q95, S56, and a water molecule. In contrast, epi-cholesterol is shifted ca. 3 Å to the mouth of the cavity and the transition to the Q95 site is prevented by an energetic barrier of 4.1 kcal·mol-1. Thus, the configuration of the hydroxyl group of cholesterol, together with the presence of a structural water molecule, is a key feature for effective absorption. Finally, whereas these findings may seemingly be challenged by single-point mutations that impair cholesterol transport but have a mild impact on the binding of cholesterol to the Q95 binding site, our results reveal that they have a drastic influence on the conformational landscape of the α8/β7 loop in the apo species compared to the wild-type protein. Overall, the results give support to the functional role played by the α8/β7 loop in regulating the access of ligands to NPC1L1, and hence to interpreting the impact of these mutations on diseases related to disruption of sterol absorption, paving the way to understanding certain physiological dysfunctions.
Collapse
Affiliation(s)
- Aitor Valdivia
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia
i Ciències de l′Alimentació—Campus Torribera,
Universitat de Barcelona, Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
- Institut
de Biomedicina (IBUB), Universitat de Barcelona, 08028 Barcelona, Spain
| | - F. Javier Luque
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia
i Ciències de l′Alimentació—Campus Torribera,
Universitat de Barcelona, Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
- Institut
de Biomedicina (IBUB), Universitat de Barcelona, 08028 Barcelona, Spain
- Institut
de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08921 Barcelona, Spain
| | - Salomé Llabrés
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia
i Ciències de l′Alimentació—Campus Torribera,
Universitat de Barcelona, Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
- Institut
de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08921 Barcelona, Spain
| |
Collapse
|
16
|
He Y, Liu K, Cao F, Song R, Liu J, Zhang Y, Li W, Han W. Using deep learning and molecular dynamics simulations to unravel the regulation mechanism of peptides as noncompetitive inhibitor of xanthine oxidase. Sci Rep 2024; 14:174. [PMID: 38168773 PMCID: PMC10761953 DOI: 10.1038/s41598-023-50686-0] [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: 10/24/2023] [Accepted: 12/22/2023] [Indexed: 01/05/2024] Open
Abstract
Xanthine oxidase (XO) is a crucial enzyme in the development of hyperuricemia and gout. This study focuses on LWM and ALPM, two food-derived inhibitors of XO. We used molecular docking to obtain three systems and then conducted 200 ns molecular dynamics simulations for the Apo, LWM, and ALPM systems. The results reveal a stronger binding affinity of the LWM peptide to XO, potentially due to increased hydrogen bond formation. Notable changes were observed in the XO tunnel upon inhibitor binding, particularly with LWM, which showed a thinner, longer, and more twisted configuration compared to ALPM. The study highlights the importance of residue F914 in the allosteric pathway. Methodologically, we utilized the perturbed response scan (PRS) based on Python, enhancing tools for MD analysis. These findings deepen our understanding of food-derived anti-XO inhibitors and could inform the development of food-based therapeutics for reducing uric acid levels with minimal side effects.
Collapse
Affiliation(s)
- Yi He
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Qianjin Road 2699, Changchun, 130012, China
| | - Kaifeng Liu
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Qianjin Road 2699, Changchun, 130012, China
| | - Fuyan Cao
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Qianjin Road 2699, Changchun, 130012, China
| | - Renxiu Song
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Qianjin Road 2699, Changchun, 130012, China
| | - Jianxuan Liu
- Jilin Academy of Chinese Medicine Sciences, Chuangju Road 155, Changchun, 130012, China
| | - Yinghua Zhang
- Jilin Academy of Chinese Medicine Sciences, Chuangju Road 155, Changchun, 130012, China.
| | - Wannan Li
- Edmond H. Fischer Signal Transduction Laboratory, School of Life Sciences, Jilin University, Qianjin Road 2699, Changchun, 130012, China.
| | - Weiwei Han
- Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, School of Life Sciences, Jilin University, Qianjin Road 2699, Changchun, 130012, China.
| |
Collapse
|
17
|
Liu X, Zheng L, Qin C, Li Y, Zhang JZH, Sun Z. Screening Power of End-Point Free-Energy Calculations in Cucurbituril Host-Guest Systems. J Chem Inf Model 2023; 63:6938-6946. [PMID: 37908066 DOI: 10.1021/acs.jcim.3c01356] [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: 11/02/2023]
Abstract
End-point free-energy methods as an indispensable component in virtual screening are commonly recognized as a tool with a certain level of screening power in pharmaceutical research. While a huge number of records could be found for end-point applications in protein-ligand, protein-protein, and protein-DNA complexes from academic and industrial reports, up to now, there is no large-scale benchmark in host-guest complexes supporting the screening power of end-point free-energy techniques. A good benchmark requires a data set of sufficient coverage of pharmaceutically relevant chemical space, a long-time sampling length supporting the trajectory approximation of the ensemble average, and a sufficient sample size of receptor-acceptor pairs to stabilize the performance statistics. In this work, selecting a popular family of macrocyclic hosts named cucurbiturils, we construct a large data set containing 154 host-guest pairs, perform extensive end-point sampling of several hundred nanosecond lengths for each system, and extract the free-energy estimates with a variety of end-point free-energy techniques, including the advanced three-trajectory dielectric-constant-variable regime proposed in our recent work. The best-performing end-point protocol employs GAFF2 for solute descriptions, the three-trajectory end-point sampling regime, and the MM/GBSA Hamiltonian in free-energy extraction, achieving a high ranking metrics of Kendall τ > 0.6, a Pearlman predictive index of ∼0.8, and a high scoring power of Pearson r > 0.8. The current project as the first large-scale systematic benchmark of end-point methods in host-guest complexes in academic publications provides solid evidence of the applicability of end-point techniques and direct guidance of computational setups in practical host-guest systems.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai 201620, China
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Chu Qin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai 201620, China
| | - Yang Li
- College of Information Science and Engineering, Shandong Agricultural University, Tai'an 271018, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- Faculty of Synthetic Biology and Institute of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen Guangdong 518055, China
| | - Zhaoxi Sun
- Changping Laboratory, Beijing 102206, China
| |
Collapse
|
18
|
Wu J, Zhang HX, Zhang J. The molecular mechanism of non-covalent inhibitor WU-04 targeting SARS-CoV-2 3CLpro and computational evaluation of its effectiveness against mainstream coronaviruses. Phys Chem Chem Phys 2023; 25:23555-23567. [PMID: 37655706 DOI: 10.1039/d3cp03828a] [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: 09/02/2023]
Abstract
There is an urgent need for highly effective therapeutic agents to interrupt the continued spread of SARS-CoV-2. As a pivotal protease in the replication process of coronaviruses, the 3CLpro protein is considered as a potential target of drug development to stop the spread and infection of the virus. In this work, molecular dynamics (MD) simulations were used to elucidate the molecular mechanism of a novel and highly effective non-covalent inhibitor, WU-04, targeting the SARS-CoV-2 3CLpro protein. The difference in dynamic behavior between the apo-3CLpro and the holo-3CLpro systems suggests that the presence of WU-04 inhibits the motion amplitude of the 3CLpro protein relative to the apo-3CLpro system, thus maintaining a stable conformational binding state. The energy calculations and interaction analysis show that the hot-spot residues Q189, M165, M49, E166, and H41 and the warm-spot residues H163 and C145 have a strong binding capacity to WU-04 by forming multiple hydrogen bonds and hydrophobic interactions, which stabilizes the binding of the inhibitor. After that, the resistance of WU-04 to the six SARS-CoV-2 variants (Alpha, Beta, Gamma, Delta, Lambda, and Omicron) and two other mainstream coronavirus (SARS-CoV and MERS-CoV) 3CLpro proteins was further investigated. Excitingly, the slight difference in energy values relative to the SARS-CoV-2 system indicates that WU-04 is still highly effective against the coronaviruses, which becomes crucial evidence that WU-04 is a pan-inhibitor of the 3CLpro protein in various SARS-CoV-2 variants and other mainstream coronaviruses. The study will hopefully provide theoretical insights for the future rational design and improvement of novel non-covalent inhibitors targeting the 3CLpro protein.
Collapse
Affiliation(s)
- Jianhua Wu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People's Republic of China.
| | - Hong-Xing Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People's Republic of China.
| | - Jilong Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People's Republic of China.
| |
Collapse
|
19
|
Zhang X, Li L, Zheng Q. Dissecting the Effect of Temperature on Hyperthermophilic Pf2001 Esterase Dimerization by Molecular Dynamics. J Chem Inf Model 2023; 63:4762-4771. [PMID: 37452749 DOI: 10.1021/acs.jcim.3c00415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2023]
Abstract
Pf2001 esterase (Pf2001) from Pyrococcus furiosus has hyperthermophilic properties and exerts a biocatalytic function in a dimeric state. Crystal structures revealed that the structural rearrangement of the cap domain is responsible for the Pf2001 dimer formation. However, the details of the cap domain remodeling and the effects of temperature on the dimerization process remain elusive at the molecular level, taking into account that experimental methods are difficult to capture the dynamic process of dimerization to some extent. Herein, four dimer models based on the monomeric crystal structure (PDB ID: 5G59) were constructed to investigate the conformational transition details and temperature effects in the dimerization by conventional molecular dynamics and accelerated molecular dynamics simulations. Our simulation results indicate that the monomer undergoes a conformational change into a "preparatory state" at high temperatures, which is more favorable for its transformation into a stable dimer. The subsequent free energy landscape analysis further identifies four intermediate states (from separated state to dimeric state) and discloses that a more accessible α-helix driven by stronger hydrophobic interactions induces a rearrangement of the cap domain, displaying a "tic-tac-toe" activation feature that is important for stabilizing the dimer interface and facilitating the formation of hydrophobic pockets. In addition, the electrostatic potential surface analysis illustrates that the weaker electrostatic repulsion (Lys and Arg) in the dimer interface at high temperatures is also a key factor for dimer stabilization. Altogether, our results can provide molecular-level insight into the dimer formation process of hyperthermophilic esterase and would be useful to understand the enzymatic specificity of α/β-hydrolase.
Collapse
Affiliation(s)
- Xue Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Lei Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| | - Qingchuan Zheng
- School of Pharmaceutical Sciences, Jilin University, Changchun 130021, China
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, China
| |
Collapse
|
20
|
Liao J, Shu Z, Gao J, Wu M, Chen C. SurfPB: A GPU-Accelerated Electrostatic Calculation and Visualization Tool for Biomolecules. J Chem Inf Model 2023; 63:4490-4496. [PMID: 37500509 DOI: 10.1021/acs.jcim.3c00745] [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: 07/29/2023]
Abstract
In this work, we present SurfPB as a useful tool for the study of biomolecules. It can do many typical calculations, including the molecular surface, electrostatic potential, solvation free energy, entropy, and binding free energy. Among all of the calculations, the entropy calculation is the most time-consuming one. In SurfPB, the calculation can be performed in a vacuum or implicit solvent and accelerated on GPU. The Poisson-Boltzmann equation solver is accelerated on GPU as well. Moreover, we developed a graphical user interface for SurfPB. It allows users to input the parameters and complete the whole calculation in a visual way. The calculated electrostatic potentials are shown on the molecular surface in a three-dimensional scene.
Collapse
Affiliation(s)
- Jun Liao
- Biomolecular Physics and Modeling Group, School of Physics Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Zirui Shu
- Biomolecular Physics and Modeling Group, School of Physics Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Junyong Gao
- Biomolecular Physics and Modeling Group, School of Physics Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Mincong Wu
- Biomolecular Physics and Modeling Group, School of Physics Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| | - Changjun Chen
- Biomolecular Physics and Modeling Group, School of Physics Huazhong University of Science and Technology, Wuhan 430074, Hubei, China
| |
Collapse
|
21
|
Imran M, Abida, Alotaibi NM, Thabet HK, Alruwaili JA, Eltaib L, Alshehri A, Alsaiari AA, Kamal M, Alshammari AMA. Repurposing Anti-Dengue Compounds against Monkeypox Virus Targeting Core Cysteine Protease. Biomedicines 2023; 11:2025. [PMID: 37509664 PMCID: PMC10377189 DOI: 10.3390/biomedicines11072025] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2023] [Revised: 07/04/2023] [Accepted: 07/06/2023] [Indexed: 07/30/2023] Open
Abstract
The monkeypox virus (MPXV) is an enveloped, double-stranded DNA virus belonging to the genus Orthopox viruses. In recent years, the virus has spread to countries where it was previously unknown, turning it into a worldwide emergency for public health. This study employs a structural-based drug design approach to identify potential inhibitors for the core cysteine proteinase of MPXV. During the simulations, the study identified two potential inhibitors, compound CHEMBL32926 and compound CHEMBL4861364, demonstrating strong binding affinities and drug-like properties. Their docking scores with the target protein were -10.7 and -10.9 kcal/mol, respectively. This study used ensemble-based protein-ligand docking to account for the binding site conformation variability. By examining how the identified inhibitors interact with the protein, this research sheds light on the workings of the inhibitors' mechanisms of action. Molecular dynamic simulations of protein-ligand complexes showed fluctuations from the initial docked pose, but they confirmed their binding throughout the simulation. The MMGBSA binding free energy calculations for CHEMBL32926 showed a binding free energy range of (-9.25 to -9.65) kcal/mol, while CHEMBL4861364 exhibited a range of (-41.66 to -31.47) kcal/mol. Later, analogues were searched for these compounds with 70% similarity criteria, and their IC50 was predicted using pre-trained machine learning models. This resulted in identifying two similar compounds for each hit with comparable binding affinity for cysteine proteinase. This study's structure-based drug design approach provides a promising strategy for identifying new drugs for treating MPXV infections.
Collapse
Affiliation(s)
- Mohd Imran
- Department of Pharmaceutical Chemistry, College of Pharmacy, Northern Border University, Rafha 91911, Saudi Arabia;
| | - Abida
- Department of Pharmaceutical Chemistry, College of Pharmacy, Northern Border University, Rafha 91911, Saudi Arabia;
| | - Nawaf M. Alotaibi
- Department of Clinical Pharmacy, College of Pharmacy, Northern Border University, Rafha 91911, Saudi Arabia
| | - Hamdy Khamees Thabet
- Chemistry Department, College of Arts and Sciences, Northern Border University, Rafha 91911, Saudi Arabia
| | - Jamal Alhameedi Alruwaili
- Medical Lab Technology Department, College of Applied Medical Sciences, Northern Border University, Arar 91431, Saudi Arabia
| | - Lina Eltaib
- Department of Pharmaceutics, College of Pharmacy, Northern Border University, Rafha 91911, Saudi Arabia
| | - Ahmed Alshehri
- Department of Pharmacology and Toxicology, College of Pharmacy, Northern Border University, Rafha 91911, Saudi Arabia
- Department of Pharmacology, College of Clinical Pharmacy, Imam Abdulrahman Bin Faisal University, King Faisal Road, Dammam 31441, Saudi Arabia
| | - Ahad Amer Alsaiari
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
| | - Mehnaz Kamal
- Department of Pharmaceutical Chemistry, College of Pharmacy, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia;
| | | |
Collapse
|
22
|
Jiang W, Yang Z, Chen P, Zhao M, Wang Y, Wang J, Li X, Wang M, Hou P. Tolperisone induces UPR-mediated tumor inhibition and synergizes with proteasome inhibitor and immunotherapy by targeting LSD1. Expert Opin Ther Targets 2023; 27:879-895. [PMID: 37704953 DOI: 10.1080/14728222.2023.2259097] [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: 05/27/2023] [Accepted: 09/11/2023] [Indexed: 09/15/2023]
Abstract
BACKGROUND Drug repurposing is an attractive strategy for extending the arsenal of oncology therapies. Tolperisone is an old centrally acting muscle relaxant used for treatment of chronic pain conditions. In this study, we investigated the therapeutic effect and mechanism of tolperisone in human cancers and explored the combination strategy with proteasome inhibitor and immunotherapy. RESEARCH DESIGN AND METHODS The antitumor effect of tolperisone was evaluated by measuring half maximal inhibitory concentration, cell death, and cell growth. RNA sequencing, western blotting, molecular docking, enzyme activity assay, and ChIP-qPCR were performed to reveal the underlying mechanism. Xenograft models were used to evaluate the efficacy of tolperisone alone or in combination with proteasome inhibitor or immunotherapy. RESULTS Tolperisone inhibited cell growth and induced cell death in human cancer cell lines. Unfolded protein responses (UPR) pathway was hyperactivated in tolperisone-treated cells. We further identified histone lysine-specific demethylase 1 (LSD1) as a potential target of tolperisone, which directly demethylates UPR-related genes in H3K4me2. Tolperisone synergistically improved the efficacy of MG132 by enhancing UPR and sensitized tumors to immunotherapy by reprogramming M2 macrophages into M1 phenotype. CONCLUSIONS Tolperisone inhibits human cancer by targeting LSD1. Repurposing tolperisone in cancer therapy by a combination strategy implies clinical potential.
Collapse
Affiliation(s)
- Wei Jiang
- Key Laboratory for Tumor Precision Medicine of Shaanxi Province and Department of Endocrinology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| | - Zhiwei Yang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Xi'an Jiaotong University, Xi'an, P.R. China
| | - Pu Chen
- Key Laboratory for Tumor Precision Medicine of Shaanxi Province and Department of Endocrinology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| | - Man Zhao
- Key Laboratory for Tumor Precision Medicine of Shaanxi Province and Department of Endocrinology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| | - Yubo Wang
- Key Laboratory for Tumor Precision Medicine of Shaanxi Province and Department of Endocrinology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| | - Jingyuan Wang
- Department of Clinical Lab Diagnosis, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| | - Xinru Li
- Key Laboratory for Tumor Precision Medicine of Shaanxi Province and Department of Endocrinology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| | - Meichen Wang
- Key Laboratory for Tumor Precision Medicine of Shaanxi Province and Department of Endocrinology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| | - Peng Hou
- Key Laboratory for Tumor Precision Medicine of Shaanxi Province and Department of Endocrinology, The First Affiliated Hospital of Xi'an Jiaotong University, Xi'an, P.R. China
| |
Collapse
|
23
|
Tang R, Wang Z, Xiang S, Wang L, Yu Y, Wang Q, Deng Q, Hou T, Sun H. Uncovering the Kinetic Characteristics and Degradation Preference of PROTAC Systems with Advanced Theoretical Analyses. JACS AU 2023; 3:1775-1789. [PMID: 37388700 PMCID: PMC10301679 DOI: 10.1021/jacsau.3c00195] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Revised: 05/16/2023] [Accepted: 05/16/2023] [Indexed: 07/01/2023]
Abstract
Proteolysis-targeting chimeras (PROTACs), which can selectively induce the degradation of target proteins, represent an attractive technology in drug discovery. A large number of PROTACs have been reported, but due to the complicated structural and kinetic characteristics of the target-PROTAC-E3 ligase ternary interaction process, the rational design of PROTACs is still quite challenging. Here, we characterized and analyzed the kinetic mechanism of MZ1, a PROTAC that targets the bromodomain (BD) of the bromodomain and extra terminal (BET) protein (Brd2, Brd3, or Brd4) and von Hippel-Lindau E3 ligase (VHL), from the kinetic and thermodynamic perspectives of view by using enhanced sampling simulations and free energy calculations. The simulations yielded satisfactory predictions on the relative residence time and standard binding free energy (rp > 0.9) for MZ1 in different BrdBD-MZ1-VHL ternary complexes. Interestingly, the simulation of the PROTAC ternary complex disintegration illustrates that MZ1 tends to remain on the surface of VHL with the BD proteins dissociating alone without a specific dissociation direction, indicating that the PROTAC prefers more to bind with E3 ligase at the first step in the formation of the target-PROTAC-E3 ligase ternary complex. Further exploration of the degradation difference of MZ1 in different Brd systems shows that the PROTAC with higher degradation efficiency tends to leave more lysine exposed on the target protein, which is guaranteed by the stability (binding affinity) and durability (residence time) of the target-PROTAC-E3 ligase ternary complex. It is quite possible that the underlying binding characteristics of the BrdBD-MZ1-VHL systems revealed by this study may be shared by different PROTAC systems as a general rule, which may accelerate rational PROTAC design with higher degradation efficiency.
Collapse
Affiliation(s)
- Rongfan Tang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Zhe Wang
- Innovation
Institute for Artificial Intelligence in Medicine of Zhejiang University,
College of Pharmaceutical Sciences, Zhejiang
University, Hangzhou 310058, Zhejiang, P. R. China
| | - Sutong Xiang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Lingling Wang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Yang Yu
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Qinghua Wang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Qirui Deng
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Tingjun Hou
- Innovation
Institute for Artificial Intelligence in Medicine of Zhejiang University,
College of Pharmaceutical Sciences, Zhejiang
University, Hangzhou 310058, Zhejiang, P. R. China
| | - Huiyong Sun
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| |
Collapse
|
24
|
Gu S, Shen C, Yu J, Zhao H, Liu H, Liu L, Sheng R, Xu L, Wang Z, Hou T, Kang Y. Can molecular dynamics simulations improve predictions of protein-ligand binding affinity with machine learning? Brief Bioinform 2023; 24:6995375. [PMID: 36681903 DOI: 10.1093/bib/bbad008] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2022] [Revised: 12/04/2022] [Accepted: 12/30/2023] [Indexed: 01/23/2023] Open
Abstract
Binding affinity prediction largely determines the discovery efficiency of lead compounds in drug discovery. Recently, machine learning (ML)-based approaches have attracted much attention in hopes of enhancing the predictive performance of traditional physics-based approaches. In this study, we evaluated the impact of structural dynamic information on the binding affinity prediction by comparing the models trained on different dimensional descriptors, using three targets (i.e. JAK1, TAF1-BD2 and DDR1) and their corresponding ligands as the examples. Here, 2D descriptors are traditional ECFP4 fingerprints, 3D descriptors are the energy terms of the Smina and NNscore scoring functions and 4D descriptors contain the structural dynamic information derived from the trajectories based on molecular dynamics (MD) simulations. We systematically investigate the MD-refined binding affinity prediction performance of three classical ML algorithms (i.e. RF, SVR and XGB) as well as two common virtual screening methods, namely Glide docking and MM/PBSA. The outcomes of the ML models built using various dimensional descriptors and their combinations reveal that the MD refinement with the optimized protocol can improve the predictive performance on the TAF1-BD2 target with considerable structural flexibility, but not for the less flexible JAK1 and DDR1 targets, when taking docking poses as the initial structure instead of the crystal structures. The results highlight the importance of the initial structures to the final performance of the model through conformational analysis on the three targets with different flexibility.
Collapse
Affiliation(s)
- Shukai Gu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Chao Shen
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Jiahui Yu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Hong Zhao
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Huanxiang Liu
- Faculty of Applied Science, Macao Polytechnic University, Macao, SAR, China
| | - Liwei Liu
- Advanced Computing and Storage Laboratory, Central Research Institute, 2012 Laboratories, Huawei Technologies Co., Ltd., Shenzhen 518129, Guangdong, China
| | - Rong Sheng
- Health Technology Development Dept, Huawei Device Co., Ltd., Dongguan 523808, Guangdong, China
| | - Lei Xu
- Institute of Bioinformatics and Medical Engineering, School of Electrical and Information Engineering, Jiangsu University of Technology, Changzhou 213001, China
| | - Zhe Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| | - Yu Kang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, China
| |
Collapse
|
25
|
Sun Y, He X, Hou T, Cai L, Man VH, Wang J. Development and test of highly accurate endpoint free energy methods. 1: Evaluation of ABCG2 charge model on solvation free energy prediction and optimization of atom radii suitable for more accurate solvation free energy prediction by the PBSA method. J Comput Chem 2023; 44:1334-1346. [PMID: 36807356 DOI: 10.1002/jcc.27089] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 12/22/2022] [Accepted: 01/29/2023] [Indexed: 02/23/2023]
Abstract
Accurate estimation of solvation free energy (SFE) lays the foundation for accurate prediction of binding free energy. The Poisson-Boltzmann (PB) or generalized Born (GB) combined with surface area (SA) continuum solvation method (PBSA and GBSA) have been widely used in SFE calculations because they can achieve good balance between accuracy and efficiency. However, the accuracy of these methods can be affected by several factors such as the charge models, polar and nonpolar SFE calculation methods and the atom radii used in the calculation. In this work, the performance of the ABCG2 (AM1-BCC-GAFF2) charge model as well as other two charge models, that is, RESP (Restrained Electrostatic Potential) and AM1-BCC (Austin Model 1-bond charge corrections), on the SFE prediction of 544 small molecules in water by PBSA/GBSA was evaluated. In order to improve the performance of the PBSA prediction based on the ABCG2 charge, we further explored the influence of atom radii on the prediction accuracy and yielded a set of atom radius parameters for more accurate SFE prediction using PBSA based on the ABCG2/GAFF2 by reproducing the thermodynamic integration (TI) calculation results. The PB radius parameters of carbon, oxygen, sulfur, phosphorus, chloride, bromide and iodine, were adjusted. New atom types, on, oi, hn1, hn2, hn3, were introduced to further improve the fitting performance. Then, we tuned the parameters in the nonpolar SFE model using the experimental SFE data and the PB calculation results. By adopting the new radius parameters and new nonpolar SFE model, the root mean square error (RMSE) of the SFE calculation for the 544 molecules decreased from 2.38 to 1.05 kcal/mol. Finally, the new radius parameters were applied in the prediction of protein-ligand binding free energies using the MM-PBSA method. For the eight systems tested, we could observe higher correlation between the experiment data and calculation results and smaller prediction errors for the absolute binding free energies, demonstrating that our new radius parameters can improve the free energy calculation using the MM-PBSA method.
Collapse
Affiliation(s)
- Yuchen Sun
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, China
| | - Lianjin Cai
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Viet Hoag Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| |
Collapse
|
26
|
Wu J, Zhang HX, Zhang J. Investigation on the interaction mechanism of different SARS-CoV-2 spike variants with hACE2: insights from molecular dynamics simulations. Phys Chem Chem Phys 2023; 25:2304-2319. [PMID: 36597957 DOI: 10.1039/d2cp04349a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Since the COVID-19 pandemic caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), SARS-CoV-2 has evolved by acquiring genomic mutations, resulting in the recent emergence of several SARS-CoV-2 variants with improved transmissibility and infectivity relative to the original strain. An underlying mechanism may be the increased ability of the mutants to bind the receptor proteins and infect the host cell. In this work, we implemented all-atom molecular dynamics (MD) simulations to study the binding and interaction of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein singly (D614G), doubly (D614G + L452R and D614G + N501Y), triply (D614G + N501Y + E484K), and quadruply (D614G + N501Y + E484K + K417T) mutated variants with the human angiotensin-converting enzyme 2 (hACE2) receptor protein in the host cell. A combination of multiple analysis approaches elucidated the effects of mutations and the extent of molecular divergence from multiple perspectives, including the dynamic correlated motions, interaction patterns, dominant motions, free energy landscape, and charge distribution on the electrostatic potential surface between the hACE2 and all RBD variants. Moreover, free energy calculations using the MM/PBSA method evaluated the binding affinity between these RBD variants and hACE2. The results showed that the D614G + N501Y + E484K variant possessed the lowest free energy value (highest affinity) compared to the D614G + N501Y + E484K + K417T, D614G + L452R, D614G + N501Y, and D614G mutants. The residue-based energy decomposition also indicated that the energy contribution of residues at the mutation site to the total binding energy was highly variable. The interaction mechanisms between the different RBD variants and hACE2 elucidated in this study will provide some insights into the development of drugs targeting the new SARS-CoV-2 variants.
Collapse
Affiliation(s)
- Jianhua Wu
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People's Republic of China.
| | - Hong-Xing Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People's Republic of China.
| | - Jilong Zhang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People's Republic of China.
| |
Collapse
|
27
|
Computationally guided discovery of novel non-steroidal AR-GR dual antagonists demonstrating potency against antiandrogen resistance. Acta Pharmacol Sin 2023:10.1038/s41401-022-01038-7. [PMID: 36639570 DOI: 10.1038/s41401-022-01038-7] [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: 10/28/2022] [Accepted: 12/01/2022] [Indexed: 01/14/2023] Open
Abstract
As a major class of medicine for treating the lethal type of castration-resistant prostate cancer (PCa), long-term use of androgen receptor (AR) antagonists commonly leads to antiandrogen resistance. When AR signaling pathway is blocked by AR-targeted therapy, glucocorticoid receptor (GR) could compensate for AR function especially at the late stage of PCa. AR-GR dual antagonist is expected to be a good solution for this situation. Nevertheless, no effective non-steroidal AR-GR dual antagonist has been reported so far. In this study, an AR-GR dual binder H18 was first discovered by combining structure-based virtual screening and biological evaluation. Then with the aid of computationally guided design, the AR-GR dual antagonist HD57 was finally identified with antagonistic activity towards both AR (IC50 = 0.394 μM) and GR (IC50 = 17.81 μM). Moreover, HD57 could effectively antagonize various clinically relevant AR mutants. Further molecular dynamics simulation provided more atomic insights into the mode of action of HD57. Our research presents an efficient and rational strategy for discovering novel AR-GR dual antagonists, and the new scaffold provides important clues for the development of novel therapeutics for castration-resistant PCa.
Collapse
|
28
|
High-Affinity Antibodies Designing of SARS-CoV-2 Based on Molecular Dynamics Simulations. Int J Mol Sci 2022; 24:ijms24010481. [PMID: 36613923 PMCID: PMC9820416 DOI: 10.3390/ijms24010481] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2022] [Revised: 12/22/2022] [Accepted: 12/23/2022] [Indexed: 12/29/2022] Open
Abstract
SARS-CoV-2 has led to a global pandemic of new crown pneumonia, which has had a tremendous impact on human society. Antibody drug therapy is one of the most effective way of combating SARS-CoV-2. In order to design potential antibody drugs with high affinity, we used antibody S309 from patients with SARS-CoV as the target antibody and RBD of S protein as the target antigen. Systems with RBD glycosylated and non-glycosylated were constructed to study the influence of glycosylation. From the results of molecular dynamics simulations, the steric effects of glycans on the surface of RBD plays a role of "wedge", which makes the L335-E340 region of RBD close to the CDR3 region of the heavy chain of antibody and increases the contact area between antigen and antibody. By mutating the key residues of antibody at the interaction interface, we found that the binding affinities of antibody mutants G103A, P28W and Y100W were all stronger than that of the wild-type, especially for the G103A mutant. G103A significantly reduces the distance between the binding region of L335-K356 in the antigen and P28-Y32 of heavy chain in the antibody through structural transition. Taken together, the antibody design method described in this work can provide theoretical guidance and a time-saving method for antibody drug design.
Collapse
|
29
|
Liu X, Zheng L, Cong Y, Gong Z, Yin Z, Zhang JZH, Liu Z, Sun Z. Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: II. regression and dielectric constant. J Comput Aided Mol Des 2022; 36:879-894. [PMID: 36394776 DOI: 10.1007/s10822-022-00487-w] [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/26/2022] [Accepted: 10/29/2022] [Indexed: 11/18/2022]
Abstract
End-point free energy calculations as a powerful tool have been widely applied in protein-ligand and protein-protein interactions. It is often recognized that these end-point techniques serve as an option of intermediate accuracy and computational cost compared with more rigorous statistical mechanic models (e.g., alchemical transformation) and coarser molecular docking. However, it is observed that this intermediate level of accuracy does not hold in relatively simple and prototypical host-guest systems. Specifically, in our previous work investigating a set of carboxylated-pillar[6]arene host-guest complexes, end-point methods provide free energy estimates deviating significantly from the experimental reference, and the rank of binding affinities is also incorrectly computed. These observations suggest the unsuitability and inapplicability of standard end-point free energy techniques in host-guest systems, and alteration and development are required to make them practically usable. In this work, we consider two ways to improve the performance of end-point techniques. The first one is the PBSA_E regression that varies the weights of different free energy terms in the end-point calculation procedure, while the second one is considering the interior dielectric constant as an additional variable in the end-point equation. By detailed investigation of the calculation procedure and the simulation outcome, we prove that these two treatments (i.e., regression and dielectric constant) are manipulating the end-point equation in a somehow similar way, i.e., weakening the electrostatic contribution and strengthening the non-polar terms, although there are still many detailed differences between these two methods. With the trained end-point scheme, the RMSE of the computed affinities is improved from the standard ~ 12 kcal/mol to ~ 2.4 kcal/mol, which is comparable to another altered end-point method (ELIE) trained with system-specific data. By tuning PBSA_E weighting factors with the host-specific data, it is possible to further decrease the prediction error to ~ 2.1 kcal/mol. These observations along with the extremely efficient optimized-structure computation procedure suggest the regression (i.e., PBSA_E as well as its GBSA_E extension) as a practically applicable solution that brings end-point methods back into the library of usable tools for host-guest binding. However, the dielectric-constant-variable scheme cannot effectively minimize the experiment-calculation discrepancy for absolute binding affinities, but is able to improve the calculation of affinity ranks. This phenomenon is somehow different from the protein-ligand case and suggests the difference between host-guest and biomacromolecular (protein-ligand and protein-protein) systems. Therefore, the spectrum of tools usable for protein-ligand complexes could be unsuitable for host-guest binding, and numerical validations are necessary to screen out really workable solutions in these 'prototypical' situations.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China.
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| | - Yalong Cong
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Zhihao Gong
- School of Micro-Nano Electronics, Zhejiang University, Hangzhou, 310027, China.,Hangzhou Global Scientific and Technological Innovation Center, Zhejiang University, Hangzhou, 310027, China
| | - Zhixiang Yin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China. .,School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China. .,Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China. .,Department of Chemistry, New York University, NY, NY, 10003, USA.
| | - Zhirong Liu
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
30
|
Wu J, Zhang J, Zhang HX. Computational Design of Miniprotein Inhibitors Targeting SARS-CoV-2 Spike Protein. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2022; 38:10690-10703. [PMID: 35984970 PMCID: PMC9437664 DOI: 10.1021/acs.langmuir.2c01699] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Revised: 07/31/2022] [Indexed: 05/16/2023]
Abstract
The ongoing pandemic of COVID-19 caused by SARS-CoV-2 has become a global health problem. There is an urgent need to develop therapeutic drugs, effective therapies, and vaccines to prevent the spread of the virus. The virus first enters the host cell through the interaction between the receptor binding domain (RBD) of spike protein and the peptidase domain (PD) of the angiotensin-converting enzyme 2 (ACE2). Therefore, blocking the binding of RBD and ACE2 is a promising strategy to inhibit the invasion and infection of the virus in the host cell. In the study, we designed several miniprotein inhibitors against SARS-CoV-2 by single/double/triple-point mutant, based on the initial inhibitor LCB3. Molecular dynamics (MD) simulations and trajectory analysis were performed for an in-depth analysis of the structural stability, essential protein motions, and per-residue energy decomposition involved in the interaction of inhibitors with the RBD. The results showed that the inhibitors have adapted the protein RBD in the binding interface, thereby forming stable complexes. These inhibitors display low binding free energy in the MM/PBSA calculations, substantiating their strong interaction with RBD. Moreover, the binding affinity of the best miniprotein inhibitor, H6Y-M7L-L17F mutant, to RBD was ∼45 980 times (ΔG = RT ln Ki) higher than that of the initial inhibitor LCB3. Following H6Y-M7L-L17F mutant, the inhibitors with strong binding activity are successively H6Y-L17F, L17F, H6Y, and F30Y mutants. Our research proves that the miniprotein inhibitors can maintain their secondary structure and have a highly stable blocking (binding) effect on SARS-CoV-2. This study proposes novel miniprotein mutant inhibitors with enhanced binding to spike protein and provides potential guidance for the rational design of new SARS-CoV-2 spike protein inhibitors.
Collapse
Affiliation(s)
- Jianhua Wu
- Institute
of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People’s Republic of China
| | - Jilong Zhang
- Institute
of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People’s Republic of China
| | - Hong-Xing Zhang
- Institute
of Theoretical Chemistry, College of Chemistry, Jilin University, Changchun 130023, Jilin, People’s Republic of China
| |
Collapse
|
31
|
Yu Y, Wang Z, Wang L, Tian S, Hou T, Sun H. Predicting the mutation effects of protein–ligand interactions via end-point binding free energy calculations: strategies and analyses. J Cheminform 2022; 14:56. [PMID: 35987841 PMCID: PMC9392442 DOI: 10.1186/s13321-022-00639-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Accepted: 08/08/2022] [Indexed: 12/04/2022] Open
Abstract
Protein mutations occur frequently in biological systems, which may impact, for example, the binding of drugs to their targets through impairing the critical H-bonds, changing the hydrophobic interactions, etc. Thus, accurately predicting the effects of mutations on biological systems is of great interests to various fields. Unfortunately, it is still unavailable to conduct large-scale wet-lab mutation experiments because of the unaffordable experimental time and financial costs. Alternatively, in silico computation can serve as a pioneer to guide the experiments. In fact, numerous pioneering works have been conducted from computationally cheaper machine-learning (ML) methods to the more expensive alchemical methods with the purpose to accurately predict the mutation effects. However, these methods usually either cannot result in a physically understandable model (ML-based methods) or work with huge computational resources (alchemical methods). Thus, compromised methods with good physical characteristics and high computational efficiency are expected. Therefore, here, we conducted a comprehensive investigation on the mutation issues of biological systems with the famous end-point binding free energy calculation methods represented by MM/GBSA and MM/PBSA. Different computational strategies considering different length of MD simulations, different value of dielectric constants and whether to incorporate entropy effects to the predicted total binding affinities were investigated to provide a more accurate way for predicting the energetic change upon protein mutations. Overall, our result shows that a relatively long MD simulation (e.g. 100 ns) benefits the prediction accuracy for both MM/GBSA and MM/PBSA (with the best Pearson correlation coefficient between the predicted ∆∆G and the experimental data of ~ 0.44 for a challenging dataset). Further analyses shows that systems involving large perturbations (e.g. multiple mutations and large number of atoms change in the mutation site) are much easier to be accurately predicted since the algorithm works more sensitively to the large change of the systems. Besides, system-specific investigation reveals that conformational adjustment is needed to refine the micro-environment of the manually mutated systems and thus lead one to understand why longer MD simulation is necessary to improve the predicting result. The proposed strategy is expected to be applied in large-scale mutation effects investigation with interpretation.
Collapse
|
32
|
Gutiérrez LJ, Tosso RD, Zarycz MNC, Enriz RD, Baldoni HA. Epitopes mapped onto SARS-CoV-2 receptor-binding motif by five distinct human neutralising antibodies. MOLECULAR SIMULATION 2022. [DOI: 10.1080/08927022.2022.2111421] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
Affiliation(s)
- Lucas J. Gutiérrez
- Multidisciplinary Institute of Biological Research (IMIBIO-SL. CONICET), San Luis, Argentina
- Faculty of Chemistry, Biochemistry and Pharmacy, National University of San Luis, San Luis, Argentina
| | - Rodrigo D. Tosso
- Multidisciplinary Institute of Biological Research (IMIBIO-SL. CONICET), San Luis, Argentina
- Faculty of Chemistry, Biochemistry and Pharmacy, National University of San Luis, San Luis, Argentina
| | - M. Natalia C. Zarycz
- Multidisciplinary Institute of Biological Research (IMIBIO-SL. CONICET), San Luis, Argentina
| | - Ricardo D. Enriz
- Multidisciplinary Institute of Biological Research (IMIBIO-SL. CONICET), San Luis, Argentina
- Faculty of Chemistry, Biochemistry and Pharmacy, National University of San Luis, San Luis, Argentina
| | - Héctor A. Baldoni
- Faculty of Chemistry, Biochemistry and Pharmacy, National University of San Luis, San Luis, Argentina
- Institute of Applied Mathematics of San Luis (IMASL. CONICET), San Luis, Argentina
| |
Collapse
|
33
|
Yang Z, Zang Y, Wang H, Kang Y, Zhang J, Li X, Zhang L, Zhang S. Recognition between CD147 and cyclophilin A deciphered by accelerated molecular dynamics simulations. Phys Chem Chem Phys 2022; 24:18905-18914. [PMID: 35913096 DOI: 10.1039/d2cp01975b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
CD147 functions as the receptor of extracellular cyclophilin A (CypA) in various diseases, and CD147-CypA binding ulteriorly underlies the pathological process of various viral infections including HIV-1, SARS, and SARS-CoV-2. Although CyPA has been identified as a key intermediate pro-inflammatory factor, the mechanism by which CD147 cooperates with CypA in the development of the cytokine storm remains largely unknown, and the binding profile of CD147 with CypA remains to be elucidated as well. Here, we prepared three binding models of the CD147-CypA complex, including the active site of CypA severally binding to the groove bound by the Ig1 and Ig2 domains (model-0), P180-G181 (model-1), and P211 (model-2) of CD147, as well as introducing mutations P180A-G181A and P211A individually in each model. All systems were studied using accelerated molecular dynamics simulations and the molecular mechanics generalized Born surface area (MM/GBSA) method. For model-0, CypA bound to the ectodomain of CD147 with the highest binding affinity. Moreover, mutations P180A-G181A of CD147 in model-0 decreased the binding affinity and weakened the dynamic correlation between CD147 and CypA, which resulted in CypA shifting from the initial binding location. Other residue mutations of CD147 did not significantly affect the CD147-CypA binding, as reflected by the energy and structural analyses. Compared with surface plasmon resonance results and nuclear magnetic resonance shift signals, CypA should tend to reciprocally bind to the groove of CD147, and the binding process might be modulated by P180-G181 rather than P211. Besides, residue R201 of CD147 is critical for CD147-CypA binding and needs further experimental verification. These findings further our understanding of the recruitment between CD147 and CypA and its potential role in the development of inflammation and viral infection.
Collapse
Affiliation(s)
- Zhiwei Yang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Yongjian Zang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - He Wang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Ying Kang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Jianwen Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Xuhua Li
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Lei Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| | - Shengli Zhang
- MOE Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, School of Physics, Xi'an Jiaotong University, Xi'an 710049, China.
| |
Collapse
|
34
|
Priya P, Basit A, Bandyopadhyay P. A strategy to optimize the peptide-based inhibitors against different mutants of the spike protein of SARS-CoV-2. J Biomol Struct Dyn 2022:1-12. [PMID: 35881159 DOI: 10.1080/07391102.2022.2103587] [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: 10/16/2022]
Abstract
The SARS-CoV-2 virus has caused high-priority health concerns at a global level. Vaccines have stalled the proliferation of viruses to some extent. Yet, the emergence of newer, potentially more infectious, and dangerous mutants such as Delta and Omicron are among the major challenges in finding a more permanent solution for this pandemic. The effectiveness of antivirals Molnupiravir and Paxlovid, authorized for emergency use by the FDA, are yet to be assessed on a larger population. Patients with a high risk of disease progression have received treatment with antibody-cocktail. Most of the mutations leading to the new lineage of SARS-CoV-2 are found in the spike protein of this virus that plays a key role in facilitating host entry. The current study has investigated how to modify a promising peptide-based inhibitor of spike protein, LCB3, against common mutations, N501Y and K417N in the target protein so that it retains its efficacy against the spike protein. LCB3 being a prototype for protein-based inhibitors is an ideal testing system to learn about protein-based inhibitors. This study proposes the substitutions of amino acid residues of LCB3 inhibitor using a structure-based approach that considers free energy decomposition of residues, the distance between atoms, and the interaction among amino acids. The binding free energy calculations suggest a possible improvement in the binding affinity of existing inhibitor LCB3 to the mutant forms of the S-protein using simple substitutions at specific positions of the inhibitor. This approach, being general, can be used in different inhibitors and other mutations and help in fighting against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Prerna Priya
- Department of Botany, Purnea Mahila College, Purnea University, Purnia, Bihar, India
| | - Abdul Basit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| |
Collapse
|
35
|
Zhao X, Xiong D, Luo S, Duan L. Molecular investigation of the dual inhibition mechanism for targeted P53 regulator MDM2/MDMX inhibitors. Phys Chem Chem Phys 2022; 24:16799-16815. [PMID: 35775962 DOI: 10.1039/d2cp01780f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Inhibitors that competitively bind MDM2/MDMX can block the inhibition of P53 by MDM2/MDMX and restart its tumor-suppressive effect. Molecular studies targeting MDM2/MDMX inhibitors have always been a hot topic in anticancer drug design. Although numerous inhibitors have been designed previously against MDM2/MDMX, their dual inhibition efficacy has not been demonstrated, and few studies assessed the general causes affecting the dual inhibition of MDM2/MDMX by these inhibitors. Here, molecular dynamics simulations and alanine scanning combined with the interaction entropy method were employed to precisely investigate whether 16 inhibitors could dually inhibit MDM2/MDMX and the similarities and differences in the interaction modes. Thereby addressing the key residue sites affecting dual inhibition. Residues L54/M53, I61/60, M62/61, Y67/66, and V93/92 of MDM2/MDMX, which are in corresponding positions in both protein structures, provide significant conditions for these inhibitors to bind to MDM2/MDMX tightly. In addition, most of these inhibitors prefer to bind MDM2 than MDMX, and residues H96 and I99 in MDM2 are attractive targets for inhibitors, resulting in inhibitors binding to MDM2/MDMX with different affinity. These key residues should be considered in the development of dual inhibitors. For these 16 inhibitors, most have dual inhibitory potential for MDM2/MDMX based on the binding affinity of the complexes. Still, it is questionable whether they can exert excellent dual inhibition considering the assessment of the hot-spots. At least their binding affinity for MDMX is not superior to that for MDM2 due to the difference in energy of the van der Waals interactions at the key sites. Furthermore, based on the analysis of three representative inhibitors (TUZ/HRH and HRQ with different binding preferences for MDM2/MDMX), 3-chloropyridine in TUZ leads to the differential binding affinity between the inhibitor and MDM2/MDMX. It readily forms hydrophobic interactions with the surrounding residues H96 and I99. But this phenomenon does not occur in the TUZ-MDMX system, implying the critical role of residues H96/P95 and I99/L98. And the completely different binding mechanism of HRQ binding to MDM2/MDMX explains its inability to inhibit MDM2 well. Thus, we are cautious about its dual inhibitory ability. Besides, HRH is more prone to strong van der Waals interactions with MDM2 than MDMX whereas its 2-chlorofluorobenzene is detrimental to this. We hope that these findings will provide reliable molecular insights for the screening and optimization of targeting MDM2/MDMX dual inhibitors.
Collapse
Affiliation(s)
- Xiaoyu Zhao
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Danyang Xiong
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Song Luo
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| |
Collapse
|
36
|
Zhao X, Xiong D, Luo S, Duan L. Origin of the tight binding mode to ACE2 triggered by multi-point mutations in the omicron variant: a dynamic insight. Phys Chem Chem Phys 2022; 24:8724-8737. [PMID: 35373810 DOI: 10.1039/d2cp00449f] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The continuous spread of the newly emerged SARS-CoV-2 Omicron variant (B.1.1.529) has become an important reason for the surge in COVID-19 infections. Its numerous mutated residues containing key sites on the receptor-binding domain (RBD) undoubtedly pose new challenges for epidemic control. Although the preventive measures are becoming more sophisticated, the effects of mutations on the binding of the virus to the receptor protein remain to be elucidated. Here, we used molecular dynamics (MD) simulations to investigate the differences in the binding mode between the Omicron variant and the angiotensin-converting enzyme 2 (ACE2) compared to the wild-type strain (WT). Multi-point mutations in the Omicron variant RBD could cause the conformation shift in the large Loop (where T478K and E484A are located), which makes it easier to wrap the N-terminal helix of ACE2 and form tighter contacts. The stronger electrostatic interaction was the main reason for its enhanced binding affinity as compared to WT. This was due to the large number of positively charged patches (N440K, T478K, Q493R, Q498R, and Y505H) formed by the substitution of neutral amino acids at multiple sites. The appearance of these highly polar hydrophilic amino acids may cause local perturbations and affect the electrostatic complementarity of RBD with the ACE2, and further mediate conformational changes. Thus, a more extensive interaction network was found in the mutation system and the complex interaction cluster was formed near E37@ACE2, which was essential for the stable binding of the two. In addition, we speculated that these mutations may affect the electrostatic complementarity with the four potential antibodies to reduce the sensitivity of the virus to antibodies. This study reveals the key details of the Omicron variant binding to ACE2 and provides important theoretical views for the enhanced infectivity of this variant. We hope that these observations can provide timely molecular insights for responding to the Omicron variant pandemic.
Collapse
Affiliation(s)
- Xiaoyu Zhao
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Danyang Xiong
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Song Luo
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan, 250014, China.
| |
Collapse
|
37
|
Li Z, Zhang JZH. Mutational Effect of Some Major COVID-19 Variants on Binding of the S Protein to ACE2. Biomolecules 2022; 12:biom12040572. [PMID: 35454161 PMCID: PMC9030943 DOI: 10.3390/biom12040572] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Revised: 03/06/2022] [Accepted: 03/07/2022] [Indexed: 12/29/2022] Open
Abstract
COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has many variants that accelerated the spread of the virus. In this study, we investigated the quantitative effect of some major mutants of the spike protein of SARS-CoV-2 binding to the human angiotensin-converting enzyme 2 (ACE2). These mutations are directly related to the Variant of Concern (VOC) including Alpha, Beta, Gamma, Delta and Omicron. Our calculations show that five major mutations (N501Y, E484K, L452R, T478K and K417N), first reported in Alpha, Beta, Gamma and Delta variants, all increase the binding of the S protein to ACE2 (except K417N), consistent with the experimental findings. We also studied an additional eight mutations of the Omicron variant that are located on the interface of the receptor binding domain (RDB) and have not been reported in other VOCs. Our study showed that most of these mutations (except Y505H and G446S) enhance the binding of the S protein to ACE2. The computational predictions helped explain why the Omicron variant quickly became dominant worldwide. Finally, comparison of several different computational methods for binding free energy calculation of these mutants was made. The alanine scanning method used in the current calculation helped to elucidate the residue-specific interactions responsible for the enhanced binding affinities of the mutants. The results show that the ASGB (alanine scanning with generalized Born) method is an efficient and reliable method for these binding free energy calculations due to mutations.
Collapse
Affiliation(s)
- Zhendong Li
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
| | - John Z. H. Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China;
- CAS Key Laboratory of Quantitative Engineering Biology, Shenzhen Institute of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
- Faculty of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
- NYU-ECNU Center for Computational Chemistry at New York University Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York, NY 10003, USA
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan 030006, China
- Correspondence:
| |
Collapse
|
38
|
Tang R, Chen P, Wang Z, Wang L, Hao H, Hou T, Sun H. Characterizing the stabilization effects of stabilizers in protein-protein systems with end-point binding free energy calculations. Brief Bioinform 2022; 23:6565618. [PMID: 35395683 DOI: 10.1093/bib/bbac127] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Revised: 02/10/2022] [Accepted: 03/15/2022] [Indexed: 02/06/2023] Open
Abstract
Drug design targeting protein-protein interactions (PPIs) associated with the development of diseases has been one of the most important therapeutic strategies. Besides interrupting the PPIs with PPI inhibitors/blockers, increasing evidence shows that stabilizing the interaction between two interacting proteins may also benefit the therapy, such as the development of various types of molecular glues/stabilizers that mostly work by stabilizing the two interacting proteins to regulate the downstream biological effects. However, characterizing the stabilization effect of a stabilizer is usually hard or too complicated for traditional experiments since it involves ternary interactions [protein-protein-stabilizer (PPS) interaction]. Thus, developing reliable computational strategies will facilitate the discovery/design of molecular glues or PPI stabilizers. Here, by fully analyzing the energetic features of the binary interactions in the PPS ternary complex, we systematically investigated the performance of molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) methods on characterizing the stabilization effects of stabilizers in 14-3-3 systems. The results show that both MM/PBSA and MM/GBSA are powerful tools in distinguishing the stabilizers from the decoys (with area under the curves of 0.90-0.93 for all tested cases) and are reasonable for ranking protein-peptide interactions in the presence or absence of stabilizers as well (with the average Pearson correlation coefficient of ~0.6 at a relatively high dielectric constant for both methods). Moreover, to give a detailed picture of the stabilization effects, the stabilization mechanism is also analyzed from the structural and energetic points of view for individual systems containing strong or weak stabilizers. This study demonstrates a potential strategy to accelerate the discovery of PPI stabilizers.
Collapse
Affiliation(s)
- Rongfan Tang
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009, Jiangsu, P. R. China
| | - Pengcheng Chen
- Institute of Big Data and Artificial Intelligence in Medicine, School of Electronics and Information Engineering, Taizhou University, Taizhou 318000, Zhejiang, P. R. China
| | - Zhe Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, P. R. China
| | - Lingling Wang
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009, Jiangsu, P. R. China
| | - Haiping Hao
- State Key Laboratory of Natural Medicines, Key Lab of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, 210009 Nanjing, China
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, Zhejiang, P. R. China
| | - Huiyong Sun
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009, Jiangsu, P. R. China
| |
Collapse
|
39
|
Panday S, Alexov E. Protein-Protein Binding Free Energy Predictions with the MM/PBSA Approach Complemented with the Gaussian-Based Method for Entropy Estimation. ACS OMEGA 2022; 7:11057-11067. [PMID: 35415339 PMCID: PMC8991903 DOI: 10.1021/acsomega.1c07037] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Accepted: 03/10/2022] [Indexed: 06/14/2023]
Abstract
Here, we present a Gaussian-based method for estimation of protein-protein binding entropy to augment the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) method for computational prediction of binding free energy (ΔG). The method is termed f5-MM/PBSA/E, where "E" stands for entropy and f5 for five adjustable parameters. The enthalpy components of ΔG (molecular mechanics, polar and non-polar solvation energies) are computed from a single implicit solvent generalized Born (GB) energy minimized structure of a protein-protein complex, while the binding entropy is computed using independently GB energy minimized unbound and bound structures. It should be emphasized that the f5-MM/PBSA/E method does not use snapshots, just energy minimized structures, and is thus very fast and computationally efficient. The method is trained and benchmarked in 5-fold validation test over a data set consisting of 46 protein-protein binding cases with experimentally determined dissociation constant K d values. This data set has been used for benchmarking in recently published protein-protein binding studies that apply conventional MM/PBSA and MM/PBSA with an enhanced sampling method. The f5-MM/PBSA/E tested on the same data set achieves similar or better performance than these computationally demanding approaches, making it an excellent choice for high throughput protein-protein binding affinity prediction studies.
Collapse
|
40
|
Mosquera-Yuqui F, Lopez-Guerra N, Moncayo-Palacio EA. Targeting the 3CLpro and RdRp of SARS-CoV-2 with phytochemicals from medicinal plants of the Andean Region: molecular docking and molecular dynamics simulations. J Biomol Struct Dyn 2022; 40:2010-2023. [PMID: 33084512 PMCID: PMC7657401 DOI: 10.1080/07391102.2020.1835716] [Citation(s) in RCA: 39] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Accepted: 10/05/2020] [Indexed: 02/08/2023]
Abstract
Given the highly contagious nature of SARS-CoV-2, it has resulted in an unprecedented number of COVID-19 infected and dead people worldwide. Since there is currently no vaccine available in the market, the identification of potential drugs is urgently needed to control the pandemic. In this study, 92 phytochemicals from medicinal plants growing in the Andean region were screened against SARS-CoV-2 3 C-like protease (3CLpro) and RNA-dependent RNA polymerase (RdRp) in their active sites through molecular docking. The cutoff values were set from the lowest docking scores of the FDA-approved drugs that are being used to treat COVID-19 patients (remdesivir, lopinavir, and ritonavir). Compounds with docking scores that were lower than cutoff values were validated by molecular dynamics simulation with GROMACS, using root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), and intermolecular hydrogen bonds (H-bonds). Furthermore, binding free energies were estimated using the MM-PBSA method, and ADMET profiles of potential inhibitors were assessed. Computational analyses revealed that the interaction with hesperidin (theoretical binding energies, ΔGbind = -15.18 kcal/mol to 3CLpro and ΔGbind = -9.46 kcal/mol to RdRp) remained stable in both enzymes, unveiling its remarkable potential as a possible multitarget antiviral agent to treat COVID-19. Importantly, lupinifolin with an estimated binding affinity to 3CLpro higher than hesperidin (ΔGbind = -20.93 kcal/mol) is also a potential inhibitor of the 3CLpro. These two compounds displayed suitable pharmacological and structural properties to be drug candidates, demonstrating to be worthy of further research.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Francisco Mosquera-Yuqui
- Departamento de Ciencias de La Vida y La Agricultura, Universidad de Las Fuerzas Armadas-ESPE, Sangolquí, Ecuador
- Grupo de Investigación y Desarrollo de la Biotecnología BioSin-Biociencias, Quito, Ecuador
| | - Nicolas Lopez-Guerra
- Grupo de Investigación y Desarrollo de la Biotecnología BioSin-Biociencias, Quito, Ecuador
- Oklahoma State University, Stillwater, OK, USA
| | - Eduardo A. Moncayo-Palacio
- Grupo de Investigación y Desarrollo de la Biotecnología BioSin-Biociencias, Quito, Ecuador
- Departamento de Química y Ciencias Exactas, Universidad Técnica Particular de Loja, Loja, Ecuador
| |
Collapse
|
41
|
|
42
|
Pang JP, Hu XP, Wang YX, Liao JN, Chai X, Wang XW, Shen C, Wang JJ, Zhang LL, Wang XY, Zhu F, Weng QJ, Xu L, Hou TJ, Li D. Discovery of a novel nonsteroidal selective glucocorticoid receptor modulator by virtual screening and bioassays. Acta Pharmacol Sin 2022; 43:2429-2438. [PMID: 35110698 PMCID: PMC8809242 DOI: 10.1038/s41401-021-00855-6] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Accepted: 12/27/2021] [Indexed: 12/31/2022] Open
Abstract
Synthetic glucocorticoids (GCs) have been widely used in the treatment of a broad range of inflammatory diseases, but their clinic use is limited by undesired side effects such as metabolic disorders, osteoporosis, skin and muscle atrophies, mood disorders and hypothalamic-pituitary-adrenal (HPA) axis suppression. Selective glucocorticoid receptor modulators (SGRMs) are expected to have promising anti-inflammatory efficacy but with fewer side effects caused by GCs. Here, we reported HT-15, a prospective SGRM discovered by structure-based virtual screening (VS) and bioassays. HT-15 can selectively act on the NF-κB/AP1-mediated transrepression function of glucocorticoid receptor (GR) and repress the expression of pro-inflammation cytokines (i.e., IL-1β, IL-6, COX-2, and CCL-2) as effectively as dexamethasone (Dex). Compared with Dex, HT-15 shows less transactivation potency that is associated with the main adverse effects of synthetic GCs, and no cross activities with other nuclear receptors. Furthermore, HT-15 exhibits very weak inhibition on the ratio of OPG/RANKL. Therefore, it may reduce the side effects induced by normal GCs. The bioactive compound HT-15 can serve as a starting point for the development of novel therapeutics for high dose or long-term anti-inflammatory treatment.
Collapse
|
43
|
Zhai J, He X, Man VH, Sun Y, Ji B, Cai L, Wang J. A multiple-step in silico screening protocol to identify allosteric inhibitors of Spike-hACE2 binding. Phys Chem Chem Phys 2022; 24:4305-4316. [PMID: 35107459 DOI: 10.1039/d1cp04736a] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
While the COVID-19 pandemic continues to worsen, effective medicines that target the life cycle of SARS-CoV-2 are still under development. As more highly infective and dangerous variants of the coronavirus emerge, the protective power of vaccines will decrease or vanish. Thus, the development of drugs, which are free of drug resistance is direly needed. The aim of this study is to identify allosteric binding modulators from a large compound library to inhibit the binding between the Spike protein of the SARS-CoV-2 virus and human angiotensin-converting enzyme 2 (hACE2). The binding of the Spike protein to hACE2 is the first step of the infection of host cells by the coronavirus. We first built a compound library containing 77 448 antiviral compounds. Molecular docking was then conducted to preliminarily screen compounds which can potently bind to the Spike protein at two allosteric binding sites. Next, molecular dynamics simulations were performed to accurately calculate the binding affinity between the spike protein and an identified compound from docking screening and to investigate whether the compound can interfere with the binding between the Spike protein and hACE2. We successfully identified two possible drug binding sites on the Spike protein and discovered a series of antiviral compounds which can weaken the interaction between the Spike protein and hACE2 receptor through conformational changes of the key Spike residues at the Spike-hACE2 binding interface induced by the binding of the ligand at the allosteric binding site. We also applied our screening protocol to another compound library which consists of 3407 compounds for which the inhibitory activities of Spike/hACE2 binding were measured. Encouragingly, in vitro data supports that the identified compounds can inhibit the Spike-ACE2 binding. Thus, we developed a promising computational protocol to discover allosteric inhibitors of the binding of the Spike protein of SARS-CoV-2 to the hACE2 receptor, and several promising allosteric modulators were discovered.
Collapse
Affiliation(s)
- Jingchen Zhai
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Yuchen Sun
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Beihong Ji
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Lianjin Cai
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| |
Collapse
|
44
|
Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex. Molecules 2022; 27:molecules27030799. [PMID: 35164065 PMCID: PMC8839245 DOI: 10.3390/molecules27030799] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 01/20/2022] [Accepted: 01/21/2022] [Indexed: 12/02/2022] Open
Abstract
The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.
Collapse
|
45
|
Gundelach L, Fox T, Tautermann CS, Skylaris CK. BRD4: quantum mechanical protein–ligand binding free energies using the full-protein DFT-based QM-PBSA method. Phys Chem Chem Phys 2022; 24:25240-25249. [PMID: 36222107 DOI: 10.1039/d2cp03705j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Fully quantum mechanical approaches to calculating protein–ligand free energies of binding have the potential to reduce empiricism and explicitly account for all physical interactions responsible for protein–ligand binding.
Collapse
Affiliation(s)
- Lennart Gundelach
- University of Southampton, Faculty of Engineering Science and Mathematics, Chemistry, University Road, Southampton, SO17 1BJ, UK
| | - Thomas Fox
- Boehringer Ingelheim Pharma GmbH & Co KG, Medicinal Chemistry, Birkendorfer Str 65, 88397, Biberach, Germany
| | - Christofer S. Tautermann
- Boehringer Ingelheim Pharma GmbH & Co KG, Medicinal Chemistry, Birkendorfer Str 65, 88397, Biberach, Germany
| | - Chris-Kriton Skylaris
- University of Southampton, Faculty of Engineering Science and Mathematics, Chemistry, University Road, Southampton, SO17 1BJ, UK
| |
Collapse
|
46
|
Havranek B, Chan KK, Wu A, Procko E, Islam SM. Computationally Designed ACE2 Decoy Receptor Binds SARS-CoV-2 Spike (S) Protein with Tight Nanomolar Affinity. J Chem Inf Model 2021; 61:4656-4669. [PMID: 34427448 PMCID: PMC8409145 DOI: 10.1021/acs.jcim.1c00783] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2021] [Indexed: 12/25/2022]
Abstract
Even with the availability of vaccines, therapeutic options for COVID-19 still remain highly desirable, especially in hospitalized patients with moderate or severe disease. Soluble ACE2 (sACE2) is a promising therapeutic candidate that neutralizes SARS CoV-2 infection by acting as a decoy. Using computational mutagenesis, we designed a number of sACE2 derivatives carrying three to four mutations. The top-predicted sACE2 decoy based on the in silico mutagenesis scan was subjected to molecular dynamics and free-energy calculations for further validation. After illuminating the mechanism of increased binding for our designed sACE2 derivative, the design was verified experimentally by flow cytometry and BLI-binding experiments. The computationally designed sACE2 decoy (ACE2-FFWF) bound the receptor-binding domain of SARS-CoV-2 tightly with low nanomolar affinity and ninefold affinity enhancement over the wild type. Furthermore, cell surface expression was slightly greater than wild-type ACE2, suggesting that the design is well-folded and stable. Having an arsenal of high-affinity sACE2 derivatives will help to buffer against the emergence of SARS CoV-2 variants. Here, we show that computational methods have become sufficiently accurate for the design of therapeutics for current and future viral pandemics.
Collapse
Affiliation(s)
- Brandon Havranek
- Department of Chemistry, University of
Illinois at Chicago, Chicago, Illinois 60607, United
States
| | - Kui K. Chan
- Orthogonal Biologics Inc.,
Urbana, Illinois 61801, United States
| | - Austin Wu
- Department of Computer Science,
Northwestern University, Evanston, Illinois 60208,
United States
| | - Erik Procko
- Department of Biochemistry and Cancer Center at
Illinois, University of Illinois, Urbana, Illinois 61801,
United States
| | - Shahidul M. Islam
- Department of Chemistry, University of
Illinois at Chicago, Chicago, Illinois 60607, United
States
| |
Collapse
|
47
|
Discovery of potential inhibitors targeting the kinase domain of polynucleotide kinase/phosphatase (PNKP): Homology modeling, virtual screening based on multiple conformations, and molecular dynamics simulation. Comput Biol Chem 2021; 94:107517. [PMID: 34456161 DOI: 10.1016/j.compbiolchem.2021.107517] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2021] [Revised: 05/11/2021] [Accepted: 05/16/2021] [Indexed: 12/15/2022]
Abstract
In recent years, the level of interest has been increased in developing the DNA-repair inhibitors, to enhance the cytotoxic effects in the treatment of cancers. Polynucleotide kinase/phosphatase (PNKP) is a critical human DNA repair enzyme that repairs DNA strand breaks by catalyzing the restoration of 5'-phosphate and 3'-hydroxyl termini that are required for subsequent processing by DNA ligases and polymerases. PNKP is the only protein that repairs the 3'-hydroxyl group and 5'-phosphate group, which depicts PNKP as a potential therapeutic target. Besides, PNKP is the only DNA-repair enzyme that contains the 5'-kinase activity, therefore, targeting this kinase domain would motivate the development of novel PNKP-specific inhibitors. However, there are neither crystal structures of human PNKP nor the kinase inhibitors reported so far. Thus, in this present study, a sequential molecular docking-based virtual screening with multiple PNKP conformations integrating homology modeling, molecular dynamics simulation, and binding free energy calculation was developed to discover novel PNKP kinase inhibitors, and the top-scored molecule was finally submitted to molecular dynamics simulation to reveal the binding mechanism between the inhibitor and PNKP. Taken together, the current study could provide some guidance for the molecular docking based-virtual screening of novel PNKP kinase inhibitors.
Collapse
|
48
|
King E, Aitchison E, Li H, Luo R. Recent Developments in Free Energy Calculations for Drug Discovery. Front Mol Biosci 2021; 8:712085. [PMID: 34458321 PMCID: PMC8387144 DOI: 10.3389/fmolb.2021.712085] [Citation(s) in RCA: 49] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 07/27/2021] [Indexed: 01/11/2023] Open
Abstract
The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.
Collapse
Affiliation(s)
- Edward King
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Erick Aitchison
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Han Li
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
| | - Ray Luo
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
- Department of Materials Science and Engineering, University of California, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, CA, United States
| |
Collapse
|
49
|
Li H, Yang W, Wang Z, Wang X, Hao Y, Xi J, Lu H, Du Z, Feng J, Zhang B, Ma D. Computational research of mTORC1 inhibitor on cerebral ischemia-reperfusion injury. Aging (Albany NY) 2021; 13:19598-19613. [PMID: 34343111 PMCID: PMC8386574 DOI: 10.18632/aging.203371] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2021] [Accepted: 07/08/2021] [Indexed: 12/16/2022]
Abstract
Ischemic stroke contributes to more than 80% of all strokes and has the four characteristics of high prevalence, high disability, high mortality, and high recurrence. Stroke is a preventable and controllable disease. In addition to controlling the primary disease, effective prevention and control measures need to be given to the occurrence and development of stroke. With the development and progress of modern treatment methods for ischemic stroke, the mortality and disability rate have decreased significantly. At present, the main treatment methods for ischemic stroke include thrombolysis, thrombus removal at the ultra-early stage, and treatment of improving collateral circulation in the acute phase. However, the ultra-early and early blood reperfusion involves reperfusion injury, which will cause secondary nerve damage, which is called cerebral ischemia/reperfusion injury (CIRI). Studies have found that autophagy is involved in the entire process of CIRI and can reduce the damage of CIRI. The mammalian target of Rapamycin (mTORC1) is the primary signal pathway regulating autophagy. And the mTORC1 inhibitor, Rapamycin, has been proved to exert neuroprotective effects in the ultra-early and early cerebral ischemia-reperfusion. Therefore, screening and designing mTORC1 inhibitors is very important to control reperfusion injury and reduce neuronal death and apoptosis. In this research, plenty of computer-assisted was applied to virtually screen and select potential mTORC1's inhibitors. We used Libdock to screen the structure and performed toxicity predictions, ADME (absorption, distribution, metabolism, excretion) to predict small molecules' pharmacological and toxicological properties. To assess the binding mechanism and affinity between the mTORC1 dimer and the ligand, molecular docking was performed. Then, the pharmacophore of small molecules in the docking conformation with the protein was supplemented by Schrodinger. Additionally, molecular dynamics simulations were conducted to assess if the ligand-receptor complex was stable in a natural environment. Furthermore, an experiment was performed to verify the inhibitory effect of compound 1 and compound 2 on mTOR protein. All in all, the study provides a hand of candidate drugs as well as pharmacological properties, which can play an essential role in mTORC1 inhibitors.
Collapse
Affiliation(s)
- Hui Li
- Department of Neurology and Neuroscience Center, The First Hospital of Jilin University, Changchun, China
| | - Wenzhuo Yang
- Clinical College, Jilin University, Changchun, China
| | - Zhenhua Wang
- Clinical College, Jilin University, Changchun, China
| | - Xu Wang
- Department of Neurology and Neuroscience Center, The First Hospital of Jilin University, Changchun, China
| | - Yulei Hao
- Department of Neurology and Neuroscience Center, The First Hospital of Jilin University, Changchun, China
| | - Jianxin Xi
- Clinical College, Jilin University, Changchun, China
| | - Han Lu
- Clinical College, Jilin University, Changchun, China
| | - Zhishan Du
- Clinical College, Jilin University, Changchun, China
| | - Jiachun Feng
- Department of Neurology and Neuroscience Center, The First Hospital of Jilin University, Changchun, China
| | - Bao Zhang
- Department of Urology Surgery, Aerospace Center Hospital, Beijing, China
| | - Di Ma
- Department of Neurology and Neuroscience Center, The First Hospital of Jilin University, Changchun, China
| |
Collapse
|
50
|
Chitranshi N, Kumar A, Sheriff S, Gupta V, Godinez A, Saks D, Sarkar S, Shen T, Mirzaei M, Basavarajappa D, Abyadeh M, Singh SK, Dua K, Zhang KYJ, Graham SL, Gupta V. Identification of Novel Cathepsin B Inhibitors with Implications in Alzheimer's Disease: Computational Refining and Biochemical Evaluation. Cells 2021; 10:cells10081946. [PMID: 34440715 PMCID: PMC8391575 DOI: 10.3390/cells10081946] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2021] [Revised: 07/23/2021] [Accepted: 07/27/2021] [Indexed: 02/07/2023] Open
Abstract
Amyloid precursor protein (APP), upon proteolytic degradation, forms aggregates of amyloid β (Aβ) and plaques in the brain, which are pathological hallmarks of Alzheimer’s disease (AD). Cathepsin B is a cysteine protease enzyme that catalyzes the proteolytic degradation of APP in the brain. Thus, cathepsin B inhibition is a crucial therapeutic aspect for the discovery of new anti-Alzheimer’s drugs. In this study, we have employed mixed-feature ligand-based virtual screening (LBVS) by integrating pharmacophore mapping, docking, and molecular dynamics to detect small, potent molecules that act as cathepsin B inhibitors. The LBVS model was generated by using hydrophobic (HY), hydrogen bond acceptor (HBA), and hydrogen bond donor (HBD) features, using a dataset of 24 known cathepsin B inhibitors of both natural and synthetic origins. A validated eight-feature pharmacophore hypothesis (Hypo III) was utilized to screen the Maybridge chemical database. The docking score, MM-PBSA, and MM-GBSA methodology was applied to prioritize the lead compounds as virtual screening hits. These compounds share a common amide scaffold, and showed important interactions with Gln23, Cys29, His110, His111, Glu122, His199, and Trp221. The identified inhibitors were further evaluated for cathepsin-B-inhibitory activity. Our study suggests that pyridine, acetamide, and benzohydrazide compounds could be used as a starting point for the development of novel therapeutics.
Collapse
Affiliation(s)
- Nitin Chitranshi
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
- Correspondence: (N.C.); (V.G.); Tel.: +61-(02)-9850-2804 (N.C.)
| | - Ashutosh Kumar
- Center for Biosystems Dynamics Research, Laboratory for Structural Bioinformatics, RIKEN, 1-7-22 Suehiro, Tsurumi, Yokohama 230-0045, Kanagawa, Japan; (A.K.); (K.Y.J.Z.)
| | - Samran Sheriff
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Veer Gupta
- School of Medicine, Faculty of Health, Deakin University, Geelong, VIC 3220, Australia;
| | - Angela Godinez
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Danit Saks
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Soumalya Sarkar
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Ting Shen
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Mehdi Mirzaei
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Devaraj Basavarajappa
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Morteza Abyadeh
- Cell Science Research Center, Department of Molecular Systems Biology, Royan Institute for Stem Cell Biology and Technology, ACECR, Tehran 1665659911, Iran;
| | - Sachin K. Singh
- School of Pharmaceutical Sciences, Lovely Professional University, Phagwara 144411, Punjab, India;
| | - Kamal Dua
- Discipline of Pharmacy, Graduate School of Health, University of Technology Sydney, Ultimo, NSW 2007, Australia;
- Australian Research Centre in Complementary and Integrative Medicine, Faculty of Health, University of Technology Sydney, Ultimo, NSW 2007, Australia
| | - Kam Y. J. Zhang
- Center for Biosystems Dynamics Research, Laboratory for Structural Bioinformatics, RIKEN, 1-7-22 Suehiro, Tsurumi, Yokohama 230-0045, Kanagawa, Japan; (A.K.); (K.Y.J.Z.)
| | - Stuart L. Graham
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
| | - Vivek Gupta
- Faculty of Medicine, Health and Human Sciences, Macquarie University, F10A, 2 Technology Place, North Ryde, NSW 2109, Australia; (S.S.); (A.G.); (D.S.); (S.S.); (T.S.); (M.M.); (D.B.); (S.L.G.)
- Correspondence: (N.C.); (V.G.); Tel.: +61-(02)-9850-2804 (N.C.)
| |
Collapse
|