An analytical approach to computing biomolecular electrostatic potential. II. Validation and applications.

Area of Science:

• Computational chemistry
• Biophysics
• Molecular modeling

Background:

• Efficient computation of electrostatic potential is crucial for molecular simulations.
• Realistic solvation conditions and mobile ions (Debye-Hückel) are essential for accuracy.
• Existing methods can be computationally intensive.

Purpose of the Study:

• To test a closed-form analytical approximation of the Poisson equation on realistic biomolecular shapes.
• To assess the accuracy of this approximation against numerical Poisson-Boltzmann (NPB) solutions.
• To evaluate the computational efficiency and applicability of the model.

Main Methods:

• Analytical approximation derived from Poisson equation for spherical geometry, extended to realistic shapes.
• Inclusion of mobile ions using Debye-Hückel theory.
• Validation against NPB reference solutions on 580 biomolecular structures.

Main Results:

• The approximation shows high accuracy, with average deviations within 0.6 kcal/mol/(kT/e) for most structures.
• Accuracy decreases in deep, narrow indentations on molecular surfaces but improves with smoothed boundaries.
• The method is computationally inexpensive, enabling analysis of large systems like viral capsids.

Conclusions:

• The analytical approximation provides a computationally efficient and accurate method for calculating electrostatic potential in biomolecules.
• The model is suitable for large-scale simulations and can reveal potential binding sites, as demonstrated with a viral capsid example.
• Further investigation into continuum dielectric models for highly curved regions is warranted.