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 pK(a) calculations using homology modeled three-dimensional (3D) structures and analyze the possibility of using the calculated pK(a) 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 pK(a) 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 pK(a) 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 pK(a) calculations. In addition, it was found that some highly conserved titratable groups either exhibit large pK(a) fluctuations among the models or have pK(a) values shifted by several pH units with respect to the pK(a) calculated for the X-ray structure. We demonstrated that such case usually indicates structural errors associated with the model. Thus, we argue that pK(a) calculations can be used for assessing the quality of the 3D models by monitoring fluctuations of the pK(a) values for highly conserved titratable residues within large sets of homologous proteins.

[1]  Philip E. Bourne,et al.  The RCSB PDB information portal for structural genomics , 2005, Nucleic Acids Res..

[2]  Lei Xie,et al.  Using multiple structure alignments, fast model building, and energetic analysis in fold recognition and homology modeling , 2003, Proteins.

[3]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[4]  G Klebe,et al.  Improving macromolecular electrostatics calculations. , 1999, Protein engineering.

[5]  M C Peitsch,et al.  ProMod and Swiss-Model: Internet-based tools for automated comparative protein modelling. , 1996, Biochemical Society transactions.

[6]  K. Sharp,et al.  Linkage of thioredoxin stability to titration of ionizable groups with perturbed pKa. , 1991, Biochemistry.

[7]  A. Warshel,et al.  Electrostatic effects in macromolecules: fundamental concepts and practical modeling. , 1998, Current opinion in structural biology.

[8]  A. Warshel,et al.  Free energy of charges in solvated proteins: microscopic calculations using a reversible charging process. , 1986, Biochemistry.

[9]  Emil Alexov,et al.  Calculating proton uptake/release and binding free energy taking into account ionization and conformation changes induced by protein–inhibitor association: Application to plasmepsin, cathepsin D and endothiapepsin–pepstatin complexes , 2004, Proteins.

[10]  Thomas L. Madden,et al.  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. , 1997, Nucleic acids research.

[11]  A. Warshel.  Computer simulations of enzyme catalysis: methods, progress, and insights. , 2003, Annual review of biophysics and biomolecular structure.

[12]  M. Ondrechen,et al.  THEMATICS: A simple computational predictor of enzyme function from structure , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[13]  E. Alexov,et al.  Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins. , 2002, Biophysical journal.

[14]  Ying Wei,et al.  Active Site Prediction for Comparative Model Structures with Thematics , 2005, J. Bioinform. Comput. Biol..

[15]  J. Mccammon,et al.  On the evaluation and optimization of protein X‐ray structures for pKa calculations , 2003, Protein science : a publication of the Protein Society.

[16]  B. Tidor,et al.  Do salt bridges stabilize proteins? A continuum electrostatic analysis , 1994, Protein science : a publication of the Protein Society.

[17]  A. Warshel.  Calculations of enzymatic reactions: calculations of pKa, proton transfer reactions, and general acid catalysis reactions in enzymes. , 1981, Biochemistry.

[18]  E. Alexov,et al.  Comparative study of the stability of poplar plastocyanin isoforms. , 2005, Biochimica et biophysica acta.

[19]  G. Labesse,et al.  Comparative modelling and immunochemical reactivity of Escherichia coli UMP kinase. , 2002, Biochemical and biophysical research communications.

[20]  M. Karplus,et al.  pKa's of ionizable groups in proteins: atomic detail from a continuum electrostatic model. , 1990, Biochemistry.

[21]  T. Blundell,et al.  Comparative protein modelling by satisfaction of spatial restraints. , 1993, Journal of molecular biology.

[22]  Ihsan A. Shehadi,et al.  Future directions in protein function prediction , 2002, Molecular Biology Reports.

[23]  Charles L. Brooks,et al.  CHARGE SCREENING AND THE DIELECTRIC CONSTANT OF PROTEINS : INSIGHTS FROM MOLECULAR DYNAMICS , 1996 .

[24]  E. Alexov,et al.  Calculated protein and proton motions coupled to electron transfer: electron transfer from QA- to QB in bacterial photosynthetic reaction centers. , 1999, Biochemistry.

[25]  Ron D. Appel,et al.  ExPASy: the proteomics server for in-depth protein knowledge and analysis , 2003, Nucleic Acids Res..

[26]  B. Honig,et al.  On the role of structural information in remote homology detection and sequence alignment: new methods using hybrid sequence profiles. , 2003, Journal of molecular biology.

[27]  M. Gilson,et al.  Prediction of pH-dependent properties of proteins. , 1994, Journal of molecular biology.

[28]  C. Brooks,et al.  Constant‐pH molecular dynamics using continuous titration coordinates , 2004, Proteins.

[29]  G. Vriend,et al.  Optimizing the hydrogen‐bond network in Poisson–Boltzmann equation‐based pKa calculations , 2001, Proteins.

[30]  Emil Alexov,et al.  Electrostatic properties of protein-protein complexes. , 2006, Biophysical journal.

[31]  Arieh Warshel,et al.  Consistent Calculations of pKa's of Ionizable Residues in Proteins: Semi-microscopic and Microscopic Approaches , 1997 .

[32]  K. Sharp,et al.  On the calculation of pKas in proteins , 1993, Proteins.

[33]  C. Marco-Marín,et al.  First-time crystallization and preliminary X-ray crystallographic analysis of a bacterial-archaeal type UMP kinase, a key enzyme in microbial pyrimidine biosynthesis. , 2005, Biochimica et biophysica acta.

[34]  M. Gilson,et al.  Computing ionization states of proteins with a detailed charge model , 1996, J. Comput. Chem..

[35]  B. Honig,et al.  Classical electrostatics in biology and chemistry. , 1995, Science.

[36]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[37]  J. Mccammon,et al.  Calculating pKa values in enzyme active sites , 2003, Protein science : a publication of the Protein Society.

[38]  B. Tidor,et al.  Surface salt bridges, double-mutant cycles, and protein stability: an experimental and computational analysis of the interaction of the Asp 23 side chain with the N-terminus of the N-terminal domain of the ribosomal protein l9. , 2003, Biochemistry.

[39]  Ronald J. Williams,et al.  Statistical criteria for the identification of protein active sites using theoretical microscopic titration curves , 2005, Proteins.

[40]  Emil Alexov,et al.  Role of the protein side‐chain fluctuations on the strength of pair‐wise electrostatic interactions: Comparing experimental with computed pKas , 2002, Proteins.

[41]  Honggao Yan,et al.  Nucleoside monophosphate kinases: structure, mechanism, and substrate specificity. , 1999, Advances in enzymology and related areas of molecular biology.

[42]  Yifan Song,et al.  Calculated proton uptake on anaerobic reduction of cytochrome C oxidase: is the reaction electroneutral? , 2006, Biochemistry.