Evaluation of poisson solvation models using a hybrid explicit/implicit solvent method

Implicit solvent methods have become popular tools in the field of protein dynamics simulations, yet evaluation of their validity has been primarily limited to comparisons with experimental and theoretical data for small molecules. In this paper, we use a recently developed hybrid explicit/implicit solvent methodology to evaluate the accuracy of several Poisson-based implicit solvent models. Specifically, we focus on the calculation of electrostatic solvation free energies of various fixed conformations for two proteins. We show that, among various dielectric boundary definitions, the Lee-Richards molecular surface has the best agreement with hybrid solvent results. Furthermore, certain modifications of the molecular surface Poisson protocol provide varied results. For instance, simple modifications of atomic radii on charged residues generally improve absolute errors but do not significantly reduce relative errors among conformations. On the other hand, using a waterprobe radius of 1.0 angstrom, as opposed to the standard value of 1.4 angstrom, to generate the molecular surface, moderately improves both absolute and relative results.