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

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].

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.

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.

Atomic Force Microscopy Study of the Temperature and Storage Duration Dependencies of Horseradish Peroxidase Oligomeric State

This paper presents an investigation of the temperature dependence of the oligomeric state of the horseradish peroxidase (HRP) enzyme on the temperature of its solution, and on the solution storage time, at the single-molecule level. Atomic force microscopy has been employed to determine how the temperature and the storage time of the HRP solution influence its aggregation upon direct adsorption of the enzyme from the solution onto bare mica substrates. In parallel, spectrophotometric measurements have been performed in order to estimate whether the HRP enzymatic activity changes over time upon the storage of the enzyme solution. The temperature dependence of the HRP oligomeric state has been studied within a broad (15–40 °C) temperature range. It has been demonstrated that the storage of the HRP solution for 14 days does not have any considerable effect on the oligomeric state of the enzyme, neither does it affect its activity. At longer storage times, AFM has allowed us to reveal a tendency of HRP to oligomerization during the storage of its buffered solution, while the enzymatic activity remains virtually unchanged even after a 1-month-long storage. By AFM, it has been revealed that after the incubation of a mica substrate in the HRP solution at various temperatures, the content of the mica-adsorbed oligomers increases insignificantly owing to a high-temperature stability of the enzyme.

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.

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.

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.

Orotidine 5'-Monophosphate Decarboxylase: The Operation of Active Site Chains Within and Across Protein Subunits

The D37 and T100′ side chains of orotidine 5′-monophosphate decarboxylase (OMPDC) interact with the C-3′ and C-2′ ribosyl hydroxyl groups, respectively, of the bound substrate. We compare the intra-subunit interactions of D37 with the inter-subunit interactions of T100′ by determining the effects of the D37G, D37A, T100′G, and T100′A substitutions on the following: (a) k_cat and k_cat/K_m values for the OMPDC-catalyzed decarboxylations of OMP and 5-fluoroorotidine 5′-monophosphate (FOMP) and (b) the stability of dimeric OMPDC relative to the monomer. The D37G and T100′A substitutions resulted in 2 kcal mol^–1 increases in ΔG^† for k_cat/K_m for the decarboxylation of OMP, while the D37A and T100′G substitutions resulted in larger 4 and 5 kcal mol^–1 increases, respectively, in ΔG^†. The D37G and T100′A substitutions both resulted in smaller 2 kcal mol^–1 decreases in ΔG^† for the decarboxylation of FOMP compared to that of OMP. These results show that the D37G and T100′A substitutions affect the barrier to the chemical decarboxylation step while the D37A and T100′G substitutions also affect the barrier to a slow, ligand-driven enzyme conformational change. Substrate binding induces the movement of an α-helix (G′98–S′106) toward the substrate C-2′ ribosyl hydroxy bound at the main subunit. The T100′G substitution destabilizes the enzyme dimer by 3.5 kcal mol^–1 compared to the monomer, which is consistent with the known destabilization of α-helices by the internal Gly side chains [Serrano, L., et al. (1992) Nature, 356, 453–455]. We propose that the T100′G substitution weakens the α-helical contacts at the dimer interface, which results in a decrease in the dimer stability and an increase in the barrier to the ligand-driven conformational change.

Role of the Carboxylate in Enzyme-Catalyzed Decarboxylation of Orotidine 5'-Monophosphate: Transition State Stabilization Dominates Over Ground State Destabilization

Kinetic parameters k_ex (s^–1) and k_ex/K_d (M^–1 s^–1) are reported for exchange for deuterium in D₂O of the C-6 hydrogen of 5-fluororotidine 5′-monophosphate (FUMP) catalyzed by the Q215A, Y217F, and Q215A/Y217F variants of yeast orotidine 5′-monophosphate decarboxylase (ScOMPDC) at pD 8.1, and by the Q215A variant at pD 7.1–9.3. The pD rate profiles for wildtype ScOMPDC and the Q215A variant are identical, except for a 2.5 log unit downward displacement in the profile for the Q215A variant. The Q215A, Y217F and Q215A/Y217F substitutions cause 1.3–2.0 kcal/mol larger increases in the activation barrier for wildtype ScOMPDC-catalyzed deuterium exchange compared with decarboxylation, because of the stronger apparent side chain interaction with the transition state for the deuterium exchange reaction. The stabilization of the transition state for the OMPDC-catalyzed deuterium exchange reaction of FUMP is ca. 19 kcal/mol smaller than the transition state for decarboxylation of OMP, and ca. 8 kcal/mol smaller than for OMPDC-catalyzed deprotonation of FUMP to form the vinyl carbanion intermediate common to OMPDC-catalyzed reactions OMP/FOMP and UMP/FUMP. We propose that ScOMPDC shows similar stabilizing interactions with the common portions of decarboxylation and deprotonation transition states that lead to formation of this vinyl carbanion intermediate, and that there is a large ca. (19–8) = 11 kcal/mol stabilization of the former transition state from interactions with the nascent CO₂ of product. The effects of Q215A and Y217F substitutions on k_cat/K_m for decarboxylation of OMP are expressed mainly as an increase in K_m for the reactions catalyzed by the variant enzymes, while the effects on k_ex/K_d for deuterium exchange are expressed mainly as an increase in k_ex. This shows that the Q215 and Y217 side chains stabilize the Michaelis complex to OMP for the decarboxylation reaction, compared with the complex to FUMP for the deuterium exchange reaction. These results provide strong support for the conclusion that interactions which stabilize the transition state for ScOMPDC-catalyzed decarboxylation at a nonpolar enzyme active site dominate over interactions that destabilize the ground-state Michaelis complex.

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.

Enzyme Architecture: Breaking Down the Catalytic Cage that Activates Orotidine 5'-Monophosphate Decarboxylase for Catalysis

We report the results of a study of the catalytic role of a network of four interacting amino acid side chains at yeast orotidine 5'-monophosphate decarboxylase ( ScOMPDC), by the stepwise replacement of all four side chains. The H-bond, which links the -CH₂OH side chain of S154 from the pyrimidine umbrella loop of ScOMPDC to the amide side chain of Q215 in the phosphodianion gripper loop, creates a protein cage for the substrate OMP. The role of this interaction in optimizing transition state stabilization from the dianion gripper side chains Q215, Y217, and R235 was probed by determining the kinetic parameter k_cat/ K_m for 16 enzyme variants, which include all combinations of single, double, triple, and quadruple S154A, Q215A, Y217F, and R235A mutations. The effects of consecutive Q215A, Y217F, and R235A mutations on Δ G^⧧ for wild-type enzyme-catalyzed decarboxylation sum to 11.6 kcal/mol, but to only 7.6 kcal/mol when starting from S154A mutant. This shows that the S154A mutation results in a (11.6-7.6) = 4.0 kcal/mol decrease in transition state stabilization from interactions with Q215, Y217, and R235. Mutant cycles show that ca. 2 kcal/mol of this 4 kcal/mol effect is from the direct interaction between the S154 and Q215 side chains and that ca. 2 kcal/mol is from a tightening in the stabilizing interactions of the Y217 and R235 side chains. The sum of the effects of individual A154S, A215Q, F217Y and A235R substitutions at the quadruple mutant of ScOMPDC to give the corresponding triple mutants, 5.6 kcal/mol, is much smaller than 16.0 kcal/mol, the sum of the effects of the related four substitutions in wild-type ScOMPDC to give the respective single mutants. The small effect of substitutions at the quadruple mutant is consistent with a large entropic cost to holding the flexible loops of ScOMPDC in the active closed conformation.

Orotidine 5'-Monophosphate Decarboxylase: Probing the Limits of the Possible for Enzyme Catalysis

The mystery associated with catalysis by what were once regarded as protein black boxes, diminished with the X-ray crystallographic determination of the three-dimensional structures of enzyme–substrate complexes. The report that several high-resolution X-ray crystal structures of orotidine 5′-monophosphate decarboxylase (OMPDC) failed to provide a consensus mechanism for enzyme-catalyzed decarboxylation of OMP to form uridine 5′-monophosphate, therefore, provoked a flurry of controversy. This controversy was fueled by the enormous 10²³-fold rate acceleration for this enzyme, which had “jolted many biochemists’ assumptions about the catalytic potential of enzymes.” Our studies on the mechanism of action of OMPDC provide strong evidence that catalysis by this enzyme is not fundamentally different from less proficient catalysts, while highlighting important architectural elements that enable a peak level of performance. Many enzymes undergo substrate-induced protein conformational changes that trap their substrates in solvent occluded protein cages, but the conformational change induced by ligand binding to OMPDC is incredibly complex, as required to enable the development of 22 kcal/mol of stabilizing binding interactions with the phosphodianion and ribosyl substrate fragments of OMP. The binding energy from these fragments is utilized to activate OMPDC for catalysis of decarboxylation at the orotate fragment of OMP, through the creation of a tight, catalytically active, protein cage from the floppy, open, unliganded form of OMPDC. Such utilization of binding energy for ligand-driven conformational changes provides a general mechanism to obtain specificity in transition state binding. The rate enhancement that results from the binding of carbon acid substrates to enzymes is partly due to a reduction in the carbon acid pK_a that is associated with ligand binding. The binding of UMP to OMPDC results in an unusually large >12 unit decrease in the pK_a = 29 for abstraction of the C-6 substrate hydrogen, due to stabilization of an enzyme-bound vinyl carbanion, which is also an intermediate of OMPDC-catalyzed decarboxylation. The protein–ligand interactions operate to stabilize the vinyl carbanion at the enzyme active site compared to aqueous solution, rather than to stabilize the transition state for the concerted electrophilic displacement of CO₂ by H⁺ that avoids formation of this reaction intermediate. There is evidence that OMPDC induces strain into the bound substrate. The interaction between the amide side chain of Gln-215 from the phosphodianion gripper loop and the hydroxymethylene side chain of Ser-154 from the pyrimidine umbrella of ScOMPDC position the amide side chain to interact with the phosphodianion of OMP. There are no direct stabilizing interactions between dianion gripper protein side chains Gln-215, Tyr-217, and Arg-235 and the pyrimidine ring at the decarboxylation transition state. Rather these side chains function solely to hold OMPDC in the catalytically active closed conformation. The hydrophobic side chains that line the active site of OMPDC in the region of the departing CO₂ product may function to stabilize the decarboxylation transition state by providing hydrophobic solvation of this product.

Enzyme Architecture: Erection of Active Orotidine 5'-Monophosphate Decarboxylase by Substrate-Induced Conformational Changes

Orotidine 5′-monophosphate decarboxylase (OMPDC) catalyzes the decarboxylation of 5-fluoroorotate (FO) with k_cat/K_m = 1.4 × 10^–7 M^–1 s^–1. Combining this and related kinetic parameters shows that the 31 kcal/mol stabilization of the transition state for decarboxylation of OMP provided by OMPDC represents the sum of 11.8 and 10.6 kcal/mol stabilization by the substrate phosphodianion and the ribosyl ring, respectively, and an 8.6 kcal/mol stabilization from the orotate ring. The transition state for OMPDC-catalyzed decarboxylation of FO is stabilized by 5.2, 7.2, and 9.0 kcal/mol, respectively, by 1.0 M phosphite dianion, d-glycerol 3-phosphate and d-erythritol 4-phosphate. The stabilization is due to the utilization of binding interactions of the substrate fragments to drive an enzyme conformational change, which locks the orotate ring of the whole substrate, or the substrate pieces in a caged complex. We propose that enzyme-activation is a possible, and perhaps probable, consequence of any substrate-induced enzyme conformational change.

Utilization of Substrate Intrinsic Binding Energy for Conformational Change and Catalytic Function in Phosphoenolpyruvate Carboxykinase

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Previous work has demonstrated that the enzyme cycles between a catalytically inactive open state and a catalytically active closed state. The transition of the enzyme between these states requires the transition of several active site loops to shift from mobile, disordered structural elements to stable ordered states. The mechanism by which these disorder-order transitions are coupled to the ligation state of the active site however is not fully understood. To further investigate the mechanisms by which the mobility of the active site loops is coupled to enzymatic function and the transitioning of the enzyme between the two conformational states, we have conducted structural and functional studies of point mutants of E89. E89 is a proposed key member of the interaction network of mobile elements as it resides in the R-loop region of the enzyme active site. These new data demonstrate the importance of the R-loop in coordinating interactions between substrates at the OAA/PEP binding site and the mobile R- and Ω-loop domains. In turn, the studies more generally demonstrate the mechanisms by which the intrinsic ligand binding energy can be utilized in catalysis to drive unfavorable conformational changes, changes that are subsequently required for both optimal catalytic activity and fidelity.

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.

Linking coupled motions and entropic effects to the catalytic activity of 2-deoxyribose-5-phosphate aldolase (DERA)

Local mutations in the phosphate binding group of DERA alter global conformation dynamics, catalytic activities and reaction entropies.

DERA, 2-deoxyribose-5-phosphate aldolase, catalyzes the retro-aldol cleavage of 2-deoxy-ribose-5-phosphate (dR5P) into glyceraldehyde-3-phosphate (G3P) and acetaldehyde in a branch of the pentose phosphate pathway. In addition to the physiological reaction, DERA also catalyzes the reverse addition reaction and, hence, is an interesting candidate for bio-catalysis of carbo-ligation reactions, which are central to synthetic chemistry. An obstacle to overcome for this enzyme to become a truly useful biocatalyst, however, is to relax the very strict dependency of this enzyme on phosphorylated substrates. We have studied herein the role of the non-canonical phosphate-binding site of this enzyme, consisting of Ser²³⁸ and Ser²³⁹, by site-directed and site-saturation mutagenesis, coupled to kinetic analysis of mutants. In addition, we have performed molecular dynamics simulations on the wild-type and four mutant enzymes, to analyse how mutations at this phosphate-binding site may affect the protein structure and dynamics. Further examination of the S239P mutant revealed that this variant increases the enthalpy change at the transition state, relative to the wild-type enzyme, but concomitant loss in entropy causes an overall relative loss in the TS free energy change. This entropy loss, as measured by the temperature dependence of catalysed rates, was mirrored in both a drastic loss in dynamics of the enzyme, which contributes to phosphate binding, as well as an overall loss in anti-correlated motions distributed over the entire protein. Our combined data suggests that the degree of anticorrelated motions within the DERA structure is coupled to catalytic efficiency in the DERA-catalyzed retro-aldol cleavage reaction, and can be manipulated for engineering purposes.