A model for dinitrogen binding in the E4 state of nitrogenase† †Electronic supplementary information (ESI) available: Further information on all broken-symmetry solutions calculated for all models. Details of the QM/MM model preparation. Spin populations of all models. Localized orbital analysis of selected models. Geometry analysis of E0 state calculated with different functionals and electronic structure analysis. Cartesian coordinates for the QM regions of all optimized structures available as XYZ files. See DOI: 10.1039/c9sc03610e
Autor: | Thorhallsson, Albert Th., Benediktsson, Bardi, Bjornsson, Ragnar |
---|---|
Jazyk: | angličtina |
Rok vydání: | 2019 |
Předmět: | |
Zdroj: | Chemical Science |
ISSN: | 2041-6539 2041-6520 |
Popis: | QM/MM calculations are used to propose a new model for the E4 state of FeMoco and how N2 binding to this state may occur. Molybdenum nitrogenase is one of the most intriguing metalloenzymes in nature, featuring an exotic iron–molybdenum–sulfur cofactor, FeMoco, whose mode of action remains elusive. In particular, the molecular and electronic structure of the N2-binding E4 state is not known. In this study we present theoretical QM/MM calculations of new structural models of the E4 state of molybdenum-dependent nitrogenase and compare to previously suggested models for this enigmatic redox state. We propose two models as possible candidates for the E4 state. Both models feature two hydrides on the FeMo cofactor, bridging atoms Fe2 and Fe6 with a terminal sulfhydryl group on either Fe2 or Fe6 (derived from the S2B bridge) and the change in coordination results in local lower-spin electronic structure at Fe2 and Fe6. These structures appear consistent with the bridging hydride proposal put forward from ENDOR studies and are calculated to be lower in energy than other proposed models for E4 at the TPSSh-QM/MM level of theory. We critically analyze the DFT method dependency in calculations of FeMoco that has resulted in strikingly different proposals for this state. Importantly, dinitrogen binds exothermically to either Fe2 or Fe6 in our models, contrary to others, an effect rationalized via the unique ligand field (from the hydrides) at the Fe with an empty coordination site. A low-spin Fe site is proposed as being important to N2 binding. Furthermore, the geometries of these states suggest a feasible reductive elimination step that could follow, as experiments indicate. Via this step, two electrons are released, reducing the cofactor to yield a distorted 4-coordinate Fe2 or Fe6 that partially activates N2. We speculate that stabilization of an N2-bound Fe(i) at Fe6 (not found for Fe2 model) via reductive elimination is a crucial part of N2 activation in nitrogenases, possibly aided by the apical heterometal ion (Mo or V). By using protons from the sulfhydryl group (to regenerate the sulfide bridge between Fe2 and Fe6) and the nearby homocitrate hydroxy group, we calculate a plausible route to yield a diazene intermediate. This is found to be more favorable with the Fe6-bound model than the Fe2-bound model; however, this protonation is uphill in energy, suggesting protonation of N2 might occur later in the catalytic cycle or via another mechanism. |
Databáze: | OpenAIRE |
Externí odkaz: |