An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations

Area of Science:

• Computational Chemistry
• Biophysics
• Theoretical Chemistry

Background:

• Accurate calculation of long-range electrostatic forces is crucial for molecular simulations of periodic systems.
• Existing methods often face computational challenges with large systems and complex interactions.

Purpose of the Study:

• To present an efficient method for evaluating long-range electrostatic forces in combined quantum mechanical and molecular mechanical (QM/MM) calculations of periodic systems.
• To enable accurate and scalable simulations of biological processes.

Main Methods:

• Developed the QM/MM-Ewald method, a linear-scaling electrostatic approach.
• Utilizes particle mesh Ewald for MM atoms and real-space multipolar expansion for QM terms.
• Incorporates a pairwise periodic correction factor for QM and QM/MM interactions.

Main Results:

• Tested the QM/MM-Ewald method in molecular dynamics simulations of ion association and phosphoryl transfer reactions.
• Compared results with traditional methods using electrostatic cutoffs and full electrostatics.
• Demonstrated the method's efficiency and accuracy in periodic boundary molecular dynamics (PBMD) simulations.

Conclusions:

• The QM/MM-Ewald method provides an efficient and accurate way to handle electrostatic interactions in periodic QM/MM calculations.
• This method facilitates the extension of linear-scaling Ewald methods to complex molecular simulations, including enzyme and ribozyme reactions.