Geometrical Analysis of Cytochrome c Unfolding

We have developed a geometrical model 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. Our calculations were carried out assuming that different regions of the polypeptide chain unfold synchronously. To test this assumption, we performed lattice Monte Carlo simulations to study systematically the possible importance of asynchronicity. Our 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.