Treatment of geometric singularities in implicit solvent models

Area of Science:

• Computational chemistry
• Biophysics
• Numerical analysis

Background:

• Geometric singularities in molecular surfaces hinder accurate Poisson-Boltzmann (PB) equation solutions.
• Existing PB solvers, like MIBPB-I, struggle with convergence and stability at these singularities.
• The matched interface and boundary (MIB) method offers accurate biomolecular electrostatics but faces limitations with complex geometries.

Purpose of the Study:

• To develop a robust Poisson-Boltzmann solver (MIBPB-II) that accurately handles geometric singularities.
• To achieve second-order convergence for arbitrarily complex molecular surfaces.
• To improve the symmetry and diagonal dominance of the solver's matrix for enhanced stability.

Main Methods:

• Introduced a new interface-based PB solver, MIBPB-II.
• Developed a systematic approach to treat geometric singularities like cusps and self-intersections.
• Implemented a procedure to ensure optimal matrix symmetry and diagonal dominance.

Main Results:

• MIBPB-II demonstrates robustness and achieves second-order convergence for complex protein surfaces.
• The solver yields accurate electrostatic potentials and solvation free energies at a coarse 0.5 Å grid spacing.
• Results show superior accuracy compared to PBEQ and APBS at finer resolutions.

Conclusions:

• MIBPB-II effectively addresses limitations of previous PB solvers concerning geometric singularities.
• The new method enables highly accurate and stable biomolecular electrostatic calculations for complex molecular structures.
• This advancement is crucial for precise modeling in computational biology and drug discovery.