Calculation of Relative Solvation Free Energy Differences by Thermodynamic Perturbation Method: Dependence of Free Energy Results on Simulation Length

摘要

Molecular dynamics (MD) simulations in conjunction with the thermodynamic cycle perturbation approach has been used to calculate relative solvation free energies for acetone to acetaldehyde, acetone to pyruvic acid, acetone to 1,1,1-trifluoroacetone, acetone to 1,1,1-trichloroacetone, acetone to 2,3-butanedione, acetone to cyclopropanone, and formaldehyde hydrate to formaldehyde. To evaluate the dependence of relative solvation free energy convergence on MD simulation length and starting configuration two studies were performed. In the first study, each simulation started from the same well-equilibrated configuration and the length was varied from 153 to 1530 ps. In the second study, the relative solvation free energy differences were calculated starting from three different configurations and using 510 ps of MD simulation for each mutation. These results clearly indicate that, even for molecules with limited conformational flexibility, a simulation length of 510 ps or greater is required to obtain satisfactory convergence and, for the mutations of large structural changes between reactant and product, such as cyclopropanone to acetone, require much longer simulation lengths to achieve satisfactory convergence. These results also show that performing one long simulation is better than averaging results from three shortest simulations of the same length using different starting conformations.