Incorporating dipolar solvents with variable density in Poisson-Boltzmann electrostatics.

Area of Science:

• Computational chemistry
• Biophysics
• Molecular modeling

Background:

• Classical Poisson-Boltzmann models treat solvent as a homogeneous dielectric.
• Advanced methods are needed for accurate electrostatic property calculations of macromolecules.
• Existing models lack detailed solvent behavior representation.

Purpose of the Study:

• To develop a new computational method for calculating electrostatic properties of macromolecules.
• To model solvent as self-orienting, interacting dipoles with variable density.
• To provide a more accurate and computationally efficient alternative to current methods.

Main Methods:

• Developed a model treating solvent as an assembly of self-orienting interacting dipoles.
• Introduced a variable dielectric constant (ε(r)) and solvent density (ρ(r)).
• Calibrated the model using small molecule and ion solvation data with two adjustable parameters.

Main Results:

• The model accurately predicts hydrophobicity scales at atomic and residue levels.
• Identifies poorly solvated patches at protein dimerization interfaces and lipid-binding sites.
• Shows qualitative agreement with molecular dynamics simulations for solvent density.

Conclusions:

• The new method offers a computationally feasible approach for analyzing macromolecular electrostatic properties.
• It provides insights into protein structure, function, and solvation.
• Potential for applications in structure-based drug design and protein engineering.