New spherical‐cutoff methods for long‐range forces in macromolecular simulation

作者：Peter J. Steinbach, Bernard R. Brooks · 发表于：Journal of Computational Chemistry · 年份：1994 · DOI：10.1002/jcc.540150702 · 被引用次数：1177 · 研究领域：Protein Structure and Dynamics、DNA and Nucleic Acid Chemistry、Advanced NMR Techniques and Applications

Abstract New atom‐ and group‐based spherical‐cutoff methods have been developed for the treatment of nonbonded interactions in molecular dynamics (MD) simulation. A new atom‐based method, force switching, leaves short‐range forces unaltered by adding a constant to the potential energy, switching forces smoothly to zero over a specified range. A simple improvement to group‐based cutoffs is presented: Switched group‐shifting shifts the group–group potential energy by a constant before being switched smoothly to zero. Also introduced are generalizations of atom‐based force shifting, which adds a constant to the Coulomb force between two charges. These new approaches are compared to existing methods by evaluating the energy of a model hydrogen‐bonding system consisting of two N ‐methyl acetamide molecules and by full MD simulation. Thirty‐five 150 ps simulations of carboxymyoglobin (MbCO) hydrated by 350 water molecules indicate that the new methods and atom‐based shifting are each able to approximate no‐cutoff results when a cutoff at or beyond 12 Å is used. However, atom‐based potential‐energy switching and truncation unacceptably contaminate group–group electrostatic interactions. Group‐based potential truncation should not be used in the presence of explicit water or other mobile electrostatic dipoles because energy is not a state function with this method, resulting in severe heating (about 4 K/ps in the simulations of hydrated MbCO). The distance‐dependent dielectric (ϵ ∝︁ ...