The discovery of Mo(III) in FeMoco: reuniting enzyme and model chemistry

Biological nitrogen fixation is enabled by molybdenum-dependent nitrogenase enzymes, which effect the reduction of dinitrogen to ammonia using an Fe7MoS9C active site, referred to as the iron molybdenum cofactor or FeMoco. In this mini-review, we summarize the current understanding of the molecular and electronic structure of FeMoco. The advances in our understanding of the active site structure are placed in context with the parallel evolution of synthetic model studies. The recent discovery of Mo(III) in the FeMoco active site is highlighted with an emphasis placed on the important role that model studies have played in this finding. In addition, the reactivities of synthetic models are discussed in terms of their relevance to the enzymatic system.


Introduction
Molybdenum-dependent nitrogenase is a complex enzyme that catalyzes the formation of ammonia through the electrochemical reduction of dinitrogen at ambient temperature and pressure, utilizing eight electrons, eight protons and 16 MgATP molecules [1,2]. Unraveling the structural and mechanistic details of how nature activates the strong bond of dinitrogen is of fundamental importance, potentially aiding in catalyst design.
The molybdenum-iron component (MoFe) of nitrogenase has two types of complex metal clusters, the "P-cluster" and the FeMo cofactor (FeMoco) (Fig. 1). The P-cluster is an eight-iron cluster that appears to serve as an electron-transfer site, while FeMoco is generally agreed to be the site of dinitrogen reduction. The FeMo cofactor consists of 7 irons, 1 molybdenum, 9 sulfides and an interstitial light atom that was recently identified as carbon by XES, ESEEM and high-resolution crystallography [3,4].
FeMoco has fascinated chemists for decades. This complex cofactor, whose structural details have been elusive for so long, is capable of chemistry that is still a challenge for synthetic chemists to mimic. While the molecular structure details of FeMoco are finally clear, many uncertainties remain about the S = 3/2 resting form (the dithionite-isolated form, also called E 0 ) of the enzyme. The total charge of the cofactor, the metal oxidation states and the electronic structure, all remain open questions-the answers to which are essential for any informed discussion of the molecular level mechanism of dinitrogen reduction.
Abstract Biological nitrogen fixation is enabled by molybdenum-dependent nitrogenase enzymes, which effect the reduction of dinitrogen to ammonia using an Fe 7 MoS 9 C active site, referred to as the iron molybdenum cofactor or FeMoco. In this mini-review, we summarize the current understanding of the molecular and electronic structure of FeMoco. The advances in our understanding of the active site structure are placed in context with the parallel evolution of synthetic model studies. The recent discovery of Mo(III) in the FeMoco active site is highlighted with an emphasis placed on the important role that model studies have played in this finding. In addition, the reactivities of synthetic models are discussed in terms of their relevance to the enzymatic system.
In this mini-review, we give a brief overview of the structural characterization of the enigmatic FeMoco cluster, including both molecular and electronic structure aspects. The evolution of our structural understanding of FeMoco is complemented by a review of FeMoco-inspired synthetic model chemistry. The strong synergy between model chemistry and the advances in our understanding of FeMoco active site structure are highlighted. Furthermore, the reactivity of the model complexes is discussed in terms of its potential relevance to the enzyme chemistry.

First structural models
Prior to the first crystal structure of MoFe protein, structural modeling of FeMoco (as well as the P-cluster) relied primarily on information from X-ray absorption spectroscopy (XAS) with a focus on the extended X-ray absorption fine structure (EXAFS) region to acquire basic molecular structural information. In 1978, Hodgson and coworkers [5] published the first structural models based on EXAFS data from Mo XAS experiments on MoFe protein (Fig. 2). The data indicated 3-4 sulfur atoms at a 2.36 Å distance from the molybdenum and 2-3 iron atoms at 2.72 Å. The models proposed, included a cubane-type structure and a linear Fe-Mo-Fe compound.
The initial EXAFS structural models for FeMoco served as inspiration for the synthetic Mo-Fe-S based chemistry that several research groups explored in detail for the next decades, particularly the groups of Holm [6] and Coucouvanis [7]. Holm and coworkers explored a self-assembly route to Mo-Fe-S clusters based on the observed spontaneous formation of Fe 4 S 4 cubane clusters from simple reactants such as FeCl 3 and NaSH and NaSR. Adding [MoS 4 ] 2− to the reaction mixture afforded the first Mo-Fe-S double cubanes, synthesized by the groups of Holm and Garner [8][9][10][11][12] in the period of 1978-1980. These compounds were the first heterometal weak field cubane clusters to be synthesized and sparked the beginning of nitrogenase-inspired model compound chemistry that continues to be explored today (Fig. 3).
As the EXAFS structural data on FeMoco suggested a single molybdenum atom per cluster (and a single cubane), single cubane Mo-Fe-S compounds were the next synthetic goal. Cleavage of double cubanes resulted in the first single [MoFe 3 S 4 ] cubanes reported between 1981 and 1983, which showed close similarity to the original single cubane FeMoco EXAFS proposal (Fig. 2) [13][14][15]. In 1982, however, Fe EXAFS data indicated secondshell Fe-Fe interactions [16], indicating that the cofactor was a larger cluster rather than a single cubane. The complexity of the FeMoco active site was confirmed in 1992, when the first crystal structure of MoFe protein became available from Rees et al. [17]. The 2.7 Å structure revealed a complicated MoFe 7 S 9 structure that was neither a cubane nor a double cubane, as a corner atom of each cube appeared to be missing (see Fig. 2).
The structure had surprising features: six three-coordinate iron atoms and a vacant cavity, and it was the source of confusion for many years.
For example, 57 Fe Mössbauer studies of FeMoco were interpreted in view of the trigonal Fe ions that had no synthetic analogs, and both isomer shifts and quadrupole splittings were found to be unusual and puzzling [18,19]. This finding, however, did inspire research focused on low-coordinate iron complexes as synthetic models of nitrogenase [20][21][22][23]. While the direct correlation to the resting form of nitrogenase no longer held upon discovery of the central atom, it is still speculated that 3-coordinate iron sites could have relevance to intermediates in the N 2 reduction process [23,24].
In fact, it was proposed early on that the trigonal irons in nitrogenase would be highly reactive and likely to bind substrate (even within the cavity). Multiple reaction mechanisms were proposed from quantum chemical calculations [25][26][27][28][29]. The vacant cavity was even proposed to become occupied by a hydride upon reduction [30]. The structure also further inspired Mo-Fe-S cluster chemistry and new cubane clusters void of a corner sulfide were synthesized [31][32][33].
The confusing aspects of the FeMoco structure were not resolved until 2002 [34], when a more highly resolved structure revealed an interstitial light atom that could be assigned as either an oxygen, a nitrogen or a carbon. In 2011, the light atom was confidently demonstrated to be carbon by both Fe X-ray emission [4] and X-ray crystallography techniques [3] (see JBIC review by Einsle [35]). The discovery of the interstitial atom confirmed the structure of FeMoco to be cubane-like after all (a fused double cubane structure with four-coordinate distorted tetrahedral irons), which re-established the biological relevance of the synthetic single and double Mo-Fe-S cubanes.
The synthetic Mo-Fe-S cubane structures are highly reminiscent of half of the full FeMoco cluster (see Fig. 4 Fig. 3) appears to be rather different than the molybdenum coordination [MoFe 3 S 4 ] cubanes with even more protein-like ligands have also been synthesized, including polycarboxylate ligands such as citrate and methyliminodiacetate [38,39].
Synthetically coupling a MoFe 3 S 4 cubane with a Fe 4 S 3 fragment has so far not succeeded although various monosulfido and edge-bridged Mo-Fe-S double cubanes have come close to a FeMoco-like structure (the [(Tp) 2 Mo 2 Fe 6 S 9 (SH) 2 ] 2 cluster in Fig. 3 mimics almost perfectly the P-cluster topology [40]). Heteroatoms have also been incorporated in Fe-S based clusters (see recent review by Lee et al. [41]). Tatsumi et al. [42] were able to incorporate an oxygen atom in an [Fe 8 S 6 O] cluster, which is reasonably close to mimicking the heteroatomic FeMoco environment, and Lee et al. [43,44] succeeded in creating iron-sulfur-imido complexes, [Fe 4 S 4-n N n ]. An [Fe 4 S 3 N t Bu] complex (see Figs. 3,4) is an analog of an [Fe 4 S 4 ] 2+ cluster with a corner nitrogen instead of sulfur that was found to be a very good structural model for the Fe cubane-half of FeMoco (RMSD of 0.05 Å for Fe-S-N/C) with the corner nitrogen taking the role of the interstitial atom which results in a similar distortion of the cubane as the FeMoco carbide.
Incorporating a μ 6 -carbide in a synthetic cluster remains a challenge. Nevertheless, it is clear that the Mo-Fe-S cubane structures synthesized, appear to be reasonable structural models for the molybdenum half of the FeMoco (sulfide instead of carbide at the cubane corner), while the well-known synthetic Fe-S cubane structures ([Fe 4 S 4 ] or [Fe 4 S 4-n X n ] structures) are quite reasonable models for the other half. Figure 4   An obvious next question then is whether the synthetic cubanes have similarities at the electronic structure level as well.

Electronic structure and oxidation states
The resting state of the FeMoco cluster has long been known to have a spin state of S = 3/2 [45,46]. The charge and metal oxidation states have not been confidently assigned, but there are generally three oxidation state assignments discussed in the literature: (1) 6Fe(II)1Fe(III) Mo(IV) [47]; (2) 4Fe(II)3Fe(III)Mo(IV) [18]; and (3) 2Fe(II)5Fe(III)Mo(IV) [48]. All 3 assignments are consistent with an S = 3/2 spin state and assume a closed-shell, diamagnetic, S = 0 Mo(IV) ion (see Fig. 5). A recent reassignment of the Mo oxidation state [49] is discussed in "Comparing the molybdenum oxidation state in MoFe protein vs Mo-Fe-S clusters". When written in terms of total charge on the cluster (an unambiguous way of presenting different charge models), with sulfur and carbon taken in their usual closed-shell forms, S 2− and C 4− (for electroncounting purposes only), the charges are [MoFe 7 S 9 C] 3− , [MoFe 7 S 9 C] 1− and [MoFe 7 S 9 C] 1+ , respectively. In view of biological iron-sulfur compounds being predominantly Fe(II) or Fe(III), it seems likely that the same applies to FeMoco.
The Mo(IV) assignment in FeMoco has been present in the nitrogenase literature since the 1980s. The assignment dates back to early Mo K-edge XAS studies by Hodgson states are shown. The actual physical oxidation state distribution is most likely more complicated due to electron delocalization [as seen in the Fe(2.5) pairs of the cubanes for example] et al. and 95 Mo ENDOR studies by Hoffman et al., but the assignments were made prior to the first crystal structure of FeMoco was published. Mo-S derived bond lengths from the Mo EXAFS region suggested an oxidation state of either Mo(III) or Mo(IV) [5,50]. Due to the significant core hole lifetime broadening at high energies, absorption edge positions did not give a clear indication of the Mo oxidation state. Later 95 Mo ENDOR experiments [51][52][53] indicated a small Mo hyperfine coupling in the protein and the data were interpreted as the molybdenum most plausibly being a closed-shell Mo(IV), rather than a S = 3/2 Mo(III) or S = 1/2 Mo(V). The authors acknowledged that the spin coupling schemes utilized were greatly simplified and that this assignment was tentative [52]. In 1988, a Mo L-edge XAS study [54] assigned the Mo atom in FeMoco as Mo(IV), when taking the ENDOR data into account, although a comparison of Mo-Fe-S based double cubanes and isolated FeMoco showed a strong similarity of the L 3 edges. However, a quantitative analysis of Mo L-edge data was prohibitive, due to both the proximity of the sulfur K-edge and the inadequacy of the available theoretical tools at that time. It seems that the initial Mo(IV) assignment of FeMoco influenced several studies on the electronic structure of FeMoco, both using experimental magnetic Mössbauer [18] and theoretical methods [48,[55][56][57][58].
The electronic structure of FeMoco is revisited in "Comparing the molybdenum oxidation state in MoFe protein vs Mo-Fe-S clusters".
Mo-Fe-S compounds with [MoFe 3 S 4 ] n cores have been prepared with oxidation states of n = 5+, 4+, 3+, 2+ and 1+, either as single or double cubane compounds (see Fig. 5). The [MoFe 3 S 4 ] 3+ oxidation level is the most easily accessible oxidation state and such cubanes are always found to have a total spin of S = 3/2 [59], while [MoFe 3 S 4 ] 2+ cubanes (accessible by chemical reduction or a separate synthetic pathway) are S = 2 [60,61]. Mössbauer spectroscopy has primarily been used to assign oxidation states in these cubanes, where the Mo oxidation state is assigned based on the individual Fe assignments, or sometimes only the average Fe oxidation state assignment using an empirical equation relating Fe oxidation state, isomer shift, and the total charge [6] [63]. Due to the potential ligand non-innocence, it is not clear whether their electronic structure is analogous to the more reduced cubanes [6].
Finally, we should mention that [MoFe 3 S 4 ] 0 cores have also been prepared [64][65][66]. We mention them separately, as these compounds are not directly related to the more oxidized cores; they have only been prepared with CO ligation on the molybdenum and formally contain a Mo(0 It is worth noting that analogous heterometal cubanes with vanadium and tungsten have also been synthesized [6]. The [VFe 3 S 4 ] 2+ core is isoelectronic to [MoFe 3 S 4 ] 3+ , has the same spin S = 3/2 ground state and is approximately isostructural, but has been found to have a more reduced Fe-part [67], suggesting an oxidation state assignment of V(III)Fe(II) 2 Fe(III). Other [VFe 3 S 4 ] n cubanes have also been synthesized and characterized [68] and found to have Fe-based reductions with the vanadium atom assigned as V(III) (d 2 configuration). These compounds are potentially relevant structural models for the iron-vanadium cofactor (FeVco) that has been suggested to be analogous to FeMoco [69], albeit this still awaits crystallographic verification.
Besides standard Mössbauer spectroscopy as a means to gain insight into the Fe electronic structure of these cubanes, various magnetic spectroscopies have been employed. EPR spectroscopy and magnetic susceptibility measurements were used to determine a ground spin state of S = 3/2 in the [MoFe 3 S 4 ] 3+ cubanes [59]. Both [MoFe 3 S 4 ] 3+ and MoFe protein were in fact found to give rise to very similar EPR spectra [60]. A study by Mascharak et al. [60] [60], two magnetic subsites were revealed, again with an intensity ratio of 2:1. Analyzing the splittings as a function of field strength was indicative of antiparallel spin coupling of the Fe ions in the cubane with the more intense subsite having a negative hyperfine field, hence being the majority spin. The majority spin species can nowadays plausibly be interpreted as a mixed-valence delocalized Fe(2.5)-Fe(2.5) pair, with the minority spin being a single Fe(III) ion. We note that while the authors of the magnetic Mössbauer studies also considered a Mo(IV)1Fe(III)2Fe(II) interpretation as a possibility, this seems unlikely in view of the cases where well-resolved Mössbauer spectra of [MoFe 3 S 4 ] 3+ cubanes show the more intense Fe species with a higher isomer shift. Further, such an assignment is also inconsistent with theoretical calculations [70,71], as noted below.
Mascharak et al. highlighted in their study the similarities in electronic structure, as determined by the magnetic spectroscopies, of synthetic cubanes vs. MoFe protein, ranging from a similar electron delocalization, antiparallel spin coupling, hyperfine couplings and similar EPR spectra.
Cook and Karplus [70,71] were the first to attempt a detailed electronic structure analysis of a [MoFe 3 S 4 ] 3+ cubane by theoretical calculations. Broken-symmetry MS-Xα calculations were performed which resulted in an electronic structure picture consisting of antiferromagnetically coupled Fe atoms, mixed-valence delocalization and a molybdenum(III) assignment with an unusual (↑↑↓) configuration and low spin population.
As the discussion above highlights, the comparison of FeMoco and Mo-Fe-S cubanes goes beyond mere molecular structure similarity. Mixed-valence delocalization, as well as antiferromagnetic spin coupling of Fe, is apparent in the electronic structure for both FeMoco and the Mo-Fe-S cubanes. Further, the EPR spectra and Mössbauer properties are similar. Yet, the oxidation state of molybdenum is almost exclusively assigned as Mo(III) in the Mo-Fe-S cubane literature, regardless of cubane oxidation state (indicating redox processes taking place at Fe), while the molybdenum in FeMoco has (until very recently) always been assigned as a closed-shell S = 0 Mo(IV) in the nitrogenase literature. This apparent discrepancy prompted us to probe directly the molybdenum oxidation state in both MoFe protein and Mo-Fe-S clusters, using a combined Mo XAS spectroscopy and theoretical approach.

Comparing the molybdenum oxidation state in MoFe protein vs Mo-Fe-S clusters
In a recent article by some of us [49], we revisited the molybdenum oxidation state assignment in FeMoco by a combined experimental and theoretical approach. By performing high-energy resolution fluorescence detected (HERFD) Mo X-ray absorption spectroscopy on MoFe protein and selected model compounds, we obtained higher resolution Mo XAS spectra than had been possible before. Two of the model compounds measured were Mo-Fe-S cubane compounds, originally described by Holm et al. One was a [MoFe 3 S 4 ] 3+ cubane with a trispyrazolylborate (Tp) ligand on the molybdenum and the other a ([MoFe 3 S 4 ] 3+ ) 2 double cubane with Tp on each molybdenum.
Since the molybdenum coordination environment in MoFe protein and the two cubane complexes is very similar, the absorption edge position should be a good indicator of a difference in oxidation state [72]. As seen in Fig. 6, both the pre-edge and the main absorption edge are very similar in MoFe protein and the 2 cubane complexes, with MoFe protein having a rising edge if anything, to slightly lower energy than either cubane model. This is in contrast to the absorption edges of the Mo(V) compounds, which have absorption edges at much higher energies. We note that an appropriate Mo(IV) reference model complex was not available for the initial study. However, we note that the Mo(IV) possibility for FeMoco can be ruled out due to the fact that the rising edge of MoFe protein appears to slightly lower energy than the Mo(III) cubane references. A Mo(IV) complex should have a rising edge between that of the Mo(III) and Mo(V) model complexes. Hence based on these data, a Mo(III) assignment was found to be most appropriate for FeMoco. Further, we note that in the original study, careful radiation damage studies were carried out, and the published data correspond to undamaged MoFe protein.
Due to the similarities in the Mo HERFD XAS data for the [MoFe 3 S 4 ] 3+ cubanes and FeMoco, a Mo(III) assignment for FeMoco appears most appropriate based only on empirical considerations. However, we note that, in the initial study, we went beyond the simple fingerprint approach and computed the Mo XAS pre-edges using time-dependent density functional theory (TDDFT) from broken-symmetry SCF solutions of model compounds and a 225-atom FeMoco model. Computed pre-edges of model compounds and our cofactor model were found to be in excellent agreement with the experimental pre-edges and a rigorous orbital analysis of broken-symmetry SCF solutions of cubane model compounds and FeMoco revealed a d 3 configuration at the molybdenum, consistent with Mo(III) for both cubanes and FeMoco. Curiously, we found the molybdenum configuration to have an unusual low spin configuration (↑↑↓) (apparently analogous to the one found by Cook and Karplus earlier for a similar model system [71]), that we have referred to as a non-Hund configuration. This configuration, combined with covalency effects, appears to be the reason for a small spin population at molybdenum (ranging from −0.1 to −0.6, depending on the DFT method) compared to an idealized value of 3.0 for an ionic Mo(III). This unusual feature may have been missed in previous DFT studies of FeMoco and appears to have not been explored further for Mo-Fe-S systems since the study by Cook and Karplus. Due to limitations of a single-reference quantum chemistry approach (whether DFT or MS-Xα), it is unfortunately not possible to tell whether the non-Hund configuration in FeMoco or Mo-Fe-S cubanes represents a genuine excited spin state of the molybdenum or whether the configuration arises due to spin-canting effects in the cubane system. A multiconfigurational ab initio approach may help to shed some light on the matter. Such approaches have recently been applied to the exchange-coupled manganese-based oxygen-evolving complex [73] as well as Fe 4 S 4 cubanes [74]. It is clear that a magnetic spectroscopic approach will be needed for direct experimental evidence of the spin structure at the molybdenum. What the DFT calculations reveal, however, is that the unusual molybdenum configuration can only be stabilized due to spin coupling with the iron atoms, and a localized orbital analysis suggests that the molybdenum electrons reside in partially bonding orbitals with Fe, indicative of a double-exchange effect between Mo and Fe. A spin-coupled Mo(III) is probably the best description of molybdenum in FeMoco for now, while the details of the Mo-Fe interaction are not understood.
The Mo HERFD XAS study was unfortunately not able to establish the total charge on the FeMoco cluster; however, it was established that for all 3 charges considered for FeMoco, the unusual Mo(III) configuration remained. For the [MoFe 7 S 9 C] 1− model, the electronic structure in the Mo-cubane fragment was found to be analogous to the electronic structure of the [MoFe 3 S 4 ] 3+ cubane, in either the [(Tp)MoFe 3 S 4 Cl 3 ] 1− cluster or the double cubane [(Tp) MoFe 3 S 4 Cl 3 ] 2 (μ-S)]). This is illustrated in Fig. 7, which shows a preliminary generalized electronic structure diagram for a Mo-Fe-S cubane, with either a sulfide or carbide corner atom: i.e., [MoFe 3 S 3 C] 1+ or [MoFe 3 S 4 ] 3+ . The two cubanes share the following features: the unusual non-Hund configuration at molybdenum, a mixed-valence Fe(2.5)-Fe(2.5) pair and an antiferromagnetically coupled Fe(III). Computed Mössbauer parameters of [(Tp) MoFe 3 S 4 Cl 3 ] 1− are in good agreement with experimental values and thus support the basic features of this diagram [49].
Finally, a recent Mo L-edge study by us reexamined previous Mo L-edge data [54] on isolated FeMoco and double cubane compounds and presented new data on MoFe protein [75]. This study provides further support for the Mo(III) assignment in FeMoco; the Mo(III) assignment holds for both the NMF-isolated cofactor and the proteinbound cofactor.

Reactivity of structural model compounds: relation to biological reactivity
The Mo-Fe-S cubane compounds were synthesized to serve both as biomimetic structural models of FeMoco and as part of exploratory synthetic chemistry. It became clear, however, that a subset of these models also exhibited interesting reactivity, possibly relevant to the reactivity of FeMoco.
Nitrogenase is known to act as a proton-reducing catalyst when no N 2 substrate is available and this seems to be nature's way of terminating the electron-transfer pathway when no substrate is present. What is also known, but not understood, is the obligatory H 2 formation accompanying NH 3 formation. Dean, Hoffman, Seefeldt et al. [76,77] have recently proposed reductive elimination of hydrides occurring at the FeMoco as an explanation of this phenomenon.
Already in 1980 [78], it was found that some reduced Mo-Fe-S cubanes were capable of H 2 evolution when protonated, and in 1983, a very detailed kinetic analysis was performed [79]. The double cubane cluster [Mo 2 Fe 6 S 8 (SPh) 9 ] 5− was found to evolve dihydrogen slowly from protonation agents like PhSH and Et 3 NH + with a yield of 100 % (based on electron count). This reactivity was connected to the property of the double cubane being able to undergo 2 low-potential 1-electron reductions [potentials given vs. saturated calomel electrode (SCE)]: [Mo 2 Fe 6 S 8 (SR) 9 ] n is thus a rare example of a single molecule that can act as a carrier of 2 electrons to be delivered to a substrate (here proton). H 2 evolution from a [Fe 4 S 4 ] cubane compound, [Fe 4 S 4 (SPh) 4 ] 3− , was also found, albeit with much lower yields (11-31 %). This was attributed to the reaction of 2 cubanes per H 2 formed, as [Fe 4 S 4 (SPh) 4 ] 3− is only a 1-electron carrier [80].
A detailed kinetic study [79] of [Mo 2 Fe 6 S 8 (SPh) 9 ] n gave experimental rate laws and proposed a preliminary Since FeMoco can also be described as a double cubane (although an unusual carbon-ligated one, [MoFe 3 S 3 ]-C-[Fe 4 S 3 ]), it seems plausible that this is an important property of the cofactor, allowing it to perform multi-electron chemistry (for example H 2 evolution, and N 2 reduction).
The H 2 evolution chemistry of synthetic Mo-Fe-S cubanes has not been explored much more since then, although we note in this context the isolated FeMo cofactor studies by Pickett et al. [80][81][82]. These studies demonstrated the H 2 evolution capabilities of the NMF-isolated cofactor (85 % conversion of electrons) and proposed hydride intermediates. Also noteworthy is recent work of Kanatzidis et al. [83] where [Mo 2 Fe 6 S 8 X 9 ] double cubanes were incorporated into a chalcogenide framework. Utilizing a [Ru(bpy) 3 ] 2+ photosensitizer, the chalcogel was capable of photocatalytic H 2 production.
The catalytic 2e − reduction of hydrazine, N 2 H 4 , to ammonia was accomplished by different [MoFe 3 S 4 ] 3+ cubanes in very detailed studies where [Fe 4 S 4 Cl 4 ] 2− clusters were found to be inactive. Through explorations of various ligands to the molybdenum atom, the most effective catalysts were found to be compounds with polycarboxylate ligands (such as citrate) on molybdenum. It was hypothesized that this ligand preference could derive from citrate protonation opening a free coordination site for hydrazine and the polycarboxylate acting as a proton shuttle. This observation is of potential relevance to the enzymatic system as molybdenum-bound homocitrate in FeMoco has been found to be crucial for dinitrogen reduction. Studies where homocitrate was replaced with other polycarboxylate ligands dramatically reduced FeMoco activity [91,92].
In the [MoFe 3 S 4 ] 3+ cubane studies, strong binding ligands at molybdenum, like PEt 3 were found to inhibit the reaction, presumably by blocking all six coordination sites, strongly suggesting molybdenum as the site of binding for hydrazine. Phenylhydrazine ligated to the molybdenum of a [MoFe 3 S 4 ] 3+ was also found by X-ray crystallography [88].
CO was found to inhibit hydrazine reduction, although CO has only been found to bind to the 1-electron reduced cubane [MoFe 3 S 4 ] 2+ , with the site of binding tentatively assigned to molybdenum [15]. In this context, we note that very recently a crystal structure of CO-inhibited MoFe protein was published. The structure shows CO replacing a sulfur-bridge in FeMoco, bound to 2 Fe atoms, i.e., a rather different type of inhibition than proposed for synthetic Mo-Fe-S cubanes.
A later study by Coucouvanis et al. [90] demonstrated the catalytic 4e − reduction of cis-dimethyldiazene by [MoFe 3 S 4 ] 3+ to methylamine, the first example of reducing a N=N bond. Again strong evidence was presented for the binding and activation of the substrate at molybdenum by phosphine inhibition studies. [Fe 4 S 4 Cl 4 ] 2− showed essentially no reactivity. The sole product, methylamine, bound at the molybdenum, was found by X-ray crystallography. A difference between the synthetic system and the enzyme was also revealed as the enzyme reduces cis-dimethyldiazene to a mixture of products: ammonia, methane and methylamine.
It should be made clear that [MoFe 3 S 4 ] (or [VFe 3 S 4 ]) cores have not been found capable of 6e − /8e − N 2 reduction, which obviously is a much more difficult reaction. Coucouvanis has suggested that the catalytic activity of [MoFe 3 S 4 ] cores could be an indicator of the direct involvement of the molybdenum atom in the mechanism at the hydrazine level, by substrate migration [93]. It is possible that activation of N 2 by these types of clusters in fact requires a large cluster such as FeMoco (i.e., a double cubane fused by an interstitial carbide) not present in any of the synthetic single or double cubanes discussed here.
Functional models for N 2 reduction have a longer history than the structural Mo-Fe-S systems. The first nitrogen complex reported by Allen and Senoff in 1965, [Ru(NH 3 ) 5 (N 2 )] + , was prepared by treating RuCl 3 (H 2 O) 3 with excess hydrazine; hydrazine was the reducing agent and (presumably) the source of ammonia and nitrogen through disproportionation [94]. This discovery initiated a search not only for more nitrogen complexes, but in particular nitrogen complexes that would generate ammonia under some set of conditions, preferably catalytically. Among the groups active in the search included groups based in England, Japan, the United States, and the Soviet Union [95][96][97][98][99][100][101]. The first dinitrogen complexes of molybdenum (Mo(0)) were prepared by Hidai and his group in 1969. The first (stoichiometric) generation of ammonia from Mo(0) and W(0) complexes through addition of acid was published by Chatt and Richards in 1975. The first catalytic reduction of dinitrogen published by Shilov and his group employed a source of molybdenum and a strong reducing agent in methanol; however, the primary reduction product in this system is hydrazine, which is then disproportionated to nitrogen and ammonia. Typically a 10:1 mixture of hydrazine and ammonia is reportedly formed.
A breakthrough molybdenum-based catalyst for N 2 reduction was published by Schrock et al. in 2003 [102-104]. The relevant Mo complexes contain the [HIPTN 3 N] 3− ligand, where HIPT = 3,5-(2,4,6-i-Pr 3 C 6 H 2 ) 2 C 6 H 3 (HexaI-soPropylTerphenyl; Fig. 9), a ligand that was designed to prevent formation of relatively stable and unreactive bimetallic (Mo-N=N-Mo) complexes, maximize steric protection of a metal coordination site in a monometallic species, and provide increased solubility in nonpolar solvents. Eight of the proposed intermediates in a hypothetical "Chatt-like" reduction of dinitrogen were prepared and characterized. These include paramagnetic The proposed catalytic reduction mechanism has been fully confirmed through computational mechanistic studies that involve the full ligand [109].
A second Mo(0)-based bimetallic catalyst was found by Nishibayashi and coworkers to reduce dinitrogen to ammonia with Co(η 5 -C 5 H 5 ) 2 as the reducing agent and lutidinium triflate as the proton source to give a maximum of ~12 equivalents of ammonia per Mo [110]. No catalytically viable intermediates have been isolated and the detailed mechanism of the reduction remains obscure. It seems unlikely that Mo(0) is the lowest oxidation state, since Co(η 5 -C 5 H 5 ) 2 is not a strong enough reducing agent to form Mo(0).
Fe-based catalysts for N 2 reduction were recently reported by Peters et al. [111,112]. The compounds are phosphine-supported trigonal geometries with axial coordination to iron by boron, silicon or carbon. The compounds bind N 2 in the other axial position trans to B, Si or C; the B-and C-containing systems catalyze reduction to NH 3 . The conditions are relatively severe with the reducing agent being KC 8 and the proton source being H(Et 2 O) 2 BAr′ 4 in THF at −78 °C. Approximately 7 equivalents of ammonia and ~15 equivalents of hydrogen are formed per Fe in the system in which the ligand contains boron trans to the nitrogen-binding site. It has been suggested that the atom trans to nitrogen in the TBP complex provides a "flexible" interaction with the metal trans to the N 2 binding site [113], and therefore is crucial to catalytic activity. However, the ligand in which a carbon atom is trans to the nitrogen-binding site gives only ~5 equivalents of ammonia per Fe, i.e., less than the system that contains boron.
None of the functional Fe or Mo systems has any structural or mechanistic feature that could be said to be "biomimetic" beyond the obvious, i.e., the presence of Fe or Mo. The only distant possibility is the trigonal arrangement of sulfur atoms around both Fe and Mo in FeMoco. Molybdenum and iron-based homogenous systems for N 2 reduction were recently reviewed by Holland et al. [114].

Uncovering the role of the molybdenum
The molecular structure of resting state FeMoco can be considered complete and recent studies in our group are moving toward a clear assignment of the charge. Combining this information with the Mo(III) assignment and electronic structure calculations gives us a simplified, but useful picture of the metal oxidation states in the cofactor. Perhaps the most important result of recent studies is that the cofactor can now be connected with the rich structural information of other Mo-Fe-S clusters. Due to complexity of the cofactor, however, the structural information that we now possess, provides only the first step toward understanding how FeMoco operates.
Molybdenum is clearly an essential part of FeMoco as is known from enzymes formed in molybdenum-deficient environments. What is still not understood is the role of the molybdenum in the mechanism of nitrogen reduction. It is important to reiterate that there is still no direct evidence for the binding site of N 2 . Current biochemical data favor one of the irons in the Mo cubane, Fe6 (X-ray structure labeling), as the binding site, mainly due to loss of activity upon increasing the bulkiness (mutating α-70 Val into Ile) of the residue directly above Fe6 [115,116]. If molybdenum is not the binding site then the atom may play a role in tuning the redox potential of the cofactor as a whole (or even a specific iron atom). Alternatively the role of the Mo may be to modify the electronic structure of the cofactor to facilitate binding or protonation of substrate. As suggested by studies by Coucouvanis et al. on the catalytic activity of [MoFe 3 S 4 ] 3+ cubanes for hydrazine and cis-dimethyldiazene reduction, the molybdenum might also play a direct role in later stages of nitrogen reduction.
The observation that molybdenum (or vanadium) appears to never change oxidation state upon reduction of Mo-Fe-S systems, as revealed by the electron-transfer chemistry of synthetic cubanes and also suggested for FeMoco, may even be an important property of the heterometal in the cluster. This would be in sharp contrast to the Schrock/Chatt cycles for nitrogen reduction, however, where molybdenum-based redox is crucial to catalytic activity.
The next challenges in FeMoco research will be to gain a better understanding of the electronic structure of FeMoco and particularly to understand its relation to reactivity. This will involve a characterization of the redox states of FeMoco (as FeMoco requires at least 3 reduction events prior to N 2 binding), where Fe-based hydrides have been proposed to be an important feature of the redox process [76,77], as well as gathering direct experimental evidence of the N 2 binding site and the nature of binding and activation. Overcoming these challenges will signify a real step toward an understanding of biological dinitrogen reduction.