Free Energy Landscapes of Alanine Dipeptide in Explicit Water Reproduced by the Force-Switching Wolf Method

J Chem Theory Comput. 2011 May 10;7(5):1484-93. doi: 10.1021/ct100357p. Epub 2011 Mar 30.

Abstract

Precise and rapid calculation of long-range interactions is of crucial importance for molecular dynamics (MD) and Monte Carlo simulations. Instead of the Ewald method or its high speed variant, PME, we applied our novel method, called the force-switching Wolf method, to computation of the free energy landscapes of a short peptide in explicit water. Wolf and co-workers showed that long-range electrostatic energy under a periodic boundary condition can be well reproduced even by truncating the contribution from the distant charges, when the charge neutrality is taken into account. We recently applied the procedure proposed by Wolf and co-workers to a mathematically consistent MD theory by means of a force-switching scheme, and we show that the total electrostatic energy for sodium chloride liquid was well conserved and stable during the MD simulation with the force-switching Wolf method. Our current results for an aqueous peptide solution with a series of canonical and multicanonical molecular dynamics simulations show that the force-switching Wolf method is not only in good accordance with the energies and forces calculated by the conventional PME method but also properly reproduces the solvation and the free energy landscapes of the peptide at 300 K.