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

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.

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.

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.

An Enzyme with High Catalytic Proficiency Utilizes Distal Site Substrate Binding Energy to Stabilize the Closed State but at the Expense of Substrate Inhibition

Understanding the factors that underpin the enormous catalytic proficiencies of enzymes is fundamental to catalysis and enzyme design. Enzymes are, in part, able to achieve high catalytic proficiencies by utilizing the binding energy derived from nonreacting portions of the substrate. In particular, enzymes with substrates containing a nonreacting phosphodianion group coordinated in a distal site have been suggested to exploit this binding energy primarily to facilitate a conformational change from an open inactive form to a closed active form, rather than to either induce ground state destabilization or stabilize the transition state. However, detailed structural evidence for the model is limited. Here, we use β-phosphoglucomutase (βPGM) to investigate the relationship between binding a phosphodianion group in a distal site, the adoption of a closed enzyme form, and catalytic proficiency. βPGM catalyzes the isomerization of β-glucose 1-phosphate to glucose 6-phosphate via phosphoryl transfer reactions in the proximal site, while coordinating a phosphodianion group of the substrate(s) in a distal site. βPGM has one of the largest catalytic proficiencies measured and undergoes significant domain closure during its catalytic cycle. We find that side chain substitution at the distal site results in decreased substrate binding that destabilizes the closed active form but is not sufficient to preclude the adoption of a fully closed, near-transition state conformation. Furthermore, we reveal that binding of a phosphodianion group in the distal site stimulates domain closure even in the absence of a transferring phosphoryl group in the proximal site, explaining the previously reported β-glucose 1-phosphate inhibition. Finally, our results support a trend whereby enzymes with high catalytic proficiencies involving phosphorylated substrates exhibit a greater requirement to stabilize the closed active form.

One-Pot and Two-Chamber Methodologies for Using Acetylene Surrogates in the Synthesis of Pyridazines and Their D-Labeled Derivatives

Acetylene surrogates are efficient tools in modern organic chemistry with largely unexplored potential in the construction of heterocyclic cores. Two novel synthetic paths to 3,6-disubstituted pyridazines were proposed using readily available acetylene surrogates through flexible C₂ unit installation procedures in a common reaction space mode (one-pot) and distributed reaction space mode (two-chamber): (1) an interaction of 1,2,4,5-tetrazine and its acceptor-functionalized derivatives with a CaC₂ -H₂ O mixture performed in a two-chamber reactor led to the corresponding pyridazines in quantitative yields; (2) [4+2] cycloaddition of 1,2,4,5-tetrazines to benzyl vinyl ether can be considered a universal synthetic path to a wide range of pyridazines. Replacing water with D₂ O and vinyl ether with its trideuterated analog in the developed procedures, a range of 4,5-dideuteropyridazines of 95-99% deuteration degree was synthesized for the first time. Quantum chemical modeling allowed to quantify the substituent effect in both synthetic pathways.

Hydride Transfer Catalyzed by Glycerol Phosphate Dehydrogenase: Recruitment of an Acidic Amino Acid Side Chain to Rescue a Damaged Enzyme

K120 of glycerol 3-phosphate dehydrogenase (GPDH) lies close to the carbonyl group of the bound dihydroxyacetone phosphate (DHAP) dianion. pH rate (pH 4.6–9.0) profiles are reported for k_cat and (k_cat/K_m)_dianion for wild type and K120A GPDH-catalyzed reduction of DHAP by NADH, and for (k_cat/K_dK_am) for activation of the variant-catalyzed reduction by CH₃CH₂NH₃⁺, where K_am and K_d are apparent dissociation constants for CH₃CH₂NH₃⁺ and DHAP, respectively. These profiles provide evidence that the K120 side chain cation, which is stabilized by an ion-pairing interaction with the D260 side chain, remains protonated between pH 4.6 and 9.0. The profiles for wild type and K120A variant GPDH show downward breaks at a similar pH value (7.6) that are attributed to protonation of the K204 side chain, which also lies close to the substrate carbonyl oxygen. The pH profiles for (k_cat/K_m)_dianion and (k_cat/K_dK_am) for the K120A variant show that the monoprotonated form of the variant is activated for catalysis by CH₃CH₂NH₃⁺ but has no detectable activity, compared to the diprotonated variant, for unactivated reduction of DHAP. The pH profile for k_cat shows that the monoprotonated K120A variant is active toward reduction of enzyme-bound DHAP, because of activation by a ligand-driven conformational change. Upward breaks in the pH profiles for k_cat and (k_cat/K_m)_dianion for K120A GPDH are attributed to protonation of D260. These breaks are consistent with the functional replacement of K120 by D260, and a plasticity in the catalytic roles of the active site side chains.

Modeling the Role of a Flexible Loop and Active Site Side Chains in Hydride Transfer Catalyzed by Glycerol-3-phosphate Dehydrogenase

Glycerol-3-phosphate dehydrogenase is a biomedically important enzyme that plays a crucial role in lipid biosynthesis. It is activated by a ligand-gated conformational change that is necessary for the enzyme to reach a catalytically competent conformation capable of efficient transition-state stabilization. While the human form (hlGPDH) has been the subject of extensive structural and biochemical studies, corresponding computational studies to support and extend experimental observations have been lacking. We perform here detailed empirical valence bond and Hamiltonian replica exchange molecular dynamics simulations of wild-type hlGPDH and its variants, as well as providing a crystal structure of the binary hlGPDH·NAD R269A variant where the enzyme is present in the open conformation. We estimated the activation free energies for the hydride transfer reaction in wild-type and substituted hlGPDH and investigated the effect of mutations on catalysis from a detailed structural study. In particular, the K120A and R269A variants increase both the volume and solvent exposure of the active site, with concomitant loss of catalytic activity. In addition, the R269 side chain interacts with both the Q295 side chain on the catalytic loop, and the substrate phosphodianion. Our structural data and simulations illustrate the critical role of this side chain in facilitating the closure of hlGPDH into a catalytically competent conformation, through modulating the flexibility of a key catalytic loop (292-LNGQKL-297). This, in turn, rationalizes a tremendous 41,000 fold decrease experimentally in the turnover number, k_cat, upon truncating this residue, as loop closure is essential for both correct positioning of key catalytic residues in the active site, as well as sequestering the active site from the solvent. Taken together, our data highlight the importance of this ligand-gated conformational change in catalysis, a feature that can be exploited both for protein engineering and for the design of allosteric inhibitors targeting this biomedically important enzyme.

A robust and tunable halogen bond organocatalyzed 2-deoxyglycosylation involving quantum tunneling

The development of noncovalent halogen bonding (XB) catalysis is rapidly gaining traction, as isolated reports documented better performance than the well-established hydrogen bonding thiourea catalysis. However, convincing cases allowing XB activation to be competitive in challenging bond formations are lacking. Herein, we report a robust XB catalyzed 2-deoxyglycosylation, featuring a biomimetic reaction network indicative of dynamic XB activation. Benchmarking studies uncovered an improved substrate tolerance compared to thiourea-catalyzed protocols. Kinetic investigations reveal an autoinductive sigmoidal kinetic profile, supporting an in situ amplification of a XB dependent active catalytic species. Kinetic isotopic effect measurements further support quantum tunneling in the rate determining step. Furthermore, we demonstrate XB catalysis tunability via a halogen swapping strategy, facilitating 2-deoxyribosylations of D-ribals. This protocol showcases the clear emergence of XB catalysis as a versatile activation mode in noncovalent organocatalysis, and as an important addition to the catalytic toolbox of chemical glycosylations.

Halogen bonding (HB) catalysis is rapidly gaining momentum, however, cases of XB activation for challenging bonds formation are rare. Here, the authors show a robust XB catalyzed 2-deoxyglycosylation with broad scope and featuring a quantum tunneling phenomenon in the proton transfer rate determining step.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg²⁺ ions, and recent advances in their transient involvement in the catalytic mechanisms.

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.