Accurate, robust and reliable calculations of Poisson-Boltzmann binding energies
This work provides practical grid spacing guidelines for researchers using Poisson-Boltzmann models, enabling faster computations without sacrificing accuracy.
The authors demonstrate that their MIBPB solver achieves accurate and reliable Poisson-Boltzmann electrostatic binding energy calculations with coarse grid spacings (e.g., 0.6 Å for binding, 1.1 Å for screening), contrary to prior warnings that 0.5 Å spacing produces unacceptable errors. Relative absolute errors in solvation free energy at 1.0 Å spacing are less than 0.2% over 153 molecules.
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, $ΔG_{\text{el}}$, and binding free energy, $ΔΔ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 $ΔΔ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 $ΔG_{\text{el}}$ obtained at the grid spacing of $1.0$ Å$ $ compared to $Δ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 $ΔΔG_{\text{el}}$ calculation. In fact, the grid spacing of $1.1$ Å$ $ appears to deliver adequate accuracy for high throughput screening.