A Systematic DFT Approach for Studying Mechanisms of Redox Active Enzymes

Unusual mechanism of aziridine biosynthesis catalysed by the αKG-dependent non-heme enzyme TqaL

Aziridines are present in many synthetic pharmaceuticals. The synthesis of the aziridine ring remains challenging due to its highly strained three-membered ring structure. Recently, a non-heme αKG-dependent enzyme, TqaL, has been demonstrated to catalyze the synthesis of aziridines from L-Val. However, the detailed reaction mechanism of the enzyme remains elusive. Herein, we reported, for the first time, the mechanism of oxidative cyclisation for aziridine synthesis catalyzed by TqaL. Following the HAA step, the reaction proceeded via a unique concerted process with a single electron transfer from the isopropyl radical to the Fe(III)-OH motif, which was coupled with the electrophilic attack of the primary amine substrate on the tertiary isopropyl radical and simultaneous proton transfer from the substrate amine to the hydroxyl group of the Fe(III)-OH to give the aziridine. This research would provide a valuable structural basis for tailoring the non-heme αKG-dependent enzyme for the biosynthesis of highly active aziridine derivatives as pharmaceuticals.

Reduction of O2 and NO in flavodiiron proteins - Tuning the energy landscape by second sphere ligation variations

Flavodiiron proteins (FDPs) constitute a large family of non-heme iron enzymes present in all domains of life. They play important roles as scavengers and detoxifiers by efficiently reducing both O2 and NO. The primary ligands of the diiron active site in all FDPs are highly conserved, indicating that the basic reaction mechanisms for O2 and NO reduction, respectively, are the same. However, the reduction activity varies significantly between different FDPs. By comparing FDPs from two different species, Thermotoga maritima and Desulfovibrio gigas, we investigate to what extent variations in the second sphere ligation can explain differences in reduction activities. Comparisons are also made between wildtype and two variants of Thermotoga maritima FDP. We use Density functional theory (DFT) calculations on a number of FDP active site models to study the reaction mechanisms for both O2 and NO reduction. For reduction of O2 we conclude that differences in activity cannot be explained by differences in the first or second active site coordination spheres, which is mainly due to a low barrier for OO bond cleavage after one proton-coupled reduction step. For NO reduction however, the rate-limiting barrier for N2O formation, a hyponitrite rotation, is high enough to be involved in the overall rate limitation. We show that second sphere residues, such as Tyr26 in Desulfovibrio gigas FDP, that can form hydrogen bonds to the rotating hyponitrite, decrease the barrier. Differences in NO reduction rate among different FDPs are most likely determined by the variation in such second sphere residues.

The Mechanism of Nitrite Reductase

Cytochrome c nitrite reductase (CcNiR) activates nitrite and produces ammonia. It is one of several enzymes that use a redox-active cofactor to perform its reaction. In this case, the cofactor has a heme with a lysine as the proximal ligand and a charged nearby arginine. The role of a tyrosine, which is also close, has been less clear. There are also four bis-histidine-ligated hemes involved in the electron transfers. CcNiR has been studied before, using essentially the same methods as here. However, the mechanism is very complicated, involving six reductions, and quite different results for the mechanism have been obtained here. For example, the tyrosine has here been found to be redox active in the final step when ammonia is produced. Also, the arginine has here been found to stay protonated throughout the mechanism, which is different from what was found in the previous study. The present results are in very good agreement with experimental findings and are, therefore, another case where the methodology has been shown to work very well. Previous examples include Photosystem II and Nitrogenase, normally considered to be the most important enzymes in nature for the development of life.

Computational Modeling of the Enzymatic Achmatowicz Rearrangement by Heme-Dependent Chloroperoxidase: Reaction Mechanism, Enantiopreference, Regioselectivity, and Substrate Specificity

The chloroperoxidase from Caldariomyces fumago (CfCPO) catalyzes the oxidative ring expansion of α-heterofunctionalized furans via the Achmatowicz rearrangement, providing an elegant tool to convert furan rings into complex-prefunctionalized scaffolds. However, the mechanism of this transformation remains unclear. Herein, the CfCPO-catalyzed reaction of rac-1-(2-furyl)ethanol (1a) is studied by quantum chemical calculations and molecular dynamics simulations. The calculations reveal that the conversion follows the general mechanism of the Achmatowicz reaction. Notably, the binding of 1a to the enzyme's active site influences the Compound I (Cpd I) formation, and the (R)-1a enantiomer binding results in a lower barrier compared to (S)-1a, explaining the observed (R)-enantiopreference toward a racemic substrate. Additionally, due to the weaker steric hindrance between the porphyrin ring and substrate, the nucleophilic attack of Cpd I on the furan core of 1a is preferred at the less-substituted C4=C5 bond, providing a rationale for the experimentally observed regioselectivity. Finally, the bottleneck residues in the substrate delivery channel and also the active site surroundings are proposed to be responsible for the substrate specificity of CfCPO. This study lays a theoretical foundation for the rational design of new CPOs that catalyze the Achmatowicz rearrangement with a broader substrate spectrum or specific stereopreference.

Nitrification Mechanisms for the P460 Enzymes

The oxidation of hydroxylamine was studied by quantum chemical modeling. Hydroxylamine is the product of ammonia oxidation in ammonia monooxygenase. That mechanism has been studied recently by quantum chemical modeling as here. Only two enzymes can oxidize hydroxylamine, hydroxylamine oxidase and cytochrome-P460. Both employ the unusual P460-heme cofactor. In hydroxylamine oxidase, there is a covalently linked tyrosine, while in cytochrome-P460, there is a covalently linked lysine. The calculations give explanations for the experimental findings that NO is the final product in hydroxylamine oxidase, while N2O is the final product in cytochrome-P460. The effect of the covalent attachments has been investigated, and reasons for their presence have been given. The methodology used, which was proven to give very good agreement with experiments for several redox enzymes, again leads to excellent agreement with experimental findings.

Sulfide release and rebinding in the mechanism for nitrogenase

Nitrogenases are the only enzymes that activate the strong triple bond in N2. The mechanism for the activation has been very difficult to determine in spite of decades of work. In previous modeling studies it has been suggested that the mechanism for nitrogen activation starts out by four pre-activation steps (A0-A4) before catalysis. That suggestion led to excellent agreement with experimental Elecrtron Paramagnetic Resonance (EPR) observations in the step where N2 becomes protonated (E4). An important part of the pre-activation is that a sulfide is released. In the present paper, the details of the pre-activation are modeled, including the release of the sulfide. Several possible transition states for the release have been obtained. An A4(E0) state is reached which is very similar to the E4 state. For completeness, the steps going back from A4(E0) to A0 after catalysis are also modeled, including the insertion of a sulfide.

Gaining Insight into the Catalytic Mechanism of the R132H IDH1 Mutant: A Synergistic DFT Cluster and Experimental Investigation

Human isocitrate dehydrogenase 1 (IDH1) is an enzyme that is found in humans that plays a critical role in aerobic metabolism. As a part of the citric acid cycle, IDH1 becomes responsible for catalyzing the oxidative decarboxylation of isocitrate to form α-ketoglutarate (αKG), with nicotinamide adenine dinucleotide phosphate (NADP+) as a cofactor. Strikingly, mutations of the IDH1 enzyme have been discovered in several cancers including glioblastoma multiforme (GBM), a highly aggressive form of brain cancer. It has been experimentally determined that single-residue IDH1 mutations occur at a very high frequency in GBM. Specifically, the IDH1 R132H mutation is known to produce (D)2-hydroxyglutarate (2HG), a recognized oncometabolite. Using the previously determined catalytic mechanism of IDH1, a DFT QM model was developed to study the mechanistic properties of IDH1 R132H compared to wild type enzyme. Validating these insights, biochemical in vitro assays of metabolites produced by mutant vs wild type enzymes were measured and compared. From the results discussed herein, we discuss the mechanistic impact of mutations in IDH1 on its ability to catalyze the formation of αKG and 2HG.

Final E5 to E8 Steps in the Nitrogenase Mechanism for Nitrogen Fixation

Nitrogenase converts nitrogen in the air to ammonia. It is often regarded as the second most important enzyme in nature after photosystem II. The mechanism for how nitrogenase is able to perform the difficult task of cleaving the strong bond in N2 is debated. It is known that for every electron that is donated to N2, two ATP are hydrolyzed. In the experimentally suggested mechanism, the activation occurs after four reductions of the ground state, but there is no suggestion for how the enzyme uses the hydrolysis energy to perform catalysis. In the theoretical mechanism, it is suggested that hydrolysis is used to reduce the electron donor. In previous papers, the steps leading to the activation of N2 in the so-called E4 state has been investigated, using both the experimental and theoretical mechanism, showing that only the theoretical one leads to agreement with EPR observations for E4. In the present paper, the four steps following E4, leading to the release of two ammonia molecules, are described using the same methodology as used in the previous studies.

Mechanisms for Methane and Ammonia Oxidation by Particulate Methane Monooxygenase

Particulate MMO (pMMO) catalyzes the oxidation of methane to methanol and also ammonia to hydroxylamine. Experimental characterization of the active site has been very difficult partly because the enzyme is membrane-bound. However, recently, there has been major progress mainly through the use of cryogenic electron microscopy (cryoEM). Electron paramagnetic resonance (EPR) and X-ray spectroscopy have also been employed. Surprisingly, the active site has only one copper. There are two histidine ligands and one asparagine ligand, and the active site is surrounded by phenyl alanines but no charged amino acids in the close surrounding. The present study is the first quantum chemical study using a model of that active site (CuD). Low barrier mechanisms have been found, where an important part is that there are two initial proton-coupled electron transfer steps to a bound O2 ligand before the substrate enters. Surprisingly, this leads to large radical character for the oxygens even though they are protonated. That result is very important for the ability to accept a proton from the substrates. Methods have been used which have been thoroughly tested for redox enzyme mechanisms.

Reduction of molecular oxygen in flavodiiron proteins - Catalytic mechanism and comparison to heme-copper oxidases

The family of flavodiiron proteins (FDPs) plays an important role in the scavenging and detoxification of both molecular oxygen and nitric oxide. Using electrons from a flavin mononucleotide cofactor molecular oxygen is reduced to water and nitric oxide is reduced to nitrous oxide and water. While the mechanism for NO reduction in FDPs has been studied extensively, there is very little information available about O2 reduction. Here we use hybrid density functional theory (DFT) to study the mechanism for O2 reduction in FDPs. An important finding is that a proton coupled reduction is needed after the O2 molecule has bound to the diferrous diiron active site and before the OO bond can be cleaved. This is in contrast to the mechanism for NO reduction, where both NN bond formation and NO bond cleavage occurs from the same starting structure without any further reduction, according to both experimental and computational results. This computational result for the O2 reduction mechanism should be possible to evaluate experimentally. Another difference between the two substrates is that the actual OO bond cleavage barrier is low, and not involved in rate-limiting the reduction process, while the barrier connected with bond cleavage/formation in the NO reduction process is of similar height as the rate-limiting steps. We suggest that these results may be part of the explanation for the generally higher activity for O2 reduction as compared to NO reduction in most FDPs. Comparisons are also made to the O2 reduction reaction in the family of heme‑copper oxidases.

DFT Mechanistic Investigation into Ni(II)-Catalyzed Hydroxylation of Benzene to Phenol by H2O2

Introduction of oxygen into aromatic C-H bonds is intriguing from both fundamental and practical perspectives. Although the 3d metal-catalyzed hydroxylation of arenes by H2O2 has been developed by several prominent researchers, a definitive mechanism for these crucial transformations remains elusive. Herein, density functional theory calculations were used to shed light on the mechanism of the established hydroxylation reaction of benzene with H2O2, catalyzed by [NiII(tepa)]2+ (tepa = tris[2-(pyridin-2-yl)ethyl]amine). Dinickel(III) bis(μ-oxo) species have been proposed as the key intermediate responsible for the benzene hydroxylation reaction. Our findings indicate that while the dinickel dioxygen species can be generated as a stable structure, it cannot serve as an active catalyst in this transformation. The calculations allowed us to unveil an unprecedented mechanism composed of six main steps as follows: (i) deprotonation of coordinated H2O2, (ii) oxidative addition, (iii) water elimination, (iv) benzene addition, (v) ketone generation, and (vi) tautomerization and regeneration of the active catalyst. Addition of benzene to oxygen, which occurs via a radical mechanism, turns out to be the rate-determining step in the overall reaction. This study demonstrates the critical role of Ni-oxyl species in such transformations, highlighting how the unpaired spin density value on oxygen and positive charges on the Ni-O• complex affect the activation barrier for benzene addition.

Scalar Relativistic All-Electron and Pseudopotential Ab Initio Study of a Minimal Nitrogenase [Fe(SH)4H]- Model Employing Coupled-Cluster and Auxiliary-Field Quantum Monte Carlo Many-Body Methods

Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to organisms. The detailed mechanism of this enzyme is currently not known, and computational studies are complicated by the fact that different density functional theory (DFT) methods give very different energetic results for calculations involving nitrogenase models. Recently, we designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S to mimic different possible protonation states of the nitrogenase active site. We showed that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical methods. Based on nonrelativistic single-reference coupled-cluster (CC) calculations, we estimated that the ΔE is 101 kJ/mol. In this study, we demonstrate that scalar relativistic effects play an important role and significantly affect ΔE. Our best revised single-reference CC estimates for ΔE are 85-91 kJ/mol, including energy corrections to account for contributions beyond triples, core-valence correlation, and basis-set incompleteness error. Among coupled-cluster approaches with approximate triples, the canonical CCSD(T) exhibits the largest error for this problem. Complementary to CC, we also used phaseless auxiliary-field quantum Monte Carlo calculations (ph-AFQMC). We show that with a Hartree-Fock (HF) trial wave function, ph-AFQMC reproduces the CC results within 5 ± 1 kJ/mol. With multi-Slater-determinant (MSD) trials, the results are 82-84 ± 2 kJ/mol, indicating that multireference effects may be rather modest. Among the DFT methods tested, τ-HCTH, r2SCAN with 10-13% HF exchange with and without dispersion, and O3LYP/O3LYP-D4, and B3LYP*/B3LYP*-D4 generally perform the best. The r2SCAN12 (with 12% HF exchange) functional mimics both the best reference MSD ph-AFQMC and CC ΔE results within 2 kJ/mol.

Computational Model Study of the Experimentally Suggested Mechanism for Nitrogenase

The mechanism for N2 activation in the E4 state of nitrogenase was investigated by model calculations. In the experimentally suggested mechanism, the E4 state is obtained after four reductions to the ground state. In a recent theoretical study, results for a different mechanism have been found in excellent agreement with available Electron Paramagnetic Resonance (EPR) experiments for E4. The two hydrides in E4 leave as H2 concertedly with the binding of N2. The mechanism suggested differs from the experimentally suggested one by a requirement for four activation steps prior to catalysis. In the present study, the experimentally suggested mechanism is studied using the same methods as those used in the previous study on the theoretical mechanism. The computed results make it very unlikely that a structure obtained after four reductions from the ground state has two hydrides, and the experimentally suggested mechanism does therefore not agree with the EPR experiments for E4. Another structure with only one hydride is here suggested to be the one that has been observed to bind N2 after only four reductions of the ground state.

The energetics of N2 reduction by vanadium containing nitrogenase

The main class of nitrogenases has a molybdenum in its cofactor. A mechanism for Mo-nitrogenase has recently been described. In the present study, another class of nitrogenases has been studied, the one with a vanadium instead of a molybdenum in its cofactor. It is generally believed that these classes use the same general mechanism to activate nitrogen. The same methodology has been used here as the one used for Mo-nitrogenase. N2 activation is known to occur after four reductions in the catalytic cycle, in the E4 state. The main features of the mechanism for Mo-nitrogenase found in the previous study are an activation process in four steps prior to catalysis, the release of a sulfide during the activation steps and the formation of H2 from two hydrides in E4, just before N2 is activated. The same features have been found here for V-nitrogenase. A difference is that five steps are needed in the activation process, which explains why the ground state of V-nitrogenase is a triplet (even number) and the one for Mo-nitrogenase is a quartet (odd number). The reason an additional step is needed for V-nitrogenase is that V3+ can be reduced to V2+, in contrast to the case for Mo3+ in Mo-nitrogenase. The fact that V3+ is Jahn-Teller active has important consequences. N2H2 is formed in E4 with reasonably small barriers.

The mechanism for N2 activation in the E4 - state of nitrogenase

Nitrogenases take nitrogen from the air and reduce it to ammonia. It has long been known that N2 becomes activated after four reductions in the catalytic cycle, in the E4 state. Several mechanisms for the activation have been suggested. In the present study a previous mechanism has been revised based on recent experimental findings. In the present mechanism N2H2 is formed in E4. As in the previously suggested mechanism, there are four initial reductions before catalysis (the A-states), after which a sulfide is released and the first state in catalysis (E0) is formed. In E4, N2 becomes bound and protonated in the Fe1, Fe2, Fe4 region, in which the hydrides have left two electrons. The rate-limiting step is the formation of N2H by a hydrogen atom transfer from Cys275 to N2 bound to Fe4, concerted with an additional electron transfer from the cofactor. The mechanism fulfills all requirements set by experiments. The activation of N2 is preceded by a formation of H2 from two hydrides, the carbide is kinetically hindered from being protonated, the E4 state is reversible. An important aspect is the presence of a water molecule in the Fe2, Fe6 region. The non-allowed formations of H2 from a hydride and a proton have been investigated and found to have higher barriers than the allowed formation of H2 from two hydrides.

Is There a Different Mechanism for Water Oxidation in Higher Plants?

The leading mechanism for the formation of O2 in photosystem II (PSII) has, during the past decade, been established as the so-called oxyl-oxo mechanism. In that mechanism, O2 is formed from a binding between an oxygen radical (oxyl) and a bridging oxo group. For the case of higher plants, that mechanism has recently been criticized. Instead, a nucleophilic attack of an oxo group on a five-coordinated Mn(V)═O group forming O2 has been suggested in a so-called water-unbound (WU) mechanism. In the present study, the WU mechanism has been investigated. It is found that the WU mechanism is just a variant of a previously suggested mechanism but with a reactant and a transition state that have much higher energies. The addition of a water molecule on the empty site of the Mn(V)═O center is very exergonic and leads back to the previously suggested oxyl-oxo mechanism.

The Quantum Chemical Cluster Approach in Biocatalysis

ConspectusThe quantum chemical cluster approach has been used for modeling enzyme active sites and reaction mechanisms for more than two decades. In this methodology, a relatively small part of the enzyme around the active site is selected as a model, and quantum chemical methods, typically density functional theory, are used to calculate energies and other properties. The surrounding enzyme is modeled using implicit solvation and atom fixing techniques. Over the years, a large number of enzyme mechanisms have been solved using this method. The models have gradually become larger as a result of the faster computers, and new kinds of questions have been addressed. In this Account, we review how the cluster approach can be utilized in the field of biocatalysis. Examples from our recent work are chosen to illustrate various aspects of the methodology. The use of the cluster model to explore substrate binding is discussed first. It is emphasized that a comprehensive search is necessary in order to identify the lowest-energy binding mode(s). It is also argued that the best binding mode might not be the productive one, and the full reactions for a number of enzyme-substrate complexes have therefore to be considered to find the lowest-energy reaction pathway. Next, examples are given of how the cluster approach can help in the elucidation of detailed reaction mechanisms of biocatalytically interesting enzymes, and how this knowledge can be exploited to develop enzymes with new functions or to understand the reasons for lack of activity toward non-natural substrates. The enzymes discussed in this context are phenolic acid decarboxylase and metal-dependent decarboxylases from the amidohydrolase superfamily. Next, the application of the cluster approach in the investigation of enzymatic enantioselectivity is discussed. The reaction of strictosidine synthase is selected as a case study, where the cluster calculations could reproduce and rationalize the selectivities of both the natural and non-natural substrates. Finally, we discuss how the cluster approach can be used to guide the rational design of enzyme variants with improved activity and selectivity. Acyl transferase from Mycobacterium smegmatis serves as an instructive example here, for which the calculations could pinpoint the factors controlling the reaction specificity and enantioselectivity. The cases discussed in this Account highlight thus the value of the cluster approach as a tool in biocatalysis. It complements experiments and other computational techniques in this field and provides insights that can be used to understand existing enzymes and to develop new variants with tailored properties.

How Protons Move in Enzymes─The Case of Nitrogenase

When moving protons in enzymes, water molecules are often used as intermediates. The water molecules used are not necessarily seen in the crystal structures if they move around at high rates. In a different situation, for metal containing cofactors in enzymes, it is sometimes necessary to move protons on the cofactor from the position they enter the cofactor to another position where the energy is lower. That is, for example, the situation in nitrogenase. In recent studies on that enzyme, prohibitively high barriers were sometimes found for transferring protons, and that was used as a strong argument against mechanisms where a sulfide is lost in the mechanism. A high barrier could be due to nonoptimal distances and angles at the transition state. In the present study, possibilities are investigated to use water molecules to reduce these barriers. The study is very general and could have been done for many other enzymes. The effect of water was found to be very large in the case of nitrogenase with a lowering of one barrier from 15.6 kcal/mol down to essentially zero. It is concluded that the effect of water molecules must be taken into account for meaningful results.

The Perfluoro Cage Effect: A Search for Electron-Encapsulating Molecules

Quantum chemical calculations have for some time predicted that perfluorinated polyhedral organic molecules should exhibit a low-energy LUMO consisting of the overlapping inward-pointing lobes of the C-F σ* orbitals. Accordingly, these molecules should be able to encapsulate an electron within the interior of their cavities. Inspired by the recent confirmation of this prediction for perfluorocubane, we have sought to identify additional perfluorinated cage molecules capable of this remarkable behavior, which we refer to as the perfluoro cage effect (PCE). Using DFT calculations with multiple well-tested exchange-correlation functionals and large STO-QZ4P basis sets, we have identified several systems including [n]prismanes (n = 3-6), [n]asteranes (n = 3-5), twistane, and two norbornadiene dimer cages that clearly exhibit the PCE. In other words, they exhibit a low-energy LUMO belonging to the total symmetric irreducible representation of the point group in question and adiabatic electron affinities ranging from somewhat under 1 eV to over 2 eV. A pronounced size effect appears to hold, with larger cages exhibiting higher electron affinities (EAs). The largest adiabatic EAs, well over 3 eV, are predicted for perfluorinated dodecahedrane and C₆₀. In contrast, the PCE is barely discernible for perfluorinated tetrahedrane and bicyclo[1.1.1]pentane.

Quantum Mechanical Study of Oxygen Ligands Protonation for the Stable States of the Laccase Active Site

Laccases are enzymes catalyzing the oxidation of a wide range of organic and inorganic substrates accompanied by molecular oxygen reduction to water. Recently, oxygen reduction by laccases has been studied by single-crystal serial X-ray crystallography with increasing absorption doses at subatomic resolution. There were two determined structures corresponding to the reduced and oxidized stable states of the laccase active site. However, the protonation of the oxygen ligands involved cannot be determined even at subatomic resolution. In the present work, the protonation of oxygen ligands in the active site of laccase for the two stable states determined in the X-ray study was explored using quantum mechanical and continuum-electrostatics calculations. This is important for understanding the reaction of the oxygen reduction mechanism in laccases. The high precision of X-ray data at subatomic resolutions allowed us to optimize the quantum mechanical calculations.

Can the E1 state in nitrogenase tell if there is an activation process prior to catalysis?

Model calculations have been performed for the singly reduced ground state of Mo-nitrogenase, usually termed E₁. Contradictory conclusions have been reached in two recent experimental studies. In a study based on EPR, it was concluded that there is a bridging hydride in E₁, while in an X-ray study it was concluded that there is no hydride in E₁. Therefore, the EPR study implies that there is an oxidation of the cofactor going from E₀ to E₁, the X-ray study implies a reduction. DFT methods have here been used, which have previously been benchmarked on a set of redox enzymes that led to the conclusion that the accuracy is about 3 kcal mol^-1 in all cases, even for redox transitions. The methodology should therefore be adequate for resolving the question of the hydride presence in E₁. As a comparison, calculations are performed on both Mo- and V-nitrogenase with the same conclusion. The conclusion from the calculations has far reaching consequences for the mechanism of nitrogenase.

Computational modeling of redox enzymes

A computational methodology is briefly described, which appears to be able to accurately describe the mechanisms of redox active enzymes. The method is built on hybrid density functional theory where the inclusion of a fraction of exact exchange is critical. Two examples of where the methodology has been applied are described. The first example is the mechanism for water oxidation in photosystem II, and the second one is the mechanism for N₂ activation by nitrogenase. The mechanism for PSII has obtained very strong support from subsequent experiments. For nitrogenase, the calculations suggest that there should be an activation process prior to catalysis, which is still strongly debated.

Porphyryne

Density functional theory calculations with the B3LYP*-D3 method with large STO-QZ4P basis sets unambiguously predict a singlet ground state for Zn-porphyryne. However, the calculations also predict a low singlet-triplet gap of about 0.4 eV and a high adiabatic electron affinity of 2.4 eV. Accordingly, the reactivity of porphyryne species may be dominated by electron transfer, hydrogen abstraction, and proton-coupled electron transfer processes.

Radical in the Peroxide-Produced F-Type Ferryl Form of Bovine Cytochrome c Oxidase

The reduction of O₂ in respiratory cytochrome c oxidases (CcO) is associated with the generation of the transmembrane proton gradient by two mechanisms. In one of them, the proton pumping, two different types of the ferryl intermediates of the catalytic heme a₃-Cu_B center P and F forms, participate. Equivalent ferryl states can be also formed by the reaction of the oxidized CcO (O) with H₂O₂. Interestingly, in acidic solutions a single molecule of H₂O₂ can generate from the O an additional F-type ferryl form (F^•) that should contain, in contrast to the catalytic F intermediate, a free radical at the heme a₃-Cu_B center. In this work, the formation and the endogenous decay of both the ferryl iron of heme a₃ and the radical in F^• intermediate were examined by the combination of four experimental approaches, isothermal titration calorimetry, electron paramagnetic resonance, and electronic absorption spectroscopy together with the reduction of this form by the defined number of electrons. The results are consistent with the generation of radicals in F^• form. However, the radical at the catalytic center is more rapidly quenched than the accompanying ferryl state of heme a₃, very likely by the intrinsic oxidation of the enzyme itself.

Energetics for CO2 Reduction by Molybdenum-Containing Formate Dehydrogenase

The level of carbon dioxide in the atmosphere has increased in a dangerous way during the past century. Methods to decrease this level are therefore of high interest at present. Inspiration to do so in an efficient way could come from biological systems. Molybdenum-containing formate dehydrogenase (Mo-FDH) is one of the most interesting enzymes in this respect. For example, the reduction potential required is not very low. The normal reaction catalyzed by Mo-FDH is actually the opposite one of oxidizing formate to CO2. However, recent electrochemical studies have shown that the reaction can be reversed by a moderate lowering of the reduction potential. The goal of the present study has been to study the full mechanism of Mo-FDH, particularly in the most interesting direction of reducing CO2, which has not been done before. The methods used are the same as those that have been shown to give excellent results for redox enzymes in all cases they have been tested. The results obtained for Mo-FDH are also in excellent agreement with the experimental results.

A Mechanism for Nitrogenase Including Loss of a Sulfide

Nitrogenase is the only enzyme in nature that can fix N₂ from the air. The active cofactor of the leading form of this enzyme contains seven irons and one molybdenum connected by sulfide bridges. In several recent experimental studies, it has been suggested that the cofactor is very flexible, and might lose one of its sulfides during catalysis. In this study, the possible loss of a sulfide has been investigated by model calculations. In previous studies, we have shown that there should be four activation steps before catalysis starts, and this study is based on that finding. It was found here that, after the four reductions in the activation steps, a sulfide will become very loosely bound and can be released in a quite exergonic step with a low barrier. The binding of N₂ has no part in that release. In our previous studies, we suggested that the central carbide should be protonated three times after the four activation steps. With the new finding, there will instead be a loss of a sulfide, as the barrier for the loss is much lower than the ones for protonating the carbide. Still, it is suggested here that the carbide will be protonated anyway, but only with one proton, in the E₃ to E₄ step. A very complicated transition state for H₂ formation involving a large structural change was obtained. The combined step, with a loss of H₂ and binding of N₂ , is calculated to be endergonic by +2.3 kcal mol^-1 ; this is in excellent agreement with experiments in which an easily reversible step has been found.

Structural Predictive Model of Presenilin-2 Protein and Analysis of Structural Effects of Familial Alzheimer's Disease Mutations

Alzheimer's disease manifests itself in brain tissue by neuronal death, due to aggregation of β-amyloid, produced by senile plaques, and hyperphosphorylation of the tau protein, which produces neurofibrillary tangles. One of the genetic markers of the disease is the gene that translates the presenilin-2 protein, which has mutations that favor the appearance of the disease and has no reported crystallographic structure. In view of this, protein modeling is performed using prediction and structural refinement tools followed by an energetic and stereochemical characterization for its validation. For the simulation, four reported mutations are chosen, which are Met239Ile, Met239Val, Ser130Leu, and Thr122Arg, all associated with various functional responses. From a theoretical analysis, a preliminary bioinformatic study is made to find the phosphorylation patterns in the protein and the hydropathic index according to the polarity and chemical environment. Molecular visualization was carried out with the Chimera 1.14 software, and the theoretical calculation with the hybrid quantum mechanics/molecular mechanics system from the semi-empirical method, with Spartan18 software and an AustinModel1 basis. These relationships allow for studying the system from a structural approach with the determination of small distance changes, potential surfaces, electrostatic maps, and angle changes, which favor the comparison between wild-type and mutant systems. With the results obtained, it is expected to complement experimental data reported in the literature from models that would allow us to understand the effects of the selected mutations.

The importance of exact exchange-A methodological investigation of NO reduction in heme-copper oxidases

Significant improvements of the density functional theory (DFT) methodology during the past few decades have made DFT calculations a powerful tool in studies of enzymatic reaction mechanisms. For metalloenzymes, however, there are still concerns about the reliability in the DFT-results. Therefore, a systematic study is performed where the fraction of exact exchange in a hybrid DFT functional is used as a parameter. By varying this parameter, a set of different but related functionals are obtained. The various functionals are applied to one of the reactions occurring in the enzyme family heme-copper oxidases, the reduction of nitric oxide (NO) to nitrous oxide (N₂O) and water. The results show that, even though certain parts of the calculated energetics exhibit large variations, the qualitative pictures of the reaction mechanisms are quite stable. Furthermore, it is found that the functional with 15% exact exchange (B3LYP^*) gives the best agreement with experimental data for the particular reactions studied. An important aspect of the procedure used is that the computational results are carefully combined with a few more general experimental data to obtain a complete description of the entire catalytic cycle of the reactions studied.

On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT

In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.

Protein Predictive Modeling and Simulation of Mutations of Presenilin-1 Familial Alzheimer's Disease on the Orthosteric Site

Alzheimer's disease pathology is characterized by β-amyloid plaques and neurofibrillary tangles. Amyloid precursor protein is processed by β and γ secretase, resulting in the production of β-amyloid peptides with a length ranging from 38 to 43 amino acids. Presenilin 1 (PS1) is the catalytic unit of γ-secretase, and more than 200 PS1 pathogenic mutations have been identified as causative for Alzheimer's disease. A complete monocrystal structure of PS1 has not been determined so far due to the presence of two flexible domains. We have developed a complete structural model of PS1 using a computational approach with structure prediction software. Missing fragments Met1-Glut72 and Ser290-Glu375 were modeled and validated by their energetic and stereochemical characteristics. Then, with the complete structure of PS1, we defined that these fragments do not have a direct effect in the structure of the pore. Next, we used our hypothetical model for the analysis of the functional effects of PS1 mutations Ala246GLu, Leu248Pro, Leu248Arg, Leu250Val, Tyr256Ser, Ala260Val, and Val261Phe, localized in the catalytic pore. For this, we used a quantum mechanics/molecular mechanics (QM/MM) hybrid method, evaluating modifications in the topology, potential surface density, and electrostatic potential map of mutated PS1 proteins. We found that each mutation exerts changes resulting in structural modifications of the active site and in the shape of the pore. We suggest this as a valid approach for functional studies of PS1 in view of the possible impact in substrate processing and for the design of targeted therapeutic strategies.

Activation of O2 and NO in heme-copper oxidases - mechanistic insights from computational modelling

Heme-copper oxidases are transmembrane enzymes involved in aerobic and anaerobic respiration. The largest subgroup contains the cytochrome c oxidases (CcO), which reduce molecular oxygen to water. A significant part of the free energy released in this exergonic process is conserved as an electrochemical gradient across the membrane, via two processes, electrogenic chemistry and proton pumping. A deviant subgroup is the cytochrome c dependent NO reductases (cNOR), which reduce nitric oxide to nitrous oxide and water. This is also an exergonic reaction, but in this case none of the released free energy is conserved. Computational studies applying hybrid density functional theory to cluster models of the bimetallic active sites in the heme-copper oxidases are reviewed. To obtain a reliable description of the reaction mechanisms, energy profiles of the entire catalytic cycles, including the reduction steps have to be constructed. This requires a careful combination of computational results with certain experimental data. Computational studies have elucidated mechanistic details of the chemical parts of the reactions, involving cleavage and formation of covalent bonds, which have not been obtainable from pure experimental investigations. Important insights regarding the mechanisms of energy conservation have also been gained. The computational studies show that the reduction potentials of the active site cofactors in the CcOs are large enough to afford electrogenic chemistry and proton pumping, i.e. efficient energy conservation. These results solve a conflict between different types of experimental data. A mechanism for the proton pumping, involving a specific and crucial role for the active site tyrosine, conserved in all CcOs, is suggested. For the cNORs, the calculations show that the low reduction potentials of the active site cofactors are optimized for fast elimination of the toxic NO molecules. At the same time, the low reduction potentials lead to endergonic reduction steps with high barriers. To prevent even higher barriers, which would lead to a too slow reaction, when the electrochemical gradient across the membrane is present, the chemistry must occur in a non-electrogenic manner. This explains why there is no energy conservation in cNOR.

The Redox-Active Tyrosine Is Essential for Proton Pumping in Cytochrome c Oxidase

Cellular respiration involves electron transport via a number of enzyme complexes to the terminal Cytochrome c oxidase (CcO), in which molecular oxygen is reduced to water. The free energy released in the reduction process is used to establish a transmembrane electrochemical gradient, via two processes, both corresponding to charge transport across the membrane in which the enzymes are embedded. First, the reduction chemistry occurring in the active site of CcO is electrogenic, which means that the electrons and protons are delivered from opposite sides of the membrane. Second, the exergonic chemistry is coupled to translocation of protons across the entire membrane, referred to as proton pumping. In the largest subfamily of the CcO enzymes, the A-family, one proton is pumped for every electron needed for the chemistry, making the energy conservation particularly efficient. In the present study, hybrid density functional calculations are performed on a model of the A-family CcOs. The calculations show that the redox-active tyrosine, conserved in all types of CcOs, plays an essential role for the energy conservation. Based on the calculations a reaction mechanism is suggested involving a tyrosyl radical (possibly mixed with tyrosinate character) in all reduction steps. The result is that the free energy released in each reduction step is large enough to allow proton pumping in all reduction steps without prohibitively high barriers when the gradient is present. Furthermore, the unprotonated tyrosine provides a mechanism for coupling the uptake of two protons per electron in every reduction step, i.e. for a secure proton pumping.

The active E4 structure of nitrogenase studied with different DFT functionals

The present study concerns the technical aspects of obtaining the energetics for the E4 state of nitrogenase, the enzyme that fixes N₂ in nature. EPR experiments have shown that the critical E4 structure that activates N₂ should contain two bridging hydrides in the FeMo-cofactor. It is furthermore in equilibrium with a structure where the two hydrides have been released and N₂ binds. These observations led to the suggestion that E4 should have two bridging hydrides and two protonated sulfides. It is important to note that the structure for E4 has not been determined, but only suggested. For a long time, no DFT study led to the suggested structure, independent of which functional was used. However, in two recent DFT studies a good agreement with the experimental suggestion was claimed to have been obtained. In one of them the TPSS functional was used. That was the first out of 11 functionals tried that led to the experimentally suggested structure. In the second of the recent DFT studies, a similar conclusion was reached using the TPSSh functional. The conclusions in the recent studies have here been studied in detail, by calculating a critical energetic value strongly implied by the same EPR experiments. Both the TPSS and TPSSh functionals have been used. The present calculations suggest that those DFT functionals would not lead to agreement with the experimental EPR results either.

A quantum chemical approach for the mechanisms of redox-active metalloenzymes

During the past 20 years, quantum chemistry has grown to be a significant part in the investigation of mechanisms for redox-active enzymes. In our group we have developed an approach that has been applied to a large number of such systems. Hybrid density functional theory (hybrid DFT) has from the start of these investigations been the leading electronic structure tool. An understanding of how the method works in practice has significantly improved the accuracy and applicability. During the past ten years, it has been found that the results for redox enzymes mainly depend on the chosen fraction of exact exchange in the functional, and that a choice of 15% has worked best. The idea has therefore been to vary that fraction over a reasonable range and study the relative energy dependence. For modeling the enzymes, a cluster approach has been developed. In the present review the development of the method we used is described from its start in work on photosystem II, fifteen years ago. Examples from a few recent applications are described, where the metals have been iron, nickel, copper, cobalt or manganese. The results are in excellent agreement with available experiments, and a large number of new predictions have been made.

During the past 20 years, quantum chemistry has grown to be a significant part in the investigation of mechanisms for redox-active enzymes.

Energetics for Proton Reduction in FeFe Hydrogenase

The energetics for proton reduction in FeFe-hydrogenase has been reinvestigated by theoretical modeling, in light of recent experiments. Two different mechanisms have been considered. In the first one, the bridging hydride position was blocked by the enzyme, which is the mechanism that has been supported by a recent spectroscopic study by Cramer et al. A major difficulty in the present study to agree with experimental energetics was to find the right position for the added proton in the first reduction step. It was eventually found that the best position was as a terminal hydride on the distal iron, which has not been suggested in any of the recent, experimentally based mechanisms. The lowest transition state was surprisingly found to be a bond formation between a proton on a cysteine and the terminal hydride. This type of TS is similar to the one for heterolytic H₂ cleavage in NiFe hydrogenase. The second mechanism investigated here is not supported by the present calculations or the recent experiments by Cramer et al., but was still studied as an interesting comparison. In that mechanism, the formation of the bridging hydride was allowed. The H–H formation barrier is only 3.6 kcal/mol higher than for the first mechanism, but there are severe problems concerning the motion of the protons.

From Synthetic to Biological Fe4 S4 Complexes: Redox Properties Correlated to Function of Radical S-Adenosylmethionine Enzymes

By employing the computational protocol for calculation of reduction potentials of the Fe₄ S₄ -containing species validated using a representative series of well-defined synthetic complexes, we focused on redox properties of two prototypical radical SAM enzymes to reveal how they transform SAM into the reactive 5'-deoxyadenosyl radical, and how they tune this radical for its proper biological function. We found the reduction potential of SAM is indeed elevated by 0.3-0.4 V upon coordination to Fe₄ S₄ , which was previously speculated in the literature. This makes a generation of 5'-deoxyadenosyl radical from SAM less endergonic (by ca. 7-9 kcal mol^-1 ) and hence more feasible in both enzymes as compared to the identical process in water. Furthermore, our calculations indicate that the enzyme-bound 5'-deoxyadenosyl radical has a significantly lower reduction potential than in referential aqueous solution, which may help the enzymes to suppress potential side redox reactions and simultaneously elevate its proton-philic character, which may, in turn, promote the radical hydrogen-atom abstraction ability.

A Theoretical Study of the Recently Suggested MnVII Mechanism for O-O Bond Formation in Photosystem II

The mechanism for water oxidation in photosystem II has been a major topic for several decades. The active catalyst has four manganese atoms connected by bridging oxo bonds, in a complex termed the oxygen-evolving complex (OEC), which also includes a calcium atom. The O-O bond of oxygen is formed after absorption of four photons in a state of the OEC termed S₄. There has been essential consensus that in the S₄ state, all manganese atoms are in the Mn(IV) oxidation state. However, recently there has been a suggestion that one of the atoms is in the Mn(VII) state. In the present computational study, the feasibility of that proposal has been investigated. It is here shown that the mechanism involving Mn(VII) has a much higher barrier for forming O₂ than the previous proposal with four Mn(IV) atoms.

Theoretical Study of O2 Reduction and Water Oxidation in Multicopper Oxidases

Recent electrochemical experiments have shown that the reduction of O₂ can be driven backward to water oxidation, which is the first case that has been successfully demonstrated for an enzyme. To understand this ability of the enzyme, both the forward reduction and backward oxidation have been studied here. For the forward reaction, a mechanism similar to earlier studies was obtained. All steps of the full catalytic cycle were obtained for the first time, and it was shown that the explicit reduction steps contribute significantly to the rate-limiting step of the O–O bond cleavage. For the backward oxidation reaction, it was found that the mechanism of the O–O bond formation is not just the reverse of the reduction step where the O–O bond is cleaved for a protonated peroxide. The formation of two fully deprotonated oxo groups was found to be important, which leads to a large radical character for one of the oxo groups. For this possibility, it is important that the pK_a of the water bound to the cofactor is quite high.

Theoretical rationalization for the equilibrium between (μ-η2:η2-peroxido)CuIICuII and bis(μ-oxido)CuIIICuIII complexes: perturbational effects from ligand frameworks

DFT calculations are carried out to investigate the geometric effects of the supporting ligands in the relative energies of the (μ-η2:η2-peroxido)CuIICuII complex 1 and the bis(μ-oxido)CuIIICuIII complex 2. The N3-tridentate ligand bearing acyclic propane diamine framework La preferentially provided 1, whereas the N3-tridentate ligand with cyclic diamine framework such as 1,4-diazacycloheptane Lb gave 2 after the oxygenation of the corresponding CuI complexes as reported previously [S. Itoh, et al., Inorg. Chem., 2014, 53, 8786-8794]. Calculations at the B3LYP*-D3 level of theory can reasonably explain the experimental results in relative energies, structures and harmonic frequencies of 1 and 2. Perturbational effects of the diamine chelates of La and Lb especially on the equilibrium of 1 and 2 are investigated in detail. In the range from 2.30 Å to 3.40 Å of the N-N distance in the diamine moiety, 1 is more stable than 2 by 8.4 kcal mol-1 at the distance of 3.40 Å. Calculated potential energies indicate that the decrease in the N-N distance is associated with a decrease in energy of 2, leading that 2 can be most stabilized at the N-N distance of 2.60 Å. Furthermore, molecular orbitals analyses are performed to explain that the energy gaps between the σ* orbital of the O-O bond and the dx2-y2 orbitals of the CuII ions of 1 get small as the diamine moiety is shrunk, leading to facilitate the O-O bond cleavage from 1 to 2.

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.

Characterisation of the Cyanate Inhibited State of Cytochrome c Oxidase

Heme-copper oxygen reductases are terminal respiratory enzymes, catalyzing the reduction of dioxygen to water and the translocation of protons across the membrane. Oxygen consumption is inhibited by various substances. Here we tested the relatively unknown inhibition of cytochrome c oxidase (CcO) with isocyanate. In contrast to other more common inhibitors like cyanide, inhibition with cyanate was accompanied with the rise of a metal to ligand charge transfer (MLCT) band around 638 nm. Increasing the cyanate concentration furthermore caused selective reduction of heme a. The presence of the CT band allowed for the first time to directly monitor the nature of the ligand via surface-enhanced resonance Raman (SERR) spectroscopy. Analysis of isotope sensitive SERR spectra in comparison with Density Functional Theory (DFT) calculations identified not only the cyanate monomer as an inhibiting ligand but suggested also presence of an uretdion ligand formed upon dimerization of two cyanate ions. It is therefore proposed that under high cyanate concentrations the catalytic site of CcO promotes cyanate dimerization. The two excess electrons that are supplied from the uretdion ligand lead to the observed physiologically inverse electron transfer from heme a₃ to heme a.

What Is the Structure of the E4 Intermediate in Nitrogenase?

Nitrogenase is the only enzyme that can cleave the strong triple bond in N₂. The active site contains a complicated MoFe₇S₉C cluster. It is believed that it needs to accept four protons and electrons, forming the E₄ state, before it can bind N₂. However, there is no consensus on the atomic structure of the E₄ state. Experimental studies indicate that it should contain two hydride ions bridging two pairs of Fe ions, and it has been suggested that both hydride ions as well as the two protons bind on the same face of the cluster. On the other hand, density functional theory (DFT) studies have indicated that it is energetically more favorable with either three hydride ions or with a triply protonated carbide ion, depending on the DFT functional. We have performed a systematic combined quantum mechanical and molecular mechanical (QM/MM) study of possible E₄ states with two bridging hydride ions. Our calculations suggest that the most favorable structure has hydride ions bridging the Fe2/6 and Fe3/7 ion pairs. In fact, such structures are 14 kJ/mol more stable than structures with three hydride ions, showing that pure DFT functionals give energetically most favorable structures in agreement with experiments. An important reason for this finding is that we have identified a new type of broken-symmetry state that involves only two Fe ions with minority spin, in contrast to the previously studied states with three Fe ions with minority spin. The energetically best structures have the two hydride ions on different faces of the FeMo cluster, whereas better agreement with ENDOR data is obtained if they are on the same face; such structures are only 6–22 kJ/mol less stable.

The mechanism for oxygen reduction in the C family cbb3 cytochrome c oxidases - Implications for the proton pumping stoichiometry

Cytochrome c oxidases (CcOs) couple the exergonic reduction of molecular oxygen to proton pumping across the membrane in which they are embedded, thereby conserving a significant part of the free energy. The A family CcOs are known to pump four protons per oxygen molecule, while there is no consensus regarding the proton pumping stoichiometry for the C family cbb₃ oxidases. Hybrid density functional theory is used here to investigate the catalytic mechanism for oxygen reduction in cbb₃ oxidases. A surprising result is that the barrier for O O bond cleavage at the mixed valence reduction level seems to be too high compared to the overall reaction rate of the enzyme. It is therefore suggested that the O O bond is cleaved only after the first proton coupled reduction step, and that this reduction step most likely is not coupled to proton pumping. Furthermore, since the cbb₃ oxidases have only one proton channel leading to the active site, it is proposed that the activated E_H intermediate, suggested to be responsible for proton pumping in one of the reduction steps in the A family, cannot be involved in the catalytic cycle for cbb₃, which results in the lack of proton pumping also in the E to R reduction step. In summary, the calculations indicate that only two protons are pumped per oxygen molecule in cbb₃ oxidases. However, more experimental information on this divergent enzyme is needed, e.g. whether the flow of electrons resembles that in the other more well-studied CcO families.

Energetics for the Mechanism of Nickel-Containing Carbon Monoxide Dehydrogenase

Nickel-containing carbon monoxide (CO) dehydrogenase is an enzyme that catalyzes the important reversible carbon dioxide reduction. Several high-resolution structures have been determined at various stages of the reduction, which can be used as good starting points for the present computational study. The cluster model is used in combination with a systematic application of the density functional theory as recently described. The results are in very good agreement with experimental evidence. There are a few important results. To explain why the X-ray structure for the reduced C_red1 state has an empty site on nickel, it is here suggested that the cluster has been over-reduced by X-rays and is therefore not the desired reduced state, which instead contains a bound CO on nickel. After an additional reduction, a hydride bound to nickel is suggested to play a role. In order to obtain energetics in agreement with experiments, it is concluded that one sulfide bridge in the Ni-Fe cluster should be protonated. The best test of the accuracy obtained is to compare the computed rate for reduction using -0.6 V with that for oxidation using -0.3 V, where good agreement was obtained. Obtaining a mechanism that is easily reversible is another demanding aspect of the modeling. Nickel oscillates between nickel(II) and nickel(I), while nickel(0) never comes in.

The mechanism for nitrogenase including all steps

Nitrogen in the air is turned into biologically useful ammonia by the nitrogenase enzyme. The leading member of this group has a cofactor with one molybdenum and seven irons linked together by sulfurs. The structure that binds N₂ has a triply protonated carbide and a rotated homocitrate. Both these structural changes are necessary for the activation.