Stability and fluctuations of amide hydrogen bonds in a bacterial cytochrome c: a molecular dynamics study

Molecular dynamics (MD) simulations on a bacterial cytochrome c were performed to investigate the lifetime and fluctuations of backbone hydrogen bonds and to correlate these data with protection factors for hydrogen exchange measured by NMR spectroscopy (Bartalesi et al. in Biochemistry, 42:10923–10930, 2003). The MD simulations provide a consistent pattern in that long lifetimes of hydrogen bonds go along with small amplitude fluctuations. In agreement with experiments, differences in stability were found with a rather flexible N-terminal segment as compared with a more rigid C-terminal part. Protection factors of backbone hydrogen exchange correlate strongly with the number of contacts but also with hydrogen-bond occupancy, hydrogen-bond survival times, as well as the inverse of fluctuations of backbone atoms and hydrogen-bond lengths derived from MD simulation data. We observed a conformational transition in the C-terminal loop, and significant motion in the N-terminal loop, which can be interpreted as being the structural units involved in the onset of the protein unfolding process in agreement with experimental evidence on mitochondrial cytochrome c.

[1]  L Mayne,et al.  Determinants of protein hydrogen exchange studied in equine cytochrome c , 1998, Protein science : a publication of the Protein Society.

[2]  S. Englander,et al.  Cooperative omega loops in cytochrome c: role in folding and function. , 2003, Journal of molecular biology.

[3]  S. Englander,et al.  Cytochrome c folding pathway: Kinetic native-state hydrogen exchange , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[4]  I. Muegge,et al.  Residence Times and Lateral Diffusion of Water at Protein Surfaces: Application to BPTI , 1995 .

[5]  E. Knapp,et al.  A consistent picture of protein dynamics. , 1984, Proceedings of the National Academy of Sciences of the United States of America.

[6]  A. D. Robertson,et al.  EX1 hydrogen exchange and protein folding. , 2004, Biochemistry.

[7]  Sozanne R. Solmaz,et al.  Electron transfer between yeast cytochrome bc(1) complex and cytochrome c: a structural analysis. , 2002, Biochimica et biophysica acta.

[8]  Antonio Rosato,et al.  Structure and dynamics of reduced Bacillus pasteurii cytochrome c: oxidation state dependent properties and implications for electron transfer processes. , 2003, Biochemistry.

[9]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[10]  L Mayne,et al.  Protein folding studied using hydrogen-exchange labeling and two-dimensional NMR. , 1992, Annual review of biophysics and biomolecular structure.

[11]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[12]  J. V. Van Beeumen,et al.  Cytochrome c-553 from the alkalophilic bacterium Bacillus pasteurii has the primary structure characteristics of a lipoprotein. , 1999, Biochemical and biophysical research communications.

[13]  A M Lesk,et al.  Helix movements and the reconstruction of the haem pocket during the evolution of the cytochrome c family. , 1985, Journal of molecular biology.

[14]  S. Marqusee,et al.  Molten globule unfolding monitored by hydrogen exchange in urea. , 1998, Biochemistry.

[15]  A J Wand,et al.  Local stability and dynamics of apocytochrome b562 examined by the dependence of hydrogen exchange on hydrostatic pressure. , 1998, Biochemistry.

[16]  R. L. Baldwin,et al.  Kinetics of hydrogen bond breakage in the process of unfolding of ribonuclease A measured by pulsed hydrogen exchange. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[17]  W. Kabsch.  A solution for the best rotation to relate two sets of vectors , 1976 .

[18]  S. Englander,et al.  Folding units govern the cytochrome c alkaline transition. , 2003, Journal of molecular biology.

[19]  Keith S. Wilson,et al.  Crystal structure of oxidized Bacillus pasteurii cytochrome c553 at 0.97-A resolution. , 2000, Biochemistry.

[20]  C. Woodward,et al.  Hydrogen exchange identifies native-state motional domains important in protein folding. , 1993, Biochemistry.

[21]  Yawen Bai,et al.  A protein folding pathway with multiple folding intermediates at atomic resolution. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[22]  Michele Vendruscolo,et al.  Rare fluctuations of native proteins sampled by equilibrium hydrogen exchange. , 2003, Journal of the American Chemical Society.

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

[24]  R. Li,et al.  The hydrogen exchange core and protein folding , 1999, Protein science : a publication of the Protein Society.

[25]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[26]  T. Sosnick,et al.  Protein folding intermediates: native-state hydrogen exchange. , 1995, Science.

[27]  G. Hummer,et al.  Conformational dynamics of cytochrome c: Correlation to hydrogen exchange , 1999, Proteins.

[28]  Tracy M. Handel,et al.  Detection of rare partially folded molecules in equilibrium with the native conformation of RNaseH , 1996, Nature Structural Biology.

[29]  C L Brooks,et al.  Characterization of "native" apomyoglobin by molecular dynamics simulation. , 1992, Journal of molecular biology.

[30]  F. Guerlesquin,et al.  Mapping the cytochrome c 553 interacting site using 1H and 15N NMR , 1999, FEBS letters.

[31]  S. Englander,et al.  How cytochrome c folds, and why: submolecular foldon units and their stepwise sequential stabilization. , 2004, Journal of molecular biology.

[32]  V. Daggett,et al.  Molecular dynamics simulations of protein unfolding and limited refolding: characterization of partially unfolded states of ubiquitin in 60% methanol and in water. , 1995, Journal of molecular biology.

[33]  G. Sciara,et al.  NMR Solution Structure, Backbone Mobility, and Homology Modeling of c‐Type Cytochromes from Gram‐Positive Bacteria , 2002, Chembiochem : a European journal of chemical biology.

[34]  Milan Hodoscek,et al.  Unfolding of the cold shock protein studied with biased molecular dynamics , 2003, Proteins.

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

[36]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[37]  S Walter Englander,et al.  Protein hydrogen exchange mechanism: Local fluctuations , 2003, Protein science : a publication of the Protein Society.

[38]  A. Fersht.  Structure and mechanism in protein science , 1998 .

[39]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[40]  E. Knapp,et al.  Tuning Heme Redox Potentials in the Cytochrome c Subunit of Photosynthetic Reaction Centers* , 2003, Journal of Biological Chemistry.