Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Steady state optimization of an ammonia reactor

A hybrid simulation technique for identification and steady state optimization of a tubular reactor used in ammonia synthesis is presented. The parameter identification program finds the catalyst activity factor and certain heat transfer coefficients that minimize the sum of squares of deviation from simulated and actual temperature measurements obtained from an operating plant. The optimization program finds the values of three flows to the reactor to maximize the ammonia yield using the estimated parameter values. Powell's direct method of optimization is used in both cases. The results obtained here are compared with the plant data.

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.

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.