Transition State Analysis Using Multiple Kinetic Isotope Effects: Mechanisms of Enzymatic and Non-enzymatic Glycoside Hydrolysis and Transfer

Clostridioides difficile TcdB Toxin Glucosylates Rho GTPase by an SNi Mechanism and Ion Pair Transition State

Toxins TcdA and TcdB from Clostridioides difficile glucosylate human colon Rho GTPases. TcdA and TcdB glucosylation of RhoGTPases results in cytoskeletal changes, causing cell rounding and loss of intestinal integrity. Clostridial toxins TcdA and TcdB are proposed to catalyze glucosylation of Rho GTPases with retention of stereochemistry from UDP-glucose. We used kinetic isotope effects to analyze the mechanisms and transition-state structures of the glucohydrolase and glucosyltransferase activities of TcdB. TcdB catalyzes Rho GTPase glucosylation with retention of stereochemistry, while hydrolysis of UDP-glucose by TcdB causes inversion of stereochemistry. Kinetic analysis revealed TcdB glucosylation via the formation of a ternary complex with no intermediate, supporting an S_Ni mechanism with nucleophilic attack and leaving group departure occurring on the same face of the glucose ring. Kinetic isotope effects combined with quantum mechanical calculations revealed that the transition states of both glucohydrolase and glucosyltransferase activities of TcdB are highly dissociative. Specifically, the TcdB glucosyltransferase reaction proceeds via an S_Ni mechanism with the formation of a distinct oxocarbenium phosphate ion pair transition state where the glycosidic bond to the UDP leaving group breaks prior to attack of the threonine nucleophile from Rho GTPase.

Matching Glycosyl Donor Reactivity to Sulfonate Leaving Group Ability Permits SN2 Glycosylations

Here we demonstrate that highly β-selective glycosylation reactions can be achieved when the electronics of a sulfonyl chloride activator and the reactivity of a glycosyl donor hemiacetal are matched. While these reactions are compatible with the acid- and base-sensitive protecting groups that are commonly used in oligosaccharide synthesis, these protecting groups are not relied upon to control selectivity. Instead, β-selectivity arises from the stereoinversion of an α-glycosyl arylsulfonate in an S_N2-like mechanism. Our mechanistic proposal is supported by NMR studies, kinetic isotope effect (KIE) measurements, and DFT calculations.

Computational simulation of mechanism and isotope effects on acetal heterolysis as a model for glycoside hydrolysis

DFT calculations for the equilibrium isotope effect for deuterium substitution at the anomeric centre Cα in 2-(p-nitrophenoxy)tetrahydropyran with continuum solvation show significant variation in the range of relative permittivity 2 ≤ ε ≤ 10. One-dimensional scans of potential energy (with implicit solvation by water) or of free energy (from QM/MM potentials of mean force with explicit aqueous solvation with a hybrid AM1/OPLS method) for heterolysis of the bond between Cα and the nucleofuge do not show a transition state. A two-dimensional free-energy surface that considers also the distance between Cα and a nucleophilic water indicates a pre-association DN*ANint ‡ mechanism with a transition state involving nucleophilic attack upon an ion-pair intermediate, and this is supported by good agreement between the mean values of the calculated and experimental α-D KIEs. However, the magnitudes of the standard deviations about the mean values for the making and breaking C–O bonds suggest that the transition state is rather plastic, with Cα–Onu≈2 ± 0.4 Å and Cα–Olg≈3 ± 0.5 Å. Not only is nucleophilic solvent assistance necessary, but there is also evidence for electrophilic assistance through specific hydrogen bonding to the nucleofuge.

C-H Bond Cleavage Is Rate-Limiting for Oxidative C-P Bond Cleavage by the Mixed Valence Diiron-Dependent Oxygenase PhnZ

PhnZ utilizes a mixed valence diiron(II/III) cofactor and O₂ to oxidatively cleave the carbon-phosphorus bond of (R)-2-amino-1-hydroxyethylphosphonic acid to form glycine and orthophosphate. The active site residues Y24 and E27 are proposed to mediate induced-fit recognition of the substrate and access of O₂ to one of the active site Fe ions. H62 is proposed to deprotonate the C1-hydroxyl of the substrate during catalysis. Kinetic isotope effects (KIEs), pH-rate dependence, and site-directed mutagenesis were used to probe the rate-determining transition state and the roles of these three active site residues. Primary deuterium KIE values of 5.5 ± 0.3 for ^D(V) and 2.2 ± 0.4 for ^D(V/K) were measured with (R)-2-amino[1-²H₁]-1-hydroxyethylphosphonic acid, indicating that cleavage of the C1-H bond of the substrate is rate-limiting. This step is also rate-limiting for PhnZ Y24F, as shown by a significant deuterium KIE value of 2.3 ± 0.1 for ^D(V). In contrast, a different reaction step appears to be rate-limiting for the PhnZ E27A and H62A variants, which exhibited ^D(V) values near unity. A solvent KIE of 2.2 ± 0.3 for ^D₂O(V) is observed for PhnZ. Significant solvent KIE values are also observed for the PhnZ Y24F and E27A variants. In contrast, the PhnZ H62A variant does not show a significant solvent KIE, suggesting that H62 is mediating proton transfer in the transition state. A proton inventory study with PhnZ indicates that 1.5 ± 0.6 protons are in flight in the rate-determining step. Overall, the rate-determining transition state for oxidative C-P bond cleavage by PhnZ is proposed to involve C-H bond cleavage that is coupled to deprotonation of the substrate C1-hydroxyl by H62.

Enzymatic Transition States and Drug Design

Transition state theory teaches that chemically stable mimics of enzymatic transition states will bind tightly to their cognate enzymes. Kinetic isotope effects combined with computational quantum chemistry provides enzymatic transition state information with sufficient fidelity to design transition state analogues. Examples are selected from various stages of drug development to demonstrate the application of transition state theory, inhibitor design, physicochemical characterization of transition state analogues, and their progress in drug development.

The Experimental Evidence in Support of Glycosylation Mechanisms at the SN1-SN2 Interface

A critical review of the state-of-the-art evidence in support of the mechanisms of glycosylation reactions is provided. Factors affecting the stability of putative oxocarbenium ions as intermediates at the SN1 end of the mechanistic continuum are first surveyed before the evidence, spectroscopic and indirect, for the existence of such species on the time scale of glycosylation reactions is presented. Current models for diastereoselectivity in nucleophilic attack on oxocarbenium ions are then described. Evidence in support of the intermediacy of activated covalent glycosyl donors is reviewed, before the influences of the structure of the nucleophile, of the solvent, of temperature, and of donor-acceptor hydrogen bonding on the mechanism of glycosylation reactions are surveyed. Studies on the kinetics of glycosylation reactions and the use of kinetic isotope effects for the determination of transition-state structure are presented, before computational models are finally surveyed. The review concludes with a critical appraisal of the state of the art.

Time-dependent view of an isotope effect in electron-nuclear nonequilibrium dynamics with applications to N2

Isotopic fractionation in the photodissociation of N₂ could explain the considerable variation in the ¹⁴N/¹⁵N ratio in different regions of our galaxy. We previously proposed that such an isotope effect is due to coupling of photoexcited bound valence and Rydberg electronic states in the frequency range where there is strong state mixing. We here identify features of the role of the mass in the dynamics through a time-dependent quantum-mechanical simulation. The photoexcitation of N₂ is by an ultrashort pulse so that the process has a sharply defined origin in time and so that we can monitor the isolated molecule dynamics in time. An ultrafast pulse is necessarily broad in frequency and spans several excited electronic states. Each excited molecule is therefore not in a given electronic state but in a superposition state. A short time after excitation, there is a fairly sharp onset of a mass-dependent large population transfer when wave packets on two different electronic states in the same molecule overlap. This coherent overlap of the wave packets on different electronic states in the region of strong coupling allows an effective transfer of population that is very mass dependent. The extent of the transfer depends on the product of the populations on the two different electronic states and on their relative phase. It is as if two molecules collide but the process occurs within one molecule, a molecule that is simultaneously in both states. An analytical toy model recovers the (strong) mass and energy dependence.

Transition State Structure and Inhibition of Rv0091, a 5'-Deoxyadenosine/5'-methylthioadenosine Nucleosidase from Mycobacterium tuberculosis

5'-Methylthioadenosine/S-adenosylhomocysteine nucleosidase (MTAN) is a bacterial enzyme that catalyzes the hydrolysis of the N-ribosidic bond in 5'-methylthioadenosine (MTA) and S-adenosylhomocysteine (SAH). MTAN activity has been linked to quorum sensing pathways, polyamine biosynthesis, and adenine salvage. Previously, the coding sequence of Rv0091 was annotated as a putative MTAN in Mycobacterium tuberculosis. Rv0091 was expressed in Escherichia coli, purified to homogeneity, and shown to be a homodimer, consistent with MTANs from other microorganisms. Substrate specificity for Rv0091 gave a preference for 5'-deoxyadenosine relative to MTA or SAH. Intrinsic kinetic isotope effects (KIEs) for the hydrolysis of [1'-(3)H], [1'-(14)C], [5'-(3)H2], [9-(15)N], and [7-(15)N]MTA were determined to be 1.207, 1.038, 0.998, 1.021, and 0.998, respectively. A model for the transition state structure of Rv0091 was determined by matching KIE values predicted via quantum chemical calculations to the intrinsic KIEs. The transition state shows a substantial loss of C1'-N9 bond order, well-developed oxocarbenium character of the ribosyl ring, and weak participation of the water nucleophile. Electrostatic potential surface maps for the Rv0091 transition state structure show similarity to DADMe-immucillin transition state analogues. DADMe-immucillin transition state analogues showed strong inhibition of Rv0091, with the most potent inhibitor (5'-hexylthio-DADMe-immucillinA) displaying a Ki value of 87 pM.

Transition State Structure of RNA Depurination by Saporin L3

Saporin L3 from the leaves of the common soapwort is a catalyst for hydrolytic depurination of adenine from RNA. Saporin L3 is a type 1 ribosome inactivating protein (RIP) composed only of a catalytic domain. Other RIPs have been used in immunotoxin cancer therapy, but off-target effects have limited their development. In the current study, we use transition state theory to understand the chemical mechanism and transition state structure of saporin L3. In favorable cases, transition state structures guide the design of transition state analogues as inhibitors. Kinetic isotope effects (KIEs) were determined for an A14C mutant of saporin L3. To permit KIE measurements, small stem-loop RNAs that contain an AGGG tetraloop structure were enzymatically synthesized with the single adenylate bearing specific isotopic substitutions. KIEs were measured and corrected for forward commitment to obtain intrinsic values. A model of the transition state structure for depurination of stem-loop RNA (5'-GGGAGGGCCC-3') by saporin L3 was determined by matching KIE values predicted via quantum chemical calculations to a family of intrinsic KIEs. This model indicates saporin L3 displays a late transition state with the N-ribosidic bond to the adenine nearly cleaved, and the attacking water nucleophile weakly bonded to the ribosyl anomeric carbon. The transition state retains partial ribocation character, a feature common to most N-ribosyl transferases. However, the transition state geometry for saporin L3 is distinct from ricin A-chain, the only other RIP whose transition state is known.

Transition state analysis of thymidine hydrolysis by human thymidine phosphorylase

Human thymidine phosphorylase (hTP) is responsible for thymidine (dT) homeostasis, and its action promotes angiogenesis. In the absence of phosphate, hTP catalyzes a slow hydrolytic depyrimidination of dT yielding thymine and 2-deoxyribose (dRib). Its transition state was characterized using multiple kinetic isotope effect (KIE) measurements. Isotopically enriched thymidines were synthesized enzymatically from glucose or (deoxy)ribose, and intrinsic KIEs were used to interpret the transition state structure. KIEs from [1'-(14)C]-, [1-(15)N]-, [1'-(3)H]-, [2'R-(3)H]-, [2'S-(3)H]-, [4'-(3)H]-, and [5'-(3)H]dTs provided values of 1.033 ± 0.002, 1.004 ± 0.002, 1.325 ± 0.003, 1.101 ± 0.004, 1.087 ± 0.005, 1.040 ± 0.003, and 1.033 ± 0.003, respectively. Transition state analysis revealed a stepwise mechanism with a 2-deoxyribocation formed early and a higher energetic barrier for nucleophilic attack of a water molecule on the high energy intermediate. An equilibrium exists between the deoxyribocation and reactants prior to the irreversible nucleophilic attack by water. The results establish activation of the thymine leaving group without requirement for phosphate. A transition state constrained to match the intrinsic KIEs was found using density functional theory. An active site histidine (His116) is implicated as the catalytic base for activation of the water nucleophile at the rate-limiting transition state. The distance between the water nucleophile and the anomeric carbon (r(C-O)) is predicted to be 2.3 A at the transition state. The transition state model predicts that deoxyribose adopts a mild 3'-endo conformation during nucleophilic capture. These results differ from the concerted bimolecular mechanism reported for the arsenolytic reaction (Birck, M. R.; Schramm, V. L. J. Am. Chem. Soc. 2004, 126, 2447-2453).

Mechanistic insights into glycosidase chemistry

The enzymatic hydrolysis of the glycosidic bond continues to gain importance, reflecting the critically important roles complex glycans play in health and disease as well as the rekindled interest in enzymatic biomass conversion. Recent advances include the broadening of our understanding of enzyme reaction coordinates, through both computational and structural studies, improved understanding of enzyme inhibition through transition state mimicry and fascinating insights into mechanism yielded by physical organic chemistry approaches.

Glycosyltransferases: structures, functions, and mechanisms

Glycosyltransferases catalyze glycosidic bond formation using sugar donors containing a nucleoside phosphate or a lipid phosphate leaving group. Only two structural folds, GT-A and GT-B, have been identified for the nucleotide sugar-dependent enzymes, but other folds are now appearing for the soluble domains of lipid phosphosugar-dependent glycosyl transferases. Structural and kinetic studies have provided new insights. Inverting glycosyltransferases utilize a direct displacement S(N)2-like mechanism involving an enzymatic base catalyst. Leaving group departure in GT-A fold enzymes is typically facilitated via a coordinated divalent cation, whereas GT-B fold enzymes instead use positively charged side chains and/or hydroxyls and helix dipoles. The mechanism of retaining glycosyltransferases is less clear. The expected two-step double-displacement mechanism is rendered less likely by the lack of conserved architecture in the region where a catalytic nucleophile would be expected. A mechanism involving a short-lived oxocarbenium ion intermediate now seems the most likely, with the leaving phosphate serving as the base.