Crystal structure of the H256A mutant of rat testis fructose-6-phosphate,2-kinase/fructose-2,6-bisphosphatase. Fructose 6-phosphate in the active site leads to mechanisms for both mutant and wild type bisphosphatase activities

PFKFB4 Drives the Oncogenicity in TP53-Mutated Hepatocellular Carcinoma in a Phosphatase-Dependent Manner

BACKGROUND & AIMS

Metabolic reprogramming is recognized as a cancer hallmark intimately linked to tumor hypoxia, which supports rapid tumor growth and mitigates the consequential oxidative stress. Phosphofructokinase-fructose bisphosphatase (PFKFB) is a family of bidirectional glycolytic enzymes possessing both kinase and phosphatase functions and has emerged as important oncogenes in multiple types of cancer. However, its clinical relevance, functional significance, and underlying mechanistic insights in hepatocellular carcinoma (HCC), the primary malignancy that develops in the most important metabolic organ, has never been addressed.

METHODS

PFKFB4 expression was examined by RNA sequencing in The Cancer Genome Atlas and our in-house HCC cohort. The up-regulation of PFKFB4 expression was confirmed further by quantitative polymerase chain reaction in an expanded hepatitis B virus-associated HCC cohort followed by clinicopathologic correlation analysis. CRISPR/Cas9-mediated PFKFB4 knockout cells were generated for functional characterization in vivo, targeted metabolomic profiling, as well as RNA sequencing analysis to comprehensively examine the impact of PFKFB4 loss in HCC.

RESULTS

PFKFB4 expression was up-regulated significantly in HCC and correlated positively with TP53 and TSC2 loss-of-function mutations. In silico transcriptome-based analysis further revealed PFKFB4 functions as a critical hypoxia-inducible gene. Clinically, PFKFB4 up-regulation was associated with more aggressive tumor behavior. Functionally, CRISPR/Cas9-mediated PFKFB4 knockout significantly impaired in vivo HCC development. Targeted metabolomic profiling revealed that PFKFB4 functions as a phosphatase in HCC and its ablation caused an accumulation of metabolites in downstream glycolysis and the pentose phosphate pathway. In addition, PFKFB4 loss induced hypoxia-responsive genes in glycolysis and reactive oxygen species detoxification. Conversely, ectopic PFKFB4 expression conferred sorafenib resistance.

CONCLUSIONS

PFKFB4 up-regulation supports HCC development and posed therapeutic implications.

Role of PFKFB3 and PFKFB4 in Cancer: Genetic Basis, Impact on Disease Development/Progression, and Potential as Therapeutic Targets

Glycolysis is a crucial metabolic process in rapidly proliferating cells such as cancer cells. Phosphofructokinase-1 (PFK-1) is a key rate-limiting enzyme of glycolysis. Its efficiency is allosterically regulated by numerous substances occurring in the cytoplasm. However, the most potent regulator of PFK-1 is fructose-2,6-bisphosphate (F-2,6-BP), the level of which is strongly associated with 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase activity (PFK-2/FBPase-2, PFKFB). PFK-2/FBPase-2 is a bifunctional enzyme responsible for F-2,6-BP synthesis and degradation. Four isozymes of PFKFB (PFKFB1, PFKFB2, PFKFB3, and PFKFB4) have been identified. Alterations in the levels of all PFK-2/FBPase-2 isozymes have been reported in different diseases. However, most recent studies have focused on an increased expression of PFKFB3 and PFKFB4 in cancer tissues and their role in carcinogenesis. In this review, we summarize our current knowledge on all PFKFB genes and protein structures, and emphasize important differences between the isoenzymes, which likely affect their kinase/phosphatase activities. The main focus is on the latest reports in this field of cancer research, and in particular the impact of PFKFB3 and PFKFB4 on tumor progression, metastasis, angiogenesis, and autophagy. We also present the most recent achievements in the development of new drugs targeting these isozymes. Finally, we discuss potential combination therapies using PFKFB3 inhibitors, which may represent important future cancer treatment options.

Mutations in the phosphatase domain of the 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase result in the transcriptional activation of the alternative oxidase and gluconeogenic pathways in Podospora anserina

By screening suppressors of a respiratory mutant lacking a functional cytochrome pathway in the filamentous fungus Podospora anserina, we isolated a mutation located in the phosphatase domain of the bi-functional enzyme 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase (PFK-2/FBPase-2). We show that the inactivation of the phosphatase but not of the kinase domain is responsible for the suppressor effect that results from the activation of the RSEs transcription factors that control expression of AOX, an alternative oxidase able to bypass the mitochondria cytochrome pathway of respiration. Remarkably, activation of the RSEs also stimulates the expression of the gluconeogenic enzymes, fructose-1,6 bi-phosphatase (FBPase-1) and phosphoenolpyruvate carboxykinase (PCK-1). We thus reveal in P. anserina an apparently paradoxical situation where the inactivation of the phosphatase domain of PFK-2/FBPase-2, supposed to stimulate glycolysis, is correlated with the transcriptional induction of the gluconeogenic enzymes. Phylogenic analysis revealed the presence of multiple presumed PFK-2/FBPase-2 isoforms in all the species of tested Ascomycetes.

Caught before Released: Structural Mapping of the Reaction Trajectory for the Sofosbuvir Activating Enzyme, Human Histidine Triad Nucleotide Binding Protein 1 (hHint1)

Human histidine triad nucleotide binding protein 1 (hHint1) is classified as an efficient nucleoside phosphoramidase and acyl-adenosine monophosphate hydrolase. Human Hint1 has been shown to be essential for the metabolic activation of nucleotide antiviral pronucleotides (i.e., proTides), such as the FDA approved hepatitis C drug, sofosbuvir. The active site of hHint1 comprises an ensemble of strictly conserved histidines, including nucleophilic His112. To structurally investigate the mechanism of hHint1 catalysis, we have designed and prepared nucleoside thiophosphoramidate substrates that are able to capture the transiently formed nucleotidylated-His112 intermediate (E*) using time-dependent crystallography. Utilizing a catalytically inactive hHint1 His112Asn enzyme variant and wild-type enzyme, the enzyme-substrate (ES¹) and product (EP²) complexes were also cocrystallized, respectively, thus providing a structural map of the reaction trajectory. On the basis of these observations and the mechanistic necessity of proton transfers, proton inventory studies were carried out. Although we cannot completely exclude the possibility of more than one proton in flight, the results of these studies were consistent with the transfer of a single proton during the formation of the intermediate. Interestingly, structural analysis revealed that the critical proton transfers required for intermediate formation and hydrolysis may be mediated by a conserved active site water channel. Taken together, our results provide mechanistic insights underpinning histidine nucleophilic catalysis in general and hHint1 catalysis, in particular, thus aiding the design of future proTides and the elucidation of the natural function of the Hint family of enzymes.

Crystal Structure of Xanthomonas AvrRxo1-ORF1, a Type III Effector with a Polynucleotide Kinase Domain, and Its Interactor AvrRxo1-ORF2

Xanthomonas oryzae pv. oryzicola (Xoc) causes bacterial leaf streak (BLS) disease on rice plants. Xoc delivers a type III effector AvrRxo1-ORF1 into rice plant cells that can be recognized by disease resistance (R) protein Rxo1, and triggers resistance to BLS disease. However, the mechanism and virulence role of AvrRxo1 is not known. In the genome of Xoc, AvrRxo1-ORF1 is adjacent to another gene AvrRxo1-ORF2, which was predicted to encode a molecular chaperone of AvrRxo1-ORF1. We report the co-purification and crystallization of the AvrRxo1-ORF1:AvrRxo1-ORF2 tetramer complex at 1.64 Å resolution. AvrRxo1-ORF1 has a T4 polynucleotide kinase domain, and expression of AvrRxo1-ORF1 suppresses bacterial growth in a manner dependent on the kinase motif. Although AvrRxo1-ORF2 binds AvrRxo1-ORF1, it is structurally different from typical effector-binding chaperones, in that it has a distinct fold containing a novel kinase-binding domain. AvrRxo1-ORF2 functions to suppress the bacteriostatic activity of AvrRxo1-ORF1 in bacterial cells.

Mechanism of dephosphorylation of glucosyl-3-phosphoglycerate by a histidine phosphatase

Mycobacterium tuberculosis (Mtb) synthesizes polymethylated polysaccharides that form complexes with long chain fatty acids. These complexes, referred to as methylglucose lipopolysaccharides (MGLPs), regulate fatty acid biosynthesis in vivo, including biosynthesis of mycolic acids that are essential for building the cell wall. Glucosyl-3-phosphoglycerate phosphatase (GpgP, EC 5.4.2.1), encoded by Rv2419c gene, catalyzes the second step of the pathway for the biosynthesis of MGLPs. The molecular basis for this dephosphorylation is currently not understood. Here, we describe the crystal structures of apo-, vanadate-bound, and phosphate-bound MtbGpgP, depicting unliganded, reaction intermediate mimic, and product-bound views of MtbGpgP, respectively. The enzyme consists of a single domain made up of a central β-sheet flanked by α-helices on either side. The active site is located in a positively charged cleft situated above the central β-sheet. Unambiguous electron density for vanadate covalently bound to His(11), mimicking the phosphohistidine intermediate, was observed. The role of residues interacting with the ligands in catalysis was probed by site-directed mutagenesis. Arg(10), His(11), Asn(17), Gln(23), Arg(60), Glu(84), His(159), and Leu(209) are important for enzymatic activity. Comparison of the structures of MtbGpgP revealed conformational changes in a key loop region connecting β1 with α1. This loop regulates access to the active site. MtbGpgP functions as dimer. L209E mutation resulted in monomeric GpgP, rendering the enzyme incapable of dephosphorylation. The structures of GpgP reported here are the first crystal structures for histidine-phosphatase-type GpgPs. These structures shed light on a key step in biosynthesis of MGLPs that could be targeted for development of anti-tuberculosis therapeutics.

Phosphofructo-2-kinase/fructose-2,6-bisphosphatase modulates oscillations of pancreatic islet metabolism

Pulses of insulin from pancreatic beta-cells help maintain blood glucose in a narrow range, although the source of these pulses is unclear. It has been proposed that a positive feedback circuit exists within the glycolytic pathway, the autocatalytic activation of phosphofructokinase-1 (PFK1), which endows pancreatic beta-cells with the ability to generate oscillations in metabolism. Flux through PFK1 is controlled by the bifunctional enzyme PFK2/FBPase2 (6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase) in two ways: via (1) production/degradation of fructose-2,6-bisphosphate (Fru2,6-BP), a potent allosteric activator of PFK1, as well as (2) direct activation of glucokinase due to a protein-protein interaction. In this study, we used a combination of live-cell imaging and mathematical modeling to examine the effects of inducibly-expressed PFK2/FBPase2 mutants on glucose-induced Ca²⁺ pulsatility in mouse islets. Irrespective of the ability to bind glucokinase, mutants of PFK2/FBPase2 that increased the kinase:phosphatase ratio reduced the period and amplitude of Ca²⁺ oscillations. Mutants which reduced the kinase:phosphatase ratio had the opposite effect. These results indicate that the main effect of the bifunctional enzyme on islet pulsatility is due to Fru2,6-BP alteration of the threshold for autocatalytic activation of PFK1 by Fru1,6-BP. Using computational models based on PFK1-generated islet oscillations, we then illustrated how moderate elevation of Fru-2,6-BP can increase the frequency of glycolytic oscillations while reducing their amplitude, with sufficiently high activation resulting in termination of slow oscillations. The concordance we observed between PFK2/FBPase2-induced modulation of islet oscillations and the models of PFK1-driven oscillations furthermore suggests that metabolic oscillations, like those found in yeast and skeletal muscle, are shaped early in glycolysis.

Molecular basis of the fructose-2,6-bisphosphatase reaction of PFKFB3: transition state and the C-terminal function

The molecular basis of fructose-2,6-bisphosphatase (F-2,6-P(2)ase) of 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFKFB) was investigated using the crystal structures of the human inducible form (PFKFB3) in a phospho-enzyme intermediate state (PFKFB3-P•F-6-P), in a transition state-analogous complex (PFKFB3•AlF(4)), and in a complex with pyrophosphate (PFKFB3•PP(i)) at resolutions of 2.45, 2.2, and 2.3 Å, respectively. Trapping the PFKFB3-P•F-6-P intermediate was achieved by flash cooling the crystal during the reaction, and the PFKFB3•AlF(4) and PFKFB3•PP(i) complexes were obtained by soaking. The PFKFB3•AlF(4) and PFKFB3•PP(i) complexes resulted in removing F-6-P from the catalytic pocket. With these structures, the structures of the Michaelis complex and the transition state were extrapolated. For both the PFKFB3-P formation and break down, the phosphoryl donor and the acceptor are located within ~5.1 Å, and the pivotal point 2-P is on the same line, suggesting an "in-line" transfer with a direct inversion of phosphate configuration. The geometry suggests that NE2 of His253 undergoes a nucleophilic attack to form a covalent N-P bond, breaking the 2O-P bond in the substrate. The resulting high reactivity of the leaving group, 2O of F-6-P, is neutralized by a proton donated by Glu322. Negative charges on the equatorial oxygen of the transient bipyramidal phosphorane formed during the transfer are stabilized by Arg252, His387, and Asn259. The C-terminal domain (residues 440-446) was rearranged in PFKFB3•PP(i), implying that this domain plays a critical role in binding of substrate to and release of product from the F-2,6-P(2) ase catalytic pocket. These findings provide a new insight into the understanding of the phosphoryl transfer reaction.

Structural and biochemical studies of TIGAR (TP53-induced glycolysis and apoptosis regulator)

Activation of the p53 tumor suppressor by cellular stress leads to variable responses ranging from growth inhibition to apoptosis. TIGAR is a novel p53-inducible gene that inhibits glycolysis by reducing cellular levels of fructose-2,6-bisphosphate, an activator of glycolysis and inhibitor of gluconeogenesis. Here we describe structural and biochemical studies of TIGAR from Danio rerio. The overall structure forms a histidine phosphatase fold with a phosphate molecule coordinated to the catalytic histidine residue and a second phosphate molecule in a position not observed in other phosphatases. The recombinant human and zebra fish enzymes hydrolyze fructose-2,6-bisphosphate as well as fructose-1,6-bisphosphate but not fructose 6-phosphate in vitro. The TIGAR active site is open and positively charged, consistent with its enzymatic function as bisphosphatase. The closest related structures are the bacterial broad specificity phosphatase PhoE and the fructose-2,6-bisphosphatase domain of the bifunctional 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase. The structural comparison shows that TIGAR combines an accessible active site as observed in PhoE with a charged substrate-binding pocket as seen in the fructose-2,6-bisphosphatase domain of the bifunctional enzyme.

Emerging roles for protein histidine phosphorylation in cellular signal transduction: lessons from the islet beta-cell

Protein phosphorylation represents one of the key regulatory events in physiological insulin secretion from the islet β-cell. In this context, several classes of protein kinases (e.g. calcium-, cyclic nucleotide- and phospholipid-dependent protein kinases and tyrosine kinases) have been characterized in the β-cell. The majority of phosphorylated amino acids identified include phosphoserine, phosphothreonine and phosphotyrosine. Protein histidine phosphorylation has been implicated in the prokaryotic and eukaryotic cellular signal transduction. Most notably, phoshohistidine accounts for 6% of total protein phosphorylation in eukaryotes, which makes it nearly 100-fold more abundant than phosphotyrosine, but less abundant than phosphoserine and phosphothreonine. However, very little is known about the number of proteins with phosphohistidines, since they are highly labile and are rapidly lost during phosphoamino acid identification under standard experimental conditions. The overall objectives of this review are to: (i) summarize the existing evidence indicating the subcellular distribution and characterization of various histidine kinases in the islet β-cell, (ii) describe evidence for functional regulation of these kinases by agonists of insulin secretion, (iii) present a working model to implicate novel regulatory roles for histidine kinases in the receptor-independent activation, by glucose, of G-proteins endogenous to the β-cell, (iv) summarize evidence supporting the localization of protein histidine phosphatases in the islet β-cell and (v) highlight experimental evidence suggesting potential defects in the histidine kinase signalling cascade in islets derived from the Goto-Kakizaki (GK) rat, a model for type 2 diabetes. Potential avenues for future research to further decipher regulatory roles for protein histidine phosphorylation in physiological insulin secretion are also discussed.

The histidine phosphatase superfamily: structure and function

The histidine phosphatase superfamily is a large functionally diverse group of proteins. They share a conserved catalytic core centred on a histidine which becomes phosphorylated during the course of the reaction. Although the superfamily is overwhelmingly composed of phosphatases, the earliest known and arguably best-studied member is dPGM (cofactor-dependent phosphoglycerate mutase). The superfamily contains two branches sharing very limited sequence similarity: the first containing dPGM, fructose-2,6-bisphosphatase, PhoE, SixA, TIGAR [TP53 (tumour protein 53)-induced glycolysis and apoptosis regulator], Sts-1 and many other activities, and the second, smaller, branch composed mainly of acid phosphatases and phytases. Human representatives of both branches are of considerable medical interest, and various parasites contain superfamily members whose inhibition might have therapeutic value. Additionally, several phosphatases, notably the phytases, have current or potential applications in agriculture. The present review aims to draw together what is known about structure and function in the superfamily. With the benefit of an expanding set of histidine phosphatase superfamily structures, a clearer picture of the conserved elements is obtained, along with, conversely, a view of the sometimes surprising variation in substrate-binding and proton donor residues across the superfamily. This analysis should contribute to correcting a history of over- and mis-annotation in the superfamily, but also suggests that structural knowledge, from models or experimental structures, in conjunction with experimental assays, will prove vital for the future description of function in the superfamily.

Design of fructose-2,6-bisphosphatase inhibitors: A novel virtual screening approach

Fructose-2,6-bisphosphatase (FBPase-2) is a switch between gluconeogenesis and glycolysis in the hepatic cells. The structural features required for inhibitory activity of FBPase-2 were unidentified; no leads are available for inhibiting this important enzyme. In this paper pharmacophore mapping, molecular docking methods were employed in a virtual screening strategy to identify leads for FBPase-2. A receptor based pharmacophore map was modeled which comprised of important interactions as observed in co-crystal of rat liver isozyme with the product inhibitor fructose-6-phosphate. The pharmacophore model was validated against two databases of best docked structural analogues of fructose-2,6-bisphosphate and fructose-6-phosphate. The query generated was submitted for flexible search of ligands in chemical databases, namely LeadQuest, Maybridge and NCI. The hits obtained were further screened by molecular docking using FlexX.

6-Phosphofructo-2-kinase and fructose-2,6-bisphosphatase in Trypanosomatidae. Molecular characterization, database searches, modelling studies and evolutionary analysis

Fructose 2,6-bisphosphate is a potent allosteric activator of trypanosomatid pyruvate kinase and thus represents an important regulator of energy metabolism in these protozoan parasites. A 6-phosphofructo-2-kinase, responsible for the synthesis of this regulator, was highly purified from the bloodstream form of Trypanosoma brucei and kinetically characterized. By searching trypanosomatid genome databases, four genes encoding proteins homologous to the mammalian bifunctional enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK-2/FBPase-2) were found for both T. brucei and the related parasite Leishmania major and four pairs in Trypanosoma cruzi. These genes were predicted to each encode a protein in which, at most, only a single domain would be active. Two of the T. brucei proteins showed most conservation in the PFK-2 domain, although one of them was predicted to be inactive due to substitution of residues responsible for ligating the catalytically essential divalent metal cation; the two other proteins were most conserved in the FBPase-2 domain. The two PFK-2-like proteins were expressed in Escherichia coli. Indeed, the first displayed PFK-2 activity with similar kinetic properties to that of the enzyme purified from T. brucei, whereas no activity was found for the second. Interestingly, several of the predicted trypanosomatid PFK-2/FBPase-2 proteins have long N-terminal extensions. The N-terminal domains of the two polypeptides with most similarity to mammalian PFK-2s contain a series of tandem repeat ankyrin motifs. In other proteins such motifs are known to mediate protein-protein interactions. Phylogenetic analysis suggests that the four different PFK-2/FBPase-2 isoenzymes found in Trypanosoma and Leishmania evolved from a single ancestral bifunctional enzyme within the trypanosomatid lineage. A possible explanation for the evolution of multiple monofunctional enzymes and for the presence of the ankyrin-motif repeats in the PFK-2 isoenzymes is presented.

The kinase activity of human brain 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase is regulated via inhibition by phosphoenolpyruvate

The two enzymatic activities of the highly conserved catalytic core of 6PF2K/Fru-2,6-P(2)ase are thought to be reciprocally regulated by the amino- and carboxy-terminal regions unique to each isoform. In this study, we describe the recombinant expression, purification, and kinetic characterization of two human brain 6PF2K/Fru-2,6-P(2)ase splice variants, HBP1 and HBP2. Interestingly, both lack an arginine which is highly conserved among other tissue isoforms, and which is understood to be critical to the fructose-2,6-bisphosphatase mechanism. As a result, the phosphatase activity of both HBP isoforms is negligible, but we found that it could be recovered by restoration of the arginine by site directed mutagenesis. We also found that AMP activated protein kinase and protein kinases A, B, and C catalyzed the phosphorylation of Ser-460 of HBP1, and that in addition both isoforms are phosphorylated at a second, as yet undetermined site by protein kinase C. However, none of the phosphorylations had any effect on the intrinsic kinetic characteristics of either enzymatic activity, and neither did point mutation (mimicking phosphorylation), deletion, and alternative-splice modification of the HBP1 carboxy-terminal region. Instead, these phosphorylations and mutations decreased the sensitivity of the 6PF2K to a potent allosteric inhibitor, phosphoenolpyruvate, which appears to be the major regulatory mechanism.

Insights into the catalytic mechanism of cofactor-independent phosphoglycerate mutase from X-ray crystallography, simulated dynamics and molecular modeling

Phosphoglycerate mutases catalyze the isomerization of 2 and 3-phosphoglycerates, and are essential for glucose metabolism in most organisms. Here, we further characterize the 2,3-bisphosphoglycerate-independent phosphoglycerate mutase (iPGM) from Bacillus stearothermophilus by determination of a high-resolution (1.4A) crystal structure of the wild-type enzyme and the crystal structure of its S62A mutant. The mutant structure surprisingly showed the replacement of one of the two catalytically essential manganese ions with a water molecule, offering an additional possible explanation for its lack of catalytic activity. Crystal structures invariably show substrate phosphoglycerate to be entirely buried in a deep cleft between the two iPGM domains. Flexibility analyses were therefore employed to reveal the likely route of substrate access to the catalytic site through an aperture created in the enzyme's surface during certain stages of the catalytic process. Several conserved residues lining this aperture may contribute to orientation of the substrate as it enters. Factors responsible for the retention of glycerate within the phosphoenzyme structure in the proposed mechanism are identified by molecular modeling of the glycerate complex of the phosphoenzyme. Taken together, these results allow for a better understanding of the mechanism of action of iPGMs. Many of the results are relevant to a series of evolutionarily related enzymes. These studies will facilitate the development of iPGM inhibitors which, due to the demonstrated importance of this enzyme in many bacteria, would be of great potential clinical significance.

Unexpected catalytic site variation in phosphoprotein phosphatase homologues of cofactor-dependent phosphoglycerate mutase

The cofactor-dependent phosphoglycerate mutase (dPGM) superfamily contains, besides mutases, a variety of phosphatases, both broadly and narrowly substrate-specific. Distant dPGM homologues, conspicuously abundant in microbial genomes, represent a challenge for functional annotation based on sequence comparison alone. Here we carry out sequence analysis and molecular modelling of two families of bacterial dPGM homologues, one the SixA phosphoprotein phosphatases, the other containing various proteins of no known molecular function. The models show how SixA proteins have adapted to phosphoprotein substrate and suggest that the second family may also encode phosphoprotein phosphatases. Unexpected variation in catalytic and substrate-binding residues is observed in the models.

Structures of phosphate and trivanadate complexes of Bacillus stearothermophilus phosphatase PhoE: structural and functional analysis in the cofactor-dependent phosphoglycerate mutase superfamily

Bacillus stearothermophilus phosphatase PhoE is a member of the cofactor-dependent phosphoglycerate mutase superfamily possessing broad specificity phosphatase activity. Its previous structural determination in complex with glycerol revealed probable bases for its efficient hydrolysis of both large, hydrophobic, and smaller, hydrophilic substrates. Here we report two further structures of PhoE complexes, to higher resolution of diffraction, which yield a better and thorough understanding of its catalytic mechanism. The environment of the phosphate ion in the catalytic site of the first complex strongly suggests an acid-base catalytic function for Glu83. It also reveals how the C-terminal tail ordering is linked to enzyme activation on phosphate binding by a different mechanism to that seen in Escherichia coli phosphoglycerate mutase. The second complex structure with an unusual doubly covalently bound trivanadate shows how covalent modification of the phosphorylable His10 is accompanied by small structural changes, presumably to catalytic advantage. When compared with structures of related proteins in the cofactor-dependent phosphoglycerate mutase superfamily, an additional phosphate ligand, Gln22, is observed in PhoE. Functional constraints lead to the corresponding residue being conserved as Gly in fructose-2,6-bisphosphatases and Thr/Ser/Cys in phosphoglycerate mutases. A number of sequence annotation errors in databases are highlighted by this analysis. B. stearothermophilus PhoE is evolutionarily related to a group of enzymes primarily present in Gram-positive bacilli. Even within this group substrate specificity is clearly variable highlighting the difficulties of computational functional annotation in the cofactor-dependent phosphoglycerate mutase superfamily.

Tissue-specific structure/function differentiation of the liver isoform of 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase

The crystal structures of the human liver 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase in three different liganding states were determined and compared with those of the rat testis isozyme. A set of amino acid sequence heterogeneity from the two distinct genes encoding the two different tissue isozymes leads to both global and local conformational differences that may cause the differences in catalytic properties of the two isozymes. The sequence differences in a beta-hairpin loop in the kinase domain causes a translational shift of several hydrophobic interactions in the dimeric contact region, and its propagation to the domains interface results in a 5 degrees twist of the entire bisphosphatase domain relative to the kinase domain. The bisphosphatase domain twist allows the dimeric interactions between the bisphosphatase domains, which are negligible in the testis enzyme, and as a result, the conformational stability of the domain is increased. Sequence polymorphisms also confer small but significant structural dissimilarities in the substrate-binding loops, allowing the differentiated catalytic properties between the two different tissue-type isozymes. Whereas the polymorphic sequence at the bisphosphatase-active pocket suggests a more suitable substrate binding, a similar extent of sequence differences at the kinase-active pocket confers a different mechanism of substrates bindings to the kinase-active pocket. It includes the ATP-sensitive unwinding of the switch helix alpha5, which is a characteristic ATP-dependent conformational change in the testis form. The sequence-dependent structural difference disallows the liver kinase to follow the ATP-switch mechanism. Altogether these suggest that the liver isoform has structural features more appropriate for an elevated bisphosphatase activity, compared with that of the testis form. The structural predisposition for bisphosphatase activity in the liver isozyme is consistent with the liver-unique glucose metabolic pathway, gluconeogenesis.

Critical switch of the metabolic fluxes by phosphofructo-2-kinase:fructose-2,6-bisphosphatase. A kinetic model

A kinetic model for the bifunctional enzyme, phosphofructo-2-kinase:fructose-2,6-bisphosphatase, is analysed by application of the graph-theoretical method, considering comparable levels for all participants. Certain elementary reactions, distributed on the enzyme surface, are considered to be co-ordinated in a single conformational transition (a model of parallel molecular operations). The method allows us to identify in the kinetic scheme its destabilising sub-scheme as a branched cycle of elementary reactions. Under certain conditions this sub-scheme induces critical phenomena (bistability or oscillations). The computer calculations for the estimated parameter values fit well the experimental observations for this system. The model explains the periodic or bistable counterphase changes of the two opposing activities of this enzyme, observed after glucose perfusion of rat hepatic enzyme samples, and predicts drastic critical changes in kinetic behaviour induced by small external signals. The model also shows the necessity of the phosphoryl intermediate in the mechanism of the bisphosphatase for the critical kind of kinetic behaviour.

Structure and mechanism of action of a cofactor-dependent phosphoglycerate mutase homolog from Bacillus stearothermophilus with broad specificity phosphatase activity

The crystal structure of Bacillus stearothermophilus PhoE (originally termed YhfR), a broad specificity monomeric phosphatase with a molecular mass of approximately 24 kDa, has been solved at 2.3 A resolution in order to investigate its structure and function. PhoE, already identified as a homolog of a cofactor-dependent phosphoglycerate mutase, shares with the latter an alpha/beta/alpha sandwich structure spanning, as a structural excursion, a smaller subdomain composed of two alpha-helices and one short beta-strand. The active site contains residues from both the alpha/beta/alpha sandwich and the sub-domain. With the exception of the hydrophilic catalytic machinery conserved throughout the cofactor-dependent phosphoglycerate mutase family, the active-site cleft is strikingly hydrophobic. Docking studies with two diverse, favored substrates show that 3-phosphoglycerate may bind to the catalytic core, while alpha-napthylphosphate binding also involves the hydrophobic portion of the active-site cleft. Combining a highly favorable phospho group binding site common to these substrate binding modes and data from related enzymes, a catalytic mechanism can be proposed that involves formation of a phosphohistidine intermediate on His10 and likely acid-base behavior of Glu83. Other structural factors contributing to the broad substrate specificity of PhoE can be identified. The dynamic independence of the subdomain may enable the active-site cleft to accommodate substrates of different sizes, although similar motions are present in simulations of cofactor-dependent phosphoglycerate mutases, perhaps favoring a more general functional role. A significant number of entries in protein sequence databases, particularly from unfinished microbial genomes, are more similar to PhoE than to cofactor-dependent phosphoglycerate mutases or to fructose-2,6-bisphosphatases. This PhoE structure will therefore serve as a valuable basis for inference of structural and functional characteristics of these proteins.

A cofactor-dependent phosphoglycerate mutase homolog from Bacillus stearothermophilus is actually a broad specificity phosphatase

The distribution of phosphoglycerate mutase (PGM) activity in bacteria is complex, with some organisms possessing both a cofactor-dependent and a cofactor-independent PGM and others having only one of these enzymes. Although Bacillus species contain only a cofactor-independent PGM, genes homologous to those encoding cofactor-dependent PGMs have been detected in this group of bacteria, but in at least one case the encoded protein lacks significant PGM activity. Here we apply sequence analysis, molecular modeling, and enzymatic assays to the cofactor-dependent PGM homologs from B. stearothermophilus and B. subtilis, and show that these enzymes are phosphatases with broad substrate specificity. Homologs from other gram-positive bacteria are also likely to possess phosphatase activity. These studies clearly show that the exploration of genomic sequences through three-dimensional modeling is capable of producing useful predictions regarding function. However, significant methodological improvements will be needed before such analysis can be carried out automatically.

Enzyme catalysis: removing chemically 'essential' residues by site-directed mutagenesis

Enzymatic catalysis relies on the action of the amino acid side chains arrayed in the enzyme active sites. Usually, only two or three 'essential' residues are directly involved in the bond making and breaking steps leading to product formation. For the past 20 years, enzymologists have been addressing the role of such residues by changing them into chemically inert side chains. Removal of an 'essential' group often does not abolish activity, but can significantly alter the catalytic mechanism. Such results underscore the sophistication of enzyme catalysis and the functional plasticity of enzyme active sites.

Binding of ATP to the fructose-2,6-bisphosphatase domain of chicken liver 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase leads to activation of its 6-phosphofructo-2-kinase

To understand the mechanism by which the activity of the 6-phosphofructo-2-kinase (6PF-2K) of chicken liver 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase is stimulated by its substrate ATP, we studied two mutants of the enzyme. Mutation of either Arg-279, the penultimate basic residue within the Walker A nucleotide-binding fold in the bisphosphatase domain, or Arg-359 to Ala eliminated the activation of the chicken 6PF-2K by ATP. Binding analysis by fluorescence spectroscopy using 2'(3')-O-(N-methylanthraniloyl)-ATP revealed that the kinase domains of these two mutants, unlike that of the wild type enzyme, showed no cooperativity in ATP binding and that the mutant enzymes possess only the high affinity ATP binding site, suggesting that the ATP binding site on the bisphosphatase domain represents the low affinity site. This conclusion was supported by the result that the affinity of ATP for the isolated bisphosphatase domain is similar to that for the low affinity site in the wild type enzyme. In addition, we found that the 6PF-2K of a chimeric enzyme, in which the last 25 residues of chicken enzyme were replaced with those of the rat enzyme, could not be activated by ATP, despite the fact that the ATP-binding properties of this chimeric enzyme were not different from those of the wild type chicken enzyme. These results demonstrate that activation of the chicken 6PF-2K by ATP may result from allosteric binding of ATP to the bisphosphatase domain where residues Arg-279 and Arg-359 are critically involved and require specific C-terminal sequences.

Characterization of active-site mutants of Schizosaccharomyces pombe phosphoglycerate mutase. Elucidation of the roles of amino acids involved in substrate binding and catalysis

The roles of a number of amino acids present at the active site of the monomeric phosphoglycerate mutase from the fission yeast Schizosaccharomyces pombe have been explored by site-directed mutagenesis. The amino acids examined could be divided broadly into those presumed from previous related structural studies to be important in the catalytic process (R14, S62 and E93) and those thought to be important in substrate binding (R94, R120 and R121). Most of these residues have not previously been studied by site-directed mutagenesis. All the mutants except R14 were expressed in an engineered null strain of Saccharomyces cerevisiae (S150-gpm:HIS) in good yield. The R14Q mutant was expressed in good yield in the transformed AH22 strain of S. cerevisiae. The S62A mutant was markedly unstable, preventing purification. The various mutants were purified to homogeneity and characterized in terms of kinetic parameters, CD and fluorescence spectra, stability towards denaturation by guanidinium chloride, and stability of phosphorylated enzyme intermediate. In addition, the binding of substrate (3-phosphoglycerate) to wild-type, E93D and R120,121Q enzymes was measured by isothermal titration calorimetry. The results provide evidence for the proposed roles of each of these amino acids in the catalytic cycle and in substrate binding, and will support the current investigation of the structure and dynamics of the enzyme using multidimensional NMR techniques.

Reaction mechanism of fructose-2,6-bisphosphatase. A mutation of nucleophilic catalyst, histidine 256, induces an alteration in the reaction pathway

A bifunctional enzyme, fructose-6-phosphate,2-kinase/fructose 2, 6-bisphosphatase (Fru-6-P,2-kinase/Fru-2,6-Pase), catalyzes synthesis and degradation of fructose 2,6-bisphosphate (Fru-2,6-P2). Previously, the rat liver Fru-2,6-Pase reaction (Fru-2,6-P2 --> Fru-6-P + Pi) has been shown to proceed via a phosphoenzyme intermediate with His258 phosphorylated, and mutation of the histidine to alanine resulted in complete loss of activity (Tauler, A., Lin, K., and Pilkis, S. J. (1990) J. Biol. Chem. 265, 15617-15622). In the present study, it is shown that mutation of the corresponding histidine (His256) of the rat testis enzyme decreases activity by less than a factor of 10 with a kcat of 17% compared with the wild type enzyme. Mutation of His390 (in close proximity to His256) to Ala results in a kcat of 12.5% compared with the wild type enzyme. Attempts to detect a phosphohistidine intermediate with the H256A mutant enzyme were unsuccessful, but the phosphoenzyme is detected in the wild type, H390A, R255A, R305S, and E325A mutant enzymes. Data demonstrate that the mutation of His256 induces a change in the phosphatase hydrolytic reaction mechanism. Elimination of the nucleophilic catalyst, H256A, results in a change in mechanism. In the H256A mutant enzyme, His390 likely acts as a general base to activate water for direct hydrolysis of the 2-phosphate of Fru-2,6-P2. Mutation of Arg255 and Arg305 suggests that the arginines probably have a role in neutralizing excess charge on the 2-phosphate and polarizing the phosphoryl for subsequent transfer to either His256 or water. The role of Glu325 is less certain, but it may serve as a general acid, protonating the leaving 2-hydroxyl of Fru-2,6-P2.