Low-temperature UV-vis absorption and Stark-effect hole-burning spectra of Zn substituted cytochrome c are studied experimentally and theoretically using quantum mechanical and Poisson-Boltzmann electrostatics models. Both the Q and Soret bands show resolved splitting at temperatures below similar to 180 K. The trend observed in the splittings when comparing cytochromes from different species is found to be the same as that observed for the Q(0,0) band of ferrous cytochrome c. The relative magnitudes of the Q and Soret splittings are found to be consistent with predictions based on Goutennan's four orbital model. For horse heart and yeast cytochrome c, which show the greatest difference in the UV-visible band splittings, Stark effect measurements on persistent spectral holes in the Q(0,0) band indicate that the protein-induced polarization is distinctly different for these two species. Incorporation of the protein electrostatic field as virtual point charges into quantum mechanical calculations utilizing the INDO/s semiempirical Hamiltonian is used to demonstrate that the effects of the protein on the heme electronic structure can be considerably different for the two proteins, consistent with the experimental observations.