The effect of heme environment on the hydrogen abstraction reaction of camphor in P450cam catalysis: a QM/MM study

The discrepancies between the published QM/MM studies (Sch?neboom, J. C.; Cohen, S.; Lin, H.; Shaik, S.; Thiel, W. J. Am. Chem. Soc. 2004, 126, 4017; Guallar, V.; Friesner, R. A. J. Am. Chem. Soc. 2004, 126, 8501) on H-abstraction of camphor in P450cam have largely been resolved. The crystallographic water molecule 903 situated near the oxo atom of Compound I acts as a catalyst for H-abstraction, lowering the barrier by about 4 kcal/mol. Spin density at the A-propionate side chain of heme can occur in the case of incomplete screening but has no major effect on the computed barrier.     

Structural studies of nonclassical cyclobutylmethyl cations by the ab initio method

Ab initio calculations at the MP2/cc-pVTZ level show that the cyclobutylmethyl cation is a nonclassical sigma-delocalized species, which is distinct from the global minimum C2-symmetric cyclopentyl cation (Schleyer, P. v. R.; Carneiro, J. W. de M.; Koch, W.; Raghavachari, K. J. Am. Chem. Soc. 1989, 111, 5475). Relatively lower level DFT calculations, on the other hand, show that the primary cyclobutylmethyl cation spontaneously collapses into the cyclopentyl cation (Prakash, G. K. S.; Reddy, V. P.; Rasul, G.; Casanova, J.; Olah, G. A. J. Am. Chem. Soc. 1998, 120, 13362). The secondary 1-cyclobutylethyl cation is also a nonclassical carbocation, as shown by calculations at the MP2/cc-pVTZ level. Two structures having energy minima are identified for the latter cation on the potential energy surface. The conformer in which the methyl group is in the exo orientation is a global minimum and is favored over the corresponding endo conformer by 1.2 kcal/mol at the MP2/cc-pVTZ//MP2/cc-pVTZ +ZPE level of calculations. The tertiary 1-cyclobutyl-1-methylethyl cation, at this level of calculations, also involves substantial nonclassical sigma-delocalization, showing that the nonclassical stabilization is more important for cyclobutylmethyl cations relative to the cyclopropylmethyl cations. The 13C NMR chemical shifts obtained from GIAO-CCSD(T)/tzp/dz calculations further substantiate the nonclassical structures for these carbocations.     

Radical clock substrates, their C-H hydroxylation mechanism by cytochrome P450, and other reactivity patterns: what does theory reveal about the clocks' behavior?

There is an ongoing and tantalizing controversy regarding the mechanism of a key process in nature, C-H hydroxylation, by the enzyme cytochrome P450 (Auclaire, K.; Hu, Z.; Little, D. M.; Ortiz de Montellano, P. R.; Groves, J. T. J. Am. Chem. Soc. 2002, 124, 6020-6027. Newcomb, M.; Aebisher, D.; Shen, R.; Esala, R.; Chandrasena, P.; Hollenberg, P. F.; Coon, M. J. J. Am. Chem. Soc. 2003, 125, 6064-6065). To definitely resolve this controversy, theory must first address the actual systems that have been used by experiment, and that generated the controversy. This is done in the present paper, which constitutes the first extensive theoretical study of such two experimental systems, trans-2-phenylmethyl-cyclopropane (1) and trans-2-phenyl-iso-propylcyclopropane (4). The theoretical study of these substrates reveals that the only low energy pathway for C-H hydroxylation is the two-state rebound mechanism described originally for methane hydroxylation (Ogliaro, F.; Harris, N.; Cohen, S.; Filatov, M.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 8977-8989). The paper shows that the scenario of a two-state rebound mechanism accommodates much of the experimental data. The computational results provide a good match to experimental results concerning the very different extents of rearrangement for 1 (20-30%) vs 4 (virtually none), lead to product isotope effect for the reaction of 1, in the direction of the experimental result, and predict as well the observed metabolic switching from methyl to phenyl hydroxylation, which occurs upon deuteration of the methyl group. Furthermore, the study reveals that an intimate ion pair species involving an alkyl carbocation derived from 4 gives no rearranged products, again in accord with experiment. This coherent match between theory and experiment cannot be merely accidental; it comes close to being aproof that the actual mechanism of C-H hydroxylation involves the two-state reactivity revealed by theory. Analysis of the rearrangement modes of the carbocations derived from 1 and 4 excludes the participation of free carbocations during the hydroxylation of these substrates. Finally, the mechanistic significance of product isotope effect (different isotope effects for the rearranged and unrearranged alcohol products) is analyzed. It is shown to be a sensitive probe of two-state reactivity; the size of the intrinsic product isotope effect and its direction reveal the structural differences of the hydrogen abstraction transition states in the low-spin vs high-spin reaction manifolds.     

O-H bond dissociation enthalpies in oximes: order restored

The O-H bond dissociation enthalpies (BDEs) of 13 oximes, RR'C=NOH, having R and/or R' = H, alkyl, and aryl are reported. Experimental anchor points used to validate the results of theoretical calculations include (1) the O-H BDEs of (t-Bu)2C=NOH, t-Bu(i-Pr)C=NOH, and t-Bu(1-Ad)C=NOH determined earlier from the heat released in the reaction of (t-Bu)2C=NO* with (PhNH)2 in benzene and EPR spectroscopy (Mahoney, L. R.; Mendenhall, G. D.; Ingold, K. U. J. Am. Chem. Soc. 1973, 95, 8610), all of which were decreased by 1.7 kcal/mol to reflect a revision to the heat of formation of (E)-azobenzene (which has significant ramifications for other BDEs) and to correct for the heat of hydrogen bonding of (t-Bu)2C=NOH (alphaH2 = 0.43 measured in this work) to benzene, and (2) the measured rates of thermal decomposition of six RR'C=NOCH2Ph at 423 or 443 K, which were used to derive O-H BDEs for the corresponding RR'C=NOH. Claims (Bordwell, F. G.; Ji, G. Z. J. Org. Chem. 1992, 57, 3019; Bordwell, F. G.; Zhang, S. J. Am. Chem. Soc. 1995, 117, 4858; and Bordwell, F. G.; Liu, W.-Z. J. Am. Chem. Soc. 1996, 118, 10819) that the O-H BDEs in mono- and diaryloximes are significantly lower than those for alkyloximes due to delocalization of the unpaired electron into the aromatic ring have always been inconsistent with the known structures of iminoxyl radicals as are the purported perpendicular structures, i.e., phi(Calpha-C=N-O*) = 90 degrees, for sterically hindered dialkyl iminoxyl radicals. The present results confirm the 1973 conclusion that simple steric effects, not electron delocalization or dramatic geometric changes, are responsible for the rather small differences in oxime O-H BDEs.     

Dinuclear Zn(II) complex catalyzed phosphodiester cleavage proceeds via a concerted mechanism: a density functional theory study

Density functional theory (DFT) calculations were used to study the mechanism for the cleavage reaction of the RNA analogue HpPNP (HpPNP = 2-hydroxypropyl-4-nitrophenyl phosphate) catalyzed by the dinuclear Zn(II) complex of 1,3-bis(1,4,7-triazacyclonon-1-yl)-2-hydroxypropane (Zn(2)(L(2)O)). We present a binding mode in which each terminal phosphoryl oxygen atom binds to one zinc center, respectively, and the nucleophilic 2-hydroxypropyl group coordinates to one of the zinc ions, while the hydroxide from deprotonation of a water molecule coordinates to the other zinc ion. Our calculations found a concerted mechanism for the HpPNP cleavage with a 16.5 kcal/mol reaction barrier. An alternative proposed stepwise mechanism through a pentavalent oxyphosphorane dianion reaction intermediate for the HpPNP cleavage was found to be less feasible with a significantly higher energy barrier. In this stepwise mechanism, the deprotonation of the nucleophilic 2-hydroxypropyl group is accompanied with nucleophilic attack in the rate-determining step. Calculations of the nucleophile (18)O kinetic isotope effect (KIE) and leaving (18)O KIE for the concerted mechanism are in reasonably good agreement with the experimental values. Our results indicate a specific-base catalysis mechanism takes place in which the deprotonation of the nucleophilic 2-hydroxypropyl group occurs in a pre-equilibrium step followed by a nucleophilic attack on the phosphorus center. Detailed comparison of the geometric and electronic structure for the HpPNP cleavage reaction mechanisms in the presence/absence of catalyst revealed that the catalyst significantly altered the determining-step transition state to become far more associative or tight, that is, bond formation to the nucleophile was remarkably more advanced than leaving group bond fission in the catalyzed mechanism. Our results are consistent with and provide a reliable interpretation for the experimental observations that suggest the reaction occurs by a concerted mechanism (see Humphry, T.; Iyer, S.; Iranzo, O.; Morrow, J. R.; Richard, J. P.; Paneth, P.; Hengge, A. C. J. Am. Chem. Soc. 2008, 130, 17858-17866) and has a specific-base catalysis character (see Yang, M.-Y.; Iranzo, O.; Richard, J. P.; Morrow, J. R. J. Am. Chem. Soc. 2005, 127, 1064-1065).     

Determinants of cyanuric acid and melamine assembly in water

While the recognition of cyanuric acid (CA) by melamine (M) and their derivatives has been known to occur in both water and organic solvents for some time, analysis of CA/M assembly in water has not been reported (Ranganathan, A.; Pedireddi, V. R.; Rao, C. N. R. J. Am. Chem. Soc.1999, 121, 1752-1753; Mathias, J. P.; Simanek, E. E.; Seto, C. T.; Whitesides, G. M. Macromol. Symp.1994, 77, 157-166; Zerkowski, J. A.; MacDonald, J. C.; Seto, C. T.; Wierda, D. A.; Whitesides, G. M. J. Am. Chem. Soc.1994, 116, 2382-2391; Mathias, J. P.; Seto, C. T.; Whitesides, G. M. Polym. Prepr.1993, 34, 92-93; Seto, C. T.; Whitesides, G. M. J. Am. Chem. Soc.1993, 115, 905-916; Zerkowski, J. A.; Seto, C. T.; Whitesides, G. M. J. Am. Chem. Soc.1992, 114, 5473-5475; Seto, C. T.; Whitesides, G. M. J. Am. Chem. Soc.1990, 112, 6409-6411; Wang, Y.; Wei, B.; Wang, Q. J. Chem. Cryst.1990, 20, 79-84; ten Cate, M. G. J.; Huskens, J.; Crego-Calama, M.; Reinhoudt, D. N. Chem.-Eur. J.2004, 10, 3632-3639). We have examined assembly of CA/M, as well as assembly of soluble trivalent CA and M derivatives (TCA/TM), in aqueous solvent, using a combination of solution phase NMR, isothermal titration and differential scanning calorimetry (ITC/DSC), cryo-transmission electron microscopy (cryo-TEM), and synthetic chemistry. While the parent heterocycles coprecipitate in water, the trivalent system displays more controlled and cooperative assembly that occurs at lower concentrations than the parent and yields a stable nanoparticle suspension. The assembly of both parent and trivalent systems is rigorously 1:1 and proceeds as an exothermic, proton-transfer coupled process in neutral pH water. Though CA and M are considered canonical hydrogen-bonding motifs in organic solvents, we find that their assembly in water is driven in large part by enthalpically favorable surface-area burial, similar to what is observed with nucleic acid recognition. There are currently few synthetic systems capable of robust molecular recognition in water that do not rely on native recognition motifs, possibly due to an incomplete understanding of recognition processes in water. This study establishes a detailed conceptual framework for considering CA/M heterocycle recognition in water which enables the future design of molecular recognition systems that function in water.     

Product distribution of peroxynitrite decay as a function of pH, temperature, and concentration.

The decay of peroxynitrite [O=NOO(-), oxoperoxonitrate(1-)] was examined as a function of concentration (0.050-2.5 mM), temperature (5-45 degrees C), and pH (2.2-10.0). Below 5 degrees C and pH 7, little amounts of the decomposition products nitrite and dioxygen are formed, even when the peroxynitrite concentration is high (2.5 mM). Instead, approximately > or =90% isomerizes to nitrate. At higher pH, decomposition increases at the expense of isomerization, up to nearly 80% at pH 10.0 at 5 degrees C and 90% at 45 degrees C. Much less nitrite and dioxygen per peroxynitrite are formed when the peroxynitrite concentration is lower; at 50 microM and pH 10.2, < or =40% decomposes. In contrast to two other reports (Pfeiffer, S.; Gorren, A. C. F.; Schmidt, K.; Werner, E. R.; Hansert, B.; Bohle, D. S.; Mayer, B. J. Biol. Chem. 1997, 272, 3465-3470, and Coddington, J. W.; Hurst, J. K.; Lymar, S. V. J. Am. Chem. Soc. 1999, 121, 2438-2443), we find that the extent of decomposition is dependent on the peroxynitrite concentration.     

A proton-shuttle mechanism mediated by the porphyrin in benzene hydroxylation by cytochrome p450 enzymes

Benzene hydroxylation is a fundamental process in chemical catalysis. In nature, this reaction is catalyzed by the enzyme cytochrome P450 via oxygen transfer in a still debated mechanism of considerable complexity. The paper uses hybrid density functional calculations to elucidate the mechanisms by which benzene is converted to phenol, benzene oxide, and ketone, by the active species of the enzyme, the high-valent iron-oxo porphyrin species. The effects of the protein polarity and hydrogen-bonding donation to the active species are mimicked, as before (Ogliaro, F.; Cohen, S.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 12892-12893). It is verified that the reaction does not proceed either by hydrogen abstraction or by initial electron transfer (Ortiz de Montellano, P. R. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; Chapter 8, pp 245-303). In accord with the latest experimental conclusions, the theoretical calculations show that the reactivity is an interplay of electrophilic and radicalar pathways, which involve an initial attack on the pi-system of the benzene to produce sigma-complexes (Korzekwa, K. R.; Swinney, D. C.; Trager, W. T. Biochemistry 1989, 28, 9019-9027). The dominant reaction channel is electrophilic and proceeds via the cationic sigma-complex,( 2)3, that involves an internal ion pair made from a cationic benzene moiety and an anionic iron porphyrin. The minor channel proceeds by intermediacy of the radical sigma-complex, (2)2, in which the benzene moiety is radicalar and the iron-porphyrin moiety is neutral. Ring closure in these intermediates produces the benzene oxide product ((2)4), which does not rearrange to phenol ((2)7) or cyclohexenone ((2)6). While such a rearrangement can occur post-enzymatically under physiological conditions by acid catalysis, the computations reveal a novel mechanism whereby the active species of the enzyme catalyzes directly the production of phenol and cyclohexenone. This enzymatic mechanism involves proton shuttles mediated by the porphyrin ring through the N-protonated intermediate, (2)5, which relays the proton either to the oxygen atom to form phenol ((2)7) or to the ortho-carbon atom to produce cyclohexenone product ((2)6). The formation of the phenol via this proton-shuttle mechanism will be competitive with the nonenzymatic conversion of benzene oxide to phenol by external acid catalysis. With the assumption that (2)5 is not fully thermalized, this novel mechanism would account also for the observation that there is a partial skeletal retention of the original hydrogen of the activated C-H bond, due to migration of the hydrogen from the site of hydroxylation to the adjacent carbon (so-called "NIH shift" (Jerina, D. M.; Daly, J. W. Science 1974, 185, 573-582)). Thus, in general, the computationally discovered mechanism of a porphyrin proton shuttle suggests thatthere is an enzymatic pathway that converts benzene directly to a phenol and ketone, in addition to nonenzymatic production of these species by conversion of arene oxide to phenol and ketone. The potential generality of protonated porphyrin intermediates in P450 chemistry is discussed in the light of the H/D exchange observed during some olefin epoxidation reactions (Groves, J. T.; Avaria-Neisser, G. E.; Fish, K. M.; Imachi, M.; Kuczkowski, R. J. Am. Chem. Soc. 1986, 108, 3837-3838) and the general observation of heme alkylation products (Kunze, K. L.; Mangold, B. L. K.; Wheeler, C.; Beilan, H. S.; Ortiz de Montellano, P. R. J. Biol. Chem. 1983, 258, 4202-4207). The competition, similarities, and differences between benzene oxidation viz. olefin epoxidation and alkanyl C-H hydroxylation are discussed, and comparison is made with relevant experimental and computational data. The dominance of low-spin reactivity in benzene hydroxylation viz. two-state reactivity (Shaik, S.; de Visser, S. P.; Ogliaro, F.; Schwarz, H.; Schr?der, D. Curr. Opin. Chem. Biol. 2002, 6, 556-567) in olefin epoxidation and alkane hydroxylation is traced to the loss of benzene resonance energy during the bond activation step.     

Controlled excitations of the Belousov-Zhabotinsky reaction: Experimental procedures

The purpose of this research was to explore the unstirred, ferroin-catalyzed Belousov-Zhabotinsky (BZ) reaction as an experimental model for the response of excitable media to small perturbations (slightly larger than the threshold for excitations). Following Showalter et al. (Showalter, K.; Noyes, R. M.; Turner, H. J.Am. Chem. Soc. 1979, 101, 7463-69), we used a positively biased silver electrode to release silver ions into a BZ reaction mixture, removing bromide ions and causing an excitation if sufficient bromide was removed. We found (1) a scaling region in which the delay before activation increased linearly as the size of the perturbation decreased, qualitatively consistent with but not fully explained by the Oregonator of Field et al. (Field, R. J.; K?r?s, E.; Noyes, R. M. J. Am. Chem. Soc. 1972, 94, 8649-64); (2) evidence for a 10 s oligomerization time scale; and (3) that activations were always delayed until after the end of a pulse of current, with the delay essentially constant for sufficiently long pulses, an effect not seen in simple ODE models but consistent with the anomalously large current apparently required for activation (Showalter, K.; Noyes, R. M. J. Am. Chem. Soc. 1976, 98, 3730-31) and explainable by bromide transport. Overall, the BZ system appeared to be well-suited as an experimental prototype, despite its complexity.     

Can ferric-superoxide act as a potential oxidant in P450(cam)? QM/MM investigation of hydroxylation, epoxidation, and sulfoxidation

In view of recent reports of high reactivity of ferric-superoxide species in heme and nonheme systems (Morokuma et al. J. Am. Chem. Soc. 2010, 132, 11993-12005; Que et al. Inorg. Chem. 2010, 49, 3618-3628; Nam et al. J. Am. Chem. Soc. 2010, 132, 5958-5959; J. Am. Chem. Soc. 2010, 132, 10668-10670), we use herein combined quantum mechanics/molecular mechanics (QM/MM) methods to explore the potential reactivity of P450(cam) ferric-superoxide toward hydroxylation, epoxidation, and sulfoxidation. The calculations demonstrate that P450 ferric-superoxide is a sluggish oxidant compared with the high-valent oxoiron porphyrin cation-radical species. As such, unlike heme enzymes with a histidine axial ligand, the P450 superoxo species does not function as an oxidant in P450(cam). The origin of this different behavior of the superoxo species of P450 vis-a?-vis other heme enzymes like tryptophan 2, 3-dioxygenase (TDO) is traced to the ability of the latter superoxo species to make a stronger FeOO-X (X = H,C) bond and to stabilize the corresponding bond-activation transition states by resonance with charge-transfer configurations. By contrast, the negatively charged thiolate ligand in the P450 superoxo species minimizes the mixing of charge transfer configurations in the transition state and raises the reaction barrier. However, as we demonstrate, an external electric field oriented along the Fe-O axis with a direction pointing from Fe toward O will quench Cpd I formation by slowing the reduction of ferric-superoxide and will simultaneously lower the barriers for oxidation by the latter species, thereby enabling observation of superoxo chemistry in P450. Other options for nascent superoxo reactivity in P450 are discussed.     

Do spin state changes matter in organometallic chemistry? A computational study

Spin changes occur often in organometallic chemistry, and their effect on kinetics is not well understood. We report computations on the singlet and triplet potential energy surfaces of several processes of this type and show that the topology of the individual surfaces, as well as of the crossing regions between them, can be used to rationalize the observed reactivity in all cases. In particular, the slow addition of dihydrogen to W[N(CH(2)CH(2)NSiMe(3))(3)]H (Schrock, R. R.; Shih, K. Y.; Dobbs, D. A.; Davis, W. M. J. Am. Chem. Soc. 1995, 117, 6609) is shown to be a "spin-blocked" reaction with a high barrier due to the crossing between reactant triplet and product singlet surfaces. In contrast, addition of CO to TpCo(CO) (Detrich, J. L.; Reinaud, O. M.; Rheingold, A. L.; Theopold, K. H. J. Am. Chem. Soc. 1995, 117, 11745) is fast because the triplet and singlet surfaces cross at low energy. Particular care is taken to use DFT methods which are in adequate agreement with experimental and high-level computational energetics for these systems.     

Origin of the nonplanarity of tetrafluoro cyclobutadiene, C4F4

Density functional theory as well as M?ller-Plesset investigations has been carried out on tetrafluoro cyclobutadiene, C4F4, to explore the origin of its nonplanarity. Although Petersson et al. (Petersson, E. J.; Fanuele, J. C.; Nimlos, M. R.; Lemal, D. M.; Ellison, G. B.; Radziszewski, J. G. J. Am. Chem. Soc. 1997, 119, 11122-11123) had earlier predicted a nonplanar geometry of this compound on the basis of spectral and bond orbital analysis, the explanation of the same from a more fundamental point of view is still missing. In the present study, we provide a heuristic explanation for the origin of nonplanarity of C4F4. The two major driving forces behind this nonplanar geometry are the unusual aromaticity of this cyclic homoatomic 4pi electron system and the second-order Jahn-Teller effect (SOJTE). These driving forces can well be explained by various energy and density parameters and also by nucleus-independent chemical shift (NICS) values. Aromaticity of a cyclic homoatomic 4pi electron system is quite remarkable. The enhancement of pi- delocalization as evidenced from molecular orbital analysis may be attributed to s-ppi mixing in nonplanar C4F4.     

Synthesis and properties of a fluorene-capped isotruxene: a new unsymmetrical star-shaped pi-system

An unsymmetrical star-shaped pi-system (1) with an isotruxene core and three fluorene arms has been synthesized, and its photophysical, electrochemical, and thermochemical properties have been investigated and compared with the corresponding symmetrical truxene derivatives 2a-d (Kanibolotsky, A. L.; Berridge, R.; Skabara, P. J.; Perepichka, I. F.; Bradley, D. D. C.; Koeberg, M. J. Am. Chem. Soc. 2004, 126, 13695-13702). Stronger electronic couplings between the arms in 1 vs 2a-d lead to lower optical band gap, larger splitting in oxidation potentials, and superior thermostability. [structure: see text]     

Capturing a synergistic effect of a conical push and an inward pull in fluoro derivatives of Li@B10H14 basket: toward a higher vertical ionization potential and nonlinear optical response

We present the design of fluoro derivatives of B(10)H(14) and Li@B(10)H(14) baskets. A synergistic effect of conical push and inward pull (reported independently in previous lithium nonlinear optical (NLO) complexes) has been explored in these derivatives to achieve a robustly large NLO response and a higher vertical ionization potential. Li@1,3,6,9-F(4)B(10)H(10), Li@6,9-F(2)B(10)H(12), and Li@2,4,6,9-F(4)B(10)H(10) exhibit first hyperpolarizability (β(0)) values as large as 181?624, 133?199, and 32?314 au; their vertical ionization potentials are 6.45, 6.30, and 6.78 eV, respectively. These values are significantly higher than those previously reported in Li-doped fluorocarbon chains at the same MP2/6-31+G* level of theory (Xu, H. L.; Li, Z. R.; Wu, D.; Wang, B. Q.; Li, Y.; Gu, F. L.; Aoki, Y. J. Am. Chem. Soc. 2007, 129, 2967). They also exceed those from our earlier designed Li@B(10)H(14) basket (Muhammad, S.; Xu, H. L.; Liao, Y.; Kan, Y. H.; Su , Z. M. J. Am. Chem. Soc. 2009, 131, 2967). In addition, new quantum chemical calculations of enthalpies of reaction (Δ(r)H°) at 298 K for B(10)H(14) and its lithium/fluoro derivatives highlight the changes in their thermodynamical aspects. The calculated enthalpies of lithiation reactions are -10.04, -11.29, and -13.18 kcal/mol for B(10)H(14), 6,9-F(2)B(10)H(12), and 2,4-F(2)B(10)H(12), respectively, demonstrating a higher probability of fluoro decaboranes for reaction with lithium. The obtained results not only explain the effect of position and number dependence of substituted fluoro atom(s) in B(10)H(14) and Li@B(10)H(14) but also elucidate a synergistic behavior to polarize a lithium excess electron for high NLO responses and vertical ionization potentials.     

Searching for the second oxidant in the catalytic cycle of cytochrome P450: a theoretical investigation of the iron(III)-hydroperoxo species and its epoxidation pathways

Iron(III)-hydroperoxo, [Por(CysS)Fe(III)-OOH](-), a key species in the catalytic cycle of cytochrome P450, was recently identified by EPR/ENDOR spectroscopies (Davydov, R.; Makris, T. M.; Kofman, V.; Werst, D. E.; Sligar, S. G.; Hoffman, B. M. J. Am. Chem. Soc. 2001, 123, 1403-1415). It constitutes the last station of the preparative steps of the enzyme before oxidation of an organic compound and is implicated as the second oxidant capable of olefin epoxidation (Vaz, A. D. N.; McGinnity, D. F.; Coon, M. J. Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 3555-3560), in addition to the penultimate active species, Compound I (Groves, J. T.; Han, Y.-Z. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; pp 3-48). In response, we present a density functional study of a model species and its ethylene epoxidation pathways. The study characterizes a variety of properties of iron(III)-hydroperoxo, such as the O-O bonding, the Fe-S bonding, Fe-O and Fe-S stretching frequencies, its electron attachment, and ionization energies. Wherever possible these properties are compared with those of Compound I. The proton affinities for protonation on the proximal and distal oxygen atoms of iron(III)-hydroperoxo, and the effect of the thiolate ligand thereof, are determined. In accordance with previous results (Harris, D. L.; Loew, G. H. J. Am. Chem. Soc. 1998, 120, 8941-8948), iron(III)-hydroperoxo is a strong base (as compared with water), and its distal protonation leads to a barrier-free formation of Compound I. The origins of this barrier-free process are discussed using a valence bond approach. It is shown that the presence of the thiolate is essential for this process, in line with the "push effect" deduced by experimentalists (Sono, M.; Roach, M. P.; Coulter, E. D.; Dawson, J. H. Chem. Rev. 1996, 96, 2841-2887). Finally, four epoxidation pathways of iron(III)-hydroxperoxo are located, in which the species transfers oxygen to ethylene either from the proximal or from the distal sites, in both concerted and stepwise manners. The barriers for the four mechanisms are 37-53 kcal/mol, in comparison with 14 kcal/mol for epoxidation by Compound I. It is therefore concluded that iron(III)-hydroperoxo, as such, cannot be a second oxidant, in line with its significant basicity and poor electron-accepting capability. Possible versions of a second oxidant are discussed.     

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.     

A coupled reference interaction site model/molecular dynamics study of the potential of mean force curve of the SN2 Cl- + CH3Cl reaction in water

An application of the coupled reference interaction site model (RISM)/simulation methodology to the calculation of the potential of mean force (PMF) curve in aqueous solution for the identity nucleophilic substitution reaction Cl(-) + CH(3)Cl is performed. The free energy of activation is calculated to be 27.1 kcal/mol which compares very well with the experimentally determined barrier height of 26.6 kcal/mol. Furthermore, the calculated PMF is almost superimposed with that previously calculated using the computationally rigorous Monte Carlo with importance sampling method (Chandrasekhar, J.; Smith, S. F.; Jorgensen, W. L. J. Am. Chem. Soc. 1985, 107, 154). Using the calculated PMF, a crude estimate of the solvated kinetic transmission coefficient also compares well with that of previous more accurate simulations. These results indicate that the coupled RISM/simulation method provides a cost-effective methodology for studying reactions in solution.     

Conserved mechanism of copper binding and transfer. A comparison of the copper-resistance proteins PcoC from Escherichia coli and CopC from Pseudomonas syringae

The copper-resistance proteins PcoC from Escherichia coli and CopC from Pseudomonas syringae exhibit 67% sequence identity, but the chemistry reported for PcoC (Peariso, K.; Huffman, D. L.; Penner-Hahn, J. E.; O'Halloran, T. V. J. Am. Chem. Soc. 2003, 125, 342-343) was distinctly different from that reported for CopC (Zhang, L.; Koay, M.; Maher, M. J.; Xiao, Z.; Wedd, A. G. J. Am. Chem. Soc. 2006, 128, 5834-5850). The source of the inconsistency has been identified, and His1 is confirmed as an unprecedented bidentate ligand in each protein. Access to a bona fide wild-type PcoC protein allowed unequivocal observation of intermediates involved in intermolecular redox copper transfer reactions.     

Characteristics of a miniature compartment-less glucose-O2 biofuel cell and its operation in a living plant

We report the temperature, pH, glucose concentration, NaCl concentration, and operating atmosphere dependence of the power output of a compartment-less miniature glucose-O(2) biofuel cell, comprised only of two bioelectrocatalyst-coated carbon fibers, each of 7 micro m diameter and 2 cm length (Mano, N.; Mao, F.; Heller, A. J. Am. Chem. Soc. 2002, 124, 12962). The bioelectrocatalyst of the anode consists of glucose oxidase from Aspergillus niger electrically "wired" by polymer I, having a redox potential of -0.19 V vs Ag/AgCl. That of the cathode consists of bilirubin oxidase from Trachyderma tsunodae "wired" by polymer II having a redox potential of +0.36 V vs Ag/AgCl (Mano, N.; Kim, H.-H.; Zhang, Y.; Heller, A. J. Am. Chem. Soc. 2002, 124, 6480. Mano, N.; Kim, H.-H.; Heller, A. J. Phys. Chem. B 2002, 106, 8842). Implantation of the fibers in the grape leads to an operating biofuel cell producing 2.4 micro W at 0.52 V.