Prediction of Absolute Hydroxyl pKa Values for 3-Hydroxypyridin-4-ones

TY - JOUR

T1 - Prediction of Absolute Hydroxyl pKa Values for 3-Hydroxypyridin-4-ones

AU - Chen, Yu-Lin

AU - Doltsinis, Nikos L.

AU - Hider, Robert C.

AU - Barlow, David

PY - 2012/10/18

Y1 - 2012/10/18

N2 - pK(a) values have been calculated for a series of 3-hydroxypyridin-4-one (HPO) chelators in aqueous solution using coordination constrained ab initio molecular dynamics (AIMD) in combination with thermodynamic integration. This dynamics-based methodology in which the solvent is treated explicitly at the ab initio level has been compared with more commonly used simple, static, approaches. Comparison with experimental numbers has confirmed that the AIMD-based approach predicts the correct trend in the pK(a) values and produces the lowest average error (similar to 0.3 pK(a) units). The corresponding pK(a) predictions made via static quantum mechanical calculations overestimate the pK(a) values by 0.3-7 pK(a) units, with the extent of error dependent on the choice of thermodynamic cycle employed. The use of simple quantitative structure property relationship methods gives prediction errors of 0.3-1 pK(a) units, with some values overestimated and some underestimated. Beyond merely calculating pK(a) values, the AIMD simulations provide valuable additional insight into the atomistic details of the proton transfer mechanism and the solvation structure and dynamics at all stages of the reaction. For all HPOs studied, it is seen that proton transfer takes place along a chain of three H2O molecules, although direct hydrogen bonds are seen to form transiently. Analysis of the solvation structure before and after the proton transfer event using radial pair distribution functions and integrated number densities suggests that the trends in the pK(a) values correlate with the strength of the hydrogen bond and the average number of solvent molecules in the vicinity of the donor oxygen.

AB - pK(a) values have been calculated for a series of 3-hydroxypyridin-4-one (HPO) chelators in aqueous solution using coordination constrained ab initio molecular dynamics (AIMD) in combination with thermodynamic integration. This dynamics-based methodology in which the solvent is treated explicitly at the ab initio level has been compared with more commonly used simple, static, approaches. Comparison with experimental numbers has confirmed that the AIMD-based approach predicts the correct trend in the pK(a) values and produces the lowest average error (similar to 0.3 pK(a) units). The corresponding pK(a) predictions made via static quantum mechanical calculations overestimate the pK(a) values by 0.3-7 pK(a) units, with the extent of error dependent on the choice of thermodynamic cycle employed. The use of simple quantitative structure property relationship methods gives prediction errors of 0.3-1 pK(a) units, with some values overestimated and some underestimated. Beyond merely calculating pK(a) values, the AIMD simulations provide valuable additional insight into the atomistic details of the proton transfer mechanism and the solvation structure and dynamics at all stages of the reaction. For all HPOs studied, it is seen that proton transfer takes place along a chain of three H2O molecules, although direct hydrogen bonds are seen to form transiently. Analysis of the solvation structure before and after the proton transfer event using radial pair distribution functions and integrated number densities suggests that the trends in the pK(a) values correlate with the strength of the hydrogen bond and the average number of solvent molecules in the vicinity of the donor oxygen.

U2 - 10.1021/jz301061m

DO - 10.1021/jz301061m

M3 - Article

SN - 1948-7185

VL - 3

SP - 2980

EP - 2985

JO - Journal of physical chemistry letters

JF - Journal of physical chemistry letters

IS - 20

ER -