Second half-reaction of nitric oxide synthase: computational insights into the initial step and key proposed intermediate

Density functional theory methods have been employed to investigate possible first steps in the second half-reaction of the mechanism of nitric oxide synthases (NOSs). In particular, reactions and complexes formed via transfer of either or both hydrogens of the substrates (NHA) -NHOH group to the Fe-bound O2 were considered. For each of these pathways, the effect of adding an extra electron from tetrahydrobiotperin (H4B) was also examined. The preferred initial pathway involves the simultaneous transfer of both hydrogens of the -NHOH group to the Fe(heme)-O2, without an additional electron, to give the Fe(heme)-HOOH species which lies only marginally higher in energy, 2.5 kcal mol(-1) or less, than the initial bound active site. An alternative mechanism in which only the -NH- proton of the -NHOH group is transferred to the Fe(heme)-O2 to give an Fe(heme)-OOH derivative is found to require only slightly more energy, approximately 2 kcal mol(-1). However, transfer of the proton back to the -NOH nitrogen occurs without a barrier at 298.15 K. Tetrahedral intermediates in which the Fe(heme)-O2 has attached at the guanidinium carbon (C(guan)) of NHA, that is, forms an Fe(heme)-O2-C(guan) link, have also been investigated. All examples of such species considered, that is, with or without hydrogen or electron transfers, lie significantly higher in energy by at least 29.0 kcal mol(-1) than the initial bound active site. Thus, it is suggested that such complexes are not mechanistically feasible. The implications of the present findings for the second half-reaction are also discussed.     

A density functional theory investigation on the mechanism of the second half-reaction of nitric oxide synthase

Density functional theory methods have been employed to systematically investigate the overall mechanism of the second half-reaction of nitric oxide synthases. The initial heme-bound hydrogen peroxide intermediate previously identified is found to first undergo a simple rotation about its O-O peroxide bond. Then, via a "ping-pong" peroxidase-like mechanism the -O(in)H- proton is transferred back onto the substrate's -NO oxygen then subsequently onto the outer oxygen of the resulting Fe(heme)-OOH species. As a result, O(out) is released as H2O with concomitant formation of a compound I-type (Fe(heme)-O) species. Formation of the final citrulline and NO products can then be achieved in one step via a tetrahedral transition structure resulting from direct attack of the Fe(heme)-O moiety at the substrate's guanidinium carbon center. The possible role of alternative mechanisms involving a protonated compound II-type species or an initial transfer of only the -NH- hydrogen of the =NHOH+ group of N(omega)-hydroxy-L-arginine is also discussed.     

Insights into the mechanism of O₂ formation and release from the Mn₄O₄L₆ "cubane" cluster

To probe photoinduced water oxidation catalyzed by the Mn?O?L? cubane clusters, we have computationally studied the mechanism and controlling factors of the O? formation from the [Mn?O?L?] catalyst, 6. It was demonstrated that dissociation of an L = H?PO?? ligand from 6 facilitates the direct O-O bond formation that proceeds with a 28.3 (33.4) kcal/mol rate-determining energy barrier at the transition state TS1. This step (the O-O single bond formation) of the reaction is a two-electron oxidation/reduction process, during which two oxo ligands are transformed into to μ2:η2-O?2? unit, and two ("distal") Mn centers are reduced from the 4+ to the 3+ oxidation state. Next two-electron oxidation/reduction occurs by "dancing" of the resulted O?2? fragment between the Mn1 and Mn2/Mn(2')-centers, keeping its strong coordination to the Mn(1')-center. As a result of this four-electron oxidation/reduction process Mn centers of the Mn?-core of I transform from {Mn1(III)-Mn(1')(III)-Mn2(IV)-Mn(2')(IV)} to {Mn1(II)-Mn(1')(II)-Mn2(III)-Mn(2')(III)} in IV. In other words, upon O? formation in cationic complex [Mn?O?L?](+), I, all four Mn-centers are reduced by one electron each. The overall reaction I → TS1 → II → III → TS2 → IV → TS3 → V → VI + O? is found to be exothermic by 15.4 (10.5) kcal/mol. We analyze the lowest spin states and geometries of all reactants, intermediates, transition states, and products of the targeted reaction.     

Mechanisms of the reactions of W AND W+ with H2O: computational studies

The mechanism of the reactions of W and W(+) with the water molecule have been studied for several lower-lying electronic states of tungsten centers at the CCSD(T)/6-311G(d,p)+SDD and B3LYP/6-31G(d,p)+SDD levels of theory. It is shown that these reactions are essentially multistate processes, during which lower-lying electronic states of the systems cross several times. They start with the formation of initial prereaction M(H(2)O) complexes with M-H(2)O bonding energies of 9.6 and 48.2 kcal/mol for M = W and W(+), followed by insertion of the metal center into an O-H bond with 20.0 and 53.3 kcal/mol barriers for neutral and cationic systems, respectively. The overall process of M + H(2)O --> t-HM(OH) is calculated to be highly exothermic, 48.4 and 48.8 kcal/mol for M = W and W(+). From the HM(OH) intermediate the reaction may proceed via several different channels, among which the stepwise HM(OH) --> HMO + H --> (H)(2)MO and concerted HM(OH) --> (H)(2)MO pathways are more favorable and can compete (energetically) with each other. For the neutral system (M = W), the concerted process is the most favorable, whereas for the charged system (M = W(+)), the stepwise pathway is slightly more favorable. From the energetically most favorable intermediate (H)(2)MO the reactions proceed via H(2)-molecule formation with a 53.1 kcal/mol activation barrier for the neutral system. For the cationic system, H-H formation and dissociation is an almost barrierless process. The overall reaction of W and W(+) with the water molecule leading to H(2) + MO formation is found to be exothermic by 48.2 and 39.8 kcal/mol, respectively. In the gas phase with the collision-less conditions the reactions W((7)S) + H(2)O --> H(2) + WO((3)Sigma(+)), and W(+)((6)D) + H(2)O --> H(2) + WO(+)((4)Sigma(+)) are expected to proceed via a 10.4 and 5.1 kcal/mol overall energy barrier corresponding to the first O-H dissociation at the TS1. On the basis of these PESs, we predict kinetic rate constants for the reactions of W and W(+) with H(2)O.     

Catalytic mechanism of dopamine beta-monooxygenase mediated by Cu(III)-oxo

Mechanisms of dopamine hydroxylation by the Cu(II)-superoxo species and the Cu(III)-oxo species of dopamine beta-monooxygenase (DBM) are discussed using QM/MM calculations for a whole-enzyme model of 4700 atoms. A calculated activation barrier for the hydrogen-atom abstraction by the Cu(II)-superoxo species is 23.1 kcal/mol, while that of the Cu(III)-oxo, which can be viewed as Cu(II)-O*, is 5.4 kcal/mol. Energies of the optimized radical intermediate in the superoxo- and oxo-mediated pathways are 18.4 and -14.2 kcal/mol, relative to the corresponding reactant complexes, respectively. These results demonstrate that the Cu(III)-oxo species can better mediate dopamine hydroxylation in the protein environment of DBM. The side chains of three amino acid residues (His415, His417, and Met490) coordinate to the Cu(B) atom, one of the copper sites in the catalytic core that plays a role for the catalytic function. The hydrogen-bonding network between dopamine and the three amino acid residues (Glu268, Glu369, and Tyr494) plays an essential role in substrate binding and the stereospecific hydroxylation of dopamine to norepinephrine. The dopamine hydroxylation by the Cu(III)-oxo species is a downhill and lower-barrier process toward the product direction with the aid of the protein environment of DBM. This enzyme is likely to use the high reactivity of the Cu(III)-oxo species to activate the benzylic C-H bond of dopamine; the enzymatic reaction can be explained by the so-called oxygen rebound mechanism.     

Reaction mechanism of naphthyl radicals with molecular oxygen. 1. Theoretical study of the potential energy surface

Potential energy surfaces (PESs) of the reactions of 1- and 2-naphthyl radicals with molecular oxygen have been investigated at the G3(MP2,CC)//B3LYP/6-311G** level of theory. Both reactions are shown to be initiated by barrierless addition of O(2) to the respective radical sites of C(10)H(7). The end-on O(2) addition leading to 1- and 2-naphthylperoxy radicals exothermic by 45-46 kcal/mol is found to be more preferable thermodynamically than the side-on addition. At the subsequent reaction step, the chemically activated 1- and 2-C(10)H(7)OO adducts can eliminate an oxygen atom leading to the formation of 1- and 2-naphthoxy radical products, respectively, which in turn can undergo unimolecular decomposition producing indenyl radical + CO via the barriers of 57.8 and 48.3 kcal/mol and with total reaction endothermicities of 14.5 and 10.2 kcal/mol, respectively. Alternatively, the initial reaction adducts can feature an oxygen atom insertion into the attacked C(6) ring leading to bicyclic intermediates a10 and a10' (from 1-naphthyl + O(2)) or b10 and b10' (from 2-naphthyl + O(2)) composed from two fused six-member C(6) and seven-member C(6)O rings. Next, a10 and a10' are predicted to decompose to C(9)H(7) (indenyl) + CO(2), 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H, and 1-C(9)H(7)O (1-benzopyranyl) + CO, whereas b10 and b10' would dissociate to C(9)H(7) (indenyl) + CO(2), 2-C(9)H(7)O (2-benzopyranyl) + CO, and 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H. On the basis of this, the 1-naphthyl + O(2) reaction is concluded to form the following products (with the overall reaction energies given in parentheses): 1-naphthoxy + O (-15.5 kcal/mol), indenyl + CO(2) (-123.9 kcal/mol), 1-benzopyranyl + CO (-97.2 kcal/mol), and 1,2-naphthoquinone + H (-63.5 kcal/mol). The 2-naphthyl + O(2) reaction is predicted to produce 2-naphthoxy + O (-10.9 kcal/mol), indenyl + CO(2) (-123.7 kcal/mol), 2-benzopyranyl + CO (-90.7 kcal/mol), and 1,2-naphthoquinone + H (-63.2 kcal/mol). Simplified kinetic calculations using transition-state theory computed rate constants at the high-pressure limit indicate that the C(10)H(7)O + O product channels are favored at high temperatures, while the irreversible oxygen atom insertion first leading to the a10 and a10' or b10 and b10' intermediates and then to their various decomposition products is preferable at lower temperatures. Among the decomposition products, indenyl + CO(2) are always most favorable at lower temperatures, but the others, 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H (from a10 and b10'), 1-C(9)H(7)O (1-benzopyranyl) + CO (from a10'), and 2-C(10)H(7)O (2-benzopyranyl) + O (from b10 and minor from b10'), may notably contribute or even become major products at higher temperatures.     

Theoretical study of [Ni (H2O)n]2+(H2O)m (n ≤ 6, m ≤ 18)

The [Ni-(H(2)O)(n)](2+)(H(2)O)(m) (n ≤ 6, m ≤ 18) complexes were studied by means of first-principles all-electron calculations performed with the BPW91 gradient corrected functional and the 6-311+G(d,p) basis sets for the H, O, and Ni atoms. Triplet states were found as low-lying states for each (n, m) combination. The estimated Ni(2+)-(H(2)O)(n) binding energies (112.8-57.4 kcal/mol for the first layer and 52.0-23.0 kcal/mol for the second one) decreases and the Ni(2+)-OH(2) bond lengths lengthen as n + m increases. With six H(2)O moieties the Ni(2+) ion furnishes its first coordination sphere of octahedral geometry. Further water addition renders the formation of the second layer. The effect of Ni(2+) on the (H(2)O)(n)···(H(2)O)(m) hydrogen bond formation for several "n" and "m" combinations was studied, revealing an enhancement of this kind of bonding, which is of key importance for the stabilization and growth of the clusters. For some n + m isomers the second layer appears before the first octahedral layer is fully formed. For example, the square planar Ni(2+)-(H(2)O)(4) core originates two-dimensional 4 + 2 and 4 + 4 isomers, where each outer water molecule accepts two H-bonds, lying 2.0 kcal/mol above the 6 and 6 + 2 ground states. The clusters were also studied by IR spectra; the OH stretching vibrational frequencies allowed us to identify the outer solvation shells by the presence of red-shifted hydrogen bond regions.     

Redox chemistry of the acetato-bridged clusters [M3(mu3-O)n(mu-O2CCH3)6(H2O)3]2+ (M = Mo, W, n = 1, 2): reversible redox between mono-micro3-oxo d8 M(III)2M(IV) and d9 M(III)3 forms

A cyclic voltammogram of aqueous 0.1 mol dm(-3) triflic acid solutions of the d6 bioxo-capped M-M bonded cluster [Mo3(mu3-O)2(O2CCH3)6(H2O)3]2+ at a glassy carbon electrode at 25 degrees C gives rise to an irreversible 3e- cathodic wave to a d9 Mo(III)3 species at -0.8 V vs. SCE which on the return scan gives rise to two anodic waves at +0.05 V vs. SCE (E(1/2), 1e- reversible to d8 Mo(III)2Mo(IV)) and +0.48 V vs. SCE (2e- irreversible back to d6 Mo(IV)3). The number of electrons passed at each redox wave has been confirmed by redox titration and controlled potential electrolysis which resulted in 90% recovery of [Mo3(mu3-O)2(O2CCH3)6(H2O)3]2+ following electrochemical re-oxidation at +0.8 V. A corresponding CV study of the d8 monoxo-capped W(III)2W(IV) cluster [W3(mu3-O)(O2CCH3)6(H2O)3]2+ gives rise to a reversible 1e- cathodic process at -0.92 V vs. SCE to give the d9 W(III)3 species [W3(mu3-O)(O2CCH3)6(H2O)3]+; the first authentic example of a W(III) complex with coordinated water ligands. However the cluster is too unstable (O2/water sensitive) to allow isolation. Comparisons with the cv study on [Mo3(mu3-O)2(O2CCH3)6(H2O)3]2+ suggest irreversible reduction of this complex to monoxo-capped [Mo(III)3(mu3-O)(O2CCH3)6(H2O)3]+ followed by reversible oxidation to its d8 counterpart [Mo3(mu3-O)(O2CCH3)6(H2O)3]2+ (Mo(III)2Mo(IV)) and finally irreversible oxidation back to the starting bioxo-capped cluster. Exposing the d9 Mo(III)3 cluster to air (O2) however gives a different final product with evidence of break up of the acetate bridged framework. Corresponding redox processes on d6 [W3(mu3-O)2(O2CCH3)6(H2O)3]2+ are too cathodic to allow similar generation of the monoxo-capped W(III)3 and W(III)2W(IV) clusters at the electrode surface.     

Active Thermochemical Tables: accurate enthalpy of formation of hydroperoxyl radical, HO2

Through the use of the Active Thermochemical Tables approach, the best currently available enthalpy of formation of HO2 has been obtained as delta(f)H(o)298 (HO2) = 2.94 +/- 0.06 kcal mol(-1) (3.64 +/- 0.06 kcal mol(-1) at 0 K). The related enthalpy of formation of the positive ion, HO2+, within the stationary electron convention is delta(f)H(o)298 (HO2+) = 264.71 +/- 0.14 kcal mol(-1) (265.41 +/- 0.14 kcal mol(-1) at 0 K), while that for the negative ion, HO2- (within the same convention), is delta(f)H(o)298 (HO2-) = -21.86 +/- 0.11 kcal mol(-1) (-21.22 +/- 0.11 kcal mol(-1) at 0 K). The related proton affinity of molecular oxygen is PA298(O2) = 100.98 +/- 0.14 kcal mol(-1) (99.81 +/- 0.14 kcal mol(-1) at 0 K), while the gas-phase acidity of H2O2 is delta(acid)G(o)298 (H2O2) = 369.08 +/- 0.11 kcal mol(-1), with the corresponding enthalpy of deprotonation of H2O2 of delta(acid)H(o)298 (H2O2) = 376.27 +/- 0.11 kcal mol(-1) (375.02 +/- 0.11 kcal mol(-1) at 0 K). In addition, a further improved enthalpy of formation of OH is briefly outlined, delta(f)H(o)298 (OH) = 8.93 +/- 0.03 kcal mol(-1) (8.87 +/- 0.03 kcal mol(-1) at 0 K), together with new and more accurate enthalpies of formation of NO, delta(f)H(o)298 (NO) = 21.76 +/- 0.02 kcal mol(-1) (21.64 +/- 0.02 kcal mol(-1) at 0 K) and NO2, delta(f)H(o)298 (NO2) = 8.12 +/- 0.02 kcal mol(-1) (8.79 +/- 0.02 kcal mol(-1) at 0 K), as well as H(2)O(2) in the gas phase, delta(f)H(o)298 (H2O2) = -32.45 +/- 0.04 kcal mol(-1) (-31.01 +/- 0.04 kcal mol(-1) at 0 K). The new thermochemistry of HO2, together with other arguments given in the present work, suggests that the previous equilibrium constant for NO + HO2 --> OH + NO2 was underestimated by a factor of approximately 2, implicating that the OH + NO2 rate was overestimated by the same factor. This point is experimentally explored in the companion paper of Srinivasan et al. (next paper in this issue).     

Adsorption of O, H, OH, and H2O on Ag(100)

The adsorption of H(2)O and its dissociation products, O, H, and OH, on Ag(100) has been studied using an ab initio embedding method. Results at different sites (atop, bridge, and hollow) are presented. The four-fold hollow site is found to be the most stable adsorption site for O, H, and OH, and the calculated adsorption energies are 87.1, 42.7, and 76.2 kcal mol(-1), respectively. The adsorption energy of water at the atop and bridge sites is almost identical with values of 11.1 and 12.0 kcal mol(-1), respectively. The formation of adsorbed OH species by adsorption of water on oxygen-precovered Ag(100) is predicted to be exothermic by 36 kcal mol(-1).     

A direct dynamics trajectory study of F- + CH(3)OOH reactive collisions reveals a major non-IRC reaction path

A direct dynamics simulation at the B3LYP/6-311+G(d,p) level of theory was used to study the F- + CH3OOH reaction dynamics. The simulations are in excellent agreement with a previous experimental study (J. Am. Chem. Soc. 2002, 124, 3196). Two product channels, HF + CH2O + OH- and HF + CH3OO-, are observed. The former dominates and occurs via an ECO2 mechanism in which F- attacks the CH3- group, abstracting a proton. Concertedly, a carbon-oxygen double bond is formed and OH- is eliminated. Somewhat surprisingly this is not the reaction path, predicted by the intrinsic reaction coordinate (IRC), which leads to a deep potential energy minimum for the CH2(OH)2...F- complex followed by dissociation to HF + CH2(OH)O-. None of the direct dynamics trajectories followed this path, which has an energy release of -63 kcal/mol and is considerably more exothermic than the ECO2 path whose energy release is -27 kcal/mol. Other product channels not observed, and which have a lower energy than that for the ECO2 path, are F- + CO + H2 + H2O (-43 kcal/mol), F- + CH2O + H2O (-51 kcal/mol), and F- + CH2(OH)2 (-60 kcal/mol). Formation of the CH3OOH...F- complex, with randomization of its internal energy, is important, and this complex dissociates via the ECO2 mechanism. Trajectories which form HF + CH3OO- are nonstatistical events and, for the 4 ps direct dynamics simulation, are not mediated by the CH3OOH...F- complex. Dissociation of this complex to form HF + CH3OO- may occur on longer time scales.     

Stepwise hydration of ionized aromatics. Energies, structures of the hydrated benzene cation, and the mechanism of deprotonation reactions

The stepwise binding energies (DeltaHdegree(n-1,n)) of 1-8 water molecules to benzene(.+) [Bz(.+)(H2O)n] were determined by equilibrium measurements using an ion mobility cell. The stepwise hydration energies, DeltaHdegree(n-1,n), are nearly constant at 8.5 +/- 1 kcal mol-1 from n = 1-6. Calculations show that in the n = 1-4 clusters, the benzene(.+) ion retains over 90% of the charge, and it is extremely solvated, that is, hydrogen bonded to an (H2O)n cluster. The binding energies and entropies are larger in the n = 7 and 8 clusters, suggesting cyclic or cage-like water structures. The concentration of the n = 3 cluster is always small, suggesting that deprotonation depletes this ion, consistent with the thermochemistry since associative deprotonation Bz(.+)(H2O)(n-1) + H2O-->C6H5. + (H2O)nH+ is thermoneutral or exothermic for n > or = 4. Associative intracluster proton transfer Bz(.+)(H2O)(n+1) + H2O-->C6H5.(H2O)nH+ would also be exothermic for n > or = 4, but lack of H/D exchange with D2O shows that the proton remains on C6H6(.+) in the observed Bz(.+)(H2O)n clusters. This suggests a barrier to intracluster proton transfer, and as a result, the [Bz(.+)(H2O)n]* activated complexes either undergo dissociative proton transfer, resulting in deprotonation and generation of (H2O)nH+, or become stabilized. The rate constant for the deprotonation reaction shows a uniquely large negative temperature coefficient of K = cT(-67+/-4) (or activation energy of -34+/- 1 kcal mol-1), caused by a multibody mechanism in which five or more components need to be assembled for the reaction.     

Is the protein surrounding the active site critical for hydrogen peroxide reduction by selenoprotein glutathione peroxidase? An ONIOM study

In this ONIOM(QM:MM) study, we evaluate the role of the protein surroundings in the mechanism of H2O2 reduction catalyzed by the glutathione peroxidase enzyme, using the whole monomer (3113 atoms in 196 amino acid residues) as a model. A new optimization scheme that allows the full optimization of transition states for large systems has been utilized. It was found that in the presence of the surrounding protein the optimized active site structure bears a closer resemblance to the one in the X-ray structure than that without the surrounding protein. H2O2 reduction occurs through a two-step mechanism. In the first step, the selenolate anion (E-Se(-)) formation occurs with a barrier of 16.4 kcal/mol and is endothermic by 12.0 kcal/mol. The Gln83 residue plays the key role of the proton abstractor, which is in line with the experimental suggestion. In the second step, the O-O bond is cleaved, and selenenic acid (R-Se-OH) and a water molecule are formed. The calculated barrier for this process is 6.0 kcal/mol, and it is exothermic by 80.9 kcal/mol. The overall barrier of 18.0 kcal/mol for H2O2 reduction is in reasonable agreement with the experimentally measured barrier of 14.9 kcal/mol. The protein surroundings has been calculated to exert a net effect of only 0.70 kcal/mol (in comparison to the "active site only" model including solvent effects) on the overall barrier, which is most likely due to the active site being located at the enzyme surface.     

New reaction model for O-O bond formation and O2 evolution catalyzed by dinuclear manganese complex

A new mechanism of the oxygen evolving reaction catalyzed by [H(2)O(terpy)Mn(μ-O)(2)Mn(terpy)OH(2)](3+) is proposed by using density functional theory. This proton coupled electron transfer (PCET) model shows reasonable barriers. Because in experiments excess oxidants (OCl(-) or HSO(5)(-)) are required to evolve oxygen from water, we considered the Mn(2) complex neutralized by three counterions. Structure optimization made the coordinated OCl(-) withdraw a H(+) from the water ligand and produces the reaction space for H(2)O(2) formation with the deprotonated OH(-) ligand. The reaction barrier for the H(2)O(2) formation from OH(-) and protonated OCl(-) depends significantly on the system charge and is 14.0 kcal/mol when the system is neutralized. The H(2)O(2) decomposes to O(2) during two PCET processes to the Mn(2) complex, both with barriers lower than 12.0 kcal/mol. In both PCET processes the spin moment of transferred electrons prefers to be parallel to that of Mn 3d electrons because of the exchange interaction. This model thus explains how the triplet O(2) molecule is produced.     

New insight into the gas-phase bimolecular self-reaction of the HOO radical

The singlet and triplet potential energy surfaces (PESs) for the gas-phase bimolecular self-reaction of HOO*, a key reaction in atmospheric environments, have been investigated by means of quantum-mechanical electronic structure methods (CASSCF and CASPT2). All the reaction pathways on both PESs consist of a first step involving the barrierless formation of a prereactive doubly hydrogen-bonded complex, which is a diradical species lying about 8 kcal/mol below the energy of the reactants at 0 K. The lowest energy reaction pathway on both PESs is the degenerate double hydrogen exchange between the HOO* moieties of the prereactive complex via a double proton transfer mechanism involving an energy barrier of only 1.1 kcal/mol for the singlet and 3.3 kcal/mol for the triplet at 0 K. The single H-atom transfer between the two HOO* moieties of the prereactive complex (yielding HOOH + O2) through a pathway keeping a planar arrangement of the six atoms involves a conical intersection between either two singlet or two triplet states of A' and A" symmetries. Thus, the lowest energy reaction pathway occurs via a nonplanar cisoid transition structure with an energy barrier of 5.8 kcal/mol for the triplet and 17.5 kcal/mol for the singlet at 0 K. The simple addition between the terminal oxygen atoms of the two HOO* moieties of the prereactive complex, leading to the straight chain H2O4 intermediate on the singlet PES, involves an energy barrier of 7.3 kcal/mol at 0 K. Because the decomposition of such an intermediate into HOOH + O2 entails an energy barrier of 45.2 kcal/mol at 0 K, it is concluded that the single H-atom transfer on the triplet PES is the dominant pathway leading to HOOH + O2. Finally, the strong negative temperature dependence of the rate constant observed for this reaction is attributed to the reversible formation of the prereactive complex in the entrance channel rather than to a short-lived tetraoxide intermediate.     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

Chemo- and periselectivity in the addition of [OsO2(CH2)2] to ethylene: a theoretical study

Quantum chemical calculations by using density functional theory at the B3LYP level have been carried out to elucidate the reaction course for the addition of ethylene to [OsO2(CH2)2] (1). The calculations predict that the kinetically most favorable reaction proceeds with an activation barrier of 8.1 kcal mol(-1) via [3+2] addition across the O=Os=CH2 moiety. This reaction is -42.4 kcal mol(-1) exothermic. Alternatively, the [3+2] addition to the H2C=Os=CH2 fragment of 1 leads to the most stable addition product 4 (-72.7 kcal mol(-1)), yet this process has a higher activation barrier (13.0 kcal mol(-1)). The [3+2] addition to the O=Os=O fragment yielding 2 is kinetically (27.5 kcal mol(-1)) and thermodynamically (-7.0 kcal mol(-1)) the least favorable [3+2] reaction. The formal [2+2] addition to the Os=O and Os=CH2 double bonds proceeds by initial rearrangement of 1 to the metallaoxirane 1 a. The rearrangement 1-->1 a and the following [2+2] additions have significantly higher activation barriers (>30 kcal mol(-1)) than the [3+2] reactions. Another isomer of 1 is the dioxoosmacyclopropane 1 b, which is 56.2 kcal mol(-1) lower in energy than 1. The activation barrier for the 1-->1 b isomerization is 15.7 kcal mol(-1). The calculations predict that there are no energetically favorable addition reactions of ethylene with 1 b. The isomeric form 1 c containing a peroxo group is too high in energy to be relevant for the reaction course. The accuracy of the B3LYP results is corroborated by high level post-HF CCSD(T) calculations for a subset of species.     

Proton transfer to nickel-thiolate complexes. 1. Protonation of [Ni(SC(6)H(4)R-4)(2)(Ph(2)PCH(2)CH(2)PPh(2))] (R = Me, MeO, H, Cl, or NO(2))

The kinetics of the equilibrium reaction between [Ni(SC(6)H(4)R-4)(2)(dppe)] (R= MeO, Me, H, Cl, or NO(2); dppe = Ph(2)PCH(2)CH(2)PPh(2)) and mixtures of [lutH](+) and lut (lut = 2,6-dimethylpyridine) in MeCN to form [Ni(SHC(6)H(4)R-4)(SC(6)H(4)R-4)(dppe)](+) have been studied using stopped-flow spectrophotometry. The kinetics for the reactions with R = MeO, Me, H, or Cl are consistent with a single-step equilibrium reaction. Investigation of the temperature dependence of the reactions shows that DeltaG = 13.6 +/- 0.3 kcal mol(-)(1) for all the derivatives but the values of DeltaH and DeltaS vary with R (R = MeO, DeltaH() = 8.5 kcal mol(-)(1), DeltaS = -16 cal K(-)(1) mol(-)(1); R = Me, DeltaH() = 10.8 kcal mol(-)(1), DeltaS = -9.5 cal K(-)(1) mol(-)(1); R = Cl, DeltaH = 23.7 kcal mol(-)(1), DeltaS = +33 cal K(-)(1) mol(-)(1)). With [Ni(SC(6)H(4)NO(2)-4)(2)(dppe)] a more complicated rate law is observed consistent with a mechanism in which initial hydrogen-bonding of [lutH](+) to the complex precedes intramolecular proton transfer. It seems likely that all the derivatives operate by this mechanism, but only with R = NO(2) (the most electron-withdrawing substituent) does the intramolecular proton transfer step become sufficiently slow to result in the change in kinetics. Studies with [lutD](+) show that the rates of proton transfer to [Ni(SC(6)H(4)R-4)(2)(dppe)] (R = Me or Cl) are associated with negligible kinetic isotope effect. The possible reasons for this are discussed. The rates of proton transfer to [Ni(SC(6)H(4)R-4)(2)(dppe)] vary with the 4-R-substituent, and the Hammett plot is markedly nonlinear. This unusual behavior is attributable to the electronic influence of R which affects the electron density at the sulfur.     

Decomposition kinetics of Ni(III)-peptide complexes with histidine and histamine as the third residue

The decomposition kinetics of the Ni(III) complexes of Gly(2)HisGly and Gly(2)Ha are studied from p[H(+)] 3.5 to 10, where His is l-histidine and Ha is histamine. In these redox reactions, at least two Ni(III) complexes are reduced to Ni(II) while oxidizing a single peptide ligand. The rate of Ni(III) loss is first order at low pH, mixed order from pH 7.0 to 8.5, and second order at higher pH. The transition from first- to second-order kinetics is attributed to the formation of an oxo-bridged Ni(III)-peptide dimer. The rates of decay of the Ni(III) complexes are general-base assisted with Br?nsted beta values of 0.62 and 0.59 for Ni(III)Gly(2)HisGly and Ni(III)Gly(2)Ha, respectively. The coordination of Gly(2)HisGly and Gly(2)Ha to Ni(II) are examined by UV-vis and CD spectroscopy. The square planar Ni(II)(H(-2)Gly(2)HisGly)(-) and Ni(II)(H(-2)Gly(2)Ha) complexes lose an additional proton from an imidazole nitrogen at high pH with pK(a) values of 11.74 and 11.54, respectively. The corresponding Ni(III) complexes have axially coordinated water molecules with pK(a) values of 9.37 and 9.44. At higher pH an additional proton is lost from the imidazole nitrogen with a pK(a) value of 10.50 to give Ni(III)(H(-3)Gly(2)Ha)(H(2)O)(OH)(2-).