Quantum chemical calculations of the iron electron structure and 57 Fe quadrupole splitting were made by density functional theory and Xα discrete variation method for the rough heme models for α-and β-subunits in deoxyhemoglobin and for deoxymyoglobin accounting stereochemical differences of the active sites in native proteins. The calculations revealed differences of quadrupole splitting temperature dependences for three models indicating sensitivity of quadrupole splitting and Fe(II) electronic structure to small variations of iron stereochemistry.