Intrinsic Kinetic Isotope Effects Research Articles

Theoretical Kinetic Isotope Effects in Establishing the Precise Biodegradation Mechanisms of Organic Pollutants.

Compound-specific isotope analysis (CSIA) for natural isotope ratios has been recognized as a promising tool to elucidate biodegradation pathways of organic pollutants by microbial enzymes by relating reported kinetic isotope effects (KIEs) to apparent KIEs (AKIEs) derived from bulk isotope fractionations (εbulk). However, for many environmental reactions, neither are the reference KIE ranges sufficiently narrow nor are the mechanisms elucidated to the point that rate-determining steps have been identified unequivocally. In this work, besides providing reference KIEs and rationalizing AKIEs, good relationships have been explained by DFT computations for diverse biodegradation pathways with known enzymatic models between the theoretical isotope fractionations (εbulk') from intrinsic KIEs on the rate-determining steps and the observed εbulk. (1) To confirm the mechanistic details of previously reported pathway-dependent CSIA, it includes isotope changes in MTBE biodegradation between hydroxylation by CYP450 and SN2 reaction by cobalamin-dependent methyltransferase, the regioselectivity of toluene biodegradation by CYP450, and the rate-determining step in toluene biodegradation by benzylsuccinate synthase. (2) To yield new fundamental insights into some unclear biodegradation pathways, it consists of the oxidative function of toluene dioxygenase in biodegradation of TCE, the epoxidation mode in biodegradation of TCE by toluene 4-monooxygenase, and the weighted average mechanism in biodegradation of cDCE by CYP450.

Employing deuterium kinetic isotope effects to uncover the mechanism of (R)-3-hydroxybutyrate dehydrogenase.

Short-chain dehydrogenases/reductases (SDR) form a large enzyme superfamily playing important roles in health and disease. Furthermore, they are useful tools in biocatalysis. Unveiling the nature of the transition state for hydride transfer is a crucial undertaking toward defining the physicochemical underpinnings of catalysis by SDR enzymes, including possible contributions from quantum mechanical tunneling. Primary deuterium kinetic isotope effects can uncover the contribution from chemistry to the rate-limiting step and potentially provide detailed information on the hydride-transfer transition state in SDR-catalyzed reactions. For the latter, however, one needs to determine the intrinsic isotope effect: that which would be measured if hydride transfer were rate determining. Alas, as is the case for many other enzymatic reactions, those catalyzed by SDRs are often limited by the rate of isotope-insensitive steps, such as product release and conformational changes, which masks the expression of the intrinsic isotope effect. This can be overcome by the powerful yet underexplored method of Palfey and Fagan via which intrinsic kinetic isotope effects can be extracted from pre-steady-state kinetics data. SDRs are ideal systems to which this method can be applied. We have employed this approach to elucidate the transition states for hydride transfer catalyzed by NADH-dependent cold- and warm-adapted (R)-3-hydroxybutyrate dehydrogenase. Experimental conditions which simplify the analysis are discussed.

Microbial Denitrification: Active Site and Reaction Path Models Predict New Isotopic Fingerprints.

The study of isotopic fingerprints in nitrate (δ15N, δ18O, Δ17O) has enabled pivotal insights into the global nitrogen cycle and revealed new knowledge gaps. Measuring populations of isotopic homologs of intact NO3 - ions (isotopologues) shows promise to advance the understanding of nitrogen cycling processes; however, we need new theory and predictions to guide laboratory experiments and field studies. We investigated the hypothesis that the isotopic composition of the residual nitrate pool is controlled by the N-O bond-breaking step in Nar dissimilatory nitrate reductase using molecular models of the enzyme active sites and associated kinetic isotope effects (KIEs). We integrated the molecular model results into reaction path models representing the reduction of nitrate under either closed-system or steady-state conditions. The predicted intrinsic KIE (15ε and 18ε) of the Nar active site matches observed fractionations in both culture and environmental studies. This is what would be expected if the isotopic composition of marine nitrate were controlled by dissimilatory nitrate reduction by Nar. For a closed system, the molecular models predict a pronounced negative 15N-18O clumping anomaly in residual nitrate. This signal could encode information about the amount of nitrate consumed in a closed system and thus constrain initial nitrate concentration and its isotopic composition. Similar clumped isotope anomalies can potentially be used to distinguish whether a system is open or closed to new nitrate addition. These mechanistic predictions can be tested and refined in combination with emerging ESI-Orbitrap measurements.

Using kinetic isotope effects to probe the mechanism of adenosylcobalamin-dependent enzymes.

Adenosylcobalamin- (AdoCbl) dependent enzyme reactions involved the transfer of hydrogen atoms between the 5'-carbon of the coenzyme and the substrates and products of the reaction. Tritium and deuterium kinetic isotope effect measurements are, therefore, a valuable tool to probe the mechanisms of AdoCbl-dependent enzymes, as they can provide information about the reaction pathway and the rate-determining step. Furthermore, if the intrinsic kinetic isotope effect can be isolated, information on the nature of the transition state associated with hydrogen transfer can be obtained. In this chapter I present methods for the preparation of isotopically-labeled AdoCbl and their use in rapid chemical quench experiments that allow isotope effects on specific steps in the reaction to be isolated. These techniques are illustrated with examples from my laboratory's studies on the AdoCbl dependent enzyme, glutamate mutase.

Evolution of Optimized Hydride Transfer Reaction and Overall Enzyme Turnover in Human Dihydrofolate Reductase.

Evolution of dihydrofolate reductase (DHFR) has been studied using the enzyme from Escherichia coli DHFR (ecDHFR) as a model, but less studies have used the enzyme from Homo sapiens DHFR (hsDHFR). Each enzyme maintains a short and narrow distribution of hydride donor-acceptor distances (DAD) at the tunneling ready state (TRS). Evolution of the enzyme was previously studied in ecDHFR where three key sites were identified as important to the catalyzed reaction. The corresponding sites in hsDHFR are F28, 62-PEKN, and 26-PPLR. Each of these sites was studied here through the creation of mutant variants of the enzyme and measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs) on the reaction. F28 is mutated first to M (F28M) and then to the L of the bacterial enzyme (F28L). The KIEs of the F28M variant are larger and more temperature-dependent than wild-type (WT), suggesting a broader and longer average DAD at the TRS. To more fully mimic ecDHFR, we also study a triple mutant of the human enzyme (F32L-PP26N-PEKN62G). Remarkably, the intrinsic KIEs, while larger in magnitude, are temperature-independent like the WT enzymes. We also construct deletion mutations of hsDHFR removing both the 62-PEKN and 26-PPLR sequences. The results mirror those described previously for insertion mutants of ecDHFR. Taken together, these results suggest a balancing act during DHFR evolution between achieving an optimal TRS for hydride transfer and preventing product inhibition arising from the different intercellular pools of NADPH and NADP+ in prokaryotic and eukaryotic cells.

Kinetic Studies of the Hydrogen Atom Transfer in a Hypoxia-Sensing Enzyme, FIH-1: KIE and O2 Reactivity.

Cellular hypoxia plays a crucial role in tissue development and adaptation to pO2. Central to cellular oxygen sensing is factor-inhibiting HIF-1α (FIH), an α-ketoglutarate (αKG)/non-heme iron(II)-dependent dioxygenase that hydroxylates a specific asparagine residue of hypoxia inducible factor-1α (HIF-1α). The high KM(O2) and rate-limiting decarboxylation step upon O2 activation are key features of the enzyme that classify it as an oxygen sensor and set it apart from other αKG/Fe(II)-dependent dioxygenases. Although the chemical intermediates following decarboxylation are presumed to follow the consensus mechanism of other αKG/Fe(II)-dependent dioxygenases, experiments have not previously demonstrated these canonical steps in FIH. In this work, a deuterated peptide substrate was used as a mechanistic probe for the canonical hydrogen atom transfer (HAT). Our data show a large kinetic isotope effect (KIE) in steady-state kinetics (Dkcat = 10 ± 1), revealing that the HAT occurs and is partially rate limiting on kcat. Kinetic studies showed that the deuterated peptide led FIH to uncouple O2 activation and provided the opportunity to spectroscopically observe the ferryl intermediate. This enzyme uncoupling was used as an internal competition with respect to the fate of the ferryl intermediate, demonstrating a large observed KIE on the uncoupling (Dk5 = 1.147 ± 0.005) and an intrinsic KIE on the HAT step (Dk > 15). The close energy barrier between αKG decarboxylation and HAT distinguishes FIH as an O2-sensing enzyme and is crucial for ensuring substrate specificity in the regulation of cellular O2 homeostasis.

Transition States for Psychrophilic and Mesophilic (R)-3-Hydroxybutyrate Dehydrogenase-Catalyzed Hydride Transfer at Sub-zero Temperatures.

(R)-3-Hydroxybutyrate dehydrogenase (HBDH) catalyzes the NADH-dependent reduction of 3-oxocarboxylates to (R)-3-hydroxycarboxylates. The active sites of a pair of cold- and warm-adapted HBDHs are identical except for a single residue, yet kinetics evaluated at -5, 0, and 5 °C show a much higher steady-state rate constant (kcat) for the cold-adapted than for the warm-adapted HBDH. Intriguingly, single-turnover rate constants (kSTO) are strikingly similar between the two orthologues. Psychrophilic HBDH primary deuterium kinetic isotope effects on kcat (Dkcat) and kSTO (DkSTO) decrease at lower temperatures, suggesting more efficient hydride transfer relative to other steps as the temperature decreases. However, mesophilic HBDH Dkcat and DkSTO are generally temperature-independent. The DkSTO data allowed calculation of intrinsic primary deuterium kinetic isotope effects. Intrinsic isotope effects of 4.2 and 3.9 for cold- and warm-adapted HBDH, respectively, at 5 °C, supported by quantum mechanics/molecular mechanics calculations, point to a late transition state for both orthologues. Conversely, intrinsic isotope effects of 5.7 and 3.1 for cold- and warm-adapted HBDH, respectively, at -5 °C indicate the transition state becomes nearly symmetric for the psychrophilic enzyme, but more asymmetric for the mesophilic enzyme. His-to-Asn and Asn-to-His mutations in the psychrophilic and mesophilic HBDH active sites, respectively, swap the single active-site position where these orthologues diverge. At 5 °C, the His-to-Asn mutation in psychrophilic HBDH decreases Dkcat to 3.1, suggesting a decrease in transition-state symmetry, while the His-to-Asn mutation in mesophilic HBDH increases Dkcat to 4.4, indicating an increase in transition-state symmetry. Hence, temperature adaptation and a single divergent active-site residue may influence transition-state geometry in HBDHs.

Temperature-Dependent Kinetic Isotope Effects in R67 Dihydrofolate Reductase from Path-Integral Simulations.

Calculation of temperature-dependent kinetic isotope effects (KIE) in enzymes presents a significant theoretical challenge. Additionally, it is not trivial to identify enzymes with available experimental accurate intrinsic KIEs in a range of temperatures. In the current work, we present a theoretical study of KIEs in the primitive R67 dihydrofolate reductase (DHFR) enzyme and compare with experimental work. The advantage of R67 DHFR is its significantly lower kinetic complexity compared to more evolved DHFR isoforms. We employ mass-perturbation-based path-integral simulations in conjunction with umbrella sampling and a hybrid quantum mechanics–molecular mechanics Hamiltonian. We obtain temperature-dependent KIEs in good agreement with experiments and ascribe the temperature-dependent KIEs primarily to zero-point energy effects. The active site in the primitive enzyme is found to be poorly preorganized, which allows excessive water access to the active site and results in loosely bound reacting ligands.

Isotopic Labeling of Formate Dehydrogenase Perturbs the Protein Dynamics.

Isotope substitution of enzymes has become a means of addressing the participation of protein motions in enzyme-catalyzed reactions. The idea is that only the enzyme mass will be altered and not the electrostatics, so that the protein dynamics are essentially the same but at lower frequencies because of the mass change. In this study, we variably label all carbon atoms in formate dehydrogenase (FDH) with 13C, all nitrogen atoms with 15N, and all nonexchangeable hydrogen atoms with deuterium and investigate the impact that isotopic substitution has on the dynamics at the active site by two-dimensional infrared spectroscopy and compare with the measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs). We show that 15N labeling of FDH has the largest effect and makes the active site more heterogeneous, whereas the addition of nonexchangeable deuterium appears to have the opposite effect of 15N on active-site dynamics, resulting in a behavior similar to that of native FDH. Nevertheless, the temperature dependence of the KIEs shows a monotonic trend with protein mass that does not correspond with the changes in dynamics. These results suggest that isotope labeling has more than just a mass effect on enzyme dynamics and may influence electrostatics in ways that complicate the interpretation of the protein isotope effect.

Evolution Conserves the Network of Coupled Residues in Dihydrofolate Reductase.

Understanding protein motions and their role in enzymatic reactions is an important and timely topic in enzymology. Protein motions that are involved in the chemical step of catalysis are particularly intriguing but difficult to identify. A global network of coupled residues in Escherichia coli dihydrofolate reductase (E. coli DHFR), which assists in catalyzing the chemical step, has previously been demonstrated through quantum mechanical/molecular mechanical and molecular dynamics simulations as well as bioinformatic analyses. A few specific residues (M42, G121, F125, 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. This study hypothesizes that the global network of residues involved in the chemical step is evolutionarily conserved and probes homologous residues of the potential global network in human DHFR through measurements of the temperature dependence of KIEsint and computer simulations based on the empirical valence bond method. We study mutants M53W and S145V. Both of these remote residues are homologous to network residues in E. coli DHFR. Non-additive isotope effects on activation energy are observed between M53 and S145, indicating their synergistic effect on the chemical step in human DHFR, which suggests that both of these residues are part of a network affecting the chemical step in enzyme catalysis. This finding supports the hypothesis that human and E. coli DHFR share similar networks, consistent with evolutionary preservation of such networks.

Understanding Biological Hydrogen Transfer Through the Lens of Temperature Dependent Kinetic Isotope Effects.

Hydrogen atom transfer (HAT) is a salient feature of many enzymatic C-H cleavage mechanisms. In systems where kinetic isolation of HAT is achieved, selective labeling of substrate with hydrogen isotopes, such as deuterium, enables the determination of intrinsic kinetic isotope effects (KIEs). While the magnitude of the KIE is itself informative, ultimately the size of the temperature dependence of the KIE, Δ Ea = Ea(D) - Ea(H), serves as a critical, and often misinterpreted (or even ignored) descriptor of the reaction coordinate. As will be highlighted in this Account, Δ Ea is one of the most robust parameters to emerge from studies of enzyme catalyzed hydrogen transfer. Kinetic parameters for C-H reactions via HAT can appear consistent with either classical "over-the-barrier" or "Bell-like tunneling correction" models. However, neither of these models is able to explain the observation of near-zero Δ Ea values with many native enzymes that increase upon extrinsic or intrinsic perturbations to function. Instead, a full tunneling model has been developed that can account for the aggregate trends in the temperature dependence of the KIE. This model is reminiscent of Marcus-like theory for electron tunneling, with the additional incorporation of an H atom donor-acceptor distance (DAD) sampling term for effective wave function overlap; the role of the latter term is manifested in the experimentally determined Δ Ea. Three enzyme systems from this laboratory that illustrate different aspects of HAT are presented: taurine dioxygenase, the dual copper β-monooxygenases, and soybean lipoxygenase (SLO). The latter provides a particularly compelling system for understanding the properties of hydrogen tunneling, showing systematic increases in Δ Ea upon reduction in the size of hydrophobic residues both proximal and distal from the active site iron cofactor. Of note, recent ENDOR-based studies of enzyme-substrate complexes with SLO indicate an increase in DAD for mutants with increased Δ Ea, observations that are inconsistent with "Bell-like correction" models. Overall, the surmounting kinetic and biophysical evidence corroborates a multidimensional approach for understanding HAT, offering a robust mechanistic explanation for the magnitude and trends of the KIE and Δ Ea. Recent DFT and QM/MM computations on SLO are compared to the developed nonadiabatic analytical constructs, providing considerable insight into ground state structures and reactivity. However, QM/MM is unable to readily reproduce the small Δ Ea values characteristic of native enzymes. Future theoretical developments to capture these experimental observations may necessitate a parsing of protein motions for local, substrate deuteration-sensitive modes from isotope-insensitive modes within the larger conformational landscape, in the process providing deeper understanding of how native enzymes have evolved to transiently optimize their active site configurations.

Bacterial versus human thymidylate synthase: Kinetics and functionality.

Thymidylate Synthase (TSase) is a highly conserved enzyme that catalyzes the production of the DNA building block thymidylate. Structurally, functionally and mechanistically, bacterial and mammalian TSases share remarkable similarities. Because of this closeness, bacterial enzymes have long been used as model systems for human TSase. Furthermore, while TSase inhibitors have long served as chemotherapeutic drugs, no TSase inhibitor serves as an antibiotic. Despite their high resemblance, the mammalian TSases are distinct in a few known aspects, such as having a N-terminal tail and two insertions in the primary sequence and active/inactive conformations. Here, we aim to comprehensively characterize human (hs) TSase and delineate its contrasts and the similarities to the well-studied Escherichia coli (ec) TSase. We found that, in contrast to ecTSase, Mg2+ does not enhance reaction rates for hsTSase. The temperature dependence of intrinsic kinetic isotope effects (KIEs), on the other hand, suggests that Mg2+ has little or no impact on the transition state of hydride transfer in either enzyme, and that the transition state for the hydride transfer in hsTSase is looser than in ecTSase. Additionally, the substrates’ binding order is strictly ordered for ecTSase but slightly less ordered for hsTSase. The observed kinetic and functional differences between bacterial and human enzymes may aid in the development of antibiotic drugs with reduced toxicity.

Protein Mass Effects on Formate Dehydrogenase.

Isotopically labeled enzymes (denoted as "heavy" or "Born-Oppenheimer" enzymes) have been used to test the role of protein dynamics in catalysis. The original idea was that the protein's higher mass would reduce the frequency of its normal-modes without altering its electrostatics. Heavy enzymes have been used to test if the vibrations in the native enzyme are coupled to the chemistry it catalyzes, and different studies have resulted in ambiguous findings. Here the temperature-dependence of intrinsic kinetic isotope effects of the enzyme formate dehydrogenase is used to examine the distribution of H-donor to H-acceptor distance as a function of the protein's mass. The protein dynamics are altered in the heavy enzyme to diminish motions that determine the transition state sampling in the native enzyme, in accordance with a Born-Oppenheimer-like effect on bond activation. Findings of this work suggest components related to fast frequencies that can be explained by Born-Oppenheimer enzyme hypothesis (vibrational) and also slower time scale events that are non-Born-Oppenheimer in nature (electrostatic), based on evaluations of protein mass dependence of donor-acceptor distance and forward commitment to catalysis along with steady state and single turnover measurements. Together, the findings suggest that the mass modulation affected both local, fast, protein vibrations associated with the catalyzed chemistry and the protein's macromolecular electrostatics at slower time scales; that is, both Born-Oppenheimer and non-Born-Oppenheimer effects are observed. Comparison to previous studies leads to the conclusion that isotopic labeling of the protein may have different effects on different systems, however, making heavy enzyme studies a very exciting technique for exploring the dynamics link to catalysis in proteins.

Kinetic Isotope Effects and Transition State Structure for Hypoxanthine-Guanine-Xanthine Phosphoribosyltransferase from Plasmodium falciparum.

Plasmodium falciparum parasites are purine auxotrophs that rely exclusively on the salvage of preformed purines from their human hosts to supply the requirement for purine nucleotides. Hypoxanthine-guanine-xanthine phosphoribosyltransferase (HGXPRT) catalyzes the freely reversible Mg2+-dependent conversion of 6-oxopurine bases to their respective nucleotides and inorganic pyrophosphate. The phosphoribosyl group is derived from 5-phospho-α-d-ribosyl 1-pyrophosphate (PRPP). The enzyme from malaria parasites (PfHGXPRT) is essential as hypoxanthine is the major precursor in purine metabolism. We used specific heavy atom labels in PRPP and hypoxanthine to measure primary (1-14C and 9-15N) and secondary (1-3H and 7-15N) intrinsic kinetic isotope effect (KIE) values for PfHGXPRT. Intrinsic isotope effects contain information for understanding enzymatic transition state properties. The transition state of PfHGXPRT was explored by matching KIE values predicted from quantum mechanical calculations to the intrinsic values determined experimentally. This approach provides information about PfHGXPRT transition state bond lengths, geometry, and atomic charge distribution. The transition state structure of PfHGXPRT was determined in the physiological direction of addition of ribose 5-phosphate to hypoxanthine by overcoming the chemical instability of PRPP. The transition state for PfHGXPRT forms nucleotides through a well-developed and near-symmetrical DN*AN, SN1-like transition state.

The Transition-State Structure for Human MAT2A from Isotope Effects

Human methionine S-adenosyltransferase (MAT2A) catalyzes the formation of S-adenosylmethionine (SAM) from ATP and methionine. Synthetic lethal genetic analysis has identified MAT2A as an anticancer target in tumor cells lacking expression of 5'-methylthioadenosine phosphorylase (MTAP). Approximately 15% of human cancers are MTAP-/-. The remainder can be rendered MTAP- through MTAP inhibitors. We used kinetic isotope effect (KIE), commitment factor (Cf), and binding isotope effect (BIE) measurements combined with quantum mechanical (QM) calculations to solve the transition state structure of human MAT2A. The reaction is characterized by an advanced SN2 transition state. The bond forming from the nucleophilic methionine sulfur to the 5'-C of ATP is 2.03 Å at the transition state (bond order of 0.67). Departure of the leaving group triphosphate of ATP is well advanced and forms a 2.32 Å bond between the 5'-C of ATP and the oxygen of the triphosphate (bond order of 0.23). Interaction of MAT2A with its MAT2B regulatory subunit causes no change in the intrinsic KIEs, indicating the same transition state structure. The transition state for MAT2A is more advanced along the reaction coordinate (more product-like) than that from the near-symmetrical transition state of methionine adenosyltransferase from E.coli.

Transition State Analysis of Adenosine Triphosphate Phosphoribosyltransferase.

Adenosine triphosphate phosphoribosyltransferase (ATP-PRT) catalyzes the first step in histidine biosynthesis, a pathway essential to microorganisms and a validated target for antimicrobial drug design. The ATP-PRT enzyme catalyzes the reversible substitution reaction between phosphoribosyl pyrophosphate and ATP. The enzyme exists in two structurally distinct forms, a short- and a long-form enzyme. These forms share a catalytic core dimer but bear completely different allosteric domains and thus distinct quaternary assemblies. Understanding enzymatic transition states can provide essential information on the reaction mechanisms and insight into how differences in domain structure influence the reaction chemistry, as well as providing a template for inhibitor design. In this study, the transition state structures for ATP-PRT enzymes from Campylobacter jejuni and Mycobacterium tuberculosis (long-form enzymes) and from Lactococcus lactis (short-form) were determined and compared. Intrinsic kinetic isotope effects (KIEs) were obtained at reaction sensitive positions for the reverse reaction using phosphonoacetic acid, an alternative substrate to the natural substrate pyrophosphate. The experimental KIEs demonstrated mechanistic similarities between the three enzymes and provided experimental boundaries for quantum chemical calculations to characterize the transition states. Predicted transition state structures support a dissociative reaction mechanism with a DN*AN‡ transition state. Weak interactions from the incoming nucleophile and a fully dissociated ATP adenine are predicted regardless of the difference in overall structure and quaternary assembly. These studies establish that despite significant differences in the quaternary assembly and regulatory machinery between ATP-PRT enzymes from different sources, the reaction chemistry and catalytic mechanism are conserved.

Kinetic Isotope Effects and Transition State Structure for Human Phenylethanolamine N-Methyltransferase.

Phenylethanolamine N-methyltransferase (PNMT) catalyzes the S-adenosyl-l-methionine (SAM)-dependent conversion of norepinephrine to epinephrine. Epinephrine has been associated with critical processes in humans including the control of respiration and blood pressure. Additionally, PNMT activity has been suggested to play a role in hypertension and Alzheimer's disease. In the current study, labeled SAM substrates were used to measure primary methyl-14C and 36S and secondary methyl-3H, 5'-3H, and 5'-14C intrinsic kinetic isotope effects for human PNMT. The transition state of human PNMT was modeled by matching kinetic isotope effects predicted via quantum chemical calculations to intrinsic values. The model provides information on the geometry and electrostatics of the human PNMT transition state structure and indicates that human PNMT catalyzes the formation of epinephrine through an early SN2 transition state in which methyl transfer is rate-limiting.

Kinetic isotope effects reveal early transition state of protein lysine methyltransferase SET8

Protein lysine methyltransferases (PKMTs) catalyze the methylation of protein substrates, and their dysregulation has been linked to many diseases, including cancer. Accumulated evidence suggests that the reaction path of PKMT-catalyzed methylation consists of the formation of a cofactor(cosubstrate)-PKMT-substrate complex, lysine deprotonation through dynamic water channels, and a nucleophilic substitution (SN2) transition state for transmethylation. However, the molecular characters of the proposed process remain to be elucidated experimentally. Here we developed a matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS) method and corresponding mathematic matrix to determine precisely the ratios of isotopically methylated peptides. This approach may be generally applicable for examining the kinetic isotope effects (KIEs) of posttranslational modifying enzymes. Protein lysine methyltransferase SET8 is the sole PKMT to monomethylate histone 4 lysine 20 (H4K20) and its function has been implicated in normal cell cycle progression and cancer metastasis. We therefore implemented the MS-based method to measure KIEs and binding isotope effects (BIEs) of the cofactor S-adenosyl-l-methionine (SAM) for SET8-catalyzed H4K20 monomethylation. A primary intrinsic 13C KIE of 1.04, an inverse intrinsic α-secondary CD3 KIE of 0.90, and a small but statistically significant inverse CD3 BIE of 0.96, in combination with computational modeling, revealed that SET8-catalyzed methylation proceeds through an early, asymmetrical SN2 transition state with the C-N and C-S distances of 2.35-2.40 Å and 2.00-2.05 Å, respectively. This transition state is further supported by the KIEs, BIEs, and steady-state kinetics with the SAM analog Se-adenosyl-l-selenomethionine (SeAM) as a cofactor surrogate. The distinct transition states between protein methyltransferases present the opportunity to design selective transition-state analog inhibitors.

Chemical Mechanism of the Branched-Chain Aminotransferase IlvE from Mycobacterium tuberculosis.

The biosynthetic pathway of the branched-chain amino acids is essential for Mycobacterium tuberculosis growth and survival. We report here the kinetic and chemical mechanism of the pyridoxal 5'-phosphate (PLP)-dependent branched-chain aminotransferase, IlvE, from M. tuberculosis (MtIlvE). This enzyme is responsible for the final step of the synthesis of the branched-chain amino acids isoleucine, leucine, and valine. As seen in other aminotransferases, MtIlvE displays a ping-pong kinetic mechanism. pK values were identified from the pH dependence on V as well as V/K, indicating that the phosphate ester of the PLP cofactor, and the α-amino group from l-glutamate and the active site Lys204, play roles in acid-base catalysis and binding, respectively. An intrinsic primary kinetic isotope effect was identified for the α-C-H bond cleavage of l-glutamate. Large solvent kinetic isotope effect values for the ping and pong half-reactions were also identified. The absence of a quininoid intermediate in combination with the Dkobs in our multiple kinetic isotope effects under single-turnover conditions suggests a concerted type of mechanism. The deprotonation of C2 of l-glutamate and the protonation of C4' of the PLP cofactor happen synchronously in the ping half-reaction. A chemical mechanism is proposed on the basis of the results obtained here.

Structure–Reactivity Effects on Intrinsic Primary Kinetic Isotope Effects for Hydride Transfer Catalyzed by Glycerol-3-phosphate Dehydrogenase

Primary deuterium kinetic isotope effects (1°DKIE) on (kcat/KGA, M–1 s–1) for dianion (X2–) activated hydride transfer from NADL to glycolaldehyde (GA) catalyzed by glycerol-3-phosphate dehydrogenase were determined over a 2100-fold range of enzyme reactivity: (X2–, 1°DKIE); FPO32–, 2.8 ± 0.1; HPO32–, 2.5 ± 0.1; SO42–, 2.8 ± 0.2; HOPO32–, 2.5 ± 0.1; S2O32–, 2.9 ± 0.1; unactivated; 2.4 ± 0.2. Similar 1°DKIEs were determined for kcat. The observed 1°DKIEs are essentially independent of changes in enzyme reactivity with changing dianion activator. The results are consistent with (i) fast and reversible ligand binding; (ii) the conclusion that the observed 1°DKIEs are equal to the intrinsic 1°DKIE on hydride transfer from NADL to GA; (iii) similar intrinsic 1°DKIEs on GPDH-catalyzed reduction of the substrate pieces and the whole physiological substrate dihydroxyacetone phosphate. The ground-state binding interactions for different X2– are similar, but there are large differences in the transition state interactions for different X2–. The changes in transition state binding interactions are expressed as changes in kcat and are proposed to represent changes in stabilization of the active closed form of GPDH. The 1°DKIEs are much smaller than observed for enzyme-catalyzed hydrogen transfer that occurs mainly by quantum-mechanical tunneling.