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 (N2O) 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.

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.

What can we learn from N2 O isotope data? - Analytics, processes and modelling.

The isotopic composition of nitrous oxide (N2 O) provides useful information for evaluating N2 O sources and budgets. Due to the co-occurrence of multiple N2 O transformation pathways, it is, however, challenging to use isotopic information to quantify the contribution of distinct processes across variable spatiotemporal scales. Here, we present an overview of recent progress in N2 O isotopic studies and provide suggestions for future research, mainly focusing on: analytical techniques; production and consumption processes; and interpretation and modelling approaches. Comparing isotope-ratio mass spectrometry (IRMS) with laser absorption spectroscopy (LAS), we conclude that IRMS is a precise technique for laboratory analysis of N2 O isotopes, while LAS is more suitable for in situ/inline studies and offers advantages for site-specific analyses. When reviewing the link between the N2 O isotopic composition and underlying mechanisms/processes, we find that, at the molecular scale, the specific enzymes and mechanisms involved determine isotopic fractionation effects. In contrast, at plot-to-global scales, mixing of N2 O derived from different processes and their isotopic variability must be considered. We also find that dual isotope plots are effective for semi-quantitative attribution of co-occurring N2 O production and reduction processes. More recently, process-based N2 O isotopic models have been developed for natural abundance and 15 N-tracing studies, and have been shown to be effective, particularly for data with adequate temporal resolution. Despite the significant progress made over the last decade, there is still great need and potential for future work, including development of analytical techniques, reference materials and inter-laboratory comparisons, further exploration of N2 O formation and destruction mechanisms, more observations across scales, and design and validation of interpretation and modelling approaches. Synthesizing all these efforts, we are confident that the N2 O isotope community will continue to advance our understanding of N2 O transformation processes in all spheres of the Earth, and in turn to gain improved constraints on regional and global budgets.

Role of the Two Metals in the Active Sites of Heme Copper Oxidases-A Study of NO Reduction in cbb3 Cytochrome c Oxidase.

The superfamily of heme copper oxidases reduces molecular oxygen or nitric oxide, and the active sites comprise a high-spin heme group (a3 or b3) and a non-heme metal (CuB or FeB). The cbb3 C family of cytochrome c oxidases, with the high-spin heme b3 and CuB in the active site, is a subfamily of the heme copper oxidases that can reduce both molecular oxygen, which is the main substrate, and nitric oxide. The mechanism for NO reduction in cbb3 oxidase is studied here using hybrid density functional theory and compared to other cytochrome c oxidases (A and B families), with a high-spin heme a3 and CuB in the active site, and to cytochrome c dependent NO reductase, with a high-spin heme b3 and a non-heme FeB in the active site. It is found that the reaction mechanism and the detailed reaction energetics of the cbb3 oxidases are not similar to those of cytochrome c dependent NO reductase, which has the same type of high-spin heme group but a different non-heme metal. This is in contrast to earlier expectations. Instead, the NO reduction mechanism in cbb3 oxidases is very similar to that in the other cytochrome c oxidases, with the same non-heme metal, CuB, and is independent of the type of high-spin heme group. The conclusion is that the type of non-heme metal (CuB or FeB) in the active site of the heme copper oxidases is more important for the reaction mechanisms than the type of high-spin heme, at least for the NO reduction reaction. The reason is that the proton-coupled reduction potentials of the active site cofactors determine the energetics for the NO reduction reaction, and they depend to a larger extent on the non-heme metal. Observed differences in NO reduction reactivity among the various cytochrome c oxidases may be explained by differences outside the BNC, affecting the rate of proton transfer, rather than in the BNC itself.

Active Site Midpoint Potentials in Different Cytochrome c Oxidase Families: A Computational Comparison.

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain, reducing molecular oxygen to water. The binuclear active site in C cO comprises a high-spin heme associated with a CuB complex and a redox active tyrosine. The electron transport in the respiratory chain is driven by increasing midpoint potentials of the involved cofactors, resulting in a release of free energy, which is stored by coupling the electron transfer to proton translocation across a membrane, building up an electrochemical gradient. In this context, the midpoint potentials of the active site cofactors in the C cOs are of special interest, since they determine the driving forces for the individual oxygen reduction steps and thereby affect the efficiency of the proton pumping. It has been difficult to obtain useful information on some of these midpoint potentials from experiments. However, since each of the reduction steps in the catalytic cycle of oxygen reduction to water corresponds to the formation of an O-H bond, they can be calculated with a reasonably high accuracy using quantum chemical methods. From the calculated O-H bond strengths, the proton-coupled midpoint potentials of the active site cofactors can be estimated. Using models representing the different families of C cO's (A, B, and C), the calculations give midpoint potentials that should be relevant during catalytic turnover. The calculations also suggest possible explanations for why some experimentally measured potentials deviate significantly from the calculated ones, i.e., for CuB in all oxidase families, and for heme b3 in the C family.

Can Reduction of NO to N2O in Cytochrome c Dependent Nitric Oxide Reductase Proceed through a Trans-Mechanism?

As part of microbial denitrification, NO is reduced to N2O in the membrane bound enzyme nitric oxide reductase, NOR. The N-N coupling occurs in the diiron binuclear active site, BNC, and different mechanisms for this reaction step have been suggested. Computational studies have supported a so-called cis:b3-mechanism, in which the hyponitrite product of the reductive N-N bond formation coordinates with one nitrogen to the heme iron and with both oxygens to the non-heme iron in the BNC. In contrast, experimental results have been interpreted to support a so-called trans-mechanism, in which the hyponitrite intermediate coordinates with one nitrogen atom to each of the two iron ions. Hybrid density functional theory is used here to perform an extensive search for possible intermediates of the NO reduction in the cNOR enzyme. It is found that hyponitrite structures coordinating with their negatively charged oxygens to the positively charged iron ions are the most stable ones. The hyponitrite intermediate involved in the suggested trans-mechanism, which only coordinates with the nitrogens to the iron ions, is found to be prohibitively high in energy, leading to a too slow reaction, which should rule out this mechanism. Furthermore, intermediates binding one NO molecule to each iron ion in the BNC, which have been suggested to initiate the trans-mechanism, are found to be too high in energy to be observable, indicating that the experimentally observed electron paramagnetic resonance signals, taken to support such an iron-nitrosyl dimer intermediate, should be reinterpreted.

The mechanism for oxygen reduction in cytochrome c dependent nitric oxide reductase (cNOR) as obtained from a combination of theoretical and experimental results.

Bacterial NO-reductases (NOR) belong to the heme-copper oxidase (HCuO) superfamily, in which most members are O2-reducing, proton-pumping enzymes. This study is one in a series aiming to elucidate the reaction mechanisms of the HCuOs, including the mechanisms for cellular energy conservation. One approach towards this goal is to compare the mechanisms for the different types of HCuOs, cytochrome c oxidase (CcO) and NOR, reducing the two substrates O2 and NO. Specifically in this study, we describe the mechanism for oxygen reduction in cytochrome c dependent NOR (cNOR). Hybrid density functional calculations were performed on large cluster models of the cNOR binuclear active site. Our results are used, together with published experimental information, to construct a free energy profile for the entire catalytic cycle. Although the overall reaction is quite exergonic, we show that during the reduction of molecular oxygen in cNOR, two of the reduction steps are endergonic with high barriers for proton uptake, which is in contrast to oxygen reduction in CcO, where all reduction steps are exergonic. This difference between the two enzymes is suggested to be important for their differing capabilities for energy conservation. An additional result from this study is that at least three of the four reduction steps are initiated by proton transfer to the active site, which is in contrast to CcO, where electrons always arrive before the protons to the active site. The roles of the non-heme metal ion and the redox-active tyrosine in the active site are also discussed.

Mechanism of Oxygen Reduction in Cytochrome c Oxidase and the Role of the Active Site Tyrosine.

Cytochrome c oxidase, the terminal enzyme in the respiratory chain, reduces molecular oxygen to water and stores the released energy through electrogenic chemistry and proton pumping across the membrane. Apart from the heme-copper binuclear center, there is a conserved tyrosine residue in the active site (BNC). The tyrosine delivers both an electron and a proton during the O-O bond cleavage step, forming a tyrosyl radical. The catalytic cycle then occurs in four reduction steps, each taking up one proton for the chemistry (water formation) and one proton to be pumped. It is here suggested that in three of the reduction steps the chemical proton enters the center of the BNC, leaving the tyrosine unprotonated with radical character. The reproprotonation of the tyrosine occurs first in the final reduction step before binding the next oxygen molecule. It is also suggested that this reduction mechanism and the presence of the tyrosine are essential for the proton pumping. Density functional theory calculations on large cluster models of the active site show that only the intermediates with the proton in the center of the BNC and with an unprotonated tyrosyl radical have a high electron affinity of similar size as the electron donor, which is essential for the ability to take up two protons per electron and thus for the proton pumping. This type of reduction mechanism is also the only one that gives a free energy profile in accordance with experimental observations for the amount of proton pumping in the working enzyme.

Protonation of the binuclear active site in cytochrome c oxidase decreases the reduction potential of CuB.

One of the remaining mysteries regarding the respiratory enzyme cytochrome c oxidase is how proton pumping can occur in all reduction steps in spite of the low reduction potentials observed in equilibrium titration experiments for two of the active site cofactors, CuB(II) and Fea3(III). It has been speculated that, at least the copper cofactor can acquire two different states, one metastable activated state occurring during enzyme turnover, and one relaxed state with lower energy, reached only when the supply of electrons stops. The activated state should have a transiently increased CuB(II) reduction potential, allowing proton pumping. The relaxed state should have a lower reduction potential, as measured in the titration experiments. However, the structures of these two states are not known. Quantum mechanical calculations show that the proton coupled reduction potential for CuB is inherently high in the active site as it appears after reaction with oxygen, which explains the observed proton pumping. It is suggested here that, when the flow of electrons ceases, a relaxed resting state is formed by the uptake of one extra proton, on top of the charge compensating protons delivered in each reduction step. The extra proton in the active site decreases the proton coupled reduction potential for CuB by almost half a volt, leading to agreement with titration experiments. Furthermore, the structure for the resting state with an extra proton is found to have a hydroxo-bridge between CuB(II) and Fea3(III), yielding a magnetic coupling that can explain the experimentally observed EPR silence.

How cytochrome c oxidase can pump four protons per oxygen molecule at high electrochemical gradient.

Experiments have shown that the A-family cytochrome c oxidases pump four protons per oxygen molecule, also at a high electrochemical gradient. This has been considered a puzzle, since two of the reduction potentials involved, Cu(II) and Fe(III), were estimated from experiments to be too low to afford proton pumping at a high gradient. The present quantum mechanical study (using hybrid density functional theory) suggests a solution to this puzzle. First, the calculations show that the charge compensated Cu(II) potential for CuB is actually much higher than estimated from experiment, of the same order as the reduction potentials for the tyrosyl radical and the ferryl group, which are also involved in the catalytic cycle. The reason for the discrepancy between theory and experiment is the very large uncertainty in the experimental observations used to estimate the equilibrium potentials, mainly caused by the lack of methods for direct determination of reduced CuB. Second, the calculations show that a high energy metastable state, labeled EH, is involved during catalytic turnover. The EH state mixes the low reduction potential of Fe(III) in heme a3 with another, higher potential, here suggested to be that of the tyrosyl radical, resulting in enough exergonicity to allow proton pumping at a high gradient. In contrast, the corresponding metastable oxidized state, OH, is not significantly higher in energy than the resting state, O. Finally, to secure the involvement of the high energy EH state it is suggested that only one proton is taken up via the K-channel during catalytic turnover.

Improved free energy profile for reduction of NO in cytochrome c dependent nitric oxide reductase (cNOR).

Quantum chemical calculations play an essential role in the elucidation of reaction mechanisms for redox-active metalloenzymes. For example, the cleavage and the formation of covalent bonds can usually not be described only on the basis of experimental information, but can be followed by the calculations. Conversely, there are properties, like reduction potentials, which cannot be accurately calculated. Therefore, computational and experimental data has to be carefully combined to obtain reliable descriptions of entire catalytic cycles involving electron and proton uptake from donors outside the enzyme. Such a procedure is illustrated here, for the reduction of nitric oxide (NO) to nitrous oxide and water in the membrane enzyme, cytochrome c dependent nitric oxide reductase (cNOR). A surprising experimental observation is that this reaction is nonelectrogenic, which means that no energy is conserved. On the basis of hybrid density functional calculations a free energy profile for the entire catalytic cycle is obtained, which agrees much better with experimental information on the active site reduction potentials than previous ones. Most importantly the energy profile shows that the reduction steps are endergonic and that the entire process is rate-limited by high proton uptake barriers during the reduction steps. This result implies that, if the reaction were electrogenic, it would become too slow when the gradient is present across the membrane. This explains why this enzyme does not conserve any of the free energy released. © 2016 Wiley Periodicals, Inc.

Proton pumping in cytochrome c oxidase: energetic requirements and the role of two proton channels.

Cytochrome c oxidase is a superfamily of membrane bound enzymes catalyzing the exergonic reduction of molecular oxygen to water, producing an electrochemical gradient across the membrane. The gradient is formed both by the electrogenic chemistry, taking electrons and protons from opposite sides of the membrane, and by proton pumping across the entire membrane. In the most efficient subfamily, the A-family of oxidases, one proton is pumped in each reduction step, which is surprising considering the fact that two of the reduction steps most likely are only weakly exergonic. Based on a combination of quantum chemical calculations and experimental information, it is here shown that from both a thermodynamic and a kinetic point of view, it should be possible to pump one proton per electron also with such an uneven distribution of the free energy release over the reduction steps, at least up to half the maximum gradient. A previously suggested pumping mechanism is developed further to suggest a reason for the use of two proton transfer channels in the A-family. Since the rate of proton transfer to the binuclear center through the D-channel is redox dependent, it might become too slow for the steps with low exergonicity. Therefore, a second channel, the K-channel, where the rate is redox-independent is needed. A redox-dependent leakage possibility is also suggested, which might be important for efficient energy conservation at a high gradient. A mechanism for the variation in proton pumping stoichiometry over the different subfamilies of cytochrome oxidase is also suggested. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Why is the reduction of NO in cytochrome c dependent nitric oxide reductase (cNOR) not electrogenic?

The membrane-bound enzyme cNOR (cytochrome c dependent nitric oxide reductase) catalyzes the reduction of NO in a non-electrogenic process. This is in contrast to the reduction of O2 in cytochrome c oxidase (CcO), the other member of the heme-copper oxidase family, which stores energy by the generation of a membrane gradient. This difference between the two enzymes has not been understood, but it has been speculated to be of kinetic origin, since per electron the NO reduction is more exergonic than the O2 reduction, and the energy should thus be enough for an electrogenic process. However, it has not been clear how and why electrogenicity, which mainly affects the thermodynamics, would slow down the very exergonic NO reduction. Quantum chemical calculations are used to construct a free energy profile for the catalytic reduction of NO in the active site of cNOR. The energy profile shows that the reduction of the NO molecules by the enzyme and the formation of N2O are very exergonic steps, making the rereduction of the enzyme endergonic and rate-limiting for the entire catalytic cycle. Therefore the NO reduction cannot be electrogenic, i.e. cannot take electrons and protons from the opposite sides of the membrane, since it would increase the endergonicity of the rereduction when the gradient is present, thereby increasing the rate-limiting barrier, and the reaction would become too slow. It also means that proton pumping coupled to electron transfer is not possible in cNOR. In CcO the corresponding rereduction of the enzyme is very exergonic.

An investigation of possible competing mechanisms for Ni-containing methyl-coenzyme M reductase.

Ni-containing methyl-coenzyme M reductase (MCR) is capable of catalyzing methane formation from methyl-coenzyme M (CH3-SCoM) and coenzyme B (CoB-SH), and also its reverse reaction (methane oxidation). Based on extensive experimental and theoretical investigations, it has turned out that a mechanism including an organometallic methyl-Ni(III)F430 intermediate is inaccessible, while another mechanism involving a methyl radical and a Ni(II)-SCoM species currently appears to be the most acceptable one for MCR. In the present paper, using hybrid density functional theory and an active-site model based on the X-ray crystal structure, two other mechanisms were studied and finally also ruled out. One of them, involving proton binding on the CH3-SCoM substrate, which should facilitate methyl-Ni(III)F430 formation, is demonstrated to be quite unfavorable since the substrate has a much smaller proton affinity than the F430 cofactor. Another one (oxidative addition mechanism) is also shown to be unfavorable for the MCR reaction, due to the large endothermicity for the formation of the ternary intermediate with side-on C-S (for CH3-SCoM) or C-H (for methane) coordination to Ni.

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.

The mechanism for proton pumping in cytochrome c oxidase from an electrostatic and quantum chemical perspective.

The mechanism for proton pumping in cytochrome c oxidase in the respiratory chain, has for decades been one of the main unsolved problems in biochemistry. However, even though several different suggested mechanisms exist, many of the steps in these mechanisms are quite similar and constitute a general consensus framework for discussing proton pumping. When these steps are analyzed, at least three critical gating situations are found, and these points are where the suggested mechanisms in general differ. The requirements for gating are reviewed and analyzed in detail, and a mechanism is suggested, where solutions for all the gating situations are formulated. This mechanism is based on an electrostatic analysis of a kinetic experiment fior the O to E transition. The key component of the mechanism is a positively charged transition state. An electron on heme a opens the gate for proton transfer from the N-side to a pump loading site (PLS). When the negative charge of the electron is compensated by a chemical proton, the positive transition state prevents backflow from the PLS to the N-side at the most critical stage of the pumping process. The mechanism has now been tested by large model DFT calculations, and these calculations give strong support for the suggested mechanism.

Theoretical study of the oxidation of phenolates by the [Cu2O2(N,N'-di-tert-butylethylenediamine)2]2+ complex.

Experiments have shown that the µ-η(2):η(2)-peroxodicopper(II) complex [Cu(2)O(2)(N,N'-di-tert-butylethylenediamine)(2)](2+) rapidly oxidizes 2,4-di-tert-butylphenolate into a mixture of catechol and quinone and that, at the extreme temperature of -120 °C, a bis-µ-oxodicopper(III)-phenolate intermediate, labeled complex A, can be observed. These experimental results suggest a new mechanism of action for the dinuclear copper-containing enzyme tyrosinase, involving an early O-O bond-cleavage step. However, whether phenolate binding occurs before or after the cleavage of the O-O bond has not been possible to answer. In this study, hybrid density functional theory is used to study the synthetic reaction and, based on the calculated free-energy profile, a mechanism is suggested for the entire phenolate-oxidation reaction that agrees with the experimental observations. Most importantly, the calculations show that the very first step in the reaction is the cleavage of the O-O bond in the peroxo complex and that, subsequently, the phenolate substrate coordinates to one of the copper ions in the bis-µ-oxodicopper(III) complex to yield the experimentally characterized phenolate intermediate (A). The oxidation of the phenolate substrate into a quinone then occurs in three steps: 1) C-O bond formation, 2) coupled internal proton and electron transfer, and 3) electron transfer coupled to proton transfer from an external donor (acidic workup, experimentally). The first of these steps is rate limiting for the decay of complex A, with a calculated free-energy barrier of 10.7 kcal mol(-1) and a deuterium kinetic isotope effect of 0.90, which are in good agreement with the experimental values of 11.2 kcal mol(-1) and 0.83(±0.09). The tert-butyl substituents on both the phenol substrate and the copper ligands need to be included in the calculations to give a correct description of the reaction mechanism.

Mechanism for N2O generation in bacterial nitric oxide reductase: a quantum chemical study.

The catalytic mechanism of reduction of NO to N(2)O in the bacterial enzyme nitric oxide reductase has been investigated using hybrid density functional theory and a model of the binuclear center (BNC) based on the newly determined crystal structure. The calculations strongly suggest a so-called cis:b(3) mechanism, while the commonly suggested trans mechanism is found to be energetically unfavorable. The mechanism suggested here involves a stable cis-hyponitrite, and it is shown that from this intermediate one N-O bond can be cleaved without the transfer of a proton or an electron into the binuclear active site, in agreement with experimental observations. The fully oxidized intermediate in the catalytic cycle and the resting form of the enzyme are suggested to have an oxo-bridged BNC with two high-spin ferric irons antiferromagnetically coupled. Both steps of reduction of the BNC after N(2)O formation are found to be pH-dependent, also in agreement with experiment. Finally, it is found that the oxo bridge in the oxidized BNC can react with NO to give nitrite, which explains the experimental observations that the fully oxidized enzyme reacts with NO, and most likely also the observed substrate inhibition at higher NO concentrations.

The alkenyl migration mechanism catalyzed by extradiol dioxygenases: a hybrid DFT study.

6-Hydroxymethyl-6-methylcyclohexa-2,4-dienone is a mechanistic probe which when incubated with an extradiol dioxygenase yields a 2-tropolone product. This observation was originally interpreted as evidence supporting a direct heterolytic 1,2-alkenyl migration mechanism for a ring expansion reaction catalyzed by this class of Fe(II)-dependent nonheme enzymes (Xin and Bugg in J Am Chem Soc 130:10422-10430, 2008). In the work reported in this contribution we used quantum chemical methods to test whether such a mechanism is energetically possible and we found that it is not, neither for the mechanistic probe nor for the native catalytic cycle intermediate. Models of increasing complexity were used to calculate energy barriers to the heterolytic 1,2-alkenyl migration and alternative radical mechanisms. It was found that the former involves substantially higher barriers than the latter. A tentative radical mechanism that accounts for the transformation of the probe substrate to 2-tropolone was also proposed, and it involves acceptable barriers.

How is methane formed and oxidized reversibly when catalyzed by Ni-containing methyl-coenzyme M reductase?

Ni-containing methyl-coenzyme M reductase (MCR) is capable of catalyzing methane formation and has recently been observed to also be able to catalyze the reverse reaction, the anaerobic oxidation of methane. The forward reaction has been extensively studied theoretically before and was found to consist of two steps. The first step is rate-limiting and the second step was therefore treated at a lower level. For an accurate treatment of the reverse reaction, both steps have to be studied at the same level. In the present paper, the mechanisms for the reversible formation and oxidation of methane catalyzed by MCR have been investigated using hybrid density functional theory with recent developments, in particular including dispersion effects. An active-site model was constructed based on the X-ray crystal structure. The calculations indicate that the MCR reaction is indeed reversible and proceeds via a methyl radical and a Ni-S(CoM) intermediate with reasonable reaction barriers in both directions. In a competing mechanism, the formation of the crucial Ni-methyl intermediate, was found to be strongly endergonic by over 20 kcal mol(-1) (including a barrier) with dispersion and entropy effects considered, and thus would not be reachable in a reasonable time under natural conditions.