1
|
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
|
2
|
Thapa B, Beckett D, Erickson J, Raghavachari K. Theoretical Study of Protein–Ligand Interactions Using the Molecules-in-Molecules Fragmentation-Based Method. J Chem Theory Comput 2018; 14:5143-5155. [DOI: 10.1021/acs.jctc.8b00531] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Jon Erickson
- Lilly Research Laboratories, Eli Lilly & Co., Indianapolis, Indiana 47285, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
3
|
Nedjoua D, Krallafa AM. Temperature effect on the structure and conformational fluctuations in two zinc knuckles from the mouse mammary tumor virus. Comput Biol Chem 2018; 74:86-93. [PMID: 29567490 DOI: 10.1016/j.compbiolchem.2018.03.005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2017] [Revised: 10/28/2017] [Accepted: 03/07/2018] [Indexed: 11/18/2022]
Abstract
Zinc fingers are small protein domains in which zinc plays a structural role, contributing to the stability of the zinc-peptide complex. Zinc fingers are structurally diverse and are present in proteins that perform a broad range of functions in various cellular processes, such as replication and repair, transcription and translation, metabolism and signaling, cell proliferation, and apoptosis. Zinc fingers typically function as interaction modules and bind to a wide variety of compounds, such as nucleic acids, proteins, and small molecules. In this study, we investigated the structural properties, in solution, of the proximal and distal zinc knuckles of the nucleocapsid (NC) protein from the mouse mammary tumor virus (MMTV) (MMTV NC). For this purpose, we performed a series of molecular dynamics simulations in aqueous solution at 300 K, 333 K, and 348 K. The temperature effect was evaluated in terms of root mean square deviation of the backbone atoms and root mean square fluctuation of the coordinating residue atoms. The stability of the zinc coordination sphere was analyzed based upon the time profile of the interatomic distances between the zinc ions and the chelator atoms. The results indicate that the hydrophobic character of the proximal zinc finger is dominant at 333 K. The low mobility of the coordinating residues suggests that the strong electrostatic effect exerted by the zinc ion on its coordinating residues is not influenced by the increase in temperature. The evolution of the structural parameters of the coordination sphere of the distal zinc finger at 300 K gives us a reasonable picture of the unfolding pathway, as proposed by Bombarda and coworkers (Bombarda et al., 2005), which can predict the binding order of the four conserved ligand-binding residues. Our results support the conclusion that the structural features can vary significantly between the two zinc knuckles of MMTV NC.
Collapse
Affiliation(s)
- Drici Nedjoua
- LCPM, Department of Chemistry, University of Oran 1 Ahmed Benbella, PO Box 1524, El m'naouer, Oran, 31000, Algeria.
| | - Abdelghani Mohamed Krallafa
- LCPM, Department of Chemistry, University of Oran 1 Ahmed Benbella, PO Box 1524, El m'naouer, Oran, 31000, Algeria.
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
Ueno-Noto K, Takano K. Water molecules inside protein structure affect binding of monosaccharides with HIV-1 antibody 2G12. J Comput Chem 2016; 37:2341-8. [DOI: 10.1002/jcc.24447] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2016] [Revised: 06/03/2016] [Accepted: 06/16/2016] [Indexed: 11/07/2022]
Affiliation(s)
- Kaori Ueno-Noto
- Chemistry Section, 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
|
7
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
8
|
Phipps MJS, Fox T, Tautermann CS, Skylaris CK. Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns. Chem Soc Rev 2015; 44:3177-211. [DOI: 10.1039/c4cs00375f] [Citation(s) in RCA: 226] [Impact Index Per Article: 25.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.
Collapse
Affiliation(s)
| | - Thomas Fox
- Lead Identification and Optimization Support
- Boehringer Ingelheim Pharma GmbH & Co. KG
- 88397 Biberach
- Germany
| | - Christofer S. Tautermann
- Lead Identification and Optimization Support
- Boehringer Ingelheim Pharma GmbH & Co. KG
- 88397 Biberach
- Germany
| | | |
Collapse
|