Model for Acetylene Reduction by Nitrogenase Derived from Density Functional Theory

A diazotrophy-ammoniotrophy dual growth model for the sulfate reducing bacterium Desulfovibrio vulgaris var. Hildenborough

Sulfate reducing bacteria (SRB) comprise one of the few prokaryotic groups in which biological nitrogen fixation (BNF) is common. Recent studies have highlighted SRB roles in N cycling, particularly in oligotrophic coastal and benthic environments where they could contribute significantly to N input. Most studies of SRB have focused on sulfur cycling and SRB growth models have primarily aimed at understanding the effects of electron sources, with N usually provided as fixed-N (nitrate, ammonium). Mechanistic links between SRB nitrogen-fixing metabolism and growth are not well understood, particularly in environments where fixed-N fluctuates. Here, we investigate diazotrophic growth of the model sulfate reducer Desulfovibrio vulgaris var. Hildenborough under anaerobic heterotrophic conditions and contrasting N availabilities using a simple cellular model with dual ammoniotrophic and diazotrophic modes. The model was calibrated using batch culture experiments with varying initial ammonium concentrations (0-3000 µM) and acetylene reduction assays of BNF activity. The model confirmed the preferential usage of ammonium over BNF for growth and successfully reproduces experimental data, with notably clear bi-phasic growth curves showing an initial ammoniotrophic phase followed by onset of BNF. Our model enables quantification of the energetic cost of each N acquisition strategy and indicates the existence of a BNF-specific limiting phenomenon, not directly linked to micronutrient (Mo, Fe, Ni) concentration, by-products (hydrogen, hydrogen sulfide), or fundamental model metabolic parameters (death rate, electron acceptor stoichiometry). By providing quantitative predictions of environment and metabolism, this study contributes to a better understanding of anaerobic heterotrophic diazotrophs in environments with fluctuating N conditions.

Carbon Monoxide Binding to the Iron-Molybdenum Cofactor of Nitrogenase: a Detailed Quantum Mechanics/Molecular Mechanics Investigation

Carbon monoxide (CO) is a well-known inhibitor of nitrogenase activity. Under turnover conditions, CO binds to FeMoco, the active site of Mo nitrogenase. Time-resolved IR measurements suggest an initial terminal CO at 1904 cm^-1 that converts to a bridging CO at 1715 cm^-1, and an X-ray structure shows that CO can displace one of the bridging belt sulfides of FeMoco. However, the CO-binding redox state(s) of FeMoco (E_n) and the role of the protein environment in stabilizing specific CO-bound intermediates remain elusive. In this work, we carry out an in-depth analysis of the CO-FeMoco interaction based on quantum chemical calculations addressing different aspects of the electronic structure. (1) The local electronic structure of the Fe-CO bond is studied through diamagnetically substituted FeMoco. (2) A cluster model of FeMoco within a polarizable continuum illustrates how CO binding may affect the spin-coupling between the metal centers. (3) A QM/MM model incorporates the explicit influence of the amino acid residues surrounding FeMoco in the MoFe protein. The QM/MM model predicts both a terminal and a bridging CO in the E₁ redox state. The scaled calculated CO frequencies (1922 and 1716 cm^-1, respectively) are in good agreement with the experimentally observed IR bands supporting CO binding to the E₁ state. Alternatively, an E₂ state QM/MM model, which has the same atomic structure as the CO-bound X-ray structure, features a semi-bridging CO with a scaled calculated frequency (1718 cm^-1) similar to the bridging CO in the E₁ model.

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.

Reduction of Substrates by Nitrogenases

Nitrogenase is the enzyme that catalyzes biological N₂ reduction to NH₃. This enzyme achieves an impressive rate enhancement over the uncatalyzed reaction. Given the high demand for N₂ fixation to support food and chemical production and the heavy reliance of the industrial Haber-Bosch nitrogen fixation reaction on fossil fuels, there is a strong need to elucidate how nitrogenase achieves this difficult reaction under benign conditions as a means of informing the design of next generation synthetic catalysts. This Review summarizes recent progress in addressing how nitrogenase catalyzes the reduction of an array of substrates. New insights into the mechanism of N₂ and proton reduction are first considered. This is followed by a summary of recent gains in understanding the reduction of a number of other nitrogenous compounds not considered to be physiological substrates. Progress in understanding the reduction of a wide range of C-based substrates, including CO and CO₂, is also discussed, and remaining challenges in understanding nitrogenase substrate reduction are considered.

Survey of the Geometric and Electronic Structures of the Key Hydrogenated Forms of FeMo-co, the Active Site of the Enzyme Nitrogenase: Principles of the Mechanistically Significant Coordination Chemistry

The enzyme nitrogenase naturally hydrogenates N2 to NH3, achieved through the accumulation of H atoms on FeMo-co, the Fe7MoS9C(homocitrate) cluster that is the catalytically active site. Four intermediates, E1H1, E2H2, E3H3, and E4H4, carry these hydrogen atoms. I report density functional calculations of the numerous possibilities for the geometric and electronic structures of these poly-hydrogenated forms of FeMo-co. This survey involves more than 100 structures, including those with bound H2, and assesses their relative energies and most likely electronic states. Twelve locations for bound H atoms in the active domain of FeMo-co, including Fe–H–Fe and Fe–H–S bridges, are studied. A significant result is that transverse Fe–H–Fe bridges (transverse to the pseudo-threefold axis of FeMo-co and shared with triply-bridging S) are not possible geometrically unless the S is hydrogenated to become doubly-bridging. The favourable Fe–H–Fe bridges are shared with doubly-bridging S. ENDOR data for an E4H4 intermediate trapped at low temperature, and interpretations in terms of the geometrical and electronic structure of E4H4, are assessed in conjunction with the calculated possibilities. The results reported here yield a set of 24 principles for the mechanistically significant coordination chemistry of H and H2 on FeMo-co, in the stages prior to N2 binding.

Reactivity toward Unsaturated Small Molecules of Thiolate-Bridged Diiron Hydride Complexes

In the presence of 1 equiv of ^tBuNC, the homolytic cleavage of the Fe^III-H bond in the diiron terminal hydride complex [Cp*Fe( t-H)(μ-η²:η⁴-bdt)FeCp*][BF₄] (1[BF₄]) smoothly took place to release 1/2 H₂, followed by binding of a ^tBuNC group to the unsaturated Fe^II center. Interestingly, upon exposure of 1[BF₄] to 1 atm of acetylene, the isomerization process of the hydride ligand from the terminal to bridging coordination site was unaffected. Upon treatment of the diiron hydride bridged complex 2[BF₄] with acetylene at 30 °C, two Fe^III-H bonds were broken, and then an acetylene molecule was coordinated to the diiron centers in a novel μ-η²:η² side-on fashion. In the above reaction system, the hydride ligands whether terminal or bridging all play a role as the electron donor for the reduction of the diiron centers from Fe^IIIFe^III to Fe^IIIFe^II. These reaction patterns are reminiscent of the vital E₄ state responsible for N₂ binding and H₂ liberation in the catalytic cycle of nitrogenase, which contains two {Fe-H-Fe} motifs as electron reservoirs for the reduction of the iron centers. Differently, when treating 1[BF₄] with TMSN₃, the terminal hydride ligand was inserted into the azide group to give a diiron amide complex 4[BF₄] in moderate yield.

Mechanisms of the S/CO/Se interchange reactions at FeMo-co, the active site cluster of nitrogenase

The active site of the N2 fixing enzyme nitrogenase is a C-centred Fe7MoS cluster (FeMo-co) containing a trigonal prism of six Fe atoms connected by a central belt of three doubly-bridging S atoms. The trigonal faces of the prism are capped via triply-bridging S atoms to Fe1 at one end and Mo at the other end. One of the central belt atoms, S2B, considered to be important in the chemical mechanism of the enzyme, has been shown by Spatzal, Rees et al. to undergo substitution by CO, and also substitution by Se in the presence of SeCN(-), under turnover conditions. Further, when turning over under C2H2 or N2/CO there is migration of Se to the other two belt bridging positions. These reactions are extraordinary, and unprecedented in metal chalcogenide cluster chemistry. Using density functional simulations, mechanisms for all of these reactions have been developed, involving the small molecules SCO, SeCO, C2H2S, C2H2Se, SeCN(-), SCN(-) functioning as carriers of S and Se atoms. The possibility that the S2B bridge position is vacant is discounted, because the barrier to formation of a bridge-void intermediate with two contiguous three-coordinate Fe atoms is too large. A bridging ligand is retained throughout the proposed mechanisms. Intermediates with Fe-C(O)-S/Se-Fe cycles and with SCO/SeCO C-bound to Fe are predicted. The energetics of the reaction trajectories show them to be feasible and easily reversible, consistent with experiment. Alternative mechanisms involving intramolecular differential rotatory rearrangements of the cluster to scramble the Se bridges are also examined, and shown to be very unlikely. The implications of these new facets of the reactivity of the FeMo-co cluster are discussed: it is considered that they are unlikely to be part of the mechanism of the physiological reactions of nitrogenase.

EXAFS and NRVS reveal a conformational distortion of the FeMo-cofactor in the MoFe nitrogenase propargyl alcohol complex

We have used EXAFS and NRVS spectroscopies to examine the structural changes in the FeMo-cofactor active site of the α-70(Ala) variant of Azotobacter vinelandii nitrogenase on binding and reduction of propargyl alcohol (PA). The Mo K-edge near-edge and EXAFS spectra are very similar in the presence and absence of PA, suggesting PA does not bind at Mo. By contrast, Fe EXAFS spectra show a clear and reproducible change in the long Fe-Fe interaction at ~3.7 Å on PA binding with the apparent appearance of a new Fe-Fe interaction at 3.99 Å. An analogous change in the long Mo-Fe 5.1 Å interaction is not seen. The NRVS spectra exclude the possibility of large-scale structural change of the FeMo-cofactor involving breaking the μ(2) Fe-S-Fe bonds of the Fe(6)S(9)X core. The simplest chemically consistent structural change is that the bound form of PA is coordinated at Fe atoms (Fe6 or Fe7) adjacent to the Mo terminus, with a concomitant movement of the Fe away from the central atom X and along the Fe-X bond by about 0.35 Å. This study comprises the first experimental evidence of the conformational changes of the FeMo-cofactor active site on binding a substrate or product.

Ligand-bound S = 1/2 FeMo-cofactor of nitrogenase: hyperfine interaction analysis and implication for the central ligand X identity

Broken symmetry density functional theory (BS-DFT) has been used to address the hyperfine parameters of the single atom ligand X, proposed to be coordinated by six iron ions in the center of the paramagnetic FeMo-cofactor (FeMoco) of nitrogenase. Using the X = N alternative, we recently found that any hyperfine signal from X would be small (calculated A(iso)(X = (14)N) = 0.3 MHz) due to both structural and electronic symmetry properties of the [Mo-7Fe-9S- X] FeMoco core in its resting S = 3/2 state. Here, we extend our BS-DFT approach to the 2e(-) reduced S = 1/2 FeMoco state. Alternative substrates coordinated to this FeMoco state effectively perturb the electronic and/or structural symmetry properties of the cofactor. Using an example of an allyl alcohol (H2C=CH-CH2-OH) product ligand, we consider three different binding modes at single iron site and three different BS-DFT spin state structures and show that this binding would enhance the key hyperfine signal A(iso)(X) by at least 1 order of magnitude (3.8 < or = A(iso)(X = (14)N) < or = 14.7 MHz), and this result should not depend strongly on the exact identity of X (nitrogen, carbon, or oxygen). The interstitial atom, when the nucleus has a nonzero magnetic moment, should therefore be observable by ESR methods for some ligand-bound FeMoco states. In addition, our results illustrate structural details and likely spin-coupling patterns for models for early intermediates in the catalytic cycle.

Modeling hydrogen evolution from the Fe4S4and Fe8S9X (X = N, C) clusters. Can a FeS high-spin cluster serve as a surrogate for the FeMo cofactor?

A high-spin model of nitrogenase with a Fe(8)S(9)X(+) cluster (X = nitrogen or carbon) is used to test a mechanism for molecular hydrogen production, which is known to accompany ammonia production. The reaction proceeds with a series of protonation-reduction (PR) steps which are considered to be spontaneous if the calculated hydrogen-cluster bond energy exceeds 35-40 kcal/mol. The novel features of this mechanism include the opening of the cluster when one of the bridging sulfides undergoes two PR steps and the direct participation of the central atom when it undergoes a PR step. After the sixth PR step, a cluster is formed which has a low barrier for loss of molecular hydrogen in an exothermic reaction step. The central atom (nitrogen or carbon) has only a minor effect on the reaction steps.

Modeling the nitrogenase FeMo cofactor with high-spin Fe8 S9 X+ (XN, C) clusters. Is the first step for N2 reduction to NH3 a concerted dihydrogen transfer?

A high-spin Fe(8)S(9)X(+) (X=N, C) cluster is used to model the reduction of molecular nitrogen to ammonia by the nitrogenase FeMo cofactor at the B3LYP/6-311G(d,p)/ECP(Fe,SDD) level of theory. A total of seventy-three structures were optimized (including three transition state optimizations) to explore the structure and energetic of N(2), C(2)H(2), and CO coordination to the Fe(8)S(9)X(+) cluster. After three protonation-reduction (PR) steps (modeled by addition of hydrogen atoms), N(2), C(2)H(2), and CO are predicted to bind to a Fe atom in the exo (cage does not open) position with binding energies of 7.6, 14.7, and 11.7 kcal/mol. With additional PR steps the coordination number of the core nitrogen atom is reduced from six to five and the bridging thiol group becomes a terminal SH(2) group. The fifth and sixth PR steps occur on the core nitrogen and the open Fe site. Coordination of N(2) is enhanced after six PR steps to give an intermediate ideally suited for a concerted dihydrogen transfer from the Fe and core nitrogen atoms to the coordinated N(2). The identity of the central atom (nitrogen or carbon) has only a minor effect on the reaction steps.

Ammonia Production at the FeMo Cofactor of Nitrogenase: Results from Density Functional Theory

Biological nitrogen fixation has been investigated beginning with the monoprotonated dinitrogen bound to the FeMo cofactor of nitrogenase up to the formation of the two ammonia molecules. The energy differences of the relevant intermediates, the reaction barriers, and potentially relevant side branches are presented. During the catalytic conversion, nitrogen bridges two Fe atoms of the central cage, replacing a sulfur bridge present before dinitrogen binds to the cofactor. A transformation from cis- to trans-diazene has been found. The strongly exothermic cleavage of the dinitrogen bond takes place, while the Fe atoms are bridged by a single nitrogen atom. The dissociation of the second ammonia from the cofactor is facilitated by the closing of the sulfur bridge following an intramolecular proton transfer. This closes the catalytic cycle.

Intrinsic Dinitrogen Activation at Bare Metal Atoms

There is ongoing interest in metal complexes which bind dinitrogen and facilitate either its reduction or oxidation under mild conditions. In nature, the enzyme nitrogenase catalyzes this process, and dinitrogen fixation occurs under mild and ambient conditions at a metal-sulfur cluster in the center of the MoFe protein, but the mechanism of this process remains largely unknown. In the last few years, new important discoveries have been made in this field. In this review are discussed recent findings on the interaction of N(2) with metal atoms and metal-atom dimers from all groups of the periodic table as provided by gas-phase as well as matrix-isolation experiments. Intrinsic dinitrogen activation at such bare metal atoms is then related to corresponding processes at complexes, clusters, and surfaces.

Exploring new frontiers of nitrogenase structure and mechanism

The mechanism of the complex enzyme nitrogenase has long been one of the most challenging problems in bioinorganic chemistry. The complexity of the metal centers of nitrogenase has stretched the boundaries of biochemical, physical and computational tools for providing insights into its structure and chemical function. Recently, there have been several key advances in crystallography and spectroscopy that have impacted the way the nitrogenase mechanism is approached. These advances have opened new frontiers in nitrogenase research, which has started to reveal novel details about the molecular structure, substrate binding and reduction. Here, we discuss these recent advances and their implications on the future of nitrogenase research.

Nitrogen Fixation under Mild Ambient Conditions: Part I—The Initial Dissociation/Association Step at Molybdenum Triamidoamine Complexes

In several recent studies Schrock and collaborators demonstrated for the first time how molecular dinitrogen can be catalytically transformed under mild and ambient conditions to ammonia by a molybdenum triamidoamine complex. In this work, we investigate the geometrical and electronic structures involved in this process of dinitrogen activation with quantum chemical methods. Density functional theory (DFT) has been employed to calculate the coordination energies of ammonia and dinitrogen relevant for the dissociation/association step in which ammonia is substituted by dinitrogen. In the DFT calculations the triamidoamine chelate ligand has been modeled by a systematic hierarchy of increasingly complex substituents at the amide nitrogen atoms. The most complex ligand considered is an experimentally known ligand with an HMT = 3,5-(2,4,6-Me3C6H2)2C6H3 substituent. Several assumptions by Schrock and collaborators on key reaction steps are confirmed by our calculations. Additional information is provided on many species not yet observed experimentally. Particular attention is paid to the role of the charge of the complexes. The investigation demonstrates that dinitrogen coordination is enhanced for the negatively charged metal fragment, that is, coordination is more favorable for the anionic metal fragment than for the neutral species. Coordination of N2 is least favorable for the cationic metal fragment. Furthermore, ammonia abstraction from the cationic complex is energetically unfavorable, while NH3 abstraction is less difficult from the neutral and easily feasible from the anionic low-spin complex.

Activation and protonation of dinitrogen at the FeMo cofactor of nitrogenase

The protonation of N2 bound to the active center of nitrogenase has been investigated using state-of-the-art density-functional theory calculations. Dinitrogen in the bridging mode is activated by forming two bonds to Fe sites, which results in a reduction of the energy for the first hydrogen transfer by 123 kJ/mol. The axial binding mode with open sulfur bridge is less reactive by 30 kJ/mol and the energetic ordering of the axial and bridged binding modes is reversed in favor of the bridging dinitrogen during the first protonation. Protonation of the central ligand is thermodynamically favorable but kinetically hindered. If the central ligand is protonated, the proton is transferred to dinitrogen following the second protonation. Protonation of dinitrogen at the Mo site does not lead to low-energy intermediates.