A buried lysine that titrates with a normal pKa: Role of conformational flexibility at the protein–water interface as a determinant of pKavalues

Previously we reported that Lys, Asp, and Glu residues at positions 66 and 92 in staphylococcal nuclease (SNase) titrate with pKa values shifted by up to 5 pKa units in the direction that promotes the neutral state. In contrast, the internal Lys‐38 in SNase titrates with a normal pKa. The crystal structure of the L38K variant shows that the side chain of Lys‐38 is buried. The ionizable moiety is ∼7 Å from solvent and ion paired with Glu‐122. This suggests that the pKa value of Lys‐38 is normal because the energetic penalty for dehydration is offset by a favorable Coulomb interaction. However, the pKa of Lys‐38 was also normal when Glu‐122 was replaced with Gln or with Ala. Continuum electrostatics calculations were unable to reproduce the pKa of Lys‐38 unless the protein was treated with an artificially high dielectric constant, consistent with structural reorganization being responsible for the normal pKa value of Lys‐38. This reorganization must be local because circular dichroism and NMR spectroscopy indicate that the L38K protein is native‐like under all conditions studied. In molecular dynamics simulations, the ion pair between Lys‐38 and Glu‐122 is unstable. The simulations show that a minor rearrangement of a loop is sufficient to allow penetration of water to the amino moiety of Lys‐38. This illustrates both the important roles of local flexibility and water penetration as determinants of pKa values of ionizable groups buried near the protein–water interface, and the challenges faced by structure‐based pKa calculations in reproducing these effects.

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

[2]  V. Hilser,et al.  Molecular mechanisms of pH‐driven conformational transitions of proteins: Insights from continuum electrostatics calculations of acid unfolding , 2006, Proteins.

[3]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[4]  Jinrang Kim,et al.  Are acidic and basic groups in buried proteins predicted to be ionized? , 2005, Journal of molecular biology.

[5]  C. Anfinsen,et al.  An improved method for the purification of staphylococcal nuclease. , 1967, The Journal of biological chemistry.

[6]  Steven T. Whitten and,et al.  pH dependence of stability of staphylococcal nuclease: evidence of substantial electrostatic interactions in the denatured state. , 2000 .

[7]  R. Raines,et al.  Microscopic pKa values of Escherichia coli thioredoxin. , 1997, Biochemistry.

[8]  E. Lattman,et al.  High apparent dielectric constants in the interior of a protein reflect water penetration. , 2000, Biophysical journal.

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

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

[11]  Jan H. Jensen,et al.  Very fast empirical prediction and rationalization of protein pKa values , 2005, Proteins.

[12]  Arieh Warshel,et al.  Microscopic simulations of macroscopic dielectric constants of solvated proteins , 1991 .

[13]  D. Shortle,et al.  Residual structure in large fragments of staphylococcal nuclease: effects of amino acid substitutions. , 1989, Biochemistry.

[14]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[15]  Eaton E Lattman,et al.  Experimental pK(a) values of buried residues: analysis with continuum methods and role of water penetration. , 2002, Biophysical journal.

[16]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[17]  W E Stites,et al.  In a staphylococcal nuclease mutant the side-chain of a lysine replacing valine 66 is fully buried in the hydrophobic core. , 1991, Journal of molecular biology.

[18]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[19]  L. R. Scott,et al.  Electrostatics and diffusion of molecules in solution: simulations with the University of Houston Brownian dynamics program , 1995 .

[20]  Michael K. Gilson,et al.  The determinants of pK(a)s in proteins , 1996 .

[21]  M. Karplus,et al.  Simulation of activation free energies in molecular systems , 1996 .

[22]  T. Harris,et al.  Effects of mutations of the active site arginine residues in 4-oxalocrotonate tautomerase on the pKa values of active site residues and on the pH dependence of catalysis. , 1999, Biochemistry.

[23]  W E Stites,et al.  Increasing the thermostability of staphylococcal nuclease: implications for the origin of protein thermostability. , 2000, Journal of molecular biology.

[24]  W E Stites,et al.  Contributions of the large hydrophobic amino acids to the stability of staphylococcal nuclease. , 1990, Biochemistry.

[25]  Z. Zolnai,et al.  Solution structures of staphylococcal nuclease from multidimensional, multinuclear NMR: Nuclease-H124L and its ternary complex with Ca2+ and thymidine-3′,5′-bisphosphate , 1997, Journal of biomolecular NMR.

[26]  C. Pace,et al.  Buried, charged, non-ion-paired aspartic acid 76 contributes favorably to the conformational stability of ribonuclease T1. , 1999, Biochemistry.

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

[28]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[29]  D. Pérahia,et al.  Polar fluctuations in proteins: molecular-dynamic studies of cytochrome c in aqueous solution. , 1996, Faraday discussions.

[30]  Jinrang Kim,et al.  Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculations. , 2006, Biochimica et biophysica acta.

[31]  A. Goldman,et al.  Buried charged surface in proteins. , 2000, Structure.

[32]  A. Warshel,et al.  The effect of protein relaxation on charge-charge interactions and dielectric constants of proteins. , 1998, Biophysical journal.

[33]  H. Kamikubo,et al.  Elucidation of information encoded in tryptophan 140 of staphylococcal nuclease , 2004, Proteins.

[34]  C. Tanford,et al.  Theory of Protein Titration Curves. I. General Equations for Impenetrable Spheres , 1957 .

[35]  Nathan A. Baker,et al.  PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations , 2004, Nucleic Acids Res..

[36]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

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

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

[39]  J. Deisenhofer,et al.  Structure of the protein subunits in the photosynthetic reaction centre of Rhodopseudomonas viridis at 3Å resolution , 1985, Nature.

[40]  Y. Yamagata,et al.  Completely buried, non-ion-paired glutamic acid contributes favorably to the conformational stability of pyrrolidone carboxyl peptidases from hyperthermophiles. , 2006, Biochemistry.

[41]  Youxing Jiang,et al.  The principle of gating charge movement in a voltage-dependent K+ channel , 2003, Nature.

[42]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

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

[44]  J. Zou,et al.  Improved methods for building protein models in electron density maps and the location of errors in these models. , 1991, Acta crystallographica. Section A, Foundations of crystallography.

[45]  M. Stahley,et al.  High apparent dielectric constant inside a protein reflects structural reorganization coupled to the ionization of an internal Asp. , 2007, Biophysical journal.

[46]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[47]  D. Case,et al.  Use of chemical shifts in macromolecular structure determination. , 2002, Methods in enzymology.

[48]  M Karplus,et al.  Improving the accuracy of protein pKa calculations: Conformational averaging versus the average structure , 1998, Proteins.

[49]  E. Lattman,et al.  X-ray and thermodynamic studies of staphylococcal nuclease variants I92E and I92K: insights into polarity of the protein interior. , 2004, Journal of molecular biology.

[50]  T. Harris,et al.  Structural Basis of Perturbed pKa Values of Catalytic Groups in Enzyme Active Sites , 2002, IUBMB life.

[51]  Structural and Dynamic Properties of the CAGQW Peptide in Water: A Molecular Dynamics Simulation Study Using Different Force Fields , 2004 .

[52]  Kevin Cowtan,et al.  research papers Acta Crystallographica Section D Biological , 2005 .

[53]  A. Garcia,et al.  Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins. , 2007, Biophysical journal.

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

[55]  S C Harvey,et al.  Dielectric relaxation spectra of water adsorbed on lysozyme. , 1972, The Journal of physical chemistry.

[56]  Yuko Okamoto,et al.  Secondary-structure preferences of force fields for proteins evaluated by generalized-ensemble simulations , 2004 .

[57]  Z. Otwinowski,et al.  Processing of X-ray diffraction data collected in oscillation mode. , 1997, Methods in enzymology.

[58]  T. Hynes,et al.  The crystal structure of staphylococcal nuclease refined at 1.7 Å resolution , 1991, Proteins.

[59]  P. Kollman,et al.  Continuum Solvent Studies of the Stability of DNA, RNA, and Phosphoramidate−DNA Helices , 1998 .

[60]  C. Chothia The nature of the accessible and buried surfaces in proteins. , 1976, Journal of molecular biology.

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

[62]  A. Warshel,et al.  Using a charging coordinate in studies of ionization induced partial unfolding. , 2006, The journal of physical chemistry. B.

[63]  S. Whitten,et al.  pH dependence of stability of staphylococcal nuclease: evidence of substantial electrostatic interactions in the denatured state. , 2000, Biochemistry.

[64]  M. Gilson,et al.  The determinants of pKas in proteins. , 1996, Biochemistry.

[65]  Alexander D. MacKerell Empirical force fields for biological macromolecules: Overview and issues , 2004, J. Comput. Chem..

[66]  S. Grzesiek,et al.  NMRPipe: A multidimensional spectral processing system based on UNIX pipes , 1995, Journal of biomolecular NMR.

[67]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[68]  J. Wyman,et al.  LINKED FUNCTIONS AND RECIPROCAL EFFECTS IN HEMOGLOBIN: A SECOND LOOK. , 1964, Advances in protein chemistry.

[69]  B. Honig,et al.  On the environment of ionizable groups in globular proteins. , 1984, Journal of molecular biology.

[70]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

[71]  G. Hummer,et al.  Peptide loop-closure kinetics from microsecond molecular dynamics simulations in explicit solvent. , 2002, Journal of the American Chemical Society.

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

[73]  E. Lattman,et al.  Experimental measurement of the effective dielectric in the hydrophobic core of a protein. , 1997, Biophysical chemistry.

[74]  D. Shortle,et al.  Evidence for strained interactions between side-chains and the polypeptide backbone. , 1994, Journal of molecular biology.

[75]  Nathan A. Baker,et al.  Improving implicit solvent simulations: a Poisson-centric view. , 2005, Current opinion in structural biology.

[76]  R J Read,et al.  Crystallography & NMR system: A new software suite for macromolecular structure determination. , 1998, Acta crystallographica. Section D, Biological crystallography.

[77]  A. Mildvan,et al.  Environments and mechanistic roles of the tyrosine residues of delta 5-3-ketosteroid isomerase. , 1993, Biochemistry.

[78]  R. Srinivasan,et al.  LINUS: A hierarchic procedure to predict the fold of a protein , 1995, Proteins.

[79]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[80]  J. Pitera,et al.  Dielectric properties of proteins from simulation: the effects of solvent, ligands, pH, and temperature. , 2001, Biophysical journal.

[81]  E. Lattman,et al.  The crystal structure of the ternary complex of staphylococcal nuclease, Ca2+ and the inhibitor pdTp, refined at 1.65 Å , 1989, Proteins.