The iron-molybdenum cofactor (FeMoco) is responsible for dinitrogen reduction in Mo nitrogenase. Unlike the resting state, E0 , reduced states of FeMoco are much less well characterized. The E2 state has been proposed to contain a hydride but direct spectroscopic evidence is still lacking. The E2 state can, however, relax back the E0 state via a H2 side-reaction, implying a hydride intermediate prior to H2 formation. This E2 →E0 pathway is one of the primary mechanisms for H2 formation under low-electron flux conditions. In this study we present an exploration of the energy surface of the E2 state. Utilizing both cluster-continuum and QM/MM calculations, we explore various classes of E2 models: including terminal hydrides, bridging hydrides with a closed or open sulfide-bridge, as well as models without. Importantly, we find the hemilability of a protonated belt-sulfide to strongly influence the stability of hydrides. Surprisingly, non-hydride models are found to be almost equally favorable as hydride models. While the cluster-continuum calculations suggest multiple possibilities, QM/MM suggests only two models as contenders for the E2 state. These models feature either i) a bridging hydride between Fe2 and Fe6 and an open sulfide-bridge with terminal SH on Fe6 (E2 -hyd) or ii) a double belt-sulfide protonated, reduced cofactor without a hydride (E2 -nonhyd). We suggest both models as contenders for the E2 redox state and further calculate a mechanism for H2 evolution. The changes in electronic structure of FeMoco during the proposed redox-state cycle, E0 →E1 →E2 →E0 , are discussed.
Keywords: cofactors; density functional calculations; hydrides; nitrogenases; quantum chemistry.
© 2021 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.