Protein dielectric constants determined from NMR chemical shift perturbations.

Understanding the connection between protein structure and function requires a quantitative understanding of electrostatic effects. Structure-based electrostatic calculations are essential for this purpose, but their use has been limited by a long-standing discussion on which value to use for the dielectric constants (ε(eff) and ε(p)) required in Coulombic and Poisson-Boltzmann models. The currently used values for ε(eff) and ε(p) are essentially empirical parameters calibrated against thermodynamic properties that are indirect measurements of protein electric fields. We determine optimal values for ε(eff) and ε(p) by measuring protein electric fields in solution using direct detection of NMR chemical shift perturbations (CSPs). We measured CSPs in 14 proteins to get a broad and general characterization of electric fields. Coulomb's law reproduces the measured CSPs optimally with a protein dielectric constant (ε(eff)) from 3 to 13, with an optimal value across all proteins of 6.5. However, when the water-protein interface is treated with finite difference Poisson-Boltzmann calculations, the optimal protein dielectric constant (ε(p)) ranged from 2 to 5 with an optimum of 3. It is striking how similar this value is to the dielectric constant of 2-4 measured for protein powders and how different it is from the ε(p) of 6-20 used in models based on the Poisson-Boltzmann equation when calculating thermodynamic parameters. Because the value of ε(p) = 3 is obtained by analysis of NMR chemical shift perturbations instead of thermodynamic parameters such as pK(a) values, it is likely to describe only the electric field and thus represent a more general, intrinsic, and transferable ε(p) common to most folded proteins.

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

[2]  C. Domene,et al.  Calculation of dipole-shielding polarizabilities (sigma(alphabetagamma)I): the influence of uniform electric field effects on the shielding of backbone nuclei in proteins. , 2003, Journal of the American Chemical Society.

[3]  W. G. Schneider,et al.  Proton Magnetic Resonance Chemical Shift of Free (Gaseous) and Associated (Liquid) Hydride Molecules , 1958 .

[4]  Kristine Steen Jensen,et al.  Millisecond dynamics in glutaredoxin during catalytic turnover is dependent on substrate binding and absent in the resting states. , 2011, Journal of the American Chemical Society.

[5]  H Luecke,et al.  Proton transfer pathways in bacteriorhodopsin at 2.3 angstrom resolution. , 1998, Science.

[6]  Transient intermediary states with high and low folding probabilities in the apparent two-state folding equilibrium of ACBP at low pH. , 2002, Journal of molecular biology.

[7]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

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

[9]  D. LeMaster,et al.  Tetrathiolate coordination of germanium(IV) in a protein active site. , 2006, Journal of inorganic biochemistry.

[10]  J. Nielsen,et al.  Titration_DB: Storage and analysis of NMR‐monitored protein pH titration curves , 2010, Proteins.

[11]  S. Boxer,et al.  Electrostatic fields near the active site of human aldose reductase: 1. New inhibitors and vibrational stark effect measurements. , 2008, Biochemistry.

[12]  Lenwood S. Heath,et al.  H++: a server for estimating pKas and adding missing hydrogens to macromolecules , 2005, Nucleic Acids Res..

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

[14]  E. Oldfield,et al.  Chemical-shift ranges in proteins , 1992 .

[15]  Arieh Warshel,et al.  Modeling electrostatic effects in proteins. , 2006, Biochimica et biophysica acta.

[16]  J. G. Batchelor Theory of linear electric field shifts in carbon-13 nuclear magnetic resonance , 1975 .

[17]  C. Castañeda,et al.  Molecular determinants of the pKa values of Asp and Glu residues in staphylococcal nuclease , 2009, Proteins.

[18]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[19]  M. DePristo,et al.  Relation between native ensembles and experimental structures of proteins. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[20]  D. Case,et al.  A new analysis of proton chemical shifts in proteins , 1991 .

[21]  A. Warshel,et al.  Control of the redox potential of cytochrome c and microscopic dielectric effects in proteins. , 1986, Biochemistry.

[22]  J. J. Led,et al.  Probing electric fields in proteins in solution by NMR spectroscopy , 2008, Proteins.

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

[24]  B. García-Moreno E.,et al.  Charges in the hydrophobic interior of proteins , 2010, Proceedings of the National Academy of Sciences.

[25]  M. Williamson,et al.  Empirical Comparisons of Models for Chemical-Shift Calculation in Proteins , 1993 .

[26]  A. Warshel,et al.  Computer simulations of enzyme catalysis: finding out what has been optimized by evolution. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Damien Farrell,et al.  Remeasuring HEWL pKa values by NMR spectroscopy: Methods, analysis, accuracy, and implications for theoretical pKa calculations , 2011, Proteins.

[28]  B. Chait,et al.  The structure of the potassium channel: molecular basis of K+ conduction and selectivity. , 1998, Science.

[29]  A. D. Buckingham,et al.  CHEMICAL SHIFTS IN THE NUCLEAR MAGNETIC RESONANCE SPECTRA OF MOLECULES CONTAINING POLAR GROUPS , 1960 .

[30]  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.

[31]  Rebecca C. Wade,et al.  Improving the Continuum Dielectric Approach to Calculating pKas of Ionizable Groups in Proteins , 1996 .

[32]  M. D. Joshi,et al.  Dissecting electrostatic interactions in Bacillus circulans xylanase through NMR-monitored pH titrations , 2011, Journal of biomolecular NMR.

[33]  E. Oldfield,et al.  Predicting Carbon-13 Nuclear Magnetic Resonance Chemical Shielding Tensors in Zwitterionic L-Threonine and L-Tyrosine via Quantum Chemistry , 1994 .

[34]  A. Warshel,et al.  Calculations of electrostatic interactions in biological systems and in solutions , 1984, Quarterly Reviews of Biophysics.

[35]  S. Boxer,et al.  Effects of buried ionizable amino acids on the reduction potential of recombinant myoglobin. , 1989, Science.

[36]  S. Withers,et al.  Probing electrostatic interactions along the reaction pathway of a glycoside hydrolase: histidine characterization by NMR spectroscopy. , 2007, Biochemistry.

[37]  Dedreia Tull,et al.  Crystallographic observation of a covalent catalytic intermediate in a β-glycosidase , 1996, Nature Structural Biology.

[38]  C Redfield,et al.  Measurement of the individual pKa values of acidic residues of hen and turkey lysozymes by two-dimensional 1H NMR. , 1994, Biophysical journal.

[39]  D. Svergun,et al.  Metal-mediated self-assembly of a beta-sandwich protein. , 2009, Chemistry.

[40]  E. Alexov,et al.  Incorporating protein conformational flexibility into the calculation of pH-dependent protein properties. , 1997, Biophysical journal.

[41]  Barbara M. Tynan-Connolly,et al.  Redesigning protein pKa values , 2006, Protein science : a publication of the Protein Society.

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

[43]  G. Otting,et al.  Hydrophobic interactions in a cyanobacterial plastocyanin-cytochrome f complex. , 2001, Journal of the American Chemical Society.

[44]  Michael J. Holst,et al.  The adaptive multilevel finite element solution of the Poisson-Boltzmann equation on massively parallel computers , 2001, IBM J. Res. Dev..

[45]  A. Karshikoff A simple algorithm for the calculation of multiple site titration curves. , 1995, Protein engineering.

[46]  A. Warshel,et al.  Chemical shifts in proteins: a shielding trajectory analysis of the fluorine nuclear magnetic resonance spectrum of the Escherichia coli galactose binding protein using a multipole shielding polarizability-local reaction field-molecular dynamics approach , 1993 .

[47]  D. LeMaster,et al.  Backbone conformational dependence of peptide acidity. , 2009, Biophysical chemistry.

[48]  R. Pethig,et al.  Dielectric studies of protein hydration and hydration-induced flexibility. , 1985, Journal of molecular biology.

[49]  Kazumasa Sakurai,et al.  Principal component analysis of the pH-dependent conformational transitions of bovine β-lactoglobulin monitored by heteronuclear NMR , 2007, Proceedings of the National Academy of Sciences.

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

[51]  Michael W. Mahoney,et al.  A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions , 2000 .

[52]  Nathan A. Baker,et al.  On the development of protein pKa calculation algorithms , 2011, Proteins.

[53]  L. Kay,et al.  Measurement of side-chain carboxyl pK(a) values of glutamate and aspartate residues in an unfolded protein by multinuclear NMR spectroscopy. , 2002, Journal of the American Chemical Society.

[54]  P. Baker,et al.  Structural origins of pH‐dependent chemical shifts in the B1 domain of protein G , 2010, Proteins.

[55]  V. Pande,et al.  Electric Fields at the Active Site of an Enzyme: Direct Comparison of Experiment with Theory , 2006, Science.

[56]  Chresten R. Søndergaard,et al.  Improving the analysis of NMR spectra tracking pH‐induced conformational changes: Removing artefacts of the electric field on the NMR chemical shift , 2010, Proteins.

[57]  S. Firbank,et al.  Regulation of protein function: crystal packing interfaces and conformational dimerization. , 2008, Biochemistry.

[58]  J. Warwicker,et al.  Simplified methods for pKa and acid pH‐dependent stability estimation in proteins: Removing dielectric and counterion boundaries , 2008, Protein science : a publication of the Protein Society.