Search

Ribosomes bind to many metabolic enzymes and change their activity. A general mechanism for ribosome-mediated amplification of metabolic enzyme activity, RAMBO, was formulated and elucidated for the glycolytic enzyme triosephosphate isomerase, TPI. The RAMBO effect results from a ribosome-dependent electric field-substrate dipole interaction energy that can increase or decrease the ground state of the reactant and product to regulate catalytic rates. NMR spectroscopy was used to determine the interaction surface of TPI binding to ribosomes and to measure the corresponding kinetic rates in the absence and presence of intact ribosome particles. Chemical cross-linking and mass spectrometry revealed potential ribosomal protein binding partners of TPI. Structural results and related changes in TPI energetics and activity show that the interaction between TPI and ribosomal protein L11 mediate the RAMBO effect.

Triosephosphate isomerase (TPI) performs the 5th step in glycolysis, operates near the limit of diffusion, and is involved in "moonlighting" functions. Its dimer was found singly phosphorylated at Ser20 (pSer20) in human cells, with this post-translational modification (PTM) showing context-dependent stoichiometry and loss under oxidative stress. We generated synthetic pSer20 proteoforms using cell-free protein synthesis that showed enhanced TPI activity by 4-fold relative to unmodified TPI. Molecular dynamics simulations show that the phosphorylation enables a channel to form that shuttles substrate into the active site. Refolding, kinetic, and crystallographic analyses of point mutants including S20E/G/Q indicate that hetero-dimerization and subunit asymmetry are key features of TPI. Moreover, characterization of an endogenous human TPI tetramer also implicates tetramerization in enzymatic regulation. S20 is highly conserved across eukaryotic TPI, yet most prokaryotes contain E/D at this site, suggesting that phosphorylation of human TPI evolved a new switch to optionally boost an already fast enzyme. Overall, complete characterization of TPI shows how endogenous proteoform discovery can prioritize functional versus bystander PTMs.

Homodimeric triosephosphate isomerases (TIMs) from Trypanosoma cruzi (TcTIM) and Trypanosoma brucei (TbTIM) have markedly similar amino-acid sequences and three-dimensional structures. However, several of their biophysical parameters, such as their susceptibility to sulfhydryl agents and their reactivation speed after being denatured, have significant differences. The causes of these differences were explored with microsecond-scale molecular dynamics (MD) simulations of three different TIM proteins: TcTIM, TbTIM, and a chimeric protein, Mut1. We examined their electrostatic interactions and explored the impact of simulation length on them. The same salt bridge between catalytic residues Lys 14 and Glu 98 was observed in all three proteins, but key differences were found in other interactions that the catalytic amino acids form. In particular, a cation-π interaction between catalytic amino acids Lys 14 and His 96 and both a salt bridge and a hydrogen bond between catalytic Glu 168 and residue Arg 100 were only observed in TcTIM. Furthermore, although TcTIM forms less hydrogen bonds than TbTIM and Mut1, its hydrogen bond network spans almost the entire protein, connecting the residues in both monomers. This work provides new insight into the mechanisms that give rise to the different behavior of these proteins. The results also show the importance of long simulations.

Every reader knows that an enzyme accelerates a reaction by reducing the activation-energy barrier. However, understanding how this is achieved by the structure of the enzyme and its interactions with stable complexes and transition states and, then, using this to (re)design enzymes to catalyze novel reactions remain the "holy grail" of mechanistic enzymology. The necessary foundation is the free-energy profile that specifies the energies of the bound substate, product, and intervening intermediates as well as the transition states by which they are interconverted. When this free-energy profile is compared to that for the uncatalyzed reaction, strategies for establishing and enhancing catalysis can be identified. This Perspective reminds readers that the first free-energy profile determined for an enzyme-catalyzed reaction, that for triosephosphate isomerase, was published in Biochemistry in 1976 by Jeremy R. Knowles, W. John Albery, and co-workers. They used the profile to propose three steps of increasing "subtlety" that can be influenced by evolutionary pressure to increase the flux through the reaction coordinate: (1) "uniform binding" of the substrate, product, and intermediates; (2) "differential binding" of complexes so that these are isoenergetic (to minimize the energy of the intervening transition states); and (3) "catalysis of an elementary step" in which the transition state for the kinetically significant chemical step is stabilized so that flux can be determined by the rate of substrate binding or product dissociation. These papers continue to guide mechanistic studies of enzyme-catalyzed reactions and provide principles for the (re)design of novel enzymes.

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

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

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.

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.

Simulation methods are valuable for elucidating protein conformational transitions between functionally diverse states given that transition pathways are difficult to capture experimentally. Nonetheless, specific computational algorithms are required because of the high free energy barriers between these different protein conformational states. Adaptively biased path optimization (ABPO) is an unrestrained, transition-path optimization method that works in a reduced-variable space to construct an adaptive biasing potential to aid convergence. ABPO was previously applied using a coarse-grained Go̅-model to study conformational activation of Lyn, a Src family tyrosine kinase. How effectively ABPO can be applied at the higher resolution of an all-atom model to explore protein conformational transitions is not yet known. Here, we report the all-atom conformational transition paths of three protein systems constructed using the ABPO methodology. Two systems, triose phosphate isomerase and dihydrofolate reductase, undergo local flipping of a short loop that promotes ligand binding. The third system, estrogen receptor α ligand binding domain, has a helix that adopts different conformations when the protein is bound to an agonist or an antagonist. For each protein, distance-based or torsion-angle reduced variables were identified from unbiased trajectories. ABPO was computed in this reduced variable space to obtain the transition path between the two states. The all-atom ABPO is shown to successfully converge an optimal transition path for each of the three systems.

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 3 - 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 3 - at high pH, in which cases the pH profiles define the p K a for deprotonation of EH·I 3 - . The linear free energy correlation, with slope of 0.73 ( r 2 = 0.96), between k cat / K m for TIM-catalyzed isomerization and the disassociation constant of PGA trianion for TIM shows that EH·I 3 - and the transition state are stabilized by similar interactions with the protein catalyst. Values of p K a = 10-10.5 were estimated for deprotonation of EH·I 3 - for wildtype TIM. This p K 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 p K a for EH·I 3 - . 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.

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

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

The main problem for enzymes from psychrophilic species, which need to work near the freezing point of liquid water, is the exponential decay of reaction rates as the temperature is decreased. Cold-adapted enzymes have solved this problem by shifting the activation enthalpy-entropy balance for the catalyzed reaction compared to those of their mesophilic orthologs. To understand the structural basis of this universal feature, it is necessary to examine pairs of such orthologous enzymes, with known three-dimensional structures, at the microscopic level. Here, we use molecular dynamics free energy calculations in combination with the empirical valence bond method to evaluate the temperature dependence of the activation free energy for differently adapted triosephosphate isomerases. The results show that the enzyme from the psychrophilic bacterium Vibrio marinus indeed displays the characteristic shift in enthalpy-entropy balance, compared to that of the yeast ortholog. The origin of this effect is found to be located in a few surface-exposed protein loops that show differential mobilities in the two enzymes. Key mutations render these loops more mobile in the cold-adapted triosephosphate isomerase, which explains both the reduced activation enthalpy contribution from the protein surface and the lower thermostability.

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.

Accurate sampling of conformational space and, in particular, the transitions between functional substates has been a challenge in molecular dynamic (MD) simulations of large biomolecular systems. We developed an Elastic Network Model (ENM)-based computational method, ClustENM, for sampling large conformational changes of biomolecules with various sizes and oligomerization states. ClustENM is an iterative method that combines ENM with energy minimization and clustering steps. It is an unbiased technique, which requires only an initial structure as input, and no information about the target conformation. To test the performance of ClustENM, we applied it to six biomolecular systems: adenylate kinase (AK), calmodulin, p38 MAP kinase, HIV-1 reverse transcriptase (RT), triosephosphate isomerase (TIM), and the 70S ribosomal complex. The generated ensembles of conformers determined at atomic resolution show good agreement with experimental data (979 structures resolved by X-ray and/or NMR) and encompass the subspaces covered in independent MD simulations for TIM, p38, and RT. ClustENM emerges as a computationally efficient tool for characterizing the conformational space of large systems at atomic detail, in addition to generating a representative ensemble of conformers that can be advantageously used in simulating substrate/ligand-binding events.

Kinetic parameters are reported for the reactions of whole substrates (kcat/Km, M(-1) s(-1)) (R)-glyceraldehyde 3-phosphate (GAP) and dihydroxyacetone phosphate (DHAP) and for the substrate pieces [(kcat/Km)E·HPi/Kd, M(-2) s(-1)] glycolaldehyde (GA) and phosphite dianion (HPi) 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.

Analysis of protein complexes by ion mobility-mass spectrometry is a valuable method for the rapid assessment of complex composition, binding stoichiometries, and structures. However, capturing labile, unknown protein assemblies directly from cells remains a challenge for the technology. Furthermore, ion mobility-mass spectrometry measurements of complexes, subcomplexes, and subunits are necessary to build complete models of intact assemblies, and such data can be difficult to acquire in a comprehensive fashion. Here, we present the use of novel mass spectrometry cleavable cross-linkers and tags to stabilize intact protein complexes for ion mobility-mass spectrometry. Our data reveal that tags and linkers bearing permanent charges are superior stabilizers relative to neutral cross-linkers, especially in the context of retaining compact forms of the assembly under a wide array of activating conditions. In addition, when cross-linked protein complexes are collisionally activated in the gas phase, a larger proportion of the product ions produced are often more compact and reflect native protein subcomplexes when compared with unmodified complexes activated in the same fashion, greatly enabling applications in structural biology.

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(HPi)K(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.

Quinolinate synthase (NadA) is an Fe4S4 cluster-containing dehydrating enzyme involved in the synthesis of quinolinic acid (QA), the universal precursor of the essential coenzyme nicotinamide adenine dinucleotide. The reaction catalyzed by NadA is not well understood, and two mechanisms have been proposed in the literature that differ in the nature of the molecule (DHAP or G-3P) that condenses with iminoaspartate (IA) to form QA. In this article, using biochemical approaches, we demonstrate that DHAP is the triose that condenses with IA to form QA. The capacity of NadA to use G-3P is due to its previously unknown triose phosphate isomerase activity.

A Monte Carlo linear inverse-kinematics method for the simulation of protein chains with fixed ends is introduced. It includes backbone bond-angle bending and simultaneous loop and ring closure to allow full proline ring flexibility. An obstacle to linear null-space methods is the eventual drift of the end group. Maintenance of the end group at its initial position by occasional reset is performed in a way that is consistent with the overall methodology and minimally disruptive to the current conformation. The implementation permitted multiple rigid regions within the chain, enabling the simulation of domain movements where domains are rigid bodies connected by flexible interdomain regions. The method was tested on polyalanine, polyglycine, loop 6 of triosephosphate isomerase, and glutamine binding protein. Simulations of glutamine binding protein, where only 11 of the 226 residues at the interdomain bending regions were flexible, accurately reproduced the experimentally determined domain movement.

The fully analytic first and second derivatives of the energy in the frozen domain formulation of the fragment molecular orbital (FMO) were developed and applied to locate transition states and determine vibrational contributions to free energies. The development is focused on the frozen domain with dimers (FDD) model. The intrinsic reaction coordinate method was interfaced with FMO. Simulations of IR and Raman spectra were enabled using FMO/FDD by developing the calculation of intensities. The accuracy is evaluated for S(N)2 reactions in explicit solvent, and for the free binding energies of a protein-ligand complex of the Trp cage protein (PDB: 1L2Y ). FMO/FDD is applied to study the keto-enol tautomeric reaction of phosphoglycolohydroxamic acid and the triosephosphate isomerase (PDB: 7TIM ), and the role of amino acid residue fragments in the reaction is discussed.

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.

Values of (k(cat)/K(m))GAP for triosephosphate isomerase-catalyzed reactions of (R)-glyceraldehyde 3-phosphate and k(cat)/K(HPi)K(GA) for reactions of the substrate pieces glycolaldehyde and HPO3(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.

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.

One-step purification of a multi-enzyme complex was developed based on a mixture of cell extracts containing three dockerin-containing enzymes and one family 3 cellulose-binding module (CBM3)-containing scaffoldin through high-affinity adsorption on low-cost solid regenerated amorphous cellulose (RAC). The three-enzyme complex, called synthetic metabolon, was self-assembled through the high-affinity interaction between the dockerin in each enzyme and three cohesins in the synthetic scaffoldin. The metabolons were either immobilized on the external surface of RAC or free when the scaffoldin contained an intein between the CBM3 and three cohesins. The immobilized and free metabolons containing triosephosphate isomerase, aldolase, and fructose 1,6-biphosphatase exhibited initial reaction rates 48 and 38 times, respectively, that of the non-complexed three-enzyme mixture at the same enzyme loading. Such reaction rate enhancements indicated strong substrate channeling among synthetic metabolons due to the close spatial organization among cascade enzymes. These results suggested that the construction of synthetic metabolons by using cohesins, dockerins, and cellulose-binding modules from cellulosomes not only decreased protein purification labor and cost for in vitro synthetic biology projects but also accelerated reaction rates by 1 order of magnitude compared to non-complexed enzymes. Synthetic metabolons would be an important biocatalytic module for in vitro and in vivo synthetic biology projects.

Recent molecular dynamics simulations have suggested important roles for nanoscale dewetting in the stability, function, and folding dynamics of proteins. Using a synergistic simulation-experimental approach on the αTS TIM barrel protein, we validated this hypothesis by revealing the occurrence of drying inside hydrophobic amino acid clusters and its manifestation in experimental measures of protein stability and structure. Cavities created within three clusters of branched aliphatic amino acids [isoleucine, leucine, and valine (ILV) clusters] were found to experience strong water density fluctuations or intermittent dewetting transitions in simulations. Individually substituting 10 residues in the large ILV cluster at the N-terminus with less hydrophobic alanines showed a weakening or diminishing effect on dewetting that depended on the site of the mutation. Our simulations also demonstrated that replacement of buried leucines with isosteric, polar asparagines enhanced the wetting of the N- and C-terminal clusters. The experimental results on the stability, secondary structure, and compactness of the native and intermediate states for the asparagine variants are consistent with the preferential drying of the large N-terminal cluster in the intermediate. By contrast, the region encompassing the small C-terminal cluster experiences only partial drying in the intermediate, and its structure and stability are unaffected by the asparagine substitution. Surprisingly, the structural distortions required to accommodate the replacement of leucine by asparagine in the N-terminal cluster revealed the existence of alternative stable folds in the native basin. This combined simulation-experimental study demonstrates the critical role of drying within hydrophobic ILV clusters in the folding and stability of the αTS TIM barrel.

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

Deamidation of asparagine is the spontaneous degradation of this residue into aspartic acid. The kinetics of this slow reaction is mainly dependent on the nature of the adjacent amino acid that follows asparagine in a peptide or protein primary sequence. In the homodimer triosephosphate isomerase (TPI), there are two main deamidation sites per subunit: Asn15-Gly16 and Asn71-Gly72 for which deamidation dynamics are known to be interrelated. In this study, we investigate the initiation of the deamidation reaction in TPI by means of molecular dynamics. Simulations based on classical AMBER force field are performed in a 60 to 90 ns time scale for six distinct samples. Conformational changes, desolvation effects, and hydrogen bond networks are analyzed to interpret the experimental findings and previous quantum mechanical (QM) results. Results that are based on desolvation analysis clarify the assignments in the literature about the different behaviors of two deamidating sites in TPI. Conformational analysis supports findings suggested by QM studies: the most favorable reaction mechanism is the one that yields to succinimide intermediate via one or two step routes. The mechanism leading to the succinimide intermediate most likely involves the formation of a tetrahedral intermediate that is formed either directly from asparagine or via a side chain tautomer intermediate. In all cases, surrounding water molecules are present to assist the reaction.

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.

Protein-protein interactions are crucial for many biological functions. The redox interactome encompasses numerous weak transient interactions in which thioredoxin plays a central role. Proteomic studies have shown that thioredoxin binds to numerous proteins belonging to various cellular processes, including energy metabolism. Thioredoxin has cross talk with other redox mechanisms involving glutathionylation and has functional overlap with glutaredoxin in deglutathionylation reactions. In this study, we have explored the structural and biochemical interactions of thioredoxin with the glycolytic enzyme, triosephosphate isomerase. Nuclear magnetic resonance chemical shift mapping methods and molecular dynamics-based docking have been applied in deriving a structural model of the thioredoxin-triosephosphate isomerase complex. The spatial proximity of active site cysteine residues of thioredoxin to reactive thiol groups on triosephosphate isomerase provides a direct link to the observed deglutathionylation of cysteine 217 in triosephosphate isomerase, thereby reversing the inhibitory effect of S-glutathionylation of triosephosphate isomerase.

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

Product yields for the reactions of (R)-glyceraldehyde 3-phosphate (GAP) in D2O at pD 7.9 catalyzed by wildtype triosephosphate isomerase from Trypanosoma brucei brucei (Tbb TIM) and a monomeric variant (monoTIM) of this wildtype enzyme were determined by (1)H NMR spectroscopy and were compared with the yields determined in earlier work for the reactions catalyzed by TIM from rabbit and chicken muscle [O'Donoghue, A. C., Amyes, T. L., and Richard, J. P. (2005), Biochemistry 44, 2610 - 2621]. Three products were observed from the reactions catalyzed by TIM: dihydroxyacetone phosphate (DHAP) from isomerization with intramolecular transfer of hydrogen, d-DHAP from isomerization with incorporation of deuterium from D2O into C-1 of DHAP, and d-GAP from incorporation of deuterium from D2O into C-2 of GAP. The yield of DHAP formed by intramolecular transfer of hydrogen decreases from 49% for the muscle enzymes to 40% for wildtype Tbb TIM to 34% for monoTIM. There is no significant difference in the ratio of the yields of d-DHAP and d-GAP for wildtype TIM from muscle sources and Trypanosoma brucei brucei, but partitioning of the enediolate intermediate of the monoTIM reaction to form d-DHAP is less favorable ((k(C1))(D)/(k(C2))(D) = 1.1) than for the wildtype enzyme ((k(C1))(D)/(k(C2))(D) = 1.7). Product yields for the wildtype Tbb TIM and monoTIM-catalyzed reactions of glycolaldehyde labeled with carbon-13 at the carbonyl carbon ([1-(13)C]-GA) at pD 7.0 in the presence of phosphite dianion and in its absence were determined by (1)H NMR spectroscopy [Go, M. K., Amyes, T. L., and Richard, J. P. (2009) Biochemistry 48, 5769-5778]. There is no detectable difference in the yields of the products of wildtype muscle and Tbb TIM-catalyzed reactions of [1-(13)C]-GA in D2O. The kinetic parameters for phosphite dianion activation of the reactions of [1-(13)C]-GA catalyzed by wildtype Tbb TIM are similar to those reported for the enzyme from rabbit muscle [Amyes, T. L. and Richard, J. P. (2007) Biochemistry 46, 5841-5854], but there is no detectable dianion activation of the reaction catalyzed by monoTIM. The engineered disruption of subunit contacts at monoTIM causes movement of the essential side chains of Lys-13 and His-95 away from the catalytic active positions. We suggest that this places an increased demand that the intrinsic binding energy of phosphite dianion be utilized to drive the change in the conformation of monoTIM back to the active structure for wildtype TIM.

The reactions of two bacterial TIM barrel prenyltransferases (PTs), MoeO5 and PcrB, were explored. MoeO5, the enzyme responsible for the first step in moenomycin biosynthesis, catalyzes the transfer of farnesyl to 3-phosphoglyceric acid (3PG) to give a product containing a cis-allylic double bond. We show that this reaction involves isomerization to a nerolidyl pyrophosphate intermediate followed by bond rotation prior to attack by the nucleophile. This mechanism is unprecedented for a prenyltransferase that catalyzes an intermolecular coupling. We also show that PcrB transfers geranyl and geranylgeranyl groups to glycerol-1-phosphate (G1P), making it the first known bacterial enzyme to use G1P as a substrate. Unlike MoeO5, PcrB catalyzes prenyl transfer without isomerization to give products that retain the trans-allylic bond of the prenyl donors. The TIM barrel family of PTs is unique in including enzymes that catalyze prenyl transfer by distinctly different reaction mechanisms.

We report that the K12G mutation in triosephosphate isomerase (TIM) from Saccharomyces cerevisiae results in (1) a approximately 50-fold increase in K(m) for the substrate glyceraldehyde 3-phosphate (GAP) and a 60-fold increase in K(i) for competitive inhibition by the intermediate analogue 2-phosphoglycolate, resulting from the loss of stabilizing ground state interactions between the alkylammonium side chain of Lys-12 and the ligand phosphodianion group; (2) a 12000-fold decrease in k(cat) for isomerization of GAP, suggesting a tightening of interactions between the side chain of Lys-12 and the substrate on proceeding from the Michaelis complex to the transition state; and (3) a 6 x 10(5)-fold decrease in k(cat)/K(m), corresponding to a total 7.8 kcal/mol stabilization of the transition state by the cationic side chain of Lys-12. The yields of the four products of the K12G TIM-catalyzed isomerization of GAP in D(2)O were quantified as dihydroxyacetone phosphate (DHAP) (27%), [1(R)-(2)H]DHAP (23%), [2(R)-(2)H]GAP (31%), and methylglyoxal (18%) from an enzyme-catalyzed elimination reaction. The K12G mutation has only a small effect on the relative yields of the three products of the transfer of a proton to the TIM-bound enediol(ate) intermediate in D(2)O, but it strongly favors catalysis of the elimination reaction to give methylglyoxal. The K12G mutation also results in a >or=14-fold decrease in k(cat)/K(m) for isomerization of bound glycolaldehyde (GA), although the dominant observed product of the mutant enzyme-catalyzed reaction of [1-(13)C]GA in D(2)O is [1-(13)C,2,2-di-(2)H]GA from a nonspecific protein-catalyzed reaction. The observation that the K12G mutation results in a large decrease in k(cat)/K(m) for the reactions of both GAP and the neutral truncated substrate [1-(13)C]GA provides evidence for a stabilizing interaction between the cationic side chain of Lys-12 and the negative charge that develops at the enolate-like oxygen in the transition state for deprotonation of the sugar substrate "piece".

Capsaicin is the major pungent component in hot chili pepper. It is known that capsaicin induces the increase of energy metabolism through adrenaline activation. The aim of this study was to investigate the mechanism underlying the effect of capsaicin on energy metabolism in human intestinal epithelial cells (Caco-2). To clarify the mechanism, we performed proteomics, real time-PCR, ATP measurement, and MTT in this study. Glycolytic enzymes, namely, phosphoglycerate mutase (PGM) and triosephosphate isomerase (TPI), were overexpressed in capsaicin-treated Caco-2 cells. mRNA expression levels of TPI and PGM were also increased in capsaicin-treated Caco-2 cells. Furthermore, intracellular adenosine triphosphate (ATP) content, which is the end product of glycolysis, was increased. To demonstrate the intracellular ATP production and ATP accumulation, the viability of Caco-2 cells incubated with glucose free medium was measured. The viability of capsaicin pretreated cells was higher than that of untreated cells. All these observations strongly indicate that capsaicin increases energy metabolism in human epithelial cells through the activation of glycolytic enzymes.

A systematic study of the convergence of QM/MM results with respect to the chosen size of the QM region is presented for two examples of peptidic systems. For this purpose, we increased the QM region to up to 1637 atoms at the HF/SVP and 383 atoms at the SOS-AO-MP2/6-31G** level. While the convergence behavior is almost independent of the chosen method and basis set, the study clearly shows that for the considered proton-transfer energy the QM/MM treatment leads to a significantly faster convergence than the pure QM treatment. This behavior can be rationalized by the fair description of the surrounding of the active center using MM methods, even though the MM description of the active center is not adequate in our present case. At the same time, the observed convergence is quite insensitive to a variation of charge surroundings in the chosen model peptides. Although the QM/MM results do converge much quicker with the system size than the pure QM ones, the data show that even for the chosen simple model systems about 150-300 QM atoms are needed to achieve accuracies in the order of 10 kJ/mol and about 300-1000 atoms for an accuracy of 2 kJ/mol with respect to a convergence with the QM-region size.

The effect of treatments with four herbicides and a safener on the activity of triosephosphate isomerase (TPI) extracted from shoots of Italian ryegrass was investigated. It was found that atrazine and fluorodifen, herbicides which interfere with photosynthesis, caused a decrease in measured enzyme activity. In addition, the in vitro effect of oxidized glutathione (GSSG), a compound produced in situations of oxidative stress, on TPI activity was investigated. It was shown that GSSG was a strong inhibitor of enzyme activity, at low concentrations in a dose-time-dependent manner. The enzyme extracts were submitted to chromatographic purifications and to two-dimensional electrophoresis. Some spots had molecular masses ranging between 20 and 30 kDa and were characterized and identified by LC-ESI-MS/MS as TPIs. The mass spectrometry also made it possible to identify the presence of cysteine residues that could be subjected to S-glutathionylation, which regulate the enzyme activity.

Product distributions for the reaction of glycolaldehyde labeled with carbon-13 at the carbonyl carbon ([1-(13)C]-GA) catalyzed by triosephosphate isomerase (TIM) in D(2)O at pD 7.0 in the presence of phosphite dianion and in its absence were determined by (1)H NMR spectroscopy. We observe three products for the relatively fast phosphite-activated reaction (Amyes, T. L., and Richard, J. P. (2007) Biochemistry 46, 5841-5854): [2-(13)C]-GA from isomerization with intramolecular transfer of hydrogen (12% of products), [2-(13)C,2-(2)H]-GA from isomerization with incorporation of deuterium from D(2)O at C-2 (64% of products), and [1-(13)C,2-(2)H]-GA from incorporation of deuterium from D(2)O at C-2 (23% of products). The much slower unactivated reaction in the absence of phosphite results in formation of the same three products along with the doubly deuterated product [1-(13)C,2,2-(2)H(2)]-GA. The two isomerization products ([2-(13)C]-GA and [2-(13)C,2-(2)H]-GA) are formed in the same relative yields in both the unactivated and the phosphite-activated reactions. However, the additional [1-(13)C,2-(2)H]-GA and the doubly deuterated [1-(13)C,2,2-(2)H(2)]-GA formed in the unactivated TIM-catalyzed reaction are proposed to result from nonspecific reaction(s) at the protein surface. The data provide evidence that phosphite dianion affects the rate, but not the product distribution, of the TIM-catalyzed reaction of [1-(13)C]-GA at the enzyme active site. They are consistent with the conclusion that both reactions occur at an unstable loop-closed form of TIM and that activation of the isomerization reaction by phosphite dianion results from utilization of the intrinsic binding energy of phosphite dianion to stabilize the active loop-closed enzyme.

The enzyme triosephosphate isomerase (TIM) has been used as a model system for understanding the relationship between protein sequence, structure, and biological function. The sequence of the active site loop (loop 6) in TIM is directly correlated with a conserved motif in loop 7. Replacement of loop 7 of chicken TIM with the corresponding loop 7 sequence from an archaeal homologue caused a 10(2)-fold loss in enzymatic activity, a decrease in substrate binding affinity, and a decrease in thermal stability. Isotope exchange studies performed by one-dimensional (1)H NMR showed that the substrate-derived proton in the enzyme is more susceptible to solvent exchange for DHAP formation in the loop 7 mutant than for WT TIM. TROSY-Hahn Echo and TROSY-selected R(1rho) experiments indicate that upon mutation of loop 7, the chemical exchange rate for active site loop motion is nearly doubled and that the coordinated motion of loop 6 is reduced relative to that of the WT. Temperature dependent NMR experiments show differing activation energies for the N- and C-terminal hinges in this mutant enzyme. Together, these data suggest that interactions between loop 6 and loop 7 are necessary to provide the proper chemical context for the enzymatic reaction to occur and that the interactions play a significant role in modulating the chemical dynamics near the active site.

Gluconeogenesis is blocked in a strain of Escherichia coli that is deficient in triosephosphate isomerase, but it was restored by the insertion of a plasmid coding for an L-glyceraldehyde 3-phosphate reductase (YghZ). This reductase provides a "bypass" that produces dihydroxyacetone phosphate (DHAP) by the consecutive enzyme-catalyzed reduction of L-glyceraldehyde 3-phosphate ( L-GAP) by NADPH to give L-glycerol 3-phosphate and reoxidation by NAD(+) catalyzed by endogenous L-glycerol 3-phosphate dehydrogenase to give DHAP. The origin of cellular L-GAP remains to be determined.

Triosephosphate isomerase (TIM) catalyzes the interconversion of d-glyceraldehyde 3-phosphate and dihydroxyacetone phosphate, an essential step in glycolytic and gluconeogenic metabolism. To uncover promiscuous isomerases embedded within the Escherichia coli genome, we searched for genes capable of restoring growth of a TIM-deficient bacterium under gluconeogenic conditions. Rather than discovering an isomerase, we selected yghZ, a gene encoding a member of the aldo-keto reductase superfamily. Here we show that YghZ catalyzes the stereospecific, NADPH-dependent reduction of l-glyceraldehyde 3-phosphate, the enantiomer of the TIM substrate. This transformation provides an alternate pathway to the formation of dihydroxyacetone phosphate.

Triosephosphate isomerase from Saccharomyces cerevisiae (wt-TIM) is an obligated homodimer. The interface of wt-TIM is formed by 34 residues. In the native dimer, each monomer buries nearly 2600 A(2) of accessible surface area (ASA), and 58.4% of the interface ASA is hydrophobic. We determined the thermodynamic and functional consequences of increasing the hydrophobic character of the wt-TIM interface. Mutations were restricted to a cluster of five nonconserved residues located far from the active site. Two different approaches, in silico design and directed evolution, were employed. In both methodologies, the obtained proteins were soluble, dimeric, and compact. In silico-designed proteins are very stable dimers that bind substrate with a wild-type-like K(m); albeit, they exhibited a very low k cat. Proteins obtained from directed evolution experiments show wild-type-like catalytic activity, while their stability is decreased. Hydrophobic replacements at the interface produced a remarkable shift in the dissociation step. For wt-TIM and for TIMs obtained by directed evolution, dissociation was observed in the first transition, with C(1/2) values ranging from 0.58 to 0.024 M GdnHCl, whereas for TIMs generated by in silico design, dissociation occurred in the last transition, with C(1/2) values ranging form 3.01 to 3.65 M GdnHCl. For the latter mutants, the stabilization of the interface changed the equilibrium transitions to a novel four-state process with two dimeric intermediates. The change in the intermediate nature suggests that the relative stabilities of different folding units are similar so that subtle alterations in their stability produce a total transformation of the folding pathway.

Cysteine 14 is an interface residue that is fundamental for the catalysis and stability of homodimeric triosephosphate isomerase from Trypanosoma brucei (TbTIM). Its side chain is surrounded by a deep pocket of 11 residues that are part of loop 3 of the adjacent monomer. Mutation of this residue to serine (producing single mutant C14S) yields a wild-type-like enzyme that is resistant to the action of sulfhydryl reagents methylmethane thiosulfonate (MMTS) and 5,5-dithiobis(2-nitrobenzoate) (DTNB). This mutant enzyme was a starting point for probing by cysteine scanning the role of four residues of loop 3 in the catalysis and stability of the enzyme. Considering that the conservative substitution of either serine or alanine with cysteine would minimally alter the structure and properties of the environment of the residue in position 14, we made double mutants C14S/A69C, C14S/S71C, C14S/A73C, and C14S/S79C. Three of these double mutants were similar in their kinetic parameters to wild-type TbTIM and the single mutant C14S, but double mutant C14S/A73C showed a greatly reduced k cat. All enzymes had similar CD spectra, but all mutants had thermal stabilities lower than that of wild-type TbTIM. Intrinsic fluorescence was also similar for all enzymes, but the double mutants bound up to 50 times more 1-anilino-8-naphthalene sulfonate (ANS) and were susceptible to digestion with subtilisin. The double mutants were also susceptible to inactivation by sulfhydryl reagents. Double mutant C14S/S79C exhibited the highest sensitivity to MMTS and DTNB, bound a significant amount of ANS, and had the highest sensitivity to subtilisin. Thus, the residues at positions 73 and 79 are critical for the catalysis and stability of TbTIM, respectively.

One of the most common, and yet least well understood, enzymatic transformations is proton abstraction from activated carbon acids such as carbonyls. Understanding the mechanism for these proton abstractions is basic to a good understanding of enzyme function. Significant controversy has arisen over the means by which a given enzyme might facilitate these deprotonations. Creating small molecule mimics of enzymes and physical organic studies that model enzymes are good approaches to probing mechanistic enzymology. This Account details a number of molecular recognition and physical organic studies, both from our laboratory and others, dealing with the elucidation of this quandary. Our analysis launches from an examination of the active sites and proposed mechanism of several enzyme-catalyzed deprotonations of carbon acids. This analysis highlights the geometries of the hydrogen bonds found at the enzyme active sites. We find evidence to support pi-oriented hydrogen bonding, rather than lone pair oriented hydrogen bonding. Our observations prompted us to study the stereochemistry of hydrogen bonding that activates carbonyl alpha-carbons to deprotonation. The results from our own thermodynamic, kinetics, and computational studies, all of which are reviewed herein, suggest that an unanticipated level of intermediate stabilization occurs via an electrophilic interaction through the pi-molecular orbital as opposed to traditional lone pair directed coordination. We do not postulate that hydrogen bonding to pi-systems is intrinsically stronger than to lone pairs, but rather that there is a greater change in bond strength during deprotonation when the hydrogen bonds are oriented at the pi-system. Through these studies, we conclude that many enzymes preferentially activate their carbon acid substrates through an electrophilic coordination directed towards the pi-bond of the carbonyl rather than the conventional lone pair directed model.

During cell volume regulation, intracellular concentration changes occur in both inorganic and organic osmolytes in order to balance the extracellular osmotic stress and maintain cell volume homeostasis. Generally, salt and urea increase the Km's of enzymes and trimethylamine N-oxide (TMAO) counteracts these effects by decreasing Km's. The hypothesis to account for these effects is that urea and salt shift the native state ensemble of the enzyme toward conformers that are substrate-binding incompetent (BI), while TMAO shifts the ensemble toward binding competent (BC) species. Km's are often complex assemblies of rate constants involving several elementary steps in catalysis, so to better understand osmolyte effects we have focused on a single elementary event, substrate binding. We test the conformational shift hypothesis by evaluating the effects of salt, urea, and TMAO on the mechanism of binding glycerol 3-phosphate, a substrate analogue, to yeast triosephosphate isomerase. Temperature-jump kinetic measurements promote a mechanism consistent with osmolyte-induced shifts in the [BI]/[BC] ratio of enzyme conformers. Importantly, salt significantly affects the binding constant through its effect on the activity coefficients of substrate, enzyme, and enzyme-substrate complex, and it is likely that TMAO and urea affect activity coefficients as well. Results indicate that the conformational shift hypothesis alone does not account for the effects of osmolytes on Km's.

More than 80% of the rate acceleration for enzymatic catalysis of the aldose-ketose isomerization of (R)-glyceraldehyde 3-phosphate (GAP) by triosephosphate isomerase (TIM) can be attributed to the phosphodianion group of GAP [Amyes, T. L., O'Donoghue, A. C., and Richard, J. P. (2001) J. Am. Chem. Soc. 123, 11325-11326]. We examine here the necessity of the covalent connection between the phosphodianion and triose sugar portions of the substrate by "carving up" GAP into the minimal neutral two-carbon sugar glycolaldehyde and phosphite dianion pieces. This "two-part substrate" preserves both the alpha-hydroxycarbonyl and oxydianion portions of GAP. TIM catalyzes proton transfer from glycolaldehyde in D2O, resulting in deuterium incorporation that can be monitored by 1H NMR spectroscopy, with kcat/Km = 0.26 M-1 s-1. Exogenous phosphite dianion results in a very large increase in the observed second-order rate constant (kcat/Km)obsd for turnover of glycolaldehyde, and the dependence of (kcat/Km)obsd on [HPO32-] exhibits saturation. The data give kcat/Km = 185 M-1 s-1 for turnover of glycolaldehyde by TIM that is saturated with phosphite dianion so that the separate binding of phosphite dianion to TIM results in a 700-fold acceleration of proton transfer from carbon. The binding of phosphite dianion to the free enzyme (Kd = 38 mM) is 700-fold weaker than its binding to the fleeting complex of TIM with the altered substrate in the transition state (Kd = 53 muM); the total intrinsic binding energy of phosphite dianion in the transition state is 5.8 kcal/mol. We propose a physical model for catalysis by TIM in which the intrinsic binding energy of the substrate phosphodianion group is utilized to drive closing of the "mobile loop" and a protein conformational change that leads to formation of an active site environment that is optimally organized for stabilization of the transition state for proton transfer from alpha-carbonyl carbon.

The importance of the fully conserved active site proline, Pro168, for the reaction mechanism of triosephosphate isomerase (TIM) has been investigated by studying the enzymatic and crystallographic properties of the P168A variant of trypanosomal TIM. In TIM, Pro168 follows the key catalytic residue Glu167, situated at the beginning of the flexible active site loop (loop 6). Turnover numbers of the P168A variant for its substrates are reduced approximately 50-fold, whereas the Km values are approximately 2 times lower. The affinity of the P168A variant for the transition state analogue 2-phosphoglycolate (2PG) is reduced 5-fold. The crystal structures of unliganded and liganded (2PG) P168A show that the phosphate moiety of 2PG is bound similarly as in wild-type TIM, whereas the interactions of the carboxylic acid moiety with the side chain of the catalytic Glu167 differ. The unique properties of the proline side chain at position 168 are required to transmit ligand binding to the conformational change of Glu167: the side chain of Glu167 flips from the inactive swung-out to the active swung-in conformation on ligand binding in wild-type TIM, whereas in the mutant this conformational change does not occur. Further structural comparisons show that in the wild-type enzyme the concerted movement of loop 6 and loop 7 from unliganded-open to liganded-closed appears to be facilitated by the interactions of the phosphate moiety with loop 7. Apparently, the rotation of 90 degrees of the Gly211-Gly212 peptide plane of loop 7 plays a key role in this concerted movement.

Solution NMR spin relaxation experiments and classical MD simulations are used to study the dynamics of triosephosphate isomerase (TIM) in complex with glycerol 3-phosphate (G3P). Three regions in TIM exhibit conformational transitions on the micros-ms time scale as detected by chemical exchange broadening effects in NMR spectroscopy: residue Lys 84 on helix C, located at the dimeric interface; active site loop 6; and helix G. The results indicate that the conformational exchange process affecting the residues of loop 6 is the correlated opening and closing of the loop. Distinct processes are responsible for the chemical exchange linebroadening observed in the other regions of TIM. MD simulations confirm that motions of individual residues within the active site loop are correlated and suggest that the chemical exchange processes observed for residues in helix G arise from transitions between 3(10)- and alpha-helical structures. The results of the joint NMR and MD study provide global insight into the role of conformational dynamic processes in the function of TIM.

We examined the interfaces of homodimeric triosephosphate isomerase (TIM) from eight different species. The crystal structures of the enzymes showed that a portion of the interface is markedly similar in TIMs from Trypanosoma cruzi (TcTIM), Trypanosoma brucei, and Leishmania mexicana and significantly different from that of TIMs from human, yeast, chicken, Plasmodium falciparum, and Entamoeba histolytica. Since this interfacial region is central in the stability of TcTIM, we hypothesized that it would be possible to find agents that selectively affect the stability of TIMs from the three trypanosomatids. We found that 6,6'-bisbenzothiazole-2,2' diamine in the low micromolar range causes a desirable irreversible inactivation of the enzymes from the three trypanosomatids and has no effect on the other five TIMs. Thus, the data indicate that it is possible to find compounds that induce selective inactivation of the enzymes from three different trypanosomatids.

Molecular docking uses the three-dimensional structure of a receptor to screen databases of small molecules for potential ligands, often based on energetic complementarity. For many docking scoring functions, which calculate nonbonded interactions, metalloenzymes are challenging because of the partial covalent nature of metal-ligand interactions. To investigate how well molecular docking can identify potential ligands of metalloenzymes using a "standard" scoring function, we have docked the MDL Drug Data Report (MDDR), a functionally annotated database of 95,000 small molecules, against the X-ray crystal structures of five metalloenzymes. These enzymes included three zinc proteases, the nickel analogue of an iron enzyme, and a molybdenum metalloenzyme. The ability of the docking program to retrospectively enrich the annotated ligands as high-scoring hits for each enzyme and to calculate proper geometries was evaluated. In all five systems, the annotated ligands within the MDDR were enriched at least 20 times over random. To test the approach prospectively, a sixth target, the zinc beta-lactamase from Bacteroides fragilis, was screened against the fragment-like subset of the ZINC database. We purchased and tested 15 compounds from among the top 50 top-ranked ligands from docking, and found 5 inhibitors with apparent K(i) values less than 120 microM, the best of which was 2 microM. A more ambitious test still was predicting actual substrates for a seventh target, a Zn-dependent phosphotriesterase from Pseudomonas diminuta. Screening the Available Chemicals Directory (ACD) identified 25 thiophosphate esters as potential substrates within the top 100 ranked compounds. Eight of these, all previously uncharacterized for this enzyme, were acquired and tested, and all were confirmed experimentally as substrates. These results suggest that a simple, noncovalent scoring function may be used to identify inhibitors of at least some metalloenzymes.