A Comparison of Seven Fast but Approximate Methods to Compute the Free Energy of Deprotonation for Amino Acids in Aqueous Solution

Abstract

In recent years, a variety of methods based on statistical mechanics have been successfully applied to calculate free energy differences of chemical reactions from molecular simulation. The accuracy and computational efficiency vary strongly between these methods. Seven approximate but fast methods to calculate free energy differences are compared in terms of accuracy and efficiency with the accurate but expensive thermodynamic integration method as reference, using 28 protonation and deprotonation reactions of aspartic acid in aqueous solution as test cases. At least two simulations are required to obtain an accurate free energy difference between two states of the system. Both, the averaged one-step perturbation method and the linear response method yield the most accurate results, while the latter method shows the fastest convergence.