Communication: Transition State Theory for dissipative systems without a dividing surface

Transition state theory is a central cornerstone in reaction dynamics. Its key step is the identification of a dividing surface that is crossed only once by all reactive trajectories. This assumption is often badly violated, especially when the reactive system is coupled to an environment. The calculations made in this way then overestimate the reaction rate and the results depend critically on the choice of the dividing surface. In this Communication, we study the phase space of a stochastically driven system close to an energetic barrier in order to identify the geometric structure unambiguously determining the reactive trajectories, which is then incorporated in a simple rate formula for reactions in condensed phase that is both independent of the dividing surface and exact.

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.

Quantum transition state theory and solving a first order master equation for complex reactions numerically

The field of chemical kinetics is an exciting and active field. The prevailing theories make a number of simplifying assumptions that do not always hold in actual cases. Another current problem concerns a development of efficient numerical algorithms for solving the master equations that arise in the description of complex reactions. The objective of the present work is to furnish a completely general and exact theory of reaction rates, in a form reminiscent of transition state theory, valid for all fluid phases and also to develop a computer program that can solve complex reactions by finding the concentrations of all participating substances as a function of time. To do so, the full quantum scattering theory is used for deriving the exact rate law, and then the resulting cumulative reaction probability is put into several equivalent forms that take into account all relativistic effects if applicable, including one that is strongly reminiscent of transition state theory, but includes corrections from scattering theory. Then two programs, one for solving complex reactions, the other for solving first order linear kinetic master equations to solve them, have been developed and tested for simple applications.

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.

Immucillin H, a powerful transition-state analog inhibitor of purine nucleoside phosphorylase, selectively inhibits human T lymphocytes

Transition-state theory has led to the design of Immucillin-H (Imm-H), a picomolar inhibitor of purine nucleoside phosphorylase (PNP). In humans, PNP is the only route for degradation of deoxyguanosine, and genetic deficiency of this enzyme leads to profound T cell-mediated immunosuppression. This study reports the biological effects and mechanism of action of Imm-H on malignant T cell lines and on normal activated human peripheral T cells. Imm-H inhibits the growth of malignant T cell leukemia lines with the induction of apoptosis. Imm-H also inhibits activated normal human T cells after antigenic stimulation in vitro. However, Imm-H did not inhibit malignant B cells, colon cancer cell lines, or normal human nonstimulated T cells, demonstrating the selective activity of Imm-H. The effects on leukemia cells were mediated by the cellular phosphorylation of deoxyguanosine and the accumulation of dGTP, an inhibitor of ribonucleotide diphosphate reductase. Cells were protected from the toxic effects of Imm-H when deoxyguanosine was absent or when deoxycytidine was present. Guanosine incorporation into nucleic acids was selectively blocked by Imm-H with no effect on guanine, adenine, adenosine, or deoxycytidine incorporation. Imm-H may have clinical potential for treatment of human T cell leukemia and lymphoma and for other diseases characterized by abnormal activation of T lymphocytes. The design of Imm-H from an enzymatic transition-state analysis exemplifies a powerful approach for developing high-affinity enzyme inhibitors with pharmacologic activity.

Non-enzymatic and enzymatic hydrolysis of alkyl halides: A haloalkane dehalogenation enzyme evolved to stabilize the gas-phase transition state of an SN2 displacement reaction

The semiempirical PM3 method, calibrated against ab initio HF/6–31+G(d) theory, has been used to elucidate the reaction of 1,2-dichloroethane (DCE) with the carboxylate of Asp-124 at the active site of haloalkane dehalogenase of Xanthobacter autothropicus. Asp-124 and 13 other amino acid side chains that make up the active site cavity (Glu-56, Trp-125, Phe-128, Phe-172, Trp-175, Leu-179, Val-219, Phe-222, Pro-223, Val-226, Leu-262, Leu-263, and His-289) were included in the calculations. The three most significant observations of the present study are that: (i) the DCE substrate and Asp-124 carboxylate, in the reactive ES complex, are present as an ion-molecule complex with a structure similar to that seen in the gas-phase reaction of AcO− with DCE; (ii) the structures of the transition states in the gas-phase and enzymatic reaction are much the same where the structure formed at the active site is somewhat exploded; and (iii) the enthalpies in going from ground states to transition states in the enzymatic and gas-phase reactions differ by only a couple kcal/mol. The dehalogenase derives its catalytic power from: (i) bringing the electrophile and nucleophile together in a low-dielectric environment in an orientation that allows the reaction to occur without much structural reorganization; (ii) desolvation; and (iii) stabilizing the leaving chloride anion by Trp-125 and Trp-175 through hydrogen bonding.

Partial molar volume, surface area, and hydration changes for equilibrium unfolding and formation of aggregation transition state: High-pressure and cosolute studies on recombinant human IFN-γ

The equilibrium dissociation of recombinant human IFN-γ was monitored as a function of pressure and sucrose concentration. The partial molar volume change for dissociation was −209 ± 13 ml/mol of dimer. The specific molar surface area change for dissociation was 12.7 ± 1.6 nm2/molecule of dimer. The first-order aggregation rate of recombinant human IFN-γ in 0.45 M guanidine hydrochloride was studied as a function of sucrose concentration and pressure. Aggregation proceeded through a transition-state species, N*. Sucrose reduced aggregation rate by shifting the equilibrium between native state (N) and N* toward the more compact N. Pressure increased aggregation rate through increased solvation of the protein, which exposes more surface area, thus shifting the equilibrium away from N toward N*. The changes in partial molar volume and specific molar surface area between the N* and N were −41 ± 9 ml/mol of dimer and 3.5 ± 0.2 nm2/molecule, respectively. Thus, the structural change required for the formation of the transition state for aggregation is small relative to the difference between N and the dissociated state. Changes in waters of hydration were estimated from both specific molar surface area and partial molar volume data. From partial molar volume data, estimates were 25 and 128 mol H2O/mol dimer for formation of the aggregation transition state and for dissociation, respectively. From surface area data, estimates were 27 and 98 mol H2O/mol dimer. Osmotic stress theory yielded values ≈4-fold larger for both transitions.

State-to-state reaction probabilities within the quantum transition state framework

Rigorous quantum dynamics calculations of reaction rates and initial state-selected reaction probabilities of polyatomic reactions can be efficiently performed within the quantum transition state concept employing flux correlation functions and wave packet propagation utilizing the multi-configurational time-dependent Hartree approach. Here, analytical formulas and a numerical scheme extending this approach to the calculation of state-to-state reaction probabilities are presented. The formulas derived facilitate the use of three different dividing surfaces: two dividing surfaces located in the product and reactant asymptotic region facilitate full state resolution while a third dividing surface placed in the transition state region can be used to define an additional flux operator. The eigenstates of the corresponding thermal flux operator then correspond to vibrational states of the activated complex. Transforming these states to reactant and product coordinates and propagating them into the respective asymptotic region, the full scattering matrix can be obtained. To illustrate the new approach, test calculations study the D + H2(ν, j) → HD(ν′, j′) + H reaction for J = 0.

Electrochemical study and dependence of 'transition state' in Co(II) and Ni(II) complexes with some antibiotics and cephalothin

Electrode kinetics and study of 'transition state' with applied potential in case of [M - antibiotics - cephalothin] system were reported at pH = 7.30 ± 0.01 at suitable supporting electrolyte at 25.0ºC. The M = Co or Ni and antibiotics were doxycycline, chlortetracycline, oxytetracycline, tetracycline, minocycline, amoxicillin and chloramphenicol used as primary ligands and cephalothin as secondary ligand. Kinetic parameters viz. transfer coefficient (a), degree of irreversibility (l), diffusion coefficient (D) and rate constant (k) were determined. The values of a and k varied from 0.41 to 0.59 and 2.60 X 10-3 cm s-1 to 9.67 X 10-3 cm s-1 in case of [Co - antibiotics - cephalothin] system. In case of [Ni - antibiotics - cephalothin], a and k varied from 0.41 to 0.58 and 2.34 X 10-3 cm s-1 to 9.19 X 10-3 cm s-1 respectively confirmed that transition state behaves between oxidant and reductant response to applied potential and it adjusts it self in such a way that the same is located midway between dropping mercury electrode and solution interface. The values of rate constant confirmed the quasireversible nature of electrode processes. The stability constants (logb) of complexes were also determined.

Theory of Breaking of a Potyrner Molecule under Tension

The thesis presents the dynamics of a polymer chain under tension. It includes existing theories of polymer fracture, important theories of reaction rates, the rate using multidimensional transition state theory and apply it to the case of polyethylene etc. The main findings of the study are; the life time of the bond is somewhat sensitive to the potential lead to rather different answers, for a given potential a rough estimate of the rate can be obtained by a simples approximation that considers the dynamics of only the bond that breaks and neglects the coupling to neighboring bonds. Dynamics of neighboring bonds would decrease the rate, but usually not more than by one order of magnitude, for the breaking of polyethylene, quantum effects are important only for temperatures below 150K, the lifetime strongly depends on the strain and as the strain varies over a narrow range, the life varies rapidly from 105 seconds to 10_5 seconds, if we change one unit of the polymer by a foreign atom, say by one sulphure atom, in the main chain itself, by a weaker bond, the rate is found to increase by orders of magnitude etc.

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.