Jw. Pitera et Wf. Van Gunsteren, One-step perturbation methods for solvation free energies of polar solutes, J PHYS CH B, 105(45), 2001, pp. 11264-11274
The one-step perturbation approach to free energy calculations permits the
calculation of relative free energies for a number of end states from a sin
gle simulation of a carefully chosen reference state. This approach has suc
cessfully been used for the calculation of relative and absolute solvation
free energies of nonpolar solutes in water, as well as relative binding fre
e energies for a family of nonpolar protein ligands. We have investigated s
everal possible approaches to extend one-step perturbation methods to also
obtain accurate solvation free energies of polar species. A "soft dipole" r
eference state permits the determination of accurate relative solvation fre
e energies for polar solutes (average unsigned error < 2.0 kJ/mol) but does
not yield accurate absolute solvation free energies. However, contrary to
expectations, simulations of a neutral reference state even yield accurate
(average unsigned error < 2.9 kJ/mol) absolute solvation free energies for
polar solutes when translational and rotational sampling of the solute is i
ncluded. In general, the choice of an appropriate reference state for one-s
tep perturbation is a challenge because of complex many-body effects. Nonet
heless, the onestep perturbation method combined with an unphysical referen
ce state yields accurate free energy estimates with an efficiency several t
imes greater than traditional approaches.