Effect of artificial periodicity in simulations of biomolecules under Ewald boundary conditions: a continuum electrostatics study. |
| |
Authors: | P H Hünenberger J A McCammon |
| |
Affiliation: | Department of Chemistry and Biochemistry, University of California at San Diego, La Jolla 92093-0365, USA. phunenbe@ucsd.edu |
| |
Abstract: | Ewald and related methods are nowadays routinely used in explicit-solvent simulations of biomolecules, although they impose an artificial periodicity in systems which are inherently non-periodic. The consequences of this approximation should be assessed, since they may crucially affect the reliability of computer simulations under Ewald boundary conditions. In the present study we use a method based on continuum electrostatics to investigate the nature and magnitude of possible periodicity-induced artifacts on the potentials of mean force for conformational equilibria in biomolecules. Three model systems and pathways are considered: polyalanine oligopeptides (unfolding), a DNA tetranucleotide (separation of the strands), and the protein Sac7d (conformations from a molecular dynamics simulation). Artificial periodicity may significantly affect these conformational equilibria, in each case stabilizing the most compact conformation of the biomolecule. Three factors enhance periodicity-induced artifacts: (i) a solvent of low dielectric permittivity; (ii) a solute size which is non-negligible compared to the size of the unit cell; and (iii) a non-neutral solute. Neither the neutrality of the solute nor the absence of charge pairs at distances exceeding half the edge of the unit cell do guarantee the absence of artifacts. |
| |
Keywords: | |
|
|