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     

Mechanistic pathways of the hydroxyl radical reactions of quinoline. 2. Computational analysis of hydroxyl radical attack at C atoms

Density functional theory (DFT) calculations are employed to compare the mechanism of the *OH attacks at all carbon atoms in quinoline. The computational analysis of the energy surface for the reaction of *OH with quinoline reveals that the formation of OH adducts proceeds through exothermic formation of pi-complexes/H-bonded complexes. The gas-phase reactions have activation energies ranging from <1.3 kcal/mol for the attack at positions C3 through C8 to 8.6 kcal/mol for the attack at the C2 position. Solvation, as described by the CPCM cavity model, lowers these activation barriers so that the attack at all carbon atoms except C2 is effectively barrierless. The *OH attack at C2 in solution is significantly different than at all other quinoline positions because it involves the only transition structure with energy higher than that of the starting materials and with an energetic barrier of 5.1 kcal/mol. The specific solvation approach also corroborates this finding because the attack at C2 was shown to have an energy barrier of 2.3 kcal/mol compared to the barrierless attack at C5. These results are in agreement with our recent experimental studies but differ from literature reports on the degradation of quinoline using the photo-Fenton reaction.     

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.     

A new understanding on how heme metabolism occurs in heme oxygenase: water-assisted oxo mechanism

Heme metabolism by heme oxygenase (HO) is investigated with quantum mechanical/molecular mechanical (QM/MM) calculations. A mechanism assisted by water is proposed: (1) an iron-oxo species and a water molecule are generated by the heterolytic cleavage of the O-O bond of an iron-hydroperoxo species in a similar way to P450-mediated reactions, (2) a hydrogen atom abstraction by the iron-oxo species from the generated water molecule and the C-O bond formation between the water molecule and the α-meso carbon take place simultaneously. The water molecule is hydrogen-bonded to the oxo ligand and to the water cluster in the active site of HO. The water cluster can control the position of the generated water molecule to ensure the regioselective oxidation of heme at the α-meso position, at the same time, can facilitate the oxidation by stabilizing a positive charge on the water molecule in the transition state. A key difference between HO and P450 is observed in the structure of the active site; Thr252 in P450 blocks the access of the water molecule to the α-meso position, and can thus suppress the undesired heme oxidation for P450.     

Role of the somersault rearrangement in the oxidation step for flavin monooxygenases (FMO). A comparison between FMO and conventional xenobiotic oxidation with hydroperoxides

Model quantum mechanical calculations presented for C-4a-flavin hydroperoxide (FlHOOH) at the B3LYP/6-311+G(d,p) level suggest a new mechanism for flavoprotein monooxygenase (FMO) oxidation involving a concerted homolytic O-O bond cleavage in concert with hydroxyl radical transfer from the flavin hydroperoxide rather than an S(N)2-like displacement by the substrate on the C-4a-hydroperoxide OOH group. Homolytic O-O bond cleavage in a somersault-like rearrangement of hydroperoxide C-4a-flavinhydroperoxide (1) (FLHO-OH → FLHO···HO) produces an internally hydrogen-bonded HO(?) radical intermediate with a classical activation barrier of 27.0 kcal/mol. Model hydroperoxide 1 is used to describe the transition state for the key oxidation step in the paradigm aromatic hydroxylase, p-hydroxybenzoate hydroxylase (PHBH). A comparison of the electron distribution in the transition structures for the PHBH hydroxylation of p-hydroxybenzoic acid (ΔE(?) = 23.0 kcal/mol) with that of oxidation of trimethylamine (ΔE(?) = 22.3 kcal/mol) and dimethyl sulfide (ΔE? = 14.1 kcal/mol) also suggests a mechanism involving a somersault mechanism in concert with transfer of an HO(?) radical to the nucleophilic heteroatom center with a hydrogen transfer back to the FLH-O residue after the barrier is crossed to produce the final product, FLH-OH. In each case the hydroxylation barrier was less than that of the O-O rearrangement barrier in the absence of a substrate supporting an overall concerted process. All three transition structures bear a resemblance to the TS for the comparable hydroxylation of isobutane (ΔE(?) = 29.2 kcal/mol) and for simple Fenton oxidation by aqueous iron(III) hydroperoxides. To our surprise the oxidation of N- and S-nucleophiles with conventional oxidants such as alkyl hydroperoxides and peracids also proceeds by HO(?) radical transfer in a manner quite similar to that for tricyclic hydroperoxide 1. Stabilization of the developing oxyradical produced by somersault rearrangement for concerted enzymatic oxidation with tricyclic hydroperoxide 1 results in a reduced overall activation barrier.     

Computer-assisted study on the reaction between pyruvate and ylide in the pathway leading to lactyl–ThDP

In this study the formation of the lactyl-thiamin diphosphate intermediate (L-ThDP) is addressed using density functional theory calculations at X3LYP/6-31++G(d,p) level of theory. The study includes potential energy surface scans, transition state search, and intrinsic reaction coordinate calculations. Reactivity is analyzed in terms of Fukui functions. The results allow to conclude that the reaction leading to the formation of L-ThDP occurs via a concerted mechanism, and during the nucleophilic attack on the pyruvate molecule, the ylide is in its AP form. The calculated activation barrier for the reaction is 19.2?kcal/mol, in agreement with the experimental reported value.     

Theoretical study of oxidation of cyclohexane diol to adipic anhydride by [Ru(IV)(O)(tpa)(H2O)]2+ complex (tpa ═ tris(2-pyridylmethyl)amine)

The catalytic conversion of 1,2-cyclohexanediol to adipic anhydride by Ru(IV)O(tpa) (tpa ═ tris(2-pyridylmethyl)amine) is discussed using density functional theory calculations. The whole reaction is divided into three steps: (1) formation of α-hydroxy cyclohexanone by dehydrogenation of cyclohexanediol, (2) formation of 1,2-cyclohexanedione by dehydrogenation of α-hydroxy cyclohexanone, and (3) formation of adipic anhydride by oxygenation of cyclohexanedione. In each step the two-electron oxidation is performed by Ru(IV)O(tpa) active species, which is reduced to bis-aqua Ru(II)(tpa) complex. The Ru(II) complex is reactivated using Ce(IV) and water as an oxygen source. There are two different pathways of the first two steps of the conversion depending on whether the direct H-atom abstraction occurs on a C-H bond or on its adjacent oxygen O-H. In the first step, the C-H (O-H) bond dissociation occurs in TS1 (TS2-1) with an activation barrier of 21.4 (21.6) kcal/mol, which is followed by abstraction of another hydrogen with the spin transition in both pathways. The second process also bifurcates into two reaction pathways. TS3 (TS4-1) is leading to dissociation of the C-H (O-H) bond, and the activation barrier of TS3 (TS4-1) is 20.2 (20.7) kcal/mol. In the third step, oxo ligand attack on the carbonyl carbon and hydrogen migration from the water ligand occur via TS5 with an activation barrier of 17.4 kcal/mol leading to a stable tetrahedral intermediate in a triplet state. However, the slightly higher energy singlet state of this tetrahedral intermediate is unstable; therefore, a spin crossover spontaneously transforms the tetrahedral intermediate into a dione complex by a hydrogen rebound and a C-C bond cleavage. Kinetic isotope effects (k(H)/k(D)) for the electronic processes of the C-H bond dissociations calculated to be 4.9-7.4 at 300 K are in good agreement with experiment values of 2.8-9.0.     

A theoretical examination of the acid-catalyzed and noncatalyzed ring-opening reaction of an oxirane by nucleophilic addition of acetate. Implications to epoxide hydrolases

Ab initio and density functional calculations have been performed to gain a better understanding of the epoxide ring-opening reaction catalyzed by epoxide hydrolase. The S(N)2 reaction of acetate with 1S,2S-trans-2-methylstyrene oxide to provide the corresponding diol acetate ester was studied with and without general-acid catalysis. MP2 and DFT (B3LYP) calculations predict, for the noncatalyzed reaction, a central barrier of approximately 20-21 kcal/mol separating the reactants from products depending on which carbon center in the epoxide is undergoing attack. From these gas-phase reactions the immediate alkoxide products are not energetically far below their associated transition states such that the reaction is predicted to be endothermic. Inclusion of aqueous solvation effects via a polarizable continuum model predicts the activation barrier to increase by almost 10 kcal/mol due to the solvation of the acetate ion nucleophile. The activation barrier for the epoxide ring-opening reaction is reduced to approximately 10 kcal/mol when phenol, as the general-acid catalyst, is included in the gas-phase calculations. This is due to the immediate product being the neutral ester rather than the corresponding alkoxide. The transition state in the general-acid-catalyzed reaction is earlier than that for the noncatalyzed reaction and the reaction is highly exothermic. Molecular mechanics calculations of 1S,2S-trans-2-methylstyrene oxide in the active site of murine epoxide hydrolase show two possible binding conformations. Both conformers have the epoxide oxygen forming hydrogen bonds with the acidic hydrogens of the catalytic tyrosines (Tyr381 and Tyr465). These two conformations likely lead to different products since the nucleophile (Asp333-CO(2)(-)) is positioned to react with either carbon center in the epoxide.     

A "push-pull" mechanism for heterolytic o-o bond cleavage in hydroperoxo manganese porphyrins

A water-soluble manganese porphyrin, 5,10,15,20-tetrakis-(1,3-dimethylimidazolium-2-yl)porphyrinatomanganese(III) (Mn(III)TDMImP) is shown to react with H(2)O(2) to generate a relatively stable dioxomanganese(V) porphyrin complex (a compound I analog). Stopped-flow kinetic studies revealed Michaelis Menton-type saturation kinetics for H(2)O(2). The visible spectrum of a compound 0 type intermediate, assigned as Mn(III)(OH)(OOH)TDMImP, can be directly observed under saturating H(2)O(2) conditions (Soret band at 428 nm and Q bands at 545 and 578 nm). The rate-determining O-O heterolysis step was found to have a very small activation enthalpy (ΔH(≠) = 4.2 ± 0.2 kcal mol(-1)) and a large, negative activation entropy (ΔS(≠) = -36 ± 1 cal mol(-1) K(-1)). The O-O bond cleavage reaction was pH independent at 8.8 < pH < 10.4 with a first-order rate constant of 66 ± 12 s(-1). These observations indicate that the O-O bond in Mn(III)(OH)(OOH)TDMImP is cleaved via a concerted "push-pull" mechanism. In the transition state, the axial (proximal) (-)OH is partially deprotonated ("push"), while the terminal oxygen in (-)OOH is partially protonated ("pull") as a water molecule is released to the medium. This mechanism is reminiscent of O-O bond cleavage in heme enzymes, such as peroxidases and cytochrome P450, and similar to the fast, reversible O-Br bond breaking and forming reaction mediated by similar manganese porphyrins. The small enthalpy of activation suggests that this O-O bond cleavage could also be made reversible.     

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.     

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.     

Can hydridic-to-protonic hydrogen bonds catalyze hydride transfers in biological systems?

Catalysis of hydride transfer by hydridic-to-protonic hydrogen (HHH) bonding in α-hydroxy carbonyl isomerization reactions was examined computationally in the lithium salts of 7-substituted endo-3-hydroxybicyclo[2.2.1]hept-5-en-2-ones. The barrier for intramolecular hydride transfer in the parent system was calculated to be 17.2 kcal/mol. Traditional proton donors, such as OH and NH(3)(+), stabilized the metal cation-bridged transition state by 1.4 and 3.3 kcal/mol, respectively. Moreover, among the conformers of the OH systems, the one in which the proton donor is able to interact with the migrating hydride (H(m)) has an activation barrier lower by 1.3 and 1.7 kcal/mol than the other possible OH conformers. By contrast, the presence of an electronegative group such as F, which disfavors the migration electronically by opposing development of hydridic charge, destabilizes the hydride migration by 1.5 kcal/mol relative to the epimeric exo system. In both ground and transition states the H(m)···H distance decreased with increasing acidity of the proton donor, reaching a minimum of 1.58 ? at the transition state for NH(3)(+). Both Mulliken and NPA charges show enhancement of negative character of the migrating hydride in the cases in which HHH bonding is possible.     

Theoretical mechanism study of UF6 hydrolysis in the gas phase

The mechanism of the gas-phase reaction UF 6 + H 2O --> UOF 4 + 2HF is explored using relativistic density functional theory calculations. Initially, H 2O coordinates with UF 6 to form a 1:1 complex UF 6.H 2O. Over an activation energy barrier of about 19 kcal/mol, H 2O transfers a H atom to a nearby ligand F, resulting in UF 5OH + HF. The eliminated HF or another H 2O molecule may form a hydrogen bond with UF 5OH. Starting from UF 5OH, the second HF elimination results in UOF 4. If UF 5OH is in the isolated form, UF 5OH --> UOF 4 + HF takes place over a barrier of 24 kcal/mol. If UF 5OH is hydrogen-bonded with H 2O or HF, the conversion barrier is less than 10 kcal/mol. Once formed, the unstable UOF 4 tends to associate with additional ligands and hydrogen-bonding donors. The calculated binding energies indicate the significance of such interactions, which may have profound impact on further HF eliminating reactions. The IR spectra features can be used to indicate the formation and interaction type of the intermediates and products.     

Does Hydrogen‐Bonding Donation to Manganese(IV)–Oxo and Iron(IV)–Oxo Oxidants Affect the Oxygen‐Atom Transfer Ability? A Computational Study

Iron(IV)–oxo intermediates are involved in oxidations catalyzed by heme and nonheme iron enzymes, including the cytochromes P450. At the distal site of the heme in P450 Compound I (Fe^IV–oxo bound to porphyrin radical), the oxo group is involved in several hydrogen‐bonding interactions with the protein, but their role in catalysis is currently unknown. In this work, we investigate the effects of hydrogen bonding on the reactivity of high‐valent metal–oxo moiety in a nonheme iron biomimetic model complex with trigonal bipyramidal symmetry that has three hydrogen‐bond donors directed toward a metal(IV)–oxo group. We show these interactions lower the oxidative power of the oxidant in reactions with dehydroanthracene and cyclohexadiene dramatically as they decrease the strength of the O? H bond (BDE_OH) in the resulting metal(III)–hydroxo complex. Furthermore, the distal hydrogen‐bonding effects cause stereochemical repulsions with the approaching substrate and force a sideways attack rather than a more favorable attack from the top. The calculations, therefore, give important new insights into distal hydrogen bonding, and show that in biomimetic, and, by extension, enzymatic systems, the hydrogen bond may be important for proton‐relay mechanisms involved in the formation of the metal–oxo intermediates, but the enzyme pays the price for this by reduced hydrogen atom abstraction ability of the intermediate. Indeed, in nonheme iron enzymes, where no proton relay takes place, there generally is no donating hydrogen bond to the iron(IV)–oxo moiety.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

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.     

Theory favors a stepwise mechanism of porphyrin degradation by a ferric hydroperoxide model of the active species of heme oxygenase

The report uses density functional theory to address the mechanism of heme degradation by the enzyme heme oxygenase (HO) using a model ferric hydroperoxide complex. HO is known to trap heme molecules and degrade them to maintain iron homeostasis in the biosystem. The degradation is initiated by complexation of the heme, then formation of the iron-hydroperoxo species, which subsequently oxidizes the meso position of the porphyrin by hydroxylation, thereby enabling eventually the cleavage of the porphyrin ring. Kinetic isotope effect studies indicate that the mechanism is assisted by general acid catalysis, via a chain of water molecules, and that all the events occur in concert. However, previous theoretical treatments indicated that the concerted mechanism has a high barrier, much higher than an alternative mechanism that is initiated by O-O bond homolysis of iron-hydroperoxide. The present contribution studies the stepwise and concerted acid-catalyzed mechanisms using H(3)O(+)(H(2)O)(n)(), n = 0-2. The effect of the acid strength is tested using the H(4)N(+)(H(2)O)(2) cluster and a fully protonated ferric hydroperoxide. All the calculations show that a stepwise mechanism that involves proton relay and O-O homolysis, in the rate-determining step, has a much lower barrier (>10 kcal/mol) than the corresponding fully concerted mechanism. The best fit of the calculated solvent kinetic isotope effect, to the experimental data, is obtained for the H(3)O(+)(H(2)O)(2) cluster. The calculated alpha-deuterium secondary kinetic isotope effect is inverse (0.95-0.98), but much less so than the experimental value (0.7). Possible reasons for this quantitative difference are discussed. Some probes are suggested that may enable experiment to distinguish the stepwise from the concerted mechanism.     

Model studies on p-hydroxybenzoate hydroxylase. The catalytic role of Arg-214 and Tyr-201 in the hydroxylation step

A model C-(4a)-flavinhydroperoxide (FlHOOH) is described that contains the tricyclic isoalloxazine moiety, the C-4a-hydroperoxide functionality, and a beta-hydroxyethyl group to model the effect of the 2'-OH group of the ribityl side chain of native FADHOOH. The electronic structures of this reduced flavin (H(3)()Fl(red)()), its N1 anion (H(2)()Fl(red)()(-)()), oxidized flavin (HFl(ox)()), and FlHOOH have been fully optimized at the B3LYP/ 6-31+G(d,p) level of theory. This model C-4a-flavinhydroperoxide is used to describe the transition state for the key step in the paradigm aromatic hydroxylase, p-hydroxybenzoate hydroxylase (PHBH): the oxidation of p-hydroxybenzoate (p-OHB). The Tyrosine-201 residue in PHBH is modeled by phenol, and Arginine-214 is modeled by guanidine. Electrophilic aromatic substitution proceeds by an S(N)2-like attack of the aromatic sextet of p-OHB phenolate anion on the distal oxygen of FlHOOH 3. The transition structure for oxygen atom transfer is fully optimized [B3LYP/6-31+G(d,p)] and has a classical activation barrier of 24.9 kcal/mol. These data suggest that the role of the Tyr-201 is to orient the p-OHB substrate and to properly align it for the oxygen transfer step. Although the negatively charged phenolate oxygen does activate the C-3 carbon of p-OHB phenolate anion toward oxidation relative to ortho oxidation of the carboxylate anion, it appears that H-bonding the Tyr-201 residue to this phenolic oxygen stabilizes both the ground state (GS) and the transition state (TS) approximately equally and therefore plays only a minor role, if any, in lowering the activation barrier. Complexation of p-OHB with guanidine has only a modest effect upon the oxidation barriers. When the complex is in the form of a salt-bridge (10a), the barrier is only slightly reduced. When the TSs are placed in THF solvent (COSMO) with full geometry optimization, salt-bridge TS-A is slightly favored (DeltaDeltaE() = 2.3 kcal/mol).     

On the mechanisms of oxidation of organic sulfides by H2O2 in aqueous solutions

The mechanism of oxidation of organic sulfides in aqueous solutions by hydrogen peroxide was investigated via ab initio calculations. Specifically, two reactions, hydrogen transfer of hydrogen peroxide to form water oxide and the oxidation of dimethyl sulfide (DMS) by hydrogen peroxide to form dimethyl sulfoxide, were studied as models of these processes in general. Solvent effects are included both via including explicitly water molecules and via the polarizable continuum model. The former was found to have a much more significant effect than the latter. When explicit water molecules are included, a mechanism different from those proposed in the literature was found. Specific interactions including hydrogen bonding with 2-3 water molecules can provide enough stabilization for the charge separation of the activation complex. The energy barrier of the oxidation of DMS by hydrogen peroxide was estimated to be 12.7 kcal/mol, within the experimental range of the oxidation of analogous compounds (10-20 kcal/mol). The major reaction coordinates of the reaction are the breaking of the O-O bond of H2O2 and the formation of the S-O bond, the transfer of hydrogen to the distal oxygen of hydrogen peroxide occurring after the system has passed the transition state. Reaction barriers of the hydrogen transfer of H2O2 are an average of 10 kcal/mol or higher than the reaction barriers of the oxidation of DMS. Therefore, a two-step oxidation mechanism in which, first, the transfer of a hydrogen atom occurs to form water oxide and, second, the transfer of oxygen to the substrate occurs is unlikely to be correct. Our proposed oxidation mechanism does not suggest a pH dependence of oxidation rate within a moderate range around neutral pH (i.e., under conditions in which hydronium and hydroxide ions do not participate directly in the reaction), and it agrees with experimental observations over moderate pH values. Also, without including a protonated solvent molecule, it has activation energies that correspond to measured activation energies.     

Catalytic mechanism and product specificity of the histone lysine methyltransferase SET7/9: an ab initio QM/MM-FE study with multiple initial structures

Histone lysine methylation is emerging as an important mechanism to regulate chromatin structure and gene activity. To provide theoretical understanding of its reaction mechanism and product specificity, ab initio quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations and molecular dynamics simulations have been carried out to investigate the histone lysine methyltransferase SET7/9. It is found that the methyl-transfer reaction catalyzed by SET7/9 is a typical in-line S(N)2 nucleophilic substitution reaction with a transition state of 70% dissociative character. The calculated average free energy barrier at the MP2(6-31+G) QM/MM level is 20.4 +/- 1.1 kcal/mol, consistent with the activation barrier of 20.9 kcal/mol estimated from the experimental reaction rate. The barrier fluctuation has a strong correlation with the nucleophilic attack distance and angle in the reactant complex. The calculation results show that the product specificity of SET7/9 as a monomethyltransferase is achieved by disrupting the formation of near-attack conformations for the dimethylation reaction.