Uncovering specific electrostatic interactions in the denatured states of proteins.

The stability and folding of proteins are modulated by energetically significant interactions in the denatured state that is in equilibrium with the native state. These interactions remain largely invisible to current experimental techniques, however, due to the sparse population and conformational heterogeneity of the denatured-state ensemble under folding conditions. Molecular dynamics simulations using physics-based force fields can in principle offer atomistic details of the denatured state. However, practical applications are plagued with the lack of rigorous means to validate microscopic information and deficiencies in force fields and solvent models. This study presents a method based on coupled titration and molecular dynamics sampling of the denatured state starting from the extended sequence under native conditions. The resulting denatured-state pK(a)s allow for the prediction of experimental observables such as pH- and mutation-induced stability changes. I show the capability and use of the method by investigating the electrostatic interactions in the denatured states of wild-type and K12M mutant of NTL9 protein. This study shows that the major errors in electrostatics can be identified by validating the titration properties of the fragment peptides derived from the sequence of the intact protein. Consistent with experimental evidence, our simulations show a significantly depressed pK(a) for Asp(8) in the denatured state of wild-type, which is due to a nonnative interaction between Asp(8) and Lys(12). Interestingly, the simulation also shows a nonnative interaction between Asp(8) and Glu(48) in the denatured state of the mutant. I believe the presented method is general and can be applied to extract and validate microscopic electrostatics of the entire folding energy landscape.

[1]  Charles L Brooks,et al.  Exploring atomistic details of pH-dependent peptide folding , 2006, Proceedings of the National Academy of Sciences.

[2]  G. Hummer,et al.  Are current molecular dynamics force fields too helical? , 2008, Biophysical journal.

[3]  Joseph A Marsh,et al.  Improved structural characterizations of the drkN SH3 domain unfolded state suggest a compact ensemble with native-like and non-native structure. , 2007, Journal of molecular biology.

[4]  Charles L Brooks,et al.  Linking folding with aggregation in Alzheimer's β-amyloid peptides , 2007, Proceedings of the National Academy of Sciences.

[5]  Ray Luo,et al.  Interplay of secondary structures and side-chain contacts in the denatured state of BBA1. , 2004, The Journal of chemical physics.

[6]  I. André,et al.  pK(a) values for side-chain carboxyl groups of a PGB1 variant explain salt and pH-dependent stability. , 2007, Biophysical journal.

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

[8]  D. Raleigh,et al.  Mutational analysis demonstrates that specific electrostatic interactions can play a key role in the denatured state ensemble of proteins. , 2005, Journal of molecular biology.

[9]  Charles L Brooks,et al.  Toward the accurate first-principles prediction of ionization equilibria in proteins. , 2006, Biochemistry.

[10]  A. Elcock.  Realistic modeling of the denatured states of proteins allows accurate calculations of the pH dependence of protein stability. , 1999, Journal of molecular biology.

[11]  A. R. Fresht.  Structure and Mechanism in Protein Science: A Guide to Enzyme Catalysis and Protein Folding , 1999 .

[12]  Y. Duan,et al.  Folding free-energy landscape of villin headpiece subdomain from molecular dynamics simulations , 2007, Proceedings of the National Academy of Sciences.

[13]  C. Dobson.  Protein folding and misfolding , 2003, Nature.

[14]  Lisa J Lapidus,et al.  How fast is protein hydrophobic collapse? , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[15]  Frank Küster,et al.  Single-molecule spectroscopy of the temperature-induced collapse of unfolded proteins , 2009, Proceedings of the National Academy of Sciences.

[16]  Sheena E Radford,et al.  An expanding arsenal of experimental methods yields an explosion of insights into protein folding mechanisms , 2009, Nature Structural &Molecular Biology.

[17]  D. Raleigh,et al.  The unfolded state of the villin headpiece helical subdomain: computational studies of the role of locally stabilized structure. , 2006, Journal of molecular biology.

[18]  C. J. Bond,et al.  Characterization of residual structure in the thermally denatured state of barnase by simulation and experiment: description of the folding pathway. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[19]  Kresten Lindorff-Larsen,et al.  Determination of an ensemble of structures representing the denatured state of the bovine acyl-coenzyme a binding protein. , 2004, Journal of the American Chemical Society.

[20]  Lauren Wickstrom,et al.  Secondary structure bias in generalized Born solvent models: comparison of conformational ensembles and free energy of solvent polarization from explicit and implicit solvation. , 2007, The journal of physical chemistry. B.

[21]  Terrence G. Oas,et al.  Backbone Dynamics of the Monomeric λ Repressor Denatured State Ensemble under Nondenaturing Conditions , 2007 .

[22]  H. Bosshard,et al.  Inverse electrostatic effect: electrostatic repulsion in the unfolded state stabilizes a leucine zipper. , 2004, Biochemistry.

[23]  Benjamin Schuler,et al.  Ultrafast dynamics of protein collapse from single-molecule photon statistics , 2007, Proceedings of the National Academy of Sciences.

[24]  F A Quiocho,et al.  Carboxypeptidase A: a protein and an enzyme. , 1971, Advances in protein chemistry.

[25]  Huan-Xiang Zhou,et al.  A Gaussian-chain model for treating residual charge–charge interactions in the unfolded state of proteins , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[26]  Jana K. Shen,et al.  Predicting pKa values with continuous constant pH molecular dynamics. , 2009, Methods in enzymology.

[27]  D. Raleigh,et al.  Thermodynamics and kinetics of non-native interactions in protein folding: a single point mutant significantly stabilizes the N-terminal domain of L9 by modulating non-native interactions in the denatured state. , 2004, Journal of molecular biology.

[28]  William Swope,et al.  Understanding folding and design: Replica-exchange simulations of ``Trp-cage'' miniproteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[29]  Michael R. Shirts,et al.  Native-like mean structure in the unfolded ensemble of small proteins. , 2002, Journal of molecular biology.

[30]  C. Dobson,et al.  Formation of native and non-native interactions in ensembles of denatured ACBP molecules from paramagnetic relaxation enhancement studies. , 2005, Journal of molecular biology.

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

[32]  V. Ramakrishnan,et al.  Ribosomal protein L9: a structure determination by the combined use of X-ray crystallography and NMR spectroscopy. , 1996, Journal of molecular biology.

[33]  Kevin L. Shaw,et al.  Charge–charge interactions influence the denatured state ensemble and contribute to protein stability , 2000, Protein science : a publication of the Protein Society.

[34]  C. Brooks,et al.  Constant pH molecular dynamics with proton tautomerism. , 2005, Biophysical journal.

[35]  B. Kuhlman,et al.  pKa values and the pH dependent stability of the N-terminal domain of L9 as probes of electrostatic interactions in the denatured state. Differentiation between local and nonlocal interactions. , 1999, Biochemistry.

[36]  D. Raleigh,et al.  Electrostatic interactions in the denatured state ensemble: their effect upon protein folding and protein stability. , 2008, Archives of biochemistry and biophysics.

[37]  V L Arcus,et al.  pKA values of carboxyl groups in the native and denatured states of barnase: the pKA values of the denatured state are on average 0.4 units lower than those of model compounds. , 1995, Biochemistry.

[38]  B. Kuhlman,et al.  Structure and stability of the N-terminal domain of the ribosomal protein L9: evidence for rapid two-state folding. , 1998, Biochemistry.

[39]  A. Garcia,et al.  Helix‐coil transition of alanine peptides in water: Force field dependence on the folded and unfolded structures , 2005, Proteins.

[40]  B. Brooks,et al.  An analysis of the accuracy of Langevin and molecular dynamics algorithms , 1988 .

[41]  Jana K. Shen.  A method to determine residue-specific unfolded-state pKa values from analysis of stability changes in single mutant cycles. , 2010, Journal of the American Chemical Society.

[42]  Electrostatic effects in unfolded staphylococcal nuclease , 2008, Protein science : a publication of the Protein Society.

[43]  C. Brooks,et al.  Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field. , 2006, Journal of the American Chemical Society.

[44]  C. Tanford.  Protein denaturation. , 1968, Advances in protein chemistry.

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

[46]  X. Daura,et al.  The Key to Solving the Protein-Folding Problem Lies in an Accurate Description of the Denatured State. , 2001, Angewandte Chemie.

[47]  A. Fersht,et al.  Perturbed pKA-values in the denatured states of proteins. , 1995, Journal of molecular biology.