Performance of Density Functional Models to Reproduce Observed 13Ca Chemical Shifts of Proteins in Solution

VILA, J.A., BALDONI, H.A. AND SCHERAGA, H.A.

JOURNAL OF COMPUTATIONAL CHEMISTRY

JOHN WILEY & SONS INC

Año: 2009 vol. 30 p. 884 - 892

0192-8651

The purpose of this work is to test several density functional models (namely, OPBE, 3LYP, OPW91, BPW91, OB98, BPBE, B971, OLYP, PBE1PBE, and B3LYP) to determine their accuracy and speed for computing 13Ca chemical shifts in proteins. The test is applied to 10 NMR-derived conformations of the 76-residue a/b protein ubiquitin (protein data bank id 1D3Z). With each functional, the 13Ca shielding was computed for 760 amino acid residues by using a combination of approaches that includes, but is not limited to, treating each amino acid X in the sequence as a terminally blocked tripeptide with the sequence Ac-GXG-NMe in the conformation of the regularized experimental protein structure. As computation of the 13Ca chemical shifts, not their shielding, is the main goal of this work, a computation of the 13Ca shielding of the reference, namely, tetramethylsilane, is investigated here and an effective and a computed tetramethylsilane shielding value for each of the functionals is provided. Despite observed small differences among all functionals tested, the results indicate that four of them, namely, OPBE, OPW91, OB98, and OLYP, provide the most accurate functionals with which to reproduce observed 13Ca chemical shifts of proteins in solution, and are among the faster ones. This study also provides evidence for the applicability of these functionals to proteins of any size or class, and for the validation of our previous results and conclusions, obtained from calculations with the slower B3LYP functional.