Precision of free energies calculated by molecular dynamics simulations of peptides in solution

Abstract Errors in free energies for molecular replacement and for conformation change of a small model peptide have been determined empirically by repeated simulations from different starting points. All calculations have been done using thermodynamic integration, in which the system's potential energy is coupled to a parameter λ, that is increased or decreased by a small amount at each step of the simulation. The effects of several factors that may alter the precision are evaluated. These factors include: the length of the simulation, the dependence of the potential energy on λ, the use of conformational restraints, and their magnitude and form. The methods used for restraint and conformational forcing are described in detail. The free energy change, calculated as the mean from several successive simulations with alternately increasing and decreasing λ, is found to be independent of the length of the simulations. As expected, longer simulations produce more precise results. The variation of the calculated free energies is found to consist of two parts, a random error and a systematic hysteresis, i.e., a dependence on the direction in which λ changes. The hysteresis varies as the inverse of the length of the simulation and the random error as the inverse square root The advantage of the use of a different (nonlinear) dependence of the attractive and repulsive parts of the nonbonded potential energy on the coupling parameter when “creating” particles in solution is found to be very large. This nonlinear coupling was found to be superior to the use of linear coupling and a nonlinear change of the coupling parameter with the simulation time. The hysteresis in conformational free energy calculations is found to increase markedly if too weak a forcing restraint is chosen. It is shown how to deconvolute the contribution of a torsional restraint from the dependence of the free energy on a torsion angle.