Abstract
The “embedded cluster reference interaction site model” (EC-RISM) approach combines statistical-mechanical integral equation theory and quantum-chemical calculations for predicting thermodynamic data for chemical reactions in solution. The electronic structure of the solute is determined self-consistently with the structure of the solvent that is described by 3D RISM integral equation theory. The continuous solvent-site distribution is mapped onto a set of discrete background charges (“embedded cluster”) that represent an additional contribution to the molecular Hamiltonian. The EC-RISM analysis of the SAMPL2 challenge set of tautomers proceeds in three stages. Firstly, the group of compounds for which quantitative experimental free energy data was provided was taken to determine appropriate levels of quantum-chemical theory for geometry optimization and free energy prediction. Secondly, the resulting workflow was applied to the full set, allowing for chemical interpretations of the results. Thirdly, disclosure of experimental data for parts of the compounds facilitated a detailed analysis of methodical issues and suggestions for future improvements of the model. Without specifically adjusting parameters, the EC-RISM model yields the smallest value of the root mean square error for the first set (0.6 kcal mol−1) as well as for the full set of quantitative reaction data (2.0 kcal mol−1) among the SAMPL2 participants.
Similar content being viewed by others
References
Kloss T, Heil J, Kast SM (2008) J Phys Chem B 112:4337–4343
Beglov D, Roux B (1997) J Phys Chem 101:7821–7826
Kovalenko A, Hirata F (1998) Chem Phys Lett 290:237–244
Taylor PJ (2009) Tautomeric teasers: SAMPL (to be replaced by overview paper)
Kast SM, Kloss T (2008) J Chem Phys 129:236101
Perkyns J, Pettitt BM (1992) Chem Phys Lett 190:626–630
Perkyns J, Pettitt BM (1992) J Chem Phys 97:7656–7666
Frisch MJ et al (2004) Gaussian 03, Rev D.02/E.01. Gaussian, Inc., Wallingford
Berendsen HJC, Grigera JR, Straatsma TP (1987) J Phys Chem 91:6269
Maw S, Sato H, Ten-no S, Hirata F (1997) Chem Phys Lett 276:20–25
Sato H, Hirata F (1999) J Chem Phys 111:8545–8555
Kovalenko A, Ten-no S, Hirata F (1999) J Comput Chem 20:928–936
Talman JD (1978) J Comput Phys 29:35–48
Rossky PJ, Friedman HL (1980) J Chem Phys 72:5694–5700
Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA (2004) J Comput Chem 25:1157–1174
Chirlian LE, Francl MM (1987) J Comput Chem 8:894–905
Acknowledgments
We thank the Deutsche Forschungsgemeinschaft, the Fonds der Chemischen Industrie, and the Adolf-Messer-Stiftung for financial support.
Author information
Authors and Affiliations
Corresponding author
Additional information
An erratum to this article can be found at http://dx.doi.org/10.1007/s10822-010-9360-6
Rights and permissions
About this article
Cite this article
Kast, S.M., Heil, J., Güssregen, S. et al. Prediction of tautomer ratios by embedded-cluster integral equation theory. J Comput Aided Mol Des 24, 343–353 (2010). https://doi.org/10.1007/s10822-010-9340-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10822-010-9340-x