Exploring the Interstitial Atom in the FeMo Cofactor of Nitrogenase: Insights from QM and QM/MM Calculations

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.

Radical SAM Enzymes and Metallocofactor Assembly: A Structural Point of View

This Review focuses on the structure-function relationship of radical S-adenosyl-l-methionine (SAM) enzymes involved in the assembly of metallocofactors corresponding to the active sites of [FeFe]-hydrogenase and nitrogenase [MoFe]-protein. It does not claim to correspond to an extensive review on the assembly machineries of these enzyme active sites, for which many good reviews are already available, but instead deals with the contribution of structural data to the understanding of their chemical mechanism (Buren et al. Chem. Rev.2020, 142 ( (25), ) 11006-11012; Britt et al. Chem. Sci.2020, 11 ( (38), ), 10313-10323). Hence, we will present the history and current knowledge about the radical SAM maturases HydE, HydG, and NifB as well as what, in our opinion, should be done in the near future to overcome the existing barriers in our understanding of this fascinating chemistry that intertwine organic radicals and organometallic complexes.

N2H2 binding to the nitrogenase FeMo cluster studied by QM/MM methods

We have made a systematic combined quantum mechanical and molecular mechanical (QM/MM) investigation of possible structures of the N₂ bound state of nitrogenase. We assume that N₂ is immediately protonated to a N₂H₂ state, thereby avoiding the problem of determining the position of the protons in the cluster. We have systematically studied both end-on and side-on structures, as well as both HNNH and NNH₂ states. Our results indicate that the binding of N₂H₂ is determined more by interactions and steric clashes with the surrounding protein than by the intrinsic preferences of the ligand and the cluster. The best binding mode with both the TPSS and B3LYP density-functional theory methods has trans-HNNH terminally bound to Fe2. It is stabilised by stacking of the substrate with His-195 and Ser-278. However, several other structures come rather close in energy (within 3-35 kJ/mol) at least in some calculations: The corresponding cis-HNNH structure terminally bound to Fe2 is second best with B3LYP. A structure with HNNH₂ terminally bound to Fe6 is second most stable with TPSS (where the third proton is transferred to the substrate from the homocitrate ligand). Structures with trans-HNNH, bound to Fe4 or Fe6, or cis-HNNH bound to Fe6 are also rather stable. Finally, with the TPSS functional, a structure with cis-HNNH side-on binding to the Fe3-Fe4-Fe5-Fe7 face of the cluster is also rather low in energy, but all side-on structures are strongly disfavoured by the B3LYP method.

What Is the Structure of the E4 Intermediate in Nitrogenase?

Nitrogenase is the only enzyme that can cleave the strong triple bond in N₂. The active site contains a complicated MoFe₇S₉C cluster. It is believed that it needs to accept four protons and electrons, forming the E₄ state, before it can bind N₂. However, there is no consensus on the atomic structure of the E₄ state. Experimental studies indicate that it should contain two hydride ions bridging two pairs of Fe ions, and it has been suggested that both hydride ions as well as the two protons bind on the same face of the cluster. On the other hand, density functional theory (DFT) studies have indicated that it is energetically more favorable with either three hydride ions or with a triply protonated carbide ion, depending on the DFT functional. We have performed a systematic combined quantum mechanical and molecular mechanical (QM/MM) study of possible E₄ states with two bridging hydride ions. Our calculations suggest that the most favorable structure has hydride ions bridging the Fe2/6 and Fe3/7 ion pairs. In fact, such structures are 14 kJ/mol more stable than structures with three hydride ions, showing that pure DFT functionals give energetically most favorable structures in agreement with experiments. An important reason for this finding is that we have identified a new type of broken-symmetry state that involves only two Fe ions with minority spin, in contrast to the previously studied states with three Fe ions with minority spin. The energetically best structures have the two hydride ions on different faces of the FeMo cluster, whereas better agreement with ENDOR data is obtained if they are on the same face; such structures are only 6–22 kJ/mol less stable.

Extremely large differences in DFT energies for nitrogenase models

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen avaiable for other organisms. It contains a complicated MoFe7S9C(homocitrate) cluster in its active site. Many computational studies with density-functional theory (DFT) of the nitrogenase enzyme have been presented, but they do not show any consensus - they do not even agree where the first four protons should be added, forming the central intermediate E4. We show that the prime reason for this is that different DFT methods give relative energies that differ by almost 600 kJ mol-1 for different protonation states. This is 4-30 times more than what is observed for other systems. The reason for this is that in some structures, the hydrogens bind to sulfide or carbide ions as protons, whereas in other structures they bind to the metals as hydride ions, changing the oxidation state of the metals, as well as the Fe-C, Fe-S and Fe-Fe distances. The energies correlate with the amount of Hartree-Fock exchange in the method, indicating a variation in the amount of static correlation in the structures. It is currently unclear which DFT method gives the best results for nitrogenase. We show that non-hybrid DFT functionals and TPSSh give the most accurate structures of the resting active site, whereas B3LYP and PBE0 give the best H2 dissociation energies. However, no DFT method indicates that a structure of E4 with two bridging hydride ions is lowest in energy, as spectroscopic experiments indicate.

A New Nitrogenase Mechanism Using a CFe8S9 Model: Does H2 Elimination Activate the Complex to N2 Addition to the Central Carbon Atom?

A truncated model of the FeMo cofactor is used to explore a new mechanism for the conversion of N2 to NH3 by the nitrogenase enzyme. After four initial protonation/reduction steps, the H4CFe8S9 cluster has two hydrogen atoms attached to sulfur, one hydrogen bridging two iron centers and one hydrogen bonded to carbon. The loss of the CH and FeHFe hydrogens as molecular hydrogen activates the cluster to addition of N2 to the carbon center. This unique step takes place at a nearly planar four-coordinate carbon center and leads to an intermediate with a significantly weakened N-N bond. A hydrogen attached to a sulfur atom is then transferred to the distal nitrogen atom. Additional prontonation/reduction steps are modeled by adding a hydrogen atom to sulfur and locating the transition states for transfer to nitrogen. The first NH3 is lost in a thermal neutral step, while the second step is endothermic. The loss of H2 activates the complex by reducing the barrier for N2 addition by 3.5 kcal/mol. Since this is the most difficult step in the mechanism, reducing the barrier for this step justifies the "extra expense" of H2 production.

Tracing the interstitial carbide of the nitrogenase cofactor during substrate turnover

The fate of the interstitial atom of the nitrogenase cofactor during substrate turnover has remained a topic of interest since the discovery of this atom more than a decade ago. In this study, we labeled the interstitial carbide atom with (14)C and (13)C isotopes and traced the fate of the isotope under turnover conditions. Our results show that the interstitial carbide cannot be exchanged upon turnover, nor can it be used as a substrate and incorporated into the products. These observations point to a role of the interstitial carbide in stabilizing the cofactor structure, although a function of this atom in indirectly modulating the reactivity of the cofactor or directly interacting with the substrate cannot be excluded.

QM and QM/MM simulations of proteins

Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.

Fischer-Tropsch chemistry at room temperature?

The unique catalytic activity of vanadium nitrogenase suggests a new direction for the direct production of biofuels from CO with either synthetic catalysts or nitrogenase-containing bacteria. The reduction of CO by V nitrogenase to light hydrocarbons shows striking similarities to the established Fischer-Tropsch process; however, the enzyme does not use H(2) directly for this reaction. ADP=adenosine diphosphate, ATP= adenosine triphosphate.

Nitrogenase structure and function relationships by density functional theory

Modern density functional theory has tremendous potential with matching popularity in metalloenzymology to reveal the unseen atomic and molecular details of structural data, spectroscopic measurements, and biochemical experiments by providing insights into unobservable structures and states, while also offering theoretical justifications for observed trends and differences. An often untapped potential of this theoretical approach is to bring together diverse experimental structural and reactivity information and allow for these to be critically evaluated at the same level. This is particularly applicable for the tantalizingly complex problem of the structure and molecular mechanism of biological nitrogen fixation. In this chapter we provide a review with extensive practical details of the compilation and evaluation of experimental data for an unbiased and systematic density functional theory analysis that can lead to remarkable new insights about the structure-function relationships of the iron-sulfur clusters of nitrogenase.