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.

Buckled nano rod - a two state system: quantum effects on its dynamics

We consider a suspended elastic rod under longitudinal compression. The compression can be used to adjust potential energy for transverse displacements from the harmonic to the double well regime. The two minima in potential energy curve describe two possible buckled states. Using transition state theory (TST) we have calculated the rate of conversion from one state to other. If the strain epsilon = 4 epsilon c the simple TST rate diverges. We suggest a method to correct this divergence for quantum calculations. We also find that zero point energy contributions can be quite large so that single mode calculations can lead to large errors in the rate.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.

How do proteins fold?

Understanding the mechanism by which an unfolded polypeptide chain folds to its unique, functional structure is a primary unsolved problem in biochemistry. Fundamental advances towards understanding how proteins fold have come from kinetic studies, Kinetic studies allow the dissection of the folding pathway of a protein into individual steps that are defined by partially-structured folding intermediates. Improvements in both the structural and temporal resolution of physical methods that are used to monitor the folding process, as well as the development of new methodologies, are now making it possible to obtain detailed structural information on protein folding pathways. The protein engineering methodology has been particularly useful in characterizing the structures of folding intermediates as well as the transition state of folding, Several characteristics of protein folding pathways have begun to emerge as general features for the folding of many different proteins. Progress in our understanding of how structure develops during folding is reviewed here.

A Theoretical Investigation of the Ni(II)-Catalyzed Hydrovinylation of Styrene

We report a detailed and full computational investigation on the hydrovinylation reaction of styrene with the Ni(II)-phospholane catalytic system, which was originally presumed to proceed through a cationic mechanism involving a nickel hydride intermediate. The following general features emerge from this study on a specific catalyst complex that was found to give quantitative yield and moderate selectivity: (a) the activation barrier for the initiation (18.8 kcal/mol) is higher than that for the reaction due to a low-lying square-planar pentenyl chelate intermediate originating from a Ni(II)-allyl catalyst precursor. Consequently there is an induction period for the catalysis; (b) the exit of product from the catalyst is via a β-H-transfer step instead of the usual β-H elimination pathway, which has a very high activation energy due to a trans effect of the phospholane ligand; (c) the turnover-limiting and enantio- determining transition state is also the β-H-transfer; (d) because of the absence of a hydride intermediate, the unwanted isomerization of the product is prevented; (e) since the enantio-discrimination is decided at the H-transfer stage itself, the configuration of the product in a catalytic cycle influences the enantioselectivity in the subsequent cycle; (f) the trans effect of the sole strong ligand in the d8 square-planar Ni(II), the stability of the η3-benzyl intermediate, and the availability of three coordination sites enable regioselective hydrovinylation over the possible oligomerization/polymerization of the olefin substrates and linear hydrovinylation. This work has also confirmed the previously recognized role of the hemilabile group at various stages in the mechanism.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

A lower bound to the survival probability and an approximate first passage time distribution for Markovian and non-Markovian dynamics in phase space

We derive a very general expression of the survival probability and the first passage time distribution for a particle executing Brownian motion in full phase space with an absorbing boundary condition at a point in the position space, which is valid irrespective of the statistical nature of the dynamics. The expression, together with the Jensen's inequality, naturally leads to a lower bound to the actual survival probability and an approximate first passage time distribution. These are expressed in terms of the position-position, velocity-velocity, and position-velocity variances. Knowledge of these variances enables one to compute a lower bound to the survival probability and consequently the first passage distribution function. As examples, we compute these for a Gaussian Markovian process and, in the case of non-Markovian process, with an exponentially decaying friction kernel and also with a power law friction kernel. Our analysis shows that the survival probability decays exponentially at the long time irrespective of the nature of the dynamics with an exponent equal to the transition state rate constant.

Kinetics and the Mechanism of Interaction of the Endoplasmic Reticulum Chaperone, Calreticulin, with Monoglucosylated (Glc1Man9GlcNAc2) Substrate

Calreticulin is a lectin-like molecular chaperone of the endoplasmic reticulum in eukaryotes. Its interaction with N-glycosylated polypeptides is mediated by the glycan, Glc(1)Man(9)GlcNAc(2), present on the target glycoproteins. In this work, binding of monoglucosyl IgG (chicken) substrate to calreticulin has been studied using real time association kinetics of the interaction with the biosensor based on surface plasmon resonance (SPR). By SPR, accurate association and dissociation rate constants were determined, and these yielded a micromolar association constant. The nature of reaction was unaffected by immobilization of either of the reactants. The Scatchard analysis values for K-a agreed web crith the one obtained by the ratio k(1)/k(-1). The interaction was completely inhibited by free oligosaccharide, Glc(1)Man(9)GlcNAc(2), whereas Man(9)GlcNAc(2) did not bind to the calreticulin-substrate complex, attesting to the exquisite specificity of this interaction. The binding of calreticulin to IgG was used for the development of immunoassay and the relative affinity of the lectin-substrate association was indirectly measured. The values are in agreement with those obtained with SPR. Although the reactions are several orders of magnitude slower than the diffusion controlled processes, the data are qualitatively and quantitatively consistent with single-step bimolecular association and dissociation reaction. Analyses of the activation parameters indicate that reaction is enthalpically driven and does not involve a highly ordered transition state. Based on these data, the mechanism of its chaperone activity is briefly discussed.

Stereoselectivity in the addition of hydrogen chloride to cyclohex-1-enecarbonitrile: Possible role of chlorine lone pair-nitrile π-orbital pseudoallylic A1, 3 interaction

The (overall trans) addition of hydrogen chloride to cyclohex-1- enecarbonitrile in anhydrous alcoholic media proceeds to give cis-2-chlorocyclohexanecarboxylate (together with some cis-2- chlorocyclohexanecarboxamide): no corresponding products with the trans-configuration are detectable. In anhydrous ether the addition proceeds to give a single isomer, presumably cis-, of 2-chlorocyclohexanecarbonitrile, indicating that the configuration of the products may not be equilibrium-controlled in alcoholic media. An examination of the steric factors indicates that the transition state for protonation of the presumed intermediate, 2-chlorocyclohexylidenemethylideneimine, leading to cis-product is favoured if interaction between the lateral π-orbital of the C-N double bond and the lone-pairs on the chlorine atom at the 2-position is large. Consideration of interactions in the transition states meets Zimmerman's criticism that invoking A1, 3 interaction existing in ground states to explain product configuration takes insufficient account of the Curtin-Hammett principle.

Kinetics and the mechanism of interaction of the endoplasmic reticulum chaperone, calreticulin, with monoglucosylated (Glc(1)Man(9)GlcNAc(2)) substrate

Calreticulin is a lectin-like molecular chaperone of the endoplasmic reticulum in eukaryotes. Its interaction with N-glycosylated polypeptides is mediated by the glycan, Glc(1)Man(9)GlcNAc(2), present on the target glycoproteins. In this work, binding of monoglucosyl IgG (chicken) substrate to calreticulin has been studied using real time association kinetics of the interaction with the biosensor based on surface plasmon resonance (SPR). By SPR, accurate association and dissociation rate constants were determined, and these yielded a micromolar association constant. The nature of reaction was unaffected by immobilization of either of the reactants. The Scatchard analysis values for K-a agreed web crith the one obtained by the ratio k(1)/k(-1). The interaction was completely inhibited by free oligosaccharide, Glc(1)Man(9)GlcNAc(2), whereas Man(9)GlcNAc(2) did not bind to the calreticulin-substrate complex, attesting to the exquisite specificity of this interaction. The binding of calreticulin to IgG was used for the development of immunoassay and the relative affinity of the lectin-substrate association was indirectly measured. The values are in agreement with those obtained with SPR. Although the reactions are several orders of magnitude slower than the diffusion controlled processes, the data are qualitatively and quantitatively consistent with single-step bimolecular association and dissociation reaction. Analyses of the activation parameters indicate that reaction is enthalpically driven and does not involve a highly ordered transition state. Based on these data, the mechanism of its chaperone activity is briefly discussed.

Asynchronous behaviour in some three-atom linear concerted radical-exchange reactions

A few simple three-atom thermoneutral radical exchange reactions (i.e. A + BC --> AB + C) are examined by ab initio SCF methods. Emphasis is laid on the detailed analysis of density matrices rather than on energetics. Results reveal that the sum of the bond orders of the breaking and forming bonds is not conserved to unity, due to development of free valence on the migrating atom 'B' in the transition state. Bond orders, free valence and spin densities on the atoms are calculated. The present analysis shows that the bond-cleavage process is always more advanced than the bond-formation process in the transition state. Further analysis shows a development of the negative spin density on the migrating atom 'B' in the transition state. The depletion of the alpha-spin density on the radical site "A" in the reactant during the reaction lags behind the growth of the alpha-spin density on the terminal atom "C" of the reactant bond, 'B-C' in the transition state. But all these processes are completed simultaneously at the end of the reaction. Hence, the reactions are asynchronous but kinetically concerted in most cases.

Computer modelling studies on the mechanism of action of ribonuclease T1

The mechanism of action of ribonuclease (RNase) T1 is still a matter of considerable debate as the results of x-ray, 2-D nmr and site-directed mutagenesis studies disagree regarding the role of the catalytically important residues. Hence computer modelling studies were carried out by energy minimisation of the complexes of RNase T1 and some of its mutants (His40Ala, His40Lys, and Glu58Ala) with the substrate guanyl cytosine (GpC), and of native RNase T1 with the reaction intermediate guanosine 2',3'-cyclic phosphate (G greater than p). The puckering of the guanosine ribose moiety in the minimum energy conformer of the RNase T1-GpC (substrate) complex was found to be O4'-endo and not C3'-endo as in the RNase T1-3'-guanylic acid (inhibitor/product) complex. A possible scheme for the mechanism of action of RNase T1 has been proposed on the basis of the arrangement of the catalytically important amino acid residues His40, Glu58, Arg77, and His92 around the guanosine ribose and the phosphate moiety in the RNase T1-GpC and RNase T1-G greater than p complexes. In this scheme, Glu58 serves as the general base group and His92 as the general acid group in the transphosphorylation step. His40 may be essential for stabilising the negatively charged phosphate moiety in the enzyme-transition state complex.

The Mechanism Of The Hydrogen Abstraction Process By Photo-Excited Ketones

Qualitative potential energy surfaces for hydrogen abstraction from alkanes containing primary, secondary and tertiary C-H bonds by a photo-excited ketone have been reported, The results suggest that the activation barriers for these processes decrease in the order primary > secondary > tertiary in agreement with the observed trend in the rate constants. The analysis of the electronic structure of the transition-state reveal that electron-transfer from hydrocarbon to ketone and formation of a new bond are almost synchronous in the hydrogen transfer process. The tunneling of hydrogen is not important in the normal temperature region even though the barriers are small.

pi.-Facial selectivities in cycloadditions to norbornyl- and norbornenyl-fused p-benzoquinones

The stereochemistry of the Diels-Alder cycloaddition of several dienes to the facially perturbed dienophiles 2,3-norbornenobenzoquinone (3) and 2,3-norbornanobenzoquinone (4) has been examined. Unambiguous structural proof for the adducts formed has been obtained from complementary 'H and I3C NMR spectral data and in two cases through X-ray crystal structure determination. While 1,3-~yclopentadiene1, ,3-~yclohexadienea, nd cyclooctatetraene exhibit preference for addition to 3 from the bottom side, the stereochemical outcome is reversed in their response to 4.1,3-DiphenyIisobenzofuran and 1,2,3,4-tetrachloro-5,5-dimethoxycyclopentadieenneg aged 3 from the top side with marked selectivity, which is further enhanced in their reaction with 4. The observed stereoselectivities seem to be essentially controlled by steric interactons at the transition state. Model calculations provide support for this interpretation.