A severe artifact in simulation of liquid water using a long cut-off length: Appearance of a strange layer structure

We report that a severe artifact appeared in molecular dynamics simulation of bulk water using the long cut-off length 18 A. Our result shows that increasing the cut-off length does not always improve the simulation result. Moreover, the use of the long cut-off length can lead to a spurious result. It is suggested that the simulation of solvated biomolecules using such a long cut-off length, which has been often performed, may contain an unexpected artifact.

[1]  D Peter Tieleman,et al.  Interactions of the designed antimicrobial peptide MB21 and truncated dermaseptin S3 with lipid bilayers: molecular-dynamics simulations. , 2003, The Biochemical journal.

[2]  Lennart Nilsson,et al.  Structure and dynamics of liquid water with different long‐range interaction truncation and temperature control methods in molecular dynamics simulations , 2002, J. Comput. Chem..

[3]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[4]  K. Tasaki,et al.  Observations concerning the treatment of long‐range interactions in molecular dynamics simulations , 1993, J. Comput. Chem..

[5]  W F van Gunsteren,et al.  An NMR‐based molecular dynamics simulation of the interaction of the lac repressor headpiece and its operator in aqueous solution , 1989, Proteins.

[6]  O. Steinhauser,et al.  Cutoff size does strongly influence molecular dynamics results on solvated polypeptides. , 1992, Biochemistry.

[7]  B. Brooks,et al.  Effect of Electrostatic Force Truncation on Interfacial and Transport Properties of Water , 1996 .

[8]  I. Shrivastava,et al.  K(+) versus Na(+) ions in a K channel selectivity filter: a simulation study. , 2002, Biophysical journal.

[9]  A Elofsson,et al.  Study of the electrostatics treatment in molecular dynamics simulations , 1999, Proteins.

[10]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[11]  Per Linse,et al.  Molecular dynamics simulations of polarizable water at different boundary conditions , 2000 .

[12]  B. Montgomery Pettitt,et al.  Structural and energetic effects of truncating long ranged interactions in ionic and polar fluids , 1985 .

[13]  Minoru Saito,et al.  Molecular dynamics simulations of proteins in water without the truncation of long-range Coulomb interactions , 1992 .

[14]  L. Nilsson,et al.  On the truncation of long-range electrostatic interactions in DNA. , 2000, Biophysical journal.

[15]  F. Fraternali,et al.  Hydrogen-bonding propensities of sphingomyelin in solution and in a bilayer assembly: a molecular dynamics study. , 2003, Biophysical journal.

[16]  H. Kono,et al.  Solvent density and long‐range dipole field around a DNA‐binding protein studied by molecular dynamics , 2000, Proteins.

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

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

[19]  B. Montgomery Pettitt,et al.  Effects of truncating long-range interactions in aqueous ionic solution simulations , 1988 .

[20]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[21]  H. Berendsen,et al.  A systematic study of water models for molecular simulation: Derivation of water models optimized for use with a reaction field , 1998 .

[22]  M. Patra,et al.  Molecular dynamics simulations of lipid bilayers: major artifacts due to truncating electrostatic interactions. , 2003, Biophysical journal.

[23]  R. Ornstein,et al.  An evaluation of implicit and explicit solvent model systems for the molecular dynamics simulation of bacteriophage T4 lysozyme , 1994, Proteins.

[24]  T. Darden,et al.  The effect of long‐range electrostatic interactions in simulations of macromolecular crystals: A comparison of the Ewald and truncated list methods , 1993 .

[25]  B. Brooks,et al.  The effects of truncating long‐range forces on protein dynamics , 1989, Proteins.

[26]  M. Bellissent-Funel,et al.  Hydration processes in biology : theoretical and experimental approaches , 1999 .

[27]  S. Hayward,et al.  Investigation of the mechanism of domain closure in citrate synthase by molecular dynamics simulation. , 2001, Journal of molecular biology.

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

[29]  Ronald M. Levy,et al.  Computer simulations of the dielectric properties of water: Studies of the simple point charge and transferrable intermolecular potential models , 1989 .

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

[31]  C. Soares,et al.  Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies. , 2003, Biophysical journal.

[32]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .