Mechanistic implications of methylglyoxal synthase complexed with phosphoglycolohydroxamic acid as observed by X-ray crystallography and NMR spectroscopy

The possible role of methylglyoxal metabolism in cancer

Tumours reprogram their metabolism to acquire an evolutionary advantage over normal cells. However, not all such metabolic pathways support energy production. An example of these metabolic pathways is the Methylglyoxal (MG) one. This pathway helps maintain the redox state, and it might act as a phosphate sensor that monitors the intracellular phosphate levels. In this work, we discuss the biochemical step of the MG pathway and interrelate it with cancer.

A novel hyperthermophilic methylglyoxal synthase: molecular dynamic analysis on the regional fluctuations

Two putative methylglyoxal synthases, which catalyze the conversion of dihydroxyacetone phosphate to methylglyoxal, from Oceanithermus profundus DSM 14,977 and Clostridium difficile 630 have been characterized for activity and thermal stability. The enzyme from O. profundus was found to be hyperthermophilic, with the optimum activity at 80 °C and the residual activity up to 59% after incubation of 15 min at 95 °C, whereas the enzyme from C. difficile was mesophilic with the optimum activity at 40 °C and the residual activity less than 50% after the incubation at 55 °C or higher temperatures for 15 min. The structural analysis of the enzymes with molecular dynamics simulation indicated that the hyperthermophilic methylglyoxal synthase has a rigid protein structure with a lower overall root-mean-square-deviation value compared with the mesophilic or thermophilic counterparts. In addition, the simulation results identified distinct regions with high fluctuations throughout those of the mesophilic or thermophilic counterparts via root-mean-square-fluctuation analysis. Specific molecular interactions focusing on the hydrogen bonds and salt bridges in the distinct regions were analyzed in terms of interatomic distances and positions of the individual residues with respect to the secondary structures of the enzyme. Key interactions including specific salt bridges and hydrogen bonds between a rigid beta-sheet core and surrounding alpha helices were found to contribute to the stabilisation of the hyperthermophilic enzyme by reducing the regional fluctuations in the protein structure. The structural information and analysis approach in this study can be further exploited for the engineering and industrial application of the enzyme.

Hydrogen Bonds: Simple after All?

Hydrogen bonds play integral roles in biological structure, function, and conformational dynamics and are fundamental to life as it has evolved on Earth. However, our understanding of these fundamental and ubiquitous interactions has seemed fractured and incomplete, and it has been difficult to extract generalities and principles about hydrogen bonds despite thousands of papers published on this topic, perhaps in part because of the expanse of this subject and the density of studies. Fortunately, recent hydrogen bond proposals, discussions, and debates have stimulated new tests and models and have led to a remarkably simple picture of the structure of hydrogen bonds. This knowledge also provides clarity concerning hydrogen bond energetics, limiting and simplifying the factors that need be considered. Herein we recount the advances that have led to this simpler view of hydrogen bond structure, dynamics, and energetics. A quantitative predictive model for hydrogen bond length can now be broadly and deeply applied to evaluate current proposals and to uncover structural features of proteins, their conformational restraints, and their correlated motions. In contrast, a quantitative energetic description of molecular recognition and catalysis by proteins remains an important ongoing challenge, although our improved understanding of hydrogen bonds may aid in testing predictions from current and future models. We close by codifying our current state of understanding into five "Rules for Hydrogen Bonding" that may provide a foundation for understanding and teaching about these vital interactions and for building toward a deeper understanding of hydrogen bond energetics.

Structural basis for the regulatory interaction of the methylglyoxal synthase MgsA with the carbon flux regulator Crh in

Utilization of energy-rich carbon sources such as glucose is fundamental to the evolutionary success of bacteria. Glucose can be catabolized via glycolysis for feeding the intermediary metabolism. The methylglyoxal synthase MgsA produces methylglyoxal from the glycolytic intermediate dihydroxyacetone phosphate. Methylglyoxal is toxic, requiring stringent regulation of MgsA activity. In the Gram-positive bacterium Bacillus subtilis, an interaction with the phosphoprotein Crh controls MgsA activity. In the absence of preferred carbon sources, Crh is present in the nonphosphorylated state and binds to and thereby inhibits MgsA. To better understand the mechanism of regulation of MgsA, here we performed biochemical and structural analyses of B. subtilis MgsA and of its interaction with Crh. Our results indicated that MgsA forms a hexamer (i.e. a trimer of dimers) in the crystal structure, whereas it seems to exist in an equilibrium between a dimer and hexamer in solution. In the hexamer, two alternative dimers could be distinguished, but only one appeared to prevail in solution. Further analysis strongly suggested that the hexamer is the biologically active form. In vitro cross-linking studies revealed that Crh interacts with the N-terminal helices of MgsA and that the Crh-MgsA binding inactivates MgsA by distorting and thereby blocking its active site. In summary, our results indicate that dimeric and hexameric MgsA species exist in an equilibrium in solution, that the hexameric species is the active form, and that binding to Crh deforms and blocks the active site in MgsA.

Crystal Structures of the Iron-Sulfur Cluster-Dependent Quinolinate Synthase in Complex with Dihydroxyacetone Phosphate, Iminoaspartate Analogues, and Quinolinate

The quinolinate synthase of prokaryotes and photosynthetic eukaryotes, NadA, contains a [4Fe-4S] cluster with unknown function. We report crystal structures of Pyrococcus horikoshii NadA in complex with dihydroxyacetone phosphate (DHAP), iminoaspartate analogues, and quinolinate. DHAP adopts a nearly planar conformation and chelates the [4Fe-4S] cluster via its keto and hydroxyl groups. The active site architecture suggests that the cluster acts as a Lewis acid in enediolate formation, like zinc in class II aldolases. The DHAP and putative iminoaspartate structures suggest a model for a condensed intermediate. The ensemble of structures suggests a two-state system, which may be exploited in early steps.

Automated protein motif generation in the structure-based protein function prediction tool ProMOL

ProMOL, a plugin for the PyMOL molecular graphics system, is a structure-based protein function prediction tool. ProMOL includes a set of routines for building motif templates that are used for screening query structures for enzyme active sites. Previously, each motif template was generated manually and required supervision in the optimization of parameters for sensitivity and selectivity. We developed an algorithm and workflow for the automation of motif building and testing routines in ProMOL. The algorithm uses a set of empirically derived parameters for optimization and requires little user intervention. The automated motif generation algorithm was first tested in a performance comparison with a set of manually generated motifs based on identical active sites from the same 112 PDB entries. The two sets of motifs were equally effective in identifying alignments with homologs and in rejecting alignments with unrelated structures. A second set of 296 active site motifs were generated automatically, based on Catalytic Site Atlas entries with literature citations, as an expansion of the library of existing manually generated motif templates. The new motif templates exhibited comparable performance to the existing ones in terms of hit rates against native structures, homologs with the same EC and Pfam designations, and randomly selected unrelated structures with a different EC designation at the first EC digit, as well as in terms of RMSD values obtained from local structural alignments of motifs and query structures. This research is supported by NIH grant GM078077.

Rationalization of allosteric pathway in Thermus sp. GH5 methylglyoxal synthase

A sequence of 10 amino acids at the C-terminus region of methylglyoxal synthase from Escherichia coli (EMGS) provides an arginine, which plays a crucial role in forming a salt bridge with a proximal aspartate residue in the neighboring subunit, consequently transferring the allosteric signal between subunits. In order to verify the role of arginine, the gene encoding MGS from a thermophile species, Thermus sp. GH5 (TMGS) lacking this arginine was cloned with an additional 30 bp sequence at the 3´-end and then expressed in form of a fusion TMGS with a 10 residual segment at the C-terminus (TMGS(+)). The resulting recombinant enzyme showed a significant increase in cooperativity towards phosphate, reflected by a change in the Hill coefficient (nH) from 1.5 to 1.99. Experiments including site directed mutagenesis for Asp-10 in TMGS and TMGS(+), two dimentional structural survey, fluorescence and irreversible thermoinactivation were carried out to confirm this pathway.

Synthesis and biochemical evaluation of selective inhibitors of class II fructose bisphosphate aldolases: towards new synthetic antibiotics

We report the synthesis and biochemical evaluation of selective inhibitors of class II (zinc-dependent) fructose bisphosphate aldolases. The most active compound is a simplified analogue of fructose bisphosphate, bearing a well-positioned metal chelating group. It is a powerful and highly selective competitive inhibitor of isolated class II aldolases. We report crystallographic studies of this inhibitor bound in the active site of the Helicobacter pylori enzyme. The compound also shows activity against Mycobacterium tuberculosis isolates.

A critical evaluation of different QM/MM frontier treatments with SCC-DFTB as the QM method

The performance of different link atom based frontier treatments in QM/MM simulations was evaluated critically with SCC-DFTB as the QM method. In addition to the analysis of gas-phase molecules as in previous studies, an important element of the present work is that chemical reactions in realistic enzyme systems were also examined. The schemes tested include all options available in the program CHARMM for SCC-DFTB/MM simulation, which treat electrostatic interactions due to the MM atoms close to the QM/MM boundary in different ways. In addition, a new approach, the divided frontier charge (DIV), has been implemented in which the partial charge associated with the frontier MM atom ("link host") is evenly distributed to the other MM atoms in the same group. The performance of these schemes was evaluated based on properties including proton affinities, deprotonation energies, dipole moments, and energetics of proton transfer reactions. Similar to previous work, it was found that calculated proton affinities and deprotonation energies of alcohols, carbonic acids, amino acids, and model DNA bases are very sensitive to the link atom scheme; the commonly used single link atom approach often gives error on the order of 15 to 20 kcal/mol. Other schemes give better and, on average, mutually comparable results. For proton transfer reactions, encouragingly, both activation barriers and reaction energies are fairly insensitive (within a typical range of 2-4 kcal/mol) to the link atom scheme due to error cancellation, and this was observed for both gas-phase and enzyme systems. Therefore, the effect of using different link atom schemes in QM/MM simulations is rather small for chemical reactions that conserve the total charge. Although the current study used an approximate DFT method as the QM level, the observed trends are expected to be applicable to QM/MM methods with use of other QM approaches. This observation does not mean to encourage QM/MM simulations without careful benchmark in the study of specific systems, rather it emphasizes that other technical details, such as the treatment of long-range electrostatics, tend to play a more important role and need to be handled carefully.

Selective prediction of interaction sites in protein structures with THEMATICS

BACKGROUND

Methods are now available for the prediction of interaction sites in protein 3D structures. While many of these methods report high success rates for site prediction, often these predictions are not very selective and have low precision. Precision in site prediction is addressed using Theoretical Microscopic Titration Curves (THEMATICS), a simple computational method for the identification of active sites in enzymes. Recall and precision are measured and compared with other methods for the prediction of catalytic sites.

RESULTS

Using a test set of 169 enzymes from the original Catalytic Residue Dataset (CatRes) it is shown that THEMATICS can deliver precise, localised site predictions. Furthermore, adjustment of the cut-off criteria can improve the recall rates for catalytic residues with only a small sacrifice in precision. Recall rates for CatRes/CSA annotated catalytic residues are 41.1%, 50.4%, and 54.2% for Z score cut-off values of 1.00, 0.99, and 0.98, respectively. The corresponding precision rates are 19.4%, 17.9%, and 16.4%. The success rate for catalytic sites is higher, with correct or partially correct predictions for 77.5%, 85.8%, and 88.2% of the enzymes in the test set, corresponding to the same respective Z score cut-offs, if only the CatRes annotations are used as the reference set. Incorporation of additional literature annotations into the reference set gives total success rates of 89.9%, 92.9%, and 94.1%, again for corresponding cut-off values of 1.00, 0.99, and 0.98. False positive rates for a 75-protein test set are 1.95%, 2.60%, and 3.12% for Z score cut-offs of 1.00, 0.99, and 0.98, respectively.

CONCLUSION

With a preferred cut-off value of 0.99, THEMATICS achieves a high success rate of interaction site prediction, about 86% correct or partially correct using CatRes/CSA annotations only and about 93% with an expanded reference set. Success rates for catalytic residue prediction are similar to those of other structure-based methods, but with substantially better precision and lower false positive rates. THEMATICS performs well across the spectrum of E.C. classes. The method requires only the structure of the query protein as input. THEMATICS predictions may be obtained via the web from structures in PDB format at: http://pfweb.chem.neu.edu/thematics/submit.html.

New competitive inhibitors of cytosolic (NADH-dependent) rabbit muscle glycerophosphate dehydrogenase

We report the synthesis and biochemical evaluation of new competitive inhibitors of the cytosolic (NADH-dependent) glycerophosphate dehydrogenase. The best tested compound, phosphono-propionohydroxamic acid, with a Ki of 6 microM, might be of interest as an anti-obesity drug.

The crystal structure of E.coli 1-deoxy-D-xylulose-5-phosphate reductoisomerase in a ternary complex with the antimalarial compound fosmidomycin and NADPH reveals a tight-binding closed enzyme conformation

The key enzyme in the non-mevalonate pathway of isoprenoid biosynthesis, 1-deoxy-D-xylulose 5-phosphate reductoisomerase (DXR) has been shown to be the target enzyme of fosmidomycin, an antimalarial, antibacterial and herbicidal compound. Here we report the crystal structure of selenomethionine-labelled Escherichia coli DXR in a ternary complex with NADPH and fosmidomycin at 2.2 A resolution. The structure reveals a considerable conformational rearrangement upon fosmidomycin binding and provides insights into the slow, tight binding inhibition mode of the inhibitor. Although the inhibitor displays an unusual non-metal mediated mode of inhibition, which is an artefact most likely due to the low metal affinity of DXR at the pH used for crystallization, the structural data add valuable information for the rational design of novel DXR inhibitors. Using this structure together with the published structural data and the 1.9 A crystal structure of DXR in a ternary complex with NADPH and the substrate 1-deoxy-D-xylulose 5-phosphate, a model for the physiologically relevant tight-binding mode of inhibition is proposed. The structure of the substrate complex must be interpreted with caution due to the presence of a second diastereomer in the active site.

Optimal alignment for enzymatic proton transfer: structure of the Michaelis complex of triosephosphate isomerase at 1.2-A resolution

In enzyme catalysis, where exquisitely positioned functionality is the sine qua non, atomic coordinates for a Michaelis complex can provide powerful insights into activation of the substrate. We focus here on the initial proton transfer of the isomerization reaction catalyzed by triosephosphate isomerase and present the crystal structure of its Michaelis complex with the substrate dihydroxyacetone phosphate at near-atomic resolution. The active site is highly compact, with unusually short and bifurcated hydrogen bonds for both catalytic Glu-165 and His-95 residues. The carboxylate oxygen of the catalytic base Glu-165 is positioned in an unprecedented close interaction with the ketone and the alpha-hydroxy carbons of the substrate (C em leader O approximately 3.0 A), which is optimal for the proton transfer involving these centers. The electrophile that polarizes the substrate, His-95, has close contacts to the substrate's O1 and O2 (N em leader O < or = 3.0 and 2.6 A, respectively). The substrate is conformationally relaxed in the Michaelis complex: the phosphate group is out of the plane of the ketone group, and the hydroxy and ketone oxygen atoms are not in the cisoid configuration. The epsilon ammonium group of the electrophilic Lys-12 is within hydrogen-bonding distance of the substrate's ketone oxygen, the bridging oxygen, and a terminal phosphate's oxygen, suggesting a role for this residue in both catalysis and in controlling the flexibility of active-site loop.

Functional specificities of methylglyoxal synthase and triosephosphate isomerase: a combined QM/MM analysis

Combined SCC-DFTB/CHARMM calculations were carried out to analyze the origin for the functional specificities of triosephosphate isomerase (TIM) and methylglyoxal synthase (MGS). The two enzymes bind to the same substrate, dihydroxyacetone phosphate (DHAP), and have rather similar active sites. However, they catalyze different reactions; TIM catalyzes the isomerization of DHAP to glyceraldehyde 3-phosphate (GAP), while MGS catalyzes the elimination of phosphate from DHAP. Similar to previous suggestions, the calculations confirmed that GAP formation is prohibited in MGS due primarily to the reduced flexibility of the catalytic base (Asp 71) compared to that in TIM (Glu 165). For the suppression of phosphate elimination in TIM, the calculations show that the widely accepted stereoelectronic argument that invokes the different phosphoryl torsion angles observed in the X-ray structures of inhibitor complexes of the two enzymes is not as important as electrostatic contributions from the protein and water molecules surrounding the phosphoryl.

Hydroxamates as substrates and inhibitors for FMN-dependent 2-hydroxy acid dehydrogenases

Long-chain hydroxy acid oxydase (HAO) is a member of a flavoenzyme family with significant amino acid sequence similarity and strongly conserved three-dimensional structure; in particular, active-site amino acids involved in catalysis are invariant, with one exception, and numerous enzymatic studies suggest an identical chemical mechanism involving an intermediate carbanion for all family members. Known physiological substrates are a variety of L-2-hydroxy acids. Peroxisomal HAO differs from the other family members in that its actual physiological substrate is not known; it was first described as an L-amino acid oxidase, and recently was identified as an enzyme that converts creatol (hydroxycreatinine) to methylguanidine (a metabolite involved in a variety of uremic syndromes). Creatol (2-amino-5-hydroxy-1-methyl-4(5H)imidazolone) is not a 2-hydroxy acid. We show in this work that 2-hydroxyphenyl acetohydroxamate (HYPAH, the hydroxamate of mandelic acid), a compound that bears similarity both to mandelate (one of the best substrates known) and to creatol, is turned over by HAO, but between 10- and 100-fold less efficiently than mandelate itself. The compound also binds to the active site of homologous flavocytochrome b(2) (L-lactate dehydrogenase). Comparative pH-rate studies for mandelate and its hydroxamate suggest that HYPAH may bind in its ionized form. Both pH-rate profiles are bell-shaped curves, as are those determined for two other family members, flavocytochrome b(2) and mandelate dehydrogenase; while the group with an acid pK(a) between 5 and 6 is most likely the active-site histidine (the residue which abstracts the substrate C2 proton), the identity of the basic group is less clear. It has been proposed to be one of the active site arginines (Lehoux, I., and Mitra, B. (1999) Biochemistry38, 5836-5848); we suggest as an alternative that it could be the lysine residue that interacts with the flavin N1 and O2 positions and stabilizes the negative charge of reduced flavin. In addition to these studies, we have found that HAO is competitively inhibited by benzohydroxamate, which is one atom shorter than HYPAH; its affinity is nearly 100-fold lower than that of the substrate, in contrast to the strong inhibition it exerts on mandelate racemase (Maurice, St. M., and Bearne, S. L. (2000) Biochemistry39, 13324-13335). In the latter case, the 100-fold higher affinity compared to mandelate was proposed to arise from the fact that the hydroxamate can mimic the enolic intermediate which lies on the reaction pathway after C2 proton abstraction. Thus our results do not support the existence of a similar enolic intermediate for HAO (and probably its homologues), although they do not disprove it.