A New and Efficient Poisson-Boltzmann Solver for Interaction of Multiple Proteins

Area of Science:

• Computational chemistry
• Biophysics
• Applied mathematics

Background:

• The Poisson-Boltzmann equation (PBE) is crucial for modeling electrostatic interactions in biological systems.
• Existing numerical methods for solving the PBE, such as finite difference and boundary element methods, face limitations in accuracy, memory management, and computational cost.
• Accurate electrostatic calculations are essential for understanding protein function and interactions.

Purpose of the Study:

• To introduce a novel semi-analytical method, PB-SAM, for solving the linearized Poisson-Boltzmann equation.
• To demonstrate the superior accuracy, memory efficiency, and cost-effectiveness of PB-SAM compared to existing numerical PBE solvers.
• To validate the PB-SAM approach using new benchmarks and a complex biological system.

Main Methods:

• Representing protein surfaces as collections of spheres.
• Utilizing analytical multipole methods for iterative charge solving.
• Developing the Poisson Boltzmann semi-analytical method (PB-SAM).
• Establishing new benchmarks involving arrays of spherical geometries and overlapping spheres with asymmetric charges.
• Applying PB-SAM to compute the potential profile of a viral protein assembly.

Main Results:

• PB-SAM achieves higher accuracy in solving the linearized PBE.
• The method offers more flexible memory management and reduced computational cost.
• New benchmarks confirm the accuracy of PB-SAM, handling mutual polarization exactly.
• PB-SAM successfully computed the electrostatic potential for a 60-monomer viral protein array.

Conclusions:

• PB-SAM presents a significant advancement in solving the linearized PBE for biomolecular systems.
• The method provides a more efficient and accurate alternative to traditional numerical approaches.
• PB-SAM's capabilities are demonstrated for complex biological structures, paving the way for enhanced molecular modeling.