1
|
Pan X, Van R, Pu J, Nam K, Mao Y, Shao Y. Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions. J Chem Theory Comput 2023; 19:8234-8244. [PMID: 37943896 PMCID: PMC10835707 DOI: 10.1021/acs.jctc.3c00973] [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] [Indexed: 11/12/2023]
Abstract
In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.
Collapse
Affiliation(s)
- Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Richard Van
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
- Laboratory of Computational Biology, National, Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| |
Collapse
|
2
|
Georgiev DD. Quantum information theoretic approach to the mind–brain problem. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2020; 158:16-32. [DOI: 10.1016/j.pbiomolbio.2020.08.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Revised: 08/02/2020] [Accepted: 08/05/2020] [Indexed: 12/25/2022]
|
3
|
Ishida T. Computational analysis of carbohydrate recognition based on hybrid QM/MM modeling: a case study of norovirus capsid protein in complex with Lewis antigen. Phys Chem Chem Phys 2018; 20:4652-4665. [PMID: 29372731 DOI: 10.1039/c7cp07701g] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Norovirus is a major pathogen of nonbacterial acute gastroenteritis in humans and animals. Carbohydrate recognition between norovirus capsid proteins and Lewis antigens is considered to play a critical role in initiating infection of eukaryotic cells. In this article, we first report a detailed atomistic simulation study of the norovirus capsid protein in complex with the Lewis antigen based on ab initio QM/MM combined with MD-FEP simulations. To understand the mechanistic details of ligand binding, we analyzed and compared the carbohydrate recognition mechanism of the wild-type P domain protein with a mutant protein. Small structural differences between two capsid proteins are observed on the weak interaction site of residue 389, which is located on the solvent exposed surface of the P domain. To further clarify affinity differences in ligand binding, we directly evaluated free energy changes of the ligand binding process. Although the mutant protein loses its interaction energy with the Lewis antigen, this small amount of energy penalty is compensated for by an increase in the solvation stability, which is induced by structural reorganization at the ligand binding site on the protein surface. As a sum of these opposite energy components, the mutant P domain obtains a slightly enhanced binding affinity for the Lewis antigen. The present computational study clearly demonstrated that a detailed free energy balance of the interaction energy between the capsid protein and the surrounding aqueous solvent is the mechanistic basis of carbohydrate recognition in the norovirus capsid protein.
Collapse
Affiliation(s)
- Toyokazu Ishida
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba Central 2, 1-1-1 Umezono, Tsukuba, 305-8568, Japan.
| |
Collapse
|
4
|
Freindorf M, Tao Y, Sethio D, Cremer D, Kraka E. New mechanistic insights into the Claisen rearrangement of chorismate – a Unified Reaction Valley Approach study. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1530464] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Marek Freindorf
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Yunwen Tao
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Daniel Sethio
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Dieter Cremer
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| | - Elfi Kraka
- Computational and Theoretical Chemistry Group (CATCO), Department of Chemistry, Southern Methodist University, Dallas, TX, USA
| |
Collapse
|
5
|
Thapa B, Beckett D, Jovan Jose KV, Raghavachari K. Assessment of Fragmentation Strategies for Large Proteins Using the Multilayer Molecules-in-Molecules Approach. J Chem Theory Comput 2018; 14:1383-1394. [DOI: 10.1021/acs.jctc.7b01198] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| |
Collapse
|
6
|
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
|
7
|
Pruitt SR, Steinmann C. Mapping Interaction Energies in Chorismate Mutase with the Fragment Molecular Orbital Method. J Phys Chem A 2017; 121:1797-1807. [DOI: 10.1021/acs.jpca.6b12830] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Spencer R. Pruitt
- Academic & Research Computing, Worcester Polytechnic Institute, Worcester, Massachusetts 01602, United States
| | - Casper Steinmann
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
8
|
Saha A, Raghavachari K. Analysis of Different Fragmentation Strategies on a Variety of Large Peptides: Implementation of a Low Level of Theory in Fragment-Based Methods Can Be a Crucial Factor. J Chem Theory Comput 2016; 11:2012-23. [PMID: 26574406 DOI: 10.1021/ct501045s] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We have investigated the performance of two classes of fragmentation methods developed in our group (Molecules-in-Molecules (MIM) and Many-Overlapping-Body (MOB) expansion), to reproduce the unfragmented MP2 energies on a test set composed of 10 small to large biomolecules. They have also been assessed to recover the relative energies of different motifs of the acetyl(ala)18NH2 system. Performance of different bond-cutting environments and the use of Hartree-Fock and different density functionals (as a low level of theory) in conjunction with the fragmentation strategies have been analyzed. Our investigation shows that while a low level of theory (for recovering long-range interactions) may not be necessary for small peptides, it provides a very effective strategy to accurately reproduce the total and relative energies of larger peptides such as the different motifs of the acetyl(ala)18NH2 system. Employing M06-2X as the low level of theory, the calculated mean total energy deviation (maximum deviation) in the total MP2 energies for the 10 molecules in the test set at MIM(d=3.5Å), MIM(η=9), and MOB(d=5Å) are 1.16 (2.31), 0.72 (1.87), and 0.43 (2.02) kcal/mol, respectively. The excellent performance suggests that such fragment-based methods should be of general use for the computation of accurate energies of large biomolecular systems.
Collapse
Affiliation(s)
- Arjun Saha
- Department of Chemistry, Indiana University , Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University , Bloomington, Indiana 47405, United States
| |
Collapse
|
9
|
Vasilevskaya T, Thiel W. Periodic Boundary Conditions in QM/MM Calculations: Implementation and Tests. J Chem Theory Comput 2016; 12:3561-70. [DOI: 10.1021/acs.jctc.6b00269] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Walter Thiel
- Max-Planck-Institut
für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
10
|
Nakata H, Fedorov DG, Nagata T, Kitaura K, Nakamura S. Simulations of Chemical Reactions with the Frozen Domain Formulation of the Fragment Molecular Orbital Method. J Chem Theory Comput 2015; 11:3053-64. [DOI: 10.1021/acs.jctc.5b00277] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hiroya Nakata
- Department
of Biomolecular Engineering, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8501, Japan
- Research
Cluster for Innovation, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Japan Society for the Promotion of Science, Kojimachi Business Center Building, 5-3-1 Kojimachi, Chiyoda-ku, Tokyo 102-0083, Japan
| | - Dmitri G. Fedorov
- Nanosystem
Research Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Takeshi Nagata
- Nanosystem
Research Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| | - Kazuo Kitaura
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| | - Shinichiro Nakamura
- Research
Cluster for Innovation, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
11
|
Christensen AS, Steinmann C, Fedorov DG, Jensen JH. Hybrid RHF/MP2 geometry optimizations with the effective fragment molecular orbital method. PLoS One 2014; 9:e88800. [PMID: 24558430 PMCID: PMC3928295 DOI: 10.1371/journal.pone.0088800] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2013] [Accepted: 01/15/2014] [Indexed: 11/18/2022] Open
Abstract
The frozen domain effective fragment molecular orbital method is extended to allow for the treatment of a single fragment at the MP2 level of theory. The approach is applied to the conversion of chorismate to prephenate by Chorismate Mutase, where the substrate is treated at the MP2 level of theory while the rest of the system is treated at the RHF level. MP2 geometry optimization is found to lower the barrier by up to 3.5 kcal/mol compared to RHF optimzations and ONIOM energy refinement and leads to a smoother convergence with respect to the basis set for the reaction profile. For double zeta basis sets the increase in CPU time relative to RHF is roughly a factor of two.
Collapse
Affiliation(s)
| | - Casper Steinmann
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense, Denmark
| | - Dmitri G. Fedorov
- Nanosystem Research Institute (NRI), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
| | - Jan H. Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
- * E-mail:
| |
Collapse
|
12
|
Choutko A, Eichenberger AP, van Gunsteren WF, Dolenc J. Exploration of swapping enzymatic function between two proteins: a simulation study of chorismate mutase and isochorismate pyruvate lyase. Protein Sci 2013; 22:809-22. [PMID: 23595942 DOI: 10.1002/pro.2264] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2013] [Accepted: 03/31/2013] [Indexed: 01/09/2023]
Abstract
The enzyme chorismate mutase EcCM from Escherichia coli catalyzes one of the few pericyclic reactions in biology, the transformation of chorismate to prephenate. The isochorismate pyruvate lyase PchB from Pseudomonas aeroginosa catalyzes another pericyclic reaction, the isochorismate to salicylate transformation. Interestingly, PchB possesses weak chorismate mutase activity as well thus being able to catalyze two distinct pericyclic reactions in a single active site. EcCM and PchB possess very similar folds, despite their low sequence identity. Using molecular dynamics simulations of four combinations of the two enzymes (EcCM and PchB) with the two substrates (chorismate and isochorismate) we show that the electrostatic field due to EcCM at atoms of chorismate favors the chorismate to prephenate transition and that, analogously, the electrostatic field due to PchB at atoms of isochorismate favors the isochorismate to salicylate transition. The largest differences between EcCM and PchB in electrostatic field strengths at atoms of the substrates are found to be due to residue side chains at distances between 0.6 and 0.8 nm from particular substrate atoms. Both enzymes tend to bring their non-native substrate in the same conformation as their native substrate. EcCM and to a lower extent PchB fail in influencing the forces on and conformations of the substrate such as to favor the other chemical reaction (isochorismate pyruvate lyase activity for EcCM and chorismate mutase activity for PchB). These observations might explain the difficulty of engineering isochorismate pyruvate lyase activity in EcCM by solely mutating active site residues.
Collapse
Affiliation(s)
- Alexandra Choutko
- Physical Chemistry, Swiss Federal Institute of Technology, ETH, Zürich, Switzerland
| | | | | | | |
Collapse
|
13
|
Nakata H, Nagata T, Fedorov DG, Yokojima S, Kitaura K, Nakamura S. Analytic second derivatives of the energy in the fragment molecular orbital method. J Chem Phys 2013; 138:164103. [DOI: 10.1063/1.4800990] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Steinmann C, Fedorov DG, Jensen JH. Mapping enzymatic catalysis using the effective fragment molecular orbital method: towards all ab initio biochemistry. PLoS One 2013; 8:e60602. [PMID: 23593259 PMCID: PMC3625203 DOI: 10.1371/journal.pone.0060602] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2012] [Accepted: 02/28/2013] [Indexed: 12/02/2022] Open
Abstract
We extend the Effective Fragment Molecular Orbital (EFMO) method to the frozen domain approach where only the geometry of an active part is optimized, while the many-body polarization effects are considered for the whole system. The new approach efficiently mapped out the entire reaction path of chorismate mutase in less than four days using 80 cores on 20 nodes, where the whole system containing 2398 atoms is treated in the ab initio fashion without using any force fields. The reaction path is constructed automatically with the only assumption of defining the reaction coordinate a priori. We determine the reaction barrier of chorismate mutase to be kcal mol−1 for MP2/cc-pVDZ and for MP2/cc-pVTZ in an ONIOM approach using EFMO-RHF/6-31G(d) for the high and low layers, respectively.
Collapse
Affiliation(s)
- Casper Steinmann
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmitri G. Fedorov
- Nanosystems Research Institute (NRI), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
| | - Jan H. Jensen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
- * E-mail:
| |
Collapse
|
15
|
Fedorov DG, Nagata T, Kitaura K. Exploring chemistry with the fragment molecular orbital method. Phys Chem Chem Phys 2012; 14:7562-77. [DOI: 10.1039/c2cp23784a] [Citation(s) in RCA: 290] [Impact Index Per Article: 24.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
16
|
Gordon MS, Fedorov DG, Pruitt SR, Slipchenko LV. Fragmentation Methods: A Route to Accurate Calculations on Large Systems. Chem Rev 2011; 112:632-72. [DOI: 10.1021/cr200093j] [Citation(s) in RCA: 836] [Impact Index Per Article: 64.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Mark S. Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Dmitri G. Fedorov
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Spencer R. Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Lyudmila V. Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
17
|
Brorsen K, Fedorov DG. Fully analytic energy gradient in the fragment molecular orbital method. J Chem Phys 2011; 134:124115. [DOI: 10.1063/1.3568010] [Citation(s) in RCA: 88] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
Claeyssens F, Ranaghan KE, Lawan N, Macrae SJ, Manby FR, Harvey JN, Mulholland AJ. Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions. Org Biomol Chem 2011; 9:1578-90. [PMID: 21243152 DOI: 10.1039/c0ob00691b] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024]
Abstract
Chorismate mutase is at the centre of current controversy about fundamental features of biological catalysts. Some recent studies have proposed that catalysis in this enzyme does not involve transition state (TS) stabilization but instead is due largely to the formation of a reactive conformation of the substrate. To understand the origins of catalysis, it is necessary to compare equivalent reactions in different environments. The pericyclic conversion of chorismate to prephenate catalysed by chorismate mutase also occurs (much more slowly) in aqueous solution. In this study we analyse the origins of catalysis by comparison of multiple quantum mechanics/molecular mechanics (QM/MM) reaction pathways at a reliable, well tested level of theory (B3LYP/6-31G(d)/CHARMM27) for the reaction (i) in Bacillus subtilis chorismate mutase (BsCM) and (ii) in aqueous solvent. The average calculated reaction (potential energy) barriers are 11.3 kcal mol(-1) in the enzyme and 17.4 kcal mol(-1) in water, both of which are in good agreement with experiment. Comparison of the two sets of reaction pathways shows that the reaction follows a slightly different reaction pathway in the enzyme than in it does in solution, because of a destabilization, or strain, of the substrate in the enzyme. The substrate strain energy within the enzyme remains constant throughout the reaction. There is no unique reactive conformation of the substrate common to both environments, and the transition state structures are also different in the enzyme and in water. Analysis of the barrier heights in each environment shows a clear correlation between TS stabilization and the barrier height. The average differential TS stabilization is 7.3 kcal mol(-1) in the enzyme. This is significantly higher than the small amount of TS stabilization in water (on average only 1.0 kcal mol(-1) relative to the substrate). The TS is stabilized mainly by electrostatic interactions with active site residues in the enzyme, with Arg90, Arg7 and Glu78 generally the most important. Conformational effects (e.g. strain of the substrate in the enzyme) do not contribute significantly to the lower barrier observed in the enzyme. The results show that catalysis is mainly due to better TS stabilization by the enzyme.
Collapse
Affiliation(s)
- Frederik Claeyssens
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK BS8 1TS
| | | | | | | | | | | | | |
Collapse
|
19
|
McGeagh JD, Ranaghan KE, Mulholland AJ. Protein dynamics and enzyme catalysis: insights from simulations. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2010; 1814:1077-92. [PMID: 21167324 DOI: 10.1016/j.bbapap.2010.12.002] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/20/2010] [Revised: 11/25/2010] [Accepted: 12/03/2010] [Indexed: 10/18/2022]
Abstract
The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.
Collapse
Affiliation(s)
- John D McGeagh
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Cantock's Close, BS8 1TS, United Kingdom
| | | | | |
Collapse
|
20
|
Ishida T. Effects of Point Mutation on Enzymatic Activity: Correlation between Protein Electronic Structure and Motion in Chorismate Mutase Reaction. J Am Chem Soc 2010; 132:7104-18. [DOI: 10.1021/ja100744h] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Toyokazu Ishida
- Research Institute for Computational Sciences (RICS), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba Central 2, 1-1-1 Umezono, Tsukuba 305-8568, Japan
| |
Collapse
|
21
|
Gao Q, Yokojima S, Fedorov DG, Kitaura K, Sakurai M, Nakamura S. Fragment-Molecular-Orbital-Method-Based ab Initio NMR Chemical-Shift Calculations for Large Molecular Systems. J Chem Theory Comput 2010. [DOI: 10.1021/ct100006n] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Qi Gao
- Mitsubishi Chemical Group Science and Technology Research Center, Inc., 1000 Kamochida-cho, Aoba-ku, Yokohama 227-8502, Japan, Center for Biological Resources and Informatics, Tokyo Institute of Technology, 4259 Nagatsuda-cho, Midori-ku, Yokohama 226-8501, Japan, The KAITEKI Institute, Inc. 14-1, Shiba 4-chome, Minato-ku, Tokyo 108-0014, Japan, RICS, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, and Graduate School of
| | - Satoshi Yokojima
- Mitsubishi Chemical Group Science and Technology Research Center, Inc., 1000 Kamochida-cho, Aoba-ku, Yokohama 227-8502, Japan, Center for Biological Resources and Informatics, Tokyo Institute of Technology, 4259 Nagatsuda-cho, Midori-ku, Yokohama 226-8501, Japan, The KAITEKI Institute, Inc. 14-1, Shiba 4-chome, Minato-ku, Tokyo 108-0014, Japan, RICS, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, and Graduate School of
| | - Dmitri G. Fedorov
- Mitsubishi Chemical Group Science and Technology Research Center, Inc., 1000 Kamochida-cho, Aoba-ku, Yokohama 227-8502, Japan, Center for Biological Resources and Informatics, Tokyo Institute of Technology, 4259 Nagatsuda-cho, Midori-ku, Yokohama 226-8501, Japan, The KAITEKI Institute, Inc. 14-1, Shiba 4-chome, Minato-ku, Tokyo 108-0014, Japan, RICS, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, and Graduate School of
| | - Kazuo Kitaura
- Mitsubishi Chemical Group Science and Technology Research Center, Inc., 1000 Kamochida-cho, Aoba-ku, Yokohama 227-8502, Japan, Center for Biological Resources and Informatics, Tokyo Institute of Technology, 4259 Nagatsuda-cho, Midori-ku, Yokohama 226-8501, Japan, The KAITEKI Institute, Inc. 14-1, Shiba 4-chome, Minato-ku, Tokyo 108-0014, Japan, RICS, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, and Graduate School of
| | - Minoru Sakurai
- Mitsubishi Chemical Group Science and Technology Research Center, Inc., 1000 Kamochida-cho, Aoba-ku, Yokohama 227-8502, Japan, Center for Biological Resources and Informatics, Tokyo Institute of Technology, 4259 Nagatsuda-cho, Midori-ku, Yokohama 226-8501, Japan, The KAITEKI Institute, Inc. 14-1, Shiba 4-chome, Minato-ku, Tokyo 108-0014, Japan, RICS, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, and Graduate School of
| | - Shinichiro Nakamura
- Mitsubishi Chemical Group Science and Technology Research Center, Inc., 1000 Kamochida-cho, Aoba-ku, Yokohama 227-8502, Japan, Center for Biological Resources and Informatics, Tokyo Institute of Technology, 4259 Nagatsuda-cho, Midori-ku, Yokohama 226-8501, Japan, The KAITEKI Institute, Inc. 14-1, Shiba 4-chome, Minato-ku, Tokyo 108-0014, Japan, RICS, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, and Graduate School of
| |
Collapse
|
22
|
Abstract
Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.
Collapse
Affiliation(s)
- Hans Martin Senn
- Department of Chemistry, WestCHEM and University of Glasgow, Glasgow G12 8QQ, UK.
| | | |
Collapse
|
23
|
|
24
|
Li H, Fedorov DG, Nagata T, Kitaura K, Jensen JH, Gordon MS. Energy gradients in combined fragment molecular orbital and polarizable continuum model (FMO/PCM) calculation. J Comput Chem 2009; 31:778-90. [DOI: 10.1002/jcc.21363] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
25
|
Fedorov DG, Jensen JH, Deka RC, Kitaura K. Covalent Bond Fragmentation Suitable To Describe Solids in the Fragment Molecular Orbital Method. J Phys Chem A 2008; 112:11808-16. [DOI: 10.1021/jp805435n] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Institute for Computational Sciences, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark, Department of Chemical Sciences, Tezpur University, Napaam, Tezpur 784028, Assam, India, and Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Jan H. Jensen
- Research Institute for Computational Sciences, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark, Department of Chemical Sciences, Tezpur University, Napaam, Tezpur 784028, Assam, India, and Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Ramesh C. Deka
- Research Institute for Computational Sciences, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark, Department of Chemical Sciences, Tezpur University, Napaam, Tezpur 784028, Assam, India, and Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Kazuo Kitaura
- Research Institute for Computational Sciences, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark, Department of Chemical Sciences, Tezpur University, Napaam, Tezpur 784028, Assam, India, and Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| |
Collapse
|
26
|
Ishida T. Probing protein environment in an enzymatic process: All-electron quantum chemical analysis combined with ab initio quantum mechanical/molecular mechanical modeling of chorismate mutase. J Chem Phys 2008; 129:125105. [DOI: 10.1063/1.2977458] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
27
|
Jiménez A, Clapés P, Crehuet R. A dynamic view of enzyme catalysis. J Mol Model 2008; 14:735-46. [DOI: 10.1007/s00894-008-0283-2] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2007] [Accepted: 02/01/2008] [Indexed: 10/22/2022]
|
28
|
Fedorov DG, Kitaura K. Extending the Power of Quantum Chemistry to Large Systems with the Fragment Molecular Orbital Method. J Phys Chem A 2007; 111:6904-14. [PMID: 17511437 DOI: 10.1021/jp0716740] [Citation(s) in RCA: 437] [Impact Index Per Article: 25.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Following the brief review of the modern fragment-based methods and other approaches to perform quantum-mechanical calculations of large systems, the theoretical development of the fragment molecular orbital method (FMO) is covered in detail, with the emphasis on the physical properties, which can be computed with FMO. The FMO-based polarizable continuum model (PCM) for treating the solvent effects in large systems and the pair interaction energy decomposition analysis (PIEDA) are described in some detail, and a range of applications of FMO to biological studies is introduced. The factors determining the relative stability of polypeptide conformers (alpha-helix, beta-turn, and extended form) are elucidated using FMO/PCM and PIEDA, and the interactions in the Trp-cage miniprotein construct (PDB: 1L2Y) are analyzed using PIEDA.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Institute for Computational Sciences (RICS), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki, Japan 305-8568.
| | | |
Collapse
|
29
|
Time-dependent density functional theory with the multilayer fragment molecular orbital method. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.07.034] [Citation(s) in RCA: 61] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
30
|
Watanabe T, Inadomi Y, Fukuzawa K, Nakano T, Tanaka S, Nilsson L, Nagashima U. DNA and Estrogen Receptor Interaction Revealed by Fragment Molecular Orbital Calculations. J Phys Chem B 2007; 111:9621-7. [PMID: 17649990 DOI: 10.1021/jp071710v] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular orbital calculations of the complex between DNA-ERE (estrogen response element) and ER (estrogen receptor)-DBD (DNA-binding domain) were performed using the fragment molecular orbital (FMO) method, which enables large-scale MO (molecular orbital) calculations by reducing the computational cost and by significantly increasing efficiency for parallel computation. Such a large system, which contains 3354 atoms, is impractical via conventional MO methods due to the immense computational cost. Details of the interaction between DNA-ERE and ER-DBD were revealed in this study as follows by using the FMO calculations to analyze the interfragment interaction energies (IFIEs) and the electrostatic potentials (ESPs). An area with a high positive ESP is identified on the DNA-binding side of ER-DBD and is the main driving force behind access to the DNA. The position of the ER-DBD monomer can be fixed on a phosphate group of DNA-ERE by the strong electrostatic interactions, whereas the rotation cannot be fixed. In contrast, both the position and rotation of the ER-DBD dimer can be fixed and can therefore form the stable (ER-DBD)2...DNA-ERE complex. Dimerization of the ER-DBD monomers, each of which have a charge of +5 , is mainly due to large attractive interaction energies of the second Zn fragments. The base pairs in the consensus sequence of DNA-ERE interact only with the recognition helix located in the major groove due to the large shielding effect of the phosphate groups of DNA. The recognition helix has weaker interactions with the base pairs than the electrostatic interactions with the phosphate groups. Thus, the DNA-binding machinery of the ER-DBD dimer, which can secure the recognition helix in the major groove of DNA, is crucial for interactions between the recognition helix and base pairs.
Collapse
Affiliation(s)
- Toshio Watanabe
- Research Institute for Computational Science, National Institute of Advanced Industrial Science and Technology, Umezono 1-1-1, Tsukuba, Ibaraki 305-8568, Japan
| | | | | | | | | | | | | |
Collapse
|
31
|
Fedorov DG, Ishimura K, Ishida T, Kitaura K, Pulay P, Nagase S. Accuracy of the three-body fragment molecular orbital method applied to Møller-Plesset perturbation theory. J Comput Chem 2007; 28:1476-1484. [PMID: 17330884 DOI: 10.1002/jcc.20645] [Citation(s) in RCA: 76] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The three-body energy expansion in the fragment molecular orbital method (FMO) was applied to the 2nd order Møller-Plesset theory (MP2). The accuracy of both the two and three-body expansions was determined for water clusters, alanine n-mers (alpha-helices and beta-strands) and one synthetic protein, using the 6-31G* and 6-311G* basis sets. At the best level of theory (three-body, two molecules/residues per fragment), the absolute errors in energy relative to ab initio MP2 were at most 1.2 and 5.0 mhartree, for the 6-31G* and 6-311G* basis sets, respectively. The relative accuracy was at worst 99.996% and 99.96%, for 6-31G* and 6-311G*, respectively. A three-body approximation was introduced and the optimum threshold value was determined. The protein calculation (6-31G*) at the production level (FMO2/2) took 3 h on 36 3.2-GHz Pentium 4 nodes and had the absolute error in the MP2 correlation energy of only 2 kcal/mol.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Kazuya Ishimura
- Department of Theoretical Molecular Science, Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Toyokazu Ishida
- National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Kazuo Kitaura
- National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Peter Pulay
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701
| | - Shigeru Nagase
- Department of Theoretical Molecular Science, Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
32
|
Lassila JK, Keeffe JR, Kast P, Mayo SL. Exhaustive Mutagenesis of Six Secondary Active-Site Residues in Escherichia coli Chorismate Mutase Shows the Importance of Hydrophobic Side Chains and a Helix N-Capping Position for Stability and Catalysis. Biochemistry 2007; 46:6883-91. [PMID: 17506527 DOI: 10.1021/bi700215x] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Secondary active-site residues in enzymes, including hydrophobic amino acids, may contribute to catalysis through critical interactions that position the reacting molecule, organize hydrogen-bonding residues, and define the electrostatic environment of the active site. To ascertain the tolerance of an important model enzyme to mutation of active-site residues that do not directly hydrogen bond with the reacting molecule, all 19 possible amino acid substitutions were investigated in six positions of the engineered chorismate mutase domain of the Escherichia coli chorismate mutase-prephenate dehydratase. The six secondary active-site residues were selected to clarify results of a previous test of computational enzyme design procedures. Five of the positions encode hydrophobic side chains in the wild-type enzyme, and one forms a helix N-capping interaction as well as a salt bridge with a catalytically essential residue. Each mutant was evaluated for its ability to complement an auxotrophic chorismate mutase deletion strain. Kinetic parameters and thermal stabilities were measured for variants with in vivo activity. Altogether, we find that the enzyme tolerated 34% of the 114 possible substitutions, with a few mutations leading to increases in the catalytic efficiency of the enzyme. The results show the importance of secondary amino acid residues in determining enzymatic activity, and they point to strengths and weaknesses in current computational enzyme design procedures.
Collapse
Affiliation(s)
- Jonathan Kyle Lassila
- Biochemistry Option, California Institute of Technology, Pasadena, California 91125, USA
| | | | | | | |
Collapse
|
33
|
Higashi M, Hayashi S, Kato S. Transition state determination of enzyme reaction on free energy surface: Application to chorismate mutase. Chem Phys Lett 2007. [DOI: 10.1016/j.cplett.2007.02.044] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
34
|
Nakanishi I, Fedorov DG, Kitaura K. Molecular recognition mechanism of FK506 binding protein: An all-electron fragment molecular orbital study. Proteins 2007; 68:145-58. [PMID: 17387719 DOI: 10.1002/prot.21389] [Citation(s) in RCA: 60] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Abstract
The fragment molecular orbital (FMO) method has enabled electronic structure calculations and geometry optimizations of very large molecules with ab initio quality. We applied the method to four FK506 binding protein (FKBP) complexes (denoted by their PDB codes 1fkb, 1fkf, 1fkg, and 1fki) containing rapamycin, FK506, and two synthetic ligands. The geometries of reduced complex models were optimized at the restricted Hartree-Fock (FMO-RHF) level using the 3-21G basis set, and then for a better estimate of binding, the energetics were refined at a higher level of theory (2nd order Møller-Plesset perturbation theory FMO-MP2 with the 6-31G* basis set). Thus, obtained binding energies were -103.9 (-82.0), -102.2 (-69.2), -70.1 (-57.7), and -71.3 (-55.3) kcal/mol for 1fkb, 1fkf, 1fkg, and 1fki, respectively, where the correlation contribution is given in parentheses. The results show that the electron correlation contribution to binding is extremely important, and it accounts for 70-80% of the binding energy. The molecular recognition mechanism of FKBP was analyzed in detail based on the FMO-pair interactions between protein residues and the ligands. Solvation effects on the protein-ligand binding were estimated using the Poisson-Boltzmann/surface area model.
Collapse
Affiliation(s)
- Isao Nakanishi
- Department of Theoretical Drug Design, Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto, Japan.
| | | | | |
Collapse
|
35
|
Fedorov DG, Ishida T, Uebayasi M, Kitaura K. The Fragment Molecular Orbital Method for Geometry Optimizations of Polypeptides and Proteins. J Phys Chem A 2007; 111:2722-32. [PMID: 17388363 DOI: 10.1021/jp0671042] [Citation(s) in RCA: 111] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The fragment molecular orbital method (FMO) has been used with a large number of wave functions for single-point calculations, and its high accuracy in comparison to ab initio methods has been well established. We have developed the analytic derivative of the electrostatic interaction between far separated fragments and performed a number of restricted Hartree-Fock (RHF) geometry optimizations using FMO and ab initio methods. In particular, the alpha-helix, beta-turn, and extended conformers of a 10-residue polyalanine were studied and the good FMO accuracy was established (the rms deviations for the former two forms were about 0.2 A and for the latter structure about 0.001 A). Met-enkephalin dimer was used as a model for the polypeptide binding and computed at the 3-21G and 6-31G* levels with a similar accuracy achieved; the error in the binding energy predictions (FMO vs ab initio) was 1-3 kcal/mol. Chignolin (PDB: 1uao) and an agonist polypeptide of the erythropoietin receptor protein (emp1) were optimized at the 3-21(+)G level, with the rms deviation from ab initio of about 0.2 A, or 0.5 degrees in terms of bond angles. The effect of solvation on the structure optimization was studied in chignolin and the Trp-cage miniprotein construct (PDB:1l2y), by describing water with TIP3P. The computed structures in gas phase and solution are compared to each other and experiment.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan.
| | | | | | | |
Collapse
|
36
|
Komeiji Y, Ishida T, Fedorov DG, Kitaura K. Change in a protein's electronic structure induced by an explicit solvent: Anab initio fragment molecular orbital study of ubiquitin. J Comput Chem 2007; 28:1750-62. [PMID: 17340606 DOI: 10.1002/jcc.20686] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
The effect of solvation on the electronic structure of the ubiquitin protein was analyzed using the ab initio fragment molecular orbital (FMO) method. FMO calculations were performed for the protein in vacuo, and the protein was immersed in an explicit solvent shell as thick as 12 A at the HF or MP2 level by using the 6-31G* basis set. The protein's physical properties examined were the net charge, the dipole moment, the internal energy, and the solvent interaction energy. Comparison of the computational results revealed the following changes in the protein upon solvation. First, the positively charged amino acid residues on the protein surface drew electrons from the solvent, while the negatively charged ones transfer electrons to the solvent. Second, the dipole moment of the protein was enhanced as a result of the polarization. Third, the internal energy of the protein was destabilized, but the destabilization was more than compensated for by the generation of a favorable protein-solvent interaction. Finally, the energetic changes were elicited both by the electron correlation effect of the first solvent shell and by the electrostatic effect of more distant solvent molecules. These findings were consistent with the picture of the solvated protein being a polarizable molecule dissolved in a dielectric media.
Collapse
Affiliation(s)
- Yuto Komeiji
- Research Institute for Computational Sciences, AIST Tsukuba Central 2, Tsukuba 305-8568, Japan.
| | | | | | | |
Collapse
|
37
|
NAKANO T, MOCHIZUKI Y, AMARI S, KOBAYASHI M, FUKUZAWA K, TANAKA S. Application of Fragment Molecular Orbital (FMO) Method to Nano-Bio Field. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2007. [DOI: 10.2477/jccj.6.173] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Tatsuya NAKANO
- Division of Medicinal Safety Science, National Institute of Health Sciences
- CREST Project, Japan Science and Technology Agency
| | - Yuji MOCHIZUKI
- CREST Project, Japan Science and Technology Agency
- Department of Chemistry, Faculty of Science, Rikkyo University
| | - Shinji AMARI
- Institute of Industrial Science, The University of Tokyo
| | - Masato KOBAYASHI
- AdvanceSoft, Center for Collaborative Research, The University of Tokyo
| | - Kaori FUKUZAWA
- CREST Project, Japan Science and Technology Agency
- Mizuho Information and Research Institute, Inc
| | - Shigenori TANAKA
- CREST Project, Japan Science and Technology Agency
- Graduate School of Science and Technology, Kobe University
| |
Collapse
|
38
|
FUKUZAWA K, NAKANO T, KATO A, MOCHIZUKI Y, TANAKA S. Applications of the Fragment Molecular Orbital Method for Bio-Macromolecules. ACTA ACUST UNITED AC 2007. [DOI: 10.2477/jccj.6.185] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
39
|
Ryde U. Accurate metal-site structures in proteins obtained by combining experimental data and quantum chemistry. Dalton Trans 2006:607-25. [PMID: 17268593 DOI: 10.1039/b614448a] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The use of molecular mechanics calculations to supplement experimental data in standard X-ray crystallography and NMR refinements is discussed and it is shown that structures can be locally improved by the use of quantum chemical calculations. Such calculations can also be used to interpret the structures, e.g. to decide the protonation state of metal-bound ligands. They have shown that metal sites in crystal structures are frequently photoreduced or disordered, which makes the interpretation of the structures hard. Similar methods can be used for EXAFS refinements to obtain a full atomic structure, rather than a set of metal-ligand distances.
Collapse
Affiliation(s)
- Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, SE-221 00, Lund, Sweden.
| |
Collapse
|
40
|
Fedorov DG, Kitaura K. The three-body fragment molecular orbital method for accurate calculations of large systems. Chem Phys Lett 2006. [DOI: 10.1016/j.cplett.2006.10.052] [Citation(s) in RCA: 66] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
41
|
Dahlke EE, Truhlar DG. Electrostatically Embedded Many-Body Expansion for Large Systems, with Applications to Water Clusters. J Chem Theory Comput 2006; 3:46-53. [DOI: 10.1021/ct600253j] [Citation(s) in RCA: 231] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Erin E. Dahlke
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
42
|
Gunner MR, Mao J, Song Y, Kim J. Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2006; 1757:942-68. [PMID: 16905113 PMCID: PMC2760439 DOI: 10.1016/j.bbabio.2006.06.005] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/07/2006] [Revised: 06/07/2006] [Accepted: 06/13/2006] [Indexed: 11/15/2022]
Abstract
A protein structure should provide the information needed to understand its observed properties. Significant progress has been made in developing accurate calculations of acid/base and oxidation/reduction reactions in proteins. Current methods and their strengths and weaknesses are discussed. The distribution and calculated ionization states in a survey of proteins is described, showing that a significant minority of acidic and basic residues are buried in the protein and that most of these remain ionized. The electrochemistry of heme and quinones are considered. Proton transfers in bacteriorhodopsin and coupled electron and proton transfers in photosynthetic reaction centers, 5-coordinate heme binding proteins and cytochrome c oxidase are highlighted as systems where calculations have provided insight into the reaction mechanism.
Collapse
Affiliation(s)
- M R Gunner
- Physics Department City College of New York, New York, NY 10031, USA.
| | | | | | | |
Collapse
|
43
|
Abstract
The energy decomposition analysis (EDA) by Kitaura and Morokuma was redeveloped in the framework of the fragment molecular orbital method (FMO). The proposed pair interaction energy decomposition analysis (PIEDA) can treat large molecular clusters and the systems in which fragments are connected by covalent bonds, such as proteins. The interaction energy in PIEDA is divided into the same contributions as in EDA: the electrostatic, exchange-repulsion, and charge transfer energies, to which the correlation (dispersion) term was added. The careful comparison to the ab initio EDA interaction energies for water clusters with 2-16 molecules revealed that PIEDA has the error of at most 1.2 kcal/mol (or about 1%). The analysis was applied to (H2O)1024, the alpha helix, beta turn, and beta strand of polyalanine (ALA)10, as well as to the synthetic protein (PDB code 1L2Y) with 20 residues. The comparative aspects of the polypeptide isomer stability are discussed in detail.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- National Institute of Advanced Industrial Science and Technology (AIST), AIST Tsukuba Central 2, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan.
| | | |
Collapse
|