1
|
Gupta AK, Maier S, Thapa B, Raghavachari K. Toward Post-Hartree-Fock Accuracy for Protein-Ligand Affinities Using the Molecules-in-Molecules Fragmentation-Based Method. J Chem Theory Comput 2024; 20:2774-2785. [PMID: 38530869 DOI: 10.1021/acs.jctc.3c01293] [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: 03/28/2024]
Abstract
The complexity and size of large molecular systems, such as protein-ligand complexes, pose computational challenges for accurate post-Hartree-Fock calculations. This study delivers a thorough benchmarking of the Molecules-in-Molecules (MIM) method, presenting a clear and accessible strategy for layer/theory selections in post-Hartree-Fock computations on substantial molecular systems, notably protein-ligand complexes. An approach is articulated, enabling augmented computational efficiency by strategically canceling out common subsystem energy terms between complexes and proteins within the supermolecular equation. Employing DLPNO-based post-Hartree-Fock methods in conjunction with the three-layer MIM method (MIM3), this study demonstrates the achievement of protein-ligand binding energies with remarkable accuracy (errors <1 kcal mol-1), while significantly reducing computational costs. Furthermore, noteworthy correlations between theoretically computed interaction energies and their experimental equivalents were observed, with R2 values of approximately 0.90 and 0.78 for CDK2 and BZT-ITK sets, respectively, thus validating the efficacy of the MIM method in calculating binding energies. By highlighting the crucial role of diffuse or small Pople-style basis sets in the middle layer for reducing energy errors, this work provides valuable insights and practical methodologies for interaction energy computations in large molecular complexes and opens avenues for their application across a diverse range of molecular systems.
Collapse
Affiliation(s)
- Ankur K Gupta
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Sarah Maier
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
2
|
Lim H, Kang DH, Kim J, Pellow-Jarman A, McFarthing S, Pellow-Jarman R, Jeon HN, Oh B, Rhee JKK, No KT. Fragment molecular orbital-based variational quantum eigensolver for quantum chemistry in the age of quantum computing. Sci Rep 2024; 14:2422. [PMID: 38287087 PMCID: PMC10825197 DOI: 10.1038/s41598-024-52926-3] [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: 12/01/2023] [Accepted: 01/25/2024] [Indexed: 01/31/2024] Open
Abstract
Quantum computers offer significant potential for complex system analysis, yet their application in large systems is hindered by limitations such as qubit availability and quantum hardware noise. While the variational quantum eigensolver (VQE) was proposed to address these issues, its scalability remains limited. Many efforts, including new ansätze and Hamiltonian modifications, have been made to overcome these challenges. In this work, we introduced the novel Fragment Molecular Orbital/Variational Quantum Eigensolver (FMO/VQE) algorithm. This method combines the fragment molecular orbital (FMO) approach with VQE and efficiently utilizes qubits for quantum chemistry simulations. Employing the UCCSD ansatz, the FMO/VQE achieved an absolute error of just 0.053 mHa with 8 qubits in a [Formula: see text] system using the STO-3G basis set, and an error of 1.376 mHa with 16 qubits in a [Formula: see text] system with the 6-31G basis set. These results indicated a significant advancement in scalability over conventional VQE, maintaining accuracy with fewer qubits. Therefore, our FMO/VQE method exemplifies how integrating fragment-based quantum chemistry with quantum algorithms can enhance scalability, facilitating more complex molecular simulations and aligning with quantum computing advancements.
Collapse
Affiliation(s)
- Hocheol Lim
- Bioinformatics and Molecular Design Research Center (BMDRC), Incheon, Republic of Korea.
| | | | - Jeonghoon Kim
- Bioinformatics and Molecular Design Research Center (BMDRC), Incheon, Republic of Korea
| | | | | | | | | | - Byungdu Oh
- Baobab AiBIO Co., Ltd., Incheon, Republic of Korea
- SKKU Advanced Institute of Nanotechnology, Sungkyunkwan University, Suwon, Republic of Korea
| | | | - Kyoung Tai No
- Bioinformatics and Molecular Design Research Center (BMDRC), Incheon, Republic of Korea.
- Baobab AiBIO Co., Ltd., Incheon, Republic of Korea.
| |
Collapse
|
3
|
Huo Y, Zhao C, Wang Y, Wang S, Mu T, Du W. Roles of Apigenin and Nepetin in the Assembly Behavior and Cytotoxicity of Prion Neuropeptide PrP106-126. ACS Chem Neurosci 2024; 15:245-257. [PMID: 38133816 DOI: 10.1021/acschemneuro.3c00417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023] Open
Abstract
Development of potential inhibitors to prevent prion protein (PrP) fibrillation is a therapeutic strategy for prion diseases. The prion neuropeptide PrP106-126, a research model of abnormal PrP (PrPSc), presents similar physicochemical and biochemical characters to PrPSc, which is also a target of potential inhibitors against prion deposition. Many flavones have antioxidant, anti-inflammatory, and antibacterial properties, and they are applied in treating prion disorder and other amyloidosis as well. However, the inhibition mechanism of flavones on PrP106-126 fibrillation is still unclear. In the current work, apigenin and nepetin were used to suppress the aggregation of PrP106-126 and to alleviate the peptide-induced cytotoxicity. The results showed that apigenin and nepetin impeded the fibril formation of PrP106-126 and depolymerized the preformed fibrils. They were bound to PrP106-126 predominantly by hydrophobic and hydrogen bonding interactions. In addition, both flavones upregulated cell viability and decreased membrane leakage through reducing peptide oligomerization. The differences in inhibition and cell protection between the two small molecules were presumably attributed to the substitution of hydroxyl and methoxy groups in nepetin, which demonstrated the significant structure-function relationship of flavones with prion neuropeptide and the prospect of flavonoids as drug candidates against prion diseases.
Collapse
Affiliation(s)
- Yan Huo
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| | - Cong Zhao
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| | - Yanan Wang
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| | - Shao Wang
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| | - Tiancheng Mu
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| | - Weihong Du
- Department of Chemistry, Renmin University of China, Beijing 100872, China
| |
Collapse
|
4
|
Yuan Z, Chen X, Fan S, Chang L, Chu L, Zhang Y, Wang J, Li S, Xie J, Hu J, Miao R, Zhu L, Zhao Z, Li H, Li S. Binding Free Energy Calculation Based on the Fragment Molecular Orbital Method and Its Application in Designing Novel SHP-2 Allosteric Inhibitors. Int J Mol Sci 2024; 25:671. [PMID: 38203841 PMCID: PMC10779950 DOI: 10.3390/ijms25010671] [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: 11/30/2023] [Revised: 12/28/2023] [Accepted: 12/29/2023] [Indexed: 01/12/2024] Open
Abstract
The accurate prediction of binding free energy is a major challenge in structure-based drug design. Quantum mechanics (QM)-based approaches show promising potential in predicting ligand-protein binding affinity by accurately describing the behavior and structure of electrons. However, traditional QM calculations face computational limitations, hindering their practical application in drug design. Nevertheless, the fragment molecular orbital (FMO) method has gained widespread application in drug design due to its ability to reduce computational costs and achieve efficient ab initio QM calculations. Although the FMO method has demonstrated its reliability in calculating the gas phase potential energy, the binding of proteins and ligands also involves other contributing energy terms, such as solvent effects, the 'deformation energy' of a ligand's bioactive conformations, and entropy. Particularly in cases involving ionized fragments, the calculation of solvation free energy becomes particularly crucial. We conducted an evaluation of some previously reported implicit solvent methods on the same data set to assess their potential for improving the performance of the FMO method. Herein, we develop a new QM-based binding free energy calculation method called FMOScore, which enhances the performance of the FMO method. The FMOScore method incorporates linear fitting of various terms, including gas-phase potential energy, deformation energy, and solvation free energy. Compared to other widely used traditional prediction methods such as FEP+, MM/PBSA, MM/GBSA, and Autodock vina, FMOScore showed good performance in prediction accuracies. By constructing a retrospective case study, it was observed that incorporating calculations for solvation free energy and deformation energy can further enhance the precision of FMO predictions for binding affinity. Furthermore, using FMOScore-guided lead optimization against Src homology-2-containing protein tyrosine phosphatase 2 (SHP-2), we discovered a novel and potent allosteric SHP-2 inhibitor (compound 8).
Collapse
Affiliation(s)
- Zhen Yuan
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Xingyu Chen
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Sisi Fan
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Longfeng Chang
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Linna Chu
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Ying Zhang
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Jie Wang
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Shuang Li
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Jinxin Xie
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Jianguo Hu
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Runyu Miao
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Lili Zhu
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Zhenjiang Zhao
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
| | - Honglin Li
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
- Innovation Center for AI and Drug Discovery, East China Normal University, Shanghai 200062, China
- Lingang Laboratory, Shanghai 200031, China
| | - Shiliang Li
- Shanghai Key Laboratory of New Drug Design, State Key Laboratory of Bioreactor Engineering, School of Pharmacy, East China University of Science & Technology, Shanghai 200237, China; (Z.Y.); (X.C.); (S.F.); (Z.Z.)
- Innovation Center for AI and Drug Discovery, East China Normal University, Shanghai 200062, China
| |
Collapse
|
5
|
Huo Y, Huang X, Wang Y, Zhao C, Zheng T, Du W. Inhibitory effects of sesquiterpene lactones on the aggregation and cytotoxicity of prion neuropeptide. Biochimie 2023; 211:131-140. [PMID: 36963557 DOI: 10.1016/j.biochi.2023.03.012] [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: 01/16/2023] [Revised: 03/12/2023] [Accepted: 03/21/2023] [Indexed: 03/26/2023]
Abstract
The misfolding and conformational transformation of prion protein (PrP) are crucial to the progression of prion diseases. Screening for available natural inhibitors against prion proteins can contribute to the rational design and development of new anti-prion drugs and therapeutic strategies. The prion neuropeptide, PrP106-126 is commonly used as a model peptide of the abnormal PrPSc, and a number of potential inhibitors were explored against the amyloid fibril formation of PrP106-126. The well-known sesquiterpene lactone, artemisinin, shows diverse biological functions in anti-malarial, anti-cancer and lowering glucose. However, its inhibitory effect on PrP106-126 fibrillation is unclear. In this work, we selected two sesquiterpene lactones, artemisinin (1) and artesunate (2), to explore their roles in PrP106-126 aggregation by a series of physicochemical and biochemical methods. The results demonstrated that 1 and 2 could effectively impede the formation of amyloid fibrils and remodel the preformed fibrils. The binding of the small molecules to PrP106-126 was dominated by electrostatic, hydrophobic and hydrogen bonding interactions. In addition, both compounds exhibited neuroprotective effects by reducing peptide oligomerization. 2 showed better inhibition and regulation on peptide aggregation and cellular viability than 1 due to its specific succinate modification. Our study provides the information of sesquiterpene lactones to prevent PrP fibril formation and other related amyloidosis.
Collapse
Affiliation(s)
- Yan Huo
- Department of Chemistry, Renmin University of China, Beijing, 100872, China
| | - Xiangyi Huang
- Department of Chemistry, Renmin University of China, Beijing, 100872, China
| | - Yanan Wang
- Department of Chemistry, Renmin University of China, Beijing, 100872, China
| | - Cong Zhao
- Department of Chemistry, Renmin University of China, Beijing, 100872, China
| | - Ting Zheng
- Department of Chemistry, Renmin University of China, Beijing, 100872, China
| | - Weihong Du
- Department of Chemistry, Renmin University of China, Beijing, 100872, China.
| |
Collapse
|
6
|
Luo L, Yang J, Wang C, Wu J, Li Y, Zhang X, Li H, Zhang H, Zhou Y, Lu A, Chen S. Natural products for infectious microbes and diseases: an overview of sources, compounds, and chemical diversities. SCIENCE CHINA. LIFE SCIENCES 2022; 65:1123-1145. [PMID: 34705221 PMCID: PMC8548270 DOI: 10.1007/s11427-020-1959-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Accepted: 07/27/2021] [Indexed: 12/13/2022]
Abstract
As coronavirus disease 2019 (COVID-19) threatens human health globally, infectious disorders have become one of the most challenging problem for the medical community. Natural products (NP) have been a prolific source of antimicrobial agents with widely divergent structures and a range vast biological activities. A dataset comprising 618 articles, including 646 NP-based compounds from 672 species of natural sources with biological activities against 21 infectious pathogens from five categories, was assembled through manual selection of published articles. These data were used to identify 268 NP-based compounds classified into ten groups, which were used for network pharmacology analysis to capture the most promising lead-compounds such as agelasine D, dicumarol, dihydroartemisinin and pyridomycin. The distribution of maximum Tanimoto scores indicated that compounds which inhibited parasites exhibited low diversity, whereas the chemistries inhibiting bacteria, fungi, and viruses showed more structural diversity. A total of 331 species of medicinal plants with compounds exhibiting antimicrobial activities were selected to classify the family sources. The family Asteraceae possesses various compounds against C. neoformans, the family Anacardiaceae has compounds against Salmonella typhi, the family Cucurbitacea against the human immunodeficiency virus (HIV), and the family Ancistrocladaceae against Plasmodium. This review summarizes currently available data on NP-based antimicrobials against refractory infections to provide information for further discovery of drugs and synthetic strategies for anti-infectious agents.
Collapse
Affiliation(s)
- Lu Luo
- Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, 100700, China
| | - Jun Yang
- Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, 100700, China
| | - Cheng Wang
- Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, 100006, China
| | - Jie Wu
- Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, 100700, China
| | - Yafang Li
- Institute of Hematology & Blood Diseases Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College, Tianjin, 300020, China
| | - Xu Zhang
- weMED Health, Houston, 77054, USA
| | - Hui Li
- Institute of Botany, Chinese Academy of Sciences, Beijing, 100093, China
| | - Hui Zhang
- Akupunktur Akademiet, Aabyhoej, Aarhus, 8230, Denmark
| | - Yumei Zhou
- The Fourth Clinical Medical College of Guangzhou University of Chinese Medicine, Guangzhou, 518033, China
| | - Aiping Lu
- School of Chinese Medicine, Hong Kong Baptist University, Hong Kong, 999077, China
| | - Shilin Chen
- Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, 100700, China.
| |
Collapse
|
7
|
Shimazaki T, Tachikawa M. Collaborative Approach between Explainable Artificial Intelligence and Simplified Chemical Interactions to Explore Active Ligands for Cyclin-Dependent Kinase 2. ACS OMEGA 2022; 7:10372-10381. [PMID: 35382271 PMCID: PMC8973106 DOI: 10.1021/acsomega.1c06976] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 03/09/2022] [Indexed: 05/13/2023]
Abstract
To improve virtual screening for drug discovery, we present a collaborative approach between explainable artificial intelligence (AI) and simplified chemical interaction scores to efficiently search for active ligands bound to the target receptor. In particular, we focus on cyclin-dependent kinase 2 (CDK2), which is well known as a cancer target protein. Docking simulation alone is insufficient to distinguish active ligands from decoy molecules. To identify active ligands, in this paper, machine learning is employed together with scoring functions that simplify the screened Coulomb and Lennard-Jones interactions between the ligands and residues of the target receptor. We demonstrate that these simplified interaction scores can significantly improve the classification ability of machine learning models. We also demonstrate that explainable AI together with the simplified scoring method can highlight the important residues of CDK2 for recognizing active ligands.
Collapse
Affiliation(s)
- Tomomi Shimazaki
- Graduate
School of Nanobioscience, Yokohama City
University, 22-2 Seto, Yokohama, Kanagawa 236-0027, Japan
| | - Masanori Tachikawa
- Graduate
School of Data Science, Yokohama City University, 22-2, Seto, Yokohama, Kanagawa 236-0027, Japan
| |
Collapse
|
8
|
Current and Future Challenges in Modern Drug Discovery. METHODS IN MOLECULAR BIOLOGY (CLIFTON, N.J.) 2021; 2114:1-17. [PMID: 32016883 DOI: 10.1007/978-1-0716-0282-9_1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Drug discovery is an expensive, time-consuming, and risky business. To avoid late-stage failure, learnings from past projects and the development of new approaches are crucial. New modalities and emerging new target spaces allow the exploration of unprecedented indications or to address so far undrugable targets. Late-stage attrition is usually attributed to the lack of efficacy or to compound-related safety issues. Efficacy has been shown to be related to a strong genetic link to human disease, a better understanding of the target biology, and the availability of biomarkers to bridge from animals to humans. Compound safety can be improved by ligand optimization, which is becoming increasingly demanding for difficult targets. Therefore, new strategies include the design of allosteric ligands, covalent binders, and other modalities. Design methods currently heavily rely on artificial intelligence and advanced computational methods such as free energy calculations and quantum chemistry. Especially for quantum chemical methods, a more detailed overview is given in this chapter.
Collapse
|
9
|
Muronetz VI, Barinova K, Kudryavtseva S, Medvedeva M, Melnikova A, Sevostyanova I, Semenyuk P, Stroylova Y, Sova M. Natural and Synthetic Derivatives of Hydroxycinnamic Acid Modulating the Pathological Transformation of Amyloidogenic Proteins. Molecules 2020; 25:E4647. [PMID: 33053854 PMCID: PMC7594092 DOI: 10.3390/molecules25204647] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2020] [Revised: 10/08/2020] [Accepted: 10/09/2020] [Indexed: 02/06/2023] Open
Abstract
This review presents the main properties of hydroxycinnamic acid (HCA) derivatives and their potential application as agents for the prevention and treatment of neurodegenerative diseases. It is partially focused on the successful use of these compounds as inhibitors of amyloidogenic transformation of proteins. Firstly, the prerequisites for the emergence of interest in HCA derivatives, including natural compounds, are described. A separate section is devoted to synthesis and properties of HCA derivatives. Then, the results of molecular modeling of HCA derivatives with prion protein as well as with α-synuclein fibrils are summarized, followed by detailed analysis of the experiments on the effect of natural and synthetic HCA derivatives, as well as structurally similar phenylacetic and benzoic acid derivatives, on the pathological transformation of prion protein and α-synuclein. The ability of HCA derivatives to prevent amyloid transformation of some amyloidogenic proteins, and their presence not only in food products but also as natural metabolites in human blood and tissues, makes them promising for the prevention and treatment of neurodegenerative diseases of amyloid nature.
Collapse
Affiliation(s)
- Vladimir I. Muronetz
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, 119992 Moscow, Russia; (K.B.); (A.M.); (I.S.); (P.S.); (Y.S.)
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, 119234 Moscow, Russia; (S.K.); (M.M.)
| | - Kseniya Barinova
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, 119992 Moscow, Russia; (K.B.); (A.M.); (I.S.); (P.S.); (Y.S.)
| | - Sofia Kudryavtseva
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, 119234 Moscow, Russia; (S.K.); (M.M.)
| | - Maria Medvedeva
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, 119234 Moscow, Russia; (S.K.); (M.M.)
| | - Aleksandra Melnikova
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, 119992 Moscow, Russia; (K.B.); (A.M.); (I.S.); (P.S.); (Y.S.)
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, 119234 Moscow, Russia; (S.K.); (M.M.)
| | - Irina Sevostyanova
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, 119992 Moscow, Russia; (K.B.); (A.M.); (I.S.); (P.S.); (Y.S.)
| | - Pavel Semenyuk
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, 119992 Moscow, Russia; (K.B.); (A.M.); (I.S.); (P.S.); (Y.S.)
| | - Yulia Stroylova
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, 119992 Moscow, Russia; (K.B.); (A.M.); (I.S.); (P.S.); (Y.S.)
- Institute of Molecular Medicine, Sechenov First Moscow State Medical University Trubetskaya St. 8, Bldg. 2, 119991 Moscow, Russia
| | - Matej Sova
- Faculty of Pharmacy, University of Ljubljana, Aškerčeva 7, 1000 Ljubljana, Slovenia;
| |
Collapse
|
10
|
Nutho B, Mahalapbutr P, Hengphasatporn K, Pattaranggoon NC, Simanon N, Shigeta Y, Hannongbua S, Rungrotmongkol T. Why Are Lopinavir and Ritonavir Effective against the Newly Emerged Coronavirus 2019? Atomistic Insights into the Inhibitory Mechanisms. Biochemistry 2020; 59:1769-1779. [PMID: 32293875 PMCID: PMC7184878 DOI: 10.1021/acs.biochem.0c00160] [Citation(s) in RCA: 153] [Impact Index Per Article: 38.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Revised: 04/08/2020] [Indexed: 12/12/2022]
Abstract
Since the emergence of a novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first reported from Wuhan, China, neither a specific vaccine nor an antiviral drug against SARS-CoV-2 has become available. However, a combination of two HIV-1 protease inhibitors, lopinavir and ritonavir, has been found to be effective against SARS-CoV, and both drugs could bind well to the SARS-CoV 3C-like protease (SARS-CoV 3CLpro). In this work, molecular complexation between each inhibitor and SARS-CoV-2 3CLpro was studied using all-atom molecular dynamics simulations, free energy calculations, and pair interaction energy analyses based on MM/PB(GB)SA and FMO-MP2/PCM/6-31G* methods. Both anti-HIV drugs interacted well with the residues at the active site of SARS-CoV-2 3CLpro. Ritonavir showed a somewhat higher number atomic contacts, a somewhat higher binding efficiency, and a somewhat higher number of key binding residues compared to lopinavir, which correspond with the slightly lower water accessibility at the 3CLpro active site. In addition, only ritonavir could interact with the oxyanion hole residues N142 and G143 via the formation of two hydrogen bonds. The interactions in terms of electrostatics, dispersion, and charge transfer played an important role in the drug binding. The obtained results demonstrated how repurposed anti-HIV drugs could be used to combat COVID-19.
Collapse
Affiliation(s)
- Bodee Nutho
- Center of Excellence in Computational Chemistry
(CECC), Department of Chemistry, Faculty of Science, Chulalongkorn
University, Bangkok 10330, Thailand
| | - Panupong Mahalapbutr
- Structural and Computational Biology Research Unit,
Department of Biochemistry, Faculty of Science, Chulalongkorn
University, Bangkok 10330, Thailand
| | - Kowit Hengphasatporn
- Center for Computational Sciences,
University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki
305-8577, Japan
| | | | - Nattapon Simanon
- Program in Bioinformatics and Computational Biology,
Graduate School, Chulalongkorn University, Bangkok 10330,
Thailand
| | - Yasuteru Shigeta
- Center for Computational Sciences,
University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki
305-8577, Japan
| | - Supot Hannongbua
- Center of Excellence in Computational Chemistry
(CECC), Department of Chemistry, Faculty of Science, Chulalongkorn
University, Bangkok 10330, Thailand
| | - Thanyada Rungrotmongkol
- Structural and Computational Biology Research Unit,
Department of Biochemistry, Faculty of Science, Chulalongkorn
University, Bangkok 10330, Thailand
- Program in Bioinformatics and Computational Biology,
Graduate School, Chulalongkorn University, Bangkok 10330,
Thailand
| |
Collapse
|
11
|
Heifetz A, Morao I, Babu MM, James T, Southey MWY, Fedorov DG, Aldeghi M, Bodkin MJ, Townsend-Nicholson A. Characterizing Interhelical Interactions of G-Protein Coupled Receptors with the Fragment Molecular Orbital Method. J Chem Theory Comput 2020; 16:2814-2824. [PMID: 32096994 PMCID: PMC7161079 DOI: 10.1021/acs.jctc.9b01136] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
G-protein coupled receptors (GPCRs) are the largest superfamily of membrane proteins, regulating almost every aspect of cellular activity and serving as key targets for drug discovery. We have identified an accurate and reliable computational method to characterize the strength and chemical nature of the interhelical interactions between the residues of transmembrane (TM) domains during different receptor activation states, something that cannot be characterized solely by visual inspection of structural information. Using the fragment molecular orbital (FMO) quantum mechanics method to analyze 35 crystal structures representing different branches of the class A GPCR family, we have identified 69 topologically equivalent TM residues that form a consensus network of 51 inter-TM interactions, providing novel results that are consistent with and help to rationalize experimental data. This discovery establishes a comprehensive picture of how defined molecular forces govern specific interhelical interactions which, in turn, support the structural stability, ligand binding, and activation of GPCRs.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec
(U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
- Institute
of Structural & Molecular Biology, Research Department of Structural
& Molecular Biology, Division of Biosciences, University College London, London, WC1E 6BT, United Kingdom
- E-mail: (A.H.)
| | - Inaki Morao
- Evotec
(U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
- E-mail: (I.M.)
| | - M. Madan Babu
- MRC
Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge CB2 0QH, United Kingdom
| | - Tim James
- Evotec
(U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | | | - Dmitri G. Fedorov
- CD-FMat,
National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Matteo Aldeghi
- Department
of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, 37077 Göttingen, Germany
| | - Michael J. Bodkin
- Evotec
(U.K.) Ltd., 114 Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - Andrea Townsend-Nicholson
- Institute
of Structural & Molecular Biology, Research Department of Structural
& Molecular Biology, Division of Biosciences, University College London, London, WC1E 6BT, United Kingdom
| |
Collapse
|
12
|
Kaliakin DS, Nakata H, Kim Y, Chen Q, Fedorov DG, Slipchenko LV. FMOxFMO: Elucidating Excitonic Interactions in the Fenna-Matthews-Olson Complex with the Fragment Molecular Orbital Method. J Chem Theory Comput 2020; 16:1175-1187. [PMID: 31841349 DOI: 10.1021/acs.jctc.9b00621] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In order to study Förster resonance energy transfer (FRET), the fragment molecular orbital (FMO) method is extended to compute electronic couplings between local excitations via the excited state transition density model, enabling efficient calculations of nonlocal excitations in a large molecular system and overcoming the previous limitation of being able to compute only local excitations. The results of these simple but accurate models are validated against full quantum calculations without fragmentation. The developed method is applied to a very important photosynthetic pigment-protein complex, the Fenna-Matthews-Olson complex (FMOc), that is responsible for the energy transfer from a chlorosome to the reaction center in the green sulfur bacteria. Absorption and circular dichroism spectra of FMOc are simulated, and the role of the molecular environment on the excitations is revealed.
Collapse
Affiliation(s)
- Danil S Kaliakin
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Hiroya Nakata
- Research Institute for Advanced Materials and Devices , Kyocera , 5-3 Hikaridai-3 , Seika-cho Soraku-gun, Kyoto 619-0237 , Japan
| | - Yongbin Kim
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Qifeng Chen
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1 , Tsukuba 305-8568 , Japan
| | - Lyudmila V Slipchenko
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| |
Collapse
|
13
|
Fedorov DG. Solvent Screening in Zwitterions Analyzed with the Fragment Molecular Orbital Method. J Chem Theory Comput 2019; 15:5404-5416. [PMID: 31461277 DOI: 10.1021/acs.jctc.9b00715] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Based on induced solvent charges, a new model of solvent screening is developed in the framework of the fragment molecular orbital combined with the polarizable continuum model. The developed model is applied to analyze interactions in a prototypical zwitterionic system, sodium chloride in water, and it is shown that the large underestimation of the interaction in the original solvent screening based on local charges is successfully corrected. The model is also applied to a complex of the Trp-cage (PDB: 1L2Y ) miniprotein with an anionic ligand, and the physical factors determined protein-ligand binding in solution are unraveled.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1 , Tsukuba 305-8568 , Japan
| |
Collapse
|
14
|
Liu KY, Carter-Fenk K, Herbert JM. Self-consistent charge embedding at very low cost, with application to symmetry-adapted perturbation theory. J Chem Phys 2019; 151:031102. [DOI: 10.1063/1.5111869] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Kuan-Yu Liu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
15
|
Fedorov DG, Brekhov A, Mironov V, Alexeev Y. Molecular Electrostatic Potential and Electron Density of Large Systems in Solution Computed with the Fragment Molecular Orbital Method. J Phys Chem A 2019; 123:6281-6290. [DOI: 10.1021/acs.jpca.9b04936] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Anton Brekhov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Vladimir Mironov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Yuri Alexeev
- Argonne Leadership Computing Facility and Computational Science Division, Argonne National Laboratory, Argonne, Illinois, 60439, United States
| |
Collapse
|
16
|
A Strength-Weaknesses-Opportunities-Threats (SWOT) Analysis of Cheminformatics in Natural Product Research. PROGRESS IN THE CHEMISTRY OF ORGANIC NATURAL PRODUCTS 2019; 110:239-271. [PMID: 31621015 DOI: 10.1007/978-3-030-14632-0_7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Cheminformatics-based techniques, such as molecular modeling, docking, virtual screening, and machine learning, are well accepted for their usefulness in drug discovery and development of therapeutically relevant small molecules. Although delayed by several decades, their application in natural product research has led to outstanding findings. Combining information obtained from different sources, i.e., virtual predictions, traditional medicine, structural, biochemical, and biological data, and handling big data effectively will open up new possibilities, but also challenges in the future. Strategies and examples will be presented on how to integrate cheminformatics in pharmacognostic workflows to benefit from these two highly complementary disciplines toward streamlining experimental efforts. While considering their limits and pitfalls and by exploiting their potential, computer-aided strategies should successfully guide future studies and thereby augment our knowledge of bioactive natural lead structures.
Collapse
|
17
|
Nakata H, Fedorov DG. Simulations of infrared and Raman spectra in solution using the fragment molecular orbital method. Phys Chem Chem Phys 2019; 21:13641-13652. [DOI: 10.1039/c9cp00940j] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
Calculation of IR and Raman spectra in solution for large molecular systems made possible with analytic FMO/PCM Hessians.
Collapse
Affiliation(s)
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat)
- National Institute of Advanced Industrial Science and Technology (AIST)
- Tsukuba
- Japan
| |
Collapse
|