Transition state structure, determination

Transition State Structure Determination from Multiple Kinetic Isotope... 

Williams 1 H 1993. Interplay of Theory and Experiment in the Determination of Transition-State Structures. Chemical Society Reviews 1 277-283. 

Another use of frequency calculations is to determine the nature of a stationary point found by a geometry optimization. As we ve noted, geometry optimizations converge to a structure on the potential energy surface where the forces on the system are essentially zero. The final structure may correspond to a minimum on the potential energy surface, or it may represent a saddle point, which is a minimum with respect to some directions on the surface and a maximum in one or more others. First order saddle points—which are a maximum in exactly one direction and a minimum in all other orthogonal directions—correspond to transition state structures linking two minima. 

Facial selectivity in 1,3-dipolar cycloadditions to cis-3,4-dimethylcyclobutene (73) (Scheme 1.21) was studied. Only phenylglyoxylo- and pyruvonitrile oxides lacked facial selectivities (anti syn = 1 1). All other nitrile oxides formed preferably anti-74. The anti/syn ratio increased from 60 40 (R = P-O2NC6K4) and 65 35 (R = Ph) to 87 13 and 92 8 for bulky ten-Bu and mesityl substituents, respectively. The transition-state structure of the cycloaddition of formonitrile oxide was determined using both HF/6-31G and B3LYP/6-31G methods. The... 

In both solvents, the variational transition state (associated with the free energy maximum) corresponds, within the numerical errors, to the dividing surface located at rc = 0. It has to be underlined that this fact is not a previous hypothesis (which would rather correspond to the Conventional Transition State Theory), but it arises, in this particular case, from the Umbrella Sampling calculations. However, there is no information about which is the location of the actual transition state structure in solution. Anyway, the definition of this saddle point has no relevance at all, because the Monte Carlo simulation provides directly the free energy barrier, the determination of the transition state structure requiring additional work and being unnecessary and unuseful. 

Wolfe and Kim (1991) also reported that the magnitude of a secondary a-deuterium KIE is primarily determined by the changes that occur in the Ca—H(D) stretching vibrations when the reactant is converted into the transition state. Wolfe and Kim calculated the transition state structures and the secondary a-deuterium KIEs for a series of identity SN2 reactions of methyl substrates [reaction (8)] at various levels of theory ranging from 4-31G to MP4/6-31 + G //6-31 + G. The KIEs were partitioned into two contributions, those from the Ca—H(D) stretching vibrations and those from the Ca—H(D) bending vibrations. 

Although the relationship between the secondary a-deuterium KIEs and transition state structure is different for the two types of transition state and interpreting secondary a-deuterium KIEs is, therefore, more difficult, it appears that the change in the KIE with substituent should be a good indicator for determining whether an SN2 transition state is symmetrical or unsymmetrical. 

Finally, as is the case for the secondary a-deuterium KIEs, the /3-deuterium KIE is assumed to vary in magnitude from near unity for a reactant-like transition state to a maximal value for a transition state resembling the carbocation formed in an SN1 reaction. The experimentally determined KIE may, therefore, be used as a measure of transition state structure provided that the maximum value of the KIE, i.e. the EIE for the formation of the carbocation, is known. 

T-secondary isotope effect can be determined. As recounted in the last item of Chart 3, such effects are expected to be measures of transition-state structure. If the transition state closely resembled reactants, then no change in the force field at the isotopic center would occur as the reactant state is converted to the transition state and the -secondary kinetic isotope effect should be 1.00. If the transition state closely resembled products, then the transition-state force field at the isotopic center would be very similar to that in the product state, and the a-secondary kinetic isotope effect should be equal to the equilibrium isotope effect, shown by Cook, Blanchard, and Cleland to be 1.13. Between these limits, the kinetic isotope effect should change monotonically from 1.00 to 1.13. 

One of the important consequences of studying catalysis by mutant enzymes in comparison with wild-type enzymes is the possibility of identifying residues involved in catalysis that are not apparent from crystal structure determinations. This has been usefully applied (Fersht et al., 1988) to the tyrosine activation step in tyrosine tRNA synthetase (47) and (49). The residues Lys-82, Arg-86, Lys-230 and Lys-233 were replaced by alanine. Each mutation was studied in turn, and comparison with the wild-type enzyme revealed that each mutant was substantially less effective in catalysing formation of tyrosyl adenylate. Kinetic studies showed that these residues interact with the transition state for formation of tyrosyl adenylate and pyrophosphate from tyrosine and ATP and have relatively minor effects on the binding of tyrosine and tyrosyl adenylate. However, the crystal structures of the tyrosine-enzyme complex (Brick and Blow, 1987) and tyrosyl adenylate complex (Rubin and Blow, 1981) show that the residues Lys-82 and Arg-86 are on one side of the substrate-binding site and Lys-230 and Lys-233 are on the opposite side. It would be concluded from the crystal structures that not all four residues could be simultaneously involved in the catalytic process. Movement of one pair of residues close to the substrate moves the other pair of residues away. It is therefore concluded from the kinetic effects observed for the mutants that, in the wild-type enzyme, formation of the transition state for the reaction involves a conformational change to a structure which differs from the enzyme structure in the complex with tyrosine or tyrosine adenylate. The induced fit to the transition-state structure must allow interaction with all four residues simultaneously. 

Both equilibrium and transition-state structure may be determined from calculation. The former requires a search for an energy minimum on a potential energy surface while the latter requires a search for an energy maximum. Lifetime or even existence is not a requirement. 

If the fragmentation of the alkoxide is rate-determining and endothermic, one would expect the transition state structure to resemble the products and thence be influenced by the relative stability of the carbonyl fragment. A transformation similar to the retroaldol cleavage of B-hydroxynitrosamines has been observed. Kruger reported that the treatment of 2-ketopropylpropylnitrosa-mine with refluxing potassium hydroxide in alcohol led to the production of methyl propylnitrosamine as is illustrated in equation 6 (10). This transformation is also under study in our laboratory. 

Figure 15.5 Transition-state structure for rate-determining hydrogen atom transfer in the methane metathesis reaction of methyllutetiocene. Note that the kinetics for this narcissistic reaction may be followed by using a label either in the reacting methane or in the methyl group of the starting organometallic...

Stable adsorption complexes are characterized by local minima on the potential energy hypersurface. The reaction pathway between two stable minima is determined by computation of a transition state structure, a saddle point on the potential energy hypersurface, characterized by a single imaginary vibrational mode. The Cartesian displacements of atoms that participate in this vibration characterize movements of these atoms along the reaction coordinate between sorption complexes. 

We discuss in this paper the unconstrained optimization of stationary points of a smooth function fix) in many variables. The emphasis is on methods useful for calculating molecular electronic energies and for determining molecular equilibrium and transition state structures. The discussion is general and practical aspects concerning computer implementations are not treated. 

---