Comparison of different schemes to treat long‐range electrostatic interactions in molecular dynamics simulations of a protein crystal

Eight molecular dynamics simulations of a ubiquitin crystal unit cell were performed to investigate the effect of different schemes to treat the long‐range electrostatic interactions as well as the need to include counter ions. A crystal system was chosen as the test system, because the higher charge density compared with a protein in solution makes it more sensitive to the way of treating the electrostatic interactions. Three different schemes of treating the long‐range interactions were compared: straight cutoff, reaction‐field approximation, and a lattice‐sum method (P3M). For each of these schemes, two simulations were performed, one with and one without the counter ions. Two additional simulations with a reaction‐field force and different initial placements of the counter ions were performed to examine the effect of the initial positions of the ions. The inclusion of long‐range electrostatic interactions using either a reaction‐field or a lattice‐sum method proved to be necessary for the simulation of crystals. These two schemes did not differ much in their ability to reproduce the crystallographic structure. The inclusion of counter ions, on the other hand, seems not necessary for obtaining a stable simulation. The initial positions of the ions have a visible but small effect on the simulation. Proteins 2001;44:509–519. © 2001 Wiley‐Liss, Inc.

[1]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[2]  E. Pietsch,et al.  Gmelin's Handbuch Der Anorganischen Chemie , 1956 .

[3]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1977, Journal of molecular biology.

[4]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[5]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1978, Archives of biochemistry and biophysics.

[6]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[7]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

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

[9]  C. Brooks.  Computer simulation of liquids , 1989 .

[10]  H. Berendsen,et al.  COMPUTER-SIMULATION OF MOLECULAR-DYNAMICS - METHODOLOGY, APPLICATIONS, AND PERSPECTIVES IN CHEMISTRY , 1990 .

[11]  H. Schreiber,et al.  Molecular dynamics studies of solvated polypeptides: Why the cut-off scheme does not work , 1992 .

[12]  Wilfred F. van Gunsteren,et al.  Consistent dielectric properties of the simple point charge and extended simple point charge water models at 277 and 300 K , 1994 .

[13]  P. Kollman,et al.  Molecular Dynamics Simulations on Solvated Biomolecular Systems: The Particle Mesh Ewald Method Leads to Stable Trajectories of DNA, RNA, and Proteins , 1995 .

[14]  Wilfred F. van Gunsteren,et al.  Determination of force field parameters for molecular simulation by molecular simulation: An application of the weak-coupling method , 1995 .

[15]  P. Argos,et al.  Knowledge‐based protein secondary structure assignment , 1995, Proteins.

[16]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[17]  P. Auffinger,et al.  A simple test for evaluating the truncation effects in simulations of systems involving charged groups , 1995 .

[18]  P. Kollman,et al.  The application of different solvation and electrostatic models in molecular dynamics simulations of ubiquitin: How well is the x‐ray structure “maintained”? , 1996, Proteins.

[19]  Alan E. Mark,et al.  The GROMOS96 Manual and User Guide , 1996 .

[20]  R. Wade,et al.  Importance of explicit salt ions for protein stability in molecular dynamics simulation. , 1998, Biophysical journal.

[21]  M. Deserno,et al.  HOW TO MESH UP EWALD SUMS. II. AN ACCURATE ERROR ESTIMATE FOR THE PARTICLE-PARTICLE-PARTICLE-MESH ALGORITHM , 1998, cond-mat/9807100.

[22]  Andrew E. Torda,et al.  The GROMOS biomolecular simulation program package , 1999 .

[23]  Gerhard Hummer,et al.  Simulation and Theory of Electrostatic Interactions in Solution: Computational Chemistry, Biophysics and Aqueous Solutions, Santa Fe, New Mexico, U. S. A., 23-25 June 1999 , 1999 .

[24]  Philippe H. Hünenberger,et al.  Lattice-sum methods for computing electrostatic interactions in molecular simulations , 1999 .

[25]  Computer simulation of biomolecules , 2000 .

[26]  J. Mccammon,et al.  Molecular Dynamics Simulations of a Polyalanine Octapeptide under Ewald Boundary Conditions: Influence of Artificial Periodicity on Peptide Conformation , 2000 .