Molecular dynamics calculations have been performed over long trajectories with the inclusion of explicit solvent molecules on the reduced and the oxidized states of yeast iso-l-cytochrome c. The resulting structures have been analyzed and compared both in terms of structural properties and dynamical behavior. The structure of the buried water molecules around the heme has been also analyzed for the two oxidation states and compared with the experimental observations on the X-ray and the solution NMR structures. From the overall analysis we learn that, as also observed experimentally through NMR, no significant differences are present between the structures of the two oxidation states beside the arrangement of a few side chains. Also the internal mobility is similar for the two oxidation states, even if interesting differences are observed for some residues, as for Tyr67, a residue present at the heme site. The location and the mobility of the ordered water molecules, observed in solution by NMR, are completely reproduced in the molecular dynamics simulations, which have been able to predict the different displacements of the catalytically relevant water molecule WAT166, similar to those observed in solution for the two oxidation states, at variance with that observed in the starting crystallographic structures. The relevance of these findings with respect to the prediction of structural and dynamical properties is discussed.
