In this study, we present an accurate protocol for the fast prediction of pKa's of carboxylic acids based on the linear relationship between computed atomic charges of the anionic form of the carboxylate fragment and their experimental pKa values. Five charge descriptors, three charge models, three solvent models, gas phase calculations and several DFT methods (combination of eight DFT functionals and fifteen basis sets) were tested. Among those, the best combination to reproduce experimental pKa's is to compute NPA atomic charge using the SMD model at the M06L/6-311G(d,p) level of theory and selecting the maximum atomic charge on carboxylic oxygen atoms (R^2 = 0.955). The applicability of the suggested protocol and its stability along geometrical changes are verified by molecular dynamics simulations performed for a set of aspartate, glutamate and alanine peptides. By reporting the calculated atomic charge of the carboxylate form into the linear relationship derived in this work, it will be possible to estimate accurately the amino acid's pKa's in protein environment.
File list (2)download file view on ChemRxiv main.pdf (733.11 KiB) download file view on ChemRxiv si-main.pdf (364.13 KiB)
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.