Dissemin is shutting down on January 1st, 2025

Published in

Royal Society of Chemistry, Physical Chemistry Chemical Physics, 4(17), p. 2859-2868, 2015

DOI: 10.1039/c4cp04538f

Links

Tools

Export citation

Search in Google Scholar

Are thermodynamic cycles necessary for continuum solvent calculation of pKas and reduction potentials?

Journal article published in 2014 by Junming Ho ORCID
This paper is available in a repository.
This paper is available in a repository.

Full text: Download

Green circle
Preprint: archiving allowed
Orange circle
Postprint: archiving restricted
Red circle
Published version: archiving forbidden
Data provided by SHERPA/RoMEO

Abstract

Continuum solvent calculations of pKas and reduction potentials usually entail the use of a thermodynamic cycle to express the reaction free energy in terms of gas phase energies and free energies of solvation. In this work, we present a systematic study comparing the solution phase free energy changes obtained in this manner with that directly computed within the SMD solvation model against a large test set of 117 pKas and 42 reduction potentials in water and DMSO. The inclusion of vibrational contributions in the free energy of solvation has negligible impact on the accuracy of thermodynamic cycle predictions of pKas and reduction potentials. Additionally, when gas phase energies in the thermodynamic cycle are computed at more accurate levels of theory, very similar results (mean unsigned difference of 0.5 kcal mol-1) can be achieved when the high level computations (MP2/GTMP2Large and G3(MP2)-RAD(+)) are directly carried out within the continuum model. Increasing the accuracy of the electronic structure theory may or may not improve the agreement with experiment suggesting that the error is largely in the solvation model. For amino acids where their gas and solution phase species exist as different tautomers, the direct approach provided a significant improvement in calculated pKas. These results demonstrate that direct calculation of solution phase pKas and reduction potentials within the SMD model provides a general and reliable approximation to corresponding thermodynamic cycle based protocols, and is recommended for systems where solvation induced changes in geometry are significant. Further studies are necessary to ascertain whether the results are generalisable to other continuum solvation models