Assessment of Methodology and Chemical Group Dependences in the Calculation of the pKa for Several Chemical Groups
Assessment of Methodology and Chemical Group Dependences in the Calculation of the pKa for Several Chemical Groups
复制标题
评估几种化学基团 pKa 计算中的方法和化学基团依赖性
DOI:
10.1021/acs.jctc.7b00587
复制
发表时间:
2017
影响因子:
5.5
通讯作者:
Morihashi Kenji
中科院分区:
文献类型:
--
作者:
Matsui Toru;Shigeta Yasuteru;Morihashi Kenji
We have investigated the dependencies of various computational methods in the calculation of acid dissociation constants (pKavalues) of certain chemical groups found in protonatable amino acids based on our previous scheme [Matsui; Phys. Chem. Chem. Phys. 2012, 14, 4181−4187]. By changing the quantum chemical (QC) method (Hartree–Fock (HF) and perturbation theory, and composite methods, or exchange–correlation functionals in density functional theory (DFT)), basis sets, solvation models, and the cavities used in the solvent models, we have exhaustively tested about 2,200 combinations to find the best combination for pKaestimation among them. Of the tested parameters, the choice of the basis set and cavity is the most crucial to reproduce experimental values compared to other factors. Concerning the basis set, the inclusion of diffuse functions is quite important for carboxyl, thiol, and phenol groups judging from the mean absolute errors (MAEs) measured from the experimental values. Of the cavity models, between the Pauling, Klamt, and the universal force field (UFF) definitions, the UFF defined cavity is the best choice, resulting in the smallest MAEs. Concerning the QC methods, hybrid DFTs and range-separated DFTs always provide better results than pure DFTs and HF. As a result, we found that LC-ϖPBE/6-31+G(d) with PCM-SMD/UFF provides the best pKaestimation with a MAE within 0.15 pKaunits.