Geometrical analysis of cytochrome c unfolding

Abstract

A geometrical model has been developed to study the unfolding of iso-1 cytochrome c. The model draws on the crystallographic data reported for this protein. These data were used to calculate the distance between specific residues in the folded state, and in a sequence of extended states defined by n = 3, 5, 7, 9, 11, 13, and 15 residue units. Exact calculations carried out for each of the 103 residues in the polypeptide chain demonstrate that different regions of the chain have different unfolding histories. Regions where there is a persistence of compact structures can be identified, and this geometrical characterization is fully consistent with analyses of time-resolved fluorescence energy-transfer (TrFET) data using dansyl-derivatized cysteine side-chain probes at positions 39, 50, 66, 85, and 99. The calculations were carried out assuming that different regions of the polypeptide chain unfold synchronously. To test this assumption, lattice Monte Carlo simulations were performed to study systematically the possible importance of asynchronicity. Calculations show that small departures from synchronous dynamics can arise if displacements of residues in the main body of the chain are much more sluggish than near-terminal residues.

Keywords:

• protein folding
• cytochrome c
• Monte Carlo simulations

Acknowledgements

Research at Caltech is supported by the National Institutes of Health (DK019038, GM068461), the National Science Foundation (CHE-0802907), and the Arnold and Mabel Beckman Foundation.