1
|
Zhou RN, Zhu ZW, Xu PY, Shen LX, Wang Z, Xue YY, Xiang YY, Cao Y, Yu XZ, Zhao J, Jin Y, Yan J, Yang Q, Fang PH, Shang WB. Rhein targets macrophage SIRT2 to promote adipose tissue thermogenesis in obesity in mice. Commun Biol 2024; 7:1003. [PMID: 39152196 PMCID: PMC11329635 DOI: 10.1038/s42003-024-06693-6] [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: 01/22/2024] [Accepted: 08/06/2024] [Indexed: 08/19/2024] Open
Abstract
Rhein, a component derived from rhubarb, has been proven to possess anti-inflammatory properties. Here, we show that rhein mitigates obesity by promoting adipose tissue thermogenesis in diet-induced obese mice. We construct a macrophage-adipocyte co-culture system and demonstrate that rhein promotes adipocyte thermogenesis through inhibiting NLRP3 inflammasome activation in macrophages. Moreover, clues from acetylome analysis identify SIRT2 as a potential drug target of rhein. We further verify that rhein directly interacts with SIRT2 and inhibits NLRP3 inflammasome activation in a SIRT2-dependent way. Myeloid knockdown of SIRT2 abrogates adipose tissue thermogenesis and metabolic benefits in obese mice induced by rhein. Together, our findings elucidate that rhein inhibits NLRP3 inflammasome activation in macrophages by regulating SIRT2, and thus promotes white adipose tissue thermogenesis during obesity. These findings uncover the molecular mechanism underlying the anti-inflammatory and anti-obesity effects of rhein, and suggest that rhein may become a potential drug for treating obesity.
Collapse
Affiliation(s)
- Ruo-Nan Zhou
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Zi-Wei Zhu
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Ping-Yuan Xu
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Li-Xuan Shen
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Ziwei Wang
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Ying-Ying Xue
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Ying-Ying Xiang
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Yue Cao
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Xi-Zhong Yu
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Juan Zhao
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Yu Jin
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Jing Yan
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China
| | - Qin Yang
- Department of Medicine, Physiology and Biophysics, University of California, Irvine, CA, 92697, USA
| | - Peng-Hua Fang
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China.
| | - Wen-Bin Shang
- Department of Endocrinology, Jiangsu Province Hospital of Chinese Medicine, the Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, 210029, China.
- Key Laboratory for Metabolic Diseases in Chinese Medicine, Nanjing University of Chinese Medicine, Nanjing, 210023, China.
| |
Collapse
|
2
|
Fan C, Basharat Z, Mah K, Wei CR. Computational approach for drug discovery against Gardnerella vaginalis in quest for safer and effective treatments for bacterial vaginosis. Sci Rep 2024; 14:17437. [PMID: 39075099 PMCID: PMC11286753 DOI: 10.1038/s41598-024-68443-2] [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: 05/20/2024] [Accepted: 07/23/2024] [Indexed: 07/31/2024] Open
Abstract
Bacterial vaginosis (BV), primarily attributed to Gardnerella vaginalis, poses significant challenges due to antibiotic resistance and suboptimal treatment outcomes. This study presents an integrated approach to identify potential drug targets and screen compounds against this bacterium by leveraging a computational methodology. Subtractive proteomics of the reference strain ASM286196v1/UMB0386 (assembly accession: GCA_002861965.1) facilitated the prioritization of proteins with essential bacterial functions and pathways as potential drug targets. We selected 3-deoxy-7-phosphoheptulonate synthase (aroG gene product; also known as DAHP synthase) for downstream analysis. Molecular docking was employed in PyRx (AutoDock Vina) to predict binding affinities between aroG inhibitors from the ZINC database and 3-deoxy-7-phosphoheptulonate synthase. Molecular dynamics simulations of 100 ns, using GROMACS, validated the stability of drug-target interactions. Additionally, ADMET profiling aided in the selection of compounds with favorable pharmacokinetic properties and safety profile for human hosts. PBPK profiling showed that ZINC98088375 had the highest bioavailability and efficient systemic circulation. Conversely, ZINC5113880 demonstrated the lowest absorption rate (39.661%). Moreover, cirrhosis, steatosis, and renal impairment appeared to influence blood concentration of the drug, impacting bioavailability. The integrative -omics approach utilized in this study underscores the potential of computer-aided drug design and offers a rational strategy for targeted inhibitor discovery against G. vaginalis. The strategy is an attempt to address the limitations of current BV treatments, including antibiotic resistance, and pave way for the development of safer and more effective therapeutics.
Collapse
Affiliation(s)
- Chenyue Fan
- Department of Research and Development, Shing Huei Group, Taipei, 10617, Taiwan
- College of Pharmacy, University of Arizona, Tuscon, AZ, 85721, USA
| | | | - Karmen Mah
- Department of Research and Development, Shing Huei Group, Taipei, 10617, Taiwan
| | - Calvin R Wei
- Department of Research and Development, Shing Huei Group, Taipei, 10617, Taiwan.
| |
Collapse
|
3
|
Pan Y, Zhao C, Fu W, Yang S, Lv S. Comparative analysis of structural dynamics and allosteric mechanisms of RecA/Rad51 family proteins: Integrated atomistic MD simulation and network-based analysis. Int J Biol Macromol 2024; 261:129843. [PMID: 38302027 DOI: 10.1016/j.ijbiomac.2024.129843] [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: 08/17/2023] [Revised: 01/26/2024] [Accepted: 01/27/2024] [Indexed: 02/03/2024]
Abstract
Homologous recombination plays a key role in double-strand break repair, stalled replication fork repair, and meiosis. The RecA/Rad51 family recombinases catalyze the DNA strand invasion reaction that occurs during homologous recombination. However, the high sequence differences between homologous groups have hindered the thoroughly studies of this ancient protein family. The dynamic mechanisms of the family, particularly at the residual level, remain poorly understood. In this work, five representative RecA/Rad51 recombinase family members from all major kingdoms of living organisms: prokaryotes, eukaryotes, archaea, and viruses, were selected to explore the molecular mechanisms behind their conserved biological significance. A variety of techniques, including all-atom molecular dynamics simulation, perturbation response scanning, and protein structure network analysis, were used to examine the flexibility and correlation of protein domains, distribution of sensors and effectors and conserved hub residues. Furthermore, the potential communication routes between the ATP-binding region and the DNA-binding region of each recombinase were identified. Our results demonstrate the conserved molecular dynamics of these recombinases in the early stage of homologous recombination, including cooperative motions between regions, conserved sensing and effecting functional residue distribution, and conserved hub residues. Meanwhile, the unique ATP-DNA communication routes of each recombinase was also revealed. These results provide new insights into the mechanism of RecA/Rad51 family proteins, and provide new theoretical guidance for the development of allosteric inhibitors and the application of RecA/Rad51 family proteins.
Collapse
Affiliation(s)
- Yue Pan
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Chong Zhao
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Wenyu Fu
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Shuo Yang
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China.
| | - Shaowu Lv
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China; Bioarchaeology Laboratory, Jilin University, 2699 Qianjin Street, Changchun 130012, China.
| |
Collapse
|
4
|
Khan H, Waqas M, Khurshid B, Ullah N, Khalid A, Abdalla AN, Alamri MA, Wadood A. Investigating the role of Sterol C24-Methyl transferase mutation on drug resistance in leishmaniasis and identifying potential inhibitors. J Biomol Struct Dyn 2023:1-14. [PMID: 37723868 DOI: 10.1080/07391102.2023.2256879] [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: 12/19/2022] [Accepted: 09/02/2023] [Indexed: 09/20/2023]
Abstract
Leishmaniasis is a fatal disease caused by the leishmania parasite. For the survival of the leishmania parasite, Sterol C24-Methyl Transferase (SMT) is essential which is an enzyme of the ergosterol pathway. SMT protein mutation is responsible for Amphotericin-B drug resistance in Leishmania, which is the main treatment for visceral leishmaniasis. Amphotericin-B resistance is caused by three mutated residues V131I, V321I and F72C. The underlying mechanisms and structural changes in SMT enzymes responsible for resistance due to mutation are still not well understood. In the current study, the potential mechanism of resistance due to these mutations and the structure variation of wild and mutant SMT proteins were investigated through molecular dynamics simulations and molecular docking analysis. The results showed that AmB established strong bonding interaction with wild SMT as compare to mutants SMT. The binding energy calculation showed that binding energy of AmB with mutants SMT increases as compare to the wild SMT. Further structural based virtual screening was carried out to design potential inhibitors for the mutant SMT. On the basis of structural-based virtual screening four inhibitors (SANC01057, SANC00882, SANC00414, SANC01047) were computationally identified as potential mutant SMT (F72C) inhibitors. This work provides valuable information for improved management of drug resistant Leishmaniasis.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Huma Khan
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Muhammad Waqas
- Natural and Medical Sciences Research Center, University of Nizwa, Birkat Al-Mouz Nizwa, Oman
| | - Beenish Khurshid
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Nazif Ullah
- Department of Biotechnology, Abdul Wali Khan University Mardan, Mardan, Pakistan
| | - Asaad Khalid
- Substance Abuse and Toxicology Research Center, Jazan University, Jazan, Saudi Arabia
| | - Ashraf N Abdalla
- Department of Pharmacology and Toxicology, College of Pharmacy, Umm Al-Qura University, Makkah, Saudi Arabia
| | - Mubarak A Alamri
- Department of Pharmaceutical Chemistry, College of Pharmacy, Prince Sattam Bin Abdulaziz University, Al-Kharj, Saudi Arabia
| | - Abdul Wadood
- Department of Biochemistry, Abdul Wali Khan University Mardan, Mardan, Pakistan
| |
Collapse
|
5
|
Bassani D, Moro S. Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies. Molecules 2023; 28:molecules28093906. [PMID: 37175316 PMCID: PMC10180087 DOI: 10.3390/molecules28093906] [Citation(s) in RCA: 24] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Revised: 04/28/2023] [Accepted: 05/02/2023] [Indexed: 05/15/2023] Open
Abstract
The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.
Collapse
Affiliation(s)
- Davide Bassani
- Pharmaceutical Research & Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., 4070 Basel, Switzerland
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131 Padova, Italy
| | - Stefano Moro
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences, University of Padova, Via Marzolo 5, 35131 Padova, Italy
| |
Collapse
|
6
|
Pan Y, Xie N, Zhang X, Yang S, Lv S. Computational Insights into the Dynamic Structural Features and Binding Characteristics of Recombinase UvsX Compared with RecA. Molecules 2023; 28:molecules28083363. [PMID: 37110596 PMCID: PMC10144138 DOI: 10.3390/molecules28083363] [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: 03/03/2023] [Revised: 04/06/2023] [Accepted: 04/07/2023] [Indexed: 04/29/2023] Open
Abstract
RecA family recombinases are the core enzymes in the process of homologous recombination, and their normal operation ensures the stability of the genome and the healthy development of organisms. The UvsX protein from bacteriophage T4 is a member of the RecA family recombinases and plays a central role in T4 phage DNA repair and replication, which provides an important model for the biochemistry and genetics of DNA metabolism. UvsX shares a high degree of structural similarity and function with RecA, which is the most deeply studied member of the RecA family. However, the detailed molecular mechanism of UvsX has not been resolved. In this study, a comprehensive all-atom molecular dynamics simulation of the UvsX protein dimer complex was carried out in order to investigate the conformational and binding properties of UvsX in combination with ATP and DNA, and the simulation of RecA was synchronized with the property comparison learning for UvsX. This study confirmed the highly conserved molecular structure characteristics and catalytic centers of RecA and UvsX, and also discovered differences in regional conformation, volatility and the ability to bind DNA between the two proteins at different temperatures, which would be helpful for the subsequent understanding and application of related recombinases.
Collapse
Affiliation(s)
- Yue Pan
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Ningkang Xie
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Xin Zhang
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Shuo Yang
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| | - Shaowu Lv
- Key Laboratory for Molecular Enzymology and Engineering of the Ministry of Education, School of Life Science, Jilin University, 2699 Qianjin Street, Changchun 130012, China
- Bioarchaeology Laboratory, Jilin University, 2699 Qianjin Street, Changchun 130012, China
| |
Collapse
|
7
|
Sun Z, He Q, Gong Z, Kalhor P, Huai Z, Liu Z. A General Picture of Cucurbit[8]uril Host–Guest Binding: Recalibrating Bonded Interactions. Molecules 2023; 28:molecules28073124. [PMID: 37049887 PMCID: PMC10095826 DOI: 10.3390/molecules28073124] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Revised: 03/15/2023] [Accepted: 03/27/2023] [Indexed: 04/03/2023] Open
Abstract
Atomic-level understanding of the dynamic feature of host–guest interactions remains a central challenge in supramolecular chemistry. The remarkable guest binding behavior of the Cucurbiturils family of supramolecular containers makes them promising drug carriers. Among Cucurbit[n]urils, Cucurbit[8]uril (CB8) has an intermediate portal size and cavity volume. It can exploit almost all host–guest recognition motifs formed by this host family. In our previous work, an extensive computational investigation of the binding of seven commonly abused and structurally diverse drugs to the CB8 host was performed, and a general dynamic binding picture of CB8-guest interactions was obtained. Further, two widely used fixed-charge models for drug-like molecules were investigated and compared in great detail, aiming at providing guidelines in choosing an appropriate charge scheme in host-guest modelling. Iterative refitting of atomic charges leads to improved binding thermodynamics and the best root-mean-squared deviation from the experimental reference is 2.6 kcal/mol. In this work, we focus on a thorough evaluation of the remaining parts of classical force fields, i.e., the bonded interactions. The widely used general Amber force fields are assessed and refitted with generalized force-matching to improve the intra-molecular conformational preference, and thus the description of inter-molecular host–guest interactions. The interaction pattern and binding thermodynamics show a significant dependence on the modelling parameters. The refitted system-specific parameter set improves the consistency of the modelling results and the experimental reference significantly. Finally, combining the previous charge-scheme comparison and the current force-field refitting, we provide general guidelines for the theoretical modelling of host–guest binding.
Collapse
|
8
|
Liu X, Zheng L, Qin C, Zhang JZH, Sun Z. Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: I. Standard procedure. J Comput Aided Mol Des 2022; 36:735-752. [PMID: 36136209 DOI: 10.1007/s10822-022-00475-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Accepted: 09/06/2022] [Indexed: 10/14/2022]
Abstract
Despite the massive application of end-point free energy methods in protein-ligand and protein-protein interactions, computational understandings about their performance in relatively simple and prototypical host-guest systems are limited. In this work, we present a comprehensive benchmark calculation with standard end-point free energy techniques in a recent host-guest dataset containing 13 host-guest pairs involving the carboxylated-pillar[6]arene host. We first assess the charge schemes for solutes by comparing the charge-produced electrostatics with many ab initio references, in order to obtain a preliminary albeit detailed view of the charge quality. Then, we focus on four modelling details of end-point free energy calculations, including the docking procedure for the generation of initial condition, the charge scheme for host and guest molecules, the water model used in explicit-solvent sampling, and the end-point methods for free energy estimation. The binding thermodynamics obtained with different modelling schemes are compared with experimental references, and some practical guidelines on maximizing the performance of end-point methods in practical host-guest systems are summarized. Further, we compare our simulation outcome with predictions in the grand challenge and discuss further developments to improve the prediction quality of end-point free energy methods. Overall, unlike the widely acknowledged applicability in protein-ligand binding, the standard end-point calculations cannot produce useful outcomes in host-guest binding and thus are not recommended unless alterations are performed.
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
| | - Chu Qin
- 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, New York, NY, 10003, USA
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
9
|
Sun S, Huggins DJ. Assessing the effect of forcefield parameter sets on the accuracy of relative binding free energy calculations. Front Mol Biosci 2022; 9:972162. [PMID: 36225254 PMCID: PMC9549959 DOI: 10.3389/fmolb.2022.972162] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Accepted: 08/10/2022] [Indexed: 11/13/2022] Open
Abstract
Software for accurate prediction of protein-ligand binding affinity can be a key enabling tool for small molecule drug discovery. Free energy perturbation (FEP) is a computational technique that can be used to compute binding affinity differences between molecules in a congeneric series. It has shown promise in reliably generating accurate predictions and is now widely used in the pharmaceutical industry. However, the high computational cost and use of commercial software, together with the technical challenges to setup, run, and analyze the simulations, limits the usage of FEP. Here, we use an automated FEP workflow which uses the open-source OpenMM package. To enable effective application of FEP, we compared the performance of different water models, partial charge assignments, and AMBER protein forcefields in eight benchmark test cases previously assembled for FEP validation studies.
Collapse
Affiliation(s)
- Shan Sun
- Tri-Institutional Therapeutics Discovery Institute, New York, NY, United States
| | - David J. Huggins
- Tri-Institutional Therapeutics Discovery Institute, New York, NY, United States
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, NY, United States
| |
Collapse
|
10
|
Sun Z, Wang M, He Q, Liu Z. Molecular Modeling of Ionic Liquids: Force‐Field Validation and Thermodynamic Perspective from Large‐Scale Fast‐Growth Solvation Free Energy Calculations. ADVANCED THEORY AND SIMULATIONS 2022. [DOI: 10.1002/adts.202200274] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| | - Mao Wang
- NCS Testing Technology Co., Ltd. No. 13, Gaoliangqiao Xiejie Beijing 100081 China
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio‐Engineering Co., Ltd. North Century Avenue 333 Ningbo 315100 China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| |
Collapse
|
11
|
Zhang Q, Yang Y, Gong X, Zhao N, Zhang Y, Liu H. Thermodynamic integration combined with molecular dynamic simulations to explore the
cross‐resistance
mechanism of isoniazid and ethionamide. Proteins 2022; 90:1142-1151. [DOI: 10.1002/prot.26295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2021] [Revised: 12/25/2021] [Accepted: 12/30/2021] [Indexed: 11/10/2022]
Affiliation(s)
| | - Yuwei Yang
- School of Pharmacy Lanzhou University Lanzhou China
| | - Xiaoqing Gong
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry Lanzhou University Lanzhou China
| | - Nannan Zhao
- School of Pharmacy Lanzhou University Lanzhou China
| | - Yang Zhang
- Supercomputing Center of Lanzhou University Lanzhou China
| | | |
Collapse
|
12
|
Sun Z, Huai Z, He Q, Liu Z. A General Picture of Cucurbit[8]uril Host-Guest Binding. J Chem Inf Model 2021; 61:6107-6134. [PMID: 34818004 DOI: 10.1021/acs.jcim.1c01208] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Describing, understanding, and designing complex interaction networks within macromolecular systems remain challenging in modern chemical research. Host-guest systems, despite their relative simplicity in both the structural feature and interaction patterns, still pose problems in theoretical modeling. The barrel-shaped supramolecular container cucurbit[8]uril (CB8) shows promising functionalities in various areas, e.g., catalysis and molecular recognition. It can stably coordinate a series of structurally diverse guests with high affinities. In this work, we examine the binding of seven commonly abused drugs to the CB8 host, aiming at providing a general picture of CB8-guest binding. Extensive sampling of the configurational space of these host-guest systems is performed, and the binding pathway and interaction patterns of CB8-guest complexes are investigated. A thorough comparison of widely used fixed-charge models for drug-like molecules is presented. Iterative refitting of the atomic charges suggests significant conformation dependence of charge generation. The initial model generated at the original conformation could be inaccurate for new conformations explored during conformational search, and the newly fitted charge set improves the prediction-experiment correlation significantly. Our investigations of the configurational space of CB8-drug complexes suggest that the host-guest interactions are more complex than expected. Despite the structural simplicities of these molecules, the conformational fluctuations of the host and the guest molecules and orientations of functional groups lead to the existence of an ensemble of binding modes. The insights of the binding thermodynamics, performance of fixed-charge models, and binding patterns of the CB8-guest systems are useful for studying and elucidating the binding mechanism of other host-guest complexes.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Zhe Huai
- XtalPi-AI Research Center (XARC), 9F, Tower A, Dongsheng Building, No. 8, Zhongguancun East Road, Haidian District, Beijing 100083, P.R. China
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd., North Century Avenue 333, Ningbo 315100, China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
13
|
Sun Z, Kalhor P, Xu Y, Liu J. Extensive numerical tests of leapfrog integrator in middle thermostat scheme in molecular simulations. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2111242] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| | - Payam Kalhor
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| | - Yang Xu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China
| |
Collapse
|
14
|
Sun Z, Liu Z. BAR‐Based Multi‐Dimensional Nonequilibrium Pulling for Indirect Construction of QM/MM Free Energy Landscapes: Varying the QM Region. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100185] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| | - Zhirong Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering Peking University Beijing 100871 China
| |
Collapse
|
15
|
|
16
|
Cao X, Tian P. "Dividing and Conquering" and "Caching" in Molecular Modeling. Int J Mol Sci 2021; 22:5053. [PMID: 34068835 PMCID: PMC8126232 DOI: 10.3390/ijms22095053] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2021] [Revised: 04/26/2021] [Accepted: 04/27/2021] [Indexed: 11/17/2022] Open
Abstract
Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, most important methodological advancements in more than half century of molecular modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but results are not transferable. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science. This framework is based on a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.
Collapse
Affiliation(s)
- Xiaoyong Cao
- School of Life Sciences, Jilin University, Changchun 130012, China;
| | - Pu Tian
- School of Life Sciences, Jilin University, Changchun 130012, China;
- School of Artificial Intelligence, Jilin University, Changchun 130012, China
| |
Collapse
|
17
|
Huai Z, Yang H, Sun Z. Binding thermodynamics and interaction patterns of human purine nucleoside phosphorylase-inhibitor complexes from extensive free energy calculations. J Comput Aided Mol Des 2021; 35:643-656. [PMID: 33759016 DOI: 10.1007/s10822-021-00382-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 03/13/2021] [Indexed: 11/29/2022]
Abstract
Human purine nucleoside phosphorylase (hPNP) plays a significant role in the catabolism of deoxyguanosine. The trimeric protein is an important target in the treatment of T-cell cancers and autoimmune disorders. Experimental studies on the inhibition of the hPNP observe that the first ligand bound to one of three subunits effectively inhibits the protein, while the binding of more ligands to the subsequent sites shows negative cooperativities. In this work, we performed extensive end-point and alchemical free energy calculations to determine the binding thermodynamics of the trimeric protein-ligand system. 13 Immucillin inhibitors with experimental results are under calculation. Two widely accepted charge schemes for small molecules including AM1-BCC and RESP are adopted for ligands. The results of RESP are in better agreement with the experimental reference. Further investigations of the interaction networks in the protein-ligand complexes reveal that several residues play significant roles in stabilizing the complex structure. The most commonly observed ones include PHE200, GLU201, MET219, and ASN243. The conformations of the protein in different protein-ligand complexes are observed to be similar. We expect these insights to aid the development of potent drugs targeting hPNP.
Collapse
Affiliation(s)
- Zhe Huai
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai, 200062, China
| | - Huaiyu Yang
- College of Engineering, Hebei Normal University, Shijiazhuang, 050024, China
| | - Zhaoxi Sun
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai, 200062, China.
| |
Collapse
|