Regioselectivity of aliphatic versus aromatic hydroxylation by a nonheme iron(II)-superoxo complex

Many enzymes in nature utilize molecular oxygen on an iron center for the catalysis of substrate hydroxylation. In recent years, great progress has been made in understanding the function and properties of iron(IV)-oxo complexes; however, little is known about the reactivity of iron(II)-superoxo intermediates in substrate activation. It has been proposed recently that iron(II)-superoxo intermediates take part as hydrogen abstraction species in the catalytic cycles of nonheme iron enzymes. To gain insight into oxygen atom transfer reactions by the nonheme iron(II)-superoxo species, we performed a density functional theory study on the aliphatic and aromatic hydroxylation reactions using a biomimetic model complex. The calculations show that nonheme iron(II)-superoxo complexes can be considered as effective oxidants in hydrogen atom abstraction reactions, for which we find a low barrier of 14.7 kcal mol(-1) on the sextet spin state surface. On the other hand, electrophilic reactions, such as aromatic hydroxylation, encounter much higher (>20 kcal mol(-1)) barrier heights and therefore are unlikely to proceed. A thermodynamic analysis puts our barrier heights into a larger context of previous studies using nonheme iron(IV)-oxo oxidants and predicts the activity of enzymatic iron(II)-superoxo intermediates.     

Computational exploration of the catalytic mechanism of dopamine beta-monooxygenase: modeling of its mononuclear copper active sites

Dopamine hydroxylation by the copper-superoxo, -hydroperoxo, and -oxo species of dopamine beta-monooxygenase (DBM) is investigated using theoretical calculations to identify the active species in its reaction and to reveal the key functions of the surrounding amino acid residues in substrate binding. A 3D model of rat DBM is constructed by homology modeling using the crystal structure of peptidylglycine alpha-hydroxylating monooxygenase (PHM) with a high sequence identity of 30% as a template. In the constructed 3D model, the CuA site in domain 1 is coordinated by three histidine residues, His265, His266, and His336, while the CuB site in domain 2 is coordinated by two histidine residues, His415 and His417, and by a methionine residue, Met490. The three Glu268, Glu369, and Tyr494 residues are suggested to play an important role in the substrate binding at the active site of DBM to enable the stereospecific hydrogen-atom abstraction. Quantum mechanical/molecular mechanical (QM/MM) calculations are performed to determine the structure of the copper-superoxo, -hydroperoxo, and -oxo species in the whole-enzyme model with about 4700 atoms. The reactivity of the three oxidants is evaluated in terms of density-functional-theory calculations with small models extracted from the QM region of the whole-enzyme model.     

Oxygen activation by the noncoupled binuclear copper site in peptidylglycine alpha-hydroxylating monooxygenase. Reaction mechanism and role of the noncoupled nature of the active site

Reaction thermodynamics and potential energy surfaces are calculated using density functional methods to investigate possible reactive Cu/O(2) species for H-atom abstraction in peptidylglycine alpha-hydroxylating monooxygenase (PHM), which has a noncoupled binuclear Cu active site. Two possible mononuclear Cu/O(2) species have been evaluated, the 2-electron reduced Cu(II)(M)-OOH intermediate and the 1-electron reduced side-on Cu(II)(M)-superoxo intermediate, which could form with comparable thermodynamics at the catalytic Cu(M) site. The substrate H-atom abstraction reaction by the Cu(II)(M)-OOH intermediate is found to be thermodynamically accessible due to the contribution of the methionine ligand, but with a high activation barrier ( approximately 37 kcal/mol, at a 3.0-A active site/substrate distance), arguing against the Cu(II)(M)-OOH species as the reactive Cu/O(2) intermediate in PHM. In contrast, H-atom abstraction from substrate by the side-on Cu(II)(M)-superoxo intermediate is a nearly isoenergetic process with a low reaction barrier at a comparable active site/substrate distance ( approximately 14 kcal/mol), suggesting that side-on Cu(II)(M)-superoxo is the reactive species in PHM. The differential reactivities of the Cu(II)(M)-OOH and Cu(II)(M)-superoxo species correlate to their different frontier molecular orbitals involved in the H-atom abstraction reaction. After the H-atom abstraction, a reasonable pathway for substrate hydroxylation involves a "water-assisted" direct OH transfer to the substrate radical, which generates a high-energy Cu(II)(M)-oxyl species. This provides the necessary driving force for intramolecular electron transfer from the Cu(H) site to complete the reaction in PHM. The differential reactivity pattern between the Cu(II)(M)-OOH and Cu(II)(M)-superoxo intermediates provides insight into the role of the noncoupled nature of PHM and dopamine beta-monooxygenase active sites, as compared to the coupled binuclear Cu active sites in hemocyanin, tyrosinase, and catechol oxidase, in O(2) activation.     

Conversion of methane to methanol at the mononuclear and dinuclear copper sites of particulate methane monooxygenase (pMMO): a DFT and QM/MM study

Methane hydroxylation at the mononuclear and dinuclear copper sites of pMMO is discussed using quantum mechanical and QM/MM calculations. Possible mechanisms are proposed with respect to the formation of reactive copper-oxo and how they activate methane. Dioxygen is incorporated into the Cu(I) species to give a Cu(II)-superoxo species, followed by an H-atom transfer from a tyrosine residue near the monocopper active site. A resultant Cu(II)-hydroperoxo species is next transformed into a Cu(III)-oxo species and a water molecule by the abstraction of an H-atom from another tyrosine residue. This process is accessible in energy under physiological conditions. Dioxygen is also incorporated into the dicopper site to form a (mu-eta(2):eta(2)-peroxo)dicopper species, which is then transformed into a bis(mu-oxo)dicopper species. The formation of this species is more favorable in energy than that of the monocopper-oxo species. The reactivity of the Cu(III)-oxo species is sufficient for the conversion of methane to methanol if it is formed in the protein environment. Since the sigma orbital localized in the Cu-O bond region is singly occupied in the triplet state, this orbital plays a role in the homolytic cleavage of a C-H bond of methane. The reactivity of the bis(mu-oxo)dicopper species is also sufficient for the conversion of methane to methanol. The mixed-valent bis(mu-oxo)Cu(II)Cu(III) species is reactive to methane because the amplitude of the sigma singly occupied MO localized on the bridging oxo moieties plays an essential role in C-H activation.     

Characterization of the structure and reactivity of monocopper-oxygen complexes supported by beta-diketiminate and anilido-imine ligands

Copper-oxygen complexes supported by beta-diketiminate and anilido-imine ligands have recently been reported (Aboelella et al., J Am Chem Soc 2004, 126, 16896; Reynolds et al., Inorg Chem 2005, 44, 6989) as potential biomimetic models for dopamine beta-monooxygenase (DbetaM) and peptidylglycine alpha-hydroxylating monooxygenase (PHM). However, in contrast to the enzymatic systems, these complexes fail to exhibit C--H hydroxylation activity (Reynolds et al., Chem Commun 2005, 2014). Quantum chemical characterization of the 1:1 Cu-O(2) model adducts and related species (Cu(III)-hydroperoxide, Cu(III)-oxo, and Cu(III)-hydroxide) indicates that the 1:1 Cu-O(2) adducts are unreactive toward substrates because of the weakness of the O--H bond that would be formed upon hydrogen-atom abstraction. This in turn is ascribed to the 1:1 adducts having both low reduction potentials and basicities. Cu(III)-oxo species on the other hand, determined to be intermediate between Cu(III)-oxo and Cu(II)-oxyl in character, are shown to be far more reactive toward substrates. Based on these results, design strategies for new DbetaM and PHM biomimetic ligands are proposed: new ligands should be made less electron rich so as to favor end-on dioxygen coordination in the 1:1 Cu-O(2) adducts. Comparison of the relative reactivities of the various copper-oxygen complexes as hydroxylating agents provides support for a Cu(II)-superoxide species as the intermediate responsible for substrate hydroxylation in DbetaM and PHM, and suggests that a Cu(III)-oxo intermediate would be competent in this process as well.     

A DFT study of the mechanism of Ni superoxide dismutase (NiSOD): role of the active site cysteine-6 residue in the oxidative half-reaction

In the present DFT study, the catalytic mechanism of H2O2 formation in the oxidative half-reaction of NiSOD, E-Ni(II) + O2- + 2H+ --> E-Ni(III) + H2O2, has been investigated. The main objective of this study is to investigate the source of two protons required in this half-reaction. The proposed mechanism consists of two steps: superoxide coordination and H2O2 formation. The effect of protonation of Cys6 and the proton donating roles of side chains (S) and backbones (B) of His1, Asp3, Cys6, and Tyr9 residues in these two steps have been studied in detail. For protonated Cys6, superoxide binding generates a Ni(III)-O2H species in a process that is exothermic by 17.4 kcal/mol (in protein environment using the continuum model). From the Ni(III)-O2H species, H2O2 formation occurs through a proton donation by His1 via Tyr9, which relative to the resting position of the enzyme is exothermic by 4.9 kcal/mol. In this pathway, a proton donating role of His1 residue is proposed. However, for unprotonated Cys6, a Ni(II)-O2- species is generated in a process that is exothermic by 11.3 kcal/mol. From the Ni(II)-O2- species, the only feasible pathway for H2O2 formation is through donation of protons by the Tyr9(S)-Asp3(S) pair. The results discussed in this study elucidate the role of the active site residues in the catalytic cycle and provide intricate details of the complex functioning of this enzyme.     

What is the active species of cytochrome P450 during camphor hydroxylation? QM/MM studies of different electronic states of compound I and of reduced and oxidized iron-oxo intermediates

We have investigated C-H hydroxylation of camphor by Compound I (Cpd I) of cytochrome P450cam in different electronic states and by its one-electron reduced and oxidized forms, using QM/MM calculations in the native protein/solvent environment. Cpd I species with five unpaired electrons (pentaradicaloids) are ca. 12 kcal/mol higher in energy than the ground state Cpd I species with three unpaired electrons (triradicaloids). The H-abstraction transition states of pentaradicaloids lie ca. 21 (9) kcal/mol above the triradicaloid (pentaradicaloid) reactants. Hydroxylation via pentaradicaloids is thus facile provided that they can react before relaxing to the ground-state triradicaloids. Excited states of Cpd I with an Fe(V)-oxo moiety lie more than 20 kcal/mol above the triradicaloid ground state in single-point gas-phase calculations, but these electronic configurations are not stable upon including the point-charge protein environment which causes SCF convergence to the triradicaloid ground state. One-electron reduced species (Cpd II) show sluggish reactivity compared with Cpd I in agreement with experimental model studies. One-electron oxidized species are more reactive than Cpd I but seem too high in energy to be accessible. The barriers to hydrogen abstraction for the various forms of Cpd I are generally not affected much by the chosen protonation states of the Asp297 and His355 residues near the propionate side chains of the heme or by the appearance of radical character at Asp297, His355, or the propionates.     

Structure,stability, and electronic and NMR properties of various oxo- and nitrido-derivatives of [L(Salen)Mn(III)]+, where L = none and imidazole. A density functional study

Structure, stability, and electronic and NMR properties of [(Salen)Mn(III)](+)-derived intermediates/reactants in the epoxidation/amination of unfunctionalized olefins, namely [(Salen)Mn(V)O](+) (1-oxo), [(Salen)Mn(IV)O] (2-oxo), and [(Salen)Mn(V)N] (3), have been studied with the B3LYP density functional method. It has been shown that the (1)A, (3)A, and (5)A states of cationic 1-oxo species are virtually degenerate, while for the neutral 2-oxo species the ground (4)A state lies 6.4 kcal/mol lower than (2)A. In the nitrido species 3, the (1)A state has been shown to be the ground state in agreement with experiment. We have investigated isomerization of 1-oxo and 2-oxo species into unusual [(OSalen)Mn(III)](+) (1-N-oxo and 1-peroxo) and [(OSalen)Mn(II)] (2-N-oxo and 2-peroxo) species, respectively. For cationic species 1, the 1-N-oxo isomers are more stable (by 10-12 kcal/mol) than the 1-oxo isomer and are separated from the latter by 21-22 kcal/mol barriers. On the other hand, 1-peroxo isomers are calculated to be 14-16 kcal/mol higher than the 1-oxo isomer. For neutral species 2, however, both 2-N-oxo and 2-peroxo isomers lie significantly higher in energy than the 2-oxo isomer. It has been shown that coordination of axial imidazole ligand alters relative energies of spin states for 1- and 2-oxo species, destabilizing low-spin states. For singlet states of H(2)Salen, 1-oxo, and 3, we have calculated (1)H, (13)C, (15)N, and (17)O NMR chemical shifts using the gauge-independent-atomic orbital (GIAO) approach.     

Computational mutation analysis of hydrogen abstraction and radical rearrangement steps in the catalysis of coenzyme B12-dependent diol dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.     

Nonheme iron-oxo and -superoxo reactivities: O2 binding and spin inversion probability matter

DFT calculated barriers for C-H activation of 1,4-cyclohexadiene by nonheme iron(IV)-oxo and iron(III)-superoxo species show that the experimental trends can be explained if the spin inversion probability of the TMC iron(IV)-oxo is assumed to be poor. Also, the TMC iron(III)-superoxo reaction proceeds with an endothermic O(2)-binding energy followed by an intrinsically reactive quintet state.     

A chromium(III)-superoxo complex in oxygen atom transfer reactions as a chemical model of cysteine dioxygenase

Metal-superoxo species are believed to play key roles in oxygenation reactions by metalloenzymes. One example is cysteine dioxygenase (CDO) that catalyzes the oxidation of cysteine with O(2), and an iron(III)-superoxo species is proposed as an intermediate that effects the sulfoxidation reaction. We now report the first biomimetic example showing that a chromium(III)-superoxo complex bearing a macrocyclic TMC ligand, [Cr(III)(O(2))(TMC)(Cl)](+), is an active oxidant in oxygen atom transfer (OAT) reactions, such as the oxidation of phosphine and sulfides. The electrophilic character of the Cr(III)-superoxo complex is demonstrated unambiguously in the sulfoxidation of para-substituted thioanisoles. A Cr(IV)-oxo complex, [Cr(IV)(O)(TMC)(Cl)](+), formed in the OAT reactions by the chromium(III)-superoxo complex, is characterized by X-ray crystallography and various spectroscopic methods. The present results support the proposed oxidant and mechanism in CDO, such as an iron(III)-superoxo species is an active oxidant that attacks the sulfur atom of the cysteine ligand by the terminal oxygen atom of the superoxo group, followed by the formation of a sulfoxide and an iron(IV)-oxo species via an O-O bond cleavage.     

Catalytic mechanism of matrix metalloproteinases: two-layered ONIOM study

The two-layered ONIOM(B3LYP:MNDO) method has been used to investigate the hydrolytical mechanism of matrix metalloproteinases (MMPs), a large family of zinc-dependent endopeptidases capable of degrading a wide range of macromolecules of the extracellular matrix. Human stromelysin-1 (MMP-3) was chosen as a physiologically important member of the MMP family. As a structural reference, X-ray data on the stromelysin-1 catalytic domain (SCD) complexed to the transition state analogue diphenyl piperidine sulfonamide inhibitor was used. The backbone spacer of 11 residues (201-211) was included in the final model, spanning the catalytic Glu202 residue and the three structural His201,205,211 zinc ligands. The polypeptide framework incorporated, partly accounting for the protein rigidity, reduces the activation free energy slightly by 1.6 kcal/mol. Essentially a single-step catalytic mechanism was obtained, generally following a classical proposal for MMPs. Glu202 here acts as a base, abstracting a proton from the metal-bound reactant water and delivering this proton to the peptide nitrogen. An auxiliary water molecule is suggested to be of crucial importance acting as an electrophilic agent to the carbonyl oxygen of the substrate. The direct inclusion of the auxiliary water molecule decreases the activation free energy by about 5 kcal/mol via donation of a strong hydrogen bond. The calculated activation barrier of 13.1 kcal/mol agrees well with experimental rates.     

QM/MM study of mechanisms for compound I formation in the catalytic cycle of cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.     

Dioxygen Oxidation Cu(II) → Cu(III) in the Copper Complex of cyclo(Lys-dHis-βAla-His): A Case Study by EXAFS and XANES Approach

A former spectroscopic study of Cu(II) coordination by the 13-membered ring cyclic tetrapeptide c(Lys-dHis-βAla-His) (DK13), revealed the presence, at alkaline pH, of a stable peptide/Cu(III) complex formed in solution by atmospheric dioxygen oxidation. To understand the nature of this coordination compound and to investigate the role of the His residues in the Cu(III) species formation, Cu K-edge XANES, and EXAFS spectra have been collected for DK13 and two other 13-membered cyclo-peptides: the diastereoisomer c(Lys-His-βAla-His) (LK13), and c(Gly-βAla-Gly-Lys) (GK13), devoid of His residues. Comparison of pre-edge peak features with those of Cu model compounds, allowed us to get information on copper oxidation state in two of the three peptides, DK13 and GK13: DK13 contains only Cu(III) ions in the experimental conditions, while GK13 binds only with Cu(II). For LK13/Cu complex, EXAFS spectrum suggested and UV-vis analysis confirmed the presence of a mixture of Cu(II) and Cu(III) coordinated species. Theoretical XANES spectra have been calculated by means of the MXAN code. The good agreement between theoretical and experimental XANES data collected for DK13, suggests that the refined structure, at least in the first coordination shell around Cu, is a good approximation of the DK13/Cu(III) coordination species present at strongly alkaline pH. All the data are consistent with a slightly distorted pyramidal CuN(4) unit, coming from the peptide bonds. Surprisingly, the His side-chains seemed not involved in the final, stable, Cu(III) scaffold.     

Computational study on mechanistic details of the aminoethanol rearrangement catalyzed by the vitamin B12-dependent ethanolamine ammonia lyase: His and Asp/Glu acting simultaneously as catalytic auxiliaries

The rearrangement of aminoethanol catalyzed by ethanolamine ammonia lyase is investigated by computational means employing DFT (B3LYP/6-31G) and ab initio molecular orbital theory (QCISD/cc-pVDZ). The study aims at providing a detailed account on various crucial aspects, in particular a distinction between a direct intramolecular migration of the partially protonated NH(2) group vs elimination of NH(4)(+). Three mechanistic scenarios were explored: (i) According to the calculations, irrespective of the nature of the protonating species, intramolecular migration of the NH(3) group is energetically less demanding than elimination of NH(4)(+). However, all computed activation enthalpies exceed the experimentally derived activation enthalpy (15 kcal/mol) associated with the rate-determining step, i.e., the hydrogen abstraction from the 5'-deoxyadenosine by the product radical. For example, when imidazole is used as a model system for His interacting with the NH(3) group of the substrate, the activation enthalpy for the migration process amounts to 27.4 kcal/mol. If acetic acid is employed to mimic Asp or Glu, the activation enthalpy is somewhat lower, being equal to 24.2 kcal/mol. (ii) For a partial deprotonation of the substrate 2 at the OH group, the rearrangement mechanism consists of the dissociation of an NH(2) radical from C(2) and its association at C(1) atom. For all investigated proton acceptors (i.e., OH(-), HCOO(-), CH(3)COO(-), CH(2)NH, imidazole), the activation enthalpy for the dissociation step also exceeds 15 kcal/mol. Typical data are 20.2 kcal/mol for Ac(-) and 23.8 kcal/mol for imidazole. (iii) However, in a synergistic action of partial protonation of the NH(2) group and partial deprotonation of the OH group by the two conceivable catalytic auxiliaries Asp/Glu and His, the activation enthalpy computed is compatible with the experimental data. For imidazole and acetate as model systems, the activation enthalpy is equal to 13.7 kcal/mol. This synergistic action of the two catalytic groups is expected to take place in a physiologically realistic pH range of 6-9.5, and the present computational findings may help to further characterize the yet unknown structural details of the ethanolamine ammonia lyase's active site.     

How heme metabolism occurs in heme oxygenase: computational study of oxygen-donation ability of the oxo and hydroperoxo species

We report a density functional theory study on the heme metabolism in heme oxygenase using iron-hydroperoxo and -oxo models. The activation energies for heme oxidation at the alpha-carbon by the iron-hydroperoxo and -oxo species are calculated to be 42.9 and 39.9 kcal/mol, respectively. These high activation barriers lead us to reconsider the catalytic mechanism of heme oxygenase     

Electron paramagnetic resonance detection of intermediates in the enzymatic cycle of an extradiol dioxygenase

Extradiol catecholic dioxygenases catalyze the cleavage of the aromatic ring of the substrate with incorporation of both oxygen atoms from O2. These enzymes are important in nature for the recovery of large amounts of carbon from aromatic compounds. The catalytic site contains either Fe or Mn coordinated by a facial triad of two His and one Glu or Asp residues. Previous studies have shown that Fe(II) and Mn(II) can be interchanged in enzymes from different organisms to catalyze similar substrate reactions. In combination, quantitative electron paramagnetic resonance spectroscopy and rapid freeze-quench experiments allow us to follow the concentrations of four different Mn species, including key metal intermediates in the catalytic cycle, as the enzyme turns over its natural substrate. Two intermediates are observed: a Mn(III)-radical species which is either Mn-superoxide or Mn-substrate radical, and a unique Mn(II) species which is involved in the rate-limiting step of the cycle and may be Mn-alkylperoxo.     

The "somersault" mechanism for the p-450 hydroxylation of hydrocarbons. The intervention of transient inverted metastable hydroperoxides

A series of model theoretical calculations are described that suggest a new mechanism for the oxidation step in enzymatic cytochrome P450 hydroxylation of saturated hydrocarbons. A new class of metastable metal hydroperoxides is described that involves the rearrangement of the ground-state metal hydroperoxide to its inverted isomeric form with a hydroxyl radical hydrogen bonded to the metal oxide (MO-OH --> MO....HO). The activation energy for this somersault motion of the FeO-OH group is 20.3 kcal/mol for the P450 model porphyrin iron(III) hydroperoxide [Por(SH)Fe(III)-OOH(-)] to produce the isomeric ferryl oxygen hydrogen bonded to an *OH radical [Por(SH)Fe(III)-O....HO(-)]. This isomeric metastable hydroperoxide, the proposed primary oxidant in the P450 hydroxylation reaction, is calculated to be 17.8 kcal/mol higher in energy than the ground-state iron(III) hydroperoxide Cpd 0. The first step of the proposed mechanism for isobutane oxidation is abstraction of a hydrogen atom from the C-H bond of isobutane by the hydrogen-bonded hydroxyl radical to produce a water molecule strongly hydrogen bonded to anionic Cpd II. The hydroxylation step involves a concerted but nonsynchronous transfer of a hydrogen atom from this newly formed, bound, water molecule to the ferryl oxygen with a concomitant rebound of the incipient *OH radical to the carbon radical of isobutane to produce the C-O bond of the final product, tert-butyl alcohol. The TS for the oxygen rebound step is 2 kcal/mol lower in energy than the hydrogen abstraction TS (DeltaE() = 19.5 kcal/mol). The overall proposed new mechanism is consistent with a lot of the ancillary experimental data for this enzymatic hydroxylation reaction.     

Basic ancillary ligands promote O-O bond formation in iridium-catalyzed water oxidation: a DFT study

The cationic iridium complex [Ir(OH(2))(2)(phpy)(2)](+) (phpy = o-phenylpyridine) is among the most efficient mononuclear catalysts for water oxidation. The postulated active species is the oxo complex [Ir(O)(X)(phpy)(2)](n), with X = OH(2) (n = +1), OH(-) (n = 0) or O(2-) (n = -1), depending on the pH. The reactivity of these species has been studied computationally at the DFT(B3LYP) level. The three [Ir(O)(X)(phpy)(2)](n) complexes have an electrophilic Ir(v)-oxo moiety, which yields an O-O bond by undergoing a nucleophilic attack of water in the critical step of the mechanism. In this step, water transfers one proton to either the Ir(V)-oxo moiety or the ancillary X ligand. Five different reaction pathways associated with this acid/base mechanism have been characterized. The calculations show that the proton is preferably accepted by the X ligand, which plays a key role in the reaction. The higher the basicity of X, the lower the energy barrier associated with O-O bond formation. The anionic species, [Ir(O)(2)(phpy)(2)](-), which has the less electrophilic Ir(V)-oxo moiety but the most basic X ligand, promotes O-O bond formation through the lowest energy barrier, 14.5 kcal mol(-1). The other two active species, [Ir(O)(OH)(phpy)(2)] and [Ir(O)(OH(2))(phpy)(2)](+), which have more electrophilic Ir(V)-oxo moieties but less basic X ligands, involve higher energy barriers, 20.2 kcal mol(-1) and 25.9 kcal mol(-1), respectively. These results are in good agreement with experiments showing important pH effects in similar catalytic systems. The theoretical insight given by the present study can be useful in the design of more efficient water oxidation catalysts. The catalytic activity may increase by using ligand scaffolds bearing internal bases.     

New model for a theoretical density functional theory investigation of the mechanism of the carbonic anhydrase: how does the internal bicarbonate rearrangement occur?

A theoretical density functional theory (DFT, B3LYP) investigation has been carried out on the catalytic cycle of the carbonic anhydrase. A model system including the Glu106 and Thr199 residues and the "deep" water molecule has been used. It has been found that the nucleophilic attack of the zinc-bound OH on the CO(2) molecule has a negligible barrier (only 1.2 kcal mol(-1)). This small value is due to a hydrogen-bond network involving Glu106, Thr199, and the deep water molecule. The two usually proposed mechanisms for the internal bicarbonate rearrangement have been carefully examined. In the presence of the two Glu106 and Thr199 residues, the direct proton transfer (Lipscomb mechanism) is a two-step process, which proceeds via a proton relay network characterized by two activation barriers of 4.4 and 9.0 kcal mol(-1). This pathway can effectively compete with a rotational mechanism (Lindskog mechanism), which has a barrier of 13.2 kcal mol(-1). The fast proton transfer found here is basically due to the effect of the Glu106 residue, which stabilizes an intermediate situation where the Glu106 fragment is protonated. In the absence of Glu106, the barrier for the proton transfer is much larger (32.3 kcal mol(-1)) and the Lindskog mechanism becomes favored.