Time-Resolved EPR Study of H2 Reductive Elimination from the Photoexcited Nitrogenase Janus E4(4H) Intermediate

QM/MM Study of Partial Dissociation of S2B for the E2 Intermediate of Nitrogenase

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available for all lifeforms. Previous computational studies have given widely diverging results regarding the reaction mechanism of the enzyme. For example, some recent studies have suggested that one of the μ2-bridging sulfide ligands (S2B) may dissociate from one of the Fe ions when protonated in the doubly reduced and protonated E2 state, whereas other studies indicated that such half-dissociated states are unfavorable. We have examined how the relative energies of 26 structures of the E2 state depend on details of combined quantum mechanical and molecular mechanical (QM/MM) calculations. We show that the selection of the broken-symmetry state, the basis set, relativistic effects, the size of the QM system, relaxation of the surroundings, and the conformations of the bound protons may affect the relative energies of the various structures by up to 12, 22, 9, 20, 37, and 33 kJ/mol, respectively. However, they do not change the preferred type of structures. On the other hand, the choice of the DFT functional strongly affects the preferences. The hybrid B3LYP functional strongly prefers doubly protonation of the central carbide ion, but such a structure is not consistent with experimental EPR data. Other functionals suggest structures with a hydride ion, in agreement with the experiments, and show that the ion bridges between Fe2 and Fe6. Moreover, there are two structures of the same type that are degenerate within 1-5 kJ/mol, in agreement with the observation of two EPR signals. However, the pure generalized gradient approximation (GGA) functional TPSS favors structures with a protonated S2B also bridging Fe2 and Fe6, whereas r2SCAN (meta-GGA) and TPSSh (hybrid) prefer structures with S2B dissociated from Fe2 (but remaining bound to Fe6). The energy difference between the two types of structure is so small (7-18 kJ/mol) that both types need to be considered in future investigations of the mechanism of nitrogenase.

13C ENDOR Characterization of the Central Carbon within the Nitrogenase Catalytic Cofactor Indicates That the CFe6 Core Is a Stabilizing "Heart of Steel"

Substrates and inhibitors of Mo-dependent nitrogenase bind and react at Fe ions of the active-site FeMo-cofactor [7Fe-9S-C-Mo-homocitrate] contained within the MoFe protein α-subunit. The cofactor contains a CFe6 core, a carbon centered within a trigonal prism of six Fe, whose role in catalysis is unknown. Targeted 13C labeling of the carbon enables electron-nuclear double resonance (ENDOR) spectroscopy to sensitively monitor the electronic properties of the Fe-C bonds and the spin-coupling scheme adopted by the FeMo-cofactor metal ions. This report compares 13CFe6 ENDOR measurements for (i) the wild-type protein resting state (E0; α-Val70) to those of (ii) α-Ile70, (iii) α-Ala70-substituted proteins; (iv) crystallographically characterized CO-inhibited "hi-CO" state; (v) E4(4H) Janus intermediate, activated for N2 binding/reduction by accumulation of 4[e-/H+]; (vi) E4(2H)* state containing a doubly reduced FeMo-cofactor without Fe-bound substrates; and (vii) propargyl alcohol reduction intermediate having allyl alcohol bound as a ferracycle to FeMo-cofactor Fe6. All states examined, both S = 1/2 and 3/2 exhibited near-zero 13C isotropic hyperfine coupling constants, Ca = [-1.3 ↔ +2.7] MHz. Density functional theory computations and natural bond orbital analysis of the Fe-C bonds show that this occurs because a (3 spin-up/3 spin-down) spin-exchange configuration of CFe6 Fe-ion spins produces cancellation of large spin-transfers to carbon in each Fe-C bond. Previous X-ray diffraction and DFT both indicate that trigonal-prismatic geometry around carbon is maintained with high precision in all these states. The persistent structure and Fe-C bonding of the CFe6 core indicate that it does not provide a functionally dynamic (hemilabile) "beating heart"─instead it acts as "a heart of steel", stabilizing the structure of the FeMo-cofactor-active site during nitrogenase catalysis.

Review on the QM/MM Methodologies and Their Application to Metalloproteins

The multiscaling quantum mechanics/molecular mechanics (QM/MM) approach was introduced in 1976, while the extensive acceptance of this methodology started in the 1990s. The combination of QM/MM approach with molecular dynamics (MD) simulation, otherwise known as the QM/MM/MD approach, is a powerful and promising tool for the investigation of chemical reactions’ mechanism of complex molecular systems, drug delivery, properties of molecular devices, organic electronics, etc. In the present review, the main methodologies in the multiscaling approaches, i.e., density functional theory (DFT), semiempirical methodologies (SE), MD simulations, MM, and their new advances are discussed in short. Then, a review on calculations and reactions on metalloproteins is presented, where particular attention is given to nitrogenase that catalyzes the conversion of atmospheric nitrogen molecules N₂ into NH₃ through the process known as nitrogen fixation and the FeMo-cofactor.

The One-Electron Reduced Active-Site FeFe-Cofactor of Fe-Nitrogenase Contains a Hydride Bound to a Formally Oxidized Metal-Ion Core

The nitrogenase active-site cofactor must accumulate 4e^-/4H⁺ (E₄(4H) state) before N₂ can bind and be reduced. Earlier studies demonstrated that this E₄(4H) state stores the reducing-equivalents as two hydrides, with the cofactor metal-ion core formally at its resting-state redox level. This led to the understanding that N₂ binding is mechanistically coupled to reductive-elimination of the two hydrides that produce H₂. The state having acquired 2e^-/2H⁺ (E₂(2H)) correspondingly contains one hydride with a resting-state core redox level. How the cofactor accommodates addition of the first e^-/H⁺ (E₁(H) state) is unknown. The Fe-nitrogenase FeFe-cofactor was used to address this question because it is EPR-active in the E₁(H) state, unlike the FeMo-cofactor of Mo-nitrogenase, thus allowing characterization by EPR spectroscopy. The freeze-trapped E₁(H) state of Fe-nitrogenase shows an S = 1/2 EPR spectrum with g = [1.965, 1.928, 1.779]. This state is photoactive, and under 12 K cryogenic intracavity, 450 nm photolysis converts to a new and likewise photoactive S = 1/2 state (denoted E₁(H)*) with g = [2.009, 1.950, 1.860], which results in a photostationary state, with E₁(H)* relaxing to E₁(H) at temperatures above 145 K. An H/D kinetic isotope effect of 2.4 accompanies the 12 K E₁(H)/E₁(H)* photointerconversion. These observations indicate that the addition of the first e^-/H⁺ to the FeFe-cofactor of Fe-nitrogenase produces an Fe-bound hydride, not a sulfur-bound proton. As a result, the cluster metal-ion core is formally one-electron oxidized relative to the resting state. It is proposed that this behavior applies to all three nitrogenase isozymes.

Thermodynamically Favourable States in the Reaction of Nitrogenase without Dissociation of any Sulfide Ligand

We have used combined quantum mechanical and molecular mechanical (QM/MM) calculations to study the reaction mechanism of nitrogenase, assuming that none of the sulfide ligands dissociates. To avoid the problem that there is no consensus regarding the structure and protonation of the E₄ state, we start from a state where N₂ is bound to the cluster and is protonated to N₂H₂, after dissociation of H₂. We show that the reaction follows an alternating mechanism with HNNH (possibly protonated to HNNH₂) and H₂NNH₂ as intermediates and the two NH₃ products dissociate at the E₇ and E₈ levels. For all intermediates, coordination to Fe6 is preferred, but for the E₄ and E₈ intermediates, binding to Fe2 is competitive. For the E₄, E₅ and E₇ intermediates we find that the substrate may abstract a proton from the hydroxy group of the homocitrate ligand of the FeMo cluster, thereby forming HNNH₂, H₂NNH₂ and NH₃ intermediates. This may explain why homocitrate is a mandatory component of nitrogenase. All steps in the suggested reaction mechanism are thermodynamically favourable compared to protonation of the nearby His‐195 group and in all cases, protonation of the NE2 atom of the latter group is preferred.

Fluorescent iron‑sulfur centers: Photochemistry of the PetA Rieske protein from Aquifex aeolicus

Cytochrome bc₁ complexes are energy-transducing enzymes and key components of respiratory electron chains. They contain Rieske 2Fe2S proteins that absorb very weakly in the visible absorption region compared to the heme cofactors of the cytochromes, but are known to yield photoproducts. Here, the photoreactions of isolated Rieske proteins from the hyperthermophilic bacterium Aquifex aeolicus are studied in two redox states using ultrafast transient fluorescence and absorption spectroscopy. We provide evidence, for the first time in iron‑sulfur proteins, of very weak fluorescence of the excited state, in the oxidized as well as the reduced state. The excited states of the oxidized and reduced forms decay in 1.5 ps and 30 ps, respectively. In both cases they give rise to product states with lifetimes beyond 1 ns, reflecting photo-reduction of oxidized centers as well as photo-oxidation of reduced centers. Potential reaction partners are discussed and studied using site-directed mutagenesis. For the reduced state, a nearby disulfide bridge is suggested as an electron acceptor. The resulting photoproducts in either state may play a role in photoactivation processes.

Cleavage of cluster iron-sulfide bonds in cyclophane-coordinated FenSm complexes

Reaction of the tri(μ-sulfido)triiron(iii) tris(β-diketiminate) cyclophane complex, Fe3S3LEt/Me (1), or of the di(μ-sulfido)diiron(iii) complex Fe2S2HLEt/Me (5), with the related tri(bromide)triiron(ii) complex Fe3Br3LEt/Me (2) results in electron and ligand redistribution to yield the mixed-ligand multiiron complexes, including Fe3Br2SLEt/Me (3) and Fe2Br2SHLEt/Me (4). The cleavage and redistribution observed in these complexes is reminiscent of necessary Fe-S bond cleavage for substrate activation in nitrogenase enzymes, and provides a new perspective on the lability of Fe-S bonds in FeS clusters.

Electron Redistribution within the Nitrogenase Active Site FeMo-Cofactor During Reductive Elimination of H2 to Achieve N≡N Triple-Bond Activation

Nitrogen fixation by nitrogenase begins with the accumulation of four reducing equivalents at the active-site FeMo-cofactor (FeMo-co), generating a state (denoted E₄(4H)) with two [Fe-H-Fe] bridging hydrides. Recently, photolytic reductive elimination (re) of the E₄(4H) hydrides showed that enzymatic re of E₄(4H) hydride yields an H₂-bound complex (E₄(H₂,2H)), in a process corresponding to a formal 2-electron reduction of the metal-ion core of FeMo-co. The resulting electron-density redistribution from Fe-H bonds to the metal ions themselves enables N₂ to bind with concomitant H₂ release, a process illuminated here by QM/MM molecular dynamics simulations. What is the nature of this redistribution? Although E₄(H₂,2H) has not been trapped, cryogenic photolysis of E₄(4H) provides a means to address this question. Photolysis of E₄(4H) causes hydride-re with release of H₂, generating doubly reduced FeMo-co (denoted E₄(2H)*), the extreme limit of the electron-density redistribution upon formation of E₄(H₂,2H). Here we examine the doubly reduced FeMo-co core of the E₄(2H)* limiting-state by ¹H, ⁵⁷Fe, and ⁹⁵Mo ENDOR to illuminate the partial electron-density redistribution upon E₄(H₂,2H) formation during catalysis, complementing these results with corresponding DFT computations. Inferences from the E₄(2H)* ENDOR results as extended by DFT computations include (i) the Mo-site participates negligibly, and overall it is unlikely that Mo changes valency throughout the catalytic cycle; and (ii) two distinctive E₄(4H) ⁵⁷Fe signals are suggested as associated with structurally identified "anchors" of one bridging hydride, two others with identified anchors of the second, with NBO-analysis further identifying one anchor of each hydride as a major recipient of electrons released upon breaking Fe-H bonds.

The Spectroscopy of Nitrogenases

Nitrogenases are responsible for biological nitrogen fixation, a crucial step in the biogeochemical nitrogen cycle. These enzymes utilize a two-component protein system and a series of iron-sulfur clusters to perform this reaction, culminating at the FeMco active site (M = Mo, V, Fe), which is capable of binding and reducing N2 to 2NH3. In this review, we summarize how different spectroscopic approaches have shed light on various aspects of these enzymes, including their structure, mechanism, alternative reactivity, and maturation. Synthetic model chemistry and theory have also played significant roles in developing our present understanding of these systems and are discussed in the context of their contributions to interpreting the nature of nitrogenases. Despite years of significant progress, there is still much to be learned from these enzymes through spectroscopic means, and we highlight where further spectroscopic investigations are needed.