Treatment of charge singularities in implicit solvent models

Area of Science:

• Computational chemistry
• Biophysics
• Applied mathematics

Background:

• Solving the Poisson-Boltzmann (PB) equation is crucial for understanding molecular electrostatics.
• Geometric singularities at dielectric interfaces pose challenges for existing PB solvers.
• Previous methods like MIBPB-II achieved good accuracy but struggled with fine meshes and charge singularities.

Purpose of the Study:

• To develop a highly accurate and efficient Poisson-Boltzmann solver.
• To rigorously handle both geometric and charge singularities in dielectric interfaces.
• To improve electrostatic potential calculations for biomolecules, especially proteins.

Main Methods:

• Utilized a Green's function formulation to address charge singularities.
• Integrated charge singularity treatment with the matched interface and boundary (MIB) method.
• Transformed charge singularities into interface flux jump conditions within the MIB framework.

Main Results:

• Developed the MIBPB-III solver, achieving high accuracy at coarse mesh sizes (1.2 Å) for proteins.
• Successfully treated geometric and charge singularities on an equal footing.
• Demonstrated MIBPB-III is approximately three times faster than APBS at comparable accuracy levels.

Conclusions:

• MIBPB-III provides reliable benchmark numerical solutions for the PB equation.
• The method offers significant computational advantages for electrostatic calculations in biomolecular systems.
• Validated extensively using diverse systems, including 24 proteins, confirming its robustness and accuracy.