Conformations generated during turnover of the Azotobacter vinelandii nitrogenase MoFe protein and their relationship to physiological function

Understanding non-reducible N2 in the mechanism of Mo-nitrogenase

In my proposed mechanism of Mo-nitrogenase there are two roles for separate N2 molecules. One N2 diffuses into the reaction zone between Fe2 and Fe6 where a strategic gallery of H atoms can capture N2 to form the Fe-bound HNNH intermediate which is then progressively hydrogenated through intermediates containing HNNH2, NH and NH2 entities and then two NH3 in sequence. The second N2 can be parked in an N2-pocket about 3.2 Å from Fe2 or bind end-on at the exo coordination site of Fe2. This second N2 is outside the reaction zone, not exposed to H atom donors, and so is 'non-reducible'. Here density functional calculations using a 485+ atom model describe the thermodynamics for non-reducible N2 moving between the N2-pocket and the exo-Fe2 position, for the resting state and 19 intermediates in the mechanism. The entropy component is estimated and included. The result is that for all intermediates with ligation by H or NHx at the endo-Fe2 position the free energy for association of non-reducible N2 at exo-Fe2 is negative. There remains some uncertainty about the status of exo-Fe2-N2 during the step in which H2 exchanges with the incoming reducible N2, where at least two unbound molecules are present. At Fe2 it is evident that attainment of octahedral coordination stereochemistry dominates the binding thermodynamics for non-reducible N2. Possibilities for experimental support of these computational conclusions are discussed.

The mechanism of Mo-nitrogenase: from N2 capture to first release of NH3

Mo-nitrogenase hydrogenates N2 to NH3. This report continues from the previous paper [I. Dance, Dalton Trans., 2024, 53, 14193-14211] that described how the active site FeMo-co of the enzyme is uniquely able to capture and activate N2, forming a key intermediate with Fe-bound HNNH. Density functional simulations with a 485+ atom model of the active site and its surroundings are used to describe here the further reactions of this HNNH intermediate. The first step is hydrogenation to form HNNH2 bridging Fe2 and Fe6. Then a single-step reaction breaks the N-N bond, generating an Fe2-NH-Fe6 bridge and forming NH3 bound to Fe6. Then NH3 dissociates from Fe6. Reaction potential energies and kinetic barriers for all steps are reported for the most favourable electronic states of the system. The steps that follow the Fe2-NH-Fe6 intermediate, forming and dissociating the second NH3, and regenerating the resting state of the enzyme, are outlined. These results provide an interpretation of the recent steady-state kinetics data and analysis by Harris et al., [Biochemistry, 2022, 61, 2131-2137] who found a slow step after the formation of the HNNH intermediate. The calculated potential energy barriers for the HNNH2 → NH + NH3 reaction (30-36 kcal mol-1) are larger than the potential energy barriers for the N2 → HNNH reaction (19-29 kcal mol-1). I propose that the post-HNNH slow step identified kinetically is the key HNNH2 → NH + NH3 reaction described here. This step and the N2-capture step are the most difficult in the conversion of N2 to 2NH3. The steps in the complete mechanism still to be computationally detailed are relatively straightforward.

The activating capture of N2 at the active site of Mo-nitrogenase

Dinitrogen is inherently inert. This report describes detailed density functional calculations (with a 485+ atom model) of mechanistic steps by which the enzyme nitrogenase activates unreactive N2 at the intact active site FeMo-co, to form a key intermediate with bound HNNH. This mechanism does not bind N2 first and then add H atoms, but rather captures N2 ('N2-ready') that diffuses in through the substrate channel and enters a strategic gallery of H atom donors in the reaction zone, between Fe2 and Fe6. This occurs at the E4 stage of the complete mechanism. Exploration of possible reactions of N2 in this space leads to the conclusion that the first reaction step is transfer of H on Fe7 to one end of N2-ready, soon followed by Fe-N bond formation, and then a second H transfer from bridging S2BH to the other N. Two H-N bonds and one or two N-Fe bonds are formed, in some cases with a single transition state. The variable positions and orientations of N2-ready lead to various reaction trajectories and products. The favourable products resulting from this capture, judged by the criteria of reaction energies, reaction barriers, and mechanistic competence for further hydrogenation reactions in the nitrogenase cycle, have Fe2-NH-NH bonding. The trajectory of one N2 capture reaction is described in detail, and calculations that separate the H atom component and the 'heavy atom' components of the classical activation energy are described, in the context of possible H atom tunneling in the activation of N2-ready. I present arguments for the activation of N2 by the pathway of concerted hydrogenation and binding of N2-ready, alternative to the commonly assumed pathway of binding N2 first, with subsequent hydrogenation. The active site of nitrogenase is well primed for the thermodynamic and kinetic advantages of N2 capture.

Analysis of early intermediate states of the nitrogenase reaction by regularization of EPR spectra

Due to the complexity of the catalytic FeMo cofactor site in nitrogenases that mediates the reduction of molecular nitrogen to ammonium, mechanistic details of this reaction remain under debate. In this study, selenium- and sulfur-incorporated FeMo cofactors of the catalytic MoFe protein component from Azotobacter vinelandii are prepared under turnover conditions and investigated by using different EPR methods. Complex signal patterns are observed in the continuous wave EPR spectra of selenium-incorporated samples, which are analyzed by Tikhonov regularization, a method that has not yet been applied to high spin systems of transition metal cofactors, and by an already established grid-of-error approach. Both methods yield similar probability distributions that reveal the presence of at least four other species with different electronic structures in addition to the ground state E0. Two of these species were preliminary assigned to hydrogenated E2 states. In addition, advanced pulsed-EPR experiments are utilized to verify the incorporation of sulfur and selenium into the FeMo cofactor, and to assign hyperfine couplings of 33S and 77Se that directly couple to the FeMo cluster. With this analysis, we report selenium incorporation under turnover conditions as a straightforward approach to stabilize and analyze early intermediate states of the FeMo cofactor.

The HD Reaction of Nitrogenase: a Detailed Mechanism

Nitrogenase is the enzyme that converts N2 to NH3 under ambient conditions. The chemical mechanism of this catalysis at the active site FeMo-co [Fe7 S9 CMo(homocitrate)] is unknown. An obligatory co-product is H2 , while exogenous H2 is a competitive inhibitor. Isotopic substitution using exogenous D2 revealed the N2 -dependent reaction D2 +2H+ +2e- →2HD (the 'HD reaction'), together with a collection of additional experimental characteristics and requirements. This paper describes a detailed mechanism for the HD reaction, developed and elaborated using density functional simulations with a 486-atom model of the active site and surrounding protein. First D2 binds at one Fe atom (endo-Fe6 coordination position), where it is flanked by H-Fe6 (exo position) and H-Fe2 (endo position). Then there is synchronous transfer of these two H atoms to bound D2 , forming one HD bound to Fe2 and a second HD bound to Fe6. These two HD dissociate sequentially. The final phase is recovery of the two flanking H atoms. These H atoms are generated, sequentially, by translocation of a proton from the protein surface to S3B of FeMo-co and combination with introduced electrons. The first H atom migrates from S3B to exo-Fe6 and the second from S3B to endo-Fe2. Reaction energies and kinetic barriers are reported for all steps. This mechanism accounts for the experimental data: (a) stoichiometry; (b) the N2 -dependence results from promotional N2 bound at exo-Fe2; (c) different N2 binding Km for the HD reaction and the NH3 formation reaction results from involvement of two different sites; (d) inhibition by CO; (e) the non-occurrence of 2HD→H2 +D2 results from the synchronicity of the two transfers of H to D2 ; (f) inhibition of HD production at high pN2 is by competitive binding of N2 at endo-Fe6; (g) the non-leakage of D to solvent follows from the hydrophobic environment and irreversibility of proton introduction.

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.

Defining Intermediates of Nitrogenase MoFe Protein during N2 Reduction under Photochemical Electron Delivery from CdS Quantum Dots

Coupling the nitrogenase MoFe protein to light-harvesting semiconductor nanomaterials replaces the natural electron transfer complex of Fe protein and ATP and provides low-potential photoexcited electrons for photocatalytic N₂ reduction. A central question is how direct photochemical electron delivery from nanocrystals to MoFe protein is able to support the multielectron ammonia production reaction. In this study, low photon flux conditions were used to identify the initial reaction intermediates of CdS quantum dot (QD):MoFe protein nitrogenase complexes under photochemical activation using EPR. Illumination of CdS QD:MoFe protein complexes led to redox changes in the MoFe protein active site FeMo-co observed as the gradual decline in the E₀ resting state intensity that was accompanied by an increase in the intensity of a new "g_eff = 4.5" EPR signal. The magnetic properties of the g_eff = 4.5 signal support assignment as a reduced S = 3/2 state, and reaction modeling was used to define it as a two-electron-reduced "E₂" intermediate. Use of a MoFe protein variant, β-188^Cys, which poises the P cluster in the oxidized P⁺ state, demonstrated that the P cluster can function as a site of photoexcited electron delivery from CdS to MoFe protein. Overall, the results establish the initial steps for how photoexcited CdS delivers electrons into the MoFe protein during reduction of N₂ to ammonia and the role of electron flux in the photochemical reaction cycle.

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.

Hydride Conformers of the Nitrogenase FeMo-cofactor Two-Electron Reduced State E2(2H), Assigned Using Cryogenic Intra Electron Paramagnetic Resonance Cavity Photolysis

Early studies in which nitrogenase was freeze-trapped during enzymatic turnover revealed the presence of high-spin ( S = ³/₂) electron paramagnetic resonance (EPR) signals from the active-site FeMo-cofactor (FeMo-co) in electron-reduced intermediates of the MoFe protein. Historically denoted as 1b and 1c, each of the signals is describable as a fictitious spin system, S' = ¹/₂, with anisotropic g' tensor, 1b with g' = [4.21, 3.76, ?] and 1c with g' = [4.69, ∼3.20, ?]. A clear discrepancy between the magnetic properties of 1b and 1c and the kinetic analysis of their appearance during pre-steady-state turnover left their identities in doubt, however. We subsequently associated 1b with the state having accumulated 2[e^-/H⁺], denoted as E₂(2H), and suggested that the reducing equivalents are stored on the catalytic FeMo-co cluster as an iron hydride, likely an [Fe-H-Fe] hydride bridge. Intra-EPR cavity photolysis (450 nm; temperature-independent from 4 to 12 K) of the E₂(2H)/1b state now corroborates the identification of this state as storing two reducing equivalents as a hydride. Photolysis converts E₂(2H)/1b to a state with the same EPR spectrum, and thus the same cofactor structure as pre-steady-state turnover 1c, but with a different active-site environment. Upon annealing of the photogenerated state at temperature T = 145 K, it relaxes back to E₂(2H)/1b. This implies that the 1c signal comes from an E₂(2H) hydride isomer of E₂(2H)/1b that stores its two reducing equivalents either as a hydride bridge between a different pair of iron atoms or an Fe-H terminal hydride.

Negative cooperativity in the nitrogenase Fe protein electron delivery cycle

Nitrogenase catalyzes the ATP-dependent reduction of dinitrogen (N₂) to two ammonia (NH₃) molecules through the participation of its two protein components, the MoFe and Fe proteins. Electron transfer (ET) from the Fe protein to the catalytic MoFe protein involves a series of synchronized events requiring the transient association of one Fe protein with each αβ half of the α₂β₂ MoFe protein. This process is referred to as the Fe protein cycle and includes binding of two ATP to an Fe protein, association of an Fe protein with the MoFe protein, ET from the Fe protein to the MoFe protein, hydrolysis of the two ATP to two ADP and two P_i for each ET, P_i release, and dissociation of oxidized Fe protein-(ADP)₂ from the MoFe protein. Because the MoFe protein tetramer has two separate αβ active units, it participates in two distinct Fe protein cycles. Quantitative kinetic measurements of ET, ATP hydrolysis, and P_i release during the presteady-state phase of electron delivery demonstrate that the two halves of the ternary complex between the MoFe protein and two reduced Fe protein-(ATP)₂ do not undergo the Fe protein cycle independently. Instead, the data are globally fit with a two-branch negative-cooperativity kinetic model in which ET in one-half of the complex partially suppresses this process in the other. A possible mechanism for communication between the two halves of the nitrogenase complex is suggested by normal-mode calculations showing correlated and anticorrelated motions between the two halves.

Reversible Protonated Resting State of the Nitrogenase Active Site

Protonated states of the nitrogenase active site are mechanistically significant since substrate reduction is invariably accompanied by proton uptake. We report the low pH characterization by X-ray crystallography and EPR spectroscopy of the nitrogenase molybdenum iron (MoFe) proteins from two phylogenetically distinct nitrogenases (Azotobacter vinelandii, Av, and Clostridium pasteurianum, Cp) at pHs between 4.5 and 8. X-ray data at pHs of 4.5–6 reveal the repositioning of side chains along one side of the FeMo-cofactor, and the corresponding EPR data shows a new S = 3/2 spin system with spectral features similar to a state previously observed during catalytic turnover. The structural changes suggest that FeMo-cofactor belt sulfurs S3A or S5A are potential protonation sites. Notably, the observed structural and electronic low pH changes are correlated and reversible. The detailed structural rearrangements differ between the two MoFe proteins, which may reflect differences in potential protonation sites at the active site among nitrogenase species. These observations emphasize the benefits of investigating multiple nitrogenase species. Our experimental data suggest that reversible protonation of the resting state is likely occurring, and we term this state “E₀H⁺”, following the Lowe–Thorneley naming scheme.

Structural characterization of CO-inhibited Mo-nitrogenase by combined application of nuclear resonance vibrational spectroscopy, extended X-ray absorption fine structure, and density functional theory: new insights into the effects of CO binding and the role of the interstitial atom

The properties of CO-inhibited Azotobacter vinelandii (Av) Mo-nitrogenase (N2ase) have been examined by the combined application of nuclear resonance vibrational spectroscopy (NRVS), extended X-ray absorption fine structure (EXAFS), and density functional theory (DFT). Dramatic changes in the NRVS are seen under high-CO conditions, especially in a 188 cm(-1) mode associated with symmetric breathing of the central cage of the FeMo-cofactor. Similar changes are reproduced with the α-H195Q N2ase variant. In the frequency region above 450 cm(-1), additional features are seen that are assigned to Fe-CO bending and stretching modes (confirmed by (13)CO isotope shifts). The EXAFS for wild-type N2ase shows evidence for a significant cluster distortion under high-CO conditions, most dramatically in the splitting of the interaction between Mo and the shell of Fe atoms originally at 5.08 Å in the resting enzyme. A DFT model with both a terminal -CO and a partially reduced -CHO ligand bound to adjacent Fe sites is consistent with both earlier FT-IR experiments, and the present EXAFS and NRVS observations for the wild-type enzyme. Another DFT model with two terminal CO ligands on the adjacent Fe atoms yields Fe-CO bands consistent with the α-H195Q variant NRVS. The calculations also shed light on the vibrational "shake" modes of the interstitial atom inside the central cage, and their interaction with the Fe-CO modes. Implications for the CO and N2 reactivity of N2ase are discussed.

A molecular pathway for the egress of ammonia produced by nitrogenase

Nitrogenase converts N2 to NH3, at one face of an Fe-Mo-S cluster (FeMo-co) buried in the protein. Through exploration of cavities in the structures of nitrogenase proteins, a pathway for the egress of ammonia from its generation site to the external medium is proposed. This pathway is conserved in the three species Azotobacter vinelandii, Klebsiella pneumoniae and Clostridium pasteurianum. A molecular mechanism for the translocation of NH3 by skipping through a sequence of hydrogen bonds involving eleven water molecules and surrounding aminoacids has been developed. The putative mechanism requires movement aside of some water molecules by up to ~ 1Å. Consistent with this, the surrounding protein is comprised of different chains and has little secondary structure: protein fluctuations are part of the mechanism. This NH3 pathway is well separated from the water chain and embedded proton wire that have been proposed for serial supply of protons to FeMo-co. Verification procedures are suggested.

Nitrogen fixation and hydrogen metabolism in cyanobacteria

This review summarizes recent aspects of (di)nitrogen fixation and (di)hydrogen metabolism, with emphasis on cyanobacteria. These organisms possess several types of the enzyme complexes catalyzing N(2) fixation and/or H(2) formation or oxidation, namely, two Mo nitrogenases, a V nitrogenase, and two hydrogenases. The two cyanobacterial Ni hydrogenases are differentiated as either uptake or bidirectional hydrogenases. The different forms of both the nitrogenases and hydrogenases are encoded by different sets of genes, and their organization on the chromosome can vary from one cyanobacterium to another. Factors regulating the expression of these genes are emerging from recent studies. New ideas on the potential physiological and ecological roles of nitrogenases and hydrogenases are presented. There is a renewed interest in exploiting cyanobacteria in solar energy conversion programs to generate H(2) as a source of combustible energy. To enhance the rates of H(2) production, the emphasis perhaps needs not to be on more efficient hydrogenases and nitrogenases or on the transfer of foreign enzymes into cyanobacteria. A likely better strategy is to exploit the use of radiant solar energy by the photosynthetic electron transport system to enhance the rates of H(2) formation and so improve the chances of utilizing cyanobacteria as a source for the generation of clean energy.

Conformational gating of electron transfer from the nitrogenase Fe protein to MoFe protein

The nitrogenase Fe protein contains a [4Fe-4S] cluster and delivers one electron at a time to the catalytic MoFe protein. During this electron delivery, the Fe protein in its [4Fe-4S](1+) reduced state (Fe(red)) binds two MgATP and forms a complex with the MoFe protein, with subsequent transfer of one electron to the MoFe protein in a reaction coupled to the hydrolysis of two ATP. Crystal structures with the nitrogenase complex in different nucleotide-bound states show major conformational changes which provide a structural underpinning to suggestions that intercomponent electron transfer (ET) is "gated" by conformational changes of the complex and/or of its component proteins. Although electron delivery is coupled to ATP hydrolysis, their connection is puzzling, for it appears that ET precedes both ATP hydrolysis and Pi release. We here test the gating hypothesis with studies of the intracomplex oxidation of Fe(red) by MoFe protein in the presence of a variety of solutes. Conformational control of this process (gating) is revealed by the finding that it responds to changes in osmotic pressure (but not viscosity), with no fewer than 80 waters being bound during the reaction. The absence of a solvent kinetic isotope effect further implies that ATP hydrolysis does not occur during the rate-limiting step of ET.