Search

EgtB is a nonheme iron enzyme catalyzing the C-S bond formation between γ-glutamyl cysteine (γGC) and N-α-trimethyl histidine (TMH) in the ergothioneine biosynthesis. Density functional calculations were performed to elucidate and delineate the reaction mechanism of this enzyme. Two different mechanisms were considered, depending on whether the sulfoxidation or the S-C bond formation takes place first. The calculations suggest that the S-O bond formation occurs first between the thiolate and the ferric superoxide, followed by homolytic O-O bond cleavage, very similar to the case of cysteine dioxygenase. Subsequently, proton transfer from a second-shell residue Tyr377 to the newly generated iron-oxo moiety takes place, which is followed by proton transfer from the TMH imidazole to Tyr377, facilitated by two crystallographically observed water molecules. Next, the S-C bond is formed between γGC and TMH, followed by proton transfer from the imidazole CH moiety to Tyr377, which was calculated to be the rate-limiting step for the whole reaction, with a barrier of 17.9 kcal/mol in the quintet state. The calculated barrier for the rate-limiting step agrees quite well with experimental kinetic data. Finally, this proton is transferred back to the imidazole nitrogen to form the product. The alternative thiyl radical attack mechanism has a very high barrier, being 25.8 kcal/mol, ruling out this possibility.

The design and synthesis of biomimetic Mn complexes to catalyze oxygen evolution is a very appealing goal because water oxidation in nature employs a Mn complex. Recently, the mononuclear Mn complex [LMn II (H 2 O) 2 ] 2+ [1, L=Py 2 N(tBu) 2 , Py=pyridyl] was reported to catalyze water oxidation electrochemically at an applied potential of 1.23 V at pH 12.2 in aqueous solution. Density functional calculations were performed to elucidate the mechanism of water oxidation promoted by this catalyst. The calculations showed that 1 can lose two protons and one electron readily to produce [LMn III (OH) 2 ] + (2), which then undergoes two sequential proton-coupled electron-transfer processes to afford [LMn V OO] + (4). The O-O bond formation can occur through direct coupling of the two oxido ligands or through nucleophilic attack of water. These two mechanisms have similar barriers of approximately 17 kcal mol -1 . The further oxidation of 4 to generate [LMn VI OO] 2+ (5), which enables O-O bond formation, has a much higher barrier. In addition, ligand degradation by C-H activation has a similar barrier to that for the O-O bond formation, and this explains the relatively low turnover number of this catalyst., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Water oxidation is a fundamental step in artificial photosynthesis for solar fuels production. In this study, we report a single-site Ru-based water oxidation catalyst, housing a dicarboxylate-benzimidazole ligand, that mediates both chemical and light-driven oxidation of water efficiently under neutral conditions. The importance of the incorporation of the negatively charged ligand framework is manifested in the low redox potentials of the developed complex, which allows water oxidation to be driven by the mild one-electron oxidant [Ru(bpy) 3 ] 3+ (bpy=2,2'-bipyridine). Furthermore, combined experimental and DFT studies provide insight into the mechanistic details of the catalytic cycle., (© 2016 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

The synthesis of two molecular iron complexes, a dinuclear iron(iii,iii) complex and a nonanuclear iron complex, based on the dinucleating ligand 2,2'-(2-hydroxy-5-methyl-1,3-phenylene)bis(1H-benzo[d]imidazole-4-carboxylic acid) is described. The two iron complexes were found to drive the oxidation of water by the one-electron oxidant [Ru(bpy)3](3+).

Nitrogen activation by nitrogenase is one of the most important enzymatic processes on earth. In spite of the determination of X-ray structures of increasingly higher resolution, the nitrogenase mechanism is still not understood. In the most recent X-ray structures it has been shown that a carbon resides in the center of the MoFe-cofactor. Its role is not known. Recent spectroscopic studies, mainly EPR, have come closest to obtaining a molecular mechanism for activating nitrogen. Two hydrides have been shown to play a key role in this context. In the present study, the mechanism for nitrogenase has been investigated by hybrid DFT using a cluster model. This approach has been shown to be very successful for predicting mechanisms for other redox-active enzymes, such as the one for photosystem II, but has so far not been used in its most recent form for nitrogenase. The mechanism obtained has large similarities to the one suggested by spectroscopy, with a reductive elimination of two hydrides just before nitrogen binding. However, a very surprising finding is that the central carbon becomes protonated and has to move out of the cavity as a methyl group before the hydrides can be formed. This has not been suggested before.

PceA is a cobalamin-dependent reductive dehalogenase that catalyzes the dechlorination of perchloroethylene to trichloroethylene and then to cis-dichloroethylene as the sole final product. The reaction mechanism and the regioselectivity of this enzyme are investigated by using density functional calculations. Four different substrates, namely, perchloroethylene, trichloroethylene, cis-dichloroethylene, and chlorotheylene, have been considered and were found to follow the same reaction mechanism pattern. The reaction starts with the reduction of Co(II) to Co(I) through a proton-coupled electron transfer process, with the proton delivered to a Tyr246 anion. This is followed by concerted C-Cl bond heterolytic cleavage and proton transfer from Tyr246 to the substrate carbon atom, generating a Co(III) -Cl intermediate. Subsequently, a one-electron transfer leads to the formation of the Co(II) -Cl product, from which the chloride and the dehalogenated product can be released from the active site. The substrate reactivity follows the trend perchloroethylene>trichloroethylene≫cis-dichloroethylene≫chlorotheylene. The barriers for the latter two substrates are significantly higher compared with those for perchloroethylene and trichloroethylene, implying that PceA does not catalyze their degradation. In addition, the formation of cis-dichloroethylene has a lower barrier by 3.8 kcal mol(-1) than the formation of trans-dichloroethylene and 1,1-dichloroethylene, reproducing the regioselectivity. These results agree quite well with the experimental findings, which show cis-dichloroethylene as the sole product in the PceA-catalyzed dechlorination of perchloethylene and trichloroethylene., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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., (© 2016 Wiley Periodicals, Inc.)

PhoX is a pentanuclear metalloenzyme that employs two ferric ions and three calcium ions to catalyze the hydrolysis of phosphomonoesters. On the basis of the X-ray structure of PhoX ( Science 2014 , 345 , 1170 - 1173 ), a model of the active site is designed, and quantum chemical calculations are used to investigate the reaction mechanism of this enzyme. The calculations support the experimental suggestion, in which the two high spin ferric ions interact in an antiferromagnetic fashion. The two step mechanism proposed by experimentalists has been investigated. The nucleophilic attack of a trinuclear bridging oxo group on the phosphorus center was calculated to be the first step, which is concomitant with the departure of the phenolate, which is stabilized by a calcium ion. The second step is a reverse attack by a water molecule activated by a calcium-bound hydroxide, leading to the regeneration of the bridging oxo group. The second step was calculated to have a barrier of 27.6 kcal/mol. The high barrier suggests that the alternative mechanism involving phosphate release directly from the active site seems to be more likely. All five metal ions are involved in the catalysis by stabilizing the pentacoordinated trigonal bipyramidal transition states.

For the main parts of the mechanism for water oxidation in photosystem II there has recently been very strong experimental support for the mechanism suggested by theoretical model studies. The question addressed in the present study is to what extent this knowledge can be used for the design of artificial catalysts. A major requirement for a useful artificial catalyst is that it is small enough to be synthesized. Small catalysts also have the big advantage that they could improve the catalysis per surface area. To make the mechanism found for PSII useful in this context, it needs to be analyzed in detail. A small model system was therefore used and the ligands were replaced one by one by water-derived ligands. Only the main chemical step of O-O bond formation was investigated in this initial study. The energetics for this small model and the larger one previously used for PSII are remarkably similar, which is the most important result of the present study. This shows that small model complexes have a potential for being very good water oxidation catalysts. It was furthermore found that there is a clear correlation between the barrier height for O-O bond formation and the type of optimal structure for the S3 state. The analysis shows that a flexible central part of the complex is the key for efficient water oxidation., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Density functional theory calculations have been used to study the reaction mechanism of water oxidation catalyzed by a tetranuclear Mn-oxo cluster Mn4O4L6 (L=(C6H4)2PO4(-)). It is proposed that the O-O bond formation mechanism is different in the gas phase and in a water solution. In the gas phase, upon phosphate ligand dissociation triggered by light absorption, the O-O bond formation starting with both the Mn4(III,III,IV,IV) and Mn4(III,IV,IV,IV) oxidation states has to take place via direct coupling of two bridging oxo groups. The calculated barriers are 42.3 and 37.1 kcal/mol, respectively, and there is an endergonicity of more than 10 kcal/mol. Additional photons are needed to overcome these large barriers. In water solution, water binding to the two vacant sites of the Mn ions, again after phosphate dissociation triggered by light absorption, is thermodynamically and kinetically very favorable. The catalytic cycle is suggested to start from the Mn4(III,III,III,IV) oxidation state. The removal of three electrons and three protons leads to the formation of a Mn4(III,IV,IV,IV)-oxyl radical complex. The O-O bond formation then proceeds via a nucleophilic attack of water on the Mn(IV)-oxyl radical assisted by a Mn-bound hydroxide that abstracts a proton during the attack. This step was calculated to be rate-limiting with a total barrier of 29.2 kcal/mol. This is followed by proton-coupled electron transfer, O2 release, and water binding to start the next catalytic cycle., (Copyright © 2014 Elsevier B.V. All rights reserved.)

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., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Insight into how H2 O is oxidized to O2 is envisioned to facilitate the rational design of artificial water oxidation catalysts, which is a vital component in solar-to-fuel conversion schemes. Herein, we report on the mechanistic features associated with a dinuclear Ru-based water oxidation catalyst. The catalytic action of the designed Ru complex was studied by the combined use of high-resolution mass spectrometry, electrochemistry, and quantum chemical calculations. Based on the obtained results, it is suggested that the designed ligand scaffold in Ru complex 1 has a non-innocent behavior, in which metal-ligand cooperation is an important part during the four-electron oxidation of H2 O. This feature is vital for the observed catalytic efficiency and highlights that the preparation of catalysts housing non-innocent molecular frameworks could be a general strategy for accessing efficient catalysts for activation of H2 O., (© 2015 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA. This is an open access article under the terms of Creative Commons Attribution NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.)

The mechanism of water reduction catalysed by a mononuclear copper complex Cu(bztpen) (bztpen = N-benzyl-N,N',N'-tris(pyridine-2-ylmethyl)ethylenediamine) has been elucidated by DFT calculations, revealing that hydrogen evolution proceeds via coupling of a Cu(ii)-hydride and a pendant pyridinium, and providing important implications for the future design of new catalytic systems for water reduction.

The increasing energy demand calls for the development of sustainable energy conversion processes. Here, the splitting of H2O to O2 and H2, or related fuels, constitutes an excellent example of solar-to-fuel conversion schemes. The critical component in such schemes has proven to be the catalyst responsible for mediating the four-electron oxidation of H2O to O2. Herein, we report on the unexpected formation of a single-site Ru complex from a ligand envisioned to accommodate two metal centers. Surprising N-N bond cleavage of the designed dinuclear ligand during metal complexation resulted in a single-site Ru complex carrying a carboxylate-amide motif. This ligand lowered the redox potential of the Ru complex sufficiently to permit H2O oxidation to be carried out by the mild one-electron oxidant [Ru(bpy)3](3+) (bpy = 2,2'-bipyridine). The work thus highlights that strongly electron-donating ligands are important elements in the design of novel, efficient H2O oxidation catalysts.

In a previous detailed study of all the steps of water oxidation in photosystem II, it was surprisingly found that O2 release is as critical for the rate as O-O bond formation. A new mechanism for O2 release has now been found, which can be described as an opening followed by a closing of the interior of the oxygen evolving complex. A transition state for peroxide rotation forming a superoxide radical, missed in the previous study, and a structural change around the outside manganese are two key steps in the new mechanism. However, O2 release may still remain rate-limiting. Additionally, for the step forming the O-O bond, an alternative, experimentally suggested, mechanism was investigated. The new model calculations can rule out the precise use of that mechanism. However, a variant with a rotation of the ligands around the outer manganese by about 30° will give a low barrier, competitive with the old DFT mechanism. Both these mechanisms use an oxyl-oxo mechanism for O-O bond formation involving the same two manganese atoms and the central oxo group (O5).

Electrocatalytic water oxidation using the oxidatively robust 2,7-[bis(2-pyridylmethyl)aminomethyl]-1,8-naphthyridine ligand (BPMAN)-based dinuclear copper(II) complex, [Cu2(BPMAN)(μ-OH)](3+), has been investigated. This catalyst exhibits high reactivity and stability towards water oxidation in neutral aqueous solutions. DFT calculations suggest that the O-O bond formation takes place by an intramolecular direct coupling mechanism rather than by a nucleophilic attack of water on the high-oxidation-state Cu(IV)=O moiety., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Most of the main features of water oxidation in photosystem II are now well understood, including the mechanism for O-O bond formation. For the intermediate S2 and S3 structures there is also nearly complete agreement between quantum chemical modeling and experiments. Given the present high degree of consensus for these structures, it is of high interest to go back to previous suggestions concerning what happens in the S2-S3 transition. Analyses of extended X-ray adsorption fine structure (EXAFS) experiments have indicated relatively large structural changes in this transition, with changes of distances sometimes larger than 0.3 Å and a change of topology. In contrast, our previous density functional theory (DFT)(B3LYP) calculations on a cluster model showed very small changes, less than 0.1 Å. It is here found that the DFT structures are also consistent with the EXAFS spectra for the S2 and S3 states within normal errors of DFT. The analysis suggests that there are severe problems in interpreting EXAFS spectra for these complicated systems.

Herein is described the preparation of a dinuclear molecular Ru catalyst for H2O oxidation. The prepared catalyst mediates the photochemical oxidation of H2O with an efficiency comparable to state-of-the-art catalysts.

The synthesis of Mn-based catalysts to mimic the structural and catalytic properties of the oxygen-evolving complex in photosystem II is a long-standing goal for researchers. An interesting result in this field came with the synthesis of a Mn complex that enables water oxidation driven by the mild single-electron oxidant [Ru(bpy)3](3+). On the basis of hybrid density functional calculations, we herein propose a water oxidation mechanism for this bioinspired Mn catalyst, where the crucial O-O bond formation proceeds from the formal Mn4(IV,IV,IV,V) state by direct coupling of a Mn(IV)-bound terminal oxyl radical and a di-Mn bridging oxo group, a mechanism quite similar to the presently leading suggestion for the natural system. Of importance here is that the designed ligand is shown to be redox-active and can therefore store redox equivalents during the catalytic transitions, thereby alleviating the redox processes at the Mn centers.

Ring hydroxylation and coupled rearrangement reactions catalyzed by 4-hydroxyphenylpyruvate dioxygenase were studied with the QM/MM method ONIOM(B3LYP:AMBER). For electrophilic attack of the ferryl species on the aromatic ring, five channels were considered: attacks on the three ring atoms closest to the oxo ligand (C1, C2, C6) and insertion of oxygen across two bonds formed by them (C1-C2, C1-C6). For the subsequent migration of the carboxymethyl substituent, two possible directions were tested (C1→C2, C1→C6), and two different mechanisms were sought (stepwise radical, single-step heterolytic). In addition, formation of an epoxide (side)product and benzylic hydroxylation, as catalyzed by the closely related hydroxymandelate synthase, were investigated. From the computed reaction free energy profiles it follows that the most likely mechanism of 4-hydroxyphenylpyruvate dioxygenase involves electrophilic attack on the C1 carbon of the ring and subsequent single-step heterolytic migration of the substituent. Computed values of the kinetic isotope effect for this step are inverse, consistent with available experimental data. Electronic structure arguments for the preferred mechanism of attack on the ring are also presented.

The electronic structure of the Mn/Fe cofactor identified in a new class of oxidases (R2lox) described by Andersson and Högbom [Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 5633] is reported. The R2lox protein is homologous to the small subunit of class Ic ribonucleotide reductase (R2c) but has a completely different in vivo function. Using multifrequency EPR and related pulse techniques, it is shown that the cofactor of R2lox represents an antiferromagnetically coupled Mn(III)/Fe(III) dimer linked by a μ-hydroxo/bis-μ-carboxylato bridging network. The Mn(III) ion is coordinated by a single water ligand. The R2lox cofactor is photoactive, converting into a second form (R2loxPhoto) upon visible illumination at cryogenic temperatures (77 K) that completely decays upon warming. This second, unstable form of the cofactor more closely resembles the Mn(III)/Fe(III) cofactor seen in R2c. It is shown that the two forms of the R2lox cofactor differ primarily in terms of the local site geometry and electronic state of the Mn(III) ion, as best evidenced by a reorientation of its unique (55)Mn hyperfine axis. Analysis of the metal hyperfine tensors in combination with density functional theory (DFT) calculations suggests that this change is triggered by deprotonation of the μ-hydroxo bridge. These results have important consequences for the mixed-metal R2c cofactor and the divergent chemistry R2lox and R2c perform.

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.

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., (Copyright © 2014 Elsevier B.V. All rights reserved.)

During recent years significant progress has been made towards the realization of a sustainable and carbon-neutral energy economy. One promising approach is photochemical splitting of H2O into O2 and solar fuels, such as H2. However, the bottleneck in such artificial photosynthetic schemes is the H2O oxidation half reaction where more efficient catalysts are required that lower the kinetic barrier for this process. In particular catalysts based on earth-abundant metals are highly attractive compared to catalysts comprised of noble metals. We have now synthesized a library of dinuclear Mn2(II,III) catalysts for H2O oxidation and studied how the incorporation of different substituents affected the electronics and catalytic efficiency. It was found that the incorporation of a distal carboxyl group into the ligand scaffold resulted in a catalyst with increased catalytic activity, most likely because of the fact that the distal group is able to promote proton-coupled electron transfer (PCET) from the high-valent Mn species, thus facilitating O-O bond formation.

The main parts of the water oxidation mechanism in photosystem II have now been established both from theory and experiments. Still, there are minor questions remaining. One of them concerns the charge and the protonation state of the oxygen evolving complex (OEC). Previously, theory and experiments have agreed that the two water derived ligands on the outer manganese should be one hydroxide and one water. In the present study it is investigated whether both of them could be water. This question is addressed by a detailed study of energy diagrams, but in this context it is more conclusive to compare the redox potential of the OEC to the one of TyrZ. Both procedures lead to the conclusion that one of the ligands is a hydroxide. Another question concerns the protonation of the second shell His337, where the results are more ambiguous. The final part of the present study describes results when calcium is removed from the OEC. Even though protons enter to compensate the charge of the missing Ca(2+), the redox potential and the pKa value of the OEC change dramatically and prevent the progress after S2.

It has experimentally been found that certain mutations close to the entry point of the proton transfer channel in cytochrome c oxidase stop proton translocation but not the oxygen reduction chemistry. This effect is termed uncoupling. Since the mutations are 20Å away from the catalytic center, this is very surprising. A new explanation for this phenomenon is suggested here, involving a local effect at the entry point of the proton channel, rather than the long range effects suggested earlier., (Copyright © 2013 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.)

The full sequence of intermediates in the water oxidation process in photosystem II has recently been characterized by model calculations, in good agreement with experiments. In the present paper, the energy diagram obtained is used as a benchmark test for several density functionals. Only the results using B3LYP with 15% or 20% show good agreement with experiments. The other functionals tried show errors for some energy levels as large as 20-30 kcal/mol. The reason for these large errors is that the error for three consecutive oxidations of Mn(III) to Mn(IV) accumulates as the cluster is oxidized.

Although metallocofactors are ubiquitous in enzyme catalysis, how metal binding specificity arises remains poorly understood, especially in the case of metals with similar primary ligand preferences such as manganese and iron. The biochemical selection of manganese over iron presents a particularly intricate problem because manganese is generally present in cells at a lower concentration than iron, while also having a lower predicted complex stability according to the Irving-Williams series (Mn(II) < Fe(II) < Ni(II) < Co(II) < Cu(II) > Zn(II)). Here we show that a heterodinuclear Mn/Fe cofactor with the same primary protein ligands in both metal sites self-assembles from Mn(II) and Fe(II) in vitro, thus diverging from the Irving-Williams series without requiring auxiliary factors such as metallochaperones. Crystallographic, spectroscopic, and computational data demonstrate that one of the two metal sites preferentially binds Fe(II) over Mn(II) as expected, whereas the other site is nonspecific, binding equal amounts of both metals in the absence of oxygen. Oxygen exposure results in further accumulation of the Mn/Fe cofactor, indicating that cofactor assembly is at least a two-step process governed by both the intrinsic metal specificity of the protein scaffold and additional effects exerted during oxygen binding or activation. We further show that the mixed-metal cofactor catalyzes a two-electron oxidation of the protein scaffold, yielding a tyrosine-valine ether cross-link. Theoretical modeling of the reaction by density functional theory suggests a multistep mechanism including a valyl radical intermediate.

Hybrid DFT model calculations have been performed for some cobalt complexes capable of oxidizing water. Since a very plausible mechanism for the oxygen-evolving complex involving the cuboidal Mn4Ca structure in photosystem II (PSII) has recently been established, the most important part of the present study concerns a detailed comparison between cobalt and manganese as water oxidation catalysts. One similarity found is that a M(IV)-O(•) state is the key precursor for O-O bond formation in both cases. This means that simply getting a M(IV) state is not enough; a formal M(V)═O state is required, with two oxidations on one center from M(III). For cobalt, not even that is enough. A singlet coupled state is required at this oxidation level, which is not the ground state. It is shown that there are also more fundamental differences between catalysts based on these metals. The favorable low-barrier direct coupling mechanism found for PSII is not possible for the corresponding cobalt complexes. The origin of this difference is explained. For the only oxygen-evolving cubic Co4O4 complex with a defined structure, described by Dismukes et al., the calculated results are in good agreement with experiments. For the Co4 models of the amorphous cobalt-oxo catalyst found by Nocera et al., higher barriers are found than the one obtained experimentally. The reasons for this are discussed.

The present status of DFT studies on water oxidation in photosystem II is described. It is argued that a full understanding of all steps is close. In each S-transition, the manganese that is oxidized and the proton released are strongly implicated, and structures of all intermediates have been determined. For the S2-state, recent important experimental findings support key elements of the structure and the mechanism. In this mechanism, the O-O bond is formed between an oxyl radical in the center of the cluster and an Mn-bridging μ-oxo ligand, which was suggested already in 2006. The DFT structure of the oxygen evolving complex, suggested in 2008, is very similar to the recent high-resolution X-ray structure. Some new aspects of the interaction between P680 and the OEC are suggested. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems., (Copyright © 2012 Elsevier B.V. All rights reserved.)

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., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Detailed mechanisms for substrate water exchange in the oxygen evolving complex in photosystem II have been determined with DFT methods for large models. Existing interpretations of the experimental water exchange results have been quite different. By many groups, these results have been the main argument against the water oxidation mechanism suggested by DFT, in which the oxygen molecule is formed between a bridging oxo and an oxyl radical ligand in the center of the OEC. That mechanism is otherwise in line with most experiments. The problem has been that the mechanism requires a rather fast exchange of a bridging oxo ligand, which is not a common finding for smaller Mn-containing model systems. However, other groups have actually favored a substrate derived oxo ligand partly based on the same experiments. In the present study, three S-states have been studied, and the rates have been well reproduced by the calculations. The surprising experimental finding that water exchange in S1 is slower than the one in S2 is reproduced and explained. The key to this rate difference is the ease by which one of the manganese centers (Mn3) is reduced. This reduction has to occur to release the substrate water from Mn3. The similar rate of the slow exchange in S2 and S3 has been rationalized on the basis of earlier experiments combined with the present calculations. The results strongly support the previous DFT-suggested water oxidation mechanism.

Activation of manganese-dependent class Ib ribonucleotide reductase by hydrogen peroxide was modeled using B3LYP* hybrid density functional theory. Class Ib ribonucleotide reductase R2 subunit (R2F) does not react with molecular oxygen. Instead R2F is proposed to react with H2O2 or HO2(-), provided by the unusual flavodoxin protein NrdI, to generate the observed manganese(III) manganese(III) tyrosyl-radical state. On the basis of the calculations, an energetically feasible reaction mechanism is suggested for activation by H2O2, which proceeds through two reductive half-reactions. In the first reductive half-reaction, H2O2 is cleaved with a barrier of 13.1 kcal mol(-1) [Mn(II)Mn(II) → Mn(III)Mn(III)], and in the second reductive half-reaction, H2O2 is cleaved with a barrier of 17.0 kcal mol(-1) [Mn(III)Mn(III) → Mn(IV)Mn(IV)]. Tyrosyl-radical formation from both the Mn(IV)Mn(IV) state and a Mn(III)Mn(IV) state, where an electron and proton have been taken up, is both kinetically and thermodynamically accessible. Hence, chemically, H2O2 is a possible oxidant for the manganese-dependent R2F. The selectivity between the second reductive half-reaction and a competing oxidative reaction, as in manganese catalase, may be the time scale for the availability of H2O2. The role of NrdI may be to provide H2O2 on the correct time scale.

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., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Hydroxymandelate synthase (HMS) and 4-hydroxyphenylpyruvate dioxygenase (HPPD) are highly related enzymes using the same substrates but catalyzing hydroxylation reactions yielding different products. The first steps of the HMS and HPPD catalytic reactions are believed to proceed in the same way and lead to an Fe(IV)═O-hydroxyphenylacetate (HPA) intermediate. Further down the catalytic cycles, HMS uses Fe(IV)═O to perform hydroxylation of the benzylic carbon, whereas in HPPD, the reactive oxoferryl intermediate attacks the aromatic ring of HPA. This study focuses on this part of the HMS catalytic cycle that starts from the oxoferryl intermediate and aims to identify interactions within the active site that are responsible for enzyme specificity. To this end, a HMS-Fe(IV)═O-HPA complex was modeled with molecular dynamics simulations. On the basis of the molecular dynamics-equilibrated structure, an active site model suitable for quantum chemical investigations was constructed and used for density functional theory (B3LYP) calculations of the mechanism of the native reaction of HMS, i.e., benzylic hydroxylation, and the alternative electrophilic attack on the ring, which is a step of the HPPD catalytic cycle. The most important result of this study is the finding that the conformation of the Ser201 side chain in the second coordination shell has a key role in directing the reaction of Fe(IV)═O into either the HMS or the HPPD channel.

The Mn 1s near-edge absorption fine structure (NEXAFS) has been computed by means of transition-state gradient-corrected density functional theory (DFT) on four Mn(4)Ca clusters modeling the successive S(0) to S(3) steps of the oxygen-evolving complex (OEC) in photosystem II (PSII). The model clusters were obtained from a previous theoretical study where they were determined by energy minimization. They are composed of Mn(III) and Mn(IV) atoms, progressing from Mn(III)(3)Mn(IV) for S(0) to Mn(III)(2)Mn(IV)(2) for S(1) to Mn(III)Mn(IV)(3) for S(2) to Mn(IV)(4) for S(3), implying an Mn-centered oxidation during each step of the photosynthetic oxygen evolution. The DFT simulations of the Mn 1s absorption edge reproduce the experimentally measured curves quite well. By the half-height method, the theoretical IPEs are shifted by 0.93 eV for the S(0) → S(1) transition, by 1.43 eV for the S(1) → S(2) transition, and by 0.63 eV for the S(2) → S(3) transition. The inflection point energy (IPE) shifts depend strongly on the method used to determine them, and the most interesting result is that the present clusters reproduce the shift in the S(2) → S(3) transition obtained by both the half-height and second-derivative methods, thus giving strong support to the previously suggested structures and assignments.

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.

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.

Quantum chemical calculations of active-site models of nitrous oxide reductase (N(2)OR) have been undertaken to elucidate the mechanism of N-O bond cleavage mediated by the supported tetranuclear Cu(4)S core (Cu(Z)) found in the enzymatic active site. Using either a minimal model previously employed by Gorelsky et al. (J. Am. Chem. Soc. 128:278-290, 2006) or a more extended model including key residue side chains in the active-site second shell, we found two distinct mechanisms. In the first model, N(2)O binds to the fully reduced Cu(Z) in a bent μ-(1,3)-O,N bridging fashion between the Cu(I) and Cu(IV) centers and subsequently extrudes N(2) while generating the corresponding bridged μ-oxo species. In the second model, substrate N(2)O binds loosely to one of the coppers of Cu(Z) in a terminal fashion, i.e., using only the oxygen atom; loss of N(2) generates the same μ-oxo copper core. The free energies of activation predicted for these two alternative pathways are sufficiently close to one another that theory does not provide decisive support for one over the other, posing an interesting problem with respect to experiments that might be designed to distinguish between the two. Effects of nearby residues and active-site water molecules are also explored.

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., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The new high-resolution X-ray structure of photosystem II has allowed more detailed studies than before of water oxidation at the oxygen evolving complex (OEC). In the present study the two final S-transitions of water oxidation are studied. The electron coupled proton transfers are followed from the center of the OEC to Asp61, which is considered as the start of the transfer chain through the protein to the lumenal side. It is found that the proton transfers occur in multiple steps. Structures of intermediates and energy diagrams are derived and compared to experimental observations. Since the new experimental structure of the OEC is very similar to the one suggested earlier by density functional calculations, the O-O bond formation step remains essentially the same as the one suggested five years ago. An interesting new result is that the barrier for proton transfer within the OEC actually competes with the O-O bond formation step of being rate-limiting.

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., (Copyright © 2011 Elsevier B.V. All rights reserved.)

Two-electron chemistry with an iron dimer, a manganese dimer, and a manganese-iron dimer as a catalyst has been modeled using B3LYP* hybrid density functional theory. The recently discovered MnFe proteins form (at least) two functionally distinct groups, performing radical generation (class Ic ribonucleotide reductase subunit II) and substrate oxidations (subunit II-like ligand-binding oxidases, R2lox), respectively. Proteins from the latter group appear to be functionally similar to the diiron carboxylate proteins that perform two-electron oxidations of substrates, such as methane monooxygenase. To qualitatively determine the potential role of a MnFe center in R2lox, methane hydroxylation with the MnFe heterodimer and with the FeFe and MnMn homodimers is studied. The redox potential of the active state of the Mn(IV)Fe(IV) heterodimer is about 7 kcal mol(-1) lower than that of the active state of the Fe(IV)Fe(IV) homodimer, leading to a high barrier for the rate-limiting hydrogen abstraction with the MnFe site. If the entropy loss is not included, the barriers are lower, and the MnFe heterodimer can therefore have a role in R2lox as an oxidase for larger substrates exergonically bound to the protein. A MnMn center has a high barrier both with and without entropy loss. The higher stability of Fe(IV) in the presence of Mn(IV) in the other site compared with a second Fe(IV) suggests an explanation for the presence of the MnFe site in R2lox: to provide a metal center that is capable of two-electron chemistry, and which is more stable and less sensitive to external reductants than an Fe(IV)Fe(IV) site.

The procedure for fixing atoms of amino acid residues in cluster model calculations on enzymes is reviewed. Examples from recent calculations on photosystem II (PSII) and Mo,Cu-dependent CO dehydrogenase are given. In this context, the cluster model work on finding a mechanism for O-O bond formation and a structure of the oxygen-evolving complex in PSII is also reviewed. In that work, fixing certain atoms played an important role. The main part of the present study concerns the mechanism in PSII using models based on the new high-resolution (1.9 Å) X-ray structure, which is compared to that using the old, theoretically suggested, structure. It is concluded that the mechanism remains the same, with a similar barrier height. Finally, a connection between the OEC structure and Mn,Ca-containing minerals is also briefly discussed., (Copyright © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Two artificial water oxidation catalysts, the blue dimer and the Llobet catalyst, have been studied using hybrid DFT methods. The results are compared to those for water oxidation in the natural photosystem II enzyme. Studies on the latter system have now reached a high level of understanding, at present much higher than the one for the artificial systems. A recent high resolution X-ray structural investigation of PSII has confirmed the main features of the structure of the oxygen evolving complex (OEC) suggested by previous DFT cluster studies. The O-O bond formation mechanism suggested is of direct coupling (DC) type between an oxygen radical and a bridging oxo ligand. A similar DC mechanism is found for the Llobet catalyst, while an acid-base (AB) mechanism is preferred for the blue dimer. All of them require at least one oxygen radical. Full energy diagrams, including both redox and chemical steps, have been constructed illustrating similarities and differences to the natural system. Unlike previous DFT studies, the results of the present study suggest that the blue dimer is rate-limited by the initial redox steps, and the Llobet catalyst by O(2) release. The results could be useful for further improvement of the artificial systems.

In the present mini-review, computational work over the past decade on water oxidation in photosystem II (PSII) is summarized. The size of the chemical model used for the oxygen evolving complex (OEC) has during this time increased from the initial 20 atoms to the present day 220 atoms. The electronic structure methods used have during the same period only undergone minor improvements. It is concluded that the results have now reached a high level of convergence and the predictions for both the structure of the OEC and the O-O bond formation mechanism are most probably of higher accuracy than presently available from experiments., (Copyright © 2011 Elsevier B.V. All rights reserved.)

A methyl-Co(cobalamin) species has been characterized to be a crucial intermediate in the last step of the de novo biosynthesis of methionine catalyzed by cobalamin-dependent methionine synthase (MetH). However, exactly how it is formed is still an open question. In the present article, the formation of the methyl-Co(cobalamin) species in MetH has been investigated with B3LYP* hybrid DFT including van der Waals (vdW) interactions (i.e., dispersion) and using a chemical model built on X-ray crystal structures. The methyl cation and radical transfer mechanisms have been examined in various protonation states. The calculations reveal that the CH(3)-Co(III)(cobalamin) formation in MetH proceeds along a stepwise pathway, where the first step is a methyl cation transfer from the protonated methyltetrahydrofolate (CH(3)-THF) substrate to the Co(I)cobalamin. The second step is a binding of His759 to the other side (α-face) of Co. The former methyl transfer is computed to be the rate-limiting step with a barrier of 18 kcal/mol, which is reduced to 13 kcal/mol when dispersion is included. For the first step, the protonation at the methyl-bound nitrogen of CH(3)-THF is very important. The methyl transfer is otherwise unreachable with a very high barrier of ~38 kcal/mol. The deprotonation of the α-face His759-Asp757-Ser810 triad is found to be much less significant but slightly facilitates the CH(3)-Co(III)Cbl formation. There has been a long-standing discrepancy of 10-20 kcal/mol between theory and experiment in previous B3LYP computations of the Co-C bond dissociation energy for the methyl-Co(cobalamin) species. The calculations indicate that the lack of dispersion (~11 kcal/mol) is the main origin of this puzzling problem. With these effects, B3LYP* gives a bond strength of 32 kcal/mol compared to the experimental value of 37 ± 3 kcal/mol. Overall, the present calculations give many examples of dispersion that makes non-negligible contributions to the energetics of enzyme reactions, especially for systems involving at least one large reacting fragment approaching or departing.

The oxygen cleavage in Chlamydia trachomatis ribonucleotide reductase (RNR) has been studied using B3LYP* hybrid density functional theory. Class Ic C. trachomatis RNR lacks the radical-bearing tyrosine, crucial for activity in conventional class I (subclass a and b) RNR. Instead of the Fe(III)Fe(III)-Tyr(rad) active state, C. trachomatis RNR has a mixed Mn(IV)Fe(III) metal center in subunit II (R2). A mixed MnFe metal center has never been observed as a radical cofactor before. The active state is generated by reductive oxygen cleavage at the metal site. On the basis of calculated barriers for oxygen cleavage in C. trachomatis R2 and R2 from Escherichia coli with a diiron, a mixed manganese-iron, and a dimanganese center, conclusions can be drawn about the effect of changing metals in R2. The oxygen cleavage is found to be governed by two factors: the redox potentials of the metals and the relative stability of the different peroxides. Mn(IV) has higher stability than Fe(IV), and the barrier is therefore lower with a mixed metal center than with a diiron center. With a dimanganese center, an asymmetric peroxide is more stable than the symmetric peroxide, and the barrier therefore becomes too high. Calculated proton-coupled redox potentials are compared to identify three possible R2 active states, the Fe(III)Fe(III)-Tyr(rad) state, the Mn(IV)Fe(III) state, and the Mn(IV)Mn(IV) state. A tentative energy profile of the thermodynamics of the radical transfer from R2 to subunit I is constructed to illustrate how the stability of the active states can be understood from a thermodynamical point of view.

Adipic acid is a key compound in the chemical industry, where it is mainly used in the production of polymers. The conventional process of its generation requires vast amounts of energy and, moreover, produces environmentally deleterious substances. Thus, there is interest in alternative ways to gain adequate amounts of adipic acid. Experimental reports on a one-pot iron-catalyzed conversion of cyclohexane to adipic acid motivated a theoretical investigation based on density functional theory calculations. The process investigated is interesting because it requires less energy than contemporary methods and does not produce environmentally harmful side products. The aim of the present contribution is to gain insight into the mechanism of the iron-catalyzed cyclohexane conversion to provide a basis for the further development of this process. The rate-limiting step of the process is discussed, but considering the accuracy of the calculations, it is difficult to ensure whether the rate-limiting step is in the substrate oxidation or in the generation of the catalytically active species. It is shown that the slowest step in the substrate oxidation is the conversion of cyclohexanol to cyclohexane-1,2-diol. Hydrogen-atom transfer from one of the OH groups of cyclohexane-1,2-diol makes the intradiol cleavage occur spontaneously.