The structural and functional properties of the active site of the bacterial hemoglobin from Thermobifida fusca are largely determined by three polar 'amino acids: TrpG8, TyrCD1, and TyrB10. We have exploited the availability of a combinatorial set of mutants, in each of which these three amino acids have been singly, doubly, or triply replaced by a Phe residue, to perform a detailed study on H-bonding interactions between the protein and heme-bound fluoride. By appropriate choice of the excitation conditions, v(Fe-F) stretching bands have been detected in the resonance Raman spectra. In the wildtype type protein and one of the mutants, two v(Fe-F) bands have been observed and assigned to the presence of two protein conformers where fluoride is singly or doubly H-bonded. Furthermore, by plotting the CT1 charge-transfer transition energy vs the v(Fe-F) wavenumbers, an empirical correlation has been found. The data are well fitted by a straight line with a positive slope. The position along the correlation line can be considered as a novel, general spectroscopic indicator of the extent of H-bonding in the active site of heme proteins. In agreement with the spectroscopic results, we have observed that the rate of ligand dissociation in stopped-flow kinetic measurements progressively increases upon substitution of the H-bonding amino acids. Molecular dynamics simulations have been performed on the fluoride complexes of native and mutated forms, indicating the prevalent interactions at the active site. All the techniques yield evidence that TrpG8 and TyrCD1 can form strong H bonds with fluoride, whereas TyrB10 plays only a minor role in the stabilization of the ligand.