Dynamics of dehaloperoxidase-hemoglobin A derived from NMR relaxation spectroscopy and molecular dynamics simulation

Dehaloperoxidase-hemoglobin is the first hemoglobin identified with biologically-relevant oxidative functions, which include peroxidase, peroxygenase and oxidase activities. Herein we report a study of the protein backbone dynamics of DHP using heteronuclear NMR relaxation methods and molecular dynamics (MD) simulations to address the role of protein dynamics in switching from one function to another. The results show that DHP's backbone helical regions and turns have average order parameters of S2 = 0.87 ± 0.03 and S2 = 0.76 ± 0.08, respectively. Furthermore, DHP is primarily a monomer in solution based on the overall tumbling correlation time τm is 9.49 ± 1.65 ns calculated using the prolate diffusion tensor model in the program relax. A number of amino acid residues have significant Rex using the Lipari-Szabo model-free formalism. These include Lys3, Ile6, Leu13, Gln18, Arg32, Ser48, Met49, Thr56, Phe60, Arg69, Thr71 Cys73, Ala77, Asn81, Gly95, Arg109, Phe115, Leu127 and Met136, which may experience slow conformational motions on the microseconds-milliseconds time scale according to the model. Caution should be used when the model contains >4 fitting parameters. The program caver3.0 was used to identify tunnels inside DHP obtained from MD simulation snapshots that are consistent with the importance of the Xe binding site, which is located at the central intersection of the tunnels. These tunnels provide diffusion pathways for small ligands such as O2, H2O and H2O2 to enter the distal pocket independently of the trajectory of substrates and inhibitors, both of which are aromatic molecules.