Development and parametrization of continuum solvent models. I. Models based on the boundary element method

A series of different simplifications of the boundary element method (BEM) for solving the Poisson–Boltzmann equation is investigated in an effort to obtain an accurate and fast enough treatment of electrostatic effects to be incorporated in Monte‐Carlo and molecular dynamics simulation methods. The tested simplifications include increasing the size of Boundary Elements, decreasing the surface dot density, and ignoring the interactions between the polarization charges. Combined with terms describing the nonelectrostatic solvation effects, the simplified BEM polarization terms were built into expressions for the solvation potential. The solvation potential is treated as empirical consistent force field equations. The intervening parameters, including atomic and probe radii, are derived by different fitting strategies of calculated vs experimental vacuum to water transfer energies of 173 charged, polar, and nonpolar small molecules. These fits are shown to yield very good correlations (rms ∼1.4 kcal/mol), e...