Optimized radii for Poisson Boltzmann calculations with the AMBER force field

@article{Swanson:2005, abstract = {Implicit solvent models are a standard tool for assessing the electrostatics of biomolecular systems. The accuracy of quantitative predictions, such as pKa values, transfer free energies, binding energies, and solvation forces, is strongly dependent on one's choice of continuum parameters: the solute charges, dielectric coefficient, and radii, which define the dielectric boundary. To ensure quantitative accuracy, these parameters can be benchmarked against explicit solvent simulations. Here we present two sets of optimized radii to define either abrupt or cubic-spline smoothed dielectric boundaries in Poisson-Boltzmann calculations of protein systems with AMBER (parm99) charges. Spline smoothing stabilizes the electrostatic potential at the molecular surface, allowing for continuum force calculations. Most implementations, however, require significantly different radii than the abrupt boundary surfaces. The optimal continuum radii are initially approximated from the solvent radial charge distribution surrounding each atom type. A genetic algorithm is then used to fine-tune the starting values to reproduce charging free energies measured from explicit solvent simulations. The optimized radii are tested on four protein-like polypeptides. The results show increased accuracy of molecular solvation energies and atomic forces relative to commonly used continuum parameter sets. These radii are suitable for Poisson-Boltzmann calculations with the AMBER force field and offer energetic congruence to any model that combines molecular mechanics and Poisson-Boltzmann solvation energies.}, annote = {Optimized radii for Poisson Boltzmann calculations with the AMBER force field}, author = {Swanson, J. M. J. and Adcock, S. A. and Mccammon, J. A. }, citeulike-article-id = {1329584}, journal = {Journal of Chemical Theory and Computation}, keywords = {electrostatics, molecular\_simulation, solvation}, local-url = {file:///C:/Documents\%20and\%20Settings/baker/My\%20Documents/References/Archive/ct049834o.pdf}, number = {3}, pages = {484--93}, posted-at = {2007-05-25 14:14:06}, priority = {0}, title = {Optimized radii for Poisson Boltzmann calculations with the AMBER force field}, url = {http://dx.doi.org/10.1021/ct049834o}, volume = {1}, year = {2005} }

TY - JOUR ID - Swanson:2005 L3 - citeulike-article-id:1329584 N2 - Implicit solvent models are a standard tool for assessing the electrostatics of biomolecular systems. The accuracy of quantitative predictions, such as pKa values, transfer free energies, binding energies, and solvation forces, is strongly dependent on one's choice of continuum parameters: the solute charges, dielectric coefficient, and radii, which define the dielectric boundary. To ensure quantitative accuracy, these parameters can be benchmarked against explicit solvent simulations. Here we present two sets of optimized radii to define either abrupt or cubic-spline smoothed dielectric boundaries in Poisson-Boltzmann calculations of protein systems with AMBER (parm99) charges. Spline smoothing stabilizes the electrostatic potential at the molecular surface, allowing for continuum force calculations. Most implementations, however, require significantly different radii than the abrupt boundary surfaces. The optimal continuum radii are initially approximated from the solvent radial charge distribution surrounding each atom type. A genetic algorithm is then used to fine-tune the starting values to reproduce charging free energies measured from explicit solvent simulations. The optimized radii are tested on four protein-like polypeptides. The results show increased accuracy of molecular solvation energies and atomic forces relative to commonly used continuum parameter sets. These radii are suitable for Poisson-Boltzmann calculations with the AMBER force field and offer energetic congruence to any model that combines molecular mechanics and Poisson-Boltzmann solvation energies. IS - 3 JF - Journal of Chemical Theory and Computation EP - 93 TI - Optimized radii for Poisson Boltzmann calculations with the AMBER force field VL - 1 SP - 484 KW - electrostatics KW - molecular_simulation KW - solvation AU - Swanson, JMJ AU - Adcock, SA AU - Mccammon, JA PY - 2005/// UR - http://dx.doi.org/10.1021/ct049834o ER -