Journal
JOURNAL OF PHYSICAL CHEMISTRY B
Volume 115, Issue 19, Pages 6011-6022Publisher
AMER CHEMICAL SOC
DOI: 10.1021/jp111271c
Keywords
-
Categories
Funding
- John von Neumann-Institut fur Computing (NIC)
- Juelich Supercomputing Centre (JSC)
- Forschungszentrum Juelich GmbH, Germany [ID HLZ18]
- Deutsche Forschungsgemeinschaft (DFG)-German Research Foundation [FE 1156/2-1]
- REA-Research Executive Agency [247500, FP7-PEOPLE-2009-IRSES]
Ask authors/readers for more resources
We report a method to dramatically improve the accuracy of hydration free energies (HFE) calculated by the 1D and 3D reference interaction site models (RISM) of molecular integral equation theory. It is shown that the errors in HFEs calculated by RISM approaches using the Gaussian fluctuations (GF) free energy functional are not random, but can be decomposed into linear combination of contributions from different structural elements of molecules (number of double bonds, number of OH groups, etc.). Therefore, by combining RISM/GF with cheminformatics, one can develop an accurate method for HFE prediction. We call this approach the structural description correction model (SDC) (Ratkova et al. J. Phys. Chem. B 2010, 114, 12068). In this work, we investigated the prediction quality of the SDC model combined with 1D and 3D RISM approaches. In parallel, we analyzed the computational performance of these two methods. The SD C model parameters were obtained by fitting against a training set of 53 simple organic molecules. To demonstrate that the values of these parameters were transferable between different classes of molecules, the models were tested against 98 more complex molecules (including 38 polyfragment compounds). The results show that the 3D RISM/SDC model predicts the HFEs with very good accuracy (RMSE of 0.47 kcal/mol), while the ID RISM approach provides only moderate accuracy (RMSE of 1.96 kcal/mol). However, a single ID RISM/SDC calculation takes only a few seconds on a PC, whereas a single 3D RISM/SDC HFE calculation is approximately 100 times more computationally expensive. Therefore, we suggest that one should use the 1D RISM/SDC model for large-scale high-throughput screening of molecular hydration properties, while further refinement of these properties for selected compounds should be carried out with the more computationally expensive but more accurate 3D RISM/SDC model.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available