Change in a protein's electronic structure induced by an explicit solvent: An ab initio fragment molecular orbital study of ubiquitin

The effect of solvation on the electronic structure of the ubiqutin 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 Å 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. © 2007 Wiley Periodicals, Inc. J Comput Chem, 2007

[1]  Noriyuki Kurita,et al.  DENSITY FUNCTIONAL CALCULATIONS ON THE INTERACTION BETWEEN CATABOLITE ACTIVATOR PROTEIN AND CYCLIC AMP USING THE FRAGMENT MOLECULAR ORBITAL METHOD , 2005 .

[2]  Kazuo Kitaura,et al.  All electron quantum chemical calculation of the entire enzyme system confirms a collective catalytic device in the chorismate mutase reaction. , 2006, The journal of physical chemistry. B.

[3]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[4]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[5]  Richard A Friesner,et al.  Structure and dynamics of the solvation of bovine pancreatic trypsin inhibitor in explicit water: a comparative study of the effects of solvent and protein polarizability. , 2005, The journal of physical chemistry. B.

[6]  Yuto Komeiji,et al.  Fragment molecular orbital method: analytical energy gradients , 2001 .

[7]  Alexander D. MacKerell,et al.  CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model , 2004, J. Comput. Chem..

[8]  M. Boero,et al.  Catalytic role of metal ion in the selection of competing reaction paths: a first principles molecular dynamics study of the enzymatic reaction in ribozyme. , 2002, Journal of the American Chemical Society.

[9]  Insights into the interplay between electronic structure and protein dynamics: The case of ubiquitin , 2005 .

[10]  Makoto Taiji,et al.  Fast and accurate molecular dynamics simulation of a protein using a special‐purpose computer , 1997 .

[11]  Kazuo Kitaura,et al.  Second order Møller-Plesset perturbation theory based upon the fragment molecular orbital method. , 2004, The Journal of chemical physics.

[12]  Yutaka Akiyama,et al.  Fragment molecular orbital method: application to polypeptides , 2000 .

[13]  Yuji Mochizuki,et al.  Large scale MP2 calculations with fragment molecular orbital scheme , 2004 .

[14]  Y. Komeiji,et al.  Molecular dynamics simulation of trp-aporepressor in a solvent. , 1991, Protein engineering.

[15]  Matteo Dal Peraro,et al.  Solute-solvent charge transfer in aqueous solution. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[16]  Kenneth M Merz,et al.  Quantum mechanics in structure-based drug design. , 2006, Current opinion in drug discovery & development.

[17]  Umpei Nagashima,et al.  A parallelized integral-direct second-order Møller–Plesset perturbation theory method with a fragment molecular orbital scheme , 2004 .

[18]  Mark S. Gordon,et al.  A new hierarchical parallelization scheme: Generalized distributed data interface (GDDI), and an application to the fragment molecular orbital method (FMO) , 2004, J. Comput. Chem..

[19]  H. Weingärtner,et al.  The Dielectric Spectrum of Ubiquitin in Aqueous Solution , 2001 .

[20]  J. Gao,et al.  A priori evaluation of aqueous polarization effects through Monte Carlo QM-MM simulations. , 1992, Science.

[21]  F. Martin,et al.  Charge distribution in the water molecule—A comparison of methods , 2005, J. Comput. Chem..

[22]  Y. Komeiji Ewald summation and multiple time step methods for molecular dynamics simulation of biological molecules , 2000 .

[23]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[24]  Kari Laasonen,et al.  ‘‘Ab initio’’ liquid water , 1993 .

[25]  V. Barone,et al.  Reliable molecular simulations of solute-solvent systems with a minimum number of solvent shells. , 2006, The Journal of chemical physics.

[26]  Fumio Hirata,et al.  The effects of solvent on the conformation and the collective motions of protein: normal mode analysis and molecular dynamics simulations of melittin in water and in vacuum , 1991 .

[27]  B. Montgomery Pettitt,et al.  Simple intramolecular model potentials for water , 1987 .

[28]  Yuichi Inadomi,et al.  Fragment molecular orbital method: application to molecular dynamics simulation, ‘ab initio FMO-MD’ , 2003 .

[29]  Y. Komeiji,et al.  A molecular dynamics study of solvent behavior around a protein , 1993, Proteins.

[30]  J. Andrew McCammon,et al.  Computer Simulation and the Design of New Biological Molecules , 1986 .

[31]  Hui Li,et al.  The polarizable continuum model (PCM) interfaced with the fragment molecular orbital method (FMO) , 2006, J. Comput. Chem..

[32]  Kazuo Kitaura,et al.  Pair interaction energy decomposition analysis , 2007, J. Comput. Chem..

[33]  Double-Metal-Ion/Single-Metal-Ion Mechanisms of the Cleavage Reaction of Ribozymes:  First-Principles Molecular Dynamics Simulations of a Fully Hydrated Model System. , 2005, Journal of chemical theory and computation.

[34]  R. Sharon,et al.  Accurate simulation of protein dynamics in solution. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[35]  Takeshi Ishikawa,et al.  Fragment molecular orbital calculations on large scale systems containing heavy metal atom , 2006 .

[36]  H. Nymeyer,et al.  Simulation of the folding equilibrium of α-helical peptides: A comparison of the generalized Born approximation with explicit solvent , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[37]  C. Bugg,et al.  Structure of ubiquitin refined at 1.8 A resolution. , 1987, Journal of molecular biology.

[38]  Petr Bour,et al.  Ab initio modeling of amide I coupling in antiparallel beta-sheets and the effect of 13C isotopic labeling on infrared spectra. , 2005, The journal of physical chemistry. B.

[39]  S. Tsuzuki,et al.  Comparison of atomic charge distributions obtained from different procedures: basis set and electron correlation effects , 1996 .

[40]  M. Klein,et al.  Reorganization free energies for long-range electron transfer in a porphyrin-binding four-helix bundle protein. , 2006, Journal of the American Chemical Society.

[41]  John Z. H. Zhang,et al.  FULL AB INITIO COMPUTATION OF PROTEIN-WATER INTERACTION ENERGIES , 2004 .

[42]  Shigenori Tanaka,et al.  Intra‐ and intermolecular interactions between cyclic‐AMP receptor protein and DNA: Ab initio fragment molecular orbital study , 2006, J. Comput. Chem..

[43]  R A Sayle,et al.  RASMOL: biomolecular graphics for all. , 1995, Trends in biochemical sciences.

[44]  M. Aida,et al.  An Ab Initio MO Study on Orbital Interaction and Charge Distribution in Alkali Metal Aqueous Solution: Li+, Na+, and K+ , 2004 .

[45]  Haibo Yu,et al.  Accounting for polarization in molecular simulation , 2005, Comput. Phys. Commun..

[46]  M Elstner,et al.  Quantum mechanics simulation of protein dynamics on long timescale , 2001, Proteins.

[47]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[48]  Salvatore Cannistraro,et al.  Molecular Dynamics of Water at the Protein-Solvent Interface , 2002 .

[49]  Yuto Komeiji,et al.  Ab initio fragment molecular orbital (FMO) method applied to analysis of the ligand-protein interaction in a pheromone-binding protein , 2005, Comput. Biol. Chem..

[50]  A. Ciechanover,et al.  Mechanisms of intracellular protein breakdown. , 1982, Annual review of biochemistry.

[51]  T. Ikeda,et al.  Hydration of Y3+ ion: a Car-Parrinello molecular dynamics study. , 2005, The Journal of chemical physics.

[52]  Kotoko Nakata,et al.  Ab initio quantum mechanical study of the binding energies of human estrogen receptor α with its ligands: An application of fragment molecular orbital method , 2005, J. Comput. Chem..

[53]  Kenneth B. Wiberg,et al.  Comparison of atomic charges derived via different procedures , 1993, J. Comput. Chem..

[54]  Alessandro Laio,et al.  A comparative theoretical study of dipeptide solvation in water , 2006, J. Comput. Chem..

[55]  Andrew S. Dixon,et al.  Charge-Transfer Interactions in Macromolecular Systems: A New View of the Protein/Water Interface , 1998 .

[56]  Kaori Fukuzawa,et al.  Fragment molecular orbital method: use of approximate electrostatic potential , 2002 .

[57]  Dean R. Haeffner,et al.  Electron distribution in water , 2000 .

[58]  Yuichi Inadomi,et al.  PEACH 4 with ABINIT-MP: a general platform for classical and quantum simulations of biological molecules. , 2004, Computational biology and chemistry.

[59]  J. Perram,et al.  Simulation of electrostatic systems in periodic boundary conditions. I. Lattice sums and dielectric constants , 1980, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[60]  Shuichi Nosé,et al.  Constant Temperature Molecular Dynamics Methods , 1991 .

[61]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[62]  Christopher M. Hadad,et al.  Comparison of different atomic charge schemes for predicting pKa variations in substituted anilines and phenols , 2002 .

[63]  R. Zhou Free energy landscape of protein folding in water: Explicit vs. implicit solvent , 2003, Proteins.