On removal of charge singularity in Poisson-Boltzmann equation

Area of Science:

• Computational chemistry
• Biophysics
• Theoretical chemistry

Background:

• Poisson-Boltzmann theory is crucial for modeling biomolecular electrostatics.
• Atomic point charges in force fields cause singularity issues in Poisson-Boltzmann equations.
• Discretization methods reduce numerical difficulties but introduce errors.

Purpose of the Study:

• To develop an efficient method to overcome charge singularity in Poisson-Boltzmann calculations.
• To achieve higher numerical accuracy in electrostatic potential solutions.
• To improve electrostatic calculations in molecular mechanics simulations.

Main Methods:

• Solving two separate equations for reaction field and total potentials simultaneously.
• Implementing the method within finite-difference Poisson-Boltzmann solvers.
• Validating the approach on small molecules and large proteins.

Main Results:

• The proposed method effectively removes charge singularity, yielding singularity-free reaction field potentials.
• High agreement was observed between computed reaction field energies and the classical finite-difference Poisson-Boltzmann method.
• The new method demonstrates faster convergence compared to classical approaches.

Conclusions:

• The developed method offers a more accurate and efficient solution for electrostatic solvation interactions.
• It provides a robust framework for advanced electrostatic calculations in complex biomolecular systems.
• This advancement facilitates more reliable molecular mechanics simulations.