文献类型：期刊文献

英文题名：Comparison of two simulation methods to compute solvation free energies and partition coefficients

作者：Yang, Li[1,2,3];Ahmed, Alauddin[1];Sandler, Stanley I.[1]

机构：[1]Univ Delaware, Ctr Mol & Engn Thermodynam, Dept Chem & Biomol Engn, Newark, DE 19716 USA;[2]E China Univ Sci & Technol, State Key Lab Chem Engn, Shanghai 200237, Peoples R China;[3]E China Univ Sci & Technol, Dept Chem, Shanghai 200237, Peoples R China

年份：2013

卷号：34

期号：4

起止页码：284

外文期刊名：JOURNAL OF COMPUTATIONAL CHEMISTRY

收录：;EI(收录号：20130415935920);WOS:【SCI-EXPANDED(收录号:WOS:000313723500003)】；

基金：Contract/grant sponsor: National Natural Science Foundation of China; Contract/grant numbers: 21136004; Contract/grant sponsor: National Basic Research Program of China; Contract/grant numbers: 2009CB219902; Contract/grant sponsor: National Science Foundation; Contract/grant number: GOALI-0853685; Contract/grant sponsor: Catalysis Center for Energy Innovation, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, and Office of Basic Energy Sciences; Contract/grant number: DE-SC0001004.

语种：英文

外文关键词：thermodynamic integration; expanded ensemble; 1-octanol; water coefficient; solvation energy; Gibbs energy

摘要：The thermodynamic integration (TI) and expanded ensemble (EE) methods are used here to calculate the hydration free energy in water, the solvation free energy in 1-octanol, and the octanol-water partition coefficient for a six compounds of varying functionality using the optimized potentials for liquid simulations (OPLS) all-atom (AA) force field parameters and atomic charges. Both methods use the molecular dynamics algorithm as a primary component of the simulation protocol, and both have found wide applications in fields such as the calculation of activity coefficients, phase behavior, and partition coefficients. Both methods result in solvation free energies and 1-octanol/water partition coefficients with average absolute deviations (AAD) from experimental data to within 4 kJ/mol and 0.5 log units, respectively. Here, we find that in simulations the OPLS-AA force field parameters (with fixed charges) can reproduce solvation free energies of solutes in 1-octanol with AAD of about half that for the solute hydration free energies using a extended simple point charge (SPC/E) model of water. The computational efficiency of the two simulation methods are compared based on the time (in nanoseconds) required to obtain similar standard deviations in the solvation free energies and 1-octanol/water partition coefficients. By this analysis, the EE method is found to be a factor of nine more efficient than the TI algorithm. For both methods, solvation free energy calculations in 1-octanol consume roughly an order of magnitude more CPU hours than the hydration free energy calculations. (C) 2012 Wiley Periodicals, Inc.