How Does the Methodology of 3D Structure Preparation Influence the Quality of pK(a) Prediction?

Investor logo
Investor logo
Investor logo
Investor logo

Warning

This publication doesn't include Faculty of Education. It includes Central European Institute of Technology. Official publication website can be found on muni.cz.
Authors

GEIDL Stanislav SVOBODOVÁ VAŘEKOVÁ Radka BENDOVÁ Veronika PETRUSEK Lukáš IONESCU Crina-Maria JURKA Zdeněk ABAGYAN Ruben A. KOČA Jaroslav

Year of publication 2015
Type Article in Periodical
Magazine / Source Journal of Chemical Information and Modeling
MU Faculty or unit

Central European Institute of Technology

Citation
Web http://pubs.acs.org/doi/pdf/10.1021/ci500758w
Doi http://dx.doi.org/10.1021/ci500758w
Field Biochemistry
Keywords ELECTRONEGATIVITY EQUALIZATION METHOD; ATOMIC CHARGE SCHEMES; CARBOXYLIC-ACIDS; PHYSICAL-PROPERTIES; PHENOLS; VALUES; VALIDATION; QSAR; DESCRIPTORS; ALGORITHM
Attached files
Description The acid dissociation constant is an important molecular property, and it can be successfully predicted by Quantitative Structure-Property Relationship (QSPR) models, even for in silico designed molecules. We analyzed how the methodology of in silico 3D structure preparation influences the quality of QSPR models. Specifically, we evaluated and compared QSPR models based on six different 3D structure sources (DTP NCI, Pubchem, Balloon, Frog2, OpenBabel, and RDKit) combined with four different types of optimization. These analyses were performed for three classes of molecules (phenols, carboxylic acids, anilines), and the QSPR model descriptors were quantum mechanical (QM) and empirical partial atomic charges. Specifically, we developed 516 QSPR models and afterward systematically analyzed the influence of the 3D structure source and other factors on their quality. Our results confirmed that QSPR models based on partial atomic charges are able to predict pK(a) with high accuracy. We also confirmed that ab initio and semiempirical QM charges provide very accurate QSPR models and using empirical charges based on electronegativity equalization is also acceptable, as well as advantageous, because their calculation is very fast. On the other hand, Gasteiger-Marsili empirical charges are not applicable for pKa prediction. We later found that QSPR models for some classes of molecules (carboxylic acids) are less accurate. In this context, we compared the influence of different 3D structure sources. We found that an appropriate selection of 3D structure source and optimization method is essential for the successful QSPR modeling of pK(a). Specifically, the 3D structures from the DTP NCI and Pubchem databases performed the best, as they provided very accurate QSPR models for all the tested molecular classes and charge calculation approaches, and they do not require optimization. Also, Frog2 performed very well. Other 3D structure sources can also be used but are not so robust, and an unfortunate combination of molecular class and charge calculation approach can produce weak QSPR models. Additionally, these 3D structures generally need optimization in order to produce good quality QSPR models.
Related projects:

You are running an old browser version. We recommend updating your browser to its latest version.