Mycobacterium heme utilization degrader (MhuD) is a heme-degrading protein from Mycobacterium tuberculosis responsible for extracting the essential nutrient iron from host-derived heme. MhuD has been previously shown to produce unique organic products compared to those of canonical heme oxygenases (HOs) as well as those of the IsdG/I heme-degrading enzymes from Staphylococcus aureus. Here, we report the X-ray crystal structure of cyanide-inhibited MhuD (MhuD–heme–CN) as well as detailed 1H nuclear magnetic resonance (NMR), UV/vis absorption, and magnetic circular dichroism (MCD) spectroscopic characterization of this species. There is no evidence for an ordered network of water molecules on the distal side of the heme substrate in the X-ray crystal structure, as was previously reported for canonical HOs. The degree of heme ruffling in the crystal structure of MhuD is greater than that observed for HO and less than that observed for IsdI. As a consequence, the Fe 3dxz-, 3dyz-, and 3dxy-based MOs are very close in energy, and the room-temperature 1H NMR spectrum of MhuD–heme–CN is consistent with population of both a 2Eg electronic state with a (dxy)2(dxz,dyz)3 electron configuration, similar to the ground state of canonical HOs, and a 2B2g state with a (dxz,dyz)4(dxy)1 electron configuration, similar to the ground state of cyanide-inhibited IsdI. Variable temperature, variable field MCD saturation magnetization data establishes that MhuD–heme–CN has a 2B2g electronic ground state with a low-lying 2Eg excited state. Our crystallographic and spectroscopic data suggest that there are both structural and electronic contributions to the α-meso regioselectivity of MhuD-catalyzed heme cleavage. The structural distortion of the heme substrate observed in the X-ray crystal structure of MhuD–heme–CN is likely to favor cleavage at the α- and γ-meso carbons, whereas the spin density distribution may favor selective oxygenation of the α-meso carbon.

For decades it has been known that an out-of-plane ruffling distortion of heme perturbs its UV-vis absorption (Abs) spectrum, but whether increased ruffling induces a red or blue shift of the Soret band has remained a topic of debate. This debate has been resolved by the spectroscopic and computational characterization of Mycobacterium tuberculosis MhuD presented here, an enzyme that converts heme, oxygen, and reducing equivalents to nonheme iron and mycobilin. W66F and W66A MhuD have been characterized using (1)H nuclear magnetic resonance, Abs, and magnetic circular dichroism spectroscopies, and the data have been used to develop an experimentally validated theoretical model of ruffled, ferric heme. The PBE density functional theory (DFT) model that has been developed accurately reproduces the observed spectral changes from wild type enzyme, and the underlying quantum mechanical origins of these ruffling-induced changes were revealed by analyzing the PBE DFT description of the electronic structure. Small amounts of heme ruffling have no influence on the energy of the Q-band and blue-shift the Soret band due to symmetry-allowed mixing of the Fe 3dxy and porphyrin a2u orbitals. Larger amounts of ruffling red-shift both the Q and Soret bands due to disruption of π-bonding within the porphyrin ring.

Recent work by several groups has established that MhuD, IsdG, and IsdI are non-canonical heme oxygenases that induce significant out-of-plane ruffling distortions of their heme substrates enroute to mycobilin or staphylobilin formation. However, clear explanations for the observations of "nested" S = ½ VTVH MCD saturation magnetization curves at cryogenic temperatures, and exchange broadened (1)H NMR resonances at physiologically-relevant temperatures have remained elusive. Here, MCD and NMR data have been acquired for F23A and F23W MhuD-heme-CN, in addition to MCD data for IsdI-heme-CN, in order to complete assembly of a library of spectroscopic data for cyanide-inhibited ferric heme with a wide range of ruffling deformations. The spectroscopic data were used to evaluate a number of computational models for cyanide-inhibited ferric heme, which ultimately led to the development of an accurate NEVPT2/CASSCF model. The resulting model has a shallow, double-well potential along the porphyrin ruffling coordinate, which provides clear explanations for the unusual MCD and NMR data. The shallow, double-well potential also implies that MhuD-, IsdG-, and IsdI-bound heme is dynamic, and the functional implications of these dynamics are discussed.