The diffusional dynamics and vibrational spectroscopy of molecular hydrogen (H2) in myoglobin (Mb) is characterized. Hydrogen has been implicated in a number of physiologically relevant processes, including cellular aging or inflammation. Here, the internal diffusion through the protein matrix was characterized, and the vibrational spectroscopy was investigated using conventional empirical energy functions and improved models able to describe higher-order electrostatic moments of the ligand. Depending on the energy function used, H2 can occupy the same internal defects as already found for Xe or CO (Xe1 to Xe4 and B-state). Furthermore, four additional sites were found, some of which had been discovered in earlier simulation studies. Simulations using a model based on a Morse oscillator and distributed charges to correctly describe the molecular quadrupole moment of H2 indicate that the vibrational spectroscopy of the ligand depends on the docking site it occupies. This is consistent with the findings for CO in Mb from experiments and simulations. For H2, the maxima of the absorption spectra cover ∼20 cm−1 which are indicative of a pronounced Stark effect of the surrounding protein matrix on the vibrational spectroscopy of the ligand. Electronic structure calculations show that H2 forms a stable complex with the heme iron (stabilized by ∼−12 kcal/mol), but splitting of H2 is unlikely due to a high activation energy (∼50 kcal/mol).
Read full abstract