Continuum solvation model: Computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation Original Research Article

A rigorous formulation of the solvation forces (first derivatives) associated with the electrostatic free energy calculated from numerical solutions of the linearized Poisson-Boltzmann equation on a discrete grid is described. The solvation forces are obtained from the formal solution of the linearized Poisson-Boltzmann equation written in terms of the Green function. An intermediate region for the solute-solvent dielectric boundary is introduced to yield a continuous solvation free energy and accurate solvation forces. A series of numerical tests show that the calculated forces agree extremely well with finite-difference derivatives of the solvation free energy. To gain a maximum efficiency, the nonpolar contribution to the free energy is expressed in terms of the discretized grid used for the electrostatic problem. The current treatment of solvation forces can be used to introduce the influence of a continuum solvation model in molecular mechanics calculations of large biological systems.