Antibiotic deactivation by a dizinc beta-lactamase: mechanistic insights from QM/MM and DFT studies

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods and density functional theory (DFT) were used to investigate the initial ring-opening step in the hydrolysis of moxalactam catalyzed by the dizinc L1 beta-lactamase from Stenotrophomonas maltophilia. Anchored at the enzyme active site via direct metal binding as suggested by a recent X-ray structure of an enzyme-product complex (Spencer, J.; et al. J. Am. Chem. Soc. 2005, 127, 14439), the substrate is well aligned with the nucleophilic hydroxide that bridges the two zinc ions. Both QM/MM and DFT results indicate that the addition of the hydroxide nucleophile to the carbonyl carbon in the substrate lactam ring leads to a metastable intermediate via a dominant nucleophilic addition barrier. The potential of mean force obtained by SCC-DFTB/MM simulations and corrected by DFT/MM calculations yields a reaction free energy barrier of 23.5 kcal/mol, in reasonable agreement with the experimental value of 18.5 kcal/mol derived from kcat of 0.15 s(-1). It is further shown that zinc-bound Asp120 plays an important role in aligning the nucleophile, but accepts the hydroxide proton only after the nucleophilic addition. The two zinc ions are found to participate intimately in the catalysis, consistent with the proposed mechanism. In particular, the Zn(1) ion is likely to serve as an "oxyanion hole" in stabilizing the carbonyl oxygen, while the Zn(2) ion acts as an electrophilic catalyst to stabilize the anionic nitrogen leaving group.     

Mechanism of H2O2 dismutation catalyzed by a new catalase mimic (a non-heme dibenzotetraaza[14]annulene-Fe(III) complex): a density functional theory investigation

The mechanism of H(2)O(2) dismutation catalyzed by the dibenzotetraaza[14]annulene-Fe(III) complex ([Fe(C(24)H(22)N(4)O(4))](+)) which was recently reported (Paschke, J.; Kirsch, M.; Korth, H. G.; de Groot, H.; Sustmann, R. J. Am. Chem. Soc. 2001, 123, 11099) has been investigated by density functional theory using the B3LYP hybrid functional. The quartet potential energy profile of the catalytic reaction has been explored. In the whole catalytic cycle, the rate-determining step is found to be the O-O bond homolytic cleavage, without the assistance of solvent molecules in the second coordination shell. The calculated free energy barrier for this step is 10.8 kcal/mol, which is in reasonable agreement with the experimental facts. The calculations also show that the hydroxyl and hydroperoxyl radicals may be generated in the reaction processes, but they can be efficiently quenched in strongly exothermic steps. This study provides a satisfactory explanation to the observed efficiency of the H(2)O(2) dismutation catalyzed by this complex.     

Glycine enolates: the effect of formation of iminium ions to simple ketones on alpha-amino carbon acidity and a comparison with pyridoxal iminium ions

Equilibrium constants in D2O were determined by 1H NMR analyses for formation of imines/iminium ions from addition of glycine methyl ester to acetone and from addition of glycine to phenylglyoxylate. First-order rate constants, also determined by 1H NMR, are reported for deuterium exchange between solvent D2O and the alpha-amino carbon of glycine methyl ester and glycine in the presence of increasing concentrations of ketone and Br?nsted bases. These rate and equilibrium data were used to calculate second-order rate constants for deprotonation by DO- and by Br?nsted bases of the alpha-imino carbon of the ketone adducts. Formation of the iminium ion between acetone and glycine methyl ester and between phenylglyoxylate and glycine is estimated to cause 7 unit and 15 unit decreases, respectively, in the pKa's of 21 and 29 for deprotonation of the parent carbon acids. The effect of formation of iminium ions to phenylglyoxylate and to 5'-deoxypyridoxal (DPL) [Toth, K.; Richard, J. P. J. Am. Chem. Soc. 2007, 129, 3013-3021] on the carbon acidity of glycine is similar. However, DPL is a much better catalyst than phenylglyoxylate of deprotonation of glycine, because of the exceptionally large thermodynamic driving force for conversion of the amino acid and DPL to the reactive iminium ion.     

Structure-activity studies on the cleavage of an RNA analogue by a potent dinuclear metal ion catalyst: effect of changing the metal ion

Dinuclear Cd(II), Cu(II), and Zn(II) complexes of L2OH (L2OH = 1,3-bis(1,4,7-triazacyclonon-1-yl)-2-hydroxypropane) are compared as catalysts for cleavage of the RNA analogue HpPNP (HpPNP = 2-hydroxypropyl 4-nitrophenyl phosphate) at 25 degrees C, I = 0.10 M (NaNO(3)). Zn(II) and Cu(II) readily form dinuclear complexes at millimolar concentrations and a 2:1 ratio of metal ion to L2OH at neutral pH. The dinuclear Zn(2)(L2O) and Cu(2)(L2O) complexes have a bridging alkoxide group that brings together the two cations in close proximity to facilitate cooperative catalysis. Under similar conditions, the dinuclear complex of Cd(II) is a minor species in solution; only at high pH values (pH 10.4) does the Cd(2)(L2O) complex become the predominant species in solution. Analysis of the second-order rate constants for cleavage of HpPNP by Zn(2)(L2O) is straightforward because a linear dependence of pseudo-first-order rate constant on dinuclear complex is observed over a wide pH range. In contrast, plots of pseudo-first-order rate constants for cleavage of HpPNP by solutions containing a 2:1 ratio of Cd(II) to L2OH as a function of increasing L2OH are curved, and second-order rate constants are obtained by fitting the kinetic data to an equation for the formation of the dinuclear Cd(II) complex as a function of pH and [L2OH]. Second-order rate constants for cleavage of HpPNP by these dinuclear complexes at pH 9.3 and 25 degrees C vary by 3 orders of magnitude in the order Cd(2)(L2O) (2.8 M(-)(1) s(-)(1)) > Zn(2)(L2O) (0.68 M(-)(1) s(-)(1)) > Cu(2)(L2O) (0.0041 M(-1) s(-1)). The relative reactivity of these complexes is discussed in terms of the different geometric preferences and Lewis acidity of the dinuclear Zn(II), Cu(II), and Cd(II) complexes, giving insight into the importance of these catalyst properties in the cleavage of phosphate diesters resembling RNA.     

A combined ONIOM quantum chemical-molecular dynamics study of zinc-uracil bond breaking in yeast cytosine deaminase

A QM/MM method that combines ONIOM quantum chemistry and molecular dynamics is developed and applied to a step in the deamination of cytosine to uracil in yeast cytosine deaminase (yCD). A two-layer ONIOM calculation is used for the reaction complex, with an inner part treated at a high level for the chemical reaction (bond breaking) and a middle part treated at a lower level for relevant protein residues that are frozen in the quantum optimization. An outer layer (protein and solvent) is treated using MD. Configurations for the entire system are generated using MD and optimized with ONIOM. The method permits the use of high-level quantum calculations along with sufficient configurational sampling to approximate the potential of mean force for certain bond-breaking reactions. A previously proposed reaction mechanism for deamination (Sklenak, S.; Yao, L. S.; Cukier, R. I.; Yan, H. G. J. Am. Chem. Soc. 2004, 126, 14879) requires breaking the bond between a catalytic zinc and the O4 of uracil in order to permit product release. Using an ONIOM approach, direct bond cleavage was found to be energetically unfavorable. In the work presented here, the combined ONIOM MD method is used to show that the barrier for bond cleavage is small, approximately 3 kcal/mol, and, consequently, should not be the rate-limiting step in the reaction.     

Linear free energy relationships and kinetic isotope effects reveal the chemistry of the Ado 2′-OH group

Using kinetic isotope effects (KIE) and Hammett correlations, we show that the main role of the adenosine 2′-OH group on deprotonation by the non nucleophilic base DBU during external acyl group transfer is to generate enhanced electron density on the attacking nucleophile through ionization. The small primary KIEs (1.2 and 1.6) and the large Hammett reaction constants (+2.25 and +3.19) obtained for the ethanolysis of 2′/3′-O-p-substituted benzoyl 5′-O-trityl adenosines and 2′-deoxyadenosines are consistent with an A_N + D_N reaction mechanism. The implications of our results are discussed in terms of chemical contributions of the 2′-OH group in the ribosome catalysis of peptide bond formation.     

Computational study of the reaction of C6F6 with [IrMe(PEt3)3]: identification of a phosphine-assisted C-F activation pathway via a metallophosphorane intermediate

Density functional theory calculations have been used to model the reaction of C6F6 with [IrMe(PEt3)3], which proceeds with both C-F and P-C bond activation to yield trans-[Ir(C6F5)(PEt3)2(PEt2F)], C2H4, and CH4 (Blum, O.; Frolow, F.; Milstein, D. J. Chem. Soc., Chem. Commun. 1991, 258). Using a model species, trans-[IrMe(PH3)2(PH2Et)], a low-energy mechanism involving nucleophilic attack of the electron-rich Ir metal center at C6F6 with displacement of fluoride has been identified. A novel feature of this process is the capture of fluoride by a phosphine ligand to generate a metallophosphorane intermediate [Ir(C6F5)(Me)(PH3)2(PH2EtF)]. These events occur in a single step via a 4-centered transition state, in a process that we have termed "phosphine-assisted C-F activation". Alternative mechanisms based on C-F activation via concerted oxidative addition or electron-transfer processes proved less favorable. From the metallophosphorane intermediate the formation of the final products can be accounted for by facile ethyl group transfer from phosphorus to iridium followed by beta-H elimination of ethene and reductive elimination of methane. The interpretation of phosphine-assisted C-F activation in terms of nucleophilic attack is supported by the reduced activation barriers computed with the more electron-rich model reactant trans-[IrMe(PMe3)2(PMe2Et)] and the higher barriers found with lesser fluorinated arenes. Reactivity patterns for a range of fluoroarenes indicate the dominance of the presence of ortho-F substituents in promoting phosphine-assisted C-F activation, and an analysis of the charge distribution and transition state geometries indicates that this process is controlled by the strength of the Ir-aryl bond that is being formed.     

Density functional theory analysis of dimethylphosphate hydrolysis: effect of solvation and nucleophile variation

We model the hydrolytic cleavage of dimethylphosphate by hydroxide ion and water in the gas phase and in implicit water using density functional theory. In all cases the rate-determining step is the nucleophilic attack. The barrier for this nucleophilic attack in vacuum is much larger for the hydroxide than for the neutral nucleophile. However, in water the barriers are similar. The rate-determining step in the attack by the neutral nucleophile involves the concerted transfer of a proton from the water molecule to the phosphate ion and the formation of the P–O bond.     

A concerted mechanism for the transfer of the thiophosphinoyl group from aryl dimethylphosphinothioate esters to oxyanionic nucleophiles in aqueous solution

Earlier work on the hydrolysis of aryl phosphinothioate esters has led to contradictory mechanistic conclusions. To resolve this mechanistic ambiguity, we have measured linear free energy relationships (beta(nuc) and beta(lg)) and kinetic isotope effects for the reactions of oxyanions with aryl dimethylphosphinothioates. For the attack of nucleophiles on 4-nitrophenyl dimethylphosphinothioate, beta(nuc) = 0.47 +/- 0.05 for phenoxide nucleophiles (pK(a) < 11) and beta(nuc) = 0.08 +/- 0.01 for hydroxide and alkoxide nucleophiles (pK(a) >or= 11). Linearity of the plot in the range that straddles the pK(a) of the leaving group (4-nitrophenoxide, pK(a) 7.14) is indicative of a concerted mechanism. The much lower value of beta(nuc) for the more basic nucleophiles reveals the importance of a desolvation step prior to rate-limiting nucleophilic attack. The reactions of a series of substituted aryl dimethylphosphinothioate esters give the same value of beta(lg) with the nucleophiles HO(-) (beta= -0.54 +/- 0.03) and PhO(-) (beta = -0.52 +/- 0.09). A significantly better Hammett correlation is obtained with sigma(-) than with sigma or sigma degrees , as expected for a transition state involving rate-limiting cleavage of the P-OAr bond. The (18)O KIE at the position of bond fission ((18)k = 1.0124 +/- 0.0008) indicates the P-O bond is approximately 40% broken, and the (15)N KIE in the leaving group ((15)k = 1.0009 +/- 0.0003) reveals the nucleofuge carries about a third of a negative charge in the transition state. Thus, both the LFER and KIE data are consistent with a concerted reaction and disfavor a stepwise mechanism.     

Alkene cis-dihydroxylation by [(Me3tacn)(CF3CO2)RuVI O2)]ClO4(Me(3)tacn = 1,4,7-Trimethyl-1,4,7-triazacyclononane): structural characterization of [3 + 2] cycloadducts and kinetic studies

cis-Dioxoruthenium(VI) complex [(Me(3)tacn)(CF(3)CO(2))Ru(VI)O(2)]ClO(4) (1, Me(3)tacn = 1,4,7-trimethyl-1,4,7-triazacyclononane) reacted with alkenes in aqueous tert-butyl alcohol to afford cis-1,2-diols in excellent yields under ambient conditions. When the reactions of 1 with alkenes were conducted in acetonitrile, oxidative C=C cleavage reaction prevailed giving carbonyl products in >90% yields without any cis-diol formation. The alkene cis-dihydroxylation and C=C cleavage reactions proceed via the formation of a [3 + 2] cycloadduct between 1 and alkenes, analogous to the related reactions with alkynes [Che et al. J. Am. Chem. Soc. 2000, 122, 11380]. With cyclooctene and trans-beta-methylstyrene as substrates, the Ru(III) cycloadducts (4a) and (4b) [formula; see text] were isolated and structurally characterized by X-ray crystal analyses. The kinetics of the reactions of 1 with a series of p-substituted styrenes has been studied in acetonitrile by stopped-flow spectrophotometry. The second-order rate constants varied by 14-fold despite an overall span of 1.3 V for the one-electron oxidation potentials of alkenes. Secondary kinetic isotope effect (KIE) was observed for the oxidation of beta-d(2)-styrene (k(H)/k(D) = 0.83 +/- 0.04) and alpha-deuteriostyrene (k(H)/k(D) = 0.96 +/- 0.03), which, together with the stereoselectivity of cis-alkene oxidation by 1, is in favor of a concerted mechanism.     

Protonated benzofuran,anthracene, naphthalene,benzene, ethene,and ethyne: measurements and estimates of pK(a) and pK(R)

Aqueous solvolyses of acyl derivatives of hydrates (water adducts) of anthracene and benzofuran yield carbocations which undergo competitive deprotonation to form the aromatic molecules and nucleophilic reaction with water to give the aromatic hydrates. Trapping experiments with azide ions yield rate constants k(p) for the deprotonation and k(H2O) for the nucleophilic reaction based on the "azide clock". Combining these with rate constants for (a) the H(+)-catalyzed reaction of the hydrate to form the carbocation and (b) hydrogen isotope exchange of the aromatic molecule (from the literature) yields pK(R) = -6.0 and -9.4 and pK(a) = -13.5 and -16.3 for the protonated anthracene and protonated benzofuran, respectively. These pK values may be compared with pK(R) = -6.7 for naphthalene hydrate (1-hydroxy-1,2-dihydronaphthalene), extrapolated to water from measurements by Pirinccioglu and Thibblin for acetonitrile-water mixtures, and pK(a) = -20.4 for the 2-protonated naphthalene from combining k(p) with an exchange rate constant. The differences between pK(R) and pK(a) correspond to pK(H2O), the equilibrium constant for hydration of the aromatic molecule (pK(H2O) = pK(R) - pK(a)). For naphthalene and anthracene values of pK(H2O) = +13.7 and +7.5 compare with independent estimates of +14.2 and +7.4. For benzene, pK(a) = -24.3 is derived from an exchange rate constant and an assigned value for the reverse rate constant close to the limit for solvent relaxation. Combining this pK(a) with calculated values of pK(H2O) gives pK(R) = -2.4 and -2.1 for protonated benzenes forming 1,2- and 1,4-hydrates, respectively. Coincidentally, the rate constant for protonation of benzene is similar to those for protonation of ethylene and acetylene (Lucchini, V.; Modena, G. J. Am. Chem Soc. 1990, 112, 6291). Values of pK(a) for the ethyl and vinyl cations (-24.8) may thus be derived in the same way as that for the benzenonium ion. Combining these with appropriate values of pK(H2O) then yields pK(R) = -39.8 and -29.6 for the vinyl and ethyl cations, respectively.     

Probing the transition states of four glucoside hydrolyses with 13C kinetic isotope effects measured at natural abundance by NMR spectroscopy

Kinetic isotope effects (KIEs) were measured for methyl glucoside (4) hydrolysis on unlabeled material by NMR. Twenty-eight (13)C KIEs were measured on the acid-catalyzed hydrolysis of alpha-4 and beta-4, as well as enzymatic hydrolyses with yeast alpha-glucosidase and almond beta-glucosidase. The 1-(13)C KIEs on the acid-catalyzed reactions of alpha-4 and beta-4, 1.007(2) and 1.010(6), respectively, were in excellent agreement with the previously reported values (1.007(1), 1.011(2): Bennet and Sinnott, J. Am. Chem. Soc. 1986, 108, 7287). Transition state analysis of the acid-catalyzed reactions using the (13)C KIEs, along with the previously reported (2)H KIEs, confirmed that both reactions proceed with a stepwise D(N)A(N) mechanism and showed that the glucosyl oxocarbenium ion intermediate exists in an E(3) sofa or (4)H(3) half-chair conformation. (13)C KIEs showed that the alpha-glucosidase reaction also proceeded through a D(N)*A(N) mechanism, with a 1-(13)C KIE of 1.010(4). The secondary (13)C KIEs showed evidence of distortions in the glucosyl ring at the transition state. For the beta-glucosidase-catalyzed reaction, the 1-(13)C KIE of 1.032(1) demonstrated a concerted A(N)D(N) mechanism. The pattern of secondary (13)C KIEs was similar to the acid-catalyzed reaction, showing no signs of distortion. KIE measurement at natural abundance makes it possible to determine KIEs much more quickly than previously, both by increasing the speed of KIE measurement and by obviating the need for synthesis of isotopically labeled compounds.     

The coupling of isonitriles and carboxylic acids occurring by sequential concerted rearrangement mechanisms

Mechanisms for the recently described reactions of isonitriles with carboxylic acids (Li, X.; Danishefsky, S. J. J. Am. Chem. Soc. 2008, 130, 5446) are explored with the B3LYP density functional method. The mechanism involves the formation of a carboxylate mixed formimidic anhydride intermediate via a concerted mechanism. This intermediate is then transformed to an N-formylamide by a concerted pseudopericyclic [1,3]-acyl shift. Mechanisms involving zwitterions or diradicals are discounted.     

Peptide hydrolysis by the binuclear zinc enzyme aminopeptidase from Aeromonas proteolytica: a density functional theory study

Aminopeptidase from Aeromonas proteolytica (AAP) is a binuclear zinc enzyme that catalyzes the cleavage of the N-terminal amino acid residue of peptides and proteins. In this study, we used density functional methods to investigate the reaction mechanism of this enzyme. A model of the active site was constructed on the basis of the X-ray crystal structure of the native enzyme, and a model dipeptide was used as a substrate. It was concluded that the hydroxide is capable of performing a nucleophilic attack at the peptide carbonyl from its bridging position without the need to first become terminal. The two zinc ions are shown to have quite different roles. Zn2 binds the amino group of the substrate, thereby orienting it toward the nucleophile, while Zn1 stabilizes the alkoxide ion of the tetrahedral intermediate, thereby lowering the barrier for the nucleophilic attack. The rate-limiting step is suggested to be the protonation of the nitrogen of the former peptide bond, which eventually leads to the cleavage of the C-N bond.     

Theoretical study of the mechanism of oxoiron(IV) formation from H2O2 and a nonheme iron(II) complex: O-O cleavage involving proton-coupled electron transfer

It has recently been shown that the nonheme oxoiron(IV) species supported by the 1,4,8,11-tetramethyl-1,4,8,11-tetraazacyclotetradecane ligand (TMC) can be generated in near-quantitative yield by reacting [Fe(II)(TMC)(OTf)(2)] with a stoichiometric amount of H(2)O(2) in CH(3)CN in the presence of 2,6-lutidine (Li, F.; England, J.; Que, L., Jr. J. Am. Chem. Soc. 2010, 132, 2134-2135). This finding has major implications for O-O bond cleavage events in both Fenton chemistry and nonheme iron enzymes. To understand the mechanism of this process, especially the intimate details of the O-O bond cleavage step, a series of density functional theory (DFT) calculations and analyses have been carried out. Two distinct reaction paths (A and B) were identified. Path A consists of two principal steps: (1) coordination of H(2)O(2) to Fe(II) and (2) a combination of partial homolytic O-O bond cleavage and proton-coupled electron transfer (PCET). The latter combination renders the rate-limiting O-O cleavage effectively a heterolytic process. Path B proceeds via a simultaneous homolytic O-O bond cleavage of H(2)O(2) and Fe-O bond formation. This is followed by H abstraction from the resultant Fe(III)-OH species by an ?OH radical. Calculations suggest that path B is plausible in the absence of base. However, once 2,6-lutidine is added to the reacting system, the reaction barrier is lowered and more importantly the mechanistic path switches to path A, where 2,6-lutidine plays an essential role as an acid-base catalyst in a manner similar to how the distal histidine or glutamate residue assists in compound I formation in heme peroxidases. The reaction was found to proceed predominantly on the quintet spin state surface, and a transition to the triplet state, the experimentally known ground state for the TMC-oxoiron(IV) species, occurs in the last stage of the oxoiron(IV) formation process.     

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.     

Unusual concerted Lewis acid-Lewis base mechanism for hydrogen activation by a phosphine-borane compound

Density functional theory calculations have been carried out to investigate the possible reaction mechanisms for the reversible dihydrogen activation catalyzed by a phosphine-borane compound, (C6H2Me3)2P(C6F4)B(C6F5)2 (Welch, G. C.; Juan, R. R. S; Masuda, J. D.; Stephan, D. W. Science. 2006, 314, 1124-1126). The present calculations show that an unusual concerted Lewis acid-Lewis base mechanism is more favorable than the proton transfer or hydride transfer mechanisms suggested previously. In the concerted Lewis acid-Lewis base mechanism, the H-H heterolytic cleavage is achieved through the simultaneous electron transfer from the lone-pair orbital of the Lewis base P center to the sigma* orbital of H2 and from the sigma orbital of H2 to the empty orbital of the Lewis acid B center. The solvent is found to dramatically change the potential energy surface. The proposed mechanism can account well for the bimolecular H-D exchange process observed in deuteration experiments and the experimental fact that the H2 activation is reversible at mild conditions.     

Hydrogen isotope effects and mechanism of aqueous ozone and peroxone decompositions

Hydrogen peroxide exalts the reactivity of aqueous ozone by reasons that remain obscure. Should H2O2 enhance free radical production, as it is generally believed, a chain mechanism propagated by (.OH/.O2-) species would account for O3 decomposition rates in neat H2O, HR-O3, and in peroxone (O3 + H2O2) solutions, HPR-O3. We found, however, that: (1) the radical mechanism correctly predicts HR-O3 but vastly overestimates HPR-O3, (2) solvent deuteration experiments preclude radical products from the (O3 + HO2-) reaction. The modest kinetic isotope effect (KIE) we measure in H2O/D2O: HR-O3/DR-O3 = 1.5 +/- 0.3, is compatible with a chain process driven by electron- and/or O-atom transfer processes. But the large KIE found in peroxone: HPR-O3/DPR-O3 = 19.6 +/- 4.0, is due to an elementary (O3 + HO2-) reaction involving H-O2- bond cleavage. Since the KIE for the hypothetical H-atom transfer: O3 + HO2- HO3. +.O2-, would emerge as a KIE1/2 factor in the rates of the ensuing radical chain, the magnitude of the observed KIE must be associated with the hydride transfer reaction that yields a diamagnetic species: O3 + HO2- HO3- + O2. HO3-/H2O3 may be the bactericidal trioxide recently identified in the antibody-catalyzed addition of O2(1Deltag) to H2O.     

Facile S(N)2 reaction in protic solvent: quantum chemical analysis

We study the effects of protic solvent (water, methanol, ethanol, and tert-butyl alcohol) and cation (Na+, K+, Cs+) on the unsymmetrical SN2 reaction X- + RY --> RX + Y- (X = F, Br; R = CH3,C3H7;Y = Cl, OMs). We describe a series of calculations for the S(N)2 reaction mechanism under the influence of cation and protic solvent, presenting the structures of pre- and postreaction complexes and transition states and the magnitude of the activation barrier. An interesting mechanism is proposed, in which the protic solvent molecules that are shielded from the nucleophile by the intervening cation act as a Lewis base to reduce the unfavorable Coulombic influence of the cation on the nucleophile. We predict that the reaction barrier for the S(N)2 reaction is significantly lowered by the cooperative effects of cation and protic solvent. We show that the cation and protic solvent, each of which has been considered to retard the SN2 reactivity of the nucleophile, can accelerate the reaction tremendously when they interact with the fluoride ion in an intricate, combined fashion. This alternative S(N)2 mechanism is discussed in relation to the recently observed phenomenal efficiency of fluorination in tert-alcohol media [Kim, D. W.; et al. J. Am. Chem. Soc. 2006, 128, 16394].