Substrate binding modifies the hinge bending characteristics of human 3-phosphoglycerate kinase: A molecular dynamics study

Temperature-jump microscopy and interaction of Hsp70 heat shock protein with a client protein in vivo

Biomolecular processes such as protein-protein interactions can depend strongly on cell type and even vary within a single cell type. Here we develop a microscope with a Peltier-controlled temperature stage, a laser temperature jump to induce heat stress, and an autofocusing feature to mitigate temperature drift during experiments, to study a protein-protein interaction in a selected cell type within a live organism, the zebrafish larva. As an application of the instrument, we show that there is considerable cell-to-cell variation of the heat shock protein Hsp70 binding to one of its clients, phosphoglycerate kinase in vivo. We adapt a key feature from our previous folding study, rare transformation of cells within the larva, so that individual cells can be imaged and differentiated for cell-to-cell response. Our approach can be extended to other organisms and cell types than the ones demonstrated in this work.

Novel inhibitors targeting the PGK1 metabolic enzyme in glycolysis exhibit effective antitumor activity against kidney renal clear cell carcinoma in vitro and in vivo

Our previous research has revealed phosphoglycerate kinase 1 (PGK1) enhances tumorigenesis and sorafenib resistance of kidney renal clear cell carcinoma (KIRC) by regulating glycolysis, so that PGK1 is a promising drug target. Herein we performed structure-based virtual screening and series of anticancer pharmaceutical experiments in vitro and in vivo to identify novel small-molecule PGK1-targeted compounds. As results, the compounds CHR-6494 and Z57346765 were screened and confirmed to specifically bind to PGK1 and significantly reduced the metabolic enzyme activity of PGK1 in glycolysis, which inhibited KIRC cell proliferation in a dose-dependent manner. While CHR-6494 showed greater anti-KIRC efficacy and fewer side effects than Z57346765 on nude mouse xenograft model. Mechanistically, CHR-9464 impeded glycolysis by decreasing the metabolic enzyme activity of PGK1 and suppressed histone H3T3 phosphorylation to inhibit KIRC cell proliferation. Z57346765 induced expression changes of genes related to cell metabolism, DNA replication and cell cycle. Overall, we screened two novel PGK1 inhibitors, CHR-6494 and Z57346765, for the first time and discovered their potent anti-KIRC effects by suppressing PGK1 metabolic enzyme activity in glycolysis.

Metastable States in the Hinge-Bending Landscape of an Enzyme in an Atomistic Cytoplasm Simulation

Many enzymes undergo major conformational changes to function in cells, particularly when they bind to more than one substrate. We quantify the large-amplitude hinge-bending landscape of human phosphoglycerate kinase (PGK) in a human cytoplasm. Approximately 70 μs of all-atom simulations, upon coarse graining, reveal three metastable states of PGK with different hinge angle distributions and additional substates. The "open" state was more populated than the "semi-open" or "closed" states. In addition to free energies and barriers within the landscape, we characterized the average transition state passage time of ≈0.3 μs and reversible substrate and product binding. Human PGK in a dilute solution simulation shows a transition directly from the open to closed states, in agreement with previous SAXS experiments, suggesting that the cell-like model environment promotes stability of the human PGK semi-open state. Yeast PGK also sampled three metastable states within the cytoplasm model, with the closed state favored in our simulation.

Signature of functional enzyme dynamics in quasielastic neutron scattering spectra: The case of phosphoglycerate kinase

We present an analysis of high-resolution quasi-elastic neutron scattering spectra of phosphoglycerate kinase which elucidates the influence of the enzymatic activity on the dynamics of the protein. We show that in the active state the inter-domain motions are amplified and the intra-domain asymptotic power-law relaxation ∝t-α is accelerated, with a reduced coefficient α. Employing an energy landscape picture of protein dynamics, this observation can be translated into a widening of the distribution of energy barriers separating conformational substates of the protein.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Computational engineering of cellulase Cel9A-68 functional motions through mutations in its linker region

Cellulase collective motions design through linker mutations leads to the enhancement of protein flexibility and function.

Molecular cloning, purification and characterization of Brugia malayi phosphoglycerate kinase

Phosphoglycerate kinase (PGK) is a glycolytic enzyme present in many parasites. It has been reported as a candidate molecule for drug and vaccine developments. In the present study, a full-length cDNA encoding the Brugia malayi 3-phosphoglycerate kinase (BmPGK) with an open reading frame of 1.3 kb was isolated and PCR amplified and cloned. The exact size of the BmPGK's ORF is 1377 bps. The BmPGK gene was subcloned into pET-28a (+) expression vector, the expressed enzyme was purified by affinity column and characterized. The SDS-PAGE analysis revealed native molecular weight of recombinant Brugia malayi 3-phosphoglycerate kinase (rBmPGK) to be ∼45 kDa. The enzyme was found sensitive to temperature and pH, it showed maximum activity at 25 °C and pH 8.5. The K_m values for PGA and ATP were 1.77 and 0.967 mM, respectively. The PGK inhibitor, clorsulon and antifilarial drugs albendazole and ivermectin inhibited the enzyme. The specific inhibitor of PGK, clorsulon, competitively inhibited enzyme with K_i value 1.88 μM. Albendazole also inhibited PGK competitively with K_i value 35.39 μM. Further these inhibitory studies were confirmed by docking and molecular simulation of drugs with enzyme. Clorsulon interacted with substrate binding site with glutamine 37 as well as in hinge regions with aspartic acid 385 and valine 387 at ADP binding site. On the other hand albendazole interacted with asparagine 335 residues. These effects were in good association with binding interactions. Thus current study might help in designing and synthesis of effective inhibitors for this novel drug target and understanding their mode of interaction with the potent anthelmintic drugs.

Conformational state distributions and catalytically relevant dynamics of a hinge-bending enzyme studied by single-molecule FRET and a coarse-grained simulation

Over the last few decades, a view has emerged showing that multidomain enzymes are biological machines evolved to harness stochastic kicks of solvent particles into highly directional functional motions. These intrinsic motions are structurally encoded, and Nature makes use of them to catalyze chemical reactions by means of ligand-induced conformational changes and states redistribution. Such mechanisms align reactive groups for efficient chemistry and stabilize conformers most proficient for catalysis. By combining single-molecule Förster resonance energy transfer measurements with normal mode analysis and coarse-grained mesoscopic simulations, we obtained results for a hinge-bending enzyme, namely phosphoglycerate kinase (PGK), which support and extend these ideas. From single-molecule Förster resonance energy transfer, we obtained insight into the distribution of conformational states and the dynamical properties of the domains. The simulations allowed for the characterization of interdomain motions of a compact state of PGK. The data show that PGK is intrinsically a highly dynamic system sampling a wealth of conformations on timescales ranging from nanoseconds to milliseconds and above. Functional motions encoded in the fold are performed by the PGK domains already in its ligand-free form, and substrate binding is not required to enable them. Compared to other multidomain proteins, these motions are rather fast and presumably not rate-limiting in the enzymatic reaction. Ligand binding slightly readjusts the orientation of the domains and feasibly locks the protein motions along a preferential direction. In addition, the functionally relevant compact state is stabilized by the substrates, and acts as a prestate to reach active conformations by means of Brownian motions.

pH-dependent relationship between thermodynamic and kinetic stability in the denaturation of human phosphoglycerate kinase 1

Human phosphoglycerate kinase 1 (hPGK1) is a glycolytic enzyme essential for ATP synthesis, and it is implicated in different pathological conditions such as inherited diseases, oncogenesis and activation of drugs for cancer and viral treatments. Particularly, mutations in hPGK1 cause human PGK1 deficiency, a rate metabolic conformational disease. We have recently found that most of these mutations cause protein kinetic destabilization by significant changes in the structure/energetics of the transition state for irreversible denaturation. In this work, we explore the relationships between protein conformation, thermodynamic and kinetic stability in hPGK1 by performing comprehensive analyses in a wide pH range (2.5-8). hPGK1 remains in a native conformation at pH 5-8, but undergoes a conformational transition to a molten globule-like state at acidic pH. Interestingly, hPGK1 kinetic stability remains essentially constant at pH 6-8, but is significantly reduced when pH is decreased from 6 to 5. We found that this decrease in kinetic stability is caused by significant changes in the energetic/structural balance of the denaturation transition state, which diverge from those found for disease-causing mutations. We also show that protein kinetic destabilization by acidic pH is strongly linked to lower thermodynamic stability, while in disease-causing mutations seems to be linked to lower unfolding cooperativity. These results highlight the plasticity of the hPGK1 denaturation mechanism that responds differently to changes in pH and in disease-causing mutations. New insight is presented into the different factors contributing to hPGK1 thermodynamic and kinetic stability and the role of denaturation mechanisms in hPGK1 deficiency.

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.

Protein Stability, Folding and Misfolding in Human PGK1 Deficiency

Conformational diseases are often caused by mutations, altering protein folding and stability in vivo. We review here our recent work on the effects of mutations on the human phosphoglycerate kinase 1 (hPGK1), with a particular focus on thermodynamics and kinetics of protein folding and misfolding. Expression analyses and in vitro biophysical studies indicate that disease-causing mutations enhance protein aggregation propensity. We found a strong correlation among protein aggregation propensity, thermodynamic stability, cooperativity and dynamics. Comparison of folding and unfolding properties with previous reports in PGKs from other species suggests that hPGK1 is very sensitive to mutations leading to enhance protein aggregation through changes in protein folding cooperativity and the structure of the relevant denaturation transition state for aggregation. Overall, we provide a mechanistic framework for protein misfolding of hPGK1, which is insightful to develop new therapeutic strategies aimed to target native state stability and foldability in hPGK1 deficient patients.

The mechanism of allosteric coupling in choline kinase α1 revealed by the action of a rationally designed inhibitor

Applying a CHOK hold: Combined experimental and computational studies of the binding mode of a rationally designed inhibitor of the dimeric choline kinase α1 (CHOKα1) explain the molecular mechanism of negative cooperativity (see scheme) and how the monomers are connected. The results give insight into how the symmetry of the dimer can be partially conserved despite a lack of conservation in the static crystal structures.

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.

Modeling of solvent flow effects in enzyme catalysis under physiological conditions

A stochastic model for the dynamics of enzymatic catalysis in explicit, effective solvents under physiological conditions is presented. Analytically-computed first passage time densities of a diffusing particle in a spherical shell with absorbing boundaries are combined with densities obtained from explicit simulation to obtain the overall probability density for the total reaction cycle time of the enzymatic system. The method is used to investigate the catalytic transfer of a phosphoryl group in a phosphoglycerate kinase-ADP-bis phosphoglycerate system, one of the steps of glycolysis. The direct simulation of the enzyme-substrate binding and reaction is carried out using an elastic network model for the protein, and the solvent motions are described by multiparticle collision dynamics which incorporates hydrodynamic flow effects. Systems where solvent-enzyme coupling occurs through explicit intermolecular interactions, as well as systems where this coupling is taken into account by including the protein and substrate in the multiparticle collision step, are investigated and compared with simulations where hydrodynamic coupling is absent. It is demonstrated that the flow of solvent particles around the enzyme facilitates the large-scale hinge motion of the enzyme with bound substrates, and has a significant impact on the shape of the probability densities and average time scales of substrate binding for substrates near the enzyme, the closure of the enzyme after binding, and the overall time of completion of the cycle.

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.

Functional domain motions in proteins on the ~1-100 ns timescale: comparison of neutron spin-echo spectroscopy of phosphoglycerate kinase with molecular-dynamics simulation

Protein function often requires large-scale domain motion. An exciting new development in the experimental characterization of domain motions in proteins is the application of neutron spin-echo spectroscopy (NSE). NSE directly probes coherent (i.e., pair correlated) scattering on the ~1-100 ns timescale. Here, we report on all-atom molecular-dynamics (MD) simulation of a protein, phosphoglycerate kinase, from which we calculate small-angle neutron scattering (SANS) and NSE scattering properties. The simulation-derived and experimental-solution SANS results are in excellent agreement. The contributions of translational and rotational whole-molecule diffusion to the simulation-derived NSE and potential problems in their estimation are examined. Principal component analysis identifies types of domain motion that dominate the internal motion's contribution to the NSE signal, with the largest being classic hinge bending. The associated free-energy profiles are quasiharmonic and the frictional properties correspond to highly overdamped motion. The amplitudes of the motions derived by MD are smaller than those derived from the experimental analysis, and possible reasons for this difference are discussed. The MD results confirm that a significant component of the NSE arises from internal dynamics. They also demonstrate that the combination of NSE with MD is potentially useful for determining the forms, potentials of mean force, and time dependence of functional domain motions in proteins.

Molecular insights on pathogenic effects of mutations causing phosphoglycerate kinase deficiency

Phosphoglycerate kinase (PGK) catalyzes an important ATP-generating step in glycolysis. PGK1 deficiency is an uncommon X-linked inherited disorder, generally characterized by various combinations of non-spherocytic hemolytic anemia, neurological dysfunctions, and myopathies. Patients rarely exhibit all three clinical features. To provide a molecular framework to the different pathological manifestations, all known mutations were reviewed and 16 mutant enzymes, obtained as recombinant forms, were functionally and structurally characterized. Most mutations heavily affect thermal stability and to a different extent catalytic efficiency, in line with the remarkably low PGK activity clinically observed in the patients. Mutations grossly impairing protein stability, but moderately affecting kinetic properties (p.I47N, p.L89P, p.C316R, p.S320N, and p.A354P) present the most homogeneous correlation with the clinical phenotype. Patients carrying these mutations display hemolytic anemia and neurological disorders, and,except for p.A354P variant, no myopaty. Variants highly perturbed in both catalytic efficiency (p.G158V, p.D164V, p.K191del, D285V, p.D315N, and p.T378P) and heat stability (all, but p.T378P) result to be mainly associated with myopathy alone. Finally, mutations faintly affecting molecular properties (p.R206P, p.E252A, p.I253T, p.V266M, and p.D268N) correlate with a wide spectrum of clinical symptoms. These are the first studies that correlate the clinical symptoms with the molecular properties of the mutant enzymes. All findings indicate that the different clinical manifestations associated with PGK1 deficiency chiefly depend on the distinctive type of perturbations caused by mutations in the PGK1 gene, highlighting the need for determination of the molecular properties of PGK variants to assist in prognosis and genetic counseling. However, the clinical symptoms can not be understood only on the bases of molecular properties of the mutant enzyme. Different (environmental, metabolic, genetic and/or epigenetic) intervening factors can contribute toward the expression of PGK deficient clinical phenotypes.

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.