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

Journal of Computational Chemistry - Tập 15 Số 7 - Trang 667-683 - 1994
Peter Steinbach1, Bernard R. Brooks1
1Laboratory of Structural Biology, Division of Computer Research & Technology, National Institutes of Health, Bethesda, Maryland 20892

Tóm tắt

AbstractNew 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 (ϵ ∝︁ r) is found to alter the temperature dependence of protein dynamics, suppressing anharmonic motion at high temperatures. Force switching and force shifting are the best atom‐based spherical cutoffs, whereas switched group‐shifting is the preferred group‐based method. To achieve realistic simulations, increasing the cutoff distance from 7.5 to 12 Å or beyond is much more important than the differences among the three best cutoff methods. © 1994 by John Wiley & Sons, Inc.This article is a US Government work and, as such, is in the public domain in the United States of America.

Từ khóa


Tài liệu tham khảo

10.1016/0022-2836(91)80181-S

10.1063/1.449621

10.1021/j100283a017

10.1063/1.452636

10.1002/prot.340060104

10.1002/jcc.540111009

10.1021/bi00140a022

10.1002/jcc.540140304

10.1002/jcc.540140306

10.1073/pnas.90.19.9135

10.1016/0009-2614(92)90053-P

10.1002/jcc.540040211

Polygen Corporation Parameter file for CHARMm version 20 Copyright1986 released August1988.

10.1016/0009-2614(88)80404-4

10.1063/1.1676585

10.1080/00268977900101751

10.1063/1.448118

10.1016/S0006-3495(91)82292-5

10.1016/S0065-3233(08)60197-7

10.1080/08927028808080956

10.1063/1.461272

10.1073/pnas.87.4.1601

10.1016/S0022-2836(05)80196-2

Allen M. P., 1987, Computer Simulation of Liquids

10.1080/00268979000101991

10.1063/1.460259

10.1063/1.445869

Brooks B. R., 1992, Chem. Des. Auto. News, 7, 16

10.1016/0022-2836(86)90470-5

10.1016/S0022-2836(05)80363-8

10.1038/337754a0

10.1063/1.463935

10.1016/0301-0104(91)87078-A

Thông tin
Thông tin xuất bản

Journal of Computational Chemistry

Tập 15 Số 7

667-683

Thông tin tác giả