ASSESSING THE QUALITY OF THE HOMOLOGY-MODELED 3D STRUCTURES FROM ELECTROSTATIC STANDPOINT: TEST ON BACTERIAL NUCLEOSIDE MONOPHOSPHATE KINASE FAMILIES
In this study, we address the issue of performing meaningful pKa calculations using homology modeled three-dimensional (3D) structures and analyze the possibility of using the calculated pKa values to detect structural defects in the models. For this purpose, the 3D structure of each member of five large protein families of a bacterial nucleoside monophosphate kinases (NMPK) have been modeled by means of homology-based approach. Further, we performed pKa calculations for the each model and for the template X-ray structures. Each bacterial NMPK family used in the study comprised on average 100 members providing a pool of sequences and 3D models large enough for reliable statistical analysis. It was shown that pKa values of titratable groups, which are highly conserved within a family, tend to be conserved among the models too. We demonstrated that homology modeled structures with sequence identity larger than 35% and gap percentile smaller than 10% can be used for meaningful pKa calculations. In addition, it was found that some highly conserved titratable groups either exhibit large pKa fluctuations among the models or have pKa values shifted by several pH units with respect to the pKa calculated for the X-ray structure. We demonstrated that such case usually indicates structural errors associated with the model. Thus, we argue that pKa calculations can be used for assessing the quality of the 3D models by monitoring fluctuations of the pKa values for highly conserved titratable residues within large sets of homologous proteins.