Communication between the nucleotide site and the main molecular hinge of 3-phosphoglycerate kinase

The Phosphoglycerate Kinase (PGK) Gene Family of Maize (Zea mays var. B73)

Phosphoglycerate kinase (PGK, E.C. 2.7.2.3) interconverts ADP + 1,3-bisphospho-glycerate (1,3-bPGA) to ATP + 3-phosphoglycerate (3PGA). While most bacteria have a single pgk gene and mammals possess two copies, plant genomes contain three or more PGK genes. In this study, we identified five Pgk genes in the Zea mays var. B73 genome, predicted to encode proteins targeted to different subcellular compartments: ZmPgk1, ZmPgk2, and ZmPgk4 (chloroplast), ZmPgk3 (cytosol), and ZmPgk5 (nucleus). The expression of ZmPgk3 was highest in non-photosynthetic tissues (roots and cobs), where PGK activity was also greatest, consistent with a function in glycolysis. Green tissues (leaf blade and husk leaf) showed intermediate levels of PGK activity, and predominantly expressed ZmPgk1 and ZmPgk2, suggesting involvement in photosynthetic metabolism. ZmPgk5 was weakly expressed and ZmPgk4 was not detected in any tissue. Phylogenetic analysis showed that the photosynthetic and glycolytic isozymes of plants clustered together, but were distinct from PGKs of animals, fungi, protozoa, and bacteria, indicating that photosynthetic and glycolytic isozymes of plants diversified after the divergence of the plant lineage from other groups. These results show the distinct role of each PGK in maize and provide the basis for future studies into the regulation and function of this key enzyme.

Weak protein-protein interactions in live cells are quantified by cell-volume modulation

Weakly bound protein complexes play a crucial role in metabolic, regulatory, and signaling pathways, due in part to the high tunability of their bound and unbound populations. This tunability makes weak binding (micromolar to millimolar dissociation constants) difficult to quantify under biologically relevant conditions. Here, we use rapid perturbation of cell volume to modulate the concentration of weakly bound protein complexes, allowing us to detect their dissociation constant and stoichiometry directly inside the cell. We control cell volume by modulating media osmotic pressure and observe the resulting complex association and dissociation by FRET microscopy. We quantitatively examine the interaction between GAPDH and PGK, two sequential enzymes in the glycolysis catalytic cycle. GAPDH and PGK have been shown to interact weakly, but the interaction has not been quantified in vivo. A quantitative model fits our experimental results with log K_d = -9.7 ± 0.3 and a 2:1 prevalent stoichiometry of the GAPDH:PGK complex. Cellular volume perturbation is a widely applicable tool to detect transient protein interactions and other biomolecular interactions in situ. Our results also suggest that cells could use volume change (e.g., as occurs upon entry to mitosis) to regulate function by altering biomolecular complex concentrations.

MgF3- and AlF4- transition state analogue complexes of yeast phosphoglycerate kinase

The phospho-transfer mechanism of yeast phosphoglycerate kinase (PGK) has been probed through formation of trifluoromagnesate (MgF₃^-) and tetrafluoroaluminate (AlF₄^-) transition state analogue complexes and analyzed using ¹⁹F, ¹H waterLOGSY and ¹H chemical shift perturbation NMR spectroscopy. We observed the first ¹⁹F NMR spectroscopic evidence for the formation of metal fluoride transition state analogues of yeast PGK and also observed significant changes to proton chemical shifts of PGK in the presence, but not in the absence, of fluoride upon titration of ligands, providing indirect evidence of the formation of a closed ternary transition state. WaterLOGSY NMR spectroscopy experiments using an uncompetitive model were used in an attempt to measure ligand binding affinities within the transition state analogue complexes.

Characterization of 3-phosphoglycerate kinase from Corynebacterium glutamicum and its impact on amino acid production

BACKGROUND

Corynebacterium glutamicum cg1790/pgk encodes an enzyme active as a 3-phosphoglycerate kinase (PGK) (EC 2.7.2.3) catalyzing phosphoryl transfer from 1,3-biphosphoglycerate (bPG) to ADP to yield 3-phosphoglycerate (3-PG) and ATP in substrate chain phosphorylation.

RESULTS

C. glutamicum 3-phosphoglycerate kinase was purified to homogeneity from the soluble fraction of recombinant E. coli. PGK(His) was found to be active as a homodimer with molecular weight of 104 kDa. The enzyme preferred conditions of pH 7.0 to 7.4 and required Mg²⁺ for its activity. PGK(His) is thermo labile and it has shown maximal activity at 50-65°C. The maximal activity of PGK(His) was estimated to be 220 and 150 U mg-1 with KM values of 0.26 and 0.11 mM for 3-phosphoglycerate and ATP, respectively. A 3-phosphoglycerate kinase negative C. glutamicum strain ∆pgk was constructed and shown to lack the ability to grow under glycolytic or gluconeogenic conditions unless PGK was expressed from a plasmid to restore growth. When pgk was overexpressed in L-arginine and L-ornithine production strains the production increased by 8% and by 17.5%, respectively.

CONCLUSION

Unlike many bacterial PGKs, C. glutamicum PGK is active as a homodimer. PGK is essential for growth of C. glutamicum with carbon sources requiring glycolysis and gluconeogenesis. Competitive inhibition by ADP reveals the critical role of PGK in gluconeogenesis by energy charge. Pgk overexpression improved the productivity in L-arginine and L-ornithine production strains.

An allosteric signaling pathway of human 3-phosphoglycerate kinase from force distribution analysis

3-Phosphogycerate kinase (PGK) is a two domain enzyme, which transfers a phosphate group between its two substrates, 1,3-bisphosphoglycerate bound to the N-domain and ADP bound to the C-domain. Indispensable for the phosphoryl transfer reaction is a large conformational change from an inactive open to an active closed conformation via a hinge motion that should bring substrates into close proximity. The allosteric pathway resulting in the active closed conformation has only been partially uncovered. Using Molecular Dynamics simulations combined with Force Distribution Analysis (FDA), we describe an allosteric pathway, which connects the substrate binding sites to the interdomain hinge region. Glu192 of alpha-helix 7 and Gly394 of loop L14 act as hinge points, at which these two secondary structure elements straighten, thereby moving the substrate-binding domains towards each other. The long-range allosteric pathway regulating hPGK catalytic activity, which is partially validated and can be further tested by mutagenesis, highlights the virtue of monitoring internal forces to reveal signal propagation, even if only minor conformational distortions, such as helix bending, initiate the large functional rearrangement of the macromolecule.

3-Phosphoglycerate kinase (PGK) is an essential enzyme for living organisms. It catalyzes the phospho-transfer reaction between two catabolites during carbohydrate metabolism. In addition to this physiological role, human PGK has been shown to phosphorylate L-nucleoside analogues, potential drugs against viral infection and cancer. PGK is a two domain enzyme, with the two substrates bound to the two separate domains. In order to perform its function the enzyme has to undergo a large conformational change involving a hinge bending to bring the substrates into close proximity. The allosteric pathway from the open non-reactive state of PGK to the closed reactive state as triggered by substrate binding has only been partially uncovered by experimental studies. Here we describe a complete allosteric pathway, which connects the substrate binding sites to the interdomain hinge region using Molecular Dynamics simulations combined with Force Distribution Analysis (FDA). While previously identified key residues involved in PGK domain closure are part of this pathway, we here fill the numerous gaps in the pathway by identifying newly uncovered residues and interesting candidates for future mutational studies.

Analysis of Chinese hamster ovary cell metabolism through a combined computational and experimental approach

Optimization of cell culture processes can benefit from the systematic analysis of experimental data and their organization in mathematical models, which can be used to decipher the effect of individual process variables on multiple outputs of interest. Towards this goal, a kinetic model of cytosolic glucose metabolism coupled with a population-level model of Chinese hamster ovary cells was used to analyse metabolic behavior under batch and fed-batch cell culture conditions. The model was parameterized using experimental data for cell growth dynamics, extracellular and intracellular metabolite profiles. The results highlight significant differences between the two culture conditions in terms of metabolic efficiency and motivate the exploration of lactate as a secondary carbon source. Finally, the application of global sensitivity analysis to the model parameters highlights the need for additional experimental information on cell cycle distribution to complement metabolomic analyses with a view to parameterize kinetic models.

Structural characterization of the ribonuclease H-like type ASKHA superfamily kinase MK0840 from Methanopyrus kandleri

Murein recycling is a process in which microorganisms recover peptidoglycan-degradation products in order to utilize them in cell wall biosynthesis or basic metabolic pathways. Methanogens such as Methanopyrus kandleri contain pseudomurein, which differs from bacterial murein in its composition and branching. Here, four crystal structures of the putative sugar kinase MK0840 from M. kandleri in apo and nucleotide-bound states are reported. MK0840 shows high similarity to bacterial anhydro-N-acetylmuramic acid kinase, which is involved in murein recycling. The structure shares a common fold with panthothenate kinase and the 2-hydroxyglutaryl-CoA dehydratase component A, both of which are members of the ASKHA (acetate and sugar kinases/Hsc70/actin) superfamily of phosphotransferases. Local conformational changes in the nucleotide-binding site between the apo and holo forms are observed upon nucleotide binding. Further insight is given into domain movements and putative active-site residues are identified.

Mechanisms of protein sequence divergence and incompatibility

Alignments of orthologous protein sequences convey a complex picture. Some positions are utterly conserved whilst others have diverged to variable degrees. Amongst the latter, many are non-exchangeable between extant sequences. How do functionally critical and highly conserved residues diverge? Why and how did these exchanges become incompatible within contemporary sequences? Our model is phosphoglycerate kinase (PGK), where lysine 219 is an essential active-site residue completely conserved throughout Eukaryota and Bacteria, and serine is found only in archaeal PGKs. Contemporary sequences tested exhibited complete loss of function upon exchanges at 219. However, a directed evolution experiment revealed that two mutations were sufficient for human PGK to become functional with serine at position 219. These two mutations made position 219 permissive not only for serine and lysine, but also to a range of other amino acids seen in archaeal PGKs. The identified trajectories that enabled exchanges at 219 show marked sign epistasis - a relatively small loss of function with respect to one amino acid (lysine) versus a large gain with another (serine, and other amino acids). Our findings support the view that, as theoretically described, the trajectories underlining the divergence of critical positions are dominated by sign epistatic interactions. Such trajectories are an outcome of rare mutational combinations. Nonetheless, as suggested by the laboratory enabled K219S exchange, given enough time and variability in selection levels, even utterly conserved and functionally essential residues may change.

Orthologs are proteins in different species sharing the same function and structure. However, the mechanisms that underline the divergence of different sequences from a single ancestor remain unclear, particularly because many amino acid exchanges between orthologs result in loss of function (incompatibility). We aimed at disentangling an ancient divergence event within the active-site of a universally spread enzyme that mediates ATP synthesis. Using laboratory evolution experiments, we found that an exchange in a functionally critical active-site residue that is incompatible within contemporary orthologs is enabled by few mutations. These mutations lead to transition sequences in which, unlike the extant sequences, a wide range of amino acids is tolerated. Our experiment reveals the properties of these transition sequences that may resemble the historical ancestral states that underlined this divergence event, and the mechanisms that led to incompatibility within the contemporary orthologs. Our results support theoretical predictions and reshape our understanding of protein structure-function. That a given position is entirely conserved and essential for function does not indicate that it will never exchange, but rather, that the exchange may depend on changes in many other positions.

Importance of aspartate residues in balancing the flexibility and fine-tuning the catalysis of human 3-phosphoglycerate kinase

The exact role of the metal ion, usually Mg(2+), in the catalysis of human 3-phosphoglycerate kinase, a well-studied two-domain enzyme, has not been clarified. Here we have prepared single and double alanine mutants of the potential metal-binding residues, D374 and D218. While all mutations weaken the catalytic interactions with Mg(2+), they surprisingly strengthen binding of both MgADP and MgATP, and the effects are even more pronounced for ADP and ATP. Thermodynamic parameters of binding indicate an increase in the binding entropy as a reason for the strengthening. In agreement with the experimental results, computer-simulated annealing calculations for the complexes of these mutants have supported the mobility of the nucleotide phosphates and, as a consequence, formation of their new interaction(s) within the active site. A similar type of mobility is suggested to be a characteristic feature of the nucleotide site of the wild-type enzyme, too, both in its inactive open conformation and in the active closed conformation. This mobility of the nucleotide phosphates that is regulated by the aspartate side chains of D218 and D374 through the complexing Mg(2+) is suggested to be essential in enzyme function.

Selectivity of kinases on the activation of tenofovir, an anti-HIV agent

Nucleoside analogues, used in HIV-therapy, need to be phosphorylated by cellular enzymes in order to become potential substrates for HIV reverse transcriptase. After incorporation into the viral DNA chain, because of lacking of their 3'-hydroxyl groups, they stop the elongation process and lead to the death of the virus. Phosphorylation of the HIV-drug derivative, tenofovir monophosphate was tested with the recombinant mammalian nucleoside diphosphate kinase (NDPK), 3-phosphoglycerate kinase (PGK), creatine kinase (CK) and pyruvate kinase (PK). Among them, only CK was found to phosphorylate tenofovir monophosphate with a reasonable rate (about 45-fold lower than with its natural substrate, ADP), while PK exhibits even lower, but still detectable activity (about 1000-fold lower compared to the value with ADP). On the other hand, neither NDPK nor PGK has any detectable activity on tenofovir monophosphate. The absence of activity with PGK is surprising, since the drug tenofovir competitively inhibits both CK and PGK towards their nucleotide substrates, with similar inhibitory constants, K(I) of 2.9 and 4.8mM, respectively. Computer modelling (docking) of tenofovir mono- or diphosphate forms to these four kinases suggests that the requirement of large-scale domain closure for functioning (as for PGK) may largely restrict their applicability for phosphorylation/activation of pro-drugs having a structure similar to tenofovir monophosphate.

A new variant of phosphoglycerate kinase deficiency (p.I371K) with multiple tissue involvement: molecular and functional characterization

Phosphoglycerate kinase (PGK) is a key glycolytic enzyme that catalyzes the reversible phosphotransfer reaction from 1,3-bisphosphoglycerate to MgADP, to form 3-phosphoglycerate and MgATP. Two isozymes encoded by distinct genes are present in humans: PGK-1, located on Xq-13.3, encodes a ubiquitous protein of 417 amino acids, whereas PGK-2 is testis-specific. PGK1 deficiency is characterized by mild to severe hemolytic anemia, neurological dysfunctions and myopathy; patients rarely exhibit all three clinical features. Nearly 40 cases have been reported, 27 of them characterized at DNA or protein level, and 20 different mutations were described. Here we report the first Italian case of PGK deficiency characterized at a molecular and biochemical level. The patient presented during infancy with hemolytic anemia, increased CPK values, and respiratory distress; the study of red blood cell enzymes showed a drastic reduction in PGK activity. In adulthood he displayed mild hemolytic anemia, mental retardation and severe myopathy. PGK-1 gene sequencing revealed the new missense mutation c.1112T>A (p.Ile371Lys). The mutation was not found among 100 normal alleles, and even if located in the third to the last nucleotide of exon 9, it did not alter mRNA splicing. The p.Ile371Lys mutation falls in a conserved region of the enzyme, near the nucleotide binding site. The mutant enzyme shows reduced catalytic rates toward both substrates (apparent k(cat) values, 12-fold lower than wild-type) and a decreased affinity toward MgATP (apparent K(m), 6-fold higher than wild-type). Moreover, it lost half of activity after nearly 9-min incubation at 45°C, a temperature that did not affect the wild-type enzyme (t(1/2)>1 h). The possible compensatory expression of PGK2 isoenzyme was investigated in the proband and in the heterozygote healthy sisters, and found to be absent. Therefore, the highly perturbed catalytic properties of the new variant p.Ile371Lys, combined with protein instability, account for the PGK deficiency found in the patient and correlate with the clinical expression of the disease.

Ligand chirality effects on the dynamics of human 3-phosphoglycerate kinase: comparison between D- and L-nucleotides

l-nucleoside analogues are now largely used as antiviral drugs for the treatment of viral infections like HBV, HCV and HIV. However, in order to be fully active, they need to be phosphorylated by cellular or viral kinases. Human 3-phosphogycerate kinase (hPGK) was shown to catalyze the last step of activation of l-enantiomers and thus constitutes an attractive target for theoretical predictions of its phosphorylation efficiency. Molecular dynamics simulations were carried out with four different nucleotides (d-/l-ADP and d-/l-CDP) in complex with hPGK and 1,3-bisphospho-d-glycerate (bPG). The binding affinities of CDPs (both enantiomers) for hPGK were found very weak while d- and l-ADP were better substrates. Interestingly, the binding affinity of the bPG substrate was found to be lower in presence of d-ADP than l-ADP which indicates a potential antagonistic effect on one substrate to the other. A detailed analysis of the simulations unravels important dynamic conditions for efficient phosphorylation. Indeed, as previously described for the natural substrate, the hinge bending motion of the domains upon substrates binding should be more correlated and directional. Interestingly, the unforeseen finding was the larger dynamics freedoms observed for the substrates that was favored by the protein atoms flexibility around the nucleobase binding site.

A spring-loaded release mechanism regulates domain movement and catalysis in phosphoglycerate kinase

Phosphoglycerate kinase (PGK) is the enzyme responsible for the first ATP-generating step of glycolysis and has been implicated extensively in oncogenesis and its development. Solution small angle x-ray scattering (SAXS) data, in combination with crystal structures of the enzyme in complex with substrate and product analogues, reveal a new conformation for the resting state of the enzyme and demonstrate the role of substrate binding in the preparation of the enzyme for domain closure. Comparison of the x-ray scattering curves of the enzyme in different states with crystal structures has allowed the complete reaction cycle to be resolved both structurally and temporally. The enzyme appears to spend most of its time in a fully open conformation with short periods of closure and catalysis, thereby allowing the rapid diffusion of substrates and products in and out of the binding sites. Analysis of the open apoenzyme structure, defined through deformable elastic network refinement against the SAXS data, suggests that interactions in a mostly buried hydrophobic region may favor the open conformation. This patch is exposed on domain closure, making the open conformation more thermodynamically stable. Ionic interactions act to maintain the closed conformation to allow catalysis. The short time PGK spends in the closed conformation and its strong tendency to rest in an open conformation imply a spring-loaded release mechanism to regulate domain movement, catalysis, and efficient product release.

Nucleotide promiscuity of 3-phosphoglycerate kinase is in focus: implications for the design of better anti-HIV analogues

The wide specificity of 3-phosphoglycerate kinase (PGK) towards its nucleotide substrate is a property that allows contribution of this enzyme to the effective phosphorylation (i.e. activation) of nucleotide-based pro-drugs against HIV. Here, the structural basis of the nucleotide-PGK interaction is characterised in comparison to other kinases, namely pyruvate kinase (PK) and creatine kinase (CK), by enzyme kinetic analysis and structural modelling (docking) studies. The results provided evidence for favouring the purine vs. pyrimidine base containing nucleotides for PGK rather than for PK or CK. This is due to the exceptional ability of PGK in forming the hydrophobic contacts of the nucleotide rings that assures the appropriate positioning of the connected phosphate-chain for catalysis. As for the D-/L-configurations of the nucleotides, the L-forms (both purine and pyrimidine) are well accepted by PGK rather than either by PK or CK. Here again the dominance of the hydrophobic interactions of the L-form of pyrimidines with PGK is underlined in comparison with those of PK or CK. Furthermore, for the l-forms, the absence of the ribose OH-groups with PGK is better tolerated for the purine than for the pyrimidine containing compounds. On the other hand, the positioning of the phosphate-chain is an even more important term for PGK in the case of both purines and pyrimidines with an L-configuration, as deduced from the present kinetic studies with various nucleotide-site mutants of PGK. These characteristics of the kinase-nucleotide interactions can provide a guideline for designing new drugs.

Transition state analogue structures of human phosphoglycerate kinase establish the importance of charge balance in catalysis

Transition state analogue (TSA) complexes formed by phosphoglycerate kinase (PGK) have been used to test the hypothesis that balancing of charge within the transition state dominates enzyme-catalyzed phosphoryl transfer. High-resolution structures of trifluoromagnesate (MgF(3)(-)) and tetrafluoroaluminate (AlF(4)(-)) complexes of PGK have been determined using X-ray crystallography and (19)F-based NMR methods, revealing the nature of the catalytically relevant state of this archetypal metabolic kinase. Importantly, the side chain of K219, which coordinates the alpha-phosphate group in previous ground state structures, is sequestered into coordinating the metal fluoride, thereby creating a charge environment complementary to the transferring phosphoryl group. In line with the dominance of charge balance in transition state organization, the substitution K219A induces a corresponding reduction in charge in the bound aluminum fluoride species, which changes to a trifluoroaluminate (AlF(3)(0)) complex. The AlF(3)(0) moiety retains the octahedral geometry observed within AlF(4)(-) TSA complexes, which endorses the proposal that some of the widely reported trigonal AlF(3)(0) complexes of phosphoryl transfer enzymes may have been misassigned and in reality contain MgF(3)(-).

Functionally compensating coevolving positions are neither homoplasic nor conserved in clades

We demonstrated that a pair of positions in phosphoglycerate kinase that score highly by three nonparametric covariation measures are important for function even though the positions can be occupied by aliphatic, aromatic, or charged residues. Examination of these pairs suggested that the majority of the covariation scores could be explained by within-clade conservation. However, an analysis of diversity showed that the conservation within clades of covarying pairs was indistinguishable from pairs of positions that do not covary, thus ruling out both clade conservation and extensive homoplasy as means to identify covarying positions. Mutagenesis showed that the residues in the covarying pair were epistatic, with the type of epistasis being dependent on the initial pair. The results show that nonconserved covarying positions that affect protein function can be identified with high precision.

Thermodynamic analysis of substrate induced domain closure of 3-phosphoglycerate kinase

The energetic changes accompanying domain closure of 3-phosphoglycerate kinase, a typical hinge-bending enzyme, were assessed. Calorimetric titrations of the enzyme with each substrate, both in the absence and presence of the other one, provide information not only about the energetics of substrate binding, but of the associated conformational changes, including domain closure. Our results suggest that conformational rearrangements in the hinge generated by binding of both substrates provide the main driving force for domain closure overcoming the slightly unfavourable contact interactions between the domains.

Myopathic form of phosphoglycerate kinase (PGK) deficiency: a new case and pathogenic considerations

We describe an 18-year-old man with muscle cramps and recurrent exertional myoglobinuria, without hemolytic anemia or brain dysfunction. Phosphoglycerate kinase (PGK) deficiency was documented in muscle and erythrocytes and molecular analysis of the PGK1 gene identified a novel mutation, T378P. This is the ninth case presenting with isolated myopathy, whereas most other patients show hereditary non-spherocytic hemolytic anemia alone or associated with brain dysfunction, and a few patients have myopathy plus brain involvement. Although the diverse tissue involvement in PGK deficiency remains unclear, all mutations in myopathic patients tend to cluster in the C terminal domain, adjacent to the substrate-binding pocket. This may lead to a failure in the closure of the N terminal and C terminal domains and loss of stability due to lack of inter-domain communication during the catalytic process.

Molecular basis for the lack of enantioselectivity of human 3-phosphoglycerate kinase

Non-natural l-nucleoside analogues are increasingly used as therapeutic agents to treat cancer and viral infections. To be active, l-nucleosides need to be phosphorylated to their respective triphosphate metabolites. This stepwise phosphorylation relies on human enzymes capable of processing l-nucleoside enantiomers. We used crystallographic analysis to reveal the molecular basis for the low enantioselectivity and the broad specificity of human 3-phosphoglycerate kinase (hPGK), an enzyme responsible for the last step of phosphorylation of many nucleotide derivatives. Based on structures of hPGK in the absence of nucleotides, and bound to l and d forms of MgADP and MgCDP, we show that a non-specific hydrophobic clamp to the nucleotide base, as well as a water-filled cavity behind it, allows high flexibility in the interaction between PGK and the bases. This, combined with the dispensability of hydrogen bonds to the sugar moiety, and ionic interactions with the phosphate groups, results in the positioning of different nucleotides so to expose their diphosphate group in a position competent for catalysis. Since the third phosphorylation step is often rate limiting, our results are expected to alleviate in silico tailoring of l-type prodrugs to assure their efficient metabolic processing.