What External Perturbations Influence the Electronic Properties of Catalase Compound I?

Theoretical Study of Hydroxylation of α- and β-Pinene by a Cytochrome P450 Monooxygenase Model

Previous studies on biocatalytic transformations of pinenes by cytochrome P450 (CYP) enzymes reveal the formation of different oxygenated products from a single substrate due to the multistate reactivity of CYP and the many reactive sites in the pinene scaffold. Up until now, the detailed mechanism of these biocatalytic transformations of pinenes have not been reported. Hereby, we report a systematic theoretical study of the plausible hydrogen abstraction and hydroxylation reactions of α- and β-pinenes by CYP using the density functional theory (DFT) method. All DFT calculations in this study were based on B3LYP/LAN computational methodology using the Gaussian09 software. We used the B3LYP functional with corrections for dispersive forces, BSSE, and anharmonicity to study the mechanism and thermodynamic properties of these reactions using a bare model (without CYP) and a pinene-CYP model. According to the potential energy surface and Boltzmann distribution for radical conformers, the major reaction products of CYP-catalyzed hydrogen abstraction from β-pinene are the doublet trans (53.4%) and doublet cis (46.1%) radical conformer at delta site. The formation of doublet cis/trans hydroxylated products released a total Gibbs free energy of about 48 kcal/mol. As for alpha pinene, the most stable radicals were trans-doublet (86.4%) and cis-doublet (13.6%) at epsilon sites, and their hydroxylation products released a total of ~50 kcal/mol Gibbs free energy. Our results highlight the likely C-H abstraction and oxygen rebounding sites accounting for the multi-state of CYP (doublet, quartet, and sextet spin states) and the formation of different conformers due to the presence of cis/trans allylic hydrogen in α-pinene and β-pinene molecules.

Monofunctional Heme-Catalases

The review focuses on four issues that are critical for the understanding of monofunctional catalases. How hydrogen peroxide (H2O2) reaches the active site and outcompetes water molecules to be able to function at a very high rate is one of the issues examined. Part of the answer is a gate valve system that is instrumental to drive out solvent molecules from the final section of the main channel. A second issue relates to how the enzyme deals with an unproductive reactive compound I (Cpd I) intermediate. Peroxidatic two and one electron donors and the transfer of electrons to the active site from NADPH and other compounds are reviewed. The new ascribed catalase reactions are revised, indicating possible measurement pitfalls. A third issue concerns the heme b to heme d oxidation, why this reaction occurs only in some large-size subunit catalases (LSCs), and the possible role of singlet oxygen in this and other modifications. The formation of a covalent bond between the proximal tyrosine with the vicinal residue is analyzed. The last issue refers to the origin and function of the additional C-terminal domain (TD) of LSCs. The TD has a molecular chaperone activity that is traced to a gene fusion between a Hsp31-type chaperone and a small-size subunit catalase (SSC).

Regulation of antioxidant defense and glyoxalase systems in cyanobacteria

Oxidative stress is common consequence of abiotic stress in plants as well as cyanobacteria caused by generation of reactive oxygen species (ROS), an inevitable product of respiration and photosynthetic electron transport. ROS act as signalling molecule at low concentration however, when its production exceeds the endurance capacity of antioxidative defence system, the organisms suffer oxidative stress. A highly toxic metabolite, methylglyoxal (MG) is also produced in cyanobacteria in response to various abiotic stresses which consequently augment the ensuing oxidative damage. Taking recourse to the common lineage of eukaryotic plants and cyanobacteria, it would be worthwhile to explore the regulatory role of glyoxalase system and antioxidative defense mechanism in combating abiotic stress in cyanobacteria. This review provides comprehensive information on the complete glyoxalase system (GlyI, GlyII and GlyIII) in cyanobacteria. Furthermore, it elucidates the recent understanding regarding the production of ROS and MG, noteworthy link between intracellular MG and ROS and its detoxification via synchronization of antioxidants (enzymatic and non-enzymatic) and glyoxalase systems using glutathione (GSH) as common co-factor.

Modeling the hydroxylation of estragole via human liver cytochrome P450

Natural compounds derived from plants are generally regarded safe and devoid of adverse effects. However, there are individual ingredients that possess toxic, genotoxic, and carcinogenic activities. These compounds when exposed at specific level become hazardous to health. Estragole (1-allyl-4-methoxybenzene) is a common component of spice plants. Its toxicity gets activated with the hydroxylation at benzylic carbon (C1') position by P450 enzymes present in the human liver. The present study grounds to explore the reaction mechanism of conversion of estragole to hydroxylated metabolite using computational methodology. Density functional theory (DFT)-based calculations were employed to explore the cytochrome P450-catalyzed mechanism at C1 position aliphatic hydroxylation of estragole. Overall reaction energy profile, electronic configuration, and 3D structure of all intermediates, transition states, and product complexes formed during the reaction along with their free energies were tried to be investigated.

Local Electric Fields as a Natural Switch of Heme-Iron Protein Reactivity

Heme-iron oxidoreductases operating through the high-valent Fe^IVO intermediates perform crucial and complicated transformations, such as oxidations of unreactive saturated hydrocarbons. These enzymes share the same Fe coordination, only differing by the axial ligation, e.g., Cys in P450 oxygenases, Tyr in catalases, and His in peroxidases. By examining ~200 heme-iron proteins, we show that the protein hosts exert highly specific intramolecular electric fields on the active sites, and there is a strong correlation between the direction and magnitude of this field and the protein function. In all heme proteins, the field is preferentially aligned with the Fe-O bond ( F_z ). The Cys-ligated P450 oxygenases have the highest average F_z of 28.5 MV cm^-1, i.e., most enhancing the oxyl-radical character of the oxo group, and consistent with the ability of these proteins to activate strong C-H bonds. In contrast, in Tyr-ligated proteins, the average F_z is only 3.0 MV cm^-1, apparently suppressing single-electron off-pathway oxidations, and in His-ligated proteins, F_z is -8.7 MV cm^-1. The operational field range is given by the trade-off between the low reactivity of the Fe^IVO Compound I at the more negative F_z , and the low selectivity at the more positive F_z . Consequently, a heme-iron site placed in the field characteristic of another heme-iron protein class loses its canonical function, and gains an adverse one. Thus, electric fields produced by the protein scaffolds, together with the nature of the axial ligand, control all heme-iron chemistry.

DFT study on the effect of proximal residues on the Mycobacterium tuberculosis catalase-peroxidase (katG) heme compound I intermediate and its bonding interaction with isoniazid

Isoniazid (INH) is converted into isonicotinyl radical through its interaction with the catalase-peroxidase (katG) enzyme present in the cells of Mycobacterium tuberculosis (M. tb.), the bacteria that causes the tuberculosis disease. This process is important because resistance of M. tb. cells to INH treatment has been associated with the failure of completion of this process. However, this process is poorly understood and there are a variety of conflicting theories about the details of the mechanism of this process. One theory suggests that INH binds to katG and transfers a single electron to the heme while the heme is in its two electron oxidized state, compound I [Fe(iv)Por˙+] (CpdI). In this study, DFT calculations at the UB3LYP/6-31g(d)-LANL2DZ level of theory are used to study the M. tb. katG CpdI molecule. Different models of the M. tb. CpdI molecule were prepared and the calculations revealed the impact of Trp321 on the electronic properties of the heme. Without Trp321 the heme assumed a triradical state with single electrons on the πxy and πyz orbitals of Fe and another on the a2u orbital of the porphyrin ring that can either be coupled with the first two, to assume a quartet state, or decoupled to form a doublet state. With Trp321, however, a transfer of an electron from the πTrp orbital to a2u porphyrin orbital leads to loss of radical character of the porphyrin and leaves the Trp321 group with a radical character. INH was observed to have strong interaction with CpdI and the absence of Trp321 significantly decreased the binding energy by 2 kcal mol-1 explaining the importance of Trp321 in the binding of INH. The results in this study show the importance of Trp321 in the binding of INH and its effect on its electronic properties, which is consistent with previous experimental findings that mutation of Trp321 results in an increase in drug resistance.

Challenging Density Functional Theory Calculations with Hemes and Porphyrins

In this paper we review recent advances in computational chemistry and specifically focus on the chemical description of heme proteins and synthetic porphyrins that act as both mimics of natural processes and technological uses. These are challenging biochemical systems involved in electron transfer as well as biocatalysis processes. In recent years computational tools have improved considerably and now can reproduce experimental spectroscopic and reactivity studies within a reasonable error margin (several kcal·mol(-1)). This paper gives recent examples from our groups, where we investigated heme and synthetic metal-porphyrin systems. The four case studies highlight how computational modelling can correctly reproduce experimental product distributions, predicted reactivity trends and guide interpretation of electronic structures of complex systems. The case studies focus on the calculations of a variety of spectroscopic features of porphyrins and show how computational modelling gives important insight that explains the experimental spectra and can lead to the design of porphyrins with tuned properties.

A comprehensive test set of epoxidation rate constants for iron(iv)-oxo porphyrin cation radical complexes

Cytochrome P450 enzymes are heme based monoxygenases that catalyse a range of oxygen atom transfer reactions with various substrates, including aliphatic and aromatic hydroxylation as well as epoxidation reactions. The active species is short-lived and difficult to trap and characterize experimentally, moreover, it reacts in a regioselective manner with substrates leading to aliphatic hydroxylation and epoxidation products, but the origin of this regioselectivity is poorly understood. We have synthesized a model complex and studied it with low-pressure Fourier transform-ion cyclotron resonance (FT-ICR) mass spectrometry (MS). A novel approach was devised using the reaction of [FeIII(TPFPP)]+ (TPFPP = meso-tetrakis(pentafluorophenyl)porphinato dianion) with iodosylbenzene as a terminal oxidant which leads to the production of ions corresponding to [FeIV(O)(TPFPP+˙)]+. This species was isolated in the gas-phase and studied in its reactivity with a variety of olefins. Product patterns and rate constants under Ideal Gas conditions were determined by FT-ICR MS. All substrates react with [FeIV(O)(TPFPP+˙)]+ by a more or less efficient oxygen atom transfer process. In addition, substrates with low ionization energies react by a charge-transfer channel, which enabled us to determine the electron affinity of [FeIV(O)(TPFPP+˙)]+ for the first time. Interestingly, no hydrogen atom abstraction pathways are observed for the reaction of [FeIV(O)(TPFPP+˙)]+ with prototypical olefins such as propene, cyclohexene and cyclohexadiene and also no kinetic isotope effect in the reaction rate is found, which suggests that the competition between epoxidation and hydroxylation - in the gas-phase - is in favour of substrate epoxidation. This notion further implies that P450 enzymes will need to adapt their substrate binding pocket, in order to enable favourable aliphatic hydroxylation over double bond epoxidation pathways. The MS studies yield a large test-set of experimental reaction rates of iron(iv)-oxo porphyrin cation radical complexes, so far unprecedented in the gas-phase, providing a benchmark for calibration studies using computational techniques. Preliminary computational results presented here confirm the observed trends excellently and rationalize the reactivities within the framework of thermochemical considerations and valence bond schemes.

The reaction mechanisms of heme catalases: an atomistic view by ab initio molecular dynamics

Catalases are ubiquitous enzymes that prevent cell oxidative damage by degrading hydrogen peroxide to water and oxygen (2H(2)O(2) → 2H(2)O+O(2)) with high efficiency. The enzyme is first oxidized to a high-valent iron intermediate, known as Compound I (Cpd I, Por(·+)-Fe(IV)=O) which, at difference from other hydroperoxidases, is reduced back to the resting state by further reacting with H(2)O(2). The normal catalase activity is reduced if Cpd I is consumed in a competing side reaction, forming a species named Cpd I*. In recent years, Density Functional Theory (DFT) methods have unraveled the electronic configuration of these high-valent iron species, helping to assign the intermediates trapped in the crystal structures of oxidized catalases. It has been demonstrated that the a priori assumption that the H(+)/H(-) type of mechanism for Cpd I reduction leads to the generation of singlet oxygen is not justified. Moreover, it has been shown by ab initio metadynamics simulations that two pathways are operative for Cpd I reduction: a His-mediated mechanism (described as H·/H(+) + e(-)) in which the distal His acts as an acid-base catalyst and a direct mechanism (described as H·/H·) in which the distal His does not play a direct role. Independently of the mechanism, the reaction proceeds by two one-electron transfers rather than one two-electron transfer, as previously assumed. Electron transfer to Cpd I, regardless of whether the electron is exogenous or endogenous, facilitates protonation of the oxoferryl group, to the point that formation of Cpd I* may be controlled by the easiness of protonation of reduced Cpd I.

Spectroscopic and theoretical studies of Ga(III)protoporphyrin-IX and its reactions with myoglobin

Ga(III)protoporphyrin-IX (Ga-PP) has been proposed as a model for the key interporphyrin interactions in malaria pigment. Unlike the paramagnetic parent iron heme derivatives, Ga-PP is readily soluble in methanol (MeOH). We report optical, mass spectroscopic, and theoretical results for Ga-PP as well as its reactions with myoglobin. UV-visible absorption and MCD spectroscopy show that Ga-PP exhibits a typical spectrum for a main group metal: a Q-band at 539 nm and a B band at 406 nm when dissolved in MeOH. We also report optical data for Zn(II)protoporphyrin IX (Zn-PP) dissolved in MeOH, which exhibits a Q-band at 545 nm and a B band at 415 nm. ESI mass spectral data for Ga-PP dissolved in MeOH show the presence of predominantly monomers, with smaller fractions of dimers [(Ga-PP)(2)] and trimers. UV-visible and MCD absorption spectroscopy and ESI mass spectral data demonstrate the successful insertion of monomeric Ga-PP into apo-Mb. Ga-PP-Mb exhibits a B band at 417 nm and Q bands at 545 and 584 nm, which are all red-shifted from the free Ga-PP values. The calculated electronic structures and frontier molecular orbitals of Ga-PP, (Ga-PP)(2) and Zn-PP fit the previously reported trends in band energies and oscillator strengths as a function of molecular orbital energies. These new data can be applied to explain the experimentally observed optical spectroscopy. The observed Q-band energies are accounted for by calculated (HOMO-LUMO) gap of the frontier MOs, while the split in the two top occupied MOs accounts for the magnitude of the Q-band oscillator strength as well as the experimentally observed Q to B band energy separation. Although Ga-PP shares more spectroscopic properties with Zn-PP than it does with Fe(III)PPIX, the trivalent oxidation state allows this molecule to be used as a model for ferric hemes in heme proteins.

Theoretical insights into heme-catalyzed oxidation of cyclohexane to adipic acid

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.

Evolution of catalases from bacteria to humans

Excessive hydrogen peroxide is harmful for almost all cell components, so its rapid and efficient removal is of essential importance for aerobically living organisms. Conversely, hydrogen peroxide acts as a second messenger in signal-transduction pathways. H(2)O(2) is degraded by peroxidases and catalases, the latter being able both to reduce H(2)O(2) to water and to oxidize it to molecular oxygen. Nature has evolved three protein families that are able to catalyze this dismutation at reasonable rates. Two of the protein families are heme enzymes: typical catalases and catalase-peroxidases. Typical catalases comprise the most abundant group found in Eubacteria, Archaeabacteria, Protista, Fungi, Plantae, and Animalia, whereas catalase-peroxidases are not found in plants and animals and exhibit both catalatic and peroxidatic activities. The third group is a minor bacterial protein family with a dimanganese active site called manganese catalases. Although catalyzing the same reaction (2 H(2)O(2)--> 2 H(2)O+ O(2)), the three groups differ significantly in their overall and active-site architecture and the mechanism of reaction. Here, we present an overview of the distribution, phylogeny, structure, and function of these enzymes. Additionally, we report about their physiologic role, response to oxidative stress, and about diseases related to catalase deficiency in humans.

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.