Ab initio, DFT and transition state theory calculations on 1,2-HF, HCI and CIF elimination reactions from CH2F-CH2Cl

This paper reports ab intio, DFT and transition state theory (TST) calculations on HF, HCI and CIF elimination reactions from CH2Cl-CH2F molecule. Both the ground state and the transition state for HX elimination reactions have been optimized at HF, MP2 and DFT calculations with 6-31G*, 6-31G** and 6-311++G** basis sets. In addition, CCSD(T) single point calculations were carried out with MP2/6-311++G** optimized geometry for more accurate determination of the energies of the minima and transition state, compared to the other methods employed here. Classical barriers are converted to Arrhenius activation energy by TST calculations for comparisons with experimental results. The pre-exponential factors, A, calculated at all levels of theory are significantly larger than the experimental values. For activation energy, E-a DFT gives good results for HF elimination, within 4-8 W mol(-1) from experimental values. None of the methods employed, including CCSD(T), give comparable results for HCI elimination reactions. However, rate constants calculated by CCSD(T) method are in very good agreement with experiment for HCI elimination and they are in reasonable agreement for HF elimination reactions. Due to the strong correlation between A and E., the rate constants could be fit to a lower A and E-a (as given by experimental fitting, corresponding to a tight TS) or to larger A and E-a (as given by high level ab initio calculations, corresponding to a loose TS). The barrier for CIF elimination is determined to be 607 U mol(-1) at HF level and it is unlikely to be important for CH2FCH2Cl. Results for other CH2X-CH2Y (X,Y = F/Cl) are included for comparison.

Aromaticity on the fly: cyclic transition state stabilization at finite temperature.

We study the transition state of pericyclic reactions at elevated temperature with unbiased ab initio molecular dynamics. We find that the transition state of the intramolecular rearrangements for barbaralane and bullvalene remains aromatic at high temperature despite the significant thermal atomic motions. Structural, magnetic, and electronic properties of the dynamical transition state show the concertedness and aromatic character. Free-energy calculations also support the validity of the transition state theory for the present rearrangement reactions. The calculations demonstrate that cyclic delocalization represents a strong force to synchronize the thermal atomic motions even at high temperatures.

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion.Including the configurational dependence will challenge the transition state theory of protein folding.

Synthesis and enzymatic evaluation of the guanosine analogue 2-amino-6-mercapto-7-methylpurine ribonucleoside (MESG): insights into the phosphorolysis reaction mechanism based on the blueprint transition state: SN1 or SN2?

A modified experimental procedure for the synthesis of MESG (2-amino-6-mercapto-7-methylpurine ribonucleoside) 1 has been successfully performed and its full characterization is presented. High resolution ESI(+)-MSMS indicates both the nucleoside bond cleavage as the main fragmentation in the gas phase and a possible SN1 mechanism. Ab initio transition state calculations based on the blue print transition state support this mechanistic rationale and discard an alternative SN2 mechanism. Assays using purine nucleoside phosphorylase (PNP) enzyme (human and M. tuberculosis sources) indicate its efficiency in the phosphorolysis of MESG and allow the quantitative determination of inorganic phosphate in real time assay.

Synthesis of unsaturated aminopyranosides as possible transition state mimics for glycosidases

Four unsaturated aminopyranosides have been prepared as possible transition-state mimics targeted towards carbohydrate processing enzymes. The conformations of the protonated aminosugars have been investigated by molecular modelling and their ability to inhibit alpha- and beta-glucosidases and an a-mannosidase have been probed. Two targets proved moderate inhibitors of alpha-glucosidases from Brewer's yeast and Bacillus stearothennophilus.

Mode-specificity and transition state-specific energy redistribution in the chemisorption of CH4 on Ni{100}

We have investigated methane (CH4) dissociative chemisorption on the Ni{100} surface by first-principles molecular dynamics (MD) simulations. Our results show that this reaction is mode-specific, with the n1 state being the most strongly coupled to efficient energy flow into the reaction coordinate when the molecule reaches the transition state. By performing MD simulations for two different transition state (TS) structures we provide evidence of TS structure-specific energy redistribution in methane chemisorption. Our results are compared with recently reported state-resolved measurement of methane adsorption probability on nickel surfaces, and we find that a strong correlation exists between the highest vibrational efficacy measured on Ni{100} for the n1 state and the calculated highest fractional vibrational energy content in this mode.

Angiotensin I-converting enzyme transition state stabilization by HIS1089 : evidence for a catalytic mechanism distinct from other gluzincin metalloproteinases

Angiotensin (Ang) I-converting enzyme (ACE) is a member of the gluzincin family of zinc metalloproteinases that contains two homologous catalytic domains. Both the N- and C-terminal domains are peptidyl-dipeptidases that catalyze Ang II formation and bradykinin degradation. Multiple sequence alignment was used to predict His¹⁰⁸⁹ as the catalytic residue in human ACE C-domain that, by analogy with the prototypical gluzincin, thermolysin, stabilizes the scissile carbonyl bond through a hydrogen bond during transition state binding. Site-directed mutagenesis was used to change His¹⁰⁸⁹ to Ala or Leu. At pH 7.5, with Ang I as substrate, k_cat/K_m values for these Ala and Leu mutants were 430 and 4,000-fold lower, respectively, compared with wild-type enzyme and were mainly due to a decrease in catalytic rate (k_cat) with minor effects on ground state substrate binding (K_m). A 120,000-fold decrease in the binding of lisinopril, a proposed transition state mimic, was also observed with the His¹⁰⁸⁹ --> Ala mutation. ACE C-domain-dependent cleavage of AcAFAA showed a pH optimum of 8.2. H1089A has a pH optimum of 5.5 with no pH dependence of its catalytic activity in the range 6.5-10.5, indicating that the His¹⁰⁸⁹ side chain allows ACE to function as an alkaline peptidyl-dipeptidase. Since transition state mutants of other gluzincins show pH optima shifts toward the alkaline, this effect of His¹⁰⁸⁹ on the ACE pH optimum and its ability to influence transition state binding of the sulfhydryl inhibitor captopril indicate that the catalytic mechanism of ACE is distinct from that of other gluzincins.

Evolutionary specialization of a tryptophan indole group for transition-state stabilization by eukaryotic transglutaminases

Covalent posttranslational protein modifications by eukaryotic transglutaminases proceed by a kinetic pathway of acylation and deacylation. Ammonia is released as the acylenzyme is formed, whereas the cross-linked product is released later in the deacylation step. Superposition of the active sites of transglutaminase type 2 (TG2) and the structurally related cysteine protease, papain, indicates that in the formation of tetrahedral intermediates, the backbone nitrogen of the catalytic Cys-277 and the N_Ɛ1 nitrogen of Trp-241 of TG2 could contribute to transition-state stabilization. The importance of this Trp-241 side chain was demonstrated by examining the kinetics of dansylcadaverine incorporation into a model peptide. Although substitution of the Trp-241 side chain with Ala or Gly had only a small effect on the Michaelis constant K_m (1.5-fold increase), it caused a >300-fold lowering of the catalytic rate constant k_cat. The wild-type and mutant TG2-catalyzed release of ammonia showed kinetics similar to the kinetics for the formation of cross-linked product, indicating that transitionstate stabilization in the acylation step was rate-limiting. In papain, a Gln residue is at the position of TG2-Trp-241. The conservation of Trp-241 in all eukaryotic transglutaminases and the finding that W241Q-TG2 had a much lower k_cat than wild-type enzyme suggest evolutionary specialization in the use of the indole group. This notion is further supported by the observation that transitionstate- stabilizing side chains of Tyr and His that operate in some serine and metalloproteases only partially substituted for Trp.

Configuration-dependent diffusion can shift the kinetic transition state and barrier height of protein folding

We show that diffusion can play an important role in protein-folding kinetics. We explicitly calculate the diffusion coefficient of protein folding in a lattice model. We found that diffusion typically is configuration- or reaction coordinate-dependent. The diffusion coefficient is found to be decreasing with respect to the progression of folding toward the native state, which is caused by the collapse to a compact state constraining the configurational space for exploration. The configuration- or position-dependent diffusion coefficient has a significant contribution to the kinetics in addition to the thermodynamic free-energy barrier. It effectively changes (increases in this case) the kinetic barrier height as well as the position of the corresponding transition state and therefore modifies the folding kinetic rates as well as the kinetic routes. The resulting folding time, by considering both kinetic diffusion and the thermodynamic folding free-energy profile, thus is slower than the estimation from the thermodynamic free-energy barrier with constant diffusion but is consistent with the results from kinetic simulations. The configuration- or coordinate-dependent diffusion is especially important with respect to fast folding, when there is a small or no free-energy barrier and kinetics is controlled by diffusion. Including the configurational dependence will challenge the transition state theory of protein folding. The classical transition state theory will have to be modified to be consistent. The more detailed folding mechanistic studies involving phi value analysis based on the classical transition state theory also will have to be modified quantitatively.