Theoretical study of model compound I complexes of horseradish peroxidase and catalase

Ground State and Excited State Tuning in Ferric Dipyrrin Complexes Promoted by Ancillary Ligand Exchange

Three ferric dipyrromethene complexes featuring different ancillary ligands were synthesized by one electron oxidation of ferrous precursors. Four-coordinate iron complexes of the type (^ArL)FeX₂ [^ArL = 1,9-(2,4,6-Ph₃C₆H₂)₂-5-mesityldipyrromethene] with X = Cl or ^tBuO were prepared and found to be high-spin (S = ⁵/₂), as determined by superconducting quantum interference device magnetometry, electron paramagnetic resonance, and ⁵⁷Fe Mössbauer spectroscopy. The ancillary ligand substitution was found to affect both ground state and excited properties of the ferric complexes examined. While each ferric complex displays reversible reduction and oxidation events, each alkoxide for chloride substitution results in a nearly 600 mV cathodic shift of the Fe^III/II couple. The oxidation event remains largely unaffected by the ancillary ligand substitution and is likely dipyrrin-centered. While the alkoxide substituted ferric species largely retain the color of their ferrous precursors, characteristic of dipyrrin-based ligand-to-ligand charge transfer (LLCT), the dichloride ferric complex loses the prominent dipyrrin chromophore, taking on a deep green color. Time-dependent density functional theory analyses indicate the weaker-field chloride ligands allow substantial configuration mixing of ligand-to-metal charge transfer into the LLCT bands, giving rise to the color changes observed. Furthermore, the higher degree of covalency between the alkoxide ferric centers is manifest in the observed reactivity. Delocalization of spin density onto the tert-butoxide ligand in (^ArL)FeCl(O^tBu) is evidenced by hydrogen atom abstraction to yield (^ArL)FeCl and HO^tBu in the presence of substrates containing weak C-H bonds, whereas the chloride (^ArL)FeCl₂ analogue does not react under these conditions.

Symmetry states of metalloporphyrin π-cation radicals, models for peroxidase compound I

Oxoferryl porphyrin π-cation radical active sites of compound I intermediates which are found in enzymes such as peroxidases and catalases have been extensively modeled by oxidized synthetic metalloporphyrins. The electronic symmetry states of these compounds were initially assigned on the basis of electronic absorption data. In recent years new experimental and theoretical results have become available which have led to a re-evaluation and modification of the original assignments. A historical perspective of these developments is provided in the context of recent NMR, resonance Raman, and other spectroscopic data and theoretical calculations for the synthetic models and enzymatic systems.

Proximal ligand effects on electronic structure and spectra of compound I of peroxidases

Computational studies exploring the extent to which differences in proximal axial ligands modulate structure, spectra, and function of peroxidases have been performed. To this end, three heme models of compound I were characterized differing only in the axial ligand. The axial ligands considered were L = ImH , Im-, that are alternative protonation models for a typical peroxidase with an imidazole ligand such as horseradish peroxidase (HRP-I), and L = SCH - that is a model for an unsual peroxidase, chloroperoxidase (CPO-I). Density functional calculations (DFTs) were performed to determine the optimized geometries and electronic structure of each of these three species. Their electronic spectra were also calculated at the DFT optimized geometries, using the INDO/S/CI method. The results of these studies led to the following conclusions: (1) the presence of the nearby Asp in a typical peroxidase does indeed decrease the energy required to deprotonate the imidazole making the two forms essentially degenerate, (2) neither the state of protonation of the imidazole ligand nor the change in axial ligand from an imidazole in typical peroxidases such as HRP to a mercaptide in CPO significantly alters the characteristics of the lowest energy spin state or the electronic structure of compound I in a way that can obviously affect function, (3) both the Im-and ImH forms of the peroxidase compound I (HRP-I) lead to the same dramatic reduction in intensity relative to the ferric resting form observed experimentally. However, only in the ImH form of HRP-I does the calculated relative shift of one component of the Soret bands relative to CPO-I agree with that observed in the transient spectra of HRP-I compared to CPO-I. These results taken together strongly indicate that factors other than the nature of the proximal axial ligand are the main determinants of function.

Light forces the pace: optical manipulation for biophotonics

The biomedical sciences have benefited immensely from photonics technologies in the last 50 years. This includes the application of minute forces that enable the trapping and manipulation of cells and single molecules. In terms of the area of biophotonics, optical manipulation has made a seminal contribution to our understanding of the dynamics of single molecules and the microrheology of cells. Here we present a review of optical manipulation, emphasizing its impact on the areas of single-molecule studies and single-cell biology, and indicating some of the key experiments in the fields.

Statistical determination of the step size of molecular motors

Molecular motors are enzymatic proteins that couple the consumption of chemical energy to mechanical displacement. In order to elucidate the translocation mechanisms of these enzymes, it is of fundamental importance to measure the physical step size. The step size can, in certain instances, be directly measured with single-molecule techniques; however, in the majority of cases individual steps are masked by noise. The step size can nevertheless be obtained from noisy single-molecule records through statistical methods. This analysis is analogous to determining the charge of the electron from current shot noise. We review methods for obtaining the step size based on analysing, in both the time and frequency domains, the variance in position from noisy single-molecule records of motor displacement. Additionally, we demonstrate how similar methods may be applied to measure the step size in bulk kinetic experiments.

The intrinsic axial ligand effect on propene oxidation by horseradish peroxidase versus cytochrome P450 enzymes

The axial ligand effect on reactivity of heme enzymes is explored by means of density functional theoretical calculations of the oxidation reactions of propene by a model compound I species of horseradish peroxidase (HRP). The results are assessed vis-a-vis those of cytochrome P450 compound I. It is shown that the two enzymatic species perform C=C epoxidation and C-H hydroxylation in a multistate reactivity scenario with Fe(III) and Fe(IV) electromeric situations and two different spin states, doublet and quartet. However, while the HRP species preferentially keeps the iron in a low oxidation state (Fe(III)), the cytochrome P450 species prefers the higher oxidation state (Fe(IV)). It is found that HRP compound I has somewhat lower barriers than those obtained by the cytochrome P450 species. Furthermore, in agreement with experimental observations and studies on model systems, HRP prefers C=C epoxidation, whereas cytochrome P450 prefers C-H hydroxylation. Thus, had the compound I species of HRP been by itself, it would have been an epoxidizing agent, and at least as reactive as cytochrome P450. In the enzyme, HRP is much less reactive than cytochrome P450, presumably because HRP reactivity is limited by the access of the substrate to compound I.

QM/MM studies of the electronic structure of the compound I intermediate in cytochrome c peroxidase and ascorbate peroxidase

Cytochrome c peroxidase (CcP) and ascorbate peroxidase (APX) both involve reactive haem oxoferryl intermediates known as 'compound I' species. These two enzymes also have a very similar structure, especially in the vicinity of the haem group. Despite this similarity, the electronic structure of compound I in the two enzymes is known to be very different. Compound I intermediates have three unpaired electrons, two of which are always situated on the Fe-O core, whilst the third is located in a porphyrin orbital in APX and many other compound I species. In CcP, however, this third unpaired electron is positioned on a tryptophan residue lying close to the haem ring. The same residue is present in the same position in APX, yet it is not oxidized in that case. We report QM/MM calculations, using accurate B3LYP density functional theory for the QM region, on the active intermediate for both enzymes. We reproduce the observed difference in electronic structure, and show that it arises as a result of subtle electrostatic effects which affect the ionization potential of both the tryptophan and porphyrin groups. The computed structures of both enzymes do not involve deprotonation of the tryptophan group, or protonation of the oxoferryl oxygen.

On the role of the axial ligand in heme proteins: a theoretical study

We present a systematic investigation of how the axial ligand in heme proteins influences the geometry, electronic structure, and spin states of the active site, and the energies of the reaction cycles. Using the density functional B3LYP method and medium-sized basis sets, we have compared models with His, His+Asp, Cys, Tyr, and Tyr+Arg as found in myoglobin and hemoglobin, peroxidases, cytochrome P450, and heme catalases, respectively. We have studied 12 reactants and intermediates of the reaction cycles of these enzymes, including complexes with H(2)O, OH(-), O(2-), CH(3)OH, O(2), H(2)O(2), and HO(2)(-) in various formal oxidation states of the iron ion (II to V). The results show that His gives ~0.6 V higher reduction potentials than the other ligands. In particular, it is harder to reduce and protonate the O(2) complex with His than with the other ligands, in accordance with the O(2) carrier function of globins and the oxidative chemistry of the other proteins. For most properties, the trend Cys<Tyr<Tyr+Arg<His+Asp<His is found, reflecting the donor capacity of the various ligands. Thus, it is easier to reduce compound I with a His+Asp ligand than with a Cys ligand, in accordance with the one-electron chemistry of peroxidases and the hydroxylation reactions of cytochromes P450. However, the Tyr complexes have an unusually low affinity for all neutral ligands, giving them a slightly enhanced driving force in the oxidation of H(2)O(2) by compound I.

Kinesin moves by an asymmetric hand-over-hand mechanism

Kinesin is a double-headed motor protein that moves along microtubules in 8-nanometer steps. Two broad classes of model have been invoked to explain kinesin movement: hand-over-hand and inchworm. In hand-over-hand models, the heads exchange leading and trailing roles with every step, whereas no such exchange is postulated for inchworm models, where one head always leads. By measuring the stepwise motion of individual enzymes, we find that some kinesin molecules exhibit a marked alternation in the dwell times between sequential steps, causing these motors to "limp" along the microtubule. Limping implies that kinesin molecules strictly alternate between two different conformations as they step, indicative of an asymmetric, hand-over-hand mechanism.

Crystal structure of manganese catalase from Lactobacillus plantarum

BACKGROUND

Catalases are important antioxidant metalloenzymes that catalyze disproportionation of hydrogen peroxide, forming dioxygen and water. Two families of catalases are known, one having a heme cofactor, and the other, a structurally distinct family containing nonheme manganese. We have solved the structure of the mesophilic manganese catalase from Lactobacillus plantarum and its azide-inhibited complex.

RESULTS

The crystal structure of the native enzyme has been solved at 1.8 A resolution by molecular replacement, and the azide complex of the native protein has been solved at 1.4 A resolution. The hexameric structure of the holoenzyme is stabilized by extensive intersubunit contacts, including a beta zipper and a structural calcium ion crosslinking neighboring subunits. Each subunit contains a dimanganese active site, accessed by a single substrate channel lined by charged residues. The manganese ions are linked by a mu1,3-bridging glutamate carboxylate and two mu-bridging solvent oxygens that electronically couple the metal centers. The active site region includes two residues (Arg147 and Glu178) that appear to be unique to the Lactobacillus plantarum catalase.

CONCLUSIONS

A comparison of L. plantarum and T. thermophilus catalase structures reveals the existence of two distinct structural classes, differing in monomer design and the organization of their active sites, within the manganese catalase family. These differences have important implications for catalysis and may reflect distinct biological functions for the two enzymes, with the L. plantarum enzyme serving as a catalase, while the T. thermophilus enzyme may function as a catalase/peroxidase.

Calculation of the electronic structure and spectra of model cytochrome P450 compound I

The electronic structure and spectra of the oxyferryl (Fe=O) compound I P450 heme species, the transient putative active intermediate of cytochrome P450s, have been calculated employing a full protoporphyrin IX heme model representation. The principal aim of this work was to compare the computed spectra of this species with the observed transient spectra attributed to it. Computations were made using both nonlocal density functional theory (DFT) and semiempirical INDO/CI methods to characterize the electronic structure of the compound I P450 species. Both methods resulted in a similar antiferromagnetic doublet as the ground state with a ferromagnetic quartet excited state partner, slightly higher in energy. The INDO/ROHF/CI semiempirical method was used to calculate the spectrum of the protoporphyrin IX P450 compound I heme species in its lowest energy antiferromagnetic doublet state at the DFT optimized geometry. As a reference, the spectrum of the ferric resting form of the protoporphyrin IX P450 heme species was also calculated. The computed shifts in the Soret and Q bands of compound I relative to the resting state were both in good agreement with the corresponding experimentally observed shifts in the transient spectra of cytochrome P450cam (Biochem. Biophys. Res. Commun. 201 (1994) 1464) and chloroperoxidase (Biochem. Biophys. Res. Commun. 94 (1980) 1123) both ascribed to their common compound I heme site. This consistency provides additional, independent support for the assignment of compound I as the origin of the reported observed transient spectra.

Understanding the structure and function of catalases: clues from molecular evolution and in vitro mutagenesis

This review gives an overview about the structural organisation of different evolutionary lines of all enzymes capable of efficient dismutation of hydrogen peroxide. Major potential applications in biotechnology and clinical medicine justify further investigations. According to structural and functional similarities catalases can be divided in three subgroups. Typical catalases are homotetrameric haem proteins. The three-dimensional structure of six representatives has been resolved to atomic resolution. The central core of each subunit reveals a characteristic "catalase fold", extremely well conserved among this group. In the native tetramer structure pairs of subunits tightly interact via exchange of their N-terminal arms. This pseudo-knot structures implies a highly ordered assembly pathway. A minor subgroup ("large catalases") possesses an extra flavodoxin-like C-terminal domain. A > or = 25 A long channel leads from the enzyme surface to the deeply buried active site. It enables rapid and selective diffusion of the substrates to the active center. In several catalases NADPH is tightly bound close to the surface. This cofactor may prevent and reverse the formation of compound II, an inactive reaction intermediate. Bifunctional catalase-peroxidase are haem proteins which probably arose via gene duplication of an ancestral peroxidase gene. No detailed structural information is currently available. Even less is know about manganese catalases. Their di-manganese reaction centers may be evolutionary.

Theoretical study of the electrostatic and steric effects on the spectroscopic characteristics of the metal-ligand unit of heme proteins. 2. C-O vibrational frequencies, 17O isotropic chemical shifts, and nuclear quadrupole coupling constants

The quantum chemical calculations, vibronic theory of activation, and London-Pople approach are used to study the dependence of the C-O vibrational frequency, 17O isotropic chemical shift, and nuclear quadrupole coupling constant on the distortion of the porphyrin ring and geometry of the CO coordination, changes in the iron-carbon and iron-imidazole distances, magnitude of the iron displacement out of the porphyrin plane, and presence of the charged groups in the heme environment. It is shown that only the electrostatic interactions can cause the variation of all these parameters experimentally observed in different heme proteins, and the heme distortions could modulate this variation. The correlations between the theoretically calculated parameters are shown to be close to the experimentally observed ones. The study of the effect of the electric field of the distal histidine shows that the presence of the four C-O vibrational bands in the infrared absorption spectra of the carbon monoxide complexes of different myoglobins and hemoglobins can be caused by the different orientations of the different tautomeric forms of the distal histidine. The dependence of the 17O isotropic chemical shift and nuclear quadrupole coupling constant on pH and the distal histidine substitution can be also explained from the same point of view.

Resonance Raman spectra of native and mesoheme-reconstituted horseradish peroxidase and their catalytic intermediates

Resonance Raman studies of native and mesoheme-reconstituted horseradish peroxidase and their catalytic intermediates, known as Compounds I and II, have been conducted using both near UV ( approximately 350 nm) and visible (406.7 nm) excitation. Careful power studies indicate that the authentic Compound I spectra are obtainable using near UV excitation, but that use of visible excitation results in contamination of the Compound I spectrum with the spectrum of a Compound II-like photoproduct. Using H218O2, the nu(Fe=O) stretching modes for both systems are unambiguously identified, for the first time, at approximately 790 cm-1. The authentic Compound I spectra are indicative of an 2A1u-like ground state for both the native and the mesoheme-reconstituted proteins. Finally, the possible biological implications of such information are briefly discussed.

Theoretical study of the distal-side steric and electrostatic effects on the vibrational characteristics of the FeCO unit of the carbonylheme proteins and their models

The vibronic theory of activation and quantum chemical intermediate neglect of differential overlap (INDO) calculations are used to study the activation of carbon monoxide (change of the C-O bond index and force field constant) by the imidazole complex with heme in dependence on the distortion of the porphyrin ring, geometry of the CO coordination, iron-carbon and iron-imidazole distances, iron displacement out of the porphyrin plane, and presence of the charged groups in the heme environment. It is shown that the main contribution to the CO activation stems from the change in the sigma donation from the 5 sigma CO orbital to iron, and back-bonding from the iron to the 2 pi orbital of CO. It follows from the results that none of the studied distortions can explain, by itself, the wide variation of the C-O vibrational frequency in the experimentally studied model compounds and heme proteins. To study the dependence of the properties of the FeCO unit on the presence of charged groups in the heme environment, the latter are simulated by the homogeneous electric field and point charges of different magnitude and location. The results show that charged groups can strongly affect the strength of the C-O bond and its vibrational frequency. It is found that the charges located on the distal side of the heme plane can affect the Fe-C and C-O bond indexes (and, consequently, the Fe-C and C-O vibrational frequencies), both in the same and in opposite directions, depending on their position. The theoretical results allow us to understand the peculiarities of the effect of charged groups on the properties of the FeCO unit both in heme proteins and in their model compounds.