Highly accurate biomolecular electrostatics in continuum dielectric environments

Area of Science:

• Computational chemistry
• Biophysics
• Molecular modeling

Background:

• Implicit solvent models are crucial for simulating biomolecules in solution.
• The Poisson-Boltzmann (PB) equation is a standard model for these simulations.
• Existing PB solvers face challenges with accuracy and convergence.

Purpose of the Study:

• To develop a novel, highly accurate Poisson-Boltzmann solver.
• To improve the convergence rates and accuracy of electrostatic calculations for biomolecules.
• To establish new reference values for biomolecular electrostatic modeling.

Main Methods:

• Development of a Poisson-Boltzmann solver utilizing the matched interface and boundary (MIB) method.
• Rigorous enforcement of electrostatic potential and flux continuity at molecular surfaces.
• Comparison of MIB solver performance against finite difference and finite element methods.

Main Results:

• The MIB-based PB solver demonstrates significantly improved convergence rates with mesh size.
• Highly accurate electrostatic potentials and solvation energies are achieved at coarse mesh sizes.
• The MIB method yields more accurate solutions to the PB equation than established methods.

Conclusions:

• The MIB method offers a new standard for accuracy in Poisson-Boltzmann solvers.
• This approach enhances the quality of electrostatic surface potentials for biomolecular interaction studies.
• The MIB-based solver is a valuable tool for computational biophysics and molecular modeling.