Mechanism for activation of triosephosphate isomerase by phosphite dianion: the role of a ligand-driven conformational change

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.

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.

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.

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.

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.

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.

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.

'Negative' and 'positive catalysis': complementary principles that shape the catalytic landscape of enzymes

Our understanding of enzyme catalysis is dominated by transition state theory. According to this concept, an enzymatic reaction is guided along a desired reaction coordinate through the stabilization of favorable transition state. But how much is the outcome of an enzyme reaction controlled by the destabilization of unwanted transition states? Here, we revive and critically review the hypothesis that the active site of enzymes also features elements of 'negative catalysis'. We provide examples that show that enzyme catalysis can be achieved by the combined action of positive and negative constraints at the active site of an enzyme. This integrated view of enzyme catalysis has direct consequences for our studies on the catalytic landscape of enzymes, as well as current efforts in enzyme engineering and the de novo-design of enzymes.

Enzyme Architecture: Amino Acid Side-Chains That Function To Optimize the Basicity of the Active Site Glutamate of Triosephosphate Isomerase

We report pH rate profiles for k_cat and K_m for the isomerization reaction of glyceraldehyde 3-phosphate catalyzed by wildtype triosephosphate isomerase (TIM) from three organisms and by ten mutants of TIM; and, for K_i for inhibition of this reaction by phosphoglycolate trianion (I^3–). The pH profiles for K_i show that the binding of I^3– to TIM (E) to form EH·I₃^– is accompanied by uptake of a proton by the carboxylate side-chain of E165, whose function is to abstract a proton from substrate. The complexes for several mutants exist mainly as E^–·I₃^– at high pH, in which cases the pH profiles define the pK_a for deprotonation of EH·I₃^–. The linear free energy correlation, with slope of 0.73 (r² = 0.96), between k_cat/K_m for TIM-catalyzed isomerization and the disassociation constant of PGA trianion for TIM shows that EH·I₃^– and the transition state are stabilized by similar interactions with the protein catalyst. Values of pK_a = 10–10.5 were estimated for deprotonation of EH·I₃^– for wildtype TIM. This pK_a decreases to as low as 6.3 for the severely crippled Y208F mutant. There is a correlation between the effect of several mutations on k_cat/K_m and on pK_a for EH·I₃^–. The results support a model where the strong basicity of E165 at the complex to the enediolate reaction intermediate is promoted by side-chains from Y208 and S211, which serve to clamp loop 6 over the substrate; I170, which assists in the creation of a hydrophobic environment for E165; and P166, which functions in driving the carboxylate side-chain of E165 toward enzyme-bound substrate.

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.

A reevaluation of the origin of the rate acceleration for enzyme-catalyzed hydride transfer

There is no consensus of opinion on the origin of the large rate accelerations observed for enzyme-catalyzed hydride transfer. The interpretation of recent results from studies on hydride transfer reactions catalyzed by alcohol dehydrogenase (ADH) focus on the proposal that the effective barrier height is reduced by quantum-mechanical tunneling through the energy barrier. This interpretation contrasts sharply with the notion that enzymatic rate accelerations are obtained through direct stabilization of the transition state for the nonenzymatic reaction in water. The binding energy of the dianion of substrate DHAP provides 11 kcal mol^-1 stabilization of the transition state for the hydride transfer reaction catalyzed by glycerol-3-phosphate dehydrogenase (GPDH). We summarize evidence that the binding interactions between (GPDH) and dianion activators are utilized directly for stabilization of the transition state for enzyme-catalyzed hydride transfer. The possibility is considered, and then discounted, that these dianion binding interactions are utilized for the stabilization of a tunnel ready state (TRS) that enables efficient tunneling of the transferred hydride through the energy barrier, and underneath the energy maximum for the transition state. It is noted that the evidence to support the existence of a tunnel-ready state for the hydride transfer reactions catalyzed by ADH is ambiguous. We propose that the rate acceleration for ADH is due to the utilization of the binding energy of the cofactor NAD+/NADH in the stabilization of the transition state for enzyme-catalyzed hydride transfer.

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.

CADEE: Computer-Aided Directed Evolution of Enzymes

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

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.

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.

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.

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: on the importance of being in a protein cage

Substrate binding occludes water from the active sites of many enzymes. There is a correlation between the burden to enzymatic catalysis of deprotonation of carbon acids and the substrate immobilization at solvent-occluded active sites for ketosteroid isomerase (KSI--small burden, substrate pKa=13), triosephosphate isomerase (TIM, substrate pKa≈18) and diaminopimelate epimerase (DAP epimerase, large burden, substrate pKa≈29) catalyzed reaction. KSI binds substrates at a surface cleft, TIM binds substrate at an exposed 'cage' formed by closure of flexible loops; and, DAP epimerase binds substrate in a tight cage formed by an 'oyster-like' clamping motion of protein domains. Directed evolution of a solvent-occluded active site at a designed protein catalyst of the Kemp elimination reaction is 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.

Enzyme architecture: the activating oxydianion binding domain for orotidine 5'-monophophate decarboxylase

Orotidine 5'-monophosphate decarboxylase catalyzes the decarboxylation of truncated substrate (1-β-D-erythrofuranosyl)orotic acid to form (1-β-D-erythrofuranosyl)uracil. This enzyme-catalyzed reaction is activated by tetrahedral oxydianions, which bind weakly to unliganded OMPDC and tightly to the enzyme-transition state complex, with the following intrinsic oxydianion binding energies (kcal/mol): SO3(2-), -8.3; HPO3(2-), -7.7; S2O3(2-), -4.6; SO4(2-), -4.5; HOPO3(2-), -3.0; HOAsO3(2-), no activation detected. We propose that the oxydianion and orotate binding domains of OMPDC perform complementary functions in catalysis of decarboxylation reactions: (1) The orotate binding domain carries out decarboxylation of the orotate ring. (2) The activating oxydianion binding domain has the cryptic function of utilizing binding interactions with tetrahedral inorganic oxydianions to drive an enzyme conformational change that results in the stabilization of transition states at the distant orotate domain.

Structural mutations that probe the interactions between the catalytic and dianion activation sites of triosephosphate isomerase

Triosephosphate isomerase (TIM) catalyzes the isomerization of dihydroxyacetone phosphate to form d-glyceraldehyde 3-phosphate. The effects of two structural mutations in TIM on the kinetic parameters for catalysis of the reaction of the truncated substrate glycolaldehyde (GA) and the activation of this reaction by phosphite dianion are reported. The P168A mutation results in similar 50- and 80-fold decreases in (kcat/Km)E and (kcat/Km)E·HPi, respectively, for deprotonation of GA catalyzed by free TIM and by the TIM·HPO(3)(2-) complex. The mutation has little effect on the observed and intrinsic phosphite dianion binding energy or the magnitude of phosphite dianion activation of TIM for catalysis of deprotonation of GA. A loop 7 replacement mutant (L7RM) of TIM from chicken muscle was prepared by substitution of the archaeal sequence 208-TGAG with 208-YGGS. L7RM exhibits a 25-fold decrease in (kcat/Km)E and a larger 170-fold decrease in (kcat/Km)E·HPi for reactions of GA. The mutation has little effect on the observed and intrinsic phosphodianion binding energy and only a modest effect on phosphite dianion activation of TIM. The observation that both the P168A and loop 7 replacement mutations affect mainly the kinetic parameters for TIM-catalyzed deprotonation but result in much smaller changes in the parameters for enzyme activation by phosphite dianion provides support for the conclusion that catalysis of proton transfer and dianion activation of TIM take place at separate, weakly interacting, sites in the protein catalyst.

Magnitude and origin of the enhanced basicity of the catalytic glutamate of triosephosphate isomerase

Glu-167 of triosephosphate isomerase from Trypanosoma brucei brucei (TbbTIM) acts as the base to deprotonate substrate to form an enediolate phosphate trianion intermediate. We report that there is a large ~6 pK unit increase in the basicity of the carboxylate side chain of Glu-167 upon binding of the inhibitor phosphoglycolate trianion (I(3-)), an analog of the enediolate phosphate intermediate, from pKEH ≈ 4 for the protonated free enzyme EH to pK(EHI) ≈ 10 for the protonated enzyme-inhibitor complex EH•I(3-). We propose that there is a similar increase in the basicity of this side chain when the physiological substrates are deprotonated by TbbTIM to form an enediolate phosphate trianion intermediate and that it makes an important contribution to the enzymatic rate acceleration. The affinity of wildtype TbbTIM for I(3-) increases 20,000-fold upon decreasing the pH from 9.3 to 4.9, because TbbTIM exists mainly in the basic form E over this pH range, while the inhibitor binds specifically to the rare protonated enzyme EH. This reflects the large increase in the basicity of the carboxylate side chain of Glu-167 upon binding of I(3-) to EH to give EH•I(3-). The I172A mutation at TbbTIM results in an ~100-fold decrease in the affinity of TbbTIM for I(3-) at pH < 6 and an ~2 pK unit decrease in the basicity of the carboxylate side chain of Glu-167 at the EH•I(3-) complex, to pK(EHI) = 7.7. Therefore, the hydrophobic side chain of Ile-172 plays a critical role in effecting the large increase in the basicity of the catalytic base upon the binding of substrate and/or inhibitors.

Specificity in transition state binding: the Pauling model revisited

Linus Pauling proposed that the large rate accelerations for enzymes are caused by the high specificity of the protein catalyst for binding the reaction transition state. The observation that stable analogues of the transition states for enzymatic reactions often act as tight-binding inhibitors provided early support for this simple and elegant proposal. We review experimental results that support the proposal that Pauling's model provides a satisfactory explanation for the rate accelerations for many heterolytic enzymatic reactions through high-energy reaction intermediates, such as proton transfer and decarboxylation. Specificity in transition state binding is obtained when the total intrinsic binding energy of the substrate is significantly larger than the binding energy observed at the Michaelis complex. The results of recent studies that aimed to characterize the specificity in binding of the enolate oxygen at the transition state for the 1,3-isomerization reaction catalyzed by ketosteroid isomerase are reviewed. Interactions between pig heart succinyl-coenzyme A:3-oxoacid coenzyme A transferase (SCOT) and the nonreacting portions of coenzyme A (CoA) are responsible for a rate increase of 3 × 10(12)-fold, which is close to the estimated total 5 × 10(13)-fold enzymatic rate acceleration. Studies that partition the interactions between SCOT and CoA into their contributing parts are reviewed. Interactions of the protein with the substrate phosphodianion group provide an ~12 kcal/mol stabilization of the transition state for the reactions catalyzed by triosephosphate isomerase, orotidine 5'-monophosphate decarboxylase, and α-glycerol phosphate dehydrogenase. The interactions of these enzymes with the substrate piece phosphite dianion provide a 6-8 kcal/mol stabilization of the transition state for reaction of the appropriate truncated substrate. Enzyme activation by phosphite dianion reflects the higher dianion affinity for binding to the enzyme-transition state complex compared with that of the free enzyme. Evidence is presented that supports a model in which the binding energy of the phosphite dianion piece, or the phosphodianion group of the whole substrate, is utilized to drive an enzyme conformational change from an inactive open form E(O) to an active closed form E(C), by closure of a phosphodianion gripper loop. Members of the enolase and haloalkanoic acid dehalogenase superfamilies use variable capping domains to interact with nonreacting portions of the substrate and sequester the substrate from interaction with bulk solvent. Interactions of this capping domain with the phenyl group of mandelate have been shown to activate mandelate racemase for catalysis of deprotonation of α-carbonyl carbon. We propose that an important function of these capping domains is to utilize the binding interactions with nonreacting portions of the substrate to activate the enzyme for catalysis.

Long-range interactions in the α subunit of tryptophan synthase help to coordinate ligand binding, catalysis, and substrate channeling

The α-subunit of tryptophan synthase (αTS) catalyzes the conversion of indole-3-glycerol phosphate to d-glyceraldehyde-3-phosphate and indole. We propose that allosteric networks intrinsic to αTS are modulated by the binding of the β-subunit to regulate αTS function. Understanding these long-range amino acid networks in αTS thus gives insight into the coordination of the two active sites within TS. In this study, we have used Ala residues as probes for structural and dynamic changes of αTS throughout its catalytic cycle, in the absence of the β-subunit. Projection analysis of the chemical shift changes by site-specific amino acid substitutions and ligand titrations indicates that αTS has three important conformational states: ligand-free, glyceraldehyde-3-phosphate-bound(like), and the active states. The amino acid networks within these conformations are different, as suggested by chemical shift correlation analysis. In particular, there are long-range connections, only in the active state, between Ala47, which reports on structural and dynamic changes associated with the general acid/base Glu49, and residues within the β2α2 loop, which contains the catalytically important Asp60 residue. These long-range interactions are likely important for coordinating chemical catalysis. In the free state, but not in the active state, there are connections between the β2α2 and β6α6 loops that likely help to coordinate substrate binding. Changes in the allosteric networks are also accompanied by protein dynamic changes. During catalytic turnover, the protein becomes more rigid on the millisecond timescale and the active-site dynamics are driven to a faster nanosecond timescale.

The Ω-loop lid domain of phosphoenolpyruvate carboxykinase is essential for catalytic function

Phosphoenolpyruvate carboxykinase (PEPCK) is an essential metabolic enzyme operating in the gluconeogenesis and glyceroneogenesis pathways. Recent studies have demonstrated that the enzyme contains a mobile active site lid domain that undergoes a transition between an open, disorded conformation and a closed, ordered conformation as the enzyme progresses through the catalytic cycle. The understanding of how this mobile domain functions in catalysis is incomplete. Previous studies showed that the closure of the lid domain stabilizes the reaction intermediate and protects the reactive intermediate from spurious protonation and thus contributes to the fidelity of the enzyme. To more fully investigate the roles of the lid domain in PEPCK function, we introduced three mutations that replaced the 11-residue lid domain with one, two, and three glycine residues. Kinetic analysis of the mutant enzymes demonstrates that none of the enzyme constructs exhibit any measurable kinetic activity, resulting in a decrease in the catalytic parameters of at least 10(6). Structural characterization of the mutants in complexes representing the catalytic cycle suggests that the inactivity is due to a role for the lid domain in the formation of the fully closed state of the enzyme that is required for catalytic function. In the absence of the lid domain, the enzyme is unable to achieve the fully closed state and is rendered inactive despite possessing all of the residues and substrates required for catalytic function. This work demonstrates how enzyme catalytic function can be abolished through the alteration of conformational equilibria despite all the elements required for chemical conversion of substrates to products remaining intact.

Mechanism for activation of triosephosphate isomerase by phosphite dianion: the role of a hydrophobic clamp

The role of the hydrophobic side chains of Ile-172 and Leu-232 in catalysis of the reversible isomerization of R-glyceraldehyde 3-phosphate (GAP) to dihydroxyacetone phosphate (DHAP) by triosephosphate isomerase (TIM) from Trypanosoma brucei brucei (Tbb) has been investigated. The I172A and L232A mutations result in 100- and 6-fold decreases in k(cat)/K(m) for the isomerization reaction, respectively. The effect of the mutations on the product distributions for the catalyzed reactions of GAP and of [1-(13)C]-glycolaldehyde ([1-(13)C]-GA) in D(2)O is reported. The 40% yield of DHAP from wild-type Tbb TIM-catalyzed isomerization of GAP with intramolecular transfer of hydrogen is found to decrease to 13% and to 4%, respectively, for the reactions catalyzed by the I172A and L232A mutants. Likewise, the 13% yield of [2-(13)C]-GA from isomerization of [1-(13)C]-GA in D(2)O is found to decrease to 2% and to 1%, respectively, for the reactions catalyzed by the I172A and L232A mutants. The decrease in the yield of the product of intramolecular transfer of hydrogen is consistent with a repositioning of groups at the active site that favors transfer of the substrate-derived hydrogen to the protein or the oxygen anion of the bound intermediate. The I172A and L232A mutations result in (a) a >10-fold decrease (I172A) and a 17-fold increase (L232A) in the second-order rate constant for the TIM-catalyzed reaction of [1-(13)C]-GA in D(2)O, (b) a 170-fold decrease (I172A) and 25-fold increase (L232A) in the third-order rate constant for phosphite dianion (HPO(3)(2-)) activation of the TIM-catalyzed reaction of GA in D(2)O, and (c) a 1.5-fold decrease (I172A) and a larger 16-fold decrease (L232A) in K(d) for activation of TIM by HPO(3)(2-) in D(2)O. The effects of the I172A mutation on the kinetic parameters for the wild-type TIM-catalyzed reactions of the whole substrate and substrate pieces are consistent with a decrease in the basicity of the carboxylate side chain of Glu-167 for the mutant enzyme. The data provide striking evidence that the L232A mutation leads to a ca. 1.7 kcal/mol stabilization of a catalytically active loop-closed form of TIM (E(C)) relative to an inactive open form (E(O)).

A Benchmark Test Suite for Proton Transfer Energies and its Use to Test Electronic Structure Model Chemistries

We present benchmark calculations of nine selected points on potential energy surfaces describing proton transfer process in three model systems, H(5)O(2) (+), CH(3)OH…H(+)…OH(2), and CH(3)COOH…OH(2). The calculated relative energies of these geometries are compared to those calculated by various wave function and density functional methods, including the polarized molecular orbital (PMO) model recently developed in our research group and other semiempirical molecular orbital methods. We found that the SCC-DFTB and PMO methods (the latter available so far only for molecules consisting of only O and H and therefore only for the first of the three model systems) give results that are, on average, within 2 kcal/mol of the benchmark results. Other semiempirical molecular orbital methods have mean unsigned errors (MUEs) of 3 to 8 kcal/mol, local density functionals have MUEs in the range 0.7 to 3.7 kcal/mol, and hybrid density functionals have MUEs of only 0.3 to 1.0 kcal/mol, with the best density functional performance obtained by hybrid meta-GGAs, especially M06 and PW6B95.

A paradigm for enzyme-catalyzed proton transfer at carbon: triosephosphate isomerase

Triosephosphate isomerase (TIM) catalyzes the stereospecific 1,2-proton shift at dihydroxyacetone phosphate (DHAP) to give (R)-glyceraldehyde 3-phosphate through a pair of isomeric enzyme-bound cis-enediolate phosphate intermediates. The chemical transformations that occur at the active site of TIM were well understood by the early 1990s. The mechanism for enzyme-catalyzed isomerization is similar to that for the nonenzymatic reaction in water, but the origin of the catalytic rate acceleration is not understood. We review the results of experimental work that show that a substantial fraction of the large 12 kcal/mol intrinsic binding energy of the nonreacting phosphodianion fragment of TIM is utilized to activate the active site side chains for catalysis of proton transfer. Evidence is presented that this activation is due to a phosphodianion-driven conformational change, the most dramatic feature of which is closure of loop 6 over the dianion. The kinetic data are interpreted within the framework of a model in which activation is due to the stabilization by the phosphodianion of a rare, desolvated, loop-closed form of TIM. The dianion binding energy is proposed to drive the otherwise thermodynamically unfavorable desolvation of the solvent-exposed active site. This reduces the effective local dielectric constant of the active site, to enhance stabilizing electrostatic interactions between polar groups and the anionic transition state, and increases the basicity of the carboxylate side chain of Glu-165 that functions to deprotonate the bound carbon acid substrate. A rebuttal is presented to the recent proposal [Samanta, M., Murthy, M. R. N., Balaram, H., and Balaram, P. (2011) ChemBioChem 12, 1886-1895] that the cationic side chain of K12 functions as an active site electrophile to protonate the carbonyl oxygen of DHAP.