The Poulos-Kraut mechanism of Compound I formation in horseradish peroxidase: a QM/MM study

QM/MM calculations are used to elucidate the Poulos-Kraut (Poulos, T. L.; Kraut, J. J. Biol. Chem. 1980, 255, 8199-8205) mechanism of O-O bond activation and Compound I (Cpd I) formation in HRP, in conditions corresponding to neutral to basic pH. Attempts to generate Compound I directly from the Fe(H2O2) complex by migrating the proton from the proximal oxygen to the distal one (1,2- proton shift) result in high barriers. The lowest energy mechanism was found to involve initial deprotonation of ferric hydrogen peroxide complex (involving spin crossover from the quartet to the doublet state) by His42 to form ferric-hydroperoxide (Cpd 0). Subsequently, the distal OH group of Cpd 0 is pulled by Arg38 and reprotonated by His42(H+) to form Cpd I and a water molecule that bridges the two residues. The structures of the intermediate and the transition state reveal the manner by which the Arg-His residues promote cooperatively the electronic reorganization that is required to attend the heterolytic O-O cleavage.     

Theoretical analysis of the structural and electronic properties of metalloporphyrin pi-cation radicals

A method for analyzing the A(1u)/A(2u) contents of metalloporphyrin pi-cation radicals is developed and applied to a series of unsubstituted planar metalloporphines (MPs) (M=Cr, Mn, Fe, Co, Ni, Cu, and Zn). The structures and electronic properties of the MPs and their cation radicals were calculated by density functional theory (DFT) and subsequently analyzed. It was found that the MPs with small core sizes have a tendency to form A(1u)-type radicals, while the MPs with large core size have a preference for an A(2u)-type. Neither of these pure-state species, however, is stable under the D(4)(h) symmetry, and both radical cation types are subject to pseudo-Jahn-Teller (pJT) distortion. The pJT distortion leads to structures with lower symmetry and states that have mixed character with respect to the A(1u) and A(2u) components. The degree of mixing could be estimated by employing orbital projection technique or a complementary spin density decomposition. Both techniques produce very similar results, pointing out that the frontier orbital, which becomes empty upon electron removal, plays a critical role in determining electronic properties.     

Orbitals and the Interpretation of Photoelectron Spectroscopy and (e,2e) Ionization Experiments

Electron momentum spectroscopy, scanning tunneling microscopy, and photoelectron spectroscopy provide unique information about electronic structure, but their interpretation has been controversial. This essay discusses a framework for interpretation. Although this interpretation is not new, we believe it is important to present this framework in light of recent publications. The key point is that these experiments provide information about how the electron distribution changes upon ionization, not how electrons behave in the pre‐ionized state. Therefore, these experiments do not lead to a “selection of the correct orbitals” in chemistry and do not overturn the well‐known conclusion that both delocalized molecular orbitals and localized molecular orbitals are useful for interpreting chemical structure and dynamics. The two types of orbitals can produce identical total molecular electron densities and therefore molecular properties. Different types of orbitals are useful for different purposes.     

Some Useful Integral Representations for Information-Theoretic Analyses

This work is an extension of our earlier article, where a well-known integral representation of the logarithmic function was explored and was accompanied with demonstrations of its usefulness in obtaining compact, easily-calculable, exact formulas for quantities that involve expectations of the logarithm of a positive random variable. Here, in the same spirit, we derive an exact integral representation (in one or two dimensions) of the moment of a nonnegative random variable, or the sum of such independent random variables, where the moment order is a general positive non-integer real (also known as fractional moments). The proposed formula is applied to a variety of examples with an information-theoretic motivation, and it is shown how it facilitates their numerical evaluations. In particular, when applied to the calculation of a moment of the sum of a large number, n, of nonnegative random variables, it is clear that integration over one or two dimensions, as suggested by our proposed integral representation, is significantly easier than the alternative of integrating over n dimensions, as needed in the direct calculation of the desired moment.     

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.     

Mechanistic Insight into Peptidyl-Cysteine Oxidation by the Copper-Dependent Formylglycine-Generating Enzyme

The copper-dependent formylglycine-generating enzyme (FGE) catalyzes the oxygen-dependent oxidation of specific peptidyl-cysteine residues to formylglycine. Our QM/MM calculations provide a very likely mechanism for this transformation. The reaction starts with dioxygen binding to the tris-thiolate Cu^I center to form a triplet Cu^II-superoxide complex. The rate-determining hydrogen atom abstraction involves a triplet-singlet crossing to form a Cu^II−OOH species that couples with the substrate radical, leading to a Cu^I-alkylperoxo intermediate. This is accompanied by proton transfer from the hydroperoxide to the S atom of the substrate via a nearby water molecule. The subsequent O−O bond cleavage is coupled with the C−S bond breaking that generates the formylglycine and a Cu^II-oxyl complex. Moreover, our results suggest that the aldehyde oxygen of the final product originates from O₂, which will be useful for future experimental work.     

Redox and acid-base coupling in ultrathin polyelectrolyte films

A single layer of poly(allylamine) with a covalently attached osmium pyridine-bipyridine complex adsorbed onto a Au surface modified by mercaptopropanesulfonate has been studied theoretically with a molecular approach and experimentally by cyclic voltammetry. These investigations have been carried out at different pHs and ionic strengths of the electrolyte solution in contact with the redox polyelectrolyte modified electrode. The theory predicts strong coupling between the acid-base and redox equilibria, particularly for low ionic strength, pH close to the pKa, and high concentration of redox sites. The coupling leads to a decrease in the peak potential at pH values above the apparent pKa of the weak polyelectrolyte, in good agreement with the experimental pH dependence at 4 mM NaNO3. Theoretical calculations suggest that the inflection point in the peak position versus pH curves can be used to estimate the apparent pKa of the amino groups in the polymer. Comparison of the apparent pKa for PAH-Os in the film with that of poly(allylamine) reported in the literature shows that the underlying charged thiol strongly influences charge regulation in the film. A systematic study of the film thickness and the degree of protonation in sulfonate and amino groups for solutions of different pH and ionic strength shows the coupling between the different interactions. It is found that the variation of the film properties has a non-monotonic dependence on bulk pH and salt concentration. For example, the film thickness shows a maximum with electrolyte ionic strength, whose origin is attributed to the balance between electrostatic amino-amino repulsions and amino-sulfonate attractions.     

Quantum mechanical/molecular mechanical study of mechanisms of heme degradation by the enzyme heme oxygenase: the strategic function of the water cluster

Heme degradation by heme oxygenase (HO) enzymes is important in maintaining iron homeostasis and prevention of oxidative stress, etc. In response to mechanistic uncertainties, we performed quantum mechanical/molecular mechanical investigations of the heme hydroxylation by HO, in the native route and with the oxygen surrogate donor H2O2. It is demonstrated that H2O2 cannot be deprotonated to yield Fe(III)OOH, and hence the surrogate reaction starts from the FeHOOH complex. The calculations show that, when starting from either Fe(III)OOH or Fe(III)HOOH, the fully concerted mechanism involving O-O bond breakage and O-C(meso) bond formation is highly disfavored. The low-energy mechanism involves a nonsynchronous, effectively concerted pathway, in which the active species undergoes first O-O bond homolysis followed by a barrier-free (small with Fe(III)HOOH) hydroxyl radical attack on the meso position of the porphyrin. During the reaction of Fe(III)HOOH, formation of the Por+*FeIV=O species, compound I, competes with heme hydroxylation, thereby reducing the efficiency of the surrogate route. All these conclusions are in accord with experimental findings (Chu, G. C.; Katakura, K.; Zhang, X.; Yoshida, T.; Ikeda-Saito, M. J. Biol. Chem. 1999, 274, 21319). The study highlights the role of the water cluster in the distal pocket in creating "function" for the enzyme; this cluster affects the O-O cleavage and the O-Cmeso formation, but more so it is responsible for the orientation of the hydroxyl radical and for the observed alpha-meso regioselectivity of hydroxylation (Ortiz de Montellano, P. R. Acc. Chem. Res. 1998, 31, 543). Differences/similarities with P450 and HRP are discussed.     

Electron transfer activation of chromopyrrolic acid by cytochrome p450 en route to the formation of an antitumor indolocarbazole derivative: theory supports experiment

QM/MM calculations support experiment and show that StaP is a P450 that functions like a peroxidase: its active species is the one-electron-reduced Cpd II species with a radical on CPA, by analogy to cytochrome c peroxidase (CcP), and its reaction with the substrate proceeds by overall proton-coupled electron transfer (PCET), in analogy to the corresponding mechanism in horseradish peroxidase (HRP). The electron transfer is enabled by His250, the presence of carboxylate groups in CPA, and by the H-bonding network that tunes the energetic of the process. Theory supports experiment but reveals some novel aspects of this unusual P450.     

A valence bond modeling of trends in hydrogen abstraction barriers and transition states of hydroxylation reactions catalyzed by cytochrome P450 enzymes

The paper outlines the fundamental factors that govern the mechanisms of alkane hydroxylation by cytochrome P450 and the corresponding barrier heights during the hydrogen abstraction and radical rebound steps of the process. This is done by a combination of density functional theory calculations for 11 alkanes and valence bond (VB) modeling of the results. The energy profiles and transition states for the various steps are reconstructed using VB diagrams (Shaik, S. S. J. Am. Chem. Soc. 1981, 103, 3692-3701. Shaik, S.; Shurki, A. Angew. Chem. Int. Ed. 1999, 38, 586-625.) and the DFT barriers are reproduced by the VB model from raw data based on C-H bond energies. The model explains a variety of other features of P450 hydroxylations: (a) the nature of the polar effect during hydrogen abstraction, (b) the difference between the activation mechanisms leading to the Fe(IV) vs the Fe(III) electromers, (c) the difference between the gas phase and the enzymatic reaction, and (d) the dependence of the rebound barrier on the spin state. The VB mechanism shows that the active species of the enzyme performs a complex reaction that involves multiple bond making and breakage mechanisms by utilizing an intermediate VB structure that cuts through the high barrier of the principal transformation between reactants and products, thereby mediating the process at a low energy cost. The correlations derived in this paper create order and organize the data for a process of a complex and important enzyme. This treatment can be generalized to the reactivity patterns of nonheme systems and synthetic iron-oxo porphyrin reagents.