No description available.

The mechanism of three important aldol reactions and a biomimetic transamination is investigated using a combination of experimental kinetic isotope effects (KIEs), standard theoretical calculations and dynamics trajectory simulations. This powerful mechanistic probe is found to be invaluable in understanding intricate details of the mechanism of these reactions. The successful application of variational transition state theory including multidimensional tunneling to theoretically predict isotope effects, described in this dissertation, represents a significant advance in our research methodology. The role of dynamic effects in aldol reactions is examined in great detail. The study of the proline catalyzed aldol reaction has revealed an intriguing new dynamic effect - quasiclassical corner cutting - where reactive trajectories cut the corner between reactant and product valleys and avoid the saddle point. This phenomenon affects the KIEs observed in this reaction in a way that is not predictable by transition state theory. The study of the Roush allylboration of aldehydes presents an example where recrossing affects experimental observations. The comparative study of the allylboration of two electronically different aldehydes, which are predicted to have different amounts of recrossing, suggests a complex interplay of tunneling and recrossing affecting the observed KIEs. The Mukaiyama aldol reaction has been investigated and the results unequivocally rule out the key carbon-carbon bond forming step as rate-limiting. This raises several interesting mechanistic scenarios - an electron transfer mechanism with two different rate-limiting steps for the two components, emerges as the most probable possibility. Finally, labeling studies of the base catalyzed 1,3- proton transfer reaction of fluorinated imines point to a stepwise process involving an azomethine ylide intermediate. It is found that dynamic effects play a role in determining the product ratio in this reaction.

Organic reactions are systematically examined experimentally and theoretically to determine the role dynamics plays in the outcome of the reaction. It is shown that trajectory studies are of vital importance in understanding reactions influenced by dynamical motion. This dissertation discusses how a combination of kinetic isotope effects, theoretical calculations, and quasiclassical dynamics trajectories aid in the understanding of the solvolysis of p-tolyldiazonium cation in water, the cycloadditions of cyclopentadiene with diphenylketene and dichloroketene, and the cycloaddition of 2- methyl-2-butene with dichloroketene. In the solvolysis of p-tolyldiazonium cation, significant 13C kinetic isotope effects are qualitatively consistent with a transition state leading to formation of an aryl cation, but on a quantitative basis, the isotope effects are not adequately accounted for by simple SN1 heterolysis to the aryl cation. The best predictions of the 13C isotope effects for the heterolytic process arise from transition structures solvated by clusters of water molecules. Dynamic trajectories starting from these transition structures afford products very slowly. The nucleophilic displacement process for aryldiazonium ions in water is determined to be at the boundary of the SN2Ar and SN1 mechanisms. The reaction of cyclopentadiene with diphenylketene affords both [4 + 2] and [2 + 2] cycloadducts directly. This is surprising. There is only one low-energy transition structure for adduct formation. Investigation of this reaction indicates that quasiclassical trajectories started from a single transition structure afford both [4 + 2] and [2 + 2] products. Overall, an understanding of the products, rates, selectivities, isotope effects, and mechanism in these reactions requires the explicit consideration of dynamic trajectories.

Non-statistical recrossing in ketene cycloadditions with alkenes, heavy-atom tunneling and the mechanism of the decarboxylation of Mandelylthiamin is investigated in this dissertation. A combination of experimental kinetic isotope effects and theoretical models and kinetic isotope effects is utilized for this endeavor. This dissertation also describes how the use of quasiclassical dynamic trajectories, microcanonical RRKM calculations, and canonical variational transition state theory in combination with small-curvature tunneling approximations is utilized to help advance our research methodology to better understand mechanism. In the cycloaddition of dichloroketene with cis-2-butene, significant amounts of recrossing is observed using quasiclassical dynamic trajectories. An unusual inverse 13C intramolecular KIE lead us to investigate the role that heavy atoms play in non-statistical recrossing. More importantly, this discovery has uncovered a new phenomena of entropic intermediates that not only applies to ketene cycloadditions, but can also be applicable to other "concerted" reactions such as Diels-Alder reactions. The ring-opening of cyclopropylcarbinyl radical has revealed that heavy-atom tunneling plays a major role. The intramolecular 13C kinetic isotope effects for the ring-opening of cyclopropylcarbinyl radical were unprecedentedly large and in combination with theoretical predictions and multidimensional tunneling corrections, the role of tunneling in this reaction can be better understood. The mechanism decarboxylation of mandelylthiamin has been extensively studied in the literature. However, until the use of theoretically predicted KIEs and theoretical binding motifs the rate-limiting step of this reaction has been hotly debated. In this dissertation, a discussion of how the theoretical KIEs indicate the initial C-C bond as the rate-limiting step and chelating binding motifs of pyridinium and mandelylthiamin to explain the observed catalysis is given.

How the fast (femtosecond-picosecond, fs-ps) protein dynamics contribute to enzymatic function has gained popularity in modern enzymology. With multiple experimental and theoretical studies developed, the most challenging part is to assess both the chemical step kinetics and the relevant motions at the transition state (TS) on the fast time scale. Formate dehydrogenase (FDH), which catalyzes a single hydride transfer reaction, is a model system to address this specific issue. I have crystallized and solved the structure of FDH from Candida boidinii (CbFDH) in complex with NAD+ and azide. With the guidance of the structure information, two active site residues were identified, V123 and I175, which could be responsible for the narrow donor-acceptor-distance (DAD) distribution observed in the wild type CbFDH. This thesis describes studies using kinetic isotope effects (KIEs) and their temperature dependence together with two-dimensional infrared spectroscopy on the recombinant CbFDH and its V123 and I175 mutants. Those mutants were designed to systematically reduce the size of their side chain (I175V, I175A, V123A, V123G and double mutant I175V/V123A), leading to broader distribution of DADs. The kinetic experiments identified a correlation between the DAD distribution and the intrinsic KIEs. The contribution of the fs-ps dynamics was examined via two-dimensional infrared spectroscopy (2D IR) by measuring the vibrational relaxation of TS analog inhibitor, aizde, reflecting the TS environmental motions. Our results provide a test of models for the kinetics of the enzyme-catalyzed reaction that invokes motions of the enzyme at the fs-ps time scale to explain the temperature dependence of intrinsic KIEs.

Protein motions are complex, including occurring at different time scales, and their roles in enzyme-catalyzed reactions have always been of great interest among enzymologists. In order to characterize the potential factors that play a role on the chemical step of enzymatic reactions, variants of dihydrofolate reductase have been used as a benchmark system to study the motions of proteins correlated with the chemical step. A “global dynamic network” of coupled residues in Escherichia coli dihydrofolate reductase (ecDHFR), which assists in catalyzing the chemical step, has been demonstrated through quantum mechanical/molecular mechanical and molecular dynamic (QM/MM/MD) simulations, as well as bioinformatic analyses. A few specific residues — M42, G121, and I14 — were shown to function synergistically with measurements of single turnover rates and the temperature dependence of intrinsic kinetic isotope effects (KIEsint) of site-directed mutants. Although similar networks have been found in other enzymes, the general features of these networks are still unclear. This project focuses on exploring homologous residues of the proposed global network in human DHFR through computer simulations and measurements of the temperature dependence of KIEsint. The mutants M53W and S145V, both remote residues, showed significant decreases in catalytic efficiency. Non-additive isotope effects on activation energy were observed between M53 and S145, indicating their synergistic effect on hydride transfer in human DHFR. Apart from the effects of the conserved residues, we also extend our studies to exploring three potential phylogenetic events that account for the discrepancies between E. coli and human DHFR. They are L28, PP insertion and PEKN insertions by phylogenetic sequence analysis. Two of them (N23PP and G51PEKN E. coli DHFR) have been proved to be important both by MD simulation and experimental probe of KIEs measurement. The experiments have found that PP insertion itself rigidified the M20 loop and motions coupled to hydride transfer were impaired, however, loop rigidification was improved after incorporating PEKN. Furthermore, deletion of PP and PEKN of the engineered human enzyme also show a similar outcome. However, the effect of the key residue of L28 is not clear. In this project, we have step-wise engineered the human DHFR to be like hagfish (F31M) and E. coli (F32L). And it is found out that there is an increase in the temperature dependence of KIEs when the enzyme was bacterilized into a more primitive variant. This indicates that not only is residue F32 important and correlated with the chemical step as indicated by bioinformatic studies, but it is possible to trace the evolutionary trajectory. A triple mutation F32L-PP26N-PEKN62G on the human DHFR was also conducted, and it is not surprising to find out that the temperature dependence of KIEs has retained its behavior like wild-type human DHFR. These results suggest that the three predicted phylogenetically coherent events coevolved together to maintain the evolutionary preservation of the protein dynamics to enable H-tunneling from well-reorganized active sites. As has been indicated by the previous project, as the enzyme evolves, the active site of the enzyme will “reorganize” to form the optimal transition state for chemical step (from F32L-F32M-wild type DHFR). Here in this project, we aimed to systematically address this point of view through a series of cyclic permutation DHFR from directed evolutions. As this primitive enzyme is 7 orders of magnitude less efficient than the well-evolved human DHFR, together with four generations of evolved variants (cp, cp’ and cp”), this provides a good model system for explorations of the molecular basis of enzyme evolution. It is found that the organizations of transition state are improved before the catalytic efficiency is enhanced as the enzyme evolves.

The present dissertation seeks to deepen our understanding of hydrogen transfers and especially C-H bond activations in enzymes. Hydrogen transfers are ubiquitous in chemistry and biology and a thorough understanding of how they occur and what factors influence them will facilitate developments in biomimetic catalysis, rational drug design, and other fields. A particular difficulty with H-transfers is the importance of nuclear quantum effects to the reaction, particularly tunneling. The overall scope of the work here aims to examine how experimental kinetic isotope effects (KIEs) can be interpreted with a particular type of tunneling model, referred to as Marcus-like models, to yield a semi-quantitative picture of the physical mechanisms of H-transfers. Previous work had used this kind of model to qualitatively interpret experimental data using a combination of intuition and generalized theories. The work here examines these theories in quantitative detail, testing and calibrating our intuition in the context of several experimental systems. The first chapter of research (ch. II) focusses on the temperature dependence of primary KIEs and how these experiments can be quantitatively interpreted as a probe for certain kinds of enzyme or solvent dynamics. The subsequent chapters (ch. III-VI) focus on the use of secondary KIEs to determine the detailed structures of tunneling ready states (TRSs) and how the dynamics of H-tunneling affect those structures. These chapters focus primarily on the TRS of the enzyme alcohol dehydrogenase, but by examining an uncatalyzed analogue to that reaction (ch. VI), the work gains some insight about similarities and differences between catalyzed and uncatalyzed reactions. In summary, the work uncovers some principles of catalysis, not just the mechanism of a catalyzed reaction. The mechanism of C-H activation presented here provides an elegant solution to problems that have been vexing to accommodate within traditional models. This work constitutes some initial steps in making Marcus-like models quantitatively useful as a supplement or even replacement for traditional models of reactivity.

Choline oxidase (E.C. 1.1.3.17) from Arthrobacter globiformis catalyzes the four-electron oxidation of choline to glycine betaine (N,N,N-trimethylglycine) via two sequential, FAD-dependent reactions in which betaine aldehyde is formed as an enzyme-bound intermediate. In each oxidative half-reaction, molecular oxygen acts as electron acceptor and is converted into hydrogen peroxide. Biochemical, structural, and mechanistic studies on the wild-type and a number of mutant variants of choline oxidase have recently been carried out, allowing for the depiction of the mechanism of alcohol oxidation catalyzed by the enzyme. Catalysis by choline oxidase is initiated by the removal of the hydroxyl proton of alcohol substrate by a catalytic base in the enzyme-substrate complex, yielding the formation of the alkoxide species. In this dissertation, the roles of His351 and conserved His466 were investigated. The results presented demonstrate that His351 is involved in the stabilization of the transition state for the hydride transfer reaction and contributes to substrate binding. His466 is likely to be a catalytic base in choline oxidase due to its dramatic effect on enzymatic activity. Comparison of choline oxidase and other enzymes within its superfamily reveals the presence of a conserved His-Asn pair within the active site of enzymes. Therefore, the role of the conserved Asn510 in choline oxidase was examined in this study. The results presented here establish the importance of Asn510 in both the reductive and oxidative half-reactions. The lost of ability to form a hydrogen bond interaction between the side chain at position 510 with neighboring residues such as His466 resulted in a change from stepwise to concerted mechanism for the cleavages of OH and CH bonds of choline, as seen in the Asn510Ala mutant. Finally, the steady-state kinetic mechanism of pyranose 2-oxidase in the pH range from 5.5 to 8.5 was investigated. It was found that pH exerts significant effects on enzyme mechanism. This study has established the involvement of the residues in the initiation of enzyme catalysis and the stabilization of the alkoxide intermediate in choline oxidase. In addition, this work demonstrates the first instance in which the kinetic mechanism of a flavin-dependent oxidase is governed by pH.

The nucleotide 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) is phosphorylated twice to become a substrate for DNA polymerases, which copy a cell’s genetic information in advance of cell division. The main route to dTMP is mediated by the enzyme thymidylate synthase (TSase) and goes through 2'-deoxyuridine 5'-monophosphate (dUMP); dUMP’s heterocyclic aromatic pyrimidine ring loses a proton from its C5 position and gains a methylene and a hydride from the other reactant, methylene tetrahydrofolate (MTHF). In general, intricate knowledge of an enzyme’s mechanism can yield insight that leads to the development of precision-targeted inhibitors tailored exactly to thymidylate synthase. In fact, even more careful targeting could be achievable: Although E. coli TSase has served as a model system, investigators have increasingly been directing their lines of inquiry toward human TSase. A general enzymatic catalytic cascade is complex, comprising substrate binding, the chemical steps and product release; typically, the product release step is rate-limiting. TSase, however, is partially rate-limited by the chemistry portion of the process. The enzymatic mechanism has been considered for decades, yet recently has undergone a reassessment. After substrate binding – for which there is strong evidence for preference to dUMP as the first ligand in the wild-type E. coli enzyme – the important events are methylene transfer from MTHF to dUMP, proton abstraction and hydride transfer. The last of these – hydride transfer – is irreversible and rate-limiting (to a large degree without Mg2+, and to a small but noticeable degree with Mg2+). The studies described here are aimed at three therapeutically relevant questions: (a) determining the extent of negative charge accumulation at the O4 position of the hydride transfer acceptor; (b) expanding knowledge of the differential properties of E. coli and human TSase; and (c) gaining insight into the molecular origin of the drug resistance seen in a clinically relevant human TSase mutant. The properties touched on in this work include steady-state kinetics; inhibition constants toward 5-fluoro dUMP, substrate binding sequence and the temperature dependency of intrinsic hydride transfer kinetic isotope effects (KIEs). Intrinsic KIEs are a specialized measurement that permits the investigator to examine a particular hydrogen transfer step in isolation; it is achieved by labeling the bond to hydrogen broken in the reaction with protium (1H, also written as H), deuterium (2H, also written as D) or tritium (3H, also written as T). The latter is radioactive. The reaction is conducted with a mixture of two hydrogen isotopes at a time, and the extent to which the heavier isotope is disfavored against reaction is assessed; this covers multiple steps. Heavier isotopes directly participating in a chemical step react slower both because of zero-point vibrational energies if a semi-classical view is taken and because of the mass-dependence of tunneling probabilities if a quantum-mechanical view is taken. Each of the two-way isotopic comparisons mentioned above furnishes an observed KIE for that competition between two isotopes. Mathematical combination of two isotopic comparisons cancels out the effect of isotopically insensitive steps and provides rich insight into the hydride transfer alone. The ultimate result is the ratio of rate constants for the isotopologues; this ratio’s magnitude and variation with temperature report on the compactness of the active site and its resistance to thermal fluctuation, respectively. Our results reveal a possible role for E. coli asparagine 177 (N177) in the hydride transfer transition state (TS) stabilization, as revealed by its disruption in the aspartate mutant, N177D. This disruption was found to be alleviated to a high extent when the substrate was changed to dCMP, consistent with the N177 stabilizing partial negative charge at the TS for hydride transfer. This has drug design implications. Our work on human TSase underscores slightly weaker substrate binding preference, insensitivity to Mg2+ and mild alteration of hydride transfer TS when compared with E. coli TSase. Finally, analysis of the Y33H mutant of human TSase – the affected residue being remote from the active site – indicated the drug resistance was because of a higher inhibition constant for 5F-dUMP and that the hydride transfer step is disrupted, with a wider variation among donor-acceptor distances (between the two carbons involved in the hydride transfer at the TS for that step). Other researchers’ crystallographic evidence reveals greater positional uncertainty for a set of active-site side chains in the E. coli equivalent mutant. In totality, the data available implicate enzyme motions as relevant to drug binding and to catalysis for human TSase. In summary, the research described herein enriches the understanding of several aspects of the behavior of multiple TSase variants – the overall performance as seen via steady-state kinetics; the pattern of substrate binding as seen with observed KIEs for the proton abstraction step; and the efficiency of active site preparation for hydride transfer as evidenced in the temperature dependency of intrinsic hydride transfer KIEs.

Type 2 diabetes is a major health concern worldwide. One of its complications is postprandial hyperglycemia, i.e., high blood glucose concentrations, caused by glucose fast release from dietary polysaccharides into the bloodstream after meals. α-Glucosidase inhibitor drugs reduce postprandial hyperglycemia by inhibiting maltase/glucoamylase (MGAM) and sucrase/isomaltase (SI). MGAM and SI transform polysaccharides into absorbable monosaccharides, and inhibiting them delays monosaccharide release into the blood. The three commercially available α-glucosidase inhibitors are limited by their absorption abilities, inhibition efficacies, and side effects, which highlights the need for more specific α-glucosidase inhibitors. Because enzymes catalyze their reactions by tightly binding to their cognate transition states (TS), TS analogs can be powerful inhibitors and potential drugs. The measurement and interpretation of kinetic isotope effects (KIEs) is the only method that can directly determine TS structures on large molecules. In this work, methods to prepare radioisotopically labelled maltose were developed, as well as methods to measure KIEs on acid- and enzyme-catalyzed maltose hydrolysis. However, the methods developed did not achieve the required precision for TS analysis. Also, KIEs were calculated computationally for a model reaction of maltose hydrolysis. / Thesis / Master of Science (MSc)