Enzymatic catalysis of proton transfer at carbon: activation of triosephosphate isomerase by phosphite dianion

Glycerol 3-Phosphate Dehydrogenase: Role of the Protein Conformational Change in Activation of a Readily Reversible Enzyme-Catalyzed Hydride Transfer Reaction

Kinetic parameters are reported for glycerol 3-phosphate dehydrogenase (GPDH)-catalyzed hydride transfer from the whole substrate glycerol 3-phosphate (G3P) or truncated substrate ethylene glycol (EtG) to NAD, and for activation of the hydride transfer reaction of EtG by phosphite dianion. These kinetic parameters were combined with parameters for enzyme-catalyzed hydride transfer in the microscopic reverse direction to give the reaction equilibrium constants Keq. Hydride transfer from G3P is favored in comparison to EtG because the carbonyl product of the former reaction is stabilized by hyperconjugative electron donation from the -CH2R keto substituent. The kinetic data show that the phosphite dianion provides the same 7.6 ± 0.1 kcal/mol stabilization of the transition states for enzyme-catalyzed reactions in the forward [reduction of NAD by EtG] and reverse [oxidation of NADH by glycolaldehyde] directions. The experimental evidence that supports a role for phosphite dianion in stabilizing the active closed form of the GPDH (EC) relative to the ca. 6 kcal/mol more unstable open form (EO) is summarized.

Triosephosphate Isomerase: The Crippling Effect of the P168A/I172A Substitution at the Heart of an Enzyme Active Site

The P168 and I172 side chains sit at the heart of the active site of triosephosphate isomerase (TIM) and play important roles in the catalysis of the isomerization reaction. The phosphodianion of substrate glyceraldehyde 3-phosphate (GAP) drives a conformational change at the TIM that creates a steric interaction with the P168 side chain that is relieved by the movement of P168 that carries the basic E167 side chain into a clamp that consists of the hydrophobic I172 and L232 side chains. The P168A/I172A substitution at TIM from Trypanosoma brucei brucei (TbbTIM) causes a large 120,000-fold decrease in kcat for isomerization of GAP that eliminates most of the difference in the reactivity of TIM compared to the small amine base quinuclidinone for deprotonation of catalyst-bound GAP. The I172A substitution causes a > 2-unit decrease in the pKa of the E167 carboxylic acid in a complex to the intermediate analog PGA, but the P168A substitution at the I172A variant has no further effect on this pKa. The P168A/I172A substitutions cause a 5-fold decrease in Km for the isomerization of GAP from a 0.9 kcal/mol stabilization of the substrate Michaelis complexes. The results show that the P168 and I172 side chains play a dual role in destabilizing the ground-state Michaelis complex to GAP and in promoting stabilization of the transition state for substrate isomerization. This is consistent with an important role for these side chains in an induced fit reaction mechanism [Richard, J. P. (2022) Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution. Biochemistry 61, 1533-1542].

Human 2'-Deoxynucleoside 5'-Phosphate N-Hydrolase 1: Mechanism of 2'-Deoxyuridine 5'-Monophosphate Hydrolysis

The enzyme 2'-deoxynucleoside 5'-phosphate N-hydrolase 1 (DNPH1) catalyzes the N-ribosidic bond cleavage of 5-hydroxymethyl-2'-deoxyuridine 5'-monophosphate to generate 2-deoxyribose 5-phosphate and 5-hydroxymethyluracil. DNPH1 accepts other 2'-deoxynucleoside 5'-monophosphates as slow-reacting substrates. DNPH1 inhibition is a promising strategy to overcome resistance to and potentiate anticancer poly(ADP-ribose) polymerase inhibitors. We solved the crystal structure of unliganded human DNPH1 and took advantage of the slow reactivity of 2'-deoxyuridine 5'-monophosphate (dUMP) as a substrate to obtain a crystal structure of the DNPH1:dUMP Michaelis complex. In both structures, the carboxylate group of the catalytic Glu residue, proposed to act as a nucleophile in covalent catalysis, forms an apparent low-barrier hydrogen bond with the hydroxyl group of a conserved Tyr residue. The crystal structures are supported by functional data, with liquid chromatography-mass spectrometry analysis showing that DNPH1 incubation with dUMP leads to slow yet complete hydrolysis of the substrate. A direct UV-vis absorbance-based assay allowed characterization of DNPH1 kinetics at low dUMP concentrations. A bell-shaped pH-rate profile indicated that acid-base catalysis is operational and that for maximum kcat/KM, two groups with an average pKa of 6.4 must be deprotonated, while two groups with an average pKa of 8.2 must be protonated. A modestly inverse solvent viscosity effect rules out diffusional processes involved in dUMP binding to and possibly uracil release from the enzyme as rate limiting to kcat/KM. Solvent deuterium isotope effects on kcat/KM and kcat were inverse and unity, respectively. A reaction mechanism for dUMP hydrolysis is proposed.

Utilization of Cofactor Binding Energy for Enzyme Catalysis: Formate Dehydrogenase-Catalyzed Reactions of the Whole NAD Cofactor and Cofactor Pieces

The pressure to optimize enzymatic rate accelerations has driven the evolution of the induced-fit mechanism for enzyme catalysts where the binding interactions of nonreacting phosphodianion or adenosyl substrate pieces drive enzyme conformational changes to form protein substrate cages that are activated for catalysis. We report the results of experiments to test the hypothesis that utilization of the binding energy of the adenosine 5'-diphosphate ribose (ADP-ribose) fragment of the NAD cofactor to drive a protein conformational change activates Candida boidinii formate dehydrogenase (CbFDH) for catalysis of hydride transfer from formate to NAD+. The ADP-ribose fragment provides a >14 kcal/mol stabilization of the transition state for CbFDH-catalyzed hydride transfer from formate to NAD+. This is larger than the ca. 6 kcal/mol stabilization of the ground-state Michaelis complex between CbFDH and NAD+ (KNAD = 0.032 mM). The ADP, AMP, and ribose 5'-phosphate fragments of NAD+ activate CbFDH for catalysis of hydride transfer from formate to nicotinamide riboside (NR). At a 1.0 M standard state, these activators stabilize the hydride transfer transition states by ≈5.5 (ADP), 5.5 (AMP), and 4.4 (ribose 5'-phosphate) kcal/mol. We propose that activation by these cofactor fragments is partly or entirely due to the ion-pair interaction between the guanidino side chain cation of R174 and the activator phosphate anion. This substitutes for the interaction between the α-adenosyl pyrophosphate anion of the whole NAD+ cofactor that holds CbFDH in the catalytically active closed conformation.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

The Role of Asn11 in Catalysis by Triosephosphate Isomerase

Four catalytic amino acids at triosephosphate isomerase (TIM) are highly conserved: N11, K13, H95, and E167. Asparagine 11 is the last of these to be characterized in mutagenesis studies. The ND2 side chain atom of N11 is hydrogen bonded to the O-1 hydroxyl of enzyme-bound dihydroxyacetone phosphate (DHAP), and it sits in an extended chain of hydrogen-bonded side chains that includes T75' from the second subunit. The N11A variants of wild-type TIM from Trypanosoma brucei brucei (TbbTIM) and Leishmania mexicana (LmTIM) undergo dissociation from the dimer to monomer under our assay conditions. Values of Kas = 8 × 103 and 1 × 106 M-1, respectively, were determined for the conversion of monomeric N11A TbbTIM and LmTIM into their homodimers. The N11A substitution at the variant of LmTIM previously stabilized by the E65Q substitution gives the N11A/E65Q variant that is stable to dissociation under our assay conditions. The X-ray crystal structure of N11A/E65Q LmTIM shows an active site that is essentially superimposable on that for wild-type TbbTIM, which also has a glutamine at position 65. A comparison of the kinetic parameters for E65Q LmTIM and N11A/E65Q LmTIM-catalyzed reactions of (R)-glyceraldehyde 3-phosphate (GAP) and (DHAP) shows that the N11A substitution results in a (13-14)-fold decrease in kcat/Km for substrate isomerization and a similar decrease in kcat for DHAP but only a 2-fold decrease in kcat for GAP.

Kinetic characterization of methylthio-d-ribose-1-phosphate isomerase

Methylthio-d-ribose-1-phosphate (MTR1P) isomerase (MtnA) catalyzes the reversible isomerization of the aldose MTR1P into the ketose methylthio-d-ribulose 1-phosphate. It serves as a member of the methionine salvage pathway that many organisms require for recycling methylthio-d-adenosine, a byproduct of S-adenosylmethionine metabolism, back to methionine. MtnA is of mechanistic interest because unlike most other aldose-ketose isomerases, its substrate exists as an anomeric phosphate ester and therefore cannot equilibrate with a ring-opened aldehyde that is otherwise required to promote isomerization. To investigate the mechanism of MtnA, it is necessary to establish reliable methods for determining the concentration of MTR1P and to measure enzyme activity in a continuous assay. This chapter describes several such protocols needed to perform steady-state kinetics measurements. It additionally outlines the preparation of [³²P]MTR1P, its use in radioactively labeling the enzyme, and the characterization of the resulting phosphoryl adduct.

Kinetics and mechanism for enzyme-catalyzed reactions of substrate pieces

The most important difference between enzyme and small molecule catalysts is that only enzymes utilize the large intrinsic binding energies of nonreacting portions of the substrate in stabilization of the transition state for the catalyzed reaction. A general protocol is described to determine the intrinsic phosphodianion binding energy for enzymatic catalysis of reactions of phosphate monoester substrates, and the intrinsic phosphite dianion binding energy in activation of enzymes for catalysis of phosphodianion truncated substrates, from the kinetic parameters for enzyme-catalyzed reactions of whole and truncated substrates. The enzyme-catalyzed reactions so-far documented that utilize dianion binding interactions for enzyme activation; and, their phosphodianion truncated substrates are summarized. A model for the utilization of dianion binding interactions for enzyme activation is described. The methods for the determination of the kinetic parameters for enzyme-catalyzed reactions of whole and truncated substrates, from initial velocity data, are described and illustrated by graphical plots of kinetic data. The results of studies on the effect of site-directed amino acid substitutions at orotidine 5'-monophosphate decarboxylase, triosephosphate isomerase, and glycerol-3-phosphate dehydrogenase provide strong support for the proposal that these enzymes utilize binding interactions with the substrate phosphodianion to hold the protein catalysts in reactive closed conformations.

Uniform binding and negative catalysis at the origin of enzymes

Enzymes are well known for their catalytic abilities, some even reaching "catalytic perfection" in the sense that the reaction they catalyze has reached the physical bound of the diffusion rate. However, our growing understanding of enzyme superfamilies has revealed that only some share a catalytic chemistry while others share a substrate-handle binding motif, for example, for a particular phosphate group. This suggests that some families emerged through a "substrate-handle-binding-first" mechanism ("binding-first" for brevity) instead of "chemistry-first" and we are, therefore, left to wonder what the role of non-catalytic binders might have been during enzyme evolution. In the last of their eight seminal, back-to-back articles from 1976, John Albery and Jeremy Knowles addressed the question of enzyme evolution by arguing that the simplest mode of enzyme evolution is what they defined as "uniform binding" (parallel stabilization of all enzyme-bound states to the same degree). Indeed, we show that a uniform-binding proto-catalyst can accelerate a reaction, but only when catalysis is already present, that is, when the transition state is already stabilized to some degree. Thus, we sought an alternative explanation for the cases where substrate-handle-binding preceded any involvement of a catalyst. We find that evolutionary starting points that exhibit negative catalysis can redirect the reaction's course to a preferred product without need for rate acceleration or product release; that is, if they do not stabilize, or even destabilize, the transition state corresponding to an undesired product. Such a mechanism might explain the emergence of "binding-first" enzyme families like the aldolase superfamily.

Enabling Role of Ligand-Driven Conformational Changes in Enzyme Evolution

Many enzymes that show a large specificity in binding the enzymatic transition state with a higher affinity than the substrate utilize substrate binding energy to drive protein conformational changes to form caged substrate complexes. These protein cages provide strong stabilization of enzymatic transition states. Using part of the substrate binding energy to drive the protein conformational change avoids a similar strong stabilization of the Michaelis complex and irreversible ligand binding. A seminal step in the development of modern enzyme catalysts was the evolution of enzymes that couple substrate binding to a conformational change. These include enzymes that function in glycolysis (triosephosphate isomerase), the biosynthesis of lipids (glycerol phosphate dehydrogenase), the hexose monophosphate shunt (6-phosphogluconate dehydrogenase), and the mevalonate pathway (isopentenyl diphosphate isomerase), catalyze the final step in the biosynthesis of pyrimidine nucleotides (orotidine monophosphate decarboxylase), and regulate the cellular levels of adenine nucleotides (adenylate kinase). The evolution of enzymes that undergo ligand-driven conformational changes to form active protein-substrate cages is proposed to proceed by selection of variants, in which the selected side chain substitutions destabilize a second protein conformer that shows compensating enhanced binding interactions with the substrate. The advantages inherent to enzymes that incorporate a conformational change into the catalytic cycle provide a strong driving force for the evolution of flexible protein folds such as the TIM barrel. The appearance of these folds represented a watershed event in enzyme evolution that enabled the rapid propagation of enzyme activities within enzyme superfamilies.

Covalent Adduct Formation in Methylthio-d-ribose-1-phosphate Isomerase: Reaction Intermediate or Artifact?

Methylthio-d-ribose-1-phosphate (MTR1P) isomerase (MtnA) functions in the methionine salvage pathway by converting the cyclic aldose MTR1P to its open-chain ketose isomer methylthio-d-ribulose 1-phosphate (MTRu1P). What is particularly challenging for this enzyme is that the substrate's phosphate ester prevents facile equilibration to an aldehyde, which in other aldose-ketose isomerases is known to activate the α-hydrogen for proton or hydride transfer between adjacent carbons. We speculated that MtnA could use covalent catalysis via a phosphorylated residue to permit isomerization by one of the canonical mechanisms, followed by phosphoryl transfer back to form the product. In apparent support of this mechanism, [32P]MTR1P was found by SDS-PAGE and gel-filtration chromatography to radiolabel the enzyme. Susceptibility of this adduct to strongly acidic and basic pH and nucleophilic agents is consistent with an acyl phosphate. C160S and D240N, mutants of two conserved active-site residues, however, exhibited no difference in radiolabeling despite a reduction in activity of ∼107, leading to the conclusion that phosphorylation is unrelated to catalysis. Unexpectedly, prolonged incubations with C160S revealed up to 30% accumulation of radioactivity, which was identified by 31P and 13C NMR to be the result of a second adduct─a hemiketal formed between Ser160 and the carbonyl of MTRu1P. These results are interpreted as indirect support for a mechanism involving transfer of the proton from C-2 to C-1 by Cys160.

Glycerol-3-Phosphate Dehydrogenase: The K120 and K204 Side Chains Define an Oxyanion Hole at the Enzyme Active Site

The cationic K120 and K204 side chains lie close to the C-2 carbonyl group of substrate dihydroxyacetone phosphate (DHAP) at the active site of glycerol-3-phosphate dehydrogenase (GPDH), and the K120 side chain is also positioned to form a hydrogen bond to the C-1 hydroxyl of DHAP. The kinetic parameters for unactivated and phosphite dianion-activated GPDH-catalyzed reduction of glycolaldehyde and acetaldehyde (AcA) show that the transition state for the former reaction is stabilized by ca 5 kcal/mole by interactions of the C-1 hydroxyl group with the protein catalyst. The K120A and K204A substitutions at wild-type GPDH result in similar decreases in k_cat, but K_m is only affected by the K120A substitution. These results are consistent with 3 kcal/mol stabilizing interactions between the K120 or K204 side chains and a negative charge at the C-2 oxygen at the transition state for hydride transfer from NADH to DHAP. This stabilization resembles that observed at oxyanion holes for other enzymes. There is no detectable rescue of the K204A variant by ethylammonium cation (EtNH₃⁺), compared with the efficient rescue of the K120A variant. This is consistent with a difference in the accessibility of the variant enzyme active sites to exogenous EtNH₃⁺. The K120A/K204A substitutions cause a (6 × 10⁶)-fold increase in the promiscuity of wild-type hlGPDH for catalysis of the reduction of AcA compared to DHAP. This may reflect conservation of the active site for an ancestral alcohol dehydrogenase, whose relative activity for catalysis of reduction of AcA increases with substitutions that reduce the activity for reduction of the specific substrate DHAP.

The role of remote flavin adenine dinucleotide pieces in the oxidative decarboxylation catalyzed by salicylate hydroxylase

Salicylate hydroxylase (NahG) has a single redox site in which FAD is reduced by NADH, the O₂ is activated by the reduced flavin, and salicylate undergoes an oxidative decarboxylation by a C(4a)-hydroperoxyflavin intermediate to give catechol. We report experimental results that show the contribution of individual pieces of the FAD cofactor to the observed enzymatic activity for turnover of the whole cofactor. A comparison of the kinetic parameters and products for the NahG-catalyzed reactions of FMN and riboflavin cofactor fragments reveal that the adenosine monophosphate (AMP) and ribitol phosphate pieces of FAD act to anchor the flavin to the enzyme and to direct the partitioning of the C(4a)-hydroperoxyflavin reaction intermediate towards hydroxylation of salicylate. The addition of AMP or ribitol phosphate pieces to solutions of the truncated flavins results in a partial restoration of the enzymatic activity lost upon truncation of FAD, and the pieces direct the reaction of the C(4a)-hydroperoxyflavin intermediate towards hydroxylation of salicylate.

Protein-Ribofuranosyl Interactions Activate Orotidine 5'-Monophosphate Decarboxylase for Catalysis

The role of a global, substrate-driven, enzyme conformational change in enabling the extraordinarily large rate acceleration for orotidine 5'-monophosphate decarboxylase (OMPDC)-catalyzed decarboxylation of orotidine 5'-monophosphate (OMP) is examined in experiments that focus on the interactions between OMPDC and the ribosyl hydroxyl groups of OMP. The D37 and T100' side chains of OMPDC interact, respectively, with the C-3' and C-2' hydroxyl groups of enzyme-bound OMP. D37G and T100'A substitutions result in 1.4 kcal/mol increases in the activation barrier ΔG^⧧ for catalysis of decarboxylation of the phosphodianion-truncated substrate 1-(β-d-erythrofuranosyl)orotic acid (EO) but result in larger 2.1-2.9 kcal/mol increases in ΔG^⧧ for decarboxylation of OMP and for phosphite dianion-activated decarboxylation of EO. This shows that these substitutions reduce transition-state stabilization by the Q215, Y217, and R235 side chains at the dianion binding site. The D37G and T100'A substitutions result in <1.0 kcal/mol increases in ΔG^⧧ for activation of OMPDC-catalyzed decarboxylation of the phosphoribofuranosyl-truncated substrate FO by phosphite dianions. Experiments to probe the effect of D37 and T100' substitutions on the kinetic parameters for d-glycerol 3-phosphate and d-erythritol 4-phosphate activators of OMPDC-catalyzed decarboxylation of FO show that ΔG^⧧ for sugar phosphate-activated reactions is increased by ca. 2.5 kcal/mol for each -OH interaction eliminated by D37G or T100'A substitutions. We conclude that the interactions between the D37 and T100' side chains and ribosyl or ribosyl-like hydroxyl groups are utilized to activate OMPDC for catalysis of decarboxylation of OMP, EO, and FO.

Adenylate Kinase-Catalyzed Reaction of AMP in Pieces: Enzyme Activation for Phosphoryl Transfer to Phosphite Dianion

The binding of adenosine 5'-triphosphate (ATP) and adenosine 5'-monophosphate (AMP) to adenylate kinase (AdK) drives closure of lids over the substrate adenosyl groups. We test the hypothesis that this conformational change activates AdK for catalysis. The rate constants for Homo sapiens adenylate kinase 1 (HsAdK1)-catalyzed phosphoryl group transfer to AMP, k_cat/K_m = 7.0 × 10⁶ M^-1 s^-1, and phosphite dianion, (k_HPi)_obs ≤1 × 10^-4 M^-1 s^-1, show that the binding energy of the adenosyl group effects a ≥7.0 × 10¹⁰-fold rate acceleration of phosphoryl transfer from ATP. The third-order rate constant of k_cat/K_HPiK_EA = 260 M^-2 s^-1 for 1-(β-d-erythrofuranosyl)adenine (EA)-activated phosphoryl transfer to phosphite dianion was determined, and the isohypophosphate reaction product characterized by ³¹P NMR. The results demonstrate the following: (i) a ≥14.7 kcal/mol stabilization of the transition state for phosphoryl transfer by the adenosyl group of AMP and a ≥2.6 × 10⁶-fold rate acceleration from the EA-driven conformational change and (ii) the recovery of ≥8.7 kcal/mol of this transition state stabilization for EA-activated phosphoryl transfer from ATP to phosphite.

Linear Free Energy Relationships for Enzymatic Reactions: Fresh Insight from a Venerable Probe

Linear free energy relationships (LFERs) for substituent effects on reactions that proceed through similar transition states provide insight into transition state structures. A classical approach to the analysis of LFERs showed that differences in the slopes of Brønsted correlations for addition of substituted alkyl alcohols to ring-substituted 1-phenylethyl carbocations and to the β-galactopyranosyl carbocation intermediate of reactions catalyzed by β-galactosidase provide evidence that the enzyme catalyst modifies the curvature of the energy surface at the saddle point for the transition state for nucleophile addition. We have worked to generalize the use of LFERs in the determination of enzyme mechanisms. The defining property of enzyme catalysts is their specificity for binding the transition state with a much higher affinity than the substrate. Triosephosphate isomerase (TIM), orotidine 5′-monophosphate decarboxylase (OMPDC), and glycerol 3-phosphate dehydrogenase (GPDH) show effective catalysis of reactions of phosphorylated substrates and strong phosphite dianion activation of reactions of phosphodianion truncated substrates, with rate constants k_cat/K_m (M^–1 s^–1) and k_cat/K_dK_HPi (M^–2 s^–1), respectively. Good linear logarithmic correlations, with a slope of 1.1, between these kinetic parameters determined for reactions catalyzed by five or more variant forms of each catalyst are observed, where the protein substitutions are mainly at side chains which function to stabilize the cage complex between the enzyme and substrate. This shows that the enzyme-catalyzed reactions of a whole substrate and substrate pieces proceed through transition states of similar structures. It provides support for the proposal that the dianion binding energy of whole phosphodianion substrates and of phosphite dianion is used to drive the conversion of these protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes that show the same activity toward catalysis of the reactions of whole and phosphodianion truncated substrates. There is a good linear correlation, with a slope of 0.73, between values of the dissociation constants log K_i for release of the transition state analog phosphoglycolate (PGA) trianion and log k_cat/K_m for isomerization of GAP for wild-type and variants of TIM. This correlation shows that the substituted amino acid side chains act to stabilize the complex between TIM and the PGA trianion and that ca. 70% of this stabilization is observed at the transition state for substrate deprotonation. The correlation provides evidence that these side chains function to enhance the basicity of the E165 side chain of TIM, which deprotonates the bound carbon acid substrate. There is a good linear correlation, with a slope of 0.74, between the values of ΔG^‡ and ΔG° determined by electron valence bond (EVB) calculations to model deprotonation of dihydroxyacetone phosphate (DHAP) in water and when bound to wild-type and variant forms of TIM to form the enediolate reaction intermediate. This correlation provides evidence that the stabilizing interactions of the transition state for TIM-catalyzed deprotonation of DHAP are optimized by placement of amino acid side chains in positions that provide for the maximum stabilization of the charged reaction intermediate, relative to the neutral substrate.

Phosphodianion Activation of Enzymes for Catalysis of Central Metabolic Reactions

The activation barriers ΔG^⧧ for k_cat/K_m for the reactions of whole substrates catalyzed by 6-phosphogluconate dehydrogenase, glucose 6-phosphate dehydrogenase, and glucose 6-phosphate isomerase are reduced by 11–13 kcal/mol by interactions between the protein and the substrate phosphodianion. Between 4 and 6 kcal/mol of this dianion binding energy is expressed at the transition state for phosphite dianion activation of the respective enzyme-catalyzed reactions of truncated substrates d-xylonate or d-xylose. These and earlier results from studies on β-phosphoglucomutase, triosephosphate isomerase, and glycerol 3-phosphate dehydrogenase define a cluster of six enzymes that catalyze reactions in glycolysis or of glycolytic intermediates, and which utilize substrate dianion binding energy for enzyme activation. Dianion-driven conformational changes, which convert flexible open proteins to tight protein cages for the phosphorylated substrate, have been thoroughly documented for five of these six enzymes. The clustering of metabolic enzymes which couple phosphodianion-driven conformational changes to enzyme activation suggests that this catalytic motif has been widely propagated in the proteome.

The Organization of Active Site Side Chains of Glycerol-3-phosphate Dehydrogenase Promotes Efficient Enzyme Catalysis and Rescue of Variant Enzymes

A comparison of the values of k_cat/K_m for reduction of dihydroxyacetone phosphate (DHAP) by NADH catalyzed by wild type and K120A/R269A variant glycerol-3-phosphate dehydrogenase from human liver (hlGPDH) shows that the transition state for enzyme-catalyzed hydride transfer is stabilized by 12.0 kcal/mol by interactions with the cationic K120 and R269 side chains. The transition state for the K120A/R269A variant-catalyzed reduction of DHAP is stabilized by 1.0 and 3.8 kcal/mol for reactions in the presence of 1.0 M EtNH₃⁺ and guanidinium cation (Gua⁺), respectively, and by 7.5 kcal/mol for reactions in the presence of a mixture of each cation at 1.0 M, so that the transition state stabilization by the ternary E·EtNH₃⁺·Gua⁺ complex is 2.8 kcal/mol greater than the sum of stabilization by the respective binary complexes. This shows that there is cooperativity between the paired activators in transition state stabilization. The effective molarities (EMs) of ∼50 M determined for the K120A and R269A side chains are ≪10⁶ M, the EM for entropically controlled reactions. The unusually efficient rescue of the activity of hlGPDH-catalyzed reactions by the HP_i/Gua⁺ pair and by the Gua⁺/EtNH₃⁺ activator pair is due to stabilizing interactions between the protein and the activator pieces that organize the K120 and R269 side chains at the active site. This “preorganization” of side chains promotes effective catalysis by hlGPDH and many other enzymes. The role of the highly conserved network of side chains, which include Q295, R269, N270, N205, T264, K204, D260, and K120, in catalysis is discussed.

The role of ligand-gated conformational changes in enzyme catalysis

Structural and biochemical studies on diverse enzymes have highlighted the importance of ligand-gated conformational changes in enzyme catalysis, where the intrinsic binding energy of the common phosphoryl group of their substrates is used to drive energetically unfavorable conformational changes in catalytic loops, from inactive open to catalytically competent closed conformations. However, computational studies have historically been unable to capture the activating role of these conformational changes. Here, we discuss recent experimental and computational studies, which can remarkably pinpoint the role of ligand-gated conformational changes in enzyme catalysis, even when not modeling the loop dynamics explicitly. Finally, through our joint analyses of these data, we demonstrate how the synergy between theory and experiment is crucial for furthering our understanding of enzyme catalysis.

Uncovering the Role of Key Active-Site Side Chains in Catalysis: An Extended Brønsted Relationship for Substrate Deprotonation Catalyzed by Wild-Type and Variants of Triosephosphate Isomerase

We report results of detailed empirical valence bond simulations that model the effect of several amino acid substitutions on the thermodynamic (ΔG°) and kinetic activation (ΔG^⧧) barriers to deprotonation of dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) bound to wild-type triosephosphate isomerase (TIM), as well as to the K12G, E97A, E97D, E97Q, K12G/E97A, I170A, L230A, I170A/L230A, and P166A variants of this enzyme. The EVB simulations model the observed effect of the P166A mutation on protein structure. The E97A, E97Q, and E97D mutations of the conserved E97 side chain result in ≤1.0 kcal mol^-1 decreases in the activation barrier for substrate deprotonation. The agreement between experimental and computed activation barriers is within ±1 kcal mol^-1, with a strong linear correlation between ΔG^⧧ and ΔG° for all 11 variants, with slopes β = 0.73 (R² = 0.994) and β = 0.74 (R² = 0.995) for the deprotonation of DHAP and GAP, respectively. These Brønsted-type correlations show that the amino acid side chains examined in this study function to reduce the standard-state Gibbs free energy of reaction for deprotonation of the weak α-carbonyl carbon acid substrate to form the enediolate phosphate reaction intermediate. TIM utilizes the cationic side chain of K12 to provide direct electrostatic stabilization of the enolate oxyanion, and the nonpolar side chains of P166, I170, and L230 are utilized for the construction of an active-site cavity that provides optimal stabilization of the enediolate phosphate intermediate relative to the carbon acid substrate.

Active-Site Glu165 Activation in Triosephosphate Isomerase and Its Deprotonation Kinetics

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) via an enediol(ate) intermediate. The active-site residue Glu165 serves as the catalytic base during catalysis. It abstracts a proton from C1 carbon of DHAP to form the reaction intermediate and donates a proton to C2 carbon of the intermediate to form product GAP. Our difference Fourier transform infrared spectroscopy studies on the yeast TIM (YeTIM)/phosphate complex revealed a C═O stretch band at 1706 cm^-1 from the protonated Glu165 carboxyl group at pH 7.5, indicating that the p K_a of the catalytic base is increased by >3.0 pH units upon phosphate binding, and that the Glu165 carboxyl environment in the complex is still hydrophilic in spite of the increased p K_a. Hence, the results show that the binding of the phosphodianion group is part of the activation mechanism which involves the p K_a elevation of the catalytic base Glu165. The deprotonation kinetics of Glu165 in the μs to ms time range were determined via infrared (IR) T-jump studies on the YeTIM/phosphate and ("heavy enzyme") [U-¹³C,-¹⁵N]YeTIM/phosphate complexes. The slower deprotonation kinetics in the ms time scale is due to phosphate dissociation modulated by the loop motion, which slows down by enzyme mass increase to show a normal heavy enzyme kinetic isotope effect (KIE) ∼1.2 (i.e., slower rate in the heavy enzyme). The faster deprotonation kinetics in the tens of μs time scale is assigned to temperature-induced p K_a decrease, while phosphate is still bound, and it shows an inverse heavy enzyme KIE ∼0.89 (faster rate in the heavy enzyme). The IR static and T-jump spectroscopy provides atomic-level resolution of the catalytic mechanism because of its ability to directly observe the bond breaking/forming process.

Protein Flexibility and Stiffness Enable Efficient Enzymatic Catalysis

The enormous rate accelerations observed for many enzyme catalysts are due to strong stabilizing interactions between the protein and reaction transition state. The defining property of these catalysts is their specificity for binding the transition state with a much higher affinity than substrate. Experimental results are presented which show that the phosphodianion-binding energy of phosphate monoester substrates is used to drive conversion of their protein catalysts from flexible and entropically rich ground states to stiff and catalytically active Michaelis complexes. These results are generalized to other enzyme-catalyzed reactions. The existence of many enzymes in flexible, entropically rich, and inactive ground states provides a mechanism for utilization of ligand-binding energy to mold these catalysts into stiff and active forms. This reduces the substrate-binding energy expressed at the Michaelis complex, while enabling the full and specific expression of large transition-state binding energies. Evidence is presented that the complexity of enzyme conformational changes increases with increases in the enzymatic rate acceleration. The requirement that a large fraction of the total substrate-binding energy be utilized to drive conformational changes of floppy enzymes is proposed to favor the selection and evolution of protein folds with multiple flexible unstructured loops, such as the TIM-barrel fold. The effect of protein motions on the kinetic parameters for enzymes that undergo ligand-driven conformational changes is considered. The results of computational studies to model the complex ligand-driven conformational change in catalysis by triosephosphate isomerase are presented.

Human Glycerol 3-Phosphate Dehydrogenase: X-ray Crystal Structures That Guide the Interpretation of Mutagenesis Studies

Human liver glycerol 3-phosphate dehydrogenase ( hlGPDH) catalyzes the reduction of dihydroxyacetone phosphate (DHAP) to form glycerol 3-phosphate, using the binding energy associated with the nonreacting phosphodianion of the substrate to properly orient the enzyme-substrate complex within the active site. Herein, we report the crystal structures for unliganded, binary E·NAD, and ternary E·NAD·DHAP complexes of wild type hlGPDH, illustrating a new position of DHAP, and probe the kinetics of multiple mutant enzymes with natural and truncated substrates. Mutation of Lys120, which is positioned to donate a proton to the carbonyl of DHAP, results in similar increases in the activation barrier to hlGPDH-catlyzed reduction of DHAP and to phosphite dianion-activated reduction of glycolaldehyde, illustrating that these transition states show similar interactions with the cationic K120 side chain. The K120A mutation results in a 5.3 kcal/mol transition state destabilization, and 3.0 kcal/mol of the lost transition state stabilization is rescued by 1.0 M ethylammonium cation. The 6.5 kcal/mol increase in the activation barrier observed for the D260G mutant hlGPDH-catalyzed reaction represents a 3.5 kcal/mol weakening of transition state stabilization by the K120A side chain and a 3.0 kcal/mol weakening of the interactions with other residues. The interactions, at the enzyme active site, between the K120 side chain and the Q295 and R269 side chains were likewise examined by double-mutant analyses. These results provide strong evidence that the enzyme rate acceleration is due mainly or exclusively to transition state stabilization by electrostatic interactions with polar amino acid side chains.

Simultaneous biosynthesis of (R)-acetoin and ethylene glycol from D-xylose through in vitro metabolic engineering

(R)-acetoin is a four-carbon platform compound used as the precursor for synthesizing novel optically active materials. Ethylene glycol (EG) is a large-volume two-carbon commodity chemical used as the anti-freezing agent and building-block molecule for various polymers. Currently established microbial fermentation processes for converting monosaccharides to either (R)-acetoin or EG are plagued by the formation of undesirable by-products. We show here that a cell-free bioreaction scheme can generate enantiomerically pure acetoin and EG as co-products from biomass-derived D-xylose. The seven-step, ATP-free system included in situ cofactor regeneration and recruited enzymes from Escherichia coli W3110, Bacillus subtilis shaijiu 32 and Caulobacter crescentus CB 2. Optimized in vitro biocatalytic conditions generated 3.2 mM (R)-acetoin with stereoisomeric purity of 99.5% from 10 mM D-xylose at 30 °C and pH 7.5 after 24 h, with an initial (R)-acetoin productivity of 1.0 mM/h. Concomitantly, EG was produced at 5.5 mM, with an initial productivity of 1.7 mM/h. This in vitro biocatalytic platform illustrates the potential for production of multiple value-added biomolecules from biomass-based sugars with no ATP requirement.

Loop Motion in Triosephosphate Isomerase Is Not a Simple Open and Shut Case

Conformational changes are crucial for the catalytic action of many enzymes. A prototypical and well-studied example is loop opening and closure in triosephosphate isomerase (TIM), which is thought to determine the rate of catalytic turnover in many circumstances. Specifically, TIM loop 6 "grips" the phosphodianion of the substrate and, together with a change in loop 7, sets up the TIM active site for efficient catalysis. Crystal structures of TIM typically show an open or a closed conformation of loop 6, with the tip of the loop moving ∼7 Å between conformations. Many studies have interpreted this motion as a two-state, rigid-body transition. Here, we use extensive molecular dynamics simulations, with both conventional and enhanced sampling techniques, to analyze loop motion in apo and substrate-bound TIM in detail, using five crystal structures of the dimeric TIM from Saccharomyces cerevisiae. We find that loop 6 is highly flexible and samples multiple conformational states. Empirical valence bond simulations of the first reaction step show that slight displacements away from the fully closed-loop conformation can be sufficient to abolish most of the catalytic activity; full closure is required for efficient reaction. The conformational change of the loops in TIM is thus not a simple "open and shut" case and is crucial for its catalytic action. Our detailed analysis of loop motion in a highly efficient enzyme highlights the complexity of loop conformational changes and their role in biological catalysis.

Primary Deuterium Kinetic Isotope Effects: A Probe for the Origin of the Rate Acceleration for Hydride Transfer Catalyzed by Glycerol-3-Phosphate Dehydrogenase

Large primary deuterium kinetic isotope effects (1° DKIEs) on enzyme-catalyzed hydride transfer may be observed when the transferred hydride tunnels through the energy barrier. The following 1° DKIEs on k_cat/K_m and relative reaction driving force are reported for wild-type and mutant glycerol-3-phosphate dehydrogenase (GPDH)-catalyzed reactions of NADL (L = H, D): wild-type GPDH, ΔΔG^⧧ = 0 kcal/mol, 1° DKIE = 1.5; N270A, 5.6 kcal/mol, 3.1; R269A, 9.1 kcal/mol, 2.8; R269A + 1.0 M guanidine, 2.4 kcal/mol, 2.7; R269A/N270A, 11.5 kcal/mol, 2.4. Similar 1° DKIEs were observed on k_cat. The narrow range of 1° DKIEs (2.4–3.1) observed for a 9.1 kcal/mol change in reaction driving force provides strong evidence that these are intrinsic 1° DKIEs on rate-determining hydride transfer. Evidence is presented that the intrinsic DKIE on wild-type GPDH-catalyzed reduction of DHAP lies in this range. A similar range of 1° DKIEs (2.4–2.9) on (k_cat/K_GA, M^–1 s^–1) was reported for dianion-activated hydride transfer from NADL to glycolaldehyde (GA) [Reyes, A. C.; Amyes, T. L.; Richard, J. P. J. Am. Chem. Soc.2016, 138, 14526–14529]. These 1° DKIEs are much smaller than those observed for enzyme-catalyzed hydrogen transfer that occurs mainly by quantum mechanical tunneling. These results support the conclusion that the rate acceleration for GPDH-catalyzed reactions is due to the stabilization of the transition state for hydride transfer by interactions with the protein catalyst. The small 1° DKIEs reported for mutant GPDH-catalyzed and for wild-type dianion-activated reactions are inconsistent with a model where the dianion binding energy is utilized in the stabilization of a tunneling ready state.

Role of Ligand-Driven Conformational Changes in Enzyme Catalysis: Modeling the Reactivity of the Catalytic Cage of Triosephosphate Isomerase

We have previously performed empirical valence bond calculations of the kinetic activation barriers, ΔG^‡_calc, for the deprotonation of complexes between TIM and the whole substrate glyceraldehyde-3-phosphate (GAP, Kulkarni et al.J. Am. Chem. Soc.2017, 139, 10514–1052528683550). We now extend this work to also study the deprotonation of the substrate pieces glycolaldehyde (GA) and GA·HP_i [HP_i = phosphite dianion]. Our combined calculations provide activation barriers, ΔG^‡_calc, for the TIM-catalyzed deprotonation of GAP (12.9 ± 0.8 kcal·mol^–1), of the substrate piece GA (15.0 ± 2.4 kcal·mol^–1), and of the pieces GA·HP_i (15.5 ± 3.5 kcal·mol^–1). The effect of bound dianion on ΔG^‡_calc is small (≤2.6 kcal·mol^–1), in comparison to the much larger 12.0 and 5.8 kcal·mol^–1 intrinsic phosphodianion and phosphite dianion binding energy utilized to stabilize the transition states for TIM-catalyzed deprotonation of GAP and GA·HP_i, respectively. This shows that the dianion binding energy is essentially fully expressed at our protein model for the Michaelis complex, where it is utilized to drive an activating change in enzyme conformation. The results represent an example of the synergistic use of results from experiments and calculations to advance our understanding of enzymatic reaction mechanisms.

Enzyme Architecture: The Role of a Flexible Loop in Activation of Glycerol-3-phosphate Dehydrogenase for Catalysis of Hydride Transfer

The side chain of Q295 of glycerol-3-phosphate dehydrogenase from human liver (hlGPDH) lies in a flexible loop, that folds over the phosphodianion of substrate dihydroxyacetone phosphate (DHAP). Q295 interacts with the side-chain cation from R269, which is ion-paired to the substrate phosphodianion. Kinetic parameters k_cat/K_m (M^–1 s^–1) and k_cat/K_GAK_HPi (M^–2 s^–1) were determined, respectively, for catalysis of the reduction of DHAP and for dianion activation of catalysis of reduction of glycolaldehyde (GA) catalyzed by wild-type, Q295G, Q295S, Q295A, and Q295N mutants of hlGPDH. These mutations result in up to a 150-fold decrease in (k_cat/K_m)_DHAP and up to a 2.7 kcal/mol decrease in the intrinsic phosphodianion binding energy. The data define a linear correlation with slope 1.1, between the intrinsic phosphodianion binding energy and the intrinsic phosphite dianion binding energy for activation of hlGPDH-catalyzed reduction of GA, that demonstrates a role for Q295 in optimizing this dianion binding energy. The R269A mutation of wild-type GPDH results in a 9.1 kcal/mol destabilization of the transition state for reduction of DHAP, but the same R269A mutation of N270A and Q295A mutants result in smaller 5.9 and 4.9 kcal/mol transition-state destabilization. Similarly, the N270A or Q295A mutations of R269A GPDH each result in large falloffs in the efficiency of rescue of the R269A mutant by guanidine cation. We conclude that N270, which interacts for the substrate phosphodianion and Q295, which interacts with the guanidine side chain of R269, function to optimize the apparent transition-state stabilization provided by the cationic side chain of R269.

A theoretical study on the size-dependence of ground-state proton transfer in phenol-ammonia clusters

Geometries and infrared (IR) spectra in the mid-IR region of phenol-(ammonia)_n (PhOH-(NH₃)_n) (n = 0-10) clusters have been studied using density functional theory (DFT) to investigate the critical number of solvent molecules necessary to promote ground-state proton transfer (GSPT). For n ≤ 8 clusters, the most stable isomer is a non-proton-transferred (non-PT) structure, and all isomers found within 1.5 kcal mol^-1 from it are also non-PT structures. For n = 9, the most stable isomer is also a non-PT structure; however, the second stable isomer is a PT structure, whose relative energy is within the experimental criterion of population (0.7 kcal mol^-1). For n = 10, the PT structure is the most stable one. We can therefore estimate that the critical size of GSPT is n = 9. This is confirmed by the fact that these calculated IR spectra are in good accordance with our previous experimental results of mid-IR spectra. It is demonstrated that characteristic changes of the ν_9a and ν₁₂ bands in the skeletal vibrational region provide clear information that the GSPT reaction has occurred. It was also found that the shortest distance between the π-ring and the solvent moiety is a good indicator of the PT reaction.

Difference FTIR Studies of Substrate Distribution in Triosephosphate Isomerase

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP), via an enediol(ate) intermediate. Determination of substrate population distribution in the TIM/substrate reaction mixture at equilibrium and characterization of the substrate-enzyme interactions in the Michaelis complex are ongoing efforts toward the understanding of the TIM reaction mechanism. By using isotope-edited difference Fourier transform infrared studies with unlabeled and 13C-labeled substrates at specific carbon(s), we are able to show that in the reaction mixture at equilibrium the keto DHAP is the dominant species and the populations of aldehyde GAP and enediol(ate) are very low, consistent with the results from previous X-ray structural and 13C NMR studies. Furthermore, within the DHAP side of the Michaelis complex, there is a set of conformational substates that can be characterized by the different C2═O stretch frequencies. The C2═O frequency differences reflect the different degree of the C2═O bond polarization due to hydrogen bonding from active site residues. The C2═O bond polarization has been considered as an important component for substrate activation within the Michaelis complex. We have found that in the enzyme-substrate reaction mixture with TIM from different organisms the number of substates and their population distribution within the DHAP side of the Michaelis complex may be different. These discoveries provide a rare opportunity to probe the interconversion dynamics of these DHAP substates and form the bases for the future studies to determine if the TIM-catalyzed reaction follows a simple linear reaction pathway, as previously believed, or follows parallel reaction pathways, as suggested in another enzyme system that also shows a set of substates in the Michaelis complex.

Enzyme Architecture: Modeling the Operation of a Hydrophobic Clamp in Catalysis by Triosephosphate Isomerase

Triosephosphate isomerase (TIM) is a proficient catalyst of the reversible isomerization of dihydroxyacetone phosphate (DHAP) to d-glyceraldehyde phosphate (GAP), via general base catalysis by E165. Historically, this enzyme has been an extremely important model system for understanding the fundamentals of biological catalysis. TIM is activated through an energetically demanding conformational change, which helps position the side chains of two key hydrophobic residues (I170 and L230), over the carboxylate side chain of E165. This is critical both for creating a hydrophobic pocket for the catalytic base and for maintaining correct active site architecture. Truncation of these residues to alanine causes significant falloffs in TIM's catalytic activity, but experiments have failed to provide a full description of the action of this clamp in promoting substrate deprotonation. We perform here detailed empirical valence bond calculations of the TIM-catalyzed deprotonation of DHAP and GAP by both wild-type TIM and its I170A, L230A, and I170A/L230A mutants, obtaining exceptional quantitative agreement with experiment. Our calculations provide a linear free energy relationship, with slope 0.8, between the activation barriers and Gibbs free energies for these TIM-catalyzed reactions. We conclude that these clamping side chains minimize the Gibbs free energy for substrate deprotonation, and that the effects on reaction driving force are largely expressed at the transition state for proton transfer. Our combined analysis of previous experimental and current computational results allows us to provide an overview of the breakdown of ground-state and transition state effects in enzyme catalysis in unprecedented detail, providing a molecular description of the operation of a hydrophobic clamp in triosephosphate isomerase.

Enzyme activation through the utilization of intrinsic dianion binding energy

43

We consider 'the proposition that the intrinsic binding energy that results from the noncovalent interaction of a specific substrate with the active site of the enzyme is considerably larger than is generally believed. An important part of this binding energy may be utilized to provide the driving force for catalysis, so that the observed binding energy represents only what is left over after this utilization' [Jencks,W.P. (1975) Adv. Enzymol. Relat. Areas. Mol. Biol. , , 219-410]. The large ~12 kcal/mol intrinsic substrate phosphodianion binding energy for reactions catalyzed by triosephosphate isomerase (TIM), orotidine 5'-monophosphate decarboxylase and glycerol-3-phosphate dehydrogenase is divided into 4-6 kcal/mol binding energy that is expressed on the formation of the Michaelis complex in anchoring substrates to the respective enzyme, and 6-8 kcal/mol binding energy that is specifically expressed at the transition state in activating the respective enzymes for catalysis. A structure-based mechanism is described where the dianion binding energy drives a conformational change that activates these enzymes for catalysis. Phosphite dianion plays the active role of holding TIM in a high-energy closed active form, but acts as passive spectator in showing no effect on transition-state structure. The result of studies on mutant enzymes is presented, which support the proposal that the dianion-driven enzyme conformational change plays a role in enhancing the basicity of side chain of E167, the catalytic base, by clamping the base between a pair of hydrophobic side chains. The insight these results provide into the architecture of enzyme active sites and the development of strategies for the de novo design of protein catalysts is discussed.

Enzyme Architecture: Self-Assembly of Enzyme and Substrate Pieces of Glycerol-3-Phosphate Dehydrogenase into a Robust Catalyst of Hydride Transfer

The stabilization of the transition state for hlGPDH-catalyzed reduction of DHAP due to the action of the phosphodianion of DHAP and the cationic side chain of R269 is between 12.4 and 17 kcal/mol. The R269A mutation of glycerol-3-phosphate dehydrogenase (hlGPDH) results in a 9.1 kcal/mol destabilization of the transition state for enzyme-catalyzed reduction of dihydroxyacetone phosphate (DHAP) by NADH, and there is a 6.7 kcal/mol stabilization of this transition state by 1.0 M guanidine cation (Gua⁺) [J. Am. Chem. Soc.2015, 137, 5312–5315]. The R269A mutant shows no detectable activity toward reduction of glycolaldehyde (GA), or activation of this reaction by 30 mM HPO₃^2–. We report the unprecedented self-assembly of R269A hlGPDH, dianions (X^2– = FPO₃^2–, HPO₃^2–, or SO₄^2–), Gua⁺ and GA into a functioning catalyst of the reduction of GA, and fourth-order reaction rate constants k_cat/K_GAK_XK_Gua. The linear logarithmic correlation (slope = 1.0) between values of k_cat/K_GAK_X for dianion activation of wildtype hlGPDH-catalyzed reduction of GA and k_cat/K_GAK_XK_Gua shows that the electrostatic interaction between exogenous dianions and the side chain of R269 is not significantly perturbed by cutting hlGPDH into R269A and Gua⁺ pieces. The advantage for connection of hlGPDH (R269A mutant + Gua⁺) and substrate pieces (GA + HP_i) pieces, (ΔG_S^‡)_HPi+E+Gua = 5.6 kcal/mol, is nearly equal to the sum of the advantage to connection of the substrate pieces, (ΔG_S^‡)_GA+HPi = 3.3 kcal/mol, for wildtype hlGPDH-catalyzed reaction of GA + HP_i, and for connection of the enzyme pieces, (ΔG_S^‡)_E+Gua = 2.4 kcal/mol, for Gua⁺ activation of the R269A hlGPDH-catalyzed reaction of DHAP.

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 (k_cat/K_GA, M^–1 s^–1) for dianion (X^2–) activated hydride transfer from NADL to glycolaldehyde (GA) catalyzed by glycerol-3-phosphate dehydrogenase were determined over a 2100-fold range of enzyme reactivity: (X^2–, 1°DKIE); FPO₃^2–, 2.8 ± 0.1; HPO₃^2–, 2.5 ± 0.1; SO₄^2–, 2.8 ± 0.2; HOPO₃^2–, 2.5 ± 0.1; S₂O₃^2–, 2.9 ± 0.1; unactivated; 2.4 ± 0.2. Similar 1°DKIEs were determined for k_cat. 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 X^2– are similar, but there are large differences in the transition state interactions for different X^2–. The changes in transition state binding interactions are expressed as changes in k_cat 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.

Structure-Function Studies of Hydrophobic Residues That Clamp a Basic Glutamate Side Chain during Catalysis by Triosephosphate Isomerase

Kinetic parameters are reported for the reactions of whole substrates (k_cat/K_m, M^–1 s^–1) (R)-glyceraldehyde 3-phosphate (GAP) and dihydroxyacetone phosphate (DHAP) and for the substrate pieces [(k_cat/K_m)_E·HPi/K_d, M^–2 s^–1] glycolaldehyde (GA) and phosphite dianion (HP_i) catalyzed by the I172A/L232A mutant of triosephosphate isomerase from Trypanosoma brucei brucei (TbbTIM). A comparison with the corresponding parameters for wild-type, I172A, and L232A TbbTIM-catalyzed reactions shows that the effect of I172A and L232A mutations on ΔG^⧧ for the wild-type TbbTIM-catalyzed reactions of the substrate pieces is nearly the same as the effect of the same mutations on TbbTIM previously mutated at the second side chain. This provides strong evidence that mutation of the first hydrophobic side chain does not affect the functioning of the second side chain in catalysis of the reactions of the substrate pieces. By contrast, the effects of I172A and L232A mutations on ΔG^⧧ for wild-type TbbTIM-catalyzed reactions of the whole substrate are different from the effect of the same mutations on TbbTIM previously mutated at the second side chain. This is due to the change in the rate-determining step that determines the barrier to the isomerization reaction. X-ray crystal structures are reported for I172A, L232A, and I172A/L232A TIMs and for the complexes of these mutants to the intermediate analogue phosphoglycolate (PGA). The structures of the PGA complexes with wild-type and mutant enzymes are nearly superimposable, except that the space opened by replacement of the hydrophobic side chain is occupied by a water molecule that lies ∼3.5 Å from the basic side chain of Glu167. The new water at I172A mutant TbbTIM provides a simple rationalization for the increase in the activation barrier ΔG^⧧ observed for mutant enzyme-catalyzed reactions of the whole substrate and substrate pieces. By contrast, the new water at the L232A mutant does not predict the decrease in ΔG^⧧ observed for the mutant enzyme-catalyzed reactions of the substrate piece GA.

Connecting Active-Site Loop Conformations and Catalysis in Triosephosphate Isomerase: Insights from a Rare Variation at Residue 96 in the Plasmodial Enzyme

Despite extensive research into triosephosphate isomerases (TIMs), there exists a gap in understanding of the remarkable conjunction between catalytic loop-6 (residues 166-176) movement and the conformational flip of Glu165 (catalytic base) upon substrate binding that primes the active site for efficient catalysis. The overwhelming occurrence of serine at position 96 (98% of the 6277 unique TIM sequences), spatially proximal to E165 and the loop-6 residues, raises questions about its role in catalysis. Notably, Plasmodium falciparum TIM has an extremely rare residue--phenylalanine--at this position whereas, curiously, the mutant F96S was catalytically defective. We have obtained insights into the influence of residue 96 on the loop-6 conformational flip and E165 positioning by combining kinetic and structural studies on the PfTIM F96 mutants F96Y, F96A, F96S/S73A, and F96S/L167V with sequence conservation analysis and comparative analysis of the available apo and holo structures of the enzyme from diverse organisms.

Role of Loop-Clamping Side Chains in Catalysis by Triosephosphate Isomerase

The side chains of Y208 and S211 from loop 7 of triosephosphate isomerase (TIM) form hydrogen bonds to backbone amides and carbonyls from loop 6 to stabilize the caged enzyme–substrate complex. The effect of seven mutations [Y208T, Y208S, Y208A, Y208F, S211G, S211A, Y208T/S211G] on the kinetic parameters for TIM catalyzed reactions of the whole substrates dihydroxyacetone phosphate and d-glyceraldehyde 3-phosphate [(k_cat/K_m)_GAP and (k_cat/K_m)_DHAP] and of the substrate pieces glycolaldehyde and phosphite dianion (k_cat/K_HPiK_GA) are reported. The linear logarithmic correlation between these kinetic parameters, with slope of 1.04 ± 0.03, shows that most mutations of TIM result in an identical change in the activation barriers for the catalyzed reactions of whole substrate and substrate pieces, so that the transition states for these reactions are stabilized by similar interactions with the protein catalyst. The second linear logarithmic correlation [slope = 0.53 ± 0.16] between k_cat for isomerization of GAP and K_d^⧧ for phosphite dianion binding to the transition state for wildtype and many mutant TIM-catalyzed reactions of substrate pieces shows that ca. 50% of the wildtype TIM dianion binding energy, eliminated by these mutations, is expressed at the wildtype Michaelis complex, and ca. 50% is only expressed at the wildtype transition state. Negative deviations from this correlation are observed when the mutation results in a decrease in enzyme reactivity at the catalytic site. The main effect of Y208T, Y208S, and Y208A mutations is to cause a reduction in the total intrinsic dianion binding energy, but the effect of Y208F extends to the catalytic site.

Rate and Equilibrium Constants for an Enzyme Conformational Change during Catalysis by Orotidine 5'-Monophosphate Decarboxylase

The caged complex between orotidine 5'-monophosphate decarboxylase (ScOMPDC) and 5-fluoroorotidine 5'-monophosphate (FOMP) undergoes decarboxylation ∼300 times faster than the caged complex between ScOMPDC and the physiological substrate, orotidine 5'-monophosphate (OMP). Consequently, the enzyme conformational changes required to lock FOMP at a protein cage and release product 5-fluorouridine 5'-monophosphate (FUMP) are kinetically significant steps. The caged form of ScOMPDC is stabilized by interactions between the side chains from Gln215, Tyr217, and Arg235 and the substrate phosphodianion. The control of these interactions over the barrier to the binding of FOMP and the release of FUMP was probed by determining the effect of all combinations of single, double, and triple Q215A, Y217F, and R235A mutations on kcat/Km and kcat for turnover of FOMP by wild-type ScOMPDC; its values are limited by the rates of substrate binding and product release, respectively. The Q215A and Y217F mutations each result in an increase in kcat and a decrease in kcat/Km, due to a weakening of the protein-phosphodianion interactions that favor fast product release and slow substrate binding. The Q215A/R235A mutation causes a large decrease in the kinetic parameters for ScOMPDC-catalyzed decarboxylation of OMP, which are limited by the rate of the decarboxylation step, but much smaller decreases in the kinetic parameters for ScOMPDC-catalyzed decarboxylation of FOMP, which are limited by the rate of enzyme conformational changes. By contrast, the Y217A mutation results in large decreases in kcat/Km for ScOMPDC-catalyzed decarboxylation of both OMP and FOMP, because of the comparable effects of this mutation on rate-determining decarboxylation of enzyme-bound OMP and on the rate-determining enzyme conformational change for decarboxylation of FOMP. We propose that kcat = 8.2 s(-1) for decarboxylation of FOMP by the Y217A mutant is equal to the rate constant for cage formation from the complex between FOMP and the open enzyme, that the tyrosyl phenol group stabilizes the closed form of ScOMPDC by hydrogen bonding to the substrate phosphodianion, and that the phenyl group of Y217 and F217 facilitates formation of the transition state for the rate-limiting conformational change. An analysis of kinetic data for mutant enzyme-catalyzed decarboxylation of OMP and FOMP provides estimates for the rate and equilibrium constants for the conformational change that traps FOMP at the enzyme active site.

Enzyme architecture: optimization of transition state stabilization from a cation-phosphodianion pair

The side chain cation of R269 lies at the surface of l-glycerol 3-phosphate dehydrogenase (GPDH) and forms an ion pair to the phosphodianion of substrate dihydroxyacetone phosphate (DHAP), which is buried at the nonpolar protein interior. The R269A mutation of GPDH results in a 110-fold increase in K_m (2.8 kcal/mol effect) and a 41 000-fold decrease in k_cat (6.3 kcal/mol effect), which corresponds to a 9.1 kcal/mol destabilization of the transition state for GPDH-catalyzed reduction of DHAP by NADH. There is a 6.7 kcal/mol stabilization of the transition state for the R269A mutant GPDH-catalyzed reaction by 1.0 M guanidinium ion, and the transition state for the reaction of the substrate pieces is stabilized by an additional 2.4 kcal/mol by their covalent attachment at wildtype GPDH. These results provide strong support for the proposal that GPDH invests the 11 kcal/mol intrinsic phosphodianion binding energy of DHAP in trapping the substrate at a nonpolar active site, where strong electrostatic interactions are favored, and obtains a 9 kcal/mol return from stabilizing interactions between the side chain cation and transition state trianion. We propose a wide propagation for the catalytic motif examined in this work, which enables strong transition state stabilization from enzyme–phosphodianion pairs.

The role of phosphate in a multistep enzymatic reaction: reactions of the substrate and intermediate in pieces

Several mechanistically unrelated enzymes utilize the binding energy of their substrate's nonreacting phosphoryl group to accelerate catalysis. Evidence for the involvement of the phosphodianion in transition state formation has come from reactions of the substrate in pieces, in which reaction of a truncated substrate lacking its phosphorylmethyl group is activated by inorganic phosphite. What has remained unknown until now is how the phosphodianion group influences the reaction energetics at different points along the reaction coordinate. 1-Deoxy-D-xylulose-5-phosphate (DXP) reductoisomerase (DXR), which catalyzes the isomerization of DXP to 2-C-methyl-D-erythrose 4-phosphate (MEsP) and subsequent NADPH-dependent reduction, presents a unique opportunity to address this concern. Previously, we have reported the effect of covalently linked phosphate on the energetics of DXP turnover. Through the use of chemically synthesized MEsP and its phosphate-truncated analogue, 2-C-methyl-D-glyceraldehyde, the current study revealed a loss of 6.1 kcal/mol of kinetic barrier stabilization upon truncation, of which 4.4 kcal/mol was regained in the presence of phosphite dianion. The activating effect of phosphite was accompanied by apparent tightening of its interactions within the active site at the intermediate stage of the reaction, suggesting a role of the phosphodianion in disfavoring intermediate release and in modulation of the on-enzyme isomerization equilibrium. The results of kinetic isotope effect and structural studies indicate rate limitation by physical steps when the covalent linkage is severed. These striking differences in the energetics of the natural reaction and the reactions in pieces provide a deeper insight into the contribution of enzyme-phosphodianion interactions to the reaction coordinate.

The activating oxydianion binding domain for enzyme-catalyzed proton transfer, hydride transfer, and decarboxylation: specificity and enzyme architecture

The kinetic parameters for activation of yeast triosephosphate isomerase (ScTIM), yeast orotidine monophosphate decarboxylase (ScOMPDC), and human liver glycerol 3-phosphate dehydrogenase (hlGPDH) for catalysis of reactions of their respective phosphodianion truncated substrates are reported for the following oxydianions: HPO3(2-), FPO3(2-), S2O3(2-), SO4(2-) and HOPO3(2-). Oxydianions bind weakly to these unliganded enzymes and tightly to the transition state complex (E·S(‡)), with intrinsic oxydianion Gibbs binding free energies that range from -8.4 kcal/mol for activation of hlGPDH-catalyzed reduction of glycolaldehyde by FPO3(2-) to -3.0 kcal/mol for activation of ScOMPDC-catalyzed decarboxylation of 1-β-d-erythrofuranosyl)orotic acid by HOPO3(2-). Small differences in the specificity of the different oxydianion binding domains are observed. We propose that the large -8.4 kcal/mol and small -3.8 kcal/mol intrinsic oxydianion binding energy for activation of hlGPDH by FPO3(2-) and S2O3(2-), respectively, compared with activation of ScTIM and ScOMPDC reflect stabilizing and destabilizing interactions between the oxydianion -F and -S with the cationic side chain of R269 for hlGPDH. These results are consistent with a cryptic function for the similarly structured oxydianion binding domains of ScTIM, ScOMPDC and hlGPDH. Each enzyme utilizes the interactions with tetrahedral inorganic oxydianions to drive a conformational change that locks the substrate in a caged Michaelis complex that provides optimal stabilization of the different enzymatic transition states. The observation of dianion activation by stabilization of active caged Michaelis complexes may be generalized to the many other enzymes that utilize substrate binding energy to drive changes in enzyme conformation, which induce tight substrate fits.

α-Fluorophosphonates reveal how a phosphomutase conserves transition state conformation over hexose recognition in its two-step reaction

β-Phosphoglucomutase (βPGM) catalyzes isomerization of β-D-glucose 1-phosphate (βG1P) into D-glucose 6-phosphate (G6P) via sequential phosphoryl transfer steps using a β-D-glucose 1,6-bisphosphate (βG16BP) intermediate. Synthetic fluoromethylenephosphonate and methylenephosphonate analogs of βG1P deliver novel step 1 transition state analog (TSA) complexes for βPGM, incorporating trifluoromagnesate and tetrafluoroaluminate surrogates of the phosphoryl group. Within an invariant protein conformation, the β-D-glucopyranose ring in the βG1P TSA complexes (step 1) is flipped over and shifted relative to the G6P TSA complexes (step 2). Its equatorial hydroxyl groups are hydrogen-bonded directly to the enzyme rather than indirectly via water molecules as in step 2. The (C)O-P bond orientation for binding the phosphate in the inert phosphate site differs by ∼ 30° between steps 1 and 2. By contrast, the orientations for the axial O-Mg-O alignment for the TSA of the phosphoryl group in the catalytic site differ by only ∼ 5°, and the atoms representing the five phosphorus-bonded oxygens in the two transition states (TSs) are virtually superimposable. The conformation of βG16BP in step 1 does not fit into the same invariant active site for step 2 by simple positional interchange of the phosphates: the TS alignment is achieved by conformational change of the hexose rather than the protein.

Enzyme architecture: deconstruction of the enzyme-activating phosphodianion interactions of orotidine 5'-monophosphate decarboxylase

The mechanism for activation of orotidine 5'-monophosphate decarboxylase (OMPDC) by interactions of side chains from Gln215 and Try217 at a gripper loop and R235, adjacent to this loop, with the phosphodianion of OMP was probed by determining the kinetic parameters k(cat) and K(m) for all combinations of single, double, and triple Q215A, Y217F, and R235A mutations. The 12 kcal/mol intrinsic binding energy of the phosphodianion is shown to be equal to the sum of the binding energies of the side chains of R235 (6 kcal/mol), Q215 (2 kcal/mol), Y217 (2 kcal/mol), and hydrogen bonds to the G234 and R235 backbone amides (2 kcal/mol). Analysis of a triple mutant cube shows small (ca. 1 kcal/mol) interactions between phosphodianion gripper side chains, which are consistent with steric crowding of the side chains around the phosphodianion at wild-type OMPDC. These mutations result in the same change in the activation barrier to the OMPDC-catalyzed reactions of the whole substrate OMP and the substrate pieces (1-β-D-erythrofuranosyl)orotic acid (EO) and phosphite dianion. This shows that the transition states for these reactions are stabilized by similar interactions with the protein catalyst. The 12 kcal/mol intrinsic phosphodianion binding energy of OMP is divided between the 8 kcal/mol of binding energy, which is utilized to drive a thermodynamically unfavorable conformational change of the free enzyme, resulting in an increase in (k(cat))(obs) for OMPDC-catalyzed decarboxylation of OMP, and the 4 kcal/mol of binding energy, which is utilized to stabilize the Michaelis complex, resulting in a decrease in (K(m))(obs).

Reflections on the catalytic power of a TIM-barrel

The TIM-barrel fold is described and its propagation throughout the enzyme universe noted. The functions of the individual front loops of the eponymous TIM-barrel of triosephosphate isomerase are presented in a discussion of: (a) electrophilic catalysis, by amino acid side chains from loops 1 and 4, of abstraction of an α-carbonyl hydrogen from substrate dihydroxyacetone phosphate (DHAP) or d-glyceraldehyde 3-phosphate (DGAP). (b) The engineering of loop 3 to give the monomeric variant monoTIM and the structure and catalytic properties of this monomer. (c) The interaction between loops 6, 7 and 8 and the phosphodianion of DHAP or DGAP. (d) The mechanism by which a ligand-gated conformational change, dominated by motion of loops 6 and 7, activates TIM for catalysis of deprotonation of DHAP or DGAP. (e) The conformational plasticity of TIM, and the utilization of substrate binding energy to "mold" the distorted active site loops of TIM mutants into catalytically active enzymes. The features of the TIM-barrel fold that favor effective protein catalysis are discussed.

Enzyme architecture: the effect of replacement and deletion mutations of loop 6 on catalysis by triosephosphate isomerase

Two mutations of the phosphodianion gripper loop in chicken muscle triosephosphate isomerase (cTIM) were examined: (1) the loop deletion mutant (LDM) formed by removal of residues 170-173 [Pompliano, D. L., et al. (1990) Biochemistry 29, 3186-3194] and (2) the loop 6 replacement mutant (L6RM), in which the N-terminal hinge sequence of TIM from eukaryotes, 166-PXW-168 (X = L or V), is replaced by the sequence from archaea, 166-PPE-168. The X-ray crystal structure of the L6RM shows a large displacement of the side chain of E168 from that for W168 in wild-type cTIM. Solution nuclear magnetic resonance data show that the L6RM results in significant chemical shift changes in loop 6 and surrounding regions, and that the binding of glycerol 3-phosphate (G3P) results in chemical shift changes for nuclei at the active site of the L6RM that are smaller than those of wild-type cTIM. Interactions with loop 6 of the L6RM stabilize the enediolate intermediate toward the elimination reaction catalyzed by the LDM. The LDM and L6RM result in 800000- and 23000-fold decreases, respectively, in kcat/Km for isomerization of GAP. Saturation of the LDM, but not the L6RM, by substrate and inhibitor phosphoglycolate is detected by steady-state kinetic analyses. We propose, on the basis of a comparison of X-ray crystal structures for wild-type TIM and the L6RM, that ligands bind weakly to the L6RM because a large fraction of the ligand binding energy is utilized to overcome destabilizing electrostatic interactions between the side chains of E168 and E129 that are predicted to develop in the loop-closed enzyme. Similar normalized yields of DHAP, d-DHAP, and d-GAP are formed in LDM- and L6RM-catalyzed reactions of GAP in D2O. The smaller normalized 12-13% yield of DHAP and d-DHAP observed for the mutant cTIM-catalyzed reactions compared with the 79% yield of these products for wild-type cTIM suggests that these mutations impair the transfer of a proton from O-2 to O-1 at the initial enediolate phosphate intermediate. No products are detected for the LDM-catalyzed isomerization reactions in D2O of [1-(13)C]GA and HPi, but the L6RM-catalyzed reaction in the presence of 0.020 M dianion gives a 2% yield of the isomerization product [2-(13)C,2-(2)H]GA.

Mechanistic Imperatives for Deprotonation of Carbon Catalyzed by Triosephosphate Isomerase: Enzyme-Activation by Phosphite Dianion

The mechanistic imperatives for catalysis of deprotonation of α-carbonyl carbon by triosephosphate isomerase (TIM) are discussed. There is a strong imperative to reduce the large thermodynamic barrier for deprotonation of carbon to form an enediolate reaction intermediate; and, a strong imperative for specificity in the expression of the intrinsic phosphodianion binding energy at the transition state for the enzyme-catalyzed reaction. Binding energies of 2 and 6 kcal/mol, respectively, have been determined for formation of phosphite dianion complexes to TIM and to the transition state for TIM-catalyzed deprotonation of the truncated substrate glycolaldehyde [T. L. Amyes, J. P. Richard, Biochemistry2007, 46, 5841]. We propose that the phosphite dianion binding energy, which is specifically expressed at the transition state complex, is utilized to stabilize a rare catalytically active loop-closed form of TIM. The results of experiments to probe the role of the side chains of Ile172 and Leu232 in activating the loop-closed form of TIM for catalysis of substrate deprotonation are discussed. Evidence is presented that the hydrophobic side chain of Ile172 assists in activating TIM for catalysis of substrate deprotonation through an enhancement of the basicity of the carboxylate side-chain of Glu167. Our experiments link the two imperatives for TIM-catalyzed deprotonation of carbon by providing evidence that the phosphodianion binding energy is utilized to drive an enzyme conformational change, which results in a reduction in the thermodynamic barrier to deprotonation of the carbon acid substrate at TIM compared with the barrier for deprotonation in water. The effects of a P168A mutation on the kinetic parameters for the reactions of whole and truncated substrates are discussed.

Enzyme architecture: remarkably similar transition states for triosephosphate isomerase-catalyzed reactions of the whole substrate and the substrate in pieces

Values of (k_cat/K_m)_GAP for triosephosphate isomerase-catalyzed reactions of (R)-glyceraldehyde 3-phosphate and k_cat/K_HPiK_GA for reactions of the substrate pieces glycolaldehyde and HPO₃^2– have been determined for wild-type and the following TIM mutants: I172V, I172A, L232A, and P168A (TIM from Trypanosoma brucei brucei); a 208-TGAG for 208-YGGS loop 7 replacement mutant (L7RM, TIM from chicken muscle); and, Y208T, Y208S, Y208A, Y208F and S211A (yeast TIM). A superb linear logarithmic correlation, with slope of 1.04 ± 0.03, is observed between the kinetic parameters for wild-type and most mutant enzymes, with positive deviations for L232A and L7RM. The unit slope shows that most mutations result in an identical change in the activation barriers for the catalyzed reactions of whole substrate and substrate pieces, so that the two transition states are stabilized by similar interactions with the protein catalyst. This is consistent with a role for dianions as active spectators, which hold TIM in a catalytically active caged form.

Study of orotidine 5'-monophosphate decarboxylase in complex with the top three OMP, BMP, and PMP ligands by molecular dynamics simulation

Catalytic mechanism of orotidine 5'-monophosphate decarboxylase (OMPDC), one of the nature most proficient enzymes which provides large rate enhancement, has not been fully understood yet. A series of 30 ns molecular dynamics (MD) simulations were run on X-ray structure of the OMPDC from Saccharomyces cerevisiae in its free form as well as in complex with different ligands, namely 1-(5'-phospho-D-ribofuranosyl) barbituric acid (BMP), orotidine 5'-monophosphate (OMP), and 6-phosphonouridine 5'-monophosphate (PMP). The importance of this biological system is justified both by its high rate enhancement and its potential use as a target in chemotherapy. This work focuses on comparing two physicochemical states of the enzyme (protonated and deprotonated Asp91) and three ligands (substrate OMP, inhibitor, and transition state analog BMP and substrate analog PMP). Detailed analysis of the active site geometry and its interactions is properly put in context by extensive comparison with relevant experimental works. Our overall results show that in terms of hydrogen bond occupancy, electrostatic interactions, dihedral angles, active site configuration, and movement of loops, notable differences among different complexes are observed. Comparison of the results obtained from these simulations provides some detailed structural data for the complexes, the enzyme, and the ligands, as well as useful insights into the inhibition mechanism of the OMPDC enzyme. Furthermore, these simulations are applied to clarify the ambiguous mechanism of the OMPDC enzyme, and imply that the substrate destabilization and transition state stabilization contribute to the mechanism of action of the most proficient enzyme, OMPDC.