Treating ion distribution with Gaussian-based smooth dielectric function in DelPhi.

Area of Science:

• Computational chemistry
• Biophysics
• Molecular modeling

Background:

• Standard Poisson-Boltzmann equation (PBE) models rely on molecular surfaces and Stern layers to define ion presence.
• Gaussian-based smooth dielectric functions in DelPhi use smooth solute-solvent boundaries, posing challenges for modeling mobile ions.
• A clear method for handling mobile ions without distinct solute-solvent boundaries is needed.

Discussion:

• This work extends the Gaussian-based smooth dielectric function to incorporate mobile ions via Boltzmann distribution.
• A desolvation penalty is introduced, influencing ion concentration based on local electrostatics and partial desolvation.
• The model's effectiveness is validated against experimental salt dependence of binding free energies.

Key Insights:

• The extended Gaussian-based approach successfully models mobile ions by combining electrostatic potential and desolvation penalties.
• High correlations (0.95 for protein-protein, 0.88 for DNA-protein) were achieved when comparing model predictions to experimental data.
• This method offers a robust way to study ion behavior near biomacromolecules.

Outlook:

• Potential for improved accuracy in molecular dynamics simulations and drug design.
• Application to a wider range of biomolecular systems and solvation phenomena.
• Further refinement of the desolvation penalty for enhanced predictive power.