1
|
Handa Y, Okuwaki K, Kawashima Y, Hatada R, Mochizuki Y, Komeiji Y, Tanaka S, Furuishi T, Yonemochi E, Honma T, Fukuzawa K. Prediction of Binding Pose and Affinity of Nelfinavir, a SARS-CoV-2 Main Protease Repositioned Drug, by Combining Docking, Molecular Dynamics, and Fragment Molecular Orbital Calculations. J Phys Chem B 2024; 128:2249-2265. [PMID: 38437183 PMCID: PMC10946393 DOI: 10.1021/acs.jpcb.3c05564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 02/09/2024] [Accepted: 02/14/2024] [Indexed: 03/06/2024]
Abstract
A novel in silico drug design procedure is described targeting the Main protease (Mpro) of the SARS-CoV-2 virus. The procedure combines molecular docking, molecular dynamics (MD), and fragment molecular orbital (FMO) calculations. The binding structure and properties of Mpro were predicted for Nelfinavir (NFV), which had been identified as a candidate compound through drug repositioning, targeting Mpro. Several poses of the Mpro and NFV complexes were generated by docking, from which four docking poses were selected by scoring with FMO energy. Then, each pose was subjected to MD simulation, 100 snapshot structures were sampled from each of the generated MD trajectories, and the structures were evaluated by FMO calculations to rank the pose based on binding energy. Several residues were found to be important in ligand recognition, including Glu47, Asp48, Glu166, Asp187, and Gln189, all of which interacted strongly with NFV. Asn142 is presumably regarded to form hydrogen bonds or CH/π interaction with NFV; however, in the present calculation, their interactions were transient. Moreover, the tert-butyl group of NFV had no interaction with Mpro. Identifying such strong and weak interactions provides candidates for maintaining and substituting ligand functional groups and important suggestions for drug discovery using drug repositioning. Besides the interaction between NFV and the amino acid residues of Mpro, the desolvation effect of the binding pocket also affected the ranking order. A similar procedure of drug design was applied to Lopinavir, and the calculated interaction energy and experimental inhibitory activity value trends were consistent. Our approach provides a new guideline for structure-based drug design starting from a candidate compound whose complex crystal structure has not been obtained.
Collapse
Affiliation(s)
- Yuma Handa
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Koji Okuwaki
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, 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
| | - Yusuke Kawashima
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Ryo Hatada
- 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
| | - Yuji Mochizuki
- 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
- Institute
of Industrial Science, University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Yuto Komeiji
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, 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
- Health
and Medical Research Institute, AIST, Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
- RIKEN
Center
for Biosystems Dynamics Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Shigenori Tanaka
- Graduate
School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Takayuki Furuishi
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Etsuo Yonemochi
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| | - Teruki Honma
- RIKEN
Center
for Biosystems Dynamics Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Kaori Fukuzawa
- Department
of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
- Graduate
School of Pharmaceutical Sciences, Osaka
University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
- Department
of Biomolecular Engineering, Graduate School of Engineering, Tohoku University, 6-6-11 Aoba, Aramaki, Aoba-ku, Sendai 980-8579, Japan
| |
Collapse
|
2
|
Lee AJ, Rackers JA, Pathak S, Bricker WP. Building an ab initio solvated DNA model using Euclidean neural networks. PLoS One 2024; 19:e0297502. [PMID: 38358990 PMCID: PMC10868815 DOI: 10.1371/journal.pone.0297502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2023] [Accepted: 01/06/2024] [Indexed: 02/17/2024] Open
Abstract
Accurately modeling large biomolecules such as DNA from first principles is fundamentally challenging due to the steep computational scaling of ab initio quantum chemistry methods. This limitation becomes even more prominent when modeling biomolecules in solution due to the need to include large numbers of solvent molecules. We present a machine-learned electron density model based on a Euclidean neural network framework that includes a built-in understanding of equivariance to model explicitly solvated double-stranded DNA. By training the machine learning model using molecular fragments that sample the key DNA and solvent interactions, we show that the model predicts electron densities of arbitrary systems of solvated DNA accurately, resolves polarization effects that are neglected by classical force fields, and captures the physics of the DNA-solvent interaction at the ab initio level.
Collapse
Affiliation(s)
- Alex J. Lee
- Department of Chemical and Biological Engineering, University of New Mexico, Albuquerque, NM, United States of America
| | - Joshua A. Rackers
- Center for Computing Research, Sandia National Laboratories, Albuquerque, NM, United States of America
| | - Shivesh Pathak
- Center for Computing Research, Sandia National Laboratories, Albuquerque, NM, United States of America
| | - William P. Bricker
- Department of Chemical and Biological Engineering, University of New Mexico, Albuquerque, NM, United States of America
| |
Collapse
|
3
|
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
|
4
|
Takaba K, Watanabe C, Tokuhisa A, Akinaga Y, Ma B, Kanada R, Araki M, Okuno Y, Kawashima Y, Moriwaki H, Kawashita N, Honma T, Fukuzawa K, Tanaka S. Protein-ligand binding affinity prediction of cyclin-dependent kinase-2 inhibitors by dynamically averaged fragment molecular orbital-based interaction energy. J Comput Chem 2022; 43:1362-1371. [PMID: 35678372 DOI: 10.1002/jcc.26940] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2022] [Revised: 03/27/2022] [Accepted: 04/02/2022] [Indexed: 01/16/2023]
Abstract
Fragment molecular orbital (FMO) method is a powerful computational tool for structure-based drug design, in which protein-ligand interactions can be described by the inter-fragment interaction energy (IFIE) and its pair interaction energy decomposition analysis (PIEDA). Here, we introduced a dynamically averaged (DA) FMO-based approach in which molecular dynamics simulations were used to generate multiple protein-ligand complex structures for FMO calculations. To assess this approach, we examined the correlation between the experimental binding free energies and DA-IFIEs of six CDK2 inhibitors whose net charges are zero. The correlation between the experimental binding free energies and snapshot IFIEs for X-ray crystal structures was R2 = 0.75. Using the DA-IFIEs, the correlation significantly improved to 0.99. When an additional CDK2 inhibitor with net charge of -1 was added, the DA FMO-based scheme with the dispersion energies still achieved R2 = 0.99, whereas R2 decreased to 0.32 employing all the energy terms of PIEDA.
Collapse
Affiliation(s)
- Kenichiro Takaba
- Pharmaceutical Research Center, Advanced Drug Discovery, Asahi Kasei Pharma Corporation, Shizuoka, Japan
| | - Chiduru Watanabe
- RIKEN Center for Biosystems Dynamics Research, Tsurumi-ku, Yokohama, Kanagawa, Japan
| | - Atsushi Tokuhisa
- RIKEN Center for Computational Science, Chuo-ku, Kobe, Hyogo, Japan
| | - Yoshinobu Akinaga
- RIKEN Center for Computational Science, Chuo-ku, Kobe, Hyogo, Japan.,Project Development Department, VINAS Co., Ltd., Osaka, Japan
| | - Biao Ma
- RIKEN Center for Computational Science, Chuo-ku, Kobe, Hyogo, Japan
| | - Ryo Kanada
- RIKEN Center for Computational Science, Chuo-ku, Kobe, Hyogo, Japan
| | - Mitsugu Araki
- Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Yasushi Okuno
- RIKEN Center for Computational Science, Chuo-ku, Kobe, Hyogo, Japan.,Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Yusuke Kawashima
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University, Tokyo, Japan
| | - Hirotomo Moriwaki
- RIKEN Center for Biosystems Dynamics Research, Tsurumi-ku, Yokohama, Kanagawa, Japan
| | | | - Teruki Honma
- RIKEN Center for Biosystems Dynamics Research, Tsurumi-ku, Yokohama, Kanagawa, Japan
| | - Kaori Fukuzawa
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University, Tokyo, Japan.,Department of Biomolecular Engineering, Graduate School of Engineering, Tohoku University, Sendai, Japan.,Graduate School of Pharmaceutical Sciences, Osaka University, Suita, Osaka, Japan
| | - Shigenori Tanaka
- Graduate School of System Informatics, Kobe University, Kobe, Japan
| |
Collapse
|
5
|
Tanaka S, Tokutomi S, Hatada R, Okuwaki K, Akisawa K, Fukuzawa K, Komeiji Y, Okiyama Y, Mochizuki Y. Dynamic Cooperativity of Ligand-Residue Interactions Evaluated with the Fragment Molecular Orbital Method. J Phys Chem B 2021; 125:6501-6512. [PMID: 34124906 DOI: 10.1021/acs.jpcb.1c03043] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
By the splendid advance in computation power realized with the Fugaku supercomputer, it has become possible to perform ab initio fragment molecular orbital (FMO) calculations for thousands of dynamic structures of protein-ligand complexes in a parallel way. We thus carried out electron-correlated FMO calculations for a complex of the 3C-like (3CL) main protease (Mpro) of the new coronavirus (SARS-CoV-2) and its inhibitor N3 incorporating the structural fluctuations sampled by classical molecular dynamics (MD) simulation in hydrated conditions. Along with a statistical evaluation of the interfragment interaction energies (IFIEs) between the N3 ligand and the surrounding amino-acid residues for 1000 dynamic structure samples, in this study we applied a novel approach based on principal component analysis (PCA) and singular value decomposition (SVD) to the analysis of IFIE data in order to extract the dynamically cooperative interactions between the ligand and the residues. We found that the relative importance of each residue is modified via the structural fluctuations and that the ligand is bound in the pharmacophore in a dynamic manner through collective interactions formed by multiple residues, thus providing new insight into structure-based drug discovery.
Collapse
Affiliation(s)
- Shigenori Tanaka
- Graduate School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Shusuke Tokutomi
- Graduate School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Ryo Hatada
- 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
| | - Koji Okuwaki
- 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
| | - Kazuki Akisawa
- 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
| | - 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 980-8579, Japan.,Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Yuto Komeiji
- Biomedical Research Institute, AIST, Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
| | - Yoshio Okiyama
- Division of Medicinal Safety Science, National Institute of Health Sciences, 3-25-26 Tonomachi, Kawasaki-ku, Kawasaki, Kanagawa 201-9501, Japan
| | - Yuji Mochizuki
- 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.,Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
6
|
Kato K, Honma T, Fukuzawa K. Intermolecular interaction among Remdesivir, RNA and RNA-dependent RNA polymerase of SARS-CoV-2 analyzed by fragment molecular orbital calculation. J Mol Graph Model 2020; 100:107695. [PMID: 32702590 PMCID: PMC7363421 DOI: 10.1016/j.jmgm.2020.107695] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Revised: 06/30/2020] [Accepted: 07/06/2020] [Indexed: 01/18/2023]
Abstract
COVID-19, a disease caused by a new strain of coronavirus (SARS-CoV-2) originating from Wuhan, China, has now spread around the world, triggering a global pandemic, leaving the public eagerly awaiting the development of a specific medicine and vaccine. In response, aggressive efforts are underway around the world to overcome COVID-19. In this study, referencing the data published on the Protein Data Bank (PDB ID: 7BV2) on April 22, we conducted a detailed analysis of the interaction between the complex structures of the RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 and Remdesivir, an antiviral drug, from the quantum chemical perspective based on the fragment molecular orbital (FMO) method. In addition to the hydrogen bonding and intra-strand stacking between complementary strands as seen in normal base pairs, Remdesivir bound to the terminus of an primer-RNA strand was further stabilized by diagonal π-π stacking with the -1A' base of the complementary strand and an additional hydrogen bond with an intra-strand base, due to the effect of chemically modified functional group. Moreover, stable OH/π interaction is also formed with Thr687 of the RdRp. We quantitatively revealed the exhaustive interaction within the complex among Remdesivir, template-primer-RNA, RdRp and co-factors, and published the results in the FMODB database.
Collapse
Affiliation(s)
- Koichiro Kato
- Department of Applied Chemistry, Graduate School of Engineering, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka, 819-0395, Japan; Center for Molecular Systems (CMS), Kyushu University, 744 Motooka, Nishi-ku, Fukuoka, 819-0395, Japan.
| | - Teruki Honma
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama, Kanagawa, 230-0045, Japan.
| | - Kaori Fukuzawa
- Department of Physical Chemistry, 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, 980-8579, Japan.
| |
Collapse
|
7
|
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
|
8
|
Rossberg A, Abe T, Okuwaki K, Barkleit A, Fukuzawa K, Nakano T, Mochizuki Y, Tsushima S. Destabilization of DNA through interstrand crosslinking by UO 22. Chem Commun (Camb) 2019; 55:2015-2018. [PMID: 30643910 DOI: 10.1039/c8cc09329f] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
UO22+ was shown to form an interstrand crosslink between two different strands of a single DNA molecule. This crosslink hardly affected the hydrogen bonds between nucleobase pairs but destabilized the π-π stacking between the two nucleobases in the vicinity of UO22+-bound phosphate. Thereby, the fragility of the DNA backbone increased upon UO22+ binding.
Collapse
Affiliation(s)
- André Rossberg
- Institute of Resource Ecology, Helmholtz-Zentrum Dresden-Rossendorf (HZDR), Dresden, 01328, Germany.
| | | | | | | | | | | | | | | |
Collapse
|
9
|
Drobot B, Schmidt M, Mochizuki Y, Abe T, Okuwaki K, Brulfert F, Falke S, Samsonov SA, Komeiji Y, Betzel C, Stumpf T, Raff J, Tsushima S. Cm3+/Eu3+induced structural, mechanistic and functional implications for calmodulin. Phys Chem Chem Phys 2019; 21:21213-21222. [DOI: 10.1039/c9cp03750k] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Trivalent lanthanide and actinide can strongly bind to calmodulin (CaM). The global structure of Ln/An-bound CaM were found to be similar to Ca-CaM but the local environment around Ln/An is distorted giving less structural rigidity to Ln/An-CaM.
Collapse
|
10
|
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
|
11
|
Komeiji Y, Okiyama Y, Mochizuki Y, Fukuzawa K. Interaction between a Single-Stranded DNA and a Binding Protein Viewed by the Fragment Molecular Orbital Method. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2018. [DOI: 10.1246/bcsj.20180150] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Yuto Komeiji
- Biomedical Research Institute, National Institute of Industrial Science and Technology, AIST Tsukuba Central 6, Tsukuba, Ibaraki 305-8566, Japan
| | - Yoshio Okiyama
- Division of Medicinal Safety Science, National Institute of Health Sciences, 3-25-26 Tonomachi, Kawasaki-ku, Kanagawa 210-9501, Japan
| | - Yuji Mochizuki
- 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
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Kaori Fukuzawa
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa-ku, Tokyo 142-8501, Japan
| |
Collapse
|
12
|
Okiyama Y, Nakano T, Watanabe C, Fukuzawa K, Mochizuki Y, Tanaka S. Fragment Molecular Orbital Calculations with Implicit Solvent Based on the Poisson–Boltzmann Equation: Implementation and DNA Study. J Phys Chem B 2018; 122:4457-4471. [DOI: 10.1021/acs.jpcb.8b01172] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Yoshio Okiyama
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, 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
| | - Chiduru Watanabe
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, 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
| | - Shigenori Tanaka
- Graduate School of System Informatics, Kobe University, 1-1 Rokkodai, Nada-ku, Kobe, Hyogo 657-8501, Japan
| |
Collapse
|
13
|
Watanabe C, Watanabe H, Fukuzawa K, Parker LJ, Okiyama Y, Yuki H, Yokoyama S, Nakano H, Tanaka S, Honma T. Theoretical Analysis of Activity Cliffs among Benzofuranone-Class Pim1 Inhibitors Using the Fragment Molecular Orbital Method with Molecular Mechanics Poisson-Boltzmann Surface Area (FMO+MM-PBSA) Approach. J Chem Inf Model 2017; 57:2996-3010. [PMID: 29111719 DOI: 10.1021/acs.jcim.7b00110] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Significant activity changes due to small structural changes (i.e., activity cliffs) of serine/threonine kinase Pim1 inhibitors were studied theoretically using the fragment molecular orbital method with molecular mechanics Poisson-Boltzmann surface area (FMO+MM-PBSA) approach. This methodology enables quantum-chemical calculations for large biomolecules with solvation. In the course of drug discovery targeting Pim1, six benzofuranone-class inhibitors were found to differ only in the position of the indole-ring nitrogen atom. By comparing the various qualities of complex structures based on X-ray, classical molecular mechanics (MM)-optimized, and quantum/molecular mechanics (QM/MM)-optimized structures, we found that the QM/MM-optimized structures provided the best correlation (R2 = 0.85) between pIC50 and the calculated FMO+MM-PBSA binding energy. Combining the classical solvation energy with the QM binding energy was important to increase the correlation. In addition, decomposition of the interaction energy into various physicochemical components by pair interaction energy decomposition analysis suggested that CH-π and electrostatic interactions mainly caused the activity differences.
Collapse
Affiliation(s)
- Chiduru Watanabe
- RIKEN Center for Life Science Technologies , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan.,Institute of Industrial Science, The University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hirofumi Watanabe
- RIKEN Center for Life Science Technologies , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Kaori Fukuzawa
- Institute of Industrial Science, The University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan.,Department of Physical Chemistry, School of Pharmacy and Pharmaceutical Sciences, Hoshi University , 2-4-41 Ebara, Shinagawa, Tokyo 142-8501, Japan
| | - Lorien J Parker
- RIKEN Structural Biology Laboratory , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan.,Department of Structural Biology, St. Vincent's Institute , 9 Princes Street, Fitzroy, Victoria 3065, Australia
| | - Yoshio Okiyama
- RIKEN Center for Life Science Technologies , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Hitomi Yuki
- RIKEN Center for Life Science Technologies , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Shigeyuki Yokoyama
- RIKEN Structural Biology Laboratory , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Hirofumi Nakano
- Drug Discovery Initiative, The University of Tokyo , 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Shigenori Tanaka
- Graduate School of System Informatics, Kobe University , 1-1 Rokkodai, Nada-ku, Kobe 657-8501, Japan
| | - Teruki Honma
- RIKEN Center for Life Science Technologies , 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| |
Collapse
|
14
|
Yamada H, Mochizuki Y, Fukuzawa K, Okiyama Y, Komeiji Y. Fragment molecular orbital (FMO) calculations on DNA by a scaled third-order Møller-Plesset perturbation (MP2.5) scheme. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2016.12.008] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
15
|
Komeij Y, Okiyama Y, Mochizuki Y, Fukuzawa K. Explicit solvation of a single-stranded DNA, a binding protein, and their complex: a suitable protocol for fragment molecular orbital calculation. CHEM-BIO INFORMATICS JOURNAL 2017. [DOI: 10.1273/cbij.17.72] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Affiliation(s)
- Yuto Komeij
- Biomedical Research Institute, National Institute of Advanced Industrial Science and Technology
| | - Yoshio Okiyama
- Division of Medicinal Safety Science, National Institute of Health Sciences
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University
- Institute of Industrial Science, The University of Tokyo
| | - Kaori Fukuzawa
- Institute of Industrial Science, The University of Tokyo
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University
| |
Collapse
|
16
|
Hydration of ligands of influenza virus neuraminidase studied by the fragment molecular orbital method. J Mol Graph Model 2016; 69:144-53. [DOI: 10.1016/j.jmgm.2016.08.004] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2016] [Revised: 07/29/2016] [Accepted: 08/10/2016] [Indexed: 11/18/2022]
|