Propionates, as peripheral groups of the heme active center in hemeproteins have been described to contribute in the modulation of heme reactivity and ligand selection. These electronic characteristics prompted the question of whether the presence of hydrogen bonding networks between propionates and distal amino acids present in the heme ligand moiety can modulate physiological relevant events, like ligand binding association and dissociation activities. Here, the role of these networks was evaluated by NMR spectroscopy using the hemoglobin I PheB10Tyr mutant from Lucina pectinata as model for TyrB10 and GlnE7 hemeproteins. H-1-NMR results for the rHblCN PheB10Tyr derivative showed chemical shifts of TyrB10 OH eta at 31.00 ppm, GlnE7 G(epsilon 1)H/N epsilon 2H at 10.66 ppm/-3.27 ppm, and PheE11 C delta H at 11.75 ppm, indicating the presence of a crowded, collapsed, and constrained distal pocket. Strong dipolar contacts and inter-residues crosspeaks between GlnE7/6-propionate group, GlnE7/TyrB10 and TyrB10/CN suggest that this hydrogen bonding network loop between GlnE7, TyrB10, 6-propionate group, and the heme ligand contribute significantly to the modulation of the heme iron electron density as well as the ligand stabilization mechanism. Therefore, the network loop presented here support the fact that the electron withdrawing character of the hydrogen bonding is controlled by the interaction of the propionates and the nearby electronic environments contributing to the modulation of the heme electron density state. Thus, we hypothesize that in hemeproteins with similar electrostatic environment the flexibility of the heme-6-propionate promotes a hydrogen bonding network loop between the 6-propionate, the heme ligand and nearby amino acids, tailoring in this way the electron density in the heme-ligand moiety. (c) 2012 Elsevier Inc. All rights reserved.