How N2 Might Be Activated by the FeMo-Cofactor in Nitrogenase

Molybdenum-Containing Metalloenzymes and Synthetic Catalysts for Conversion of Small Molecules

The energy deficiency and environmental problems have motivated researchers to develop energy conversion systems into a sustainable pathway, and the development of catalysts holds the center of the research endeavors. Natural catalysts such as metalloenzymes have maintained energy cycles on Earth, thus proving themselves the optimal catalysts. In the previous research results, the structural and functional analogs of enzymes and nano-sized electrocatalysts have shown promising activities in energy conversion reactions. Mo ion plays essential roles in natural and artificial catalysts, and the unique electrochemical properties render its versatile utilization as an electrocatalyst. In this review paper, we show the current understandings of the Mo-enzyme active sites and the recent advances in the synthesis of Mo-catalysts aiming for high-performing catalysts.

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 Fe₇MoS₉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.

Electron paramagnetic resonance spectroscopy

EPR spectroscopy has been an important tool in nitrogenase research for the last 50 years. The three metalloclusters in nitrogenase, the Fe protein [4Fe-4S] cluster, and the MoFe protein P-cluster, and FeMo-cofactor, all have EPR spectra when poised in the appropriate paramagnetic states. EPR spectroscopy can probe changes in the electronic properties of each metal cluster, such as when substrates bind, and can provide a definitive method for observing changes in the redox states of the clusters. In this chapter, the methods for analysis of the three metal clusters of nitrogenase by EPR spectroscopy are described, along with methods for trapping substrate-derived intermediates on the active site that are amenable to characterization by EPR and other magnetic resonance spectroscopy techniques.

Mixed-valence dinitrogen-bridged Fe(0)/Fe(II) complex

The reactions of a dinitrogen-bridged Fe(II)/Fe(II) complex [(FeH(PP(3)))(2)(μ-N(2))](2+) (3) (PP(3) = P(CH(2)CH(2)PMe(2))(3)) with base were investigated using (15)N labeling techniques to enhance characterization. In the presence of base, 3 is initially deprotonated to the Fe(II)/Fe(0) dinitrogen-bridged complex [(FeH(PP(3)))(μ-N(2))(Fe(PP(3)))](+) (4) and then to the symmetrical Fe(0)/Fe(0) dinitrogen-bridged complex (Fe(PP(3)))(2)(μ-N(2)) (5). [(FeH(PP(3)))(μ-N(2))(Fe(PP(3)))](+) (4) exhibits unusual long-range (31)P-(31)P NMR coupling through the bridging dinitrogen ligand from the phosphines at the Fe(0) center and those at the Fe(II) center. Reaction of 4 with base under an atmosphere of argon resulted in the known dinitrogen Fe(0) complex Fe(N(2))(PP(3)) (6) and a solvent C-H activation product. Complexes 3, 4, and 5 were fully characterized by multinuclear NMR spectroscopy, and complexes 3 and 4 by X-ray crystallography.

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.

Reaction of H2 with a Binuclear Zirconium Dinitrogen Complex - Evaluation of Theoretical Models and Hybrid Approaches

Molecular orbital and hybrid ONIOM (both IMOMO and IMOMM) calculations have been carried out on the important reaction of H2 with a binuclear zirconium dinitrogen complex to test the efficacy of several structural models of the ancillary ligand. The complete experimental ligand, PhP(CH2SiMe2NSiMe2CH2)2PPh, in the zirconium complex has been treated at the IMOMM level, while two smaller approximations of the ligand, HP(CH2SiH2NSiH2CH2)2PH and (PH3)2(NH2)2, have received the full molecular orbital treatment. The mechanism of dihydrogen addition has been compared with our earlier study (Basch, Musaev, and Morokuma J. Am. Chem. Soc. 1999, 121, 5754-5761). We find that the substituent effects do cause some small changes in both the structures of the complexes studied and the activation energies of the transition structures. However for the most part the potential energy profiles are very similar to our earlier study and lend support to our use of simple theoretical models to represent moderately large experimental structures.

An Electronic Perspective on the Reduction of an N?N Double Bond at a Conserved Dimolybdenum Core

Density functional theory has been used to assess the role of the bimetallic core in supporting reductive cleavage of the N=N double bond in [Cp2Mo2(mu-SMe)3(mu-eta1:eta1-HN=NPh)]+. The HOMO of the complex, the Mo-Mo delta orbital, plays a key role as a source of high-energy electrons, available for transfer into the vacant orbitals of the N=N unit. As a result, the metal centres cycle between the Mo(III) and Mo(IV) oxidation states. The symmetry of the Mo-Mo delta "buffer" orbital has a profound influence on the reaction pathway, because significant overlap with the redox-active orbital on the N=N unit (pi* or sigma*) is required for efficient electron transfer. The orthogonality of the Mo-Mo delta and N-N sigma* orbitals in the eta1:eta1 coordination mode ensures that electron transfer into the N-N sigma bond is effectively blocked, and a rate-limiting eta1:eta1-->eta1 rearrangement is a necessary precursor to cleavage of the bond.

Chemical activity of the nitrogenase FeMo cofactor with a central nitrogen ligand: density functional study

We investigate the chemical consequences of a central ligand in the nitrogenase FeMo cofactor using density functional calculations. Several studies have shown that the central ligand most probably is a nitrogen atom, but the consequences for the chemical reactivity of the cofactor are unknown. We investigate several possible routes for insertion of the central nitrogen ligand and conclude that all routes involve barriers and intermediate states, which are inaccessible at ambient conditions. On this basis we suggest that the central nitrogen ligand is present at all times during the reaction. Furthermore, we investigate how the FeMoco with the central ligand can interact with N(2) and reduce it.

An Evaluation by Density Functional Theory of M−M Interactions in Organometallic Clusters with the [Fe3MoS3]2+ Cores

Density functional theory calculations were carried out on the structurally characterized [(Cl(4)-cat)Mo(py)Fe(3)S(3) (CO)(4)(P(n)Pr(3))(3)], A, and (Cl(4)-cat)Mo(py)Fe(3)S(3)(CO)(6)(PEt(3))(2), B, and also on A(2)(-) and B(2+) clusters. The Fe-Fe distances in these molecules depend on the total number of valence electrons (60 e(-) in A and B(2)(+) and 62 e(-) in A(2)(-) and B) and undergo great structural changes upon addition or removal of electrons. The changes are consistent with known electron-counting rules in organometallic chemistry. The weak nature of the Fe-Fe bonding interactions in these clusters is apparent in the very similar energies of states with widely different Fe-Fe distances.

Density Functional Theory Calculations and Exploration of a Possible Mechanism of N2 Reduction by Nitrogenase

Density functional theory (DFT) calculations have been performed on the nitrogenase cofactor, FeMoco. Issues that have been addressed concern the nature of M-M interactions and the identity and origin of the central light atom, revealed in a recent crystallographic study of the FeMo protein of nitrogenase (Einsle, O.; et al. Science 2002, 297, 871). Introduction of Se in place of the S atoms in the cofactor and energy minimization results in an optimized structure very similar to that in the native enzyme. The nearly identical, short, lengths of the Fe-Fe distances in the Se and S analogues are interpreted in terms of M-M weak bonding interactions. DFT calculations with O or N as the central atoms in the FeMoco marginally support the assignment of the central atom as N rather than O. The assumption was made that the central atom is the N atom, and steps of a catalytic cycle were calculated starting with either of two possible states for the cofactor and maintaining the same charge throughout (by addition of equal numbers of H(+) and e(-)) between steps. The states were [(Cl)Fe(II)(6)Fe(III)Mo(IV)S(9)(H(+))(3)N(3-)(Gl)(Im)](2-), [I-N-3H](2-), and [(Cl)Fe(II)(4)Fe(III)(3)Mo(IV)S(9)(H(+))(3)N(3-)(Gl)(Im)], [I-N-3H](0) (Gl = deprotonated glycol; Im = imidazole). These are the triply protonated ENDOR/ESEEM [I-N](5-) and Mössbauer [I-N](3-) models, respectively. The proposed mechanism explores the possibilities that (a) redox-induced distortions facilitate insertion of N(2) and derivative substrates into the Fe(6) central unit of the cofactor, (b) the central atom in the cofactor is an exchangeable nitrogen, and (c) the individual steps are related by H(+)/e(-) additions (and reduction of substrate) or aquation/dehydration (and distortion of the Fe(6) center). The Delta E's associated with the individual steps of the proposed mechanism are small and either positive or negative. The largest positive Delta E is +121 kJ/mol. The largest negative Delta E of -333 kJ/mol is for the FeMoco with a N(3-) in the center (the isolated form) and an intermediate in the proposed mechanism.

Density Functional Calculations on the Binding of Dinitrogen to the FeFe Cofactor in Fe-Only Nitrogenase: FeFeco(μ6-N2) as Intermediate in Nitrogen Fixation

The geometries and stabilities of the FeFe cofactor at different oxidation states and its complexes with N(2) have been determined by density functional calculations. These calculations support an EPR-inactive resting state of the FeFe cofactor with four Fe(2+) and four Fe(3+) sites (4Fe(2+)4Fe(3+)). FeFeco(mu(6)-N(2)) with a central dinitrogen ligand is predicted to be the most stable complex of the FeFe cofactor with N(2). It is easily formed by penetration of N(2) into the trigonal Fe(6) prism of the FeFe cofactor with an approximate barrier of 4 kcal mol(-1). The present DFT results suggest that an FeFeco(mu(6)-N(2)) entity is a plausible intermediate in dinitrogen fixation by nitrogenase. CO is calculated to bind even more strongly than N(2) to the FeFe cofactor so that CO may inhibit the reduction of nitrogen by Fe-only nitrogenase.

Structural, spectroscopic, and redox consequences of a central ligand in the FeMoco of nitrogenase: a density functional theoretical study

Broken symmetry density functional and electrostatics calculations have been used to shed light on which of three proposed atoms, C, N, or O, is most likely to be present in the center of the FeMoco, the active site of nitrogenase. At the Mo(4+)4Fe(2+)3Fe(3+) oxidation level, a central N(3-) anion results in (1) calculated Fe-N bond distances that are in very good agreement with the recent high-resolution X-ray data of Einsle et al.; (2) a calculated redox potential of 0.19 eV versus the standard hydrogen electrode (SHE) for FeMoco(oxidized) + e(-) --> FeMoco(resting), in good agreement with the measured value of -0.042 V in Azotobacter vinelandii; and (3) average Mössbauer isomer shift values (IS(av) = 0.48 mm s(-1)) compatible with experiment (IS(av) = 0.40 mm s(-1)). At the more reduced Mo(4+)6Fe(2+)1Fe(3+) level, the calculated geometry around a central N(3-) anion still correlates well with the X-ray data, but the average Mössbauer isomer shift value (IS(av) = 0.54 mm s(-1)) and the redox potential of -2.21 eV show a much poorer agreement with experiment. These calculated structural, spectroscopic, and redox data indicate the most likely iron oxidation state for the resting FeMoco of nitrogenase to be 4Fe(2+)3Fe(3+). At this favored oxidation state, oxygen or carbon coordination leads to (1) Fe-O distances in poor agreement and Fe-C distances in good agreement with experiment and (2) calculated redox potentials of +0.97 eV for O(2-) and -1.31 eV for C(4-). The calculated structural parameters and/or redox data suggest either O(2-) or C(4-) is unlikely as a central anion.

Metal substitution in the active site of nitrogenase MFe(7)S(9) (M = Mo(4+), V(3+), Fe(3+))

The unifying view that molybdenum is the essential component in nitrogenase has changed over the past few years with the discovery of a vanadium-containing nitrogenase and an iron-only nitrogenase. The principal question that has arisen for the alternative nitrogenases concerns the structures of their corresponding cofactors and their metal-ion valence assignments and whether there are significant differences with that of the more widely known molybdenum-iron cofactor (FeMoco). Spin-polarized broken-symmetry (BS) density functional theory (DFT) calculations are used to assess which of the two possible metal-ion valence assignments (4Fe(2+)4Fe(3+) or 6Fe(2+)2Fe(3+)) for the iron-only cofactor (FeFeco) best represents the resting state. For the 6Fe(2+)2Fe(3+) oxidation state, the spin coupling pattern for several spin state alignments compatible with S = 0 were generated and assessed by energy criteria. The most likely BS spin state is composed of a 4Fe cluster with spin S(a) = (7)/(2) antiferromagnetically coupled to a 4Fe' cluster with spin S(b) = (7)/(2). This state has the lowest DFT energy for the isolated FeFeco cluster and displays calculated Mössbauer isomer shifts consistent with experiment. Although the S = 0 resting state of FeFeco has recently been proposed to have metal-ion valencies of 4Fe(2+)4Fe(3+) (derived from experimental Mössbauer isomer shifts), our isomer shift calculations for the 4Fe(2+)4Fe(3+) oxidation state are in poorer agreement with experiment. Using the Mo(4+)6Fe(2+)Fe(3+) oxidation level of the cofactor as a starting point, the structural consequences of replacement of molybdenum (Mo(4+)) with vanadium (V(3+)) or iron (Fe(3+)) in the cofactor have been investigated. The size of the cofactor cluster shows a dependency on the nature of the heterometal and increases in the order FeMoco < FeVco < FeFeco.

Substitution reactions of [Fe(4)S(4)Cl(4)](2-) with Bu(t)NC or Et(2)NCS(2)(-): trapping intermediates and detecting new pathways

Kinetic studies on the substitution reaction between [Fe(4)S(4)Cl(4)](2-) and Bu(t)NC or Et(2)NCS(2)(-) are reported. The binding of small molecules and ions to Fe-S clusters is a fundamental step in substitution reactions but can be difficult to follow directly because these reactions are rapid and often associated with small spectroscopic changes. A novel kinetic method is reported which allows the time course of molecule and ion binding to Fe-S clusters to be followed by monitoring the lability of the cluster. Using a stopped-flow, sequential-mix apparatus, [Fe(4)S(4)Cl(4)](2-) and L (L = Et(2)NCS(2)(-) or Bu(t)NC) are rapidly mixed, and after a known time (delta) the resulting solution is mixed with a solution of PhS(-). The thiolate substitutes for the chloro ligands on the cluster, in a reaction which is easy to follow because of the large change in the visible absorption spectrum. The rate of this substitution is extremely sensitive to whether L is bound to the cluster or not. By correlation of delta with the rate of the reaction with PhS(-), the time course of the reaction between [Fe(4)S(4)Cl(4)](2-) and L can be mapped out. In studies where L = Bu(t)NC this technique has allowed the detection of an intermediate ([Fe(4)S(4)Cl(4)(CNBu(t))](2-)) which cannot be detected spectrophotometrically. In further studies, the substitution reactions of [Fe(4)S(4)Cl(4)](2-) with PhS(-), Et(2)NCS(2)(-), or Bu(t)NC are all perturbed by the addition of Cl(-). In all cases a common pathway for substitution is evident, but with Et(2)NCS(2)(-) an additional, slower pathway becomes apparent under conditions where the common pathway is completely inhibited by Cl(-).

Nitrogenase MoFe-protein at 1.16 A resolution: a central ligand in the FeMo-cofactor

A high-resolution crystallographic analysis of the nitrogenase MoFe-protein reveals a previously unrecognized ligand coordinated to six iron atoms in the center of the catalytically essential FeMo-cofactor. The electron density for this ligand is masked in structures with resolutions lower than 1.55 angstroms, owing to Fourier series termination ripples from the surrounding iron and sulfur atoms in the cofactor. The central atom completes an approximate tetrahedral coordination for the six iron atoms, instead of the trigonal coordination proposed on the basis of lower resolution structures. The crystallographic refinement at 1.16 angstrom resolution is consistent with this newly detected component being a light element, most plausibly nitrogen. The presence of a nitrogen atom in the cofactor would have important implications for the mechanism of dinitrogen reduction by nitrogenase.

Binding modes for the first coupled electron and proton addition to FeMoco of nitrogenase

A combined broken-symmetry density functional and electrostatics approach has been used to model the one-electron reduced and protonated state of the iron-molybdenum cofactor active site of nitrogenase. The active site of the protein contains Fe, Mo, S, N, and O atoms, and many possible sites for protonation have been examined. A novel hydridic proton asymmetrically located in the central cavity created by six Fe sites is most favored from the calculations. Under physiological turnover conditions of low electron flux, the formation of this iron-hydride intermediate may represent a first step towards cofactor liberation of dihydrogen in the absence of dinitrogen.

Insights into properties and energetics of iron-sulfur proteins from simple clusters to nitrogenase

Some of the principal physical features of iron-sulfur clusters in proteins are analyzed, including metal-ligand covalency, spin polarization, spin coupling, valence delocalization, valence interchange and small reorganization energies, with emphasis on recent spectroscopic and theoretical work. The current state of structural, spectroscopic, and computational knowledge for the iron-sulfur clusters in the nitrogenase iron and iron-molybdenum proteins is examined by comparison and contrast to 'simpler' ironclusters. The differing interactions of the nitrogenase iron and iron-molybdenum clusters compared with those of other iron-sulfur clusters with the protein and solvent environment are also explored.

The Fe-only nitrogenase and the Mo nitrogenase from Rhodobacter capsulatus: a comparative study on the redox properties of the metal clusters present in the dinitrogenase components

The dinitrogenase component proteins of the conventional Mo nitrogenase (MoFe protein) and of the alternative Fe-only nitrogenase (FeFe protein) were both isolated and purified from Rhodobacter capsulatus, redox-titrated according to the same procedures and subjected to an EPR spectroscopic comparison. In the course of an oxidative titration of the MoFe protein (Rc1Mo) three significant S = 1/2 EPR signals deriving from oxidized states of the P-cluster were detected: (1) a rhombic signal (g = 2.07, 1.96 and 1.83), which showed a bell-shaped redox curve with midpoint potentials (Em) of -195 mV (appearance) and -30 mV (disappearance), (2) an axial signal (g(parallel) = 2.00, g perpendicular = 1.90) with almost identical redox properties and (3) a second rhombic signal (g = 2.03, 2.00, 1.90) at higher redox potentials (> 100 mV). While the 'low-potential' rhombic signal and the axial signal have been both attributed to the one-electron-oxidized P-cluster (P1+) present in two conformationally different proteins, the 'high-potential' rhombic signal has been suggested rather to derive from the P3+ state. Upon oxidation, the FeFe protein (Rc1Fe) exhibited three significant S = 1/2 EPR signals as well. However, the Rc1Fe signals strongly deviated from the MoFe protein signals, suggesting that they cannot simply be assigned to different P-cluster states. (a) The most prominent feature is an unusually broad signal at g = 2.27 and 2.06, which proved to be fully reversible and to correlate with catalytic activity. The cluster giving rise to this signal appears to be involved in the transfer of two electrons. The midpoint potentials determined were: -80 mV (appearance) and 70 mV (disappearance). (b) Under weakly acidic conditions (pH 6.4) a slightly altered EPR signal occurred. It was characterized by a shift of the g values to 2.22 and 2.05 and by the appearance of an additional negative absorption-shaped peak at g = 1.86. (c) A very narrow rhombic EPR signal at g = 2.00, 1.98 and 1.96 appeared at positive redox potentials (Em = 80 mV, intensity maximum at 160 mV). Another novel S = 1/2 signal at g = 1.96, 1.92 and 1.77 was observed on further, enzymatic reduction of the dithionite-reduced state of Rc1Fe with the dinitrogenase reductase component (Rc2Fe) of the same enzyme system (turnover conditions in the presence of N2 and ATP). When the Rc1Mo protein was treated analogously, neither this 'turnover signal' nor any other S = 1/2 signal were detectable. All Rc1Fe-specific EPR signals detected are discussed and tentatively assigned with special consideration of the reference spectra obtained from Rc1Mo preparations.

FeMo cofactor of nitrogenase: a density functional study of states M(N), M(OX), M(R), and M(I)

The M(N) S = (3)/(2) resting state of the FeMo cofactor of nitrogenase has been proposed to have metal-ion valencies of either Mo(4+)6Fe(2+)Fe(3+) (derived from metal hyperfine interactions) or Mo(4+)4Fe(2+)3Fe(3+) (from Mössbauer isomer shifts). Spin-polarized broken-symmetry (BS) density functional theory (DFT) calculations have been undertaken to determine which oxidation level best represents the M(N) state and to provide a framework for understanding its energetics and spectroscopy. For the Mo(4+)6Fe(2+)Fe(3+) oxidation state, the spin coupling pattern for several spin state alignments compatible with S = (3)/(2) were generated and assessed by energy and geometric criteria. The most likely BS spin state is composed of a Mo3Fe cluster with spin S(a) = 2 antiferromagnetically coupled to a 4Fe' cluster with spin S(b) = (7)/(2). This state has a low DFT energy for the isolated FeMoco cluster and the lowest energy when the interaction with the protein and solvent environment is included. This spin state also displays calculated metal hyperfine and Mössbauer isomer shifts compatible with experiment, and optimized geometries that are in excellent agreement with the protein X-ray data. Our best model for the actual spin-coupled state within FeMoco alters this BS state by a slight canting of spins and is analogous in several respects to that found in the 8Fe P-cluster in the same protein. The spin-up and spin-down components of the LUMO contain atomic contributions from Mo(4+) and the homocitrate and from the central prismane Fe sites and muS(2) atoms, respectively. This qualitative picture of the accepting orbitals for M(N) is consistent with observations from Mössbauer spectra of the one-electron reduced states. Similar calculations for the Mo(4+)4Fe(2+)3Fe(3+) oxidation state yield results that are in poorer agreement with experiment. Using the Mo(4+)6Fe(2+)Fe(3+) oxidation level as the most plausible resting state, the geometric, electronic and energetic properties of the one-electron redox transition to the oxidized state, M(OX), catalytically observed M(R) and radiolytically reduced M(I) states have also been explored.

Theoretical studies of biological nitrogen fixation. I. Density functional modeling of the Mo-site of the FeMo-cofactor

The Mo-site and its ligand environment of the FeMo-cofactor (FeMo-co) were studied using the hybrid density functional method B3LYP. The structure and stability of the model complex (S-ligand)3(N-ligand)Mo[(S)-OCH(CH3)C(O)O-] along with its various protonated and reduced/oxidized forms were calculated. Several hypotheses were tested: (i) ligand environment of the Mo-site, (ii) monodentate vs bidentate coordination of the Mo-bound homocitrate ligand, (iii) substrate coordination to the Mo center, and (iv) Mo-His interaction. It was found that the decoordination of one of the homocitrate (lactate in the model) "legs", the bidentate-->monodentate rearrangement, does not occur spontaneously upon either single/double protonation or one-electron reduction. However, it could occur only upon substrate coordination to the Mo-center of the single-protonated forms of the complex. It was shown that one-electron reduction, single-protonation, and substrate coordination facilitate the bidentate<-->monodentate rearrangement of the homocitrate (lactate) ligand of FeMo-co. It was demonstrated that the smallest acceptable model of His ligand in FeMo-co is methylimidazolate (MeIm-). Our studies suggest that the epsilon-N of the FeMo-co-bound His residue is not protonated, and as a consequence the cluster is tightly bound to the protein matrix via a strong Mo-N delta bond.