A simplified confinement method for calculating absolute free energies and free energy and entropy differences

A simple and robust formulation of the path-independent confinement method for the calculation of free energies is presented. The simplified confinement method (SCM) does not require matrix diagonalization or switching off the molecular force field, and has a simple convergence criterion. The method can be readily implemented in molecular dynamics programs with minimal or no code modifications. Because the confinement method is a special case of thermodynamic integration, it is trivially parallel over the integration variable. The accuracy of the method is demonstrated using a model diatomic molecule, for which exact results can be computed analytically. The method is then applied to the alanine dipeptide in vacuum, and to the α-helix ↔ β-sheet transition in a 16-residue peptide modeled in implicit solvent. The SCM requires less effort for the calculation of free energy differences than previous formulations because it does not require computing normal modes. The SCM has a diminished advantage for determining absolute free energy values, because it requires decreasing the MD integration step to obtain accurate results. An approximate confinement procedure is introduced, which can be used to estimate directly the configurational entropy difference between two macrostates, without the need for additional computation of the difference in the free energy or enthalpy. The approximation has convergence properties similar to those of the standard confinement method for the calculation of free energies. The use of the approximation requires about 5 times less wall-clock simulation time than that needed to compute enthalpy differences to similar precision from an MD trajectory. For the biomolecular systems considered in this study, the errors in the entropy approximation are under 10%. Practical applications of the methods to proteins are currently limited to implicit solvent simulations.