1
|
Nakamura S, Akaki T, Nishiwaki K, Nakatani M, Kawase Y, Takahashi Y, Nakanishi I. System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study. J Comput Chem 2023; 44:824-831. [PMID: 36444861 DOI: 10.1002/jcc.27044] [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/07/2022] [Revised: 11/01/2022] [Accepted: 11/03/2022] [Indexed: 11/30/2022]
Abstract
The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.
Collapse
Affiliation(s)
- Shinya Nakamura
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Tatsuo Akaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan.,Chemical Research Laboratories, Central Pharmaceutical Research Institute, Japan Tobacco Inc., Osaka, Japan
| | - Keiji Nishiwaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Midori Nakatani
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuji Kawase
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuki Takahashi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Isao Nakanishi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| |
Collapse
|
2
|
Nakamura T, Yokaichiya T, Fedorov DG. Analysis of Guest Adsorption on Crystal Surfaces Based on the Fragment Molecular Orbital Method. J Phys Chem A 2022; 126:957-969. [PMID: 35080391 DOI: 10.1021/acs.jpca.1c10229] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
For gaining insights into interactions in periodic systems, an analysis is developed based on the fragment molecular orbital method combined with periodic boundary conditions. The adsorption energy is decomposed into guest and surface polarization and deformation energy, guest-surface and guest-guest interactions, and the vibrational free energy. The analysis is applied to the adsorption of guest molecules to Ih (001) ice surface. The cooperativity effects result in a non-linear change in the adsorption energy with coverage due to many-body effects. The role of dispersion is found to be dominant for guests with long hydrophobic tails. A rule is proposed relating the length of the alkyl tail with the formation of the guest layer. The computed binding enthalpies are in good agreement with experimental values. For high coverage, adsorbed molecules can form an ordered layer known as self-assembled monolayer.
Collapse
Affiliation(s)
- Taiji Nakamura
- 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
| | - Tomoko Yokaichiya
- 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
| | - 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
|
3
|
Watanabe C, Okiyama Y, Tanaka S, Fukuzawa K, Honma T. Molecular recognition of SARS-CoV-2 spike glycoprotein: quantum chemical hot spot and epitope analyses. Chem Sci 2021; 12:4722-4739. [PMID: 35355624 PMCID: PMC8892577 DOI: 10.1039/d0sc06528e] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2020] [Accepted: 02/21/2021] [Indexed: 12/18/2022] Open
Abstract
Due to the COVID-19 pandemic, researchers have attempted to identify complex structures of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S-protein) with angiotensin-converting enzyme 2 (ACE2) or a blocking antibody. However, the molecular recognition mechanism-critical information for drug and antibody design-has not been fully clarified at the amino acid residue level. Elucidating such a microscopic mechanism in detail requires a more accurate molecular interpretation that includes quantum mechanics to quantitatively evaluate hydrogen bonds, XH/π interactions (X = N, O, and C), and salt bridges. In this study, we applied the fragment molecular orbital (FMO) method to characterize the SARS-CoV-2 S-protein binding interactions with not only ACE2 but also the B38 Fab antibody involved in ACE2-inhibitory binding. By analyzing FMO-based interaction energies along a wide range of binding interfaces carefully, we identified amino acid residues critical for molecular recognition between S-protein and ACE2 or B38 Fab antibody. Importantly, hydrophobic residues that are involved in weak interactions such as CH-O hydrogen bond and XH/π interactions, as well as polar residues that construct conspicuous hydrogen bonds, play important roles in molecular recognition and binding ability. Moreover, through these FMO-based analyses, we also clarified novel hot spots and epitopes that had been overlooked in previous studies by structural and molecular mechanical approaches. Altogether, these hot spots/epitopes identified between S-protein and ACE2/B38 Fab antibody may provide useful information for future antibody design, evaluation of the binding property of the SARS-CoV-2 variants including its N501Y, and small or medium drug design against the SARS-CoV-2.
Collapse
Affiliation(s)
- Chiduru Watanabe
- Center for Biosystems Dynamics Research, RIKEN 1-7-22 Suehiro-cho, Tsurumi-ku Yokohama Kanagawa 230-0045 Japan +81-45-503-9432 +81-45-503-9551
- JST, PRESTO 4-1-8, Honcho Kawaguchi Saitama 332-0012 Japan
| | - Yoshio Okiyama
- Division of Medicinal Safety Science, National Institute of Health Sciences 3-25-26 Tonomachi, Kawasaki-ku Kawasaki Kanagawa 210-9501 Japan
| | - Shigenori Tanaka
- Department of Computational Science, Graduate School of System Informatics, Kobe University 1-1 Rokkodai, Nada-ku Kobe Hyogo 657-8501 Japan
| | - Kaori Fukuzawa
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University 2-4-41 Ebara, Shinagawa-ku Tokyo 142-8501 Japan
- Department of Biomolecular Engineering, Graduate School of Engineering, Tohoku University 6-6-11 Aoba, Aramaki, Aoba-ku Sendai Miyagi 980-8579 Japan
| | - Teruki Honma
- Center for Biosystems Dynamics Research, RIKEN 1-7-22 Suehiro-cho, Tsurumi-ku Yokohama Kanagawa 230-0045 Japan +81-45-503-9432 +81-45-503-9551
| |
Collapse
|
4
|
Abstract
The understanding of binding interactions between a protein and a small molecule plays a key role in the rationalization of potency and selectivity and in design of new ideas. However, even when a target of interest is structurally enabled, visual inspection and force field-based molecular mechanics calculations cannot always explain the full complexity of the molecular interactions that are critical in drug design. Quantum mechanical methods have the potential to address this shortcoming, but traditionally, computational expense has made the application of these calculations impractical. The fragment molecular orbital (FMO) method offers a solution that combines accuracy, speed, and the ability to characterize important interactions (i.e. its strength in kcal/mol and chemical nature: hydrophobic, electrostatic, etc) that would otherwise be hard to detect. In this chapter, we describe the FMO method and illustrate its application in the discovery of the benzothiazole (BZT) series as novel tyrosine kinase ITK inhibitors for treatment of allergic asthma.
Collapse
|
5
|
Analyzing GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method. Methods Mol Biol 2021. [PMID: 32016893 DOI: 10.1007/978-1-0716-0282-9_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
G-protein-coupled receptors (GPCRs) have enormous physiological and biomedical importance, and therefore it is not surprising that they are the targets of many prescribed drugs. Further progress in GPCR drug discovery is highly dependent on the availability of protein structural information. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanics (QM) approaches are often too computationally expensive to be of practical use in time-sensitive situations, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed, and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule toward ligand binding, including an analysis of their chemical nature. Such information is essential for an efficient structure-based drug design (SBDD) process. In this chapter, we describe how to use FMO in the characterization of GPCR-ligand interactions.
Collapse
|
6
|
Reiter-Scherer V, Cuellar-Camacho JL, Bhatia S, Haag R, Herrmann A, Lauster D, Rabe JP. Force Spectroscopy Shows Dynamic Binding of Influenza Hemagglutinin and Neuraminidase to Sialic Acid. Biophys J 2019; 116:1037-1048. [PMID: 30799074 DOI: 10.1016/j.bpj.2019.01.041] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2018] [Revised: 11/30/2018] [Accepted: 01/14/2019] [Indexed: 10/27/2022] Open
Abstract
The influenza A virus infects target cells through multivalent interactions of its major spike proteins, hemagglutinin (HA) and neuraminidase (NA), with the cellular receptor sialic acid (SA). HA is known to mediate the attachment of the virion to the cell, whereas NA enables the release of newly formed virions by cleaving SA from the cell. Because both proteins target the same receptor but have antagonistic functions, virus infection depends on a properly tuned balance of the kinetics of HA and NA activities for viral entry to and release from the host cell. Here, dynamic single-molecule force spectroscopy, based on scanning force microscopy, was employed to determine these bond-specific kinetics, characterized by the off rate koff, rupture length xβ and on rate kon, as well as the related free-energy barrier ΔG and the dissociation constant KD. Measurements were conducted using surface-immobilized HA and NA of the influenza A virus strain A/California/04/2009 and a novel, to our knowledge, synthetic SA-displaying receptor for functionalization of the force probe. Single-molecule force spectroscopy at force loading rates between 100 and 50,000 pN/s revealed most probable rupture forces of the protein-SA bond in the range of 10-100 pN. Using an extension of the widely applied Bell-Evans formalism by Friddle, De Yoreo, and co-workers, it is shown that HA features a smaller xβ, a larger koff and a smaller ΔG than NA. Measurements of the binding probability at increasing contact time between the scanning force microscopy force probe and the surface allow an estimation of KD, which is found to be three times as large for HA than for NA. This suggests a stronger interaction for NA-SA than for HA-SA. The biological implications in regard to virus binding to the host cell and the release of new virions from the host cell are discussed.
Collapse
Affiliation(s)
| | | | - Sumati Bhatia
- Department of Chemistry, Freie Universität Berlin, Berlin, Germany
| | - Rainer Haag
- Department of Chemistry, Freie Universität Berlin, Berlin, Germany
| | - Andreas Herrmann
- Department of Biology and IRI Life Sciences, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Daniel Lauster
- Department of Biology and IRI Life Sciences, Humboldt-Universität zu Berlin, Berlin, Germany.
| | - Jürgen P Rabe
- Department of Physics and IRIS Adlershof, Humboldt-Universität zu Berlin, Berlin, Germany.
| |
Collapse
|
7
|
Yamanaka M. <b>Random matrix theory for an inter-fragment interaction energy matrix in fragment molecular orbital method </b>. CHEM-BIO INFORMATICS JOURNAL 2018. [DOI: 10.1273/cbij.18.123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Affiliation(s)
- Masanori Yamanaka
- Department of Physics, College of Science and Technology, Nihon University
| |
Collapse
|
8
|
Chudyk EI, Sarrat L, Aldeghi M, Fedorov DG, Bodkin MJ, James T, Southey M, Robinson R, Morao I, Heifetz A. Exploring GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method. Methods Mol Biol 2018; 1705:179-195. [PMID: 29188563 DOI: 10.1007/978-1-4939-7465-8_8] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity. It is essential for an efficient structure-based drug design (SBDD) process. FMO enables ab initio approaches to be applied to systems that conventional quantum-mechanical (QM) methods would find challenging. The key advantage of the Fragment Molecular Orbital Method (FMO) is that it can reveal atomistic details about the individual contributions and chemical nature of each residue and water molecule toward ligand binding which would otherwise be difficult to detect without using QM methods. In this chapter, we demonstrate the typical use of FMO to analyze 19 crystal structures of β1 and β2 adrenergic receptors with their corresponding agonists and antagonists.
Collapse
Affiliation(s)
- Ewa I Chudyk
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Laurie Sarrat
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Matteo Aldeghi
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Dmitri G Fedorov
- CD-FMat, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki, 305-8568, Japan
| | - Mike J Bodkin
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Tim James
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Michelle Southey
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Roger Robinson
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Inaki Morao
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK.
| |
Collapse
|
9
|
Using the fragment molecular orbital method to investigate agonist-orexin-2 receptor interactions. Biochem Soc Trans 2016; 44:574-81. [PMID: 27068972 PMCID: PMC5264495 DOI: 10.1042/bst20150250] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2016] [Indexed: 12/11/2022]
Abstract
The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity and is essential for an efficient structure-based drug discovery (SBDD) process. Clearly, to begin SBDD, a structure is needed, and although there has been fantastic progress in solving G-protein-coupled receptor (GPCR) crystal structures, the process remains quite slow and is not currently feasible for every GPCR or GPCR-ligand complex. This situation significantly limits the ability of X-ray crystallography to impact the drug discovery process for GPCR targets in 'real-time' and hence there is still a need for other practical and cost-efficient alternatives. We present here an approach that integrates our previously described hierarchical GPCR modelling protocol (HGMP) and the fragment molecular orbital (FMO) quantum mechanics (QM) method to explore the interactions and selectivity of the human orexin-2 receptor (OX2R) and its recently discovered nonpeptidic agonists. HGMP generates a 3D model of GPCR structures and its complexes with small molecules by applying a set of computational methods. FMO allowsab initioapproaches to be applied to systems that conventional QM methods would find challenging. The key advantage of FMO is that it can reveal information on the individual contribution and chemical nature of each residue and water molecule to the ligand binding that normally would be difficult to detect without QM. We illustrate how the combination of both techniques provides a practical and efficient approach that can be used to analyse the existing structure-function relationships (SAR) and to drive forward SBDD in a real-world example for which there is no crystal structure of the complex available.
Collapse
|
10
|
Ratcliff LE, Mohr S, Huhs G, Deutsch T, Masella M, Genovese L. Challenges in large scale quantum mechanical calculations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1290] [Citation(s) in RCA: 70] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Laura E. Ratcliff
- Argonne Leadership Computing Facility Argonne National Laboratory Lemon IL USA
| | - Stephan Mohr
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Georg Huhs
- Department of Computer Applications in Science and Engineering Barcelona Supercomputing Center (BSC‐CNS) Barcelona Spain
| | - Thierry Deutsch
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| | - Michel Masella
- Laboratoire de Biologie Structurale et Radiologie, Service de Bioénergétique, Biologie Structurale et Mécanisme Institut de Biologie et de Technologie de Saclay, CEA Saclay Gif‐sur‐Yvette Cedex France
| | - Luigi Genovese
- University Grenoble Alpes INAC‐MEM Grenoble France
- CEA, INAC‐MEM Grenoble France
| |
Collapse
|
11
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
12
|
Heifetz A, Chudyk EI, Gleave L, Aldeghi M, Cherezov V, Fedorov DG, Biggin PC, Bodkin MJ. The Fragment Molecular Orbital Method Reveals New Insight into the Chemical Nature of GPCR–Ligand Interactions. J Chem Inf Model 2015; 56:159-72. [PMID: 26642258 DOI: 10.1021/acs.jcim.5b00644] [Citation(s) in RCA: 79] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Our interpretation of ligand-protein interactions is often informed by high-resolution structures, which represent the cornerstone of structure-based drug design. However, visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum Mechanics approaches are often too computationally expensive, but one method, Fragment Molecular Orbital (FMO), offers an excellent compromise and has the potential to reveal key interactions that would otherwise be hard to detect. To illustrate this, we have applied the FMO method to 18 Class A GPCR-ligand crystal structures, representing different branches of the GPCR genome. Our work reveals key interactions that are often omitted from structure-based descriptions, including hydrophobic interactions, nonclassical hydrogen bonds, and the involvement of backbone atoms. This approach provides a more comprehensive picture of receptor-ligand interactions than is currently used and should prove useful for evaluation of the chemical nature of ligand binding and to support structure-based drug design.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Ewa I. Chudyk
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Laura Gleave
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Matteo Aldeghi
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Vadim Cherezov
- Department
of Chemistry, Bridge Institute, University of Southern California, Los Angeles, California 90089, United States
- Laboratory
for Structural Biology of GPCRs, Moscow Institute of Physics and Technology, Dolgoprudny 141700, Russia
| | - Dmitri G. Fedorov
- NMRI, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Philip C. Biggin
- Department
of Biochemistry, University of Oxford, South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Mike J. Bodkin
- Evotec (U.K.) Ltd., 114 Innovation
Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| |
Collapse
|
13
|
Genheden S, Ryde U, Söderhjelm P. Binding affinities by alchemical perturbation using QM/MM with a large QM system and polarizable MM model. J Comput Chem 2015; 36:2114-24. [PMID: 26280564 DOI: 10.1002/jcc.24048] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2015] [Revised: 06/13/2015] [Accepted: 07/17/2015] [Indexed: 12/19/2022]
Abstract
The most general way to improve the accuracy of binding-affinity calculations for protein-ligand systems is to use quantum-mechanical (QM) methods together with rigorous alchemical-perturbation (AP) methods. We explore this approach by calculating the relative binding free energy of two synthetic disaccharides binding to galectin-3 at a reasonably high QM level (dispersion-corrected density functional theory with a triple-zeta basis set) and with a sufficiently large QM system to include all short-range interactions with the ligand (744-748 atoms). The rest of the protein is treated as a collection of atomic multipoles (up to quadrupoles) and polarizabilities. Several methods for evaluating the binding free energy from the 3600 QM calculations are investigated in terms of stability and accuracy. In particular, methods using QM calculations only at the endpoints of the transformation are compared with the recently proposed non-Boltzmann Bennett acceptance ratio (NBB) method that uses QM calculations at several stages of the transformation. Unfortunately, none of the rigorous approaches give sufficient statistical precision. However, a novel approximate method, involving the direct use of QM energies in the Bennett acceptance ratio method, gives similar results as NBB but with better precision, ∼3 kJ/mol. The statistical error can be further reduced by performing a greater number of QM calculations.
Collapse
Affiliation(s)
- Samuel Genheden
- School of Chemistry, University of Southampton, Highfield, Southampton, SO17 1BJ, United Kingdom
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| | - Pär Söderhjelm
- Department of Biophysical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| |
Collapse
|
14
|
Phipps MJS, Fox T, Tautermann CS, Skylaris CK. Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns. Chem Soc Rev 2015; 44:3177-211. [DOI: 10.1039/c4cs00375f] [Citation(s) in RCA: 226] [Impact Index Per Article: 25.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.
Collapse
Affiliation(s)
| | - Thomas Fox
- Lead Identification and Optimization Support
- Boehringer Ingelheim Pharma GmbH & Co. KG
- 88397 Biberach
- Germany
| | - Christofer S. Tautermann
- Lead Identification and Optimization Support
- Boehringer Ingelheim Pharma GmbH & Co. KG
- 88397 Biberach
- Germany
| | | |
Collapse
|
15
|
Nishimoto Y, Fedorov DG, Irle S. Density-Functional Tight-Binding Combined with the Fragment Molecular Orbital Method. J Chem Theory Comput 2014; 10:4801-12. [DOI: 10.1021/ct500489d] [Citation(s) in RCA: 73] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
| | - Dmitri G. Fedorov
- Nanosystem
Research Institute (NRI), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba 305-8565, Japan
| | | |
Collapse
|
16
|
Interaction energy analysis on specific binding of influenza virus hemagglutinin to avian and human sialosaccharide receptors: Importance of mutation-induced structural change. J Mol Graph Model 2014; 53:48-58. [DOI: 10.1016/j.jmgm.2014.07.004] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2014] [Revised: 07/04/2014] [Accepted: 07/07/2014] [Indexed: 11/19/2022]
|
17
|
Alexeev Y, Fedorov DG, Shvartsburg AA. Effective Ion Mobility Calculations for Macromolecules by Scattering on Electron Clouds. J Phys Chem A 2014; 118:6763-72. [DOI: 10.1021/jp505012c] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yuri Alexeev
- Argonne
Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Dmitri G. Fedorov
- Nanosystem
Research Institute, National Institute of Advanced Industrial Science and Technology, Tsukuba 305-8568, Japan
| | - Alexandre A. Shvartsburg
- Biological
Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
18
|
Richard RM, Lao KU, Herbert JM. Understanding the many-body expansion for large systems. I. Precision considerations. J Chem Phys 2014; 141:014108. [DOI: 10.1063/1.4885846] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- 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
|
19
|
Koyama Y, Ueno-Noto K, Takano K. Affinity of HIV-1 antibody 2G12 with monosaccharides: A theoretical study based on explicit and implicit water models. Comput Biol Chem 2014; 49:36-44. [DOI: 10.1016/j.compbiolchem.2014.01.013] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2013] [Revised: 01/14/2014] [Accepted: 01/14/2014] [Indexed: 11/25/2022]
|
20
|
Asada N, Fedorov DG, Kitaura K, Nakanishi I, Merz KM. An Efficient Method to Evaluate Intermolecular Interaction Energies in Large Systems Using Overlapping Multicenter ONIOM and the Fragment Molecular Orbital Method. J Phys Chem Lett 2012; 3:2604-2610. [PMID: 23050059 PMCID: PMC3462021 DOI: 10.1021/jz3010688] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
We propose an approach based on the overlapping multicenter ONIOM to evaluate intermolecular interaction energies in large systems and demonstrate its accuracy on several representative systems in the complete basis set limit at the MP2 and CCSD(T) level of theory. In the application to the intermolecular interaction energy between insulin dimer and 4'-hydroxyacetanilide at the MP2/CBS level, we use the fragment molecular orbital method for the calculation of the entire complex assigned to the lowest layer in three-layer ONIOM. The developed method is shown to be efficient and accurate in the evaluation of the protein-ligand interaction energies.
Collapse
Affiliation(s)
- Naoya Asada
- Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Dmitri G. Fedorov
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Kazuo Kitaura
- Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Isao Nakanishi
- Department of Pharmaceutical Sciences, Kinki University, 3-4-1, Kowakae, Higashi-Osaka, Osaka 577-8502, Japan
| | - Kenneth M. Merz
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, P.O. Box 118435, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
21
|
Nagata T, Fedorov DG, Sawada T, Kitaura K. Analysis of Solute–Solvent Interactions in the Fragment Molecular Orbital Method Interfaced with Effective Fragment Potentials: Theory and Application to a Solvated Griffithsin–Carbohydrate Complex. J Phys Chem A 2012; 116:9088-99. [DOI: 10.1021/jp304991a] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Takeshi Nagata
- Graduate School of Pharmaceutical
Sciences, Kyoto University, 46-29 Yoshidashimo
Adachi, Sakyo-ku, Kyoto 606-8501, Japan
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Dmitri G. Fedorov
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Toshihiko Sawada
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
- Core Research for Evolutional
Science and Technology (CREST), Japan Science and Technology Agency (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012,
Japan
| | - Kazuo Kitaura
- Graduate School of Pharmaceutical
Sciences, Kyoto University, 46-29 Yoshidashimo
Adachi, Sakyo-ku, Kyoto 606-8501, Japan
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
22
|
Nakata H, Fedorov DG, Nagata T, Yokojima S, Ogata K, Kitaura K, Nakamura S. Unrestricted Hartree-Fock based on the fragment molecular orbital method: Energy and its analytic gradient. J Chem Phys 2012; 137:044110. [DOI: 10.1063/1.4737860] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Nagata T, Fedorov DG, Li H, Kitaura K. Analytic gradient for second order Møller-Plesset perturbation theory with the polarizable continuum model based on the fragment molecular orbital method. J Chem Phys 2012; 136:204112. [DOI: 10.1063/1.4714601] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Bowler DR, Miyazaki T. O(N) methods in electronic structure calculations. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2012; 75:036503. [PMID: 22790422 DOI: 10.1088/0034-4885/75/3/036503] [Citation(s) in RCA: 210] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Linear-scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high-performance computers. The linear-scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas are then discussed. The applications of linear-scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear-scaling methods are discussed.
Collapse
Affiliation(s)
- D R Bowler
- London Centre for Nanotechnology, UCL, 17-19 Gordon St, London WC1H 0AH, UK.
| | | |
Collapse
|
25
|
Fedorov DG, Nagata T, Kitaura K. Exploring chemistry with the fragment molecular orbital method. Phys Chem Chem Phys 2012; 14:7562-77. [DOI: 10.1039/c2cp23784a] [Citation(s) in RCA: 290] [Impact Index Per Article: 24.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
26
|
Fedorov DG, Kitaura K. Energy Decomposition Analysis in Solution Based on the Fragment Molecular Orbital Method. J Phys Chem A 2011; 116:704-19. [DOI: 10.1021/jp209579w] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Dmitri G. Fedorov
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Kazuo Kitaura
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
- Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| |
Collapse
|
27
|
Prediction of probable mutations in influenza virus hemagglutinin protein based on large-scale ab initio fragment molecular orbital calculations. J Mol Graph Model 2011; 30:110-9. [DOI: 10.1016/j.jmgm.2011.06.011] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2011] [Revised: 06/24/2011] [Accepted: 06/27/2011] [Indexed: 11/22/2022]
|
28
|
Gordon MS, Fedorov DG, Pruitt SR, Slipchenko LV. Fragmentation Methods: A Route to Accurate Calculations on Large Systems. Chem Rev 2011; 112:632-72. [DOI: 10.1021/cr200093j] [Citation(s) in RCA: 836] [Impact Index Per Article: 64.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Mark S. Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Dmitri G. Fedorov
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Spencer R. Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Lyudmila V. Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
29
|
Okiyama Y, Fukuzawa K, Yamada H, Mochizuki Y, Nakano T, Tanaka S. Counterpoise-corrected interaction energy analysis based on the fragment molecular orbital scheme. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.04.070] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|