Premium
Improved precision and efficiency of free energy calculations for small systems using λ‐scaled atomic masses and separating conformational and transformational sampling
Author(s) -
Carlsson Peter,
Nilsson Lennart
Publication year - 2003
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.10221
Subject(s) - scaling , molecular dynamics , chemistry , chemical physics , umbrella sampling , energy (signal processing) , molecule , thermodynamic integration , molecular physics , computational chemistry , atomic physics , physics , quantum mechanics , mathematics , geometry , organic chemistry
We present results showing the importance of appropriate treatment of atomic masses in molecular dynamics (MD)‐based single topology free‐energy perturbations (FEPs) on small molecule systems. The reversibility of gas phase simulations is significantly improved by scaling the atomic mass of mutated atoms with the λ variable normally used for the scaling of energy terms. Because this effect is less pronounced for solvated systems, it will not cancel in estimates of the relative hydration free energy difference. The advantage of mass scaling is demonstrated by a null mutation of ethane to ethane and the calculation of the relative hydration free energy difference between ethane and n ‐propane. Furthermore, it is found that the simulation time necessary for converged MD/FEPs is prohibitively large for relative hydration free energy calculations on cyclic alkanes. Therefore, we explore an alternative free energy pathway including strongly constrained conformations to improve convergence in FEP simulations of flexible molecules. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1383–1389, 2003