A Newton-like iterative method implemented in the DelPhi for solving the nonlinear Poisson-Boltzmann equation

Area of Science:

• Computational biology
• Biophysics
• Biomolecular modeling

Background:

• DelPhi is a widely used scientific program for calculating electrostatic potentials of biomolecules.
• It solves the nonlinear Poisson-Boltzmann equation (PBE) using the finite difference method.
• The existing Successive Over Relaxation (SOR) algorithm has limitations in handling highly nonlinear PBE problems.

Purpose of the Study:

• Introduce a novel Newton-like method for solving the nonlinear PBE within the DelPhi program.
• Improve the stability and convergence of DelPhi for complex biomolecular systems.
• Provide a more robust computational tool for electrostatic analysis.

Main Methods:

• Implementation of a Newton-like iterative method to solve the nonlinear PBE.
• Comparison of the new method's performance against the Successive Over Relaxation (SOR) algorithm.
• Testing on diverse biomolecular examples to evaluate stability and convergence.

Main Results:

• The Newton-like method demonstrates superior stability compared to the SOR method for nonlinear PBE problems.
• The new method achieves convergence for cases with very strong nonlinearity, where SOR may fail.
• Enhanced reliability of electrostatic potential calculations in DelPhi.

Conclusions:

• The novel Newton-like method offers a significant advancement for DelPhi's nonlinear PBE solver.
• This improvement enhances DelPhi's capability in accurately modeling biomolecular electrostatics.
• The updated DelPhi version provides a more robust and efficient tool for biophysical research.