Calculating solvation energies by means of a fluctuating charge model combined with continuum solvent model

Continuum solvent models have shown to be very efficient for calculating solvation energy of biomolecules in solution. However, in order to produce accurate results, besides atomic radii or volumes, an appropriate set of partial charges of the molecule is needed. Here, a set of partial charges produced by a fluctuating charge model-the atom-bond electronegativity equalization method model (ABEEMσπ) fused into molecular mechanics is used to fit for the analytical continuum electrostatics model of generalized-Born calculations. Because the partial atomic charges provided by the ABEEMσπ model can well reflect the polarization effect of the solute induced by the continuum solvent in solution, accurate and rapid calculations of the solvation energies have been performed for series of compounds involving 105 small neutral molecules, twenty kinds of dipeptides and several protein fragments. The solvation energies of small neutral molecules computed with the combination of the GB model with the fluctuating charge protocol (ABEEMσπ∕GB) show remarkable agreement with the experimental results, with a correlation coefficient of 0.97, a slope of 0.95, and a bias of 0.34 kcal∕mol. Furthermore, for twenty kinds of dipeptides and several protein fragments, the results obtained from the analytical ABEEMσπ∕GB model calculations correlate well with those from ab initio and Poisson-Boltzmann calculations. The remarkable agreement between the solvation energies computed with the ABEEMσπ∕GB model and PB model provides strong motivation for the use of ABEEMσπ∕GB solvent model in the simulation of biochemical systems.