We have measured the low temperature (T) 20 K) absorption spectra of the N52A, N52V, N52I, Y67F, and N52AY67F mutants of ferrous Saccharomyces cerevisiae (baker's yeast) cytochrome c. All the bands in the Q(0)- and Q(v)-band region are split, and the intensity distributions among the split bands are highly asymmetric. The spectra were analyzed by a decomposition into Voigtian profiles. The spectral parameters thus obtained were further analyzed in terms of the vibronic coupling model of Schweitzer-Stenner and Bigman (Schweitzer-Stenner, R.; Bigman, D. J. Phys. Chem. B 2001, 7064-7073) to identify parameters related to electronic and vibronic perturbations of the heme macrocycle. We report that the electronic perturbation is of B-1g symmetry and reflects the heterogeneity of the electric field at the heme, that is, the difference between the gradients along the perpendicular N-Fe-N axis of the heme core. We found that all the investigated mutations substantially increase this electronic perturbation, so that the spectral properties become similar to those of horse heart cytochrome c. Moreover, the electronic perturbation was found to correlate nonlinearly with the enthalpy changes associated with the reduction of the heme iron. Group theoretical arguments are invoked to propose a simple model which explains how a perturbation of the obtained symmetry can stabilize the reduced state of the heme iron. Finally, vibronic coupling parameters obtained from the analysis of the Q(v)-band region suggest that the investigated mutations decrease the nonplanar deformations of the heme group. This finding was reproduced by a normal mode structural decomposition (NSD) analysis of the N52V and N52VY67F heme conformations obtained from a 1 ns molecular dynamics simulation. We argue that the reduced nonplanarity contributes to the stabilization of the reduced state.