Electronic state of the molecular oxygen released by catalase

In catalases, the high redox intermediate known as compound I (Cpd I) is reduced back to the resting state by means of hydrogen peroxide in a 2-electron reaction [Cpd I (Por(*+)-Fe(IV)O) + H(2)O(2) --> Enz (Por-Fe(III)) + H(2)O + O(2)]. It has been proposed that this reaction takes place via proton transfer toward the distal His and hydride transfer toward the oxoferryl oxygen (H(+)/H(-) scheme) and some authors have related it to singlet oxygen generation. Here, we consider the possible reaction schemes and qualitatively analyze the electronic state of the species involved to show that the commonly used association of the H(+)/H(-) scheme with singlet oxygen production is not justified. The analysis is complemented with density functional theory (DFT) calculations for a gas-phase active site model of the reactants and products.     

The structures and electronic configuration of compound I intermediates of Helicobacter pylori and Penicillium vitale catalases determined by X-ray crystallography and QM/MM density functional theory calculations

The structures of Helicobacter pylori (HPC) and Penicillium vitale (PVC) catalases, each with two subunits in the crystal asymmetric unit, oxidized with peroxoacetic acid are reported at 1.8 and 1.7 A resolution, respectively. Despite the similar oxidation conditions employed, the iron-oxygen coordination length is 1.72 A for PVC, close to what is expected for a Fe=O double bond, and 1.80 and 1.85 A for HPC, suggestive of a Fe-O single bond. The structure and electronic configuration of the oxoferryl heme and immediate protein environment is investigated further by QM/MM density functional theory calculations. Four different active site electronic configurations are considered, Por*+-FeIV=O, Por*+-FeIV=O...HisH+, Por*+-FeIV-OH+ and Por-FeIV-OH (a protein radical is assumed in the latter configuration). The electronic structure of the primary oxidized species, Por*+-FeIV=O, differs qualitatively between HPC and PVC with an A2u-like porphyrin radical delocalized on the porphyrin in HPC and a mixed A1u-like "fluctuating" radical partially delocalized over the essential distal histidine, the porphyrin, and, to a lesser extent, the proximal tyrosine residue. This difference is rationalized in terms of HPC containing heme b and PVC containing heme d. It is concluded that compound I of PVC contains an oxoferryl Por*+-FeIV=O species with partial protonation of the distal histidine and compound I of HPC contains a hydroxoferryl Por-FeIV-OH with the second oxidation equivalent delocalized as a protein radical. The findings support the idea that there is a relation between radical migration to the protein and protonation of the oxoferryl bond in catalase.     

On the functional role of a water molecule in clade 3 catalases: a proposal for the mechanism by which NADPH prevents the formation of compound II

X-ray structures of the 13 different monofunctional heme catalases published to date were scrutinized in order to gain insight in the mechanism by which NADPH in Clade 3 catalases may protect the reactive ferryloxo intermediate Compound I (Cpd I; por (*+)Fe (IV)O) against deactivation to the catalytically inactive intermediate Compound II (Cpd II; porFe (IV)O). Striking similarities in the molecular network of the protein subunits encompassing the heme center and the surface-bound NADPH were found for all of the Clade 3 catalases. Unique features in this region are the presence of a water molecule (W1) adjacent to the 4-vinyl group of heme and a serine residue or a second water molecule hydrogen-bonded to both W1 and the carbonyl group of a threonine-proline linkage, with the proline in van der Waals contact with the dihydronicotinamide group of NADPH. A mechanism is proposed in which a hydroxyl anion released from W1 undergoes reversible nucleophilic addition to the terminal carbon of the 4-vinyl group of Cpd I, thereby producing a neutral porphyrin pi-radical ferryloxo (HO-por (*)Fe (IV)O) species of reduced reactivity. This structure is suggested to be the elusive Cpd II' intermediate proposed in previous studies. An accompanying proton-shifting process along the hydrogen-bonded network is believed to facilitate the NADPH-mediated reduction of Cpd I to ferricatalase and to serve as a funnel for electron transfer from NADPH to the heme center to restore the catalase Fe (III) resting state. The proposed reaction paths were fully supported as chemically reasonable and energetically feasible by means of density functional theory calculations at the (U)B3LYP/6-31G* level. A particularly attractive feature of the present mechanism is that the previously discussed formation of protein-derived radicals is avoided.     

Versatility of the electronic structure of compound I in catalase-peroxidases

Catalase-peroxidases (KatGs) are bifunctional heme proteins, belonging to the family of class I peroxidases, that are able to catalyze both catalatic and peroxidatic reactions within a peroxidase-like structure. We investigated the electronic structure of reaction intermediates of the catalytic cycle of KatGs by means of density functional theory (DFT) QM/MM calculations. The outcome was that the ionization state of the KatG-specific covalent adduct (Met264-Tyr238-Trp111) affects the radical character of compound I (Cpd I). Specifically, in the optimized structures, substantial radical character is observed on the proximal Trp330 when Tyr238 is protonated, whereas when Tyr238 is deprotonated the radical localizes on the Met+-Tyr(O-)-Trp adduct. These findings are not affected by protein thermal fluctuations, although details of the spin density distribution are affected by the geometry of the active site. Calculations provide structures in good agreement with the crystal structure of BpKatG Cpd I. They also provide an explanation for the experimental findings of the mobile and catalatic-specific residue Arg426 being 100% in conformation R in the X-ray structure of BpKatG treated with organic peroxides. The role of different Cpd I forms in the catalase and peroxidase reaction pathways is discussed.     

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.     

Theoretical study of the mechanisms of substrate recognition by catalase

A variety of theoretical methods including classical molecular interaction potentials, classical molecular dynamics, and activated molecular dynamics have been used to analyze the substrate recognition mechanisms of peroxisomal catalase from Saccharomyces cerevisiae. Special attention is paid to the existence of channels connecting the heme group with the exterior of the protein. On the basis of these calculations a rationale is given for the unique catalytic properties of this enzyme, as well as for the change in enzyme efficiency related to key mutations. According to our calculations the water is expected to be a competitive inhibitor of the enzyme, blocking the access of hydrogen peroxide to the active site. The main channel is the preferred route for substrate access to the enzyme and shows a cooperative binding to hydrogen peroxide. However, the overall affinity of the main channel for H(2)O(2) is only slightly larger than that for H(2)O. Alternative channels connecting the heme group with the monomer interface and the NADP(H) binding site are detected. These secondary channels might be important for product release.     

Fe(III) complexes of 1,4,8,11-tetraaza[14]annulenes as catalase mimics

The development of enzyme mimics of catalase which decompose hydrogen peroxide to water and molecular oxygen according to the 2:1 stoichiometry of native catalase and in aqueous solution at pH 7 and at micromolar concentrations of the enzyme model and hydrogen peroxide is reported. For this purpose, iron(III) complexes of 1,4,8,11-tetraaza[14]annulenes are prepared by various procedures. Efficacious preparations utilize reaction of the [N4] macrocyles with FeII salts in the presence of triphenylamine, followed by gentle oxidation of the FeII complexes by molecular oxygen or by tris(4-bromophenyl)aminium hexachloroantimonate. The complexes are characterized by SQUID magnetometry and by M?ssbauer, EPR, and UV/vis spectrometry. In the solid state, the iron(III) center of the catalytically active complexes exists in the intermediate (quartet, S = 3/2) spin state. Several of these complexes decompose hydrogen peroxide in aqueous buffer solution at pH 7.2 at room temperature with turnover numbers between 40 and 80. The apparent second-order rate constant for hydrogen peroxide decomposition is in the range of 1400-2400 M(-1) s(-1), about 3 orders of magnitude lower than the value for native catalase. Besides oxygen production, a non-oxygen releasing pathway of hydrogen peroxide decomposition is unveiled.     

Quantum mechanical/molecular mechanical study on the mechanisms of compound I formation in the catalytic cycle of chloroperoxidase: an overview on heme enzymes

The formation of Compound I (Cpd I), the active species of the enzyme chloroperoxidase (CPO), was studied using QM/MM calculation. Starting from the substrate complex with hydrogen peroxide, FeIII-HOOH, we examined two alternative mechanisms on the three lowest spin-state surfaces. The calculations showed that the preferred pathway involves heterolytic O-O cleavage that proceeds via the iron hydroperoxide species, i.e., Compound 0 (Cpd 0), on the doublet-state surface. This process is effectively concerted, with a barrier of 12.4 kcal/mol, and is catalyzed by protonation of the distal OH group of Cpd 0. By comparison, the path that involves a direct O-O cleavage from FeIII-HOOH is less favored. A proton coupled electron transfer (PCET) feature was found to play an important role in the mechanism nascent from Cpd 0. Initially, the O-O cleavage progresses in a homolytic sense, but as soon as the proton is transferred to the distal OH, it triggers an electron transfer from the heme-oxo moiety to form water and Cpd I. This study enables us to generalize the mechanisms of O-O activation, elucidated so far by QM/MM calculations, for other heme enzymes, e.g., cytochrome P450cam, horseradish peroxidase (HRP), nitric oxide synthase (NOS), and heme oxygenase (HO). Much like for CPO, in the cases of P450 and HRP, the PCET lowers the barrier below the purely homolytic cleavage alternative (in our case, the homolytic mechanism is calculated directly from FeIII-HOOH). By contrast, the absence of PCET in HO, along with the robust water cluster, prefers a homolytic cleavage mechanism.     

Direct electrochemistry of catalase on glassy carbon electrodes.

Catalase was investigated as a possible catalyst of the electrochemical reduction of oxygen on glassy carbon electrodes. The presence of catalase dissolved in solution only provoked a moderate current increase, which was fully explained by the catalase-catalysed disproportionation of hydrogen peroxide (Scheme I). When catalase was adsorbed from dimethylsulfoxide on the surface of electrodes that did not undergo any electrochemical pre-treatment (EP), catalase efficiently catalysed oxygen reduction via direct electron transfer from the electrode (Scheme II). The results are discussed with respect to the electrode surface properties and the enzyme structure.     

Hydrogen peroxide decomposition by a non-heme iron(III) catalase mimic: a DFT study

Non-heme iron(III) complexes of 14-membered tetraaza macrocycles have previously been found to catalytically decompose hydrogen peroxide to water and molecular oxygen, like the native enzyme catalase. Here the mechanism of this reaction is theoretically investigated by DFT calculations at the (U)B3LYP/6-31G* level, with focus on the reactivity of the possible spin states of the FeIII complexes. The computations suggest that H2O2 decomposition follows a homolytic route with intermediate formation of an iron(IV) oxo radical cation species (L.+FeIV==O) that resembles Compound I of natural iron porphyrin systems. Along the whole catalytic cycle, no significant energetic differences were found for the reaction proceeding on the doublet (S=1/2) or on the quartet (S=3/2) hypersurface, with the single exception of the rate-determining O--O bond cleavage of the first associated hydrogen peroxide molecule, for which reaction via the doublet state is preferred. The sextet (S=5/2) state of the FeIII complexes appears to be unreactive in catalase-like reactions.     

Compound I in heme thiolate enzymes: a comparative QM/MM study

This study directly compares the active species of heme enzymes, so-called Compound I (Cpd I), across the heme-thiolate enzyme family. Thus, sixty-four different Cpd I structures are calculated by hybrid quantum mechanical/molecular mechanical (QM/MM) methods using four different cysteine-ligated heme enzymes (P450(cam), the mutant P450(cam)-L358P, CPO and NOS) with varying QM region sizes in two multiplicities each. The overall result is that these Cpd I species are similar to each other with regard to many characteristic features. Hence, using the more stable CPO Cpd I as a model for P450 Cpd I in experiments should be a reasonable approach. However, systematic differences were also observed, and it is shown that NOS stands out in most comparisons. By analyzing the electrical field generated by the enzyme on the QM region, one can see that (a) the protein exerts a large influence and modifies all the Cpd I species compared with the gas-phase situation and (b) in NOS this field is approximately planar to the heme plane, whereas it is approximately perpendicular in the other enzymes, explaining the deviating results on NOS. The calculations on the P450(cam) mutant L358P show that the effects of removing the hydrogen bond between the heme sulfur and L358 are small at the Cpd I stage. Finally, Mossbauer parameters are calculated for the different Cpd I species, enabling future comparisons with experiments. These results are discussed in the broader context of recent findings of Cpd I species that exhibit large variations in the electronic structure due to the presence of the substrate.     

QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding

Cytochrome P450 3A4 is involved in the metabolism of 50% of all swallowed drugs. The enzyme functions by means of a high-valent iron-oxo species, called compound I (Cpd I), which is formed after entrance of the substrate to the active site. We explored the features of Cpd I using hybrid quantum mechanical/molecular mechanical calculations on various models that are either substrate-free or containing one and two molecules of diazepam as a substrate. M?ssbauer parameters of Cpd I were computed. Our major finding shows that without the substrate, Cpd I tends to elongate its Fe-S bond, localize the radical on the sulfur, and form hydrogen bonds with A305 and T309, which may hypothetically lead to Cpd I consumption by H-abstraction. However, the positioning of diazepam close to Cpd I, as enforced by the effector molecule, was found to strengthen the NH...S interactions of the conserved I443 and G444 residues with the proximal cysteinate ligand. These interactions are known to stabilize the Fe-S bond, and as such, the presence of the substrate leads to a shorter Fe-S bond and it prevents the localization of the radical on the sulfur. This diazepam-Cpd I stabilization was manifested in the 1W0E conformer. The effector substrate did not influence Cpd I directly but rather by positioning the active substrate close to Cpd I, thus displacing the hydrogen bonds with A305 and T309, and thereby giving preference to substrate oxidation. It is hypothesized that these effects on Cpd I, promoted by the restrained substrate, may be behind the special metabolic behavior observed in cases of multiple substrate binding (also called cooperative binding). This restraint constitutes a mechanism whereby substrates stabilize Cpd I sufficiently long to affect monooxygenation by P450s at the expense of Cpd I destruction by the protein residues.     

Compound I of nitric oxide synthase: the active site protonation state

A quantum mechanical/molecular mechanical (QM/MM) study of the formation of the elusive active species Compound I (Cpd I) of nitric oxide synthase (NOS) from the oxyferrous intermediate shows that two protons have to be provided to produce a reaction that is reasonably exothermic and that leads to the appearance of a radical on the tetrahydrobiopterin cofactor. Molecular dynamics and energy considerations show that a possible source of proton is the water H-bond chain formed from the surface to the active site, but that a water molecule by itself cannot be the source of the proton; an H3O+ species that is propagated along the chain is more likely. The QM/MM calculations demonstrate that Cpd I and H2O are formed from the ferric-hydrogen peroxide complex in a unique heterolytic O-O cleavage mechanism. The properties of the so-formed Cpd I are compared with those of the known species of chloroperoxidase, and the geometry and spin densities are found to be compatible. The M?ssbauer parameters are calculated and may serve as experimental probes in attempts to characterize NOS Cpd I.     

微量热法研究过氧化氢酶底物抑制动力学

过氧化氢酶是需氧生物体内抗氧化酶系的重要组分。过氧化氢酶催化过氧化氢分解是一个两底物酶促反应,依照Chance提出的机理,反应速率方程具有一级反应方程的形式。此反应在高浓度底物存在的情况下,表现出明显的不可逆底物抑制。本研究用热动力学方法研究了这一反应,提出了一种不可逆底物抑制机理,并应用该机理求出了相关动力学参数。在310.15K,pH=7.0时k0=9.6×10^5L·mol^-1^·s^-1,k1/k2=2.9×10^6。实验结果证明此机理正确有效。     

(salen)MnIII compounds as nonpeptidyl mimics of catalase. Mechanism-based tuning of catalase activity: a theoretical study

We present the results of the first theoretical investigation of salen-manganese complexes as synthetic catalytic scavengers of hydrogen peroxide molecules that mimic catalase enzymes. Catalase mimics can be used as therapeutic agents against oxidative stress in treatment of many diseases, including Alzheimer's disease, stroke, heart disease, aging, and cancer. A ping-pong mechanism approach has been considered to describe the H2O2 dismutation reaction. The real compounds reacting with a peroxide molecule were utilized in our BP density functional calculations to avoid uncertainties connected with using incomplete models. Part I of the dismutation reaction-converting a peroxide molecule into a water molecule with simultaneous oxidation of the metal atom of the catalyst-can be done quite effectively at the Mn catalytic center. To act as catalytic scavengers of hydrogen peroxide, the oxomanganese salen complexes have to be deoxidized during part II of the dismutation reaction. It has been shown that there are two possible reaction routes for the second part of the dismutation reaction: the top and the side substrate approach routes. Our results suggest that the catalyst could be at least temporarily deactivated (poisoned) in the side approach reaction route due to the formation of a kinetically stable intermediate. Overall, the side approach reaction route for the catalyst recovery is the bottleneck for the whole dismutation process. On the basis of the detailed knowledge of the mode of action of the (salen)MnIII catalase mimics, we suggest and rationalize structural changes of the catalyst that should lead to better therapeutic properties. The available experimental data support our conclusions. Our findings on the reaction dismutation mechanism could be the starting point for further improvement of salen-manganese complexes as synthetic catalytic scavengers of reactive oxygen species.     

Active species of horseradish peroxidase (HRP) and cytochrome P450: two electronic chameleons

The active site of HRP Compound I (Cpd I) is modeled using hybrid density functional theory (UB3LYP). The effects of neighboring amino acids and of environmental polarity are included. The low-lying states have porphyrin radical cationic species (Por(*)(+)). However, since the Por(*)(+) species is a very good electron acceptor, other species, which can be either the ligand or side chain amino acid residues, may participate in electron donation to the Por(*)(+) moiety, thereby making Cpd I behave like a chemical chameleon. Thus, this behavior that was noted before for Cpd I of P450 is apparently much more wide ranging than initially appreciated. Since chemical chameleonic behavior property was found to be expressed not only in the properties of Cpd I itself, but also in its reactivity, the roots of this phenomenon are generalized. A comparative discussion of Cpd I species follows for the enzymes HRP, CcP, APX, CAT (catalase), and P450.     

Flow-injection assay of catalase activity.

A novel flow-injection assay (FIA) system with a double line for catalase activity was constructed in which an oxidase is immobilized and the substrate is continuously pumped to reduce the dissolved oxygen and to generate a given level of hydrogen peroxide. The catalase in a sample decomposed the hydrogen peroxide, and thus the increase in dissolved oxygen dependent on the activity was amperometrically monitored using a Clark-type oxygen electrode. Among the examined several oxidases, uricase was most suitable for the continuous formation of hydrogen peroxide from a consideration of the stability and the conversion efficiency. Under the optimum conditions, a linear calibration curve was obtained in the range from 21 to 210 units/mg and the reproducibility (CV) was better than 2% by 35 successive determinations of 210 units/ml catalase preparation. The sampling frequency was about 15 samples/h. The present FIA system was applicable to monitor the inactivation of catalase by glycation.     

DFT Calculations Suggest a New Type of Self-Protection and Self-Inhibition Mechanism in the Mammalian Heme Enzyme Myeloperoxidase: Nucleophilic Addition of a Functional Water rather than One-Electron Reduction

The mammalian heme enzyme myeloperoxidase (MPO) catalyzes the reaction of Cl(-) to the antimicrobial-effective molecule HOCl. During the catalytic cycle, a reactive intermediate "Compound?I" (Cpd?I) is generated. Cpd?I has the ability to destroy the enzyme. Indeed, in the absence of any substrate, Cpd?I decays with a half-life of 100?ms to an intermediate called Compound?II (Cpd?II), which is typically the one-electron reduced Cpd?I. However, the nature of Cpd?II, its spectroscopic properties, and the source of the additional electron are only poorly understood. On the basis of DFT and time-dependent (TD)-DFT quantum chemical calculations at the PBE0/6-31G* level, we propose an extended mechanism involving a new intermediate, which allows MPO to protect itself from self-oxidation or self-destruction during the catalytic cycle. Because of its similarity in electronic structure to Cpd?II, we named this intermediate Cpd?II'. However, the suggested mechanism and our proposed functional structure of Cpd?II' are based on the hypothesis that the heme is reduced by charge separation caused by reaction with a water molecule, and not, as is normally assumed, by the transfer of an electron. In the course of this investigation, we found a second intermediate, the reduced enzyme, towards which the new mechanism is equally transferable. In analogy to Cpd?II', we named it Fe(II') . The proposed new intermediates Cpd?II' and Fe(II') allow the experimental findings, which have been well documented in the literature for decades but not so far understood, to be explained for the first time. These encompass a)?the spontaneous decay of Cpd?I, b)?the unusual (chlorin-like) UV/Vis, circular dichroism (CD), and resonance Raman spectra, c)?the inability of reduced MPO to bind CO, d)?the fact that MPO-Cpd?II reduces SCN(-) but not Cl(-) , and e)?the experimentally observed auto-oxidation/auto-reduction features of the enzyme. Our new mechanism is also transferable to cytochromes, and could well be viable for heme enzymes in general.     

Soluble and immobilized catalase

This article examines the effect of pressure on the steady-state kinetics and long-term deactivation of the enzyme catalase supported on porous alumina. The reaction studied is the decomposition of hydrogen peroxide. The results of studies carried out in a continuous stirred-tank reactor under isothermal conditions are presented and compared with results obtained for soluble catalase. For soluble catalase, it is found that in the range of pressures studied, the oxygen flow rate increases with increase in pressure up to a certain value and then decreases. Hydrogen peroxide concentration appears to have a strong influence on pressure effects. With immobilized catalase, the pressure effects are not as prominent. Fluorescent microscopy studies of the immobilized enzyme suggest that this is probably because of pore diffusional limitations.