1
|
Fedorov DG. Partition analysis of dipole moments in solution applied to functional groups in polypeptide motifs. Phys Chem Chem Phys 2024; 26:18614-18628. [PMID: 38919134 DOI: 10.1039/d4cp01654h] [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: 06/27/2024]
Abstract
A partition analysis based on segments is developed for density functional theory defining solute dipole moments of functional groups, and the corresponding induced solvent dipoles representing solvent screening. The accuracy of dipoles from the fragment molecular orbital method is evaluated in comparison to unfragmented values. The analysis is applied to evaluate dipole moments of side chains, amino and carbonyl groups in common polypeptide motifs, α-helixes, β-turns, and random coils in solution. The membrane domain of the ATP synthase (1B9U) is analyzed, revealing the effect of the bend splitting of the α-helix into two pieces.
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
|
2
|
Guareschi R, Lukac I, Gilbert IH, Zuccotto F. SophosQM: Accurate Binding Affinity Prediction in Compound Optimization. ACS OMEGA 2023; 8:15083-15098. [PMID: 37151542 PMCID: PMC10157843 DOI: 10.1021/acsomega.2c08132] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 04/07/2023] [Indexed: 05/09/2023]
Abstract
The optimization of compounds' binding affinity for a biological target is a crucial aspect of the drug development process. Being able to accurately predict binding energies in advance of synthesizing compounds would have a massive impact on the speed of the drug discovery process. The ideal binding affinity prediction method should combine accuracy, reliability, and speed. In this paper, we present SophosQM, a quantum mechanics (QM)-based approach, which can accurately predict the binding affinities of compounds to proteins. The binding affinity predictive models generated by SophosQM are based on the fragment molecular orbital (FMO) method to estimate the enthalpic component of the binding free energy, and a macroscopic descriptor, clog P, is used as an approximation of the entropic component. The affinity prediction is performed using multilinear regression, fitting the experimental values against the FMO-computed enthalpic term and clog P. The quality of the prediction can be assessed in terms of the correlation coefficient between experimental and predicted values. In this work, the method's reliability and accuracy are exemplified by applying SophosQM to 70 compounds binding to six different targets of pharmaceutical relevance. Overall, the results show a very satisfactory performance with a global correlation coefficient in the order of 0.9. Our predictions also show a satisfactory performance compared to data based on free energy perturbation. Finally, SophosQM can also be applied in high-throughput mode by using semiempirical QM methods to evaluate large portions of chemical space, while retaining a good level of accuracy, but decreasing the computing time to just a few seconds per compound.
Collapse
|
3
|
Remington JM, McKay KT, Beckage NB, Ferrell JB, Schneebeli ST, Li J. GPCRLigNet: rapid screening for GPCR active ligands using machine learning. J Comput Aided Mol Des 2023; 37:147-156. [PMID: 36840893 PMCID: PMC10379640 DOI: 10.1007/s10822-023-00497-2] [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: 11/16/2022] [Accepted: 02/03/2023] [Indexed: 02/26/2023]
Abstract
Molecules with bioactivity towards G protein-coupled receptors represent a subset of the vast space of small drug-like molecules. Here, we compare machine learning models, including dilated graph convolutional networks, that conduct binary classification to quickly identify molecules with activity towards G protein-coupled receptors. The models are trained and validated using a large set of over 600,000 active, inactive, and decoy compounds. The best performing machine learning model, dubbed GPCRLigNet, was a surprisingly simple feedforward dense neural network mapping from Morgan fingerprints to activity. Incorporation of GPCRLigNet into a high-throughput virtual screening workflow is demonstrated with molecular docking towards a particular G protein-coupled receptor, the pituitary adenylate cyclase-activating polypeptide receptor type 1. Through rigorous comparison of docking scores for molecules selected with and without using GPCRLigNet, we demonstrate an enrichment of potentially potent molecules using GPCRLigNet. This work provides a proof of principle that GPCRLigNet can effectively hone the chemical search space towards ligands with G protein-coupled receptor activity.
Collapse
Affiliation(s)
- Jacob M Remington
- Department of Chemistry, University of Vermont, Burlington, VT, 05405, USA
| | - Kyle T McKay
- Department of Chemistry, University of Vermont, Burlington, VT, 05405, USA
| | - Noah B Beckage
- Department of Chemistry, University of Vermont, Burlington, VT, 05405, USA
| | - Jonathon B Ferrell
- Department of Chemistry, University of Vermont, Burlington, VT, 05405, USA
| | - Severin T Schneebeli
- Department of Chemistry, University of Vermont, Burlington, VT, 05405, USA.,Department of Industrial and Physical Pharmacy, Department of Chemistry, Purdue University, West Lafayette, IN, 47906, USA.,Department of Pathology, University of Vermont, Burlington, VT, 05405, USA
| | - Jianing Li
- Department of Chemistry, University of Vermont, Burlington, VT, 05405, USA. .,Department of Pathology, University of Vermont, Burlington, VT, 05405, USA. .,Department of Medicinal Chemistry and Molecular Pharmacology, Purdue University, West Lafayette, IN, 47906, USA.
| |
Collapse
|
4
|
Fedorov DG. Parametrized quantum-mechanical approaches combined with the fragment molecular orbital method. J Chem Phys 2022; 157:231001. [PMID: 36550057 DOI: 10.1063/5.0131256] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022] Open
Abstract
Fast parameterized methods such as density-functional tight-binding (DFTB) facilitate realistic calculations of large molecular systems, which can be accelerated by the fragment molecular orbital (FMO) method. Fragmentation facilitates interaction analyses between functional parts of molecular systems. In addition to DFTB, other parameterized methods combined with FMO are also described. Applications of FMO methods to biochemical and inorganic systems are reviewed.
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
|
5
|
The Importance of Charge Transfer and Solvent Screening in the Interactions of Backbones and Functional Groups in Amino Acid Residues and Nucleotides. Int J Mol Sci 2022; 23:ijms232113514. [PMID: 36362296 PMCID: PMC9654426 DOI: 10.3390/ijms232113514] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 10/26/2022] [Accepted: 10/28/2022] [Indexed: 11/06/2022] Open
Abstract
Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated.
Collapse
|
6
|
Fedorov DG. Polarization energies in the fragment molecular orbital method. J Comput Chem 2022; 43:1094-1103. [PMID: 35446441 DOI: 10.1002/jcc.26869] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Revised: 03/25/2022] [Accepted: 04/05/2022] [Indexed: 12/23/2022]
Abstract
Using isolated and polarized states of fragments, a method for computing the polarization energies in density functional theory (DFT) and density-functional tight-binding (DFTB) is developed in the framework of the fragment molecular orbital method. For DFTB, the method is extended into the use of periodic boundary conditions (PBC), for which a new component, a periodic self-polarization energy, is derived. The couplings of the polarization to other components in the pair interaction energy analysis (PIEDA) are derived for DFT and DFTB, and compared to Hartree-Fock and second-order Møller-Plesset perturbation theory (MP2). The effect of the self-consistent (DFT) and perturbative (MP2) treatment of the electron correlation on the polarization is discussed. The difference in the polarization in the bulk (PBC) and micro (cluster) solvation is elucidated.
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
|
7
|
Identification of Novel Natural Product Inhibitors against Matrix Metalloproteinase 9 Using Quantum Mechanical Fragment Molecular Orbital-Based Virtual Screening Methods. Int J Mol Sci 2022; 23:ijms23084438. [PMID: 35457257 PMCID: PMC9030947 DOI: 10.3390/ijms23084438] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 04/11/2022] [Accepted: 04/14/2022] [Indexed: 11/22/2022] Open
Abstract
Matrix metalloproteinases (MMPs) are calcium-dependent zinc-containing endopeptidases involved in multiple cellular processes. Among the MMP isoforms, MMP-9 regulates cancer invasion, rheumatoid arthritis, and osteoarthritis by degrading extracellular matrix proteins present in the tumor microenvironment and cartilage and promoting angiogenesis. Here, we identified two potent natural product inhibitors of the non-catalytic hemopexin domain of MMP-9 using a novel quantum mechanical fragment molecular orbital (FMO)-based virtual screening workflow. The workflow integrates qualitative pharmacophore modeling, quantitative binding affinity prediction, and a raw material search of natural product inhibitors with the BMDMS-NP library. In binding affinity prediction, we made a scoring function with the FMO method and applied the function to two protein targets (acetylcholinesterase and fibroblast growth factor 1 receptor) from DUD-E benchmark sets. In the two targets, the FMO method outperformed the Glide docking score and MM/PBSA methods. By applying this workflow to MMP-9, we proposed two potent natural product inhibitors (laetanine 9 and genkwanin 10) that interact with hotspot residues of the hemopexin domain of MMP-9. Laetanine 9 and genkwanin 10 bind to MMP-9 with a dissociation constant (KD) of 21.6 and 0.614 μM, respectively. Overall, we present laetanine 9 and genkwanin 10 for MMP-9 and demonstrate that the novel FMO-based workflow with a quantum mechanical approach is promising to discover potent natural product inhibitors of MMP-9, satisfying the pharmacophore model and good binding affinity.
Collapse
|
8
|
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
|
9
|
Xu Z, Wauchope OR, Frank AT. Navigating Chemical Space by Interfacing Generative Artificial Intelligence and Molecular Docking. J Chem Inf Model 2021; 61:5589-5600. [PMID: 34633194 DOI: 10.1021/acs.jcim.1c00746] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Here, we report the implementation and application of a simple, structure-aware framework to generate target-specific screening libraries. Our approach combines advances in generative artificial intelligence (AI) with conventional molecular docking to explore chemical space conditioned on the unique physicochemical properties of the active site of a biomolecular target. As a demonstration, we used our framework, which we refer to as sample-and-dock, to construct focused libraries for cyclin-dependent kinase type-2 (CDK2) and the active site of the main protease (Mpro) of the SARS-CoV-2 virus. We envision that the sample-and-dock framework could be used to generate theoretical maps of the chemical space specific to a given target and so provide information about its molecular recognition characteristics.
Collapse
Affiliation(s)
- Ziqiao Xu
- Chemistry Department, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, United States
| | - Orrette R Wauchope
- Department of Natural Sciences, City University of New York, Baruch College, New York, New York 10010, United States
| | - Aaron T Frank
- Biophysics Program, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
10
|
Nakamura T, Yokaichiya T, Fedorov DG. Quantum-Mechanical Structure Optimization of Protein Crystals and Analysis of Interactions in Periodic Systems. J Phys Chem Lett 2021; 12:8757-8762. [PMID: 34478310 DOI: 10.1021/acs.jpclett.1c02510] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
A fast quantum-mechanical approach, density-functional tight-binding combined with the fragment molecular orbital method and periodic boundary conditions, is used to optimize atomic coordinates and cell parameters for a set of protein crystals: 1ETL, 5OQZ, 3Q8J, 1CBN, and 2VB1. Good agreement between experimental and calculated structures is obtained for both atomic coordinates and cell parameters. Sterical clashes present in the experimental structures are corrected by simulations. The partition analysis is extended to treat periodic boundary conditions and applied to analyze protein-solvent interactions in crystals.
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
|
11
|
Abstract
Vibrational energies are partitioned into the contributions of molecular parts called segments, for instance, residues in proteins. The fragment molecular orbital method is used to facilitate vibrational calculations of large systems at the DFTB and HF-3c levels. The vibrational analysis is combined with the partitioning of the electronic energy, yielding free-energy contributions of segments to the binding energy, pinpointing hot spots for drug discovery and other studies. The analysis is illustrated on two protein-ligand complexes in solution.
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
|
12
|
Ma B, Terayama K, Matsumoto S, Isaka Y, Sasakura Y, Iwata H, Araki M, Okuno Y. Structure-Based de Novo Molecular Generator Combined with Artificial Intelligence and Docking Simulations. J Chem Inf Model 2021; 61:3304-3313. [PMID: 34242036 DOI: 10.1021/acs.jcim.1c00679] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Recently, molecular generation models based on deep learning have attracted significant attention in drug discovery. However, most existing molecular generation models have serious limitations in the context of drug design wherein they do not sufficiently consider the effect of the three-dimensional (3D) structure of the target protein in the generation process. In this study, we developed a new deep learning-based molecular generator, SBMolGen, that integrates a recurrent neural network, a Monte Carlo tree search, and docking simulations. The results of an evaluation using four target proteins (two kinases and two G protein-coupled receptors) showed that the generated molecules had a better binding affinity score (docking score) than the known active compounds, and the generated molecules possessed a broader chemical space distribution. SBMolGen not only generates novel binding active molecules but also presents 3D docking poses with target proteins, which will be useful in subsequent drug design. The code is available at https://github.com/clinfo/SBMolGen.
Collapse
Affiliation(s)
- Biao Ma
- Center for Cluster Development and Coordination, Foundation for Biomedical Research and Innovation at Kobe, 1-5-2, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Kei Terayama
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan.,Graduate School of Medical Life Science, Yokohama City University, 1-7-29, Suehiro-cho, Tsurumi-ku, Yohohama, Kanagawa 230-0045, Japan
| | - Shigeyuki Matsumoto
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Yuta Isaka
- Center for Cluster Development and Coordination, Foundation for Biomedical Research and Innovation at Kobe, 1-5-2, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yoko Sasakura
- Center for Cluster Development and Coordination, Foundation for Biomedical Research and Innovation at Kobe, 1-5-2, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Hiroaki Iwata
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Mitsugu Araki
- Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| | - Yasushi Okuno
- Center for Cluster Development and Coordination, Foundation for Biomedical Research and Innovation at Kobe, 1-5-2, Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,RIKEN Center for Computational Science, 7-1-26 Minatojima-minamimachi, Chuo-ku, Kobe, Hyogo 650-0047, Japan.,Graduate School of Medicine, Kyoto University, 53 Shogoin-Kawaharacho, Sakyo-ku, Kyoto 606-8507, Japan
| |
Collapse
|
13
|
Tielker N, Eberlein L, Hessler G, Schmidt KF, Güssregen S, Kast SM. Quantum-mechanical property prediction of solvated drug molecules: what have we learned from a decade of SAMPL blind prediction challenges? J Comput Aided Mol Des 2021; 35:453-472. [PMID: 33079358 PMCID: PMC8018924 DOI: 10.1007/s10822-020-00347-5] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2020] [Accepted: 09/26/2020] [Indexed: 01/26/2023]
Abstract
Joint academic-industrial projects supporting drug discovery are frequently pursued to deploy and benchmark cutting-edge methodical developments from academia in a real-world industrial environment at different scales. The dimensionality of tasks ranges from small molecule physicochemical property assessment over protein-ligand interaction up to statistical analyses of biological data. This way, method development and usability both benefit from insights gained at both ends, when predictiveness and readiness of novel approaches are confirmed, but the pharmaceutical drug makers get early access to novel tools for the quality of drug products and benefit of patients. Quantum-mechanical and simulation methods particularly fall into this group of methods, as they require skills and expense in their development but also significant resources in their application, thus are comparatively slowly dripping into the realm of industrial use. Nevertheless, these physics-based methods are becoming more and more useful. Starting with a general overview of these and in particular quantum-mechanical methods for drug discovery we review a decade-long and ongoing collaboration between Sanofi and the Kast group focused on the application of the embedded cluster reference interaction site model (EC-RISM), a solvation model for quantum chemistry, to study small molecule chemistry in the context of joint participation in several SAMPL (Statistical Assessment of Modeling of Proteins and Ligands) blind prediction challenges. Starting with early application to tautomer equilibria in water (SAMPL2) the methodology was further developed to allow for challenge contributions related to predictions of distribution coefficients (SAMPL5) and acidity constants (SAMPL6) over the years. Particular emphasis is put on a frequently overlooked aspect of measuring the quality of models, namely the retrospective analysis of earlier datasets and predictions in light of more recent and advanced developments. We therefore demonstrate the performance of the current methodical state of the art as developed and optimized for the SAMPL6 pKa and octanol-water log P challenges when re-applied to the earlier SAMPL5 cyclohexane-water log D and SAMPL2 tautomer equilibria datasets. Systematic improvement is not consistently found throughout despite the similarity of the problem class, i.e. protonation reactions and phase distribution. Hence, it is possible to learn about hidden bias in model assessment, as results derived from more elaborate methods do not necessarily improve quantitative agreement. This indicates the role of chance or coincidence for model development on the one hand which allows for the identification of systematic error and opportunities toward improvement and reveals possible sources of experimental uncertainty on the other. These insights are particularly useful for further academia-industry collaborations, as both partners are then enabled to optimize both the computational and experimental settings for data generation.
Collapse
Affiliation(s)
- Nicolas Tielker
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany
| | - Gerhard Hessler
- R&D Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany
| | - K Friedemann Schmidt
- R&D Preclinical Safety, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany
| | - Stefan Güssregen
- R&D Integrated Drug Discovery, Sanofi-Aventis Deutschland GmbH, 65926, Frankfurt am Main, Germany.
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Otto-Hahn-Str. 4a, 44227, Dortmund, Germany.
| |
Collapse
|
14
|
Nishimoto Y, Fedorov DG. The fragment molecular orbital method combined with density-functional tight-binding and periodic boundary conditions. J Chem Phys 2021; 154:111102. [PMID: 33752370 DOI: 10.1063/5.0039520] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The density-functional tight-binding (DFTB) formulation of the fragment molecular orbital method is combined with periodic boundary conditions. Long-range electrostatics and dispersion are evaluated with the Ewald summation technique. The first analytic derivatives of the energy with respect to atomic coordinates and lattice parameters are formulated. The accuracy of the method is established in comparison to numerical gradients and DFTB without fragmentation. The largest elementary cell in this work has 1631 atoms. The method is applied to elucidate the polarization, charge transfer, and interactions in the solution.
Collapse
Affiliation(s)
- Yoshio Nishimoto
- Graduate School of Science, Kyoto University, Kitashirakawa Oiwakecho, Sakyoku, Kyoto 606-8502, Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
15
|
Abstract
Computational methods for modeling biochemical processes implemented in GAMESS package are reviewed; in particular, quantum mechanics combined with molecular mechanics (QM/MM), semi-empirical, and fragmentation approaches. A detailed summary of capabilities is provided for the QM/MM implementation in QuanPol program and the fragment molecular orbital (FMO) method. Molecular modeling and visualization packages useful for biochemical simulations with GAMESS are described. GAMESS capabilities with corresponding references are tabulated for reader's convenience.
Collapse
|
16
|
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
|
17
|
Abstract
High-order charge transfer is incorporated into the fragment molecular orbital (FMO) method using a charge transfer state with fractional charges. This state is used for a partition analysis of properties based on segments that may be different from fragments in FMO. The partition analysis is also formulated for calculations without fragmentation. All development in this work is limited to density-functional tight-binding. The analysis is applied to a water cluster, crambin (PDB: 1CBN), and two complexes of Trp-cage (1L2Y) with ligands. The contributions of functional groups in ligands are obtained, providing useful information for drug discovery.
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
|
18
|
Cortés-García CJ, Chacón-García L, Mejía-Benavides JE, Díaz-Cervantes E. Tackling the SARS-CoV-2 main protease using hybrid derivatives of 1,5-disubstituted tetrazole-1,2,3-triazoles: an in silico assay. PEERJ PHYSICAL CHEMISTRY 2020. [DOI: 10.7717/peerj-pchem.10] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
In regard to the actual public health global emergency and, based on the state of the art about the ways to inhibit the SARS-CoV-2 treating the COVID19, a family of 1,5-disubstituted tetrazole-1,2,3-triazoles, previously synthesized, have been evaluated through in silico assays against the main protease of the mentioned virus (CoV-2-MPro). The results show that three of these compounds present a more favorable interaction with the selected target than the co-crystallized molecule, which is a peptide-like derivative. It was also found that also hydrophobic interactions play a key role in the ligand-target molecular couplings, due to the higher hydrophobic surfaces into the active site. Finally, a pharmacophore model has been proposed based on the results below, and a family of 1,5-DT derivatives has been designed and tested with the same methods employed in this work. It was concluded that the compound with the isatin as a substituent (P8) present the higher ligand-target interaction, which makes this a strong drug candidate against COVID19, due can inhibit the CoV-2-MProprotein.
Collapse
Affiliation(s)
- Carlos J. Cortés-García
- Laboratorio de Diseño Molecular, Instituto de Investigaciones Químico-Biológicas, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, México
| | - Luis Chacón-García
- Laboratorio de Diseño Molecular, Instituto de Investigaciones Químico-Biológicas, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, México
| | - Jorge Emmanuel Mejía-Benavides
- Departamento de Enfermería y Obstetricia, Centro Interdisciplinario del Noreste (CINUG), Universidad de Guanajuato, Tierra Blanca, Guanajuato, México
| | - Erik Díaz-Cervantes
- Departamento de Alimentos, Centro Interdisciplinario del Noreste (CINUG), Universidad de Guanajuato, Tierra Blanca, Guanajuato, México
| |
Collapse
|
19
|
Nakliang P, Lazim R, Chang H, Choi S. Multiscale Molecular Modeling in G Protein-Coupled Receptor (GPCR)-Ligand Studies. Biomolecules 2020; 10:E631. [PMID: 32325877 PMCID: PMC7226129 DOI: 10.3390/biom10040631] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Revised: 04/13/2020] [Accepted: 04/17/2020] [Indexed: 12/17/2022] Open
Abstract
G protein-coupled receptors (GPCRs) are major drug targets due to their ability to facilitate signal transduction across cell membranes, a process that is vital for many physiological functions to occur. The development of computational technology provides modern tools that permit accurate studies of the structures and properties of large chemical systems, such as enzymes and GPCRs, at the molecular level. The advent of multiscale molecular modeling permits the implementation of multiple levels of theories on a system of interest, for instance, assigning chemically relevant regions to high quantum mechanics (QM) level of theory while treating the rest of the system using classical force field (molecular mechanics (MM) potential). Multiscale QM/MM molecular modeling have far-reaching applications in the rational design of GPCR drugs/ligands by affording precise ligand binding configurations through the consideration of conformational plasticity. This enables the identification of key binding site residues that could be targeted to manipulate GPCR function. This review will focus on recent applications of multiscale QM/MM molecular simulations in GPCR studies that could boost the efficiency of future structure-based drug design (SBDD) strategies.
Collapse
Affiliation(s)
| | | | | | - Sun Choi
- College of Pharmacy and Graduate School of Pharmaceutical Sciences, Ewha Womans University, Seoul 03760, Korea; (P.N.); (R.L.); (H.C.)
| |
Collapse
|
20
|
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
|
21
|
Kusumoto M, Ueno-Noto K, Takano K. Systematic Interaction Analysis of Anti-Human Immunodeficiency Virus Type-1 Neutralizing Antibodies with High Mannose Glycans Using Fragment Molecular Orbital and Molecular Dynamics Methods. J Comput Chem 2020; 41:31-42. [PMID: 31565801 DOI: 10.1002/jcc.26073] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2019] [Revised: 08/18/2019] [Accepted: 09/01/2019] [Indexed: 11/10/2022]
Abstract
A series of broadly neutralizing antibodies called PGT have been shown to be bound directly to human immunodeficiency virus type-1 via high mannose glycans on glycoprotein gp120. Despite the sequence similarities of amino acids of the antibodies, their affinities to the glycan differ. Glycan-antibody interactions among these antibodies are systematically compared with quantum chemical fragment molecular orbital calculations and molecular dynamics simulations. The differences among structural stability of the glycan in the active site of the complexes and total interaction energies as well as binding free energies between the glycan and antibodies agree well with the experimentally shown affinities of the glycan to the antibodies. The terminal saccharide, Man D3, is structurally stable and responsible for the glycan-antibody binding through electrostatic and dispersion interactions. The structural stability of nonterminal saccharides such as Man 4 or Man C plays substantial roles in the interaction via direct hydrogen bonds. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Miyu Kusumoto
- Department of Chemistry and Biochemistry, Graduate School of Humanities and Sciences, Ochanomizu University, 2-1-1 Otsuka, Bunkyo-ku, Tokyo, 112-8610, Japan
| | - Kaori Ueno-Noto
- Center for Natural Sciences, College of Liberal Arts and Sciences, Kitasato University, 1-15-1 Kitasato, Minami-ku, Sagamihara, Kanagawa, 252-0373, Japan
| | - Keiko Takano
- Department of Chemistry and Biochemistry, Graduate School of Humanities and Sciences, Ochanomizu University, 2-1-1 Otsuka, Bunkyo-ku, Tokyo, 112-8610, Japan
| |
Collapse
|
22
|
Morao I, Heifetz A, Fedorov DG. Accurate Scoring in Seconds with the Fragment Molecular Orbital and Density-Functional Tight-Binding Methods. Methods Mol Biol 2020; 2114:143-148. [PMID: 32016891 DOI: 10.1007/978-1-0716-0282-9_9] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The accurate evaluation of receptor-ligand interactions is an essential part of rational drug design. While quantum mechanical (QM) methods have been a promising means by which to achieve this, traditional QM is not applicable for large biological systems due to its high computational cost. Here, the fragment molecular orbital (FMO) method has been combined with the density-functional tight-binding (DFTB) method to compute energy calculations of biological systems in seconds. FMO-DFTB outperformed GBVI/WSA in identifying a set of 10 binders versus a background of 500 decoys applied to human k-opioid receptor. The significant increase in the speed and the high accuracy achieved with FMO-DFTB calculations allows FMO to be applied in areas of drug discovery that were not previously accessible to traditional QM methodologies. For the first time, it is now possible to perform FMO calculations in a high-throughput manner.
Collapse
Affiliation(s)
- Inaki Morao
- Evotec (UK) Ltd., Abingdon, Oxfordshire, UK.
| | | | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
| |
Collapse
|
23
|
Abstract
Basic concepts in the analysis of binding using the fragment molecular orbital method are discussed at length: polarization, desolvation, and interaction. The components in the pair interaction energy decomposition analysis are introduced, and the analysis is illustrated for a water dimer and a protein-ligand complex.
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
|
24
|
Abstract
This chapter describes the current status of development of the fragment molecular orbital (FMO) method for analyzing the electronic state and intermolecular interactions of biomolecular systems in solvent. The orbital energies and the inter-fragment interaction energies (IFIEs) for a specific molecular structure can be obtained directly by performing FMO calculations by exposing water molecules and counterions around biomolecular systems. Then, it is necessary to pay attention to the thickness of the water shell surrounding the biomolecules. The single-point calculation for snapshots from MD trajectory does not incorporate the effects of temperature and configurational fluctuation, but the SCIFIE (statistically corrected IFIE) method is proposed as a many-body correlated method that partially compensates for this deficiency. Furthermore, implicit continuous dielectric models have been developed as effective approaches to incorporating the screening effect of the solvent in thermal equilibrium, and we illustrate their usefulness for theoretical evaluation of IFIEs and ligand-binding free energy on the basis of the FMO-PBSA (Poisson-Boltzmann surface area) method and other computational methods.
Collapse
|
25
|
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: 22] [Impact Index Per Article: 4.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
|
26
|
Townsend-Nicholson A, Altwaijry N, Potterton A, Morao I, Heifetz A. Computational prediction of GPCR oligomerization. Curr Opin Struct Biol 2019; 55:178-184. [PMID: 31170578 DOI: 10.1016/j.sbi.2019.04.005] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2019] [Revised: 04/16/2019] [Accepted: 04/17/2019] [Indexed: 01/08/2023]
Abstract
There has been a recent and prolific expansion in the number of GPCR crystal structures being solved: in both active and inactive forms and in complex with ligand, with G protein and with each other. Despite this, there is relatively little experimental information about the precise configuration of GPCR oligomers during these different biologically relevant states. While it may be possible to identify the experimental conditions necessary to crystallize a GPCR preferentially in a specific structural conformation, computational approaches afford a potentially more tractable means of describing the probability of formation of receptor dimers and higher order oligomers. Ensemble-based computational methods based on structurally determined dimers, coupled with a computational workflow that uses quantum mechanical methods to analyze the chemical nature of the molecular interactions at a GPCR dimer interface, will generate the reproducible and accurate predictions needed to predict previously unidentified GPCR dimers and to inform future advances in our ability to understand and begin to precisely manipulate GPCR oligomers in biological systems. It may also provide information needed to achieve an increase in the number of experimentally determined oligomeric GPCR structures.
Collapse
Affiliation(s)
- 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.
| | - Nojood Altwaijry
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom
| | - Andrew Potterton
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom; Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Inaki Morao
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| |
Collapse
|
27
|
Vuong VQ, Nishimoto Y, Fedorov DG, Sumpter BG, Niehaus TA, Irle S. The Fragment Molecular Orbital Method Based on Long-Range Corrected Density-Functional Tight-Binding. J Chem Theory Comput 2019; 15:3008-3020. [PMID: 30998360 DOI: 10.1021/acs.jctc.9b00108] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The presently available linear scaling approaches to density-functional tight-binding (DFTB) based on the fragment molecular orbital (FMO) method are severely impacted by the problem of artificial charge transfer due to the self-interaction error (SIE), which hampers the simulation of zwitterionic systems such as biopolymers or ionic liquids. Here we report an extension of FMO-DFTB where we included a long-range corrected (LC) functional designed to mitigate the DFTB SIE, called the FMO-LC-DFTB method, resulting in a robust method which succeeds in simulating zwitterionic systems. Both energy and analytic gradient are developed for the gas phase and the polarizable continuum model of solvation. The scaling of FMO-LC-DFTB with system size N is shown to be almost linear, O( N1.13-1.28), and its numerical accuracy is established for a variety of representative systems including neutral and charged polypeptides. It is shown that pair interaction energies between fragments for two mini-proteins are in excellent agreement with results from long-range corrected density functional theory. The new method was employed in long time scale (1 ns) molecular dynamics simulations of the tryptophan cage protein (PDB: 1L2Y ) in the gas phase for four different protonation states and in stochastic global minimum structure searches for 1-ethyl-3-methylimidazolium nitrate ionic liquid clusters containing up to 2300 atoms.
Collapse
Affiliation(s)
- Van Quan Vuong
- Bredesen Center for Interdisciplinary Research and Graduate Education , University of Tennessee , Knoxville , Tennessee 37996 , United States
| | - Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry , Kyoto University , Kyoto 606-8501 , Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Tsukuba 305-8568 , Japan
| | - Bobby G Sumpter
- Center for Nanophase Materials Sciences and Computational Sciences and Engineering Division , Oak Ridge National Laboratory , Oak Ridge , Tennessee 37831 , United States
| | - Thomas A Niehaus
- Univ Lyon, Université Claude Bernard Lyon 1 , CNRS, Institut Lumière Matière , F-69622 Villeurbanne , France
| | - Stephan Irle
- Bredesen Center for Interdisciplinary Research and Graduate Education , University of Tennessee , Knoxville , Tennessee 37996 , United States.,Center for Nanophase Materials Sciences and Computational Sciences and Engineering Division , Oak Ridge National Laboratory , Oak Ridge , Tennessee 37831 , United States.,Chemical Sciences Division , Oak Ridge National Laboratory , Oak Ridge , Tennessee 37831 , United States
| |
Collapse
|
28
|
Heifetz A, James T, Southey M, Morao I, Aldeghi M, Sarrat L, Fedorov DG, Bodkin MJ, Townsend-Nicholson A. Characterising GPCR-ligand interactions using a fragment molecular orbital-based approach. Curr Opin Struct Biol 2019; 55:85-92. [PMID: 31022570 DOI: 10.1016/j.sbi.2019.03.021] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2019] [Revised: 02/19/2019] [Accepted: 03/14/2019] [Indexed: 10/27/2022]
Abstract
There has been fantastic progress in solving GPCR crystal structures. 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 mechanical approaches (QM) are often too computationally expensive, 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 towards ligand binding, including an analysis of their chemical nature.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom.
| | - Tim James
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Michelle Southey
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Inaki Morao
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Matteo Aldeghi
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Laurie Sarrat
- Evotec (France) SAS, 195 Route d' Espagne, 31036 Toulouse, France
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (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, 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
|
29
|
Śliwa P, Kurczab R, Kafel R, Drabczyk A, Jaśkowska J. Recognition of repulsive and attractive regions of selected serotonin receptor binding site using FMO-EDA approach. J Mol Model 2019; 25:114. [PMID: 30955095 DOI: 10.1007/s00894-019-3995-6] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2018] [Accepted: 03/14/2019] [Indexed: 12/28/2022]
Abstract
The complexes of selected long-chain arylpiperazines with homology models of 5-HT1A, 5-HT2A, and 5-HT7 receptors were investigated using quantum mechanical methods. The molecular geometries of the ligand-receptor complexes were firstly optimized with the Our own N-layered Integrated molecular Orbital and molecular Mechanics (ONIOM) method. Next, the fragment molecular orbitals method with an energy decomposition analysis scheme (FMO-EDA) was employed to estimate the interaction energies in binding sites. The results clearly showed that orthosteric binding sites of studied serotonin receptors have both attractive and repulsive regions. In the case of 5-HT1A and 5-HT2A two repulsive areas, located in the lower part of the binding pocket, and one large area of attraction engaging many residues at the top of all helices were identified. Additionally, for the 5-HT7 receptor, the third area of destabilization located at the extracellular end of the helix 6 was found.
Collapse
Affiliation(s)
- Paweł Śliwa
- Faculty of Chemical Engineering and Technology, Cracow University of Technology, 24 Warszawska, 31-155, Kraków, Poland.
| | - Rafał Kurczab
- Department of Medicinal Chemistry, Institute of Pharmacology, Polish Academy of Sciences, 12 Smȩtna, 31-343, Kraków, Poland
| | - Rafał Kafel
- Department of Medicinal Chemistry, Institute of Pharmacology, Polish Academy of Sciences, 12 Smȩtna, 31-343, Kraków, Poland
| | - Anna Drabczyk
- Faculty of Chemical Engineering and Technology, Cracow University of Technology, 24 Warszawska, 31-155, Kraków, Poland
| | - Jolanta Jaśkowska
- Faculty of Chemical Engineering and Technology, Cracow University of Technology, 24 Warszawska, 31-155, Kraków, Poland
| |
Collapse
|
30
|
Watanabe C, Watanabe H, Okiyama Y, Takaya D, Fukuzawa K, Tanaka S, Honma T. Development of an automated fragment molecular orbital (FMO) calculation protocol toward construction of quantum mechanical calculation database for large biomolecules . CHEM-BIO INFORMATICS JOURNAL 2019. [DOI: 10.1273/cbij.19.5] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Affiliation(s)
| | - Hirofumi Watanabe
- Education Center on Computational Science and Engineering, Kobe University
| | - Yoshio Okiyama
- Center for Biosystems Dynamics Research, RIKEN
- National Institute of Health Sciences
| | | | - Kaori Fukuzawa
- Center for Biosystems Dynamics Research, RIKEN
- Department of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University
| | | | | |
Collapse
|
31
|
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
|
32
|
Okiyama Y, Watanabe C, Fukuzawa K, Mochizuki Y, Nakano T, Tanaka S. Fragment Molecular Orbital Calculations with Implicit Solvent Based on the Poisson-Boltzmann Equation: II. Protein and Its Ligand-Binding System Studies. J Phys Chem B 2018; 123:957-973. [PMID: 30532968 DOI: 10.1021/acs.jpcb.8b09326] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this study, the electronic properties of bioactive proteins were analyzed using an ab initio fragment molecular orbital (FMO) methodology in solution: coupling with an implicit solvent model based on the Poisson-Boltzmann surface area called as FMO-PBSA. We investigated the solvent effects on practical and heterogeneous targets with uneven exposure to solvents unlike deoxyribonucleic acid analyzed in our recent study. Interfragment interaction energy (IFIE) and its decomposition analyses by FMO-PBSA revealed solvent-screening mechanisms that affect local stability inside ubiquitin protein: the screening suppresses excessiveness in bare charge-charge interactions and enables an intuitive IFIE analysis. The electrostatic character and associated solvation free energy also give consistent results as a whole to previous studies on the explicit solvent model. Moreover, by using the estrogen receptor alpha (ERα) protein bound to ligands, we elucidated the importance of specific interactions that depend on the electric charge and activatability as agonism/antagonism of the ligand while estimating the influences of the implicit solvent on the ligand and helix-12 bindings. The predicted ligand-binding affinities of bioactive compounds to ERα also show a good correlation with their in vitro activities. The FMO-PBSA approach would thus be a promising tool both for biological and pharmaceutical research targeting proteins.
Collapse
Affiliation(s)
- Yoshio Okiyama
- Institute of Industrial Science , The University of Tokyo , 4-6-1 Komaba , Meguro-ku, Tokyo 153-8505 , Japan.,Division of Medicinal Safety Science , National Institute of Health Sciences , 3-25-26 Tonomachi , Kawasaki-ku, Kawasaki , Kanagawa 210-9501 , Japan
| | - Chiduru Watanabe
- Institute of Industrial Science , The University of Tokyo , 4-6-1 Komaba , Meguro-ku, Tokyo 153-8505 , Japan.,RIKEN Center for Biosystems Dynamics Research , 1-7-22 Suehiro-cho , Tsurumi-ku, Yokohama , Kanagawa 230-0045 , Japan
| | - Kaori Fukuzawa
- Institute of Industrial Science , The University of Tokyo , 4-6-1 Komaba , Meguro-ku, Tokyo 153-8505 , Japan.,Faculty of Pharmaceutical Sciences , Hoshi University , 2-4-41 Ebara , Shinagawa-ku, Tokyo 142-8501 , Japan
| | - Yuji Mochizuki
- Institute of Industrial Science , The University of Tokyo , 4-6-1 Komaba , Meguro-ku, Tokyo 153-8505 , Japan.,Department of Chemistry and Research Center for Smart Molecules, Faculty of Science , Rikkyo University , 3-34-1 Nishi-ikebukuro , Toshima-ku, Tokyo 171-8501 , Japan
| | - Tatsuya Nakano
- Institute of Industrial Science , The University of Tokyo , 4-6-1 Komaba , Meguro-ku, Tokyo 153-8505 , Japan.,Division of Medicinal Safety Science , National Institute of Health Sciences , 3-25-26 Tonomachi , Kawasaki-ku, Kawasaki , Kanagawa 210-9501 , Japan
| | - Shigenori Tanaka
- Graduate School of System Informatics , Kobe University , 1-1 Rokkodai, Nada-ku, Kobe , Hyogo 657-8501 , Japan
| |
Collapse
|
33
|
Nakata H, Fedorov DG. Analytic second derivatives for the efficient electrostatic embedding in the fragment molecular orbital method. J Comput Chem 2018; 39:2039-2050. [DOI: 10.1002/jcc.25360] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2018] [Revised: 04/27/2018] [Accepted: 04/29/2018] [Indexed: 01/09/2023]
Affiliation(s)
- Hiroya Nakata
- Department of Fundamental Technology Research; Research and Development Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho; Kirishima-shi Kagoshima, 899-4312 Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono; Tsukuba Ibaraki, 305-8568 Japan
| |
Collapse
|
34
|
Empirical corrections and pair interaction energies in the fragment molecular orbital method. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.06.025] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
35
|
Nishimoto Y, Fedorov DG. Adaptive frozen orbital treatment for the fragment molecular orbital method combined with density-functional tight-binding. J Chem Phys 2018; 148:064115. [DOI: 10.1063/1.5012935] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
36
|
Fedorov DG, Kitaura K. Pair Interaction Energy Decomposition Analysis for Density Functional Theory and Density-Functional Tight-Binding with an Evaluation of Energy Fluctuations in Molecular Dynamics. J Phys Chem A 2018; 122:1781-1795. [DOI: 10.1021/acs.jpca.7b12000] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/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
| | - Kazuo Kitaura
- Advanced
Institute for Computational Science (AICS), RIKEN, 7-1-26 Minatojima-Minami-Machi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Fukui
Institute for Fundamental Chemistry, Kyoto University, Takano-Nishihiraki-cho
34-4, Sakyou-ku, Kyoto 606-8103, Japan
| |
Collapse
|
37
|
Heifetz A, Southey M, Morao I, Townsend-Nicholson A, Bodkin MJ. Computational Methods Used in Hit-to-Lead and Lead Optimization Stages of Structure-Based Drug Discovery. Methods Mol Biol 2018; 1705:375-394. [PMID: 29188574 DOI: 10.1007/978-1-4939-7465-8_19] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
GPCR modeling approaches are widely used in the hit-to-lead (H2L) and lead optimization (LO) stages of drug discovery. The aims of these modeling approaches are to predict the 3D structures of the receptor-ligand complexes, to explore the key interactions between the receptor and the ligand and to utilize these insights in the design of new molecules with improved binding, selectivity or other pharmacological properties. In this book chapter, we present a brief survey of key computational approaches integrated with hierarchical GPCR modeling protocol (HGMP) used in hit-to-lead (H2L) and in lead optimization (LO) stages of structure-based drug discovery (SBDD). We outline the differences in modeling strategies used in H2L and LO of SBDD and illustrate how these tools have been applied in three drug discovery projects.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK. .,Division of Biosciences, Research Department of Structural and Molecular Biology, Institute of Structural and Molecular Biology, University College London, London, WC1E 6BT, UK.
| | - Michelle Southey
- 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
| | - Andrea Townsend-Nicholson
- Division of Biosciences, Research Department of Structural and Molecular Biology, University College London, Darwin Building, Gower Street, London, WC1E 6BT, UK
| | - Mike J Bodkin
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| |
Collapse
|
38
|
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
|
39
|
Morao I, Fedorov DG, Robinson R, Southey M, Townsend‐Nicholson A, Bodkin MJ, Heifetz A. Rapid and accurate assessment of GPCR-ligand interactions Using the fragment molecular orbital-based density-functional tight-binding method. J Comput Chem 2017; 38:1987-1990. [PMID: 28675443 PMCID: PMC5600120 DOI: 10.1002/jcc.24850] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2017] [Revised: 05/07/2017] [Accepted: 05/16/2017] [Indexed: 01/17/2023]
Abstract
The reliable and precise evaluation of receptor–ligand interactions and pair‐interaction energy is an essential element of rational drug design. While quantum mechanical (QM) methods have been a promising means by which to achieve this, traditional QM is not applicable for large biological systems due to its high computational cost. Here, the fragment molecular orbital (FMO) method has been used to accelerate QM calculations, and by combining FMO with the density‐functional tight‐binding (DFTB) method we are able to decrease computational cost 1000 times, achieving results in seconds, instead of hours. We have applied FMO‐DFTB to three different GPCR–ligand systems. Our results correlate well with site directed mutagenesis data and findings presented in the published literature, demonstrating that FMO‐DFTB is a rapid and accurate means of GPCR–ligand interactions. © 2017 Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Inaki Morao
- Computational ChemistryEvotec (UK) Ltd114 Innovation Drive, Milton ParkAbingdonOxfordshireOX14 4RZUnited Kingdom
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat)National Institute of Advanced Industrial Science and Technology (AIST)1‐1‐1 UmezonoTsukubaIbaraki305‐8568Japan
| | - Roger Robinson
- Computational ChemistryEvotec (UK) Ltd114 Innovation Drive, Milton ParkAbingdonOxfordshireOX14 4RZUnited Kingdom
| | - Michelle Southey
- Computational ChemistryEvotec (UK) Ltd114 Innovation Drive, Milton ParkAbingdonOxfordshireOX14 4RZUnited Kingdom
| | - Andrea Townsend‐Nicholson
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of BiosciencesUniversity College LondonLondonWC1E 6BTUnited Kingdom
| | - Mike J. Bodkin
- Computational ChemistryEvotec (UK) Ltd114 Innovation Drive, Milton ParkAbingdonOxfordshireOX14 4RZUnited Kingdom
| | - Alexander Heifetz
- Computational ChemistryEvotec (UK) Ltd114 Innovation Drive, Milton ParkAbingdonOxfordshireOX14 4RZUnited Kingdom
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of BiosciencesUniversity College LondonLondonWC1E 6BTUnited Kingdom
| |
Collapse
|
40
|
Fedorov DG. The fragment molecular orbital method: theoretical development, implementation in
GAMESS
, and applications. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1322] [Citation(s) in RCA: 71] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/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)TsukubaJapan
| |
Collapse
|