Modeling Charge-Sign Asymmetric Solvation Free Energies With Nonlinear Boundary Conditions

Abstract

We show that charge-sign-dependent asymmetric hydration can be modeled accurately using linear Poisson theory but replacing the standard electric-displacement boundary condition with a simple nonlinear boundary condition. Using a single multiplicative scaling factor to determine atomic radii from molecular dynamics Lennard-Jones parameters, the new model accurately reproduces MD free-energy calculations of hydration asymmetries for (i) monatomic ions, (ii) titratable amino acids in both their protonated and unprotonated states, and (iii) the Mobley "bracelet" and "rod" test problems [J. Phys. Chem. B, v. 112:2408, 2008]. Remarkably, the model also justifies the use of linear response expressions for charging free energies. Our boundary-element method implementation demonstrates the ease with which other continuum-electrostatic solvers can be extended to include asymmetry.

7 pages, 2 figures, accepted to Journal of Chemical Physics

References in corpus (1)

• Modeling Charge-Sign Asymmetric Solvation Free Energies With Nonlinear Boundary Conditions

Cited by in corpus (5)

• Dielectric continuum methods for quantum chemistry
• Modeling Charge-Sign Asymmetric Solvation Free Energies With Nonlinear Boundary Conditions
• Accurate Evaluation of Charge Asymmetry in Aqueous Solvation
• A Nonlinear Boundary Condition for Continuum Models of Biomolecular Electrostatics
• A Simple Electrostatic Model for the Hard-Sphere Solute Component of Nonpolar Solvation