Crystal structure of the mutant yeast triosephosphate isomerase in which the catalytic base glutamic acid 165 is changed to aspartic acid

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.

Therapeutic Targeting of Cancer Metabolism with Triosephosphate Isomerase

The increase in glycolytic flux in cancer, known as aerobic glycolysis, is one of the most important hallmarks of cancer. Therefore, glycolytic enzymes have importance in understanding the molecular mechanism of cancer progression. Triosephosphate isomerase (TPI) is one of the key glycolytic enzymes. Furthermore, it takes a part in gluconeogenesis, pentose phosphate pathway and fatty acid biosynthesis. To date, it has been shown altered levels of TPI in various cancer types, especially in metastatic phenotype. According to other studies, TPI might be considered as a potential therapeutic target and a cancer-related biomarker in different types of cancer. However, its function in tumor formation and development has not been fully understood. Here, we reviewed the relationship between TPI and cancer for the first time.

Integrated Use of Biochemical, Native Mass Spectrometry, Computational, and Genome-Editing Methods to Elucidate the Mechanism of a Salmonella deglycase

Salmonellais a foodborne pathogen that causes annually millions of cases of salmonellosis globally, yet Salmonella-specific antibacterials are not available. During inflammation, Salmonella utilizes the Amadori compound fructose-asparagine (F-Asn) as a nutrient through the successive action of three enzymes, including the terminal FraB deglycase. Salmonella mutants lacking FraB are highly attenuated in mouse models of inflammation due to the toxic build-up of the substrate 6-phosphofructose-aspartate (6-P-F-Asp). This toxicity makes Salmonella FraB an appealing drug target, but there is currently little experimental information about its catalytic mechanism. Therefore, we sought to test our postulated mechanism for the FraB-catalyzed deglycation of 6-P-F-Asp (via an enaminol intermediate) to glucose-6-phosphate and aspartate. A FraB homodimer model generated by RosettaCM was used to build substrate-docked structures that, coupled with sequence alignment of FraB homologs, helped map a putative active site. Five candidate active-site residues-including three expected to participate in substrate binding-were mutated individually and characterized. Native mass spectrometry and ion mobility were used to assess collision cross sections and confirm that the quaternary structure of the mutants mirrored the wild type, and that there are two active sites/homodimer. Our biochemical studies revealed that FraB Glu214Ala, Glu214Asp, and His230Ala were inactive in vitro, consistent with deprotonated-Glu214 and protonated-His230 serving as a general base and a general acid, respectively. Glu214Ala or His230Ala introduced into the Salmonella chromosome by CRISPR/Cas9-mediated genome editing abolished growth on F-Asn. Results from our computational and experimental approaches shed light on the catalytic mechanism of Salmonella FraB and of phosphosugar deglycases in general.

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.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

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.

Rescue of K12G triosephosphate isomerase by ammonium cations: the reaction of an enzyme in pieces

The K12G mutation at yeast triosephosphate isomerase (TIM) results in a 5.5 × 10(5)-fold decrease in k(cat)/K(m) for isomerization of glyceraldehyde 3-phosphate, and the activity of this mutant can be successfully "rescued" by NH(4)(+) and primary alkylammonium cations. The transition state for the K12G mutant TIM-catalyzed reaction is stabilized by 1.5 kcal/mol by interaction with NH(4)(+). The larger 3.9 kcal/mol stabilization by CH(3)CH(2)CH(2)CH(2)NH(3)(+) is due to hydrophobic interactions between the mutant enzyme and the butyl side chain of the cation activator. There is no significant transfer of a proton from alkylammonium cations to GAP at the transition state for the K12G mutant TIM-catalyzed reaction, because activation by a series of RNH(3)(+) shows little or no dependence on the pK(a) of RNH(3)(+). A comparison of k(cat)/K(m) = 6.6 × 10(6) M(-1) s(-1) for the wildtype TIM-catalyzed isomerization of GAP and the third-order rate constant of 150 M(-2) s(-1) for activation by NH(4)(+) of the K12G mutant TIM-catalyzed isomerization shows that stabilization of the bound transition state by the effectively intramolecular interaction of the cationic side chain of Lys-12 at wildtype TIM is 6.3 kcal/mol greater than that for the corresponding intermolecular interaction of NH(4)(+) at K12G mutant TIM.

Testing geometrical discrimination within an enzyme active site: constrained hydrogen bonding in the ketosteroid isomerase oxyanion hole

Enzymes are classically proposed to accelerate reactions by binding substrates within active-site environments that are structurally preorganized to optimize binding interactions with reaction transition states rather than ground states. This is a remarkably formidable task considering the limited 0.1-1 A scale of most substrate rearrangements. The flexibility of active-site functional groups along the coordinate of substrate rearrangement, the distance scale on which enzymes can distinguish structural rearrangement, and the energetic significance of discrimination on that scale remain open questions that are fundamental to a basic physical understanding of enzyme active sites and catalysis. We bring together 1.2-1.5 A resolution X-ray crystallography, (1)H and (19)F NMR spectroscopy, quantum mechanical calculations, and transition-state analogue binding measurements to test the distance scale on which noncovalent forces can constrain the structural relaxation or translation of side chains and ligands along a specific coordinate and the energetic consequences of such geometric constraints within the active site of bacterial ketosteroid isomerase (KSI). Our results strongly suggest that packing and binding interactions within the KSI active site can constrain local side-chain reorientation and prevent hydrogen bond shortening by 0.1 A or less. Further, this constraint has substantial energetic effects on ligand binding and stabilization of negative charge within the oxyanion hole. These results provide evidence that subtle geometric effects, indistinguishable in most X-ray crystallographic structures, can have significant energetic consequences and highlight the importance of using synergistic experimental approaches to dissect enzyme function.

Stress resistance and signal fidelity independent of nuclear MAPK function

Elevated external solute stimulates a conserved MAPK cascade that elicits responses that maintain osmotic balance. The yeast high-osmolarity glycerol (HOG) pathway activates Hog1 MAPK (mammalian ortholog p38alpha/SAPKalpha), which enters the nucleus and induces expression of >50 genes, implying that transcriptional up-regulation is necessary to cope with hyperosmotic stress. Contrary to this expectation, we show here that cells lacking the karyopherin required for Hog1 nuclear import or in which Hog1 is anchored at the plasma membrane (or both) can withstand long-term hyperosmotic challenge by ionic and nonionic solutes without exhibiting the normal change in transcriptional program (comparable with hog1Delta cells), as judged by mRNA hybridization and microarray analysis. For such cells to survive hyperosmotic stress, systematic genetic analysis ruled out the need for any Hog1-dependent transcription factor, the Hog1-activated MAPKAP kinases, or ion, glycerol, and water channels. By contrast, enzymes needed for glycerol production were essential for viability. Thus, control of intracellular glycerol formation by Hog1 is critical for maintenance of osmotic balance but not transcriptional induction of any gene.

Slow proton transfer from the hydrogen-labelled carboxylic acid side chain (Glu-165) of triosephosphate isomerase to imidazole buffer in D2O

The catalytic base at the active site of triosephosphate isomerase (TIM) was labelled with -H by abstraction of a proton from substrate d-glyceraldehyde 3-phosphate to form an enzyme-bound enediol(ate) in D2O solvent. The partitioning of this labelled enzyme between intramolecular transfer of -H to form dihydroxyacetone phosphate (DHAP), and irreversible exchange with -D from solvent was examined by determining the yields of H- and D-labelled products by 1H NMR spectroscopy. The yield of hydrogen-labelled product DHAP remains constant as the concentration of the basic form of imidazole buffer is increased from 0.014 to 0.56 M. This shows that the active site of free TIM, which has an open conformation needed to allow substrate binding, adopts a closed conformation at the enediolate-complex intermediate where the catalytic side chain is sequestered from interaction with imidazole dissolved in D2O.

Local encoding of computationally designed enzyme activity

One aim of computational protein design is to introduce novel enzyme activity into proteins of known structure by predicting mutations that stabilize transition states. Previously, we showed that it is possible to introduce triose phosphate isomerase activity into the ribose-binding protein of Escherichia coli by constructing 17 mutations in the first two layers of residues that surround the wild-type ligand-binding site. Here, we report that these mutations can be "transplanted" into a homologous ribose-binding protein, isolated from the hyperthermophilic bacterium Thermoanaerobacter tengcongensis, with retention of catalytic activity, substrate affinity, and reaction pH dependence. The observed 10(5)-10(6)-fold rate enhancement corresponds to 70% of the maximally known transition-state binding energy. The wild-type sequences in these two homologues are almost perfectly conserved in the vicinity of their ribose-binding sites, but diverge significantly at increasing distance from these sites. The results demonstrate that the computationally designed mutations are sufficient to encode the observed enzyme activity, that all the observed activity is encoded locally within the layer of residues directly in contact with the substrate and that, in this case, at least 70% of transition state stabilization energy can be achieved using straightforward considerations of stereochemical complementarity between enzyme and reactants.

AB INITIO QUANTUM CHEMICAL AND MIXED QUANTUM MECHANICS/MOLECULAR MECHANICS (QM/MM) METHODS FOR STUDYING ENZYMATIC CATALYSIS

We describe large scale ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic reactions. First, technical aspects of the methodology are reviewed, including the hybrid density functional theory (DFT) methods that are typically employed for the QM aspect of the calculations, and various approaches to defining the interface between the QM and MM regions in QM/MM approaches. The modeling of the enzymatic catalytic cycle for three examples--methane monooxygenase, cytochrome P450, and triose phosphate isomerase--are discussed in some depth, followed by a brief summary of other systems that have been investigated by ab initio methods over the past several years. Finally, a discussion of the qualitative and quantitative conclusions concerning enzymatic catalysis that are available from modern ab initio approaches is presented, followed by a conclusion briefly summarizing future prospects.

Structural origins of efficient proton abstraction from carbon by a catalytic antibody

Antibody 34E4 catalyzes the conversion of benzisoxazoles to salicylonitriles with high rates and multiple turnovers. The crystal structure of its complex with the benzimidazolium hapten at 2.5-angstroms resolution shows that a combination of hydrogen bonding, pi stacking, and van der Waals interactions is exploited to position both the base, Glu(H50), and the substrate for efficient proton transfer. Suboptimal placement of the catalytic carboxylate, as observed in the 2.8-angstroms structure of the Glu(H50)Asp variant, results in substantially reduced catalytic efficiency. In addition to imposing high positional order on the transition state, the antibody pocket provides a highly structured microenvironment for the reaction in which the carboxylate base is activated through partial desolvation, and the highly polarizable transition state is stabilized by dispersion interactions with the aromatic residue Trp(L91) and solvation of the leaving group oxygen by external water. The enzyme-like efficiency of general base catalysis in this system directly reflects the original hapten design, in which a charged guanidinium moiety was strategically used to elicit an accurately positioned functional group in an appropriate reaction environment and suggests that even larger catalytic effects may be achievable by extending this approach to the induction of acid-base pairs capable of bifunctional catalysis.

Potential use of sugar binding proteins in reactors for regeneration of CO2 fixation acceptor D-Ribulose-1,5-bisphosphate

Sugar binding proteins and binders of intermediate sugar metabolites derived from microbes are increasingly being used as reagents in new and expanding areas of biotechnology. The fixation of carbon dioxide at emission source has recently emerged as a technology with potentially significant implications for environmental biotechnology. Carbon dioxide is fixed onto a five carbon sugar D-ribulose-1,5-bisphosphate. We present a review of enzymatic and non-enzymatic binding proteins, for 3-phosphoglycerate (3PGA), 3-phosphoglyceraldehyde (3PGAL), dihydroxyacetone phosphate (DHAP), xylulose-5-phosphate (X5P) and ribulose-1,5-bisphosphate (RuBP) which could be potentially used in reactors regenerating RuBP from 3PGA. A series of reactors combined in a linear fashion has been previously shown to convert 3-PGA, (the product of fixed CO₂ on RuBP as starting material) into RuBP (Bhattacharya et al., 2004; Bhattacharya, 2001). This was the basis for designing reactors harboring enzyme complexes/mixtures instead of linear combination of single-enzyme reactors for conversion of 3PGA into RuBP. Specific sugars in such enzyme-complex harboring reactors requires removal at key steps and fed to different reactors necessitating reversible sugar binders. In this review we present an account of existing microbial sugar binding proteins and their potential utility in these operations.

Computational modeling of the catalytic reaction in triosephosphate isomerase

We present a comprehensive analysis of the catalytic cycle of the enzyme triosephosphate isomerase (TIM), including both the reactive chemistry and the catalytic loop and side-chain motions. Combining accurate mixed quantum mechanics/molecular mechanics (QM/MM) and protein structure prediction methods, we have modeled both the structural and chemical aspects of the reversible isomerization of dihydroxyacetone phosphate (DHAP) to d-glyceraldehyde 3-phosphate (GAP), for which there is a wealth of experimental data. The conjunction of this novel computational approach with the use of the recent near-atomic resolution TIM-DHAP Michaelis complex PDB structure, 1NEY.pdb, has enabled us to obtain robust qualitative and, where available, quantitative agreement with a wide range of experimental data. Among the principal conclusions that we are able to draw are the importance of the monoanionic (as opposed to dianioic) form of the substrate phosphate group in the catalytic cycle, detailed positioning and energetics of the key catalytic residues in the active-site, the flexible nature of Glu165, which favors its direct involvement in the formation of the enediol intermediate, energetics of the open and closed form of the catalytic loop region in the presence and absence of substrate, and quantitative reproduction of various experimentally measured reaction rates, typically to within approximately 1 kcal/mol. Our results are consistent with the available experimental data, and provide an initial picture as to why loop opening when GAP is the product has a higher barrier than when DHAP is the product.

The importance of GLU361 position in the reaction catalyzed by cholesterol oxidase

Cholesterol oxidase stereospecifically isomerizes cholest-5-en-3-one to cholest-4-en-3-one. When the base catalyst for isomerization, Glu361, is mutated to Asp, the rate of deprotonation of cholest-5-en-3-one is not affected, but protonation of the dienolic intermediate becomes rate-limiting. This may be a consequence of the large distance between the catalytic base and carbon-6 of the intermediate in the mutant enzyme.

An interface point-mutation variant of triosephosphate isomerase is compactly folded and monomeric at low protein concentrations

Wild-type trypanosomal triosephosphate isomerase (wtTIM) is a very tight dimer. The interface residue His-47 of wtTIM has been mutated into an asparagine. Ultracentrifugation data show that this variant (H47N) only dimerises at protein concentrations above 3 mg/ml. H47N has been characterised at a protein concentration where it is predominantly a monomer. Circular dichroism measurements in the near-UV and far-UV show that this monomer is a compactly folded protein with secondary structure similar as in wtTIM. The thermal stability of the monomeric H47N is decreased compared to wtTIM: temperature gradient gel electrophoresis (TGGE) measurements give Tm-values of 41 degrees C for wtTIM, whereas the Tm-value for the monomeric form of H47N is approximately 7 degrees C lower.