Long range Debye-Hückel correction for computation of grid-based electrostatic forces between biomacromolecules

Area of Science:

• Computational biophysics
• Molecular modeling

Background:

• Brownian dynamics (BD) simulations enable atomic-detail studies of large molecular systems, including intracellular environments.
• Computational cost, particularly for interaction forces and energies, is a key challenge in BD simulations.
• Grid-based methods for computing macromolecular interactions can introduce finite size errors for long-range forces like electrostatics.

Purpose of the Study:

• To implement and evaluate a Debye-Hückel correction for grid-based electrostatic potentials in BD simulations.
• To improve the accuracy of simulating large biomolecular systems, specifically protein solutions.

Main Methods:

• Implementation of a Debye-Hückel correction within the SDA BD simulation software.
• Application of the corrected method to simulate solutions of bovine serum albumin and hen egg white lysozyme.
• Utilizing precomputed interaction potentials on three-dimensional discretized grids.

Main Results:

• The Debye-Hückel correction significantly increased the accuracy of protein-protein interaction profiles.
• Improved accuracy was also observed in protein diffusion coefficients, especially at low ionic strength.
• The correction effectively addressed finite size errors in grid-based electrostatic calculations.

Conclusions:

• The Debye-Hückel correction offers a computationally efficient way to handle long-range electrostatic interactions in large biomacromolecular systems.
• This method enhances the accuracy of BD simulations for protein solutions.
• The implementation is also applicable to implicit solvent molecular dynamics simulations using gridded potentials.