Nitrogenase is the only enzyme that can cleave the triple bond in N 2 , making nitrogen available for all lifeforms. Previous computational studies have given widely diverging results regarding the reaction mechanism of the enzyme. For example, some recent studies have suggested that one of the μ 2 -bridging sulfide ligands (S2B) may dissociate from one of the Fe ions when protonated in the doubly reduced and protonated E 2 state, whereas other studies indicated that such half-dissociated states are unfavorable. We have examined how the relative energies of 26 structures of the E 2 state depend on details of combined quantum mechanical and molecular mechanical (QM/MM) calculations. We show that the selection of the broken-symmetry state, the basis set, relativistic effects, the size of the QM system, relaxation of the surroundings, and the conformations of the bound protons may affect the relative energies of the various structures by up to 12, 22, 9, 20, 37, and 33 kJ/mol, respectively. However, they do not change the preferred type of structures. On the other hand, the choice of the DFT functional strongly affects the preferences. The hybrid B3LYP functional strongly prefers doubly protonation of the central carbide ion, but such a structure is not consistent with experimental EPR data. Other functionals suggest structures with a hydride ion, in agreement with the experiments, and show that the ion bridges between Fe2 and Fe6. Moreover, there are two structures of the same type that are degenerate within 1−5 kJ/mol, in agreement with the observation of two EPR signals. However, the pure generalized gradient approximation (GGA) functional TPSS favors structures with a protonated S2B also bridging Fe2 and Fe6, whereas r 2 SCAN (meta-GGA) and TPSSh (hybrid) prefer structures with S2B dissociated from Fe2 (but remaining bound to Fe6). The energy difference between the two types of structure is so small (7−18 kJ/mol) that both types need to be considered in future investigations of the mechanism of nitrogenase.

Nitrogenase is the only enzyme that can convert N2 to NH3. Crystallographic structures have indicated that one of the sulfide ligands of the active‐site FeMo cluster, S2B, can be replaced by an inhibitor, like CO and OH−, and it has been suggested that it may be displaced also during the normal reaction. We have investigated possible proton transfer pathways within the FeMo cluster during the conversion of N2H2 to two molecules of NH3, assuming that the protons enter the cluster at the S3B, S4B or S5A sulfide ions and are then transferred to the substrate. We use combined quantum mechanical and molecular mechanical (QM/MM) calculations with the TPSS and B3LYP functionals. The calculations indicate that the barriers for these reactions are reasonable if the S2B ligand remains bound to the cluster, but they become prohibitively high if S2B has dissociated. This suggests that it is unlikely that S2B reversibly dissociates during the normal reaction cycle.

Nitrogenase is the only enzyme that can convert N 2 to NH 3 . Crystallographic structures have indicated that one of the sulfide ligands of the active-site FeMo cluster, S2B, can be replaced by an inhibitor, like CO and OH À , and it has been suggested that it may be displaced also during the normal reaction. We have investigated possible proton transfer pathways within the FeMo cluster during the conversion of N 2 H 2 to two molecules of NH 3 , assuming that the protons enter the cluster at the S3B, S4B or S5A sulfide ions and are then transferred to the substrate. We use combined quantum mechanical and molecular mechanical (QM/MM) calculations with the TPSS and B3LYP functionals. The calculations indicate that the barriers for these reactions are reasonable if the S2B ligand remains bound to the cluster, but they become prohibitively high if S2B has dissociated. This suggests that it is unlikely that S2B reversibly dissociates during the normal reaction cycle.

We have calculated redox potentials of the two metal clusters in Mo-nitrogenase with quantum mechanical (QM) calculations. We employ an approach calibrated for iron–sulfur clusters with 1–4 Fe ions, involving QM-cluster calculations in continuum solvent and large QM systems (400–500 atoms), based on structures from combined QM and molecular mechanics (QM/MM) geometry optimisations. Calculations on the P-cluster show that we can reproduce the experimental redox potentials within 0.33 V. This is similar to the accuracy obtained for the smaller clusters, although two of the redox reactions involve also proton transfer. The calculated P1+/PN redox potential is nearly the same independently of whether P1+ is protonated or deprotonated, explaining why redox titrations do not show any pH dependence. For the FeMo cluster, the calculations clearly show that the formal oxidation state of the cluster in the resting E0 state is , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E0–E4 states are nearly constant, as is expected if the electrons are delivered by the same site (the P-cluster). However, the redox potentials are insensitive to the formal oxidation states of the Fe ion (i.e., whether the added protons bind to sulfide or Fe ions). Finally, we show that the later (E4–E8) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E0/E1 couple.