Accurate, robust and reliable calculations of Poisson-Boltzmann binding energies

Abstract

Poisson-Boltzmann (PB) model is one of the most popular implicit solvent models in biophysical modeling and computation. The ability of providing accurate and reliable PB estimation of electrostatic solvation free energy, $\Delta G_{\text{el}}$, and binding free energy, $\Delta \Delta G_{\text{el}}$, is of tremendous significance to computational biophysics and biochemistry. Recently, it has been warned in the literature (Journal of Chemical Theory and Computation 2013, 9, 3677-3685) that the widely used grid spacing of $0.5$ Å produces unacceptable errors in $\Delta \Delta G_{\text{el}}$ estimation with the solvent exclude surface (SES). In this work, we investigate the grid dependence of our PB solver (MIBPB) with SESs for estimating both electrostatic solvation free energies and electrostatic binding free energies. It is found that the relative absolute error of $\Delta G_{\text{el}}$ obtained at the grid spacing of $1.0$ Å compared to $\Delta G_{\text{el}}$ at $0.2$ Å averaged over 153 molecules is less than 0.2\%. Our results indicate that the use of grid spacing $0.6$ Å ensures accuracy and reliability in $\Delta \Delta G_{\text{el}}$ calculation. In fact, the grid spacing of $1.1$ Å appears to deliver adequate accuracy for high throughput screening.

26 pages, 7 figures

References in corpus (3)

• Second order Method for Solving 3D Elasticity Equations with Complex and Sharp Interfaces
• A GPU-accelerated Direct-sum Boundary Integral Poisson-Boltzmann Solver
• Parallel Higher-order Boundary Integral Electrostatics Computation on Molecular Surfaces with Curved Triangulation

Cited by in corpus (1)

• Some challenges of diffused interfaces in implicit-solvent models