Mechanism of the Tyrosine Ammonia Lyase Reaction—Tandem Nucleophilic and Electrophilic Enhancement by a Proton Transfer

Quantum mechanics/molecular mechanics calculations in tyrosine ammonia lyase (TAL) ruled out the hypothetical Friedel–Crafts (FC) route for ammonia elimination from L ‐tyrosine due to the high energy of FC intermediates. The calculated pathway from the zwitterionic L ‐tyrosine‐binding state (0.0 kcal mol^?1) to the product‐binding state ((E)‐coumarate+H₂N? MIO; ?24.0 kcal mol^?1; MIO=3,5‐dihydro‐5‐methylidene‐4H‐imidazol‐4‐one) involves an intermediate (IS, ?19.9 kcal mol^?1), which has a covalent bond between the N atom of the substrate and MIO, as well as two transition states (TS1 and TS2). TS1 (14.4 kcal mol^?1) corresponds to a proton transfer from the substrate to the N1 atom of MIO by Tyr300? OH. Thus, a tandem nucleophilic activation of the substrate and electrophilic activation of MIO happens. TS2 (5.2 kcal mol^?1) indicates a concerted C? N bond breaking of the N‐MIO intermediate and deprotonation of the pro‐S β position by Tyr60. Calculations elucidate the role of enzymic bases (Tyr60 and Tyr300) and other catalytically relevant residues (Asn203, Arg303, and Asn333, Asn435), which are fully conserved in the amino acid sequences and in 3D structures of all known MIO‐containing ammonia lyases and 2,3‐aminomutases.     

人工设计逆醛缩酶RA95.5-8F催化β-羟基酮C—C裂解的理论研究

采用量子力学与分子力学组合(QM/MM)方法对人工设计逆醛缩酶RA95.5-8F催化β-羟基酮化合物裂解反应的机理进行了研究. 结果表明, 裂解反应主要包括赖氨酸Lys1083对底物的亲核进攻、 Schiff碱形成、 烯胺水解及C—N断裂等过程, C—N键裂解生成丙酮为整个反应的决速步骤, 能垒为106.27 kJ/mol; 活性中心的赖氨酸Lys1083、 酪氨酸Tyr1051、 天冬酰胺Asn1110和酪氨酸Tyr1180构成一个催化四联体, Lys1083通过与底物形成席夫碱对底物进行活化, Tyr1051作为催化酸碱参与质子转移过程, 催化四联体的氢键网络有利于反应过渡态的稳定并使R-构型的底物更容易结合在活性位点, 导致RA95.5-8F对R构型底物具有高的选择性和催化活性.     

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

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

Experimental and theoretical studies on the thermal decomposition of heterocyclic nitrosimines

A series of substituted 2-nitrosiminobenzothiazolines (2) were synthesized by the nitrosation of the corresponding 2-iminobenzothiazolines (6). Thermal decomposition of 2a--f and of the seleno analogue 7 in methanol and of 3-methyl-2-nitrosobenzothiazoline (2a) in acetonitrile, 1,4-dioxane, and cyclohexane followed first-order kinetics. The activation parameters for thermal deazetization of 2a were measured in cyclohexane (Delta H(++) = 25.3 +/- 0.5 kcal/mol, Delta S(++) = 1.3 +/- 1.5 eu) and in methanol (Delta H(++) = 22.5 +/- 0.7 kcal/mol, Delta S(++) = -12.9 +/- 2.1 eu). These results indicate a unimolecular decomposition and are consistent with a proposed stepwise mechanism involving cyclization of the nitrosimine followed by loss of N(2). The ground-state conformations of the parent nitrosiminothiazoline (9a) and transition states for rotation around the exocyclic C==N bond, electrocyclic ring closure, and loss of N(2) were calculated using ab initio molecular orbital theory at the MP2/6-31G* level. The calculated gas-phase barrier height for the loss of N(2) from 9a (25.2 kcal/mol, MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE) compares favorably with the experimental barrier for 2a of 25.3 kcal/mol in cyclohexane. The potential energy surface is unusual; the rotational transition state 9a-rot-ts connects directly to the orthogonal transition state for ring-closure 9aTS. The decoupling of rotational and pseudopericyclic bond-forming transition states is contrasted with the single pericyclic transition state (15TS) for the electrocyclic ring-opening of oxetene (15) to acrolein (16). For comparison, the calculated homolytic strength of the N--NO bond is 40.0 kcal/mol (MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE).     

Reactivity of alkyl versus silyl peroxides. The consequences of 1, 2-silicon bridging on the epoxidation of alkenes with silyl hydroperoxides and bis(trialkylsilyl)peroxides

The bond dissociation energies for a series of silyl peroxides have been calculated at the G2 and CBS-Q levels of theory. A comparison is made with the O-O BDE of the corresponding dialkyl peroxides, and the effect of the O-O bond strength on the activation barrier for oxygen atom transfer is discussed. The O-O bond dissociation enthalpies (DeltaH(298)) for bis (trimethylsilyl) peroxide (1) and trimethylsilyl hydroperoxide (2) are 54.8 and 53.1 kcal/mol, respectively at the G2 (MP2) and CBS-Q levels of theory. The O-O bond dissociation energies computed at G2 and G2(MP2) levels for bis(tert-butyl) peroxide and tert-butyl hydroperoxide are 45.2 and 48.3 kcal/mol, respectively. The barrier height for 1,2-methyl migration from silicon to oxygen in trimethylsilyl hydroperoxide is 47.9 kcal/mol (MP4//MP2/6-31G). The activation energy for the oxidation of trimethylamine to its N-oxide by bis(trimethylsilyl) peroxide is 28.2 kcal/mol (B3LYP/6-311+G(3df,2p)// B3LYP/6-31G(d)). 1,2-Silicon bridging in the transition state for oxygen atom transfer to a nucleophilic amine results in a significant reduction in the barrier height. The barrier for the epoxidation of E-2-butene with bis(dimethyl(trifluoromethyl))silyl peroxide is 25.8 kcal/mol; a reduction of 7.5 kcal/mol relative to epoxidation with 1. The activation energy calculated for the epoxidation of E-2-butene with F(3)SiOOSiF(3) is reduced to only 2.2 kcal/mol reflecting the inductive effect of the electronegative fluorine atoms.     

Determination of the C4-H bond dissociation energies of NADH models and their radical cations in acetonitrile

Heterolytic and homolytic bond dissociation energies of the C4-H bonds in ten NADH models (seven 1,4-dihydronicotinamide derivatives, two Hantzsch 1,4-dihydropyridine derivatives, and 9,10-dihydroacridine) and their radical cations in acetonitrile were evaluated by titration calorimetry and electrochemistry, according to the four thermodynamic cycles constructed from the reactions of the NADH models with N,N,N',N'-tetramethyl-p-phenylenediamine radical cation perchlorate in acetonitrile (note: C9-H bond rather than C4-H bond for 9,10-dihydroacridine; however, unless specified, the C9-H bond will be described as a C4-H bond for convenience). The results show that the energetic scales of the heterolytic and homolytic bond dissociation energies of the C4-H bonds cover ranges of 64.2-81.1 and 67.9-73.7 kcal mol(-1) for the neutral NADH models, respectively, and the energetic scales of the heterolytic and homolytic bond dissociation energies of the (C4-H)(.+) bonds cover ranges of 4.1-9.7 and 31.4-43.5 kcal mol(-1) for the radical cations of the NADH models, respectively. Detailed comparison of the two sets of C4-H bond dissociation energies in 1-benzyl-1,4-dihydronicotinamide (BNAH), Hantzsch 1,4-dihydropyridine (HEH), and 9,10-dihydroacridine (AcrH(2)) (as the three most typical NADH models) shows that for BNAH and AcrH(2), the heterolytic C4-H bond dissociation energies are smaller (by 3.62 kcal mol(-1)) and larger (by 7.4 kcal mol(-1)), respectively, than the corresponding homolytic C4-H bond dissociation energy. However, for HEH, the heterolytic C4-H bond dissociation energy (69.3 kcal mol(-1)) is very close to the corresponding homolytic C4-H bond dissociation energy (69.4 kcal mol(-1)). These results suggests that the hydride is released more easily than the corresponding hydrogen atom from BNAH and vice versa for AcrH(2), and that there are two almost equal possibilities for the hydride and the hydrogen atom transfers from HEH. Examination of the two sets of the (C4-H)(.+) bond dissociation energies shows that the homolytic (C4-H)(.+) bond dissociation energies are much larger than the corresponding heterolytic (C4-H)(.+) bond dissociation energies for the ten NADH models by 23.3-34.4 kcal mol(-1); this suggests that if the hydride transfer from the NADH models is initiated by a one-electron transfer, the proton transfer should be more likely to take place than the corresponding hydrogen atom transfer in the second step. In addition, some elusive structural information about the reaction intermediates of the NADH models was obtained by using Hammett-type linear free-energy analysis.     

Cheletropic decomposition of cyclic nitrosoamines revisited: the nature of the transition states and a critical role of the ring strain

The cheletropic decompositions of 1-nitrosoaziridine (1), 1-nitroso-Delta(3)-pyrroline (2), 7-nitroso-7-azabicyclo[2.2. 1]hepta-2,5-diene (3), and 6-nitroso-6-azabicyclo[2.1.1]hexa-4-ene (4) have been studied theoretically using high level ab initio computations. Activation parameters of the decomposition of nitrosoaziridine 1 were obtained experimentally in heptane (DeltaH()(298) = 18.6 kcal mol(-)(1), DeltaS()(298) = -7.6 cal mol(-)(1) K(-)(1)) and methanol (20.3 kcal mol(-)(1), 0.3 cal mol(-)(1) K(-)(1)). Among employed theoretical methods (B3LYP, MP2, CCD, CCSD(T)//CCD), the B3LYP method in conjunction with 6-31+G, 6-311+G, and 6-311++G(3df,2pd) basis sets gives the best agreement with experimental data. It was found that typical N-nitrosoheterocycles 2-4 which have high N-N bond rotation barriers (>16 kcal mol(-)(1)) extrude nitrous oxide via a highly asynchronous transition state with a planar ring nitrogen atom. Nitrosoaziridine 1, with a low rotation barrier (<9 kcal mol(-)(1)) represents a special case. This compound can eliminate N(2)O via a low energy linear synperiplanar transition state (DeltaH()(298) = 20.6 kcal mol(-)(1), DeltaS()(298) = 2.5 cal mol(-)(1) K(-)(1)). Two higher energy transition states are also available. The B3LYP activation barriers of the cheletropic fragmentation of nitrosoheterocycles 2-4 decrease in the series: 2 (58 kcal mol(-)(1)) > 3 (18 kcal mol(-)(1)) > 4 (12) kcal mol(-)(1). The relative strain energies increase in the same order: 2 (0 kcal mol(-)(1)) < 3 (39 kcal mol(-)(1)) < 4 (52 kcal mol(-)(1)). Comparison of the relative energies of 2-4 and their transition states on a common scale where the energy of nitrosopyrroline 2 is assumed as reference indicates that the thermal stability of the cyclic nitrosoamines toward cheletropic decomposition is almost entirely determined by the ring strain.     

Computational investigation of hydrogen abstraction from 2-aminoethanol by the 1,5-dideoxyribose-5-yl radical: a model study of a reaction occurring in the active site of ethanolamine ammonia lyase

Hydrogen abstraction from 2-aminoethanol by the 5'-deoxyadenosyl radical, which is formed upon Co--C bond homolysis in coenzyme B(12), was investigated by theoretical means with employment of the DFT (B3LYP) and ab initio (MP2) approaches. As a model system for the 5'-deoxyadenosyl moiety the computationally less demanding 1,5-dideoxyribose was employed; two conformers, which differ in ring conformation (C2- and C3-endo), were considered. If hydrogen is abstracted from "free" substrate by the C2-endo conformer of the 1,5-dideoxyribose-5-yl radical, the activation enthalpy is 16.7 kcal mol(-1); with the C3-endo counterpart, the value is 17.3 kcal mol(-1). These energetic requirements are slightly above the activation enthalpy limit (15 kcal mol(-1)) determined experimentally for the rate-determining step of the sequence, that is, hydrogen delivery from 5'-deoxyadenosine to the product radical. The activation enthalpy is lower when the substrate interacts with at least one amino acid from the active site. According to the computations, when a His model system partially protonates the substrate the activation enthalpy is 4.5 kcal mol(-1) for the C3-endo conformer and 5.8 kcal mol(-1) for the C2-endo counterpart. As hydrogen abstraction from the fully as well as the partially protonated substrate is preceded by the formation of quite stable encounter complexes, the actual activation barriers are around 13-15 kcal mol(-1). A synergistic interaction of 2-aminoethanol with two amino acids where His partially protonates the NH(2) group and Asp partially deprotonates the OH group of the substrate results in an activation enthalpy of 12.4 kcal mol(-1) for the C3-endo conformer and 13.2 kcal mol(-1) for the C2-endo counterpart. However, if encounter complexes exist in the active site, the actual activation barriers are much higher (>25 kcal mol(-1)) than that reported for the rate-determining step. These findings together with previous computations suggest that the energetics of the initial hydrogen abstraction decrease with an interaction of the substrate with only a protonating auxiliary, but for the rearrangement of the radical the synergistic effects of two auxiliaries are essential to pull the barrier below the limit of 15 kcal mol(-1).     

Transition-state characterization of the ammonia ionization process in aqueous solution via the free-energy gradient method

For the ionization process of ammonia in aqueous solution, the transition-state (TS) structure was fully optimized for the first time on the free-energy surface (FES) by applying the free-energy gradient (FEG) method combined with a hybrid quantum mechanical and molecular mechanical molecular dynamics (QM/MM-MD) method. In aqueous solution, the ionization process was found to proceed by way of a clear TS (R(N1-H5) = 1.512 A), which does not exist in the gas phase. The free-energy (FE) of activation for ionization obtained was 14.7 kcal/mol, within the classical approximation, via the QM/MM-MD FEG method, and is found to be in good agreement with 9.57 kcal/mol estimated from the TS theory using the experimental value of the rate constant. Apart from the dynamic correction, it is indicated that the theoretical value would be improved to be 10.28 kcal/mol if the electronic-state calculation could be executed at the B3LYP/6-31G(d) level of theory.     

Internal degrees of freedom,structural motifs,and conformational energetics of the 5'-deoxyadenosyl radical: implications for function in adenosylcobalamin-dependent enzymes. A computational study

The potential energy surface of the free 5'-deoxyadenosyl radical in the gas phase is explored using density functional and second-order M?ller-Plesset perturbation theories with 6-31G(d) and 6-31++G(d,p) basis sets and interpreted in terms of attractive and repulsive interactions. The 5',8-cyclization is found to be exothermic by approximately 20 kcal/mol but kinetically unfavorable; the lowest cyclization transition state (TS) lies about 7 kcal/mol higher than the highest TS for conversion between most of the open isomers. In open isomers, the two energetically most important attractive interactions are the hydrogen bonds (a) between the 2'-OH group and the N3 adenine center and (b) between the 2'-OH and 3'-OH groups. The relative ribose-adenine rotation about the C1'-N9 glycosyl bond in a certain range changes the energy by as much as 10-15 kcal/mol, the origin being (i) the repulsive 2'-H.H-C8 and O1'.N3 and (ii) the attractive 2'-OH.N3 ribose-adenine interactions. The hypothetical synergy between the glycosyl rotation and the Co-C bond scission may contribute to the experimentally established labilization of the Co-C bond in enzyme-bound adenosylcobalamin. The computational results are not inconsistent with the rotation about the C1'-N9 glycosyl bond being the principal coordinate for long-range radical migration in coenzyme B(12)-dependent enzymes. The effect of the protein environment on the model system results reported here remains an open question.     

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

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

Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions

Chorismate mutase is at the centre of current controversy about fundamental features of biological catalysts. Some recent studies have proposed that catalysis in this enzyme does not involve transition state (TS) stabilization but instead is due largely to the formation of a reactive conformation of the substrate. To understand the origins of catalysis, it is necessary to compare equivalent reactions in different environments. The pericyclic conversion of chorismate to prephenate catalysed by chorismate mutase also occurs (much more slowly) in aqueous solution. In this study we analyse the origins of catalysis by comparison of multiple quantum mechanics/molecular mechanics (QM/MM) reaction pathways at a reliable, well tested level of theory (B3LYP/6-31G(d)/CHARMM27) for the reaction (i) in Bacillus subtilis chorismate mutase (BsCM) and (ii) in aqueous solvent. The average calculated reaction (potential energy) barriers are 11.3 kcal mol(-1) in the enzyme and 17.4 kcal mol(-1) in water, both of which are in good agreement with experiment. Comparison of the two sets of reaction pathways shows that the reaction follows a slightly different reaction pathway in the enzyme than in it does in solution, because of a destabilization, or strain, of the substrate in the enzyme. The substrate strain energy within the enzyme remains constant throughout the reaction. There is no unique reactive conformation of the substrate common to both environments, and the transition state structures are also different in the enzyme and in water. Analysis of the barrier heights in each environment shows a clear correlation between TS stabilization and the barrier height. The average differential TS stabilization is 7.3 kcal mol(-1) in the enzyme. This is significantly higher than the small amount of TS stabilization in water (on average only 1.0 kcal mol(-1) relative to the substrate). The TS is stabilized mainly by electrostatic interactions with active site residues in the enzyme, with Arg90, Arg7 and Glu78 generally the most important. Conformational effects (e.g. strain of the substrate in the enzyme) do not contribute significantly to the lower barrier observed in the enzyme. The results show that catalysis is mainly due to better TS stabilization by the enzyme.     

Molecule-induced alkane homolysis with dioxiranes.

The mechanisms of C-H and C-C bond activations with dimethyldioxirane (DMD) were studied experimentally and computationally at the B3LYP/6-311+G**//B3LYP/6-31G* density functional theory level for the propellanes 3,6-dehydrohomoadamantane (2) and 1,3-dehydroadamantane (3). The sigma(C-C) activation of 3 with DMD (Delta G(*) = 23.9 kcal mol(-1) and Delta G(r) = -5.4 kcal mol(-1)) is the first example of a molecule-induced homolytic C-C bond cleavage. The C-H bond hydroxylation observed for 2 is highly exergonic (Delta G(r) = -74.4 kcal mol(-1)) and follows a concerted pathway (Delta G(*) = 34.8 kcal mol(-1)), in contrast to its endergonic molecule-induced homolysis (Delta G(*) = 28.8 kcal mol(-1) and Delta G(r) = +9.2 kcal mol(-1)). The reactivities of 2 and 3 with CrO(2)Cl(2), which follow a molecule-induced homolytic activation mechanism, parallel the DMD results only for highly reactive 3, but differ considerably for more stable propellanes such as 4-phenyl-3,6-dehydrohomoadamantane (1) and 2.     

Mechanism of N2O reduction by the mu4-S tetranuclear CuZ cluster of nitrous oxide reductase

Reaction thermodynamics and potential energy surfaces are calculated using density functional theory to investigate the mechanism of the reductive cleavage of the N-O bond by the mu(4)-sulfide-bridged tetranuclear Cu(Z) site of nitrous oxide reductase. The Cu(Z) cluster provides an exogenous ligand-binding site, and, in its fully reduced 4Cu(I) state, the cluster turns off binding of stronger donor ligands while enabling the formation of the Cu(Z)-N(2)O complex through enhanced Cu(Z) --> N(2)O back-donation. The two copper atoms (Cu(I) and Cu(IV)) at the ligand-binding site of the cluster play a crucial role in the enzymatic function, as these atoms are directly involved in bridged N(2)O binding, bending the ligand to a configuration that resembles the transition state (TS) and contributing the two electrons for N(2)O reduction. The other atoms of the Cu(Z) cluster are required for extensive back-bonding with minimal sigma ligand-to-metal donation for the N(2)O activation. The low reaction barrier (18 kcal mol(-)(1)) of the direct cleavage of the N-O bond in the Cu(Z)-N(2)O complex is due to the stabilization of the TS by a strong Cu(IV)(2+)-O(-) bond. Due to the charge transfer from the Cu(Z) cluster to the N(2)O ligand, noncovalent interactions with the protein environment stabilize the polar TS and reduce the activation energy to an extent dependent on the strength of proton donor. After the N-O bond cleavage, the catalytic cycle consists of a sequence of alternating protonation/one-electron reduction steps which return the Cu(Z) cluster to the fully reduced (4Cu(I)) state for future turnover.     

Theoretical studies of hetero-diels-alder reactions involving N-sulfinyl dienophiles

The gas-phase hetero-Diels-Alder reactions between butadiene and X-substituted sulfinyl dienophiles, O(-)-S(+)=N-X, are investigated theoretically at the B3LYP/6-31G level. The Z-forms of the dienophiles are found to be more stable (by 5-7 kcal mol(-)(1)) than the E-forms. Four modes of cycloadducts are considered: Z-endo; Z-exo; E(X)(-)(endo)(); E(X)(-)(exo)(). Five factors are responsible for the decreasing energetic preferences of the adducts in the order E(X)(-)(endo)() > E(X)(-)(exo)() > Z-endo > Z-exo: (i) The sigma-sigma proximate charge-transfer interactions in the TS; (ii) the relative sizes of the LUMO AO coefficients on S and N atoms; (iii) steric hindrance in the TS; (iv) the levels of the ground state and the LUMOs of the dienophile; (v) bond energies of the C-S and C-N bonds that are formed in the TS. All the reactions proceed concertedly, but the adduct formation is asynchronous. The endo-additions are favored over the exo-additions kinetically (lower DeltaG()) as well as thermodynamically (lower DeltaG degrees ). The major secondary orbital interaction determining the endo preference is that between the lone pair on N (n(N)) and the d(3) (C(3)-C(4)) sigma orbital (n(N)-sigma(d3)) interactions, whereas the larger AO lobe (LUMO) sizes on S favor a greater degree of d(5) (C-S) bond formation than d(6) (C-N) bond. The solvent, C(6)H(6), uniformly lowers the activation barriers so that the energetic preferences in the gas phase between various modes are maintained in solution.     

Bonding in ammonia borane: an analysis based on the natural orbitals for chemical valence and the extended transition state method (ETS-NOCV)

In the present study the natural orbitals for chemical valence (NOCVs) combined with the energy decomposition scheme (ETS) were used to characterize bonding in various clusters of ammonia borane (borazane): dimer D, trimer TR, tetramer TE, and the crystal based models: nonamer N and tetrakaidecamer TD. ETS-NOCV results have shown that shortening of the B-N bond (by ~0.1 ?) in ammonia borane crystal (as compared to isolated borazane molecule) is related to the enhancement of donation (by 6.5 kcal/mol) and electrostatic (by 11.3 kcal/mol) contributions. This, in turn, is caused solely by the electrostatic dipole-dipole interaction between ammonia borane units; dihydrogen bonding, BH···HN, formed between borazane units exhibits no direct impact on B-N bond contraction. On the other hand, formation of dihydrogen bonding appeared to be very important in the total stabilization of single borazane unit, namely, ETS-based data indicated that it leads to significant electronic stabilization ΔE(orb) = -17.5 kcal/mol, which is only slightly less important than the electrostatic term, ΔE(elstat) = -19.4 kcal/mol. Thus, both factors contribute to relatively high melting point of the borazane crystal. Deformation density contributions (Δρ(i)) obtained from NOCVs allowed to conclude that dihydrogen bonding is primarily based on outflow of electron density from B-H bonding orbitals to the empty σ*(N-H) (charge transfer component). Equally important is the covalent contribution resulting from the shift of the electron density from hydrogen atoms of both NH and BH groups to the interatomic regions of NH···HB. Quantitatively, averaged electronic strength of dihydrogen bond per one BH···HN link varies from 1.95 kcal/mol (for the crystal structure model, N), 2.47 kcal/mol (for trimer TR), through 2.65 kcal/mol (for tetramer TE), up to 3.95 kcal/mol (for dimer D).     

A comparative study of claisen and cope rearrangements catalyzed by chorismate mutase. An insight into enzymatic efficiency: transition state stabilization or substrate preorganization?

In this work we present a detailed analysis of the activation free energies and averaged interactions for the Claisen and Cope rearrangements of chorismate and carbachorismate catalyzed by Bacillus subtilischorismate mutase (BsCM) using quantum mechanics/molecular mechanics (QM/MM) simulation methods. In gas phase, both reactions are described as concerted processes, with the activation free energy for carbachorismate being about 10-15 kcal mol(-)(1) larger than for chorismate, at the AM1 and B3LYP/6-31G levels. Aqueous solution and BsCM active site environments reduce the free energy barriers for both reactions, due to the fact that in these media the two carboxylate groups can be approached more easily than in the gas phase. The enzyme specifically reduces the activation free energy of the Claisen rearrangement about 3 kcal mol(-)(1) more than that for the Cope reaction. This result is due to a larger transition state stabilization associated to the formation of a hydrogen bond between Arg90 and the ether oxygen. When this oxygen atom is changed by a methylene group, the interaction is lost and Arg90 moves inside the active site establishing stronger interactions with one of the carboxylate groups. This fact yields a more intense rearrangement of the substrate structure. Comparing two reactions in the same enzyme, we have been able to obtain conclusions about the relative magnitude of the substrate preorganization and transition state stabilization effects. Transition state stabilization seems to be the dominant effect in this case.     

What Is the Catalytic Mechanism of Enzymatic Histone N-Methyl Arginine Demethylation and Can It Be Influenced by an External Electric Field?

Arginine methylation is an important mechanism of epigenetic regulation. Some Fe(II) and 2-oxoglutarate dependent Jumonji-C (JmjC) Nϵ-methyl lysine histone demethylases also have N-methyl arginine demethylase activity. We report combined molecular dynamic (MD) and Quantum Mechanical/Molecular Mechanical (QM/MM) studies on the mechanism of N-methyl arginine demethylation by human KDM4E and compare the results with those reported for N-methyl lysine demethylation by KDM4A. At the KDM4E active site, Glu191, Asn291, and Ser197 form a conserved scaffold that restricts substrate dynamics; substrate binding is also mediated by an out of active site hydrogen-bond between the substrate Ser1 and Tyr178. The calculations imply that in either C−H or N−H potential bond cleaving pathways for hydrogen atom transfer (HAT) during N-methyl arginine demethylation, electron transfer occurs via a σ-channel; the transition state for the N−H pathway is ∼10 kcal/mol higher than for the C−H pathway due to the higher bond dissociation energy of the N−H bond. The results of applying external electric fields (EEFs) reveal EEFs with positive field strengths parallel to the Fe=O bond have a significant barrier-lowering effect on the C−H pathway, by contrast, such EEFs inhibit the N−H activation rate. The overall results imply that KDM4 catalyzed N-methyl arginine demethylation and N-methyl lysine demethylation occur via similar C−H abstraction and rebound mechanisms leading to methyl group hydroxylation, though there are differences in the interactions leading to productive binding of intermediates.     

Experimental determination of the nN --> sigma*P-O interaction energy of O-equatorial C-apical phosphoranes bearing a primary amino group

The reaction of a chlorophosphorane (9-Cl) with primary amines produced anti-apicophilic spirophosphoranes (5, O-equatorial phosphoranes), which violate the apicophilicity concept, having an apical carbon-equatorial oxygen configuration, along with the ordinarily expected O-apical stereoisomers (6) with the apical oxygen-equatorial carbon configuration. Although the amino group is electronegative in nature, the O-equatorial phosphoranes were found to be stable at room temperature and could still be converted to their more stable O-apical pseudorotamers (6) when they were heated in solution. X-ray analysis implied that this remarkable stability of the O-equatorial isomers could be attributed to the orbital interaction between the lone-pair electrons of the nitrogen atom (n(N)) and the antibonding sigma(P-O) orbital in the equatorial plane. A kinetic study of the isomerization of 5 to 6 and that between diastereomeric O-apical phosphoranes 13b-exo and 13b-endo revealed that 5b bearing an n-propylamino substituent at the central phosphorus atom was found to be less stable than the corresponding isomeric 6b by ca. 7.5 kcal mol(-1). This value was smaller than the difference in energy (11.9 kcal mol(-1)) between the O-equatorial (1b) and the O-apical n-butylphosphorane (2b) by 4.4 kcal mol(-1). This value of 4.4 kcal mol(-1) can be regarded as the stabilization energy induced by the n(N) --> sigma(P-O) interaction. The experimentally determined value was in excellent agreement with that derived from density functional theory (DFT) calculations at the B3PW91 level (4.0 kcal mol(-1)) between the nonsubstituted aminophosphoranes (5g is less stable than 6g by 10.1 kcal mol(-1)) and their P-methyl-substituted counterparts (1a is less stable than 2a by 14.1 kcal mol(-1)).