Molecular Dynamics Simulations of a Polyalanine Octapeptide under Ewald Boundary Conditions: Influence of Artificial Periodicity on Peptide Conformation Export

@article{citeulike:1807941, abstract = {Ewald and related mesh methods are nowadays routinely used in explicit-solvent simulations of solvated biomolecules, although they impose an artificial periodicity in systems which are inherently nonperiodic. In the present study, we investigate the consequences of this approximation for the conformational equilibrium of a polyalanine octapeptide (with charged termini) in water. We report three explicit-solvent molecular dynamics simulations of this peptide in cubic unit cells of edges L = 2, 3, and 4 nm, using the particleparticleparticlemesh (P3M) method for handling electrostatic interactions. The initial configuration of the peptide is alpha-helical. In the largest unit cell (L = 4 nm), the helix unfolds quickly toward configurations with shorter end-to-end distances. By contrast, in the two smaller unit cells (L = 2 and 3 nm), the alpha-helix remains stable during 2 ns. Backbone fluctuations are somewhat larger in the medium (L = 3 nm) compared to the smallest unit cell. These differences are rationalized using a continuum electrostatics analysis of configurations from the simulations. These calculations show that the alpha-helical conformation is stabilized by artificial periodicity relative to any other configuration sampled during the trajectories. This artificial stabilization is larger for smaller unit cells, and is responsible for the absence of unfolding in the two smaller unit cells and the reduced backbone fluctuations in the smallest unit cell. These results suggest that artificial periodicity imposed by the use of infinite periodic (Ewald) boundary conditions in explicit-solvent simulations of biomolecules may significantly perturb the potentials of mean force for conformational equilibria, and even in some cases invert the relative stabilities of the folded and unfolded states.}, address = {Department of Chemistry and Biochemistry and Department of Pharmacology, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093-0365 and Laboratorium f\"{u}r Physikalische Chemie, ETH Zentrum, Universit\"{a}tstrasse 6, CH-8092 Z\"{u}rich, Switzerland}, author = {Weber, Wolfgang and Hunenberger, Philippe H. and Mccammon, J. Andrew}, citeulike-article-id = {1807941}, citeulike-linkout-0 = {http://dx.doi.org/10.1021/jp9937757}, citeulike-linkout-1 = {http://pubs.acs.org/doi/abs/10.1021/jp9937757}, day = {1}, doi = {10.1021/jp9937757}, journal = {The Journal of Physical Chemistry B}, keywords = {ewald, protein\_electrostatics}, month = {April}, number = {15}, pages = {3668--3675}, posted-at = {2007-10-22 22:46:06}, priority = {3}, title = {Molecular Dynamics Simulations of a Polyalanine Octapeptide under Ewald Boundary Conditions: Influence of Artificial Periodicity on Peptide Conformation}, url = {http://dx.doi.org/10.1021/jp9937757}, volume = {104}, year = {2000} }

TY - JOUR ID - citeulike:1807941 L3 - citeulike-article-id:1807941 N2 - Ewald and related mesh methods are nowadays routinely used in explicit-solvent simulations of solvated biomolecules, although they impose an artificial periodicity in systems which are inherently nonperiodic. In the present study, we investigate the consequences of this approximation for the conformational equilibrium of a polyalanine octapeptide (with charged termini) in water. We report three explicit-solvent molecular dynamics simulations of this peptide in cubic unit cells of edges L = 2, 3, and 4 nm, using the particleparticleparticlemesh (P3M) method for handling electrostatic interactions. The initial configuration of the peptide is alpha-helical. In the largest unit cell (L = 4 nm), the helix unfolds quickly toward configurations with shorter end-to-end distances. By contrast, in the two smaller unit cells (L = 2 and 3 nm), the alpha-helix remains stable during 2 ns. Backbone fluctuations are somewhat larger in the medium (L = 3 nm) compared to the smallest unit cell. These differences are rationalized using a continuum electrostatics analysis of configurations from the simulations. These calculations show that the alpha-helical conformation is stabilized by artificial periodicity relative to any other configuration sampled during the trajectories. This artificial stabilization is larger for smaller unit cells, and is responsible for the absence of unfolding in the two smaller unit cells and the reduced backbone fluctuations in the smallest unit cell. These results suggest that artificial periodicity imposed by the use of infinite periodic (Ewald) boundary conditions in explicit-solvent simulations of biomolecules may significantly perturb the potentials of mean force for conformational equilibria, and even in some cases invert the relative stabilities of the folded and unfolded states. IS - 15 JF - The Journal of Physical Chemistry B AD - Department of Chemistry and Biochemistry and Department of Pharmacology, University of California at San Diego, 9500 Gilman Drive, La Jolla, California 92093-0365 and Laboratorium für Physikalische Chemie, ETH Zentrum, Universitätstrasse 6, CH-8092 Zürich, Switzerland EP - 3675 TI - Molecular Dynamics Simulations of a Polyalanine Octapeptide under Ewald Boundary Conditions: Influence of Artificial Periodicity on Peptide Conformation VL - 104 SP - 3668 KW - ewald KW - protein_electrostatics AU - Weber, Wolfgang AU - Hunenberger, Philippe AU - Mccammon, Andrew PY - 2000/04/01/ UR - http://dx.doi.org/10.1021/jp9937757 ER -