Catalytic Mechanism of a Heme and Tyrosyl Radical-Containing Fatty Acid α-(Di)oxygenase

The steady-state catalytic mechanism of a fatty acid α-(di)oxygenase is examined, revealing that a persistent tyrosyl radical (Tyr379(?)) effects O(2) insertion into C(α)-H bonds of fatty acids. The initiating C(α)-H homolysis step is characterized by apparent rate constants and deuterium kinetic isotope effects (KIEs) that increase hyperbolically upon raising the concentration of O(2). These results are consistent with H(?) tunneling, transitioning from a reversible to an irreversible regime. The limiting deuterium KIEs increase from ～30 to 120 as the fatty acid chain is shortened from that of the native substrate. In addition, activation barriers increase in a manner that reflects decreased fatty acid binding affinities. Anaerobic isotope exchange experiments provide compelling evidence that Tyr379(?) initiates catalysis by H(?) abstraction. C(α)-H homolysis is kinetically driven by O(2) trapping of the α-carbon radical and reduction of a putative peroxyl radical intermediate to a 2(R)-hydroperoxide product. These findings add to a body of work which establishes large-scale hydrogen tunneling in proteins. This particular example is novel because it involves a protein-derived amino acid radical.     

Kinetic effects of increased proton transfer distance on proton-coupled oxidations of phenol-amines

To test the effect of varying the proton donor-acceptor distance in proton-coupled electron transfer (PCET) reactions, the oxidation of a bicyclic amino-indanol (2) is compared with that of a closely related phenol with an ortho CPh(2)NH(2) substituent (1). Spectroscopic, structural, thermochemical, and computational studies show that the two amino-phenols are very similar, except that the O···N distance (d(ON)) is >0.1 ? longer in 2 than in 1. The difference in d(ON) is 0.13 ± 0.03 ? from X-ray crystallography and 0.165 ? from DFT calculations. Oxidations of these phenols by outer-sphere oxidants yield distonic radical cations (?)OAr-NH(3)(+) by concerted proton-electron transfer (CPET). Simple tunneling and classical kinetic models both predict that the longer donor-acceptor distance in 2 should lead to slower reactions, by ca. 2 orders of magnitude, as well as larger H/D kinetic isotope effects (KIEs). However, kinetic studies show that the compound with the longer proton-transfer distance, 2, exhibits smaller KIEs and has rate constants that are quite close to those of 1. For example, the oxidation of 2 by the triarylamminium radical cation N(C(6)H(4)OMe)(3)(?+) (3a(+)) occurs at (1.4 ± 0.1) × 10(4) M(-1) s(-1), only a factor of 2 slower than the closely related reaction of 1 with N(C(6)H(4)OMe)(2)(C(6)H(4)Br)(?+) (3b(+)). This difference in rate constants is well accounted for by the slightly different free energies of reaction: ΔG° (2 + 3a(+)) = +0.078 V versus ΔG° (1 + 3b(+)) = +0.04 V. The two phenol-amines do display some subtle kinetic differences: for instance, compound 2 has a shallower dependence of CPET rate constants on driving force (Br?nsted α, Δ ln(k)/Δ ln(K(eq))). These results show that the simple tunneling model is not a good predictor of the effect of proton donor-acceptor distance on concerted-electron transfer reactions involving strongly hydrogen-bonded systems. Computational analysis of the observed similarity of the two phenols emphasizes the importance of the highly anharmonic O···H···N potential energy surface and the influence of proton vibrational excited states.     

Deuterium kinetic isotope effects on the gas-phase reactions of C2H with H2(D2) and CH4(CD4)

Kinetics of the ethynyl (C(2)H) radical reactions with H(2), D(2), CH(4) and CD(4) was studied over the temperature range of 295-396 K by a pulsed laser photolysis/chemiluminescence technique. The C(2)H radicals were generated by ArF excimer-laser photolysis of C(2)H(2) or CF(3)C(2)H and were monitored by the chemiluminescence of CH(A(2)Δ) produced by their reaction with O(2) or O((3)P). The measured absolute rate constants for H(2) and CH(4) agreed well with the available literature data. The primary kinetic isotope effects (KIEs) were determined to be k(H(2))/k(D(2)) = 2.48 ± 0.14 and k(CH(4))/k(CD(4)) = 2.45 ± 0.16 at room temperature. Both of the KIEs increased as the temperature was lowered. The KIEs were analyzed by using the variational transition state theory with semiclassical small-curvature tunneling corrections. With anharmonic corrections on the loose transitional vibrational modes of the transition states, the theoretical predictions satisfactorily reproduced the experimental KIEs for both C(2)H + H(2)(D(2)) and C(2)H + CH(4)(CD(4)) reactions.     

Kinetic isotope effects for concerted multiple proton transfer: a direct dynamics study of an active-site model of carbonic anhydrase II

The rate constant of the reaction catalyzed by the enzyme carbonic anhydrase II, which removes carbon dioxide from body fluids, is calculated for a model of the active site. The rate-determining step is proton transfer from a zinc-bound water molecule to a histidine residue via a bridge of two or more water molecules. The structure of the active site is known from X-ray studies except for the number and location of the water molecules. Model calculations are reported for a system of 58 atoms including a four-coordinated zinc ion connected to a methylimidazole molecule by a chain of two waters, constrained to reproduce the size of the active site. The structure and vibrational force field are calculated by an approximate density functional treatment of the proton-transfer step at the Self-Consistent-Charge Density Functional Tight Binding (SCC-DFTB) level. A single transition state is found indicating concerted triple proton transfer. Direct-dynamics calculations for proton and deuteron transfer and combinations thereof, based on the Approximate Instanton Method and on Variational Transition State Theory with Tunneling Corrections, are in fair agreement and yield rates that are considerably higher and kinetic isotope effects (KIEs) that are somewhat higher than experiment. Classical rate constants obtained from Transition State Theory are smaller than the quantum values but the corresponding KIEs are five times larger. For multiple proton transfer along water bridges classical KIEs are shown to be generally larger than quantum KIEs, which invalidates the standard method to distinguish tunneling and over-barrier transfer. In the present case, a three-way comparison of classical and quantum results with the observed data is necessary to conclude that proton transfer along the bridge proceeds by tunneling. The results suggest that the two-water bridge is present in low concentrations but makes a substantial contribution to proton transport because of its high efficiency. Bridging structures containing more water molecules may have lower energies but are expected to be less efficient. The observed exponential dependence of the KIEs on the deuterium concentration in H(2)O/D(2)O mixtures implies concerted transfer and thus rules out substantial contributions from structures that lead to stepwise transfer via solvated hydronium ions, which presumably dominate proton transfer in less efficient carbonic anhydrase isozymes.     

Alpha-secondary isotope effects as probes of "tunneling-ready" configurations in enzymatic H-tunneling: insight from environmentally coupled tunneling models

Using alpha-secondary kinetic isotope effects (2 degrees KIEs) in conjunction with primary (1 degrees ) KIEs, we have investigated the mechanism of environmentally coupled hydrogen tunneling in the reductive half-reactions of two homologous flavoenzymes, morphinone reductase (MR) and pentaerythritol tetranitrate reductase (PETNR). We find exalted 2 degrees KIEs (1.17-1.18) for both enzymes, consistent with hydrogen tunneling. These 2 degrees KIEs, unlike 1 degrees KIEs, are independent of promoting motions-a nonequilibrium pre-organization of cofactor and active site residues that is required to bring the reactants into a "tunneling-ready" configuration. That these 2 degrees KIEs are identical suggests the geometries of the "tunneling-ready" configurations in both enzymes are indistinguishable, despite the fact that MR, but not PETNR, has a clearly temperature-dependent 1 degrees KIE. The work emphasizes the benefit of combining studies of 1 degrees and 2 degrees KIEs to report on pre-organization and local geometries within the context of contemporary environmentally coupled frameworks for H-tunneling.     

Hydrogen atom vs electron transfer in catecholase-mimetic oxidations by superoxometal complexes. Deuterium kinetic isotope effects

Dioximato-cobalt(II), -iron(II) and -manganese(II) complexes (1)-(6), acting as functional catecholase and phenoxazinone synthase models, exhibit a deuterium kinetic isotope effect predicted by theory (k4H/k4D < or = 3) in the catalytic oxidative dehydrogenation of 3,5-di-tert-butylcatechol and 2-aminophenol by O2. KIEs in the range of (k4H/k4D approximately 1.79-3.51) are observed with (1) and (2) as catalysts, pointing to hydrogen atom transfer in the rate-determining step from the substrate hydroxy group to the metal-bound superoxo ligand. Less significant KIEs (1.06-1.20) are exhibited by catalysts systems (3)-(6), indicating that proton-coupled electron transfer is the preferred route in those cases.     

Theoretical Study on the Gas‐Phase SN2 Reaction of Microhydrated Fluoride with Methyl Fluoride

The rate constants for the gas‐phase S_N2 reaction of F^?(H₂O) with CH₃F have been calculated using the dual‐level variational transition state theory including multidimensional tunneling from 50 to 500 K. Tunneling was found to dominate the reaction below 200 K. The deuterium, ¹³C, and ¹⁴C kinetic isotope effects (KIEs) and solvent (D₂O) isotope effects (SKIEs) were also calculated in the same temperature range. The results indicated that the deuterium and heavy water substitutions resulted in inverse KIEs (0.6～0.8 ) while the ¹³C and ¹⁴C substitutions resulted in normal KIEs (1.0～1.2) at room temperature. The calculated carbon KIEs increased significantly below 80 K due to the differences in the magnitude of the tunneling effects for different isotopic substitutions.     

Insights into the hydrogen-abstraction reactions of diol dehydratase: relevance to the catalytic mechanism and suicide inactivation

High-level quantum chemistry calculations have been used to examine the hydrogen-abstraction reactions of diol dehydratase (DDH) in the context of both the catalytic mechanism and the enzyme dysfunction phenomenon termed suicide inactivation. The barriers for the catalytic hydrogen-abstraction reactions of ethane-1,2-diol and propane-1,2-diol are examined in isolation, as well as in the presence of various Br?nsted acids and bases. Modest changes in the magnitudes of the initial and final abstraction barriers are seen, depending on the strength of the acid or base, and on whether these effects are considered individually or together. The most significant changes (ca. 20 kJ mol(-1)) are found for the initial abstraction barrier when the spectator OH group is partially deprotonated. Kinetic isotope effects including Eckart tunneling corrections (KIEs) have also been calculated for these model systems. We find that contributions from tunneling are of a magnitude similar to that of the contributions from semiclassical theory alone, meaning that quantum effects serve to significantly accelerate the rate of hydrogen transfer. The calculated KIEs for the partially deprotonated system are in qualitative agreement with experimentally determined values. In complementary investigations, the ability of DDH to become deactivated by certain substrate analogues is examined. In all cases, the formation of a stable radical intermediate causes the hydrogen re-abstraction step to become an extremely endothermic process. The consequent inability of 5'-deoxyadenosyl radical to be regenerated breaks the catalytic cycle, resulting in the suicide inactivation of DDH.     

18O kinetic isotope effects in non-heme iron enzymes: probing the nature of Fe/O2 intermediates

Contrasted here are the competitive 18O/16O kinetic isotope effects (18O KIEs) on kcat/Km(O2) for three non-heme iron enzymes that activate O2 at an iron center coordinated by a 2-His-1-carboxylate facial triad: taurine dioxygenase (TauD), (S)-(2)-hydroxypropylphosphonic acid epoxidase (HppE), and 1-aminocyclopropyl-1-carboxylic acid oxidase (ACCO). Measured 18O KIEs of 1.0102 +/- 0.0002 (TauD), 1.0120 +/- 0.0002 (HppE), and 1.0215 +/- 0.0005 (ACCO) suggest the formation in the rate-limiting step of O2 activation of an FeIII-peroxohemiketal, FeIII-OOH, and FeIV O species, respectively. The comparison of the measured 18O KIEs with calculated or experimental 18O equilibrium isotope effects (18O EIEs) provides new insights into the O2 activation through an inner-sphere mechanism at a non-heme iron center.     

Synthesis of isoreticular zinc(II)-phosphonocarboxylate frameworks and their application in the Friedel-Crafts benzylation reaction

Three isoreticular zinc(II)-phosphonocarboxylate frameworks, namely {[Zn(3)(pbdc)(2)]·2H(3)O}(n) (ZnPC-2), {[Zn(3)(pbdc)(2)]·Hpd·H(3)O·4H(2)O}(n) (Hpd@ZnPC-2) and {[Co(1.5)Zn(1.5)(pbdc)(2)]·2H(3)O}(n) (CoZnPC-2) (H(4)pbdc=5-phosphonobenzene-1,3-dicarboxylic acid, pd=pyrrolidine), were solvothermally synthesized. ZnPC-2 has a 3D structure based on trinuclear Zn(II) clusters (Zn(3)-SBU) showing 3D interconnected channels. Hpd@ZnPC-2 contains an isoreticular framework of ZnPC-2 with small channels blocked by Hpd molecules. In CoZnPC-2, Zn(II) ions in ZnPC-2 are partially substituted by Co(II) ions. The Friedel-Crafts benzylation reactions were carried out over these isoreticular porous materials. The catalytic results reveal that ZnPC-2 is an excellent heterogeneous Lewis acid catalyst with a high selectivity (>90%) towards less bulky para-oriented products. The catalytic reaction has been proved to occur inside the pore of ZnPC-2, and the immobilized Zn(3)-SBUs are the active sites.     

Tunneling in C-H oxidation reactions by an oxoiron(IV) porphyrin radical cation: direct measurements of very large H/D kinetic isotope effects

Rate constants for oxidations of benzyl alcohol-d0 and -d7 by oxoiron(IV) tetramesitylporphyrin radical cation perchlorate in acetonitrile were measured in single turnover kinetic studies. The kinetic isotope effect (kH/kD) increased from 28 at 23 degrees C to 360 at -30 degrees C due to extensive hydrogen atom tunneling that was analyzed in terms of a parabolic energy barrier to tunneling. Similarly, large KIE values were found for oxidations of ethylbenzene-d0 and -d10 at room temperature. The large KIE values are a function of the porphyrin identity, and porphyrins containing electron-withdrawing groups display normal KIEs. KIEs found under catalytic turnover conditions are somewhat smaller than those obtained in single turnover reactions. The results should serve as benchmarks for computational studies of C-H oxidations by porphyrin and heme-iron-oxo systems.     

星际媒介H2NCH2CN与H2O反应中的H2O分子偶合质子转移机理和氢隧道效应

H2NCH2CN+H2O→H2NCH2C(OH)NH是一个重要的反应, 涉及到星际媒介中甘氨酸的形成, 与早期地球上的氨基酸起源有关. 如果没有考虑氢隧道效应, 在MP2/6-311+G(d,p)级别上计算反应能垒是254.7 kJ·mol-1, 在星际媒介中该气相反应很难进行. 在星际媒介冰颗粒表面上, 水分子催化反应增强了该化学反应的活性. H2NCH2CN与(H2O)3反应中的两个水分子作为催化剂降低活化能77.5 kJ·mol-1和活化自由能70.9 kJ·mol-1, 并且通过氢键桥协同传递质子. 量子氢隧道对于该反应进行至关紧要,采用小弯曲隧道(SCT)近似和正则变分过渡态理论(CVT)方法研究. 温度50 K时, 速率常数kSCT/CVT为1.86×10-23 cm3·molecule-1·s-1, 表明在星际媒介中通过质子隧道机理该反应容易进行. 研究结果与地球上的氨基酸起源于地球本身物质的观点相一致.     

Investigating Inner-Sphere Reorganization via Secondary Kinetic Isotope Effects in the C-H Cleavage Reaction Catalyzed by Soybean Lipoxygenase: Tunneling in the Substrate Backbone as Well as the Transferred Hydrogen

This work describes the application of NMR to the measurement of secondary deuterium (2° (2)H) and carbon-13 ((13)C) kinetic isotope effects (KIEs) at positions 9-13 within the substrate linoleic acid (LA) of soybean lipoxygenase-1. The KIEs have been measured using LA labeled with either protium (11,11-h2-LA) or deuterium (11,11-d2-LA) at the reactive C11 position, which has been previously shown to yield a primary deuterium isotope effect of ca. 80. The conditions of measurement yield the intrinsic 2° (2)H and (13)C KIEs on k(cat)/K(m) directly for 11,11-d2-LA, whereas the values for the 2° (2)H KIEs for 11,11-h2-LA are obtained after correction for a kinetic commitment. The pattern of the resulting 2° (2)H and (13)C isotope effects reveals values that lie far above those predicted from changes in local force constants. Additionally, many of the experimental values cannot be modeled by electronic effects, torsional strain, or the simple inclusion of a tunneling correction to the rate. Although previous studies have shown the importance of extensive tunneling for cleavage of the primary hydrogen at C11 of LA, the present findings can only be interpreted by extending the conclusion of nonclassical behavior to the secondary hydrogens and carbons that flank the position undergoing C-H bond cleavage. A quantum mechanical method introduced by Buhks et al. [J. Phys. Chem. 1981, 85, 3763] to model the inner-sphere reorganization that accompanies electron transfer has been shown to be able to reproduce the scale of the 2° (2)H KIEs.     

Isotope effects and heavy-atom tunneling in the Roush allylboration of aldehydes

Intermolecular (13)C kinetic isotope effects (KIEs) for the Roush allylboration of p-anisaldehyde were determined using a novel approach. The experimental (13)C KIEs fit qualitatively with the expected rate-limiting cyclic transition state, but they are far higher than theoretical predictions based on conventional transition state theory. This discrepancy is attributed to a substantial contribution of heavy-atom tunneling to the reaction, and this is supported by multidimensional tunneling calculations that reproduce the observed KIEs.     

Kinetic isotope effects calculated with the instanton method

The ring-opening reaction of the cyclopropylcarbinyl radical proceeds via heavy-atom tunneling at low temperature. We used instanton theory to calculate tunneling rates and kinetic isotope effects with on-the-fly calculation of energies by density functional theory (B3LYP). The accuracy was verified by explicitly correlated coupled-cluster calculations (UCCSD(T)-F12). At cryogenic temperatures, we found protium/deuterium KIEs up to 13 and inverse KIEs down to 0.2. We also studied an intramolecular tautomerization reaction. A simple and computationally efficient method is proposed to calculate KIEs with the instanton method: the instanton path is assumed to be independent of the atomic masses. This results in surprisingly good estimates of the KIEs for the cyclopropylcarbinyl radical and for the secondary KIEs of the tautomerization. Challenges and capabilities of the instanton method for calculating KIEs are discussed.     

Y2O3/Ag3VO4复合可见光催化剂的制备及其光催化性能

以Na3VO4.12H2O,AgNO3和Y(NO3)3.6H2O为原料,采用浸渍法制备了Y2O3/Ag3VO4复合可见光催化剂,并用XRD,SEM,XPS,UV-Vis等测试手段表征了试样的结构和性能。结果显示,Y2 O3/Ag3VO4复合可见光催化剂为单斜结构,Y以Y2 O3的形式分散在Ag3VO4晶体的表面。UV-Vis测试结果表明,Y2O3/Ag3VO4较纯Ag3VO4吸收带边发生了红移,在可见光区的吸收增强;以金属卤灯(波长大于400 nm)为光源,研究了Y2O3/Ag3VO4催化剂对甲基橙(MO)的可见光催化降解性能。结果发现,Y2O3/Ag3VO4复合可见光催化剂的光催化活性较纯Ag3VO4均有大幅提高,其中Y掺杂量为4%时活性最高。     

Kinetic isotope effects on the rate-limiting step of heme oxygenase catalysis indicate concerted proton transfer/heme hydroxylation

Heme oxygenase (HO) catalyzes the O2 and NADPH/cytochrome P450 reductase-dependent conversion of heme to biliverdin, free iron ion, and CO through a process in which the heme participates as both dioxygen-activating prosthetic group and substrate. We earlier confirmed that the first step of HO catalysis is a monooxygenation in which the addition of one electron and two protons to the HO oxy-ferroheme produces ferric-alpha-meso-hydroxyheme (h). Cryoreduction/EPR and ENDOR measurements further showed that hydroperoxo-ferri-HO converts directly to h in a single kinetic step without formation of a Compound I. We here report details of that rate-limiting step. One-electron 77 K cryoreduction of human oxy-HO and annealing at 200 K generates a structurally relaxed hydroperoxo-ferri-HO species, denoted R. We here report the cryoreduction/annealing experiments that directly measure solvent and secondary kinetic isotope effects (KIEs) of the rate-limiting R --> h conversion, using enzyme prepared with meso-deuterated heme and in H2O/D2O buffers to measure the solvent KIE (solv-KIE), and the secondary KIE (sec-KIE) associated with the conversion. This approach is unique in that KIEs measured by monitoring the rate-limiting step are not susceptible to masking by KIEs of other processes, and these results represent the first direct measurement of the KIEs of product formation by a kinetically competent reaction intermediate in any dioxygen-activating heme enzyme.The observation of both solv-KIE(298) = 1.8 and sec-KIE(298) = 0.8 (inverse) indicates that the rate-limiting step for formation of h by HO is a concerted process: proton transfer to the hydroperoxo-ferri-heme through the distal-pocket H-bond network, likely from a carboxyl group acting as a general acid catalyst, occurring in synchrony with bond formation between the terminal hydroperoxo-oxygen atom and the alpha-meso carbon to form a tetrahedral hydroxylated-heme intermediate. Subsequent rearrangement and loss of H2O then generates h.     

Large Isotope Effects in Organometallic Chemistry

The kinetic isotope effect (KIE) is key to understanding reaction mechanisms in many areas of chemistry and chemical biology, including organometallic chemistry. This ratio of rate constants, k_H/k_D, typically falls between 1–7. However, KIEs up to 105 have been reported, and can even be so large that reactivity with deuterium is unobserved. We collect here examples of large KIEs across organometallic chemistry, in catalytic and stoichiometric reactions, along with their mechanistic interpretations. Large KIEs occur in proton transfer reactions such as protonation of organometallic complexes and clusters, protonolysis of metal–carbon bonds, and dihydrogen reactivity. C−H activation reactions with large KIEs occur with late and early transition metals, photogenerated intermediates, and abstraction by metal-oxo complexes. We categorize the mechanistic interpretations of large KIEs into the following three types: (a) proton tunneling, (b) compound effects from multiple steps, and (c) semi-classical effects on a single step. This comprehensive collection of large KIEs in organometallics provides context for future mechanistic interpretation.     

Microwave measurements of proton tunneling and structural parameters for the propiolic acid-formic acid dimer

Microwave spectra of the propiolic acid-formic acid doubly hydrogen bonded complex were measured in the 1 GHz to 21 GHz range using four different Fourier transform spectrometers. Rotational spectra for seven isotopologues were obtained. For the parent isotopologue, a total of 138 a-dipole transitions and 28 b-dipole transitions were measured for which the a-dipole transitions exhibited splittings of a few MHz into pairs of lines and the b-type dipole transitions were split by ~580 MHz. The transitions assigned to this complex were fit to obtain rotational and distortion constants for both tunneling levels: A(0+) = 6005.289(8), B(0+) = 930.553(8), C(0+) = 803.9948(6) MHz, Δ(0+)(J) = 0.075(1), Δ(0+)(JK) = 0.71(1), and δ(0+)(j) = -0.010(1) kHz and A(0-) = 6005.275(8), B(0-) = 930.546(8), C(0-) = 803.9907(5) MHz, Δ(0-)(J) = 0.076(1), Δ(0-)(JK) = 0.70(2), and δ(0-)(j) = -0.008(1) kHz. Double resonance experiments were used on some transitions to verify assignments and to obtain splittings for cases when the b-dipole transitions were difficult to measure. The experimental difference in energy between the two tunneling states is 291.428(5) MHz for proton-proton exchange and 3.35(2) MHz for the deuterium-deuterium exchange. The vibration-rotation coupling constant between the two levels, F(ab), is 120.7(2) MHz for the proton-proton exchange. With one deuterium atom substituted in either of the hydrogen-bonding protons, the tunneling splittings were not observed for a-dipole transitions, supporting the assignment of the splitting to the concerted proton tunneling motion. The spectra were obtained using three Flygare-Balle type spectrometers and one chirped-pulse machine at the University of Virginia. Rotational constants and centrifugal distortion constants were obtained for HCOOH···HOOCCCH, H(13)COOH···HOOCCCH, HCOOD···HOOCCCH, HCOOH···DOOCCCH, HCOOD···DOOCCCH, DCOOH···HOOCCCH, and DCOOD···HOOCCCH. High-level ab initio calculations provided initial rotational constants for the complex, structural parameters, and some details of the proton tunneling potential energy surface. A least squares fit to the isotopic data reveals a planar structure that is slightly asymmetric in the OH distances. The formic OH···O propiolic hydrogen bond length is 1.8 ? and the propiolic OH···O formic hydrogen bond length is 1.6 ?, for the equilibrium configuration. The magnitude of the dipole moment was experimentally determined to be 1.95(3) × 10(-30) C m (0.584(8) D) for the 0(+) states and 1.92(5) × 10(-30) C m (0.576(14) D) for the 0(-) states.     

Large tunneling effect on the hydrogen transfer in bis(μ-oxo)dicopper enzyme: a theoretical study

Type-III copper-containing enzymes have dicopper centers in their active sites and exhibit a novel capacity for activating aliphatic C-H bonds in various substrates by taking molecular oxygen. Dicopper enzyme models developed by Tolman and co-workers reveal exceptionally large kinetic isotope effects (KIEs) for the hydrogen transfer process, indicating a significant tunneling effect. In this work, we demonstrate that variational transition state theory allows accurate prediction of the KIEs and Arrhenius parameters for such model systems. This includes multidimensional tunneling based on state-of-the-art quantum-mechanical calculations of the minimum-energy path (MEP). The computational model of bis(μ-oxo)dicopper enzyme consists of 70 atoms, resulting in a 204-dimensional potential energy surface. The calculated values of E(a)(H) - E(a)(D), A(H)/A(D), and the KIE at 233 K are -1.86 kcal/mol, 0.51, and 28.1, respectively, for the isopropyl ligand system. These values agree very well with experimental values within the limits of experimental error. For the representative tunneling path (RTP) at 233 K, the pre- and post-tunneling configurations are 3.3 kcal/mol below the adiabatic energy maximum, where the hydrogen travels 0.54 ? by tunneling. We found that tunneling is very efficient for hydrogen transfer and that the RTP is very different from the MEP. It is mainly heavy atoms that move as the reaction proceeds from the reactant complex to the pretunneling configuration, and the hydrogen atom suddenly hops at that point.