Analysis of integral expressions for effective Born radii
- 8 November 2007
- journal article
- research article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 127 (18), 185101
- https://doi.org/10.1063/1.2783847
Abstract
Generalized Born (GB) models provide a computationally efficient means of representing the electrostatic effects of solvent and are widely used, especially in molecular dynamics (MD). Accurate and facile computation of the effective Born radii is a key for the performance of GB models. Here, we examine a simple integral prescription, R 6 , based on the exact solution of the Poisson-Boltzmann (PB) equation for a perfect sphere. Numerical tests on 22 molecules representing a variety of structural classes show that R 6 may be more accurate than the more complex integral-based approaches such as GBMV2. At the same time, R 6 is computationally less demanding. Fundamental limitations of current integration-based methods for calculating effective radii, including R 6 , are explored and the deviations from the numerical PB results are correlated with specific topological and geometrical features of the molecular surface. A small systematic bias observed in the R 6 -based radii can be removed with a single, transferable constant offset; when the resulting effective radii are used in the “classical” (Still et al.’s) GB formula to compute the electrostatic solvation free energy, the average deviation from the PB reference is no greater than when the “perfect” (PB-based) effective radii are used. This deviation is also appreciably smaller than the uncertainty of the PB reference itself, as estimated by comparison to explicit solvent.This publication has 52 references indexed in Scilit:
- Generalized Born Model with a Simple, Robust Molecular Volume CorrectionJournal of Chemical Theory and Computation, 2006
- Effective Born radii in the generalized Born approximation: The importance of being perfectJournal of Computational Chemistry, 2002
- Introducing an Implicit Membrane in Generalized Born/Solvent Accessibility Continuum Solvent ModelsThe Journal of Physical Chemistry B, 2002
- Electrostatics of nanosystems: Application to microtubules and the ribosomeProceedings of the National Academy of Sciences of the United States of America, 2001
- Generalized Born Models of Macromolecular Solvation EffectsAnnual Review of Physical Chemistry, 2000
- Generalized Born Model Based on a Surface Integral FormulationThe Journal of Physical Chemistry B, 1998
- Continuum solvation model: Computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equationComputer Physics Communications, 1998
- VMD: Visual molecular dynamicsJournal of Molecular Graphics, 1996
- Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian Dynamics programComputer Physics Communications, 1995
- Semianalytical treatment of solvation for molecular mechanics and dynamicsJournal of the American Chemical Society, 1990