The Effect of Long-range Electrostatic Interactions in Simulations of Macromolecular Crystals: A Comparison of the Ewald and Truncated List Methods

The Journal of Chemical Physics vol. 99  p. 8345  DOI: 10.1063/1.465608  Published: 1993-11-15 

Darrin M. York [ ] , Tom A. Darden, Lee G. Pedersen

  View Full Article
 Download PDF


Simulations of the HIV-1 protease unit cell using a 9 Å cutoff, 9/18 Å ‘‘twin-range’’ cutoff, and full Ewald sums have been carried out to 300 ps. The results indicate that long-range electrostaticinteractions are essential for proper representation of the HIV-1 protease crystal structure. The 9 Å simulation did not converge in 300 ps. Inclusion of a 9/18 Å ‘‘twin-range’’ cutoff showed significant improvement. Simulation using the Ewald summation convention gave the best overall agreement with x-ray crystallographic data, and showed the least internal differences in the time average structures of the asymmetric units. The Ewald simulation represents an efficient implementation of the Particle Mesh Ewald method [Darden et al., J. Chem. Phys. 98, 10?089 (1993)], and illustrates the importance of including long-range electrostatic forces in large macromolecular systems.