Electrostatic interactions in protein solution--a comparison between Poisson-Boltzmann and Monte Carlo calculations

Area of Science:

• Computational chemistry
• Biophysics
• Electrochemistry

Background:

• The Poisson-Boltzmann (PB) equation is a fundamental tool for modeling electrostatic interactions in biological systems.
• Understanding ion binding to proteins is crucial for predicting protein function and behavior.
• Monte Carlo (MC) simulations offer a detailed, atomistic approach to studying these interactions.

Purpose of the Study:

• To evaluate the accuracy of the Poisson-Boltzmann (PB) approximation and its linearized version.
• To compare PB and MC simulation results against experimental data for protein-ion binding.
• To assess the performance of PB models under varying salt concentrations and protein mutations.

Main Methods:

• The study employed the Poisson-Boltzmann equation solved via finite difference methods in spherical polar coordinates.
• Monte Carlo simulations were conducted for comparison.
• The calcium binding constant of the protein calbindin served as a test case, with the protein modeled as charged and neutral spheres in an electrolyte solution.

Main Results:

• Both the PB approximation and MC simulations showed quantitative agreement with experimental results for calbindin's calcium binding.
• The linearized PB equation demonstrated high accuracy but showed reduced reliability with divalent ions.
• The linearized PB model failed to accurately describe cation and anion concentration profiles, even in 1:1 salt solutions.

Conclusions:

• The PB approximation provides a reliable framework for studying protein electrostatics and ion binding.
• The linearized PB equation is a useful but limited approximation, particularly for systems with multivalent ions.
• Accurate modeling of ion concentration profiles requires careful consideration of PB model limitations.