Wildtype and engineered monomeric triosephosphate isomerase from Trypanosoma brucei: partitioning of reaction intermediates in D2O and activation by phosphite dianion

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

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.

Q-RepEx: A Python pipeline to increase the sampling of empirical valence bond simulations

The exploration of chemical systems occurs on complex energy landscapes. Comprehensively sampling rugged energy landscapes with many local minima is a common problem for molecular dynamics simulations. These multiple local minima trap the dynamic system, preventing efficient sampling. This is a particular challenge for large biochemical systems with many degrees of freedom. Replica exchange molecular dynamics (REMD) is an approach that accelerates the exploration of the conformational space of a system, and thus can be used to enhance the sampling of complex biomolecular processes. In parallel, the empirical valence bond (EVB) approach is a powerful approach for modeling chemical reactivity in biomolecular systems. Here, we present an open-source Python-based tool that interfaces with the Q simulation package, and increases the sampling efficiency of the EVB free energy perturbation/umbrella sampling approach by means of REMD. This approach, Q-RepEx, both decreases the computational cost of the associated REMD-EVB simulations, and opens the door to more efficient studies of biochemical reactivity in systems with significant conformational fluctuations along the chemical reaction coordinate.

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.

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.

Evaluating the catalytic importance of a conserved Glu97 residue in triosephosphate isomerase

Investigating enzyme activity is central to our understanding of biological function, and the design of biocatalysts continues to find applications in synthesis. While a role for active site residues can be proposed based on structure and mechanism, our understanding of the catalytic importance for residues surrounding the active site is less well understood. In triosephosphate isomerase (TIM), Glu97 is situated adjacent to the active site and is found in essentially all sequences. Prior studies reported mutation of Glu97 to Asp and Gln in TIM from Plasmodium falciparum (PfTIM) led to a 100- and 4000-fold decrease in activity, respectively, while the E97D mutation in TIM from Gallus gallus (cTIM) had no effect on activity. To investigate further the question of how mutations in essentially superimposable structures give different effects, we mutated E97 in TIM from Trypanosoma brucei brucei (TbbTIM), Saccharomyces cerevisiae (yTIM), and human (hTIM). The E97D, E97A, and E97Q mutations led to a ∼three-tenfold decrease in activity, a modest effect compared to the 10²-10³-fold effect in PfTIM. CD and fluorescence studies showed the overall structures for the mutants were essentially unchanged. Structural analysis shows that several residues surrounding E97 differ between PfTIM and TIM from the other organisms, and rearrangements or mispositioning of residues in PfTIM may lead to the different rate effects. The results illustrate the interplay of active site and surrounding residues in affecting catalysis and highlight that understanding of the role of residues surrounding the active site may aid in the incorporation of favorable or avoidance of unfavorable interactions when designing 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.

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.

Evaluation of Antigens for Development of a Serological Test for Human African Trypanosomiasis

Background

Control and elimination of human African trypanosomiasis (HAT) can be accelerated through the use of diagnostic tests that are more accurate and easier to deploy. The goal of this work was to evaluate the immuno-reactivity of antigens and identify candidates to be considered for development of a simple serological test for the detection of Trypanosoma brucei gambiense or T. b. rhodesiense infections, ideally both.

Methodology/Principal Findings

The reactivity of 35 antigens was independently evaluated by slot blot and ELISA against sera from both T. b. gambiense and T. b. rhodesiense infected patients and controls. The antigens that were most reactive by both tests to T. b. gambiense sera were the membrane proteins VSG LiTat 1.3, VSG LiTat 1.5 and ISG64. Reactivity to T. b. rhodesiense sera was highest with VSG LiTat 1.3, VSG LiTat 1.5 and SRA, although much lower than with T. b. gambiense samples. The reactivity of all possible combinations of antigens was also calculated. When the slot blot results of 2 antigens were paired, a VSG LiTat 1.3- ISG75 combination performed best on T. b. gambiense sera, while a VSG LiTat 1.3-VSG LiTat 1.5 combination was the most reactive using ELISA. A combination of SRA and either VSG LiTat 1.3 or VSG LiTat 1.5 had the highest reactivity on T. b. rhodesiense sera according to slot blot, while in ELISA, pairing SRA with either GM6 or VSG LiTat 1.3 yielded the best results.

Conclusions

This study identified antigens that were highly reactive to T. b. gambiense sera, which could be considered for developing a serological test for gambiense HAT, either individually or in combination. Antigens with potential for inclusion in a test for T. b. rhodesiense HAT were also identified, but because their reactivity was comparatively lower, a search for additional antigens would be required before developing a test for this form of the disease.

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.

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.

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.

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.

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

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

The L232A mutation in triosephosphate isomerase (TIM) from Trypanosoma brucei brucei results in a small 6-fold decrease in k(cat)/K(m) for the reversible enzyme-catalyzed isomerization of glyceraldehyde 3-phosphate to give dihydroxyacetone phosphate. In contrast, this mutation leads to a 17-fold increase in the second-order rate constant for the TIM-catalyzed proton transfer reaction of the truncated substrate piece [1-(13)C]glycolaldehyde ([1-(13)C]-GA) in D(2)O, a 25-fold increase in the third-order rate constant for the reaction of the substrate pieces GA and phosphite dianion (HPO(3)(2-)), and a 16-fold decrease in K(d) for binding of HPO(3)(2-) to the free enzyme. Most significantly, the mutation also results in an 11-fold decrease in the extent of activation of the enzyme toward turnover of GA by bound HPO(3)(2-). 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)). We propose that this is due to the relief, in L232A mutant TIM, of unfavorable steric interactions between the bulky hydrophobic side chain of Leu-232 and the basic carboxylate side chain of Glu-167, the catalytic base, which destabilize E(c) relative to E(o).