The different features of the nature of heme iron-proximal binding in methemoglobin and metmyoglobin were investigated by high resolution proton nuclear magnetic resonance spectroscopy at 220 MHz. From the comparison of the hyperfine-shifted heme methyl resonances of various methemoglobin and metmyoglobin derivatives it was revealed that the metHb derivatives exhibit smaller hyperfine shifts than the corresponding metMb derivatives. This observation was interpreted in terms of the stronger interaction between proximal histidine and ferric heme iron in methemoglobin than in metmyoglobin. The spectral contributions from the alpha and beta subunits in methemoglobin in high spin state were successfully separated by the analysis of the temperature dependent heme methyl shifts of the methemoglobin azide complex which is in the high and low thermal spin equilibrium. It was shown that the beta heme methyl resonances have smaller hyperfine shifts than the alpha heme ones in ferric high spin state. This tendency in methemoglobin derivatives was also interpreted in terms of the stronger iron-histidine binding in the beta subunits than in the alpha subunits within the tetramer. From the comparisons between methemoglobin and metmyoglobin, and between the alpha and beta subunits in methemoglobin, the order of the strength of the iron-histidine bond interaction was deduced as follows: metmyoglobin less than alpha subunits in methemoglobin less than or equal to beta units in methemoglobin.