Comparison of Proteolytic Susceptibility in Phosphoglycerate Kinases from Yeast and E. coli: Modulation of Conformational Ensembles Without Altering Structure or Stability

Synergistic effect of conformational changes in phosphoglycerate kinase 1 product release

In the glycolysis pathway, phosphoglycerate kinase 1 (PGK1) transfers one phosphoryl-group from 1,3-diphosphoglycerate (1,3BPG) to ADP to product 3-phosphoglycerate (3PG) and ATP. The catalytic process is accompanied with the conversion between the open conformation and the closed conformation of PGK1. However, the dynamic collaboration mechanism between the PGK1 conformation transition and the products releasing process remains poorly understood. Here using molecular dynamics simulations combined with molecular mechanics generalized born surface area (MM/GBSA) analysis, we demonstrated that PGK1 in the closed conformation first releases the product ATP to reach a semi-open conformation, and releases the product 3PG to achieve the full open conformation, which could accept new substrates ADP and 1,3BPG for the next cycle. It is noteworthy that the phosphorylation of PGK1 at T243 causes the loop region (residues L248-E260) flip outside the protein, and the phosphorylation of Y324 leads PGK1 become looser. Both modifications cause the exposure of the ADP/ATP binding site, which was beneficial for the substrates/products binding/releasing of PGK1. In addition, the other post translational modifications (PTMs) were also able to regulate the ligands binding/releasing with different effects. Our results revealed the dynamic cooperative molecular mechanism of PGK1 conformational transition with products releasing, as well as the influence of PTMs, which would contribute to the understanding of PGK1 substrates/products conversion process and the development of small molecule drugs targeting PGK1.Communicated by Ramaswamy H. Sarma.

Investigating the Prevalence of RNA-Binding Metabolic Enzymes in E. coli

An open research field in cellular regulation is the assumed crosstalk between RNAs, metabolic enzymes, and metabolites, also known as the REM hypothesis. High-throughput assays have produced extensive interactome data with metabolic enzymes frequently found as hits, but only a few examples have been biochemically validated, with deficits especially in prokaryotes. Therefore, we rationally selected nineteen Escherichia coli enzymes from such datasets and examined their ability to bind RNAs using two complementary methods, iCLIP and SELEX. Found interactions were validated by EMSA and other methods. For most of the candidates, we observed no RNA binding (12/19) or a rather unspecific binding (5/19). Two of the candidates, namely glutamate-5-kinase (ProB) and quinone oxidoreductase (QorA), displayed specific and previously unknown binding to distinct RNAs. We concentrated on the interaction of QorA to the mRNA of yffO, a grounded prophage gene, which could be validated by EMSA and MST. Because the physiological function of both partners is not known, the biological relevance of this interaction remains elusive. Furthermore, we found novel RNA targets for the MS2 phage coat protein that served us as control. Our results indicate that RNA binding of metabolic enzymes in procaryotes is less frequent than suggested by the results of high-throughput studies, but does occur.

Loss of stability and unfolding cooperativity in hPGK1 upon gradual structural perturbation of its N-terminal domain hydrophobic core

Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype-phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.

Co-Translational Folding of Multi-Domain Proteins

The majority of proteins in nature are composed of multiple domains connected in a single polypeptide. How these long sequences fold into functional structures without forming toxic misfolds or aggregates is poorly understood. Their folding is inextricably linked to protein synthesis and interactions with cellular machinery, making mechanistic studies challenging. Recent progress has revealed critical features of multi-domain protein folding in isolation and in the context of translation by the ribosome. In this review, we discuss challenges and progress in understanding multi-domain protein folding, and highlight how molecular interactions shape folding and misfolding pathways. With the development of new approaches and model systems, the stage is now set for mechanistically exploring the folding of large multi-domain proteins.

Facile tethering of stable and unstable proteins for optical tweezers experiments

Single-molecule force spectroscopy with optical tweezers has emerged as a powerful tool for dissecting protein folding. The requirement to stably attach "molecular handles" to specific points in the protein of interest by preparative biochemical techniques is a limiting factor in applying this methodology, especially for large or unstable proteins that are difficult to produce and isolate. Here, we present a streamlined approach for creating stable and specific attachments using autocatalytic covalent tethering. The high specificity of coupling allowed us to tether ribosome-nascent chain complexes, demonstrating its suitability for investigating complex macromolecular assemblies. We combined this approach with cell-free protein synthesis, providing a facile means of preparing samples for single-molecule force spectroscopy. The workflow eliminates the need for biochemical protein purification during sample preparation for single-molecule measurements, making structurally unstable proteins amenable to investigation by this powerful single-molecule technique. We demonstrate the capabilities of this approach by carrying out pulling experiments with an unstructured domain of elongation factor G that had previously been refractory to analysis. Our approach expands the pool of proteins amenable to folding studies, which should help to reduce existing biases in the currently available set of protein folding models.

Dynamic 3D proteomes reveal protein functional alterations at high resolution in situ

Biological processes are regulated by intermolecular interactions and chemical modifications that do not affect protein levels, thus escaping detection in classical proteomic screens. We demonstrate here that a global protein structural readout based on limited proteolysis-mass spectrometry (LiP-MS) detects many such functional alterations, simultaneously and in situ, in bacteria undergoing nutrient adaptation and in yeast responding to acute stress. The structural readout, visualized as structural barcodes, captured enzyme activity changes, phosphorylation, protein aggregation, and complex formation, with the resolution of individual regulated functional sites such as binding and active sites. Comparison with prior knowledge, including other ‘omics data, showed that LiP-MS detects many known functional alterations within well-studied pathways. It suggested distinct metabolite-protein interactions and enabled identification of a fructose-1,6-bisphosphate-based regulatory mechanism of glucose uptake in E. coli. The structural readout dramatically increases classical proteomics coverage, generates mechanistic hypotheses, and paves the way for in situ structural systems biology.

•

Dynamic structural proteomic screens detect functional changes at high resolution

•

Detect enzyme activity, phosphorylation, and molecular interactions in situ

•

Generate new molecular hypotheses and increase functional proteomics coverage

•

Enabled discovery of a regulatory mechanism of glucose uptake in E. coli

Phosphoglycerate kinase: structural aspects and functions, with special emphasis on the enzyme from Kinetoplastea

Phosphoglycerate kinase (PGK) is a glycolytic enzyme that is well conserved among the three domains of life. PGK is usually a monomeric enzyme of about 45 kDa that catalyses one of the two ATP-producing reactions in the glycolytic pathway, through the conversion of 1,3-bisphosphoglycerate (1,3BPGA) to 3-phosphoglycerate (3PGA). It also participates in gluconeogenesis, catalysing the opposite reaction to produce 1,3BPGA and ADP. Like most other glycolytic enzymes, PGK has also been catalogued as a moonlighting protein, due to its involvement in different functions not associated with energy metabolism, which include pathogenesis, interaction with nucleic acids, tumorigenesis progression, cell death and viral replication. In this review, we have highlighted the overall aspects of this enzyme, such as its structure, reaction kinetics, activity regulation and possible moonlighting functions in different protistan organisms, especially both free-living and parasitic Kinetoplastea. Our analysis of the genomes of different kinetoplastids revealed the presence of open-reading frames (ORFs) for multiple PGK isoforms in several species. Some of these ORFs code for unusually large PGKs. The products appear to contain additional structural domains fused to the PGK domain. A striking aspect is that some of these PGK isoforms are predicted to be catalytically inactive enzymes or ‘dead’ enzymes. The roles of PGKs in kinetoplastid parasites are analysed, and the apparent significance of the PGK gene duplication that gave rise to the different isoforms and their expression in Trypanosoma cruzi is discussed.

Engineering a Trypsin-Resistant Thermophilic α-Galactosidase to Enhance Pepsin Resistance, Acidic Tolerance, Catalytic Performance, and Potential in the Food and Feed Industry

α-Galactosidase has potential applications, and attempts to improve proteolytic resistance of enzymes have important values. We use a novel strategy for genetic manipulation of a pepsin-sensitive region specific for a pepsin-sensitive but trypsin-resistant high-temperature-active Gal27B from Neosartorya fischeri to screen mutants with enhanced pepsin resistance. All enzymes were produced in Pichia pastoris to identify the roles of loop 4 (Gal27B-A23) and its key residue at position 156 (Gly156Arg/Pro/His) in pepsin resistance. Gal27B-A23 and Gly156Arg/Pro/His elevated pepsin resistance, thermostability, stability at low pH, activity toward raffinose (5.3-6.9-fold) and stachyose (about 1.3-fold), and catalytic efficiencies (up to 4.9-fold). Replacing the pepsin cleavage site Glu155 with Gly improved pepsin resistance but had no effect on pepsin resistance when Arg/Pro/His was at position 156. Thus, pepsin resistance could appear to occur through steric hindrance between the residue at the altered site and neighboring pepsin active site. In the presence of pepsin or trypsin, all mutations increased the ability of Gal27B to hydrolyze galactosaccharides in soybean flour (up to 9.6- and 4.3-fold, respectively) and promoted apparent metabolizable energy and nutrient digestibility in soybean meal for broilers (1.3-1.8-fold). The high activity and tolerance to heat, low pH, and protease benefit food and feed industry in a cost-effective way.

Escherichia Coli Increases its ATP Concentration in Weakly Acidic Environments Principally through the Glycolytic Pathway

Acid resistance is an intrinsic characteristic of intestinal bacteria in order to survive passage through the stomach. Adenosine triphosphate (ATP), the ubiquitous chemical used to power metabolic reactions, activate signaling cascades, and form precursors of nucleic acids, was also found to be associated with the survival of Escherichia coli (E. coli) in acidic environments. The metabolic pathway responsible for elevating the level of ATP inside these bacteria during acid adaptation has been unclear. E. coli uses several mechanisms of ATP production, including oxidative phosphorylation, glycolysis and the oxidation of organic compounds. To uncover which is primarily used during adaptation to acidic conditions, we broadly analyzed the levels of gene transcription of multiple E. coli metabolic pathway components. Our findings confirmed that the primary producers of ATP in E. coli undergoing mild acidic stress are the glycolytic enzymes Glk, PykF and Pgk, which are also essential for survival under markedly acidic conditions. By contrast, the transcription of genes related to oxidative phosphorylation was downregulated, despite it being the major producer of ATP in neutral pH environments.

Phosphoglycerate kinase 1 (PGK1) in cancer: A promising target for diagnosis and therapy

Phosphoglycerate kinase 1 (PGK1) is the first critical enzyme to produce ATP in the glycolytic pathway. PGK1 is not only a metabolic enzyme but also a protein kinase, which mediates the tumor growth, migration and invasion through phosphorylation some important substrates. Moreover, PGK1 is associated with poor treatment and prognosis of cancers. This manuscript reviews the structure, functions, post-translational modifications (PTMs) of PGK1 and its relationship with tumors, which demonstrates that PGK1 has indispensable value in the tumor progression. The current review highlights the important role of PGK1 in anticancer treatments.

Unraveling the Mechanical Unfolding Pathways of a Multidomain Protein: Phosphoglycerate Kinase

Phosphoglycerate kinase (PGK) is a highly conserved enzyme that is crucial for glycolysis. PGK is a monomeric protein composed of two similar domains and has been the focus of many studies for investigating interdomain interactions within the native state and during folding. Previous studies used traditional biophysical methods (such as circular dichroism, tryptophan fluorescence, and NMR) to measure signals over a large ensemble of molecules, which made it difficult to observe transient changes in stability or structure during unfolding and refolding of single molecules. Here, we unfold single molecules of PGK using atomic force spectroscopy and steered molecular dynamic computer simulations to examine the conformational dynamics of PGK during its unfolding process. Our results show that after the initial forced separation of its domains, yeast PGK (yPGK) does not follow a single mechanical unfolding pathway; instead, it stochastically follows two distinct pathways: unfolding from the N-terminal domain or unfolding from the C-terminal domain. The truncated yPGK N-terminal domain unfolds via a transient intermediate, whereas the structurally similar isolated C-terminal domain has no detectable intermediates throughout its mechanical unfolding process. The N-terminal domain in the full-length yPGK displays a strong unfolding intermediate 13% of the time, whereas the truncated domain (yPGKNT) transitions through the intermediate 81% of the time. This effect indicates that the mechanical properties of yPGK cannot be simply deduced from the mechanical properties of its constituents. We also find that Escherichia coli PGK is significantly less mechanically stable as compared to yPGK, contrary to bulk unfolding measurements. Our results support the growing body of observations that the folding behavior of multidomain proteins is difficult to predict based solely on the studies of isolated domains.

Engineering of Yersinia Phytases to Improve Pepsin and Trypsin Resistance and Thermostability and Application Potential in the Food and Feed Industry

Susceptibility to proteases usually limits the application of phytase. We sought to improve the pepsin and trypsin resistance of YeAPPA from Yersinia enterocolitica and YkAPPA from Y. kristensenii by optimizing amino acid polarity and charge. The predicted pepsin/trypsin cleavage sites F89/K226 in pepsin/trypsin-sensitive YeAPPA and the corresponding sites (F89/E226) in pepsin-sensitive but trypsin-resistant YkAPPA were substituted with S and H, respectively. Six variants were produced in Pichia pastoris for catalytic and biochemical characterization. F89S, E226H, and F89S/E226H elevated pepsin resistance and thermostability and K226H and F89S/K226H improved pepsin and trypsin resistance and stability at 60 °C and low pH. All the variants increased the ability of the proteins to hydrolyze phytate in corn meal by 2.6-14.9-fold in the presence of pepsin at 37 °C and low pH. This study developed a genetic manipulation strategy specific for pepsin/trypsin-sensitive phytases that can improve enzyme tolerance against proteases and heat and benefit the food and feed industry in a cost-effective way.

The ribosome destabilizes native and non-native structures in a nascent multidomain protein

Correct folding is a prerequisite for the biological activity of most proteins. Folding has largely been studied using in vitro refolding assays with isolated small, robustly folding proteins. A substantial fraction of all cellular proteomes is composed of multidomain proteins that are often not amenable to this approach, and their folding remains poorly understood. These large proteins likely begin to fold during their synthesis by the ribosome, a large molecular machine that translates the genetic code. The ribosome affects how folding proceeds, but the underlying mechanisms remain largely obscure. We have utilized optical tweezers to study the folding of elongation factor G, a multidomain protein composed of five domains. We find that interactions among unfolded domains interfere with productive folding in the full-length protein. The N-terminal G-domain constitutes an independently folding unit that, upon in vitro refolding, adopts two similar states that correspond to the natively folded and a non-native, possibly misfolded structure. The ribosome destabilizes both of these states, suggesting a mechanism by which terminal misfolding into highly stable, non-native structures is avoided. The ribosome may thus directly contribute to efficient folding by modulating the folding of nascent multidomain proteins.

Novel Immunomodulatory Flagellin-Like Protein FlaC in Campylobacter jejuni and Other Campylobacterales

Flagellins not only are important for bacterial motility but are major bacterial proteins that can modulate host responses via Toll-like receptor 5 (TLR5) or other pattern recognition receptors. Campylobacterales colonizing the intestinal tracts of different host species harbor a gene coding for an unusual flagellin, FlaC, that is not involved in motility but is secreted and possesses a chimeric amino acid sequence composed of TLR5-activating and non-TLR5-activating flagellin sequences. Campylobacter jejuni FlaC activates cells to increase in cytokine expression in chicken and human cells, promotes cross-tolerance to TLR4 ligands, and alters chicken cecal microbiota. We propose that FlaC is a secreted effector flagellin that has specifically evolved to modulate the immune response in the intestinal tract in the presence of the resident microbiota and may contribute to bacterial persistence. The results also strengthen the role of the flagellar type III apparatus as a functional secretion system for bacterial effector proteins.

The human diarrheal pathogens Campylobacter jejuni and Campylobacter coli interfere with host innate immune signaling by different means, and their flagellins, FlaA and FlaB, have a low intrinsic property to activate the innate immune receptor Toll-like receptor 5 (TLR5). We have investigated here the hypothesis that the unusual secreted, flagellin-like molecule FlaC present in C. jejuni, C. coli, and other Campylobacterales might activate cells via TLR5 and interact with TLR5. FlaC shows striking sequence identity in its D1 domains to TLR5-activating flagellins of other bacteria, such as Salmonella, but not to nonstimulating Campylobacter flagellins. We overexpressed and purified FlaC and tested its immunostimulatory properties on cells of human and chicken origin. Treatment of cells with highly purified FlaC resulted in p38 activation. FlaC directly interacted with TLR5. Preincubation with FlaC decreased the responsiveness of chicken and human macrophage-like cells toward the bacterial TLR4 agonist lipopolysaccharide (LPS), suggesting that FlaC mediates cross-tolerance. C. jejuni flaC mutants induced an increase of cell responses in comparison to those of the wild type, which was suppressed by genetic complementation. Supplementing excess purified FlaC likewise reduced the cellular response to C. jejuni. In vivo, the administration of ultrapure FlaC led to a decrease in cecal interleukin 1β (IL-1β) expression and a significant change of the cecal microbiota in chickens. We propose that Campylobacter spp. have evolved a novel type of secreted immunostimulatory flagellin-like effector in order to specifically modulate host responses, for example toward other pattern recognition receptor (PRR) ligands, such as LPS.

IMPORTANCE Flagellins not only are important for bacterial motility but are major bacterial proteins that can modulate host responses via Toll-like receptor 5 (TLR5) or other pattern recognition receptors. Campylobacterales colonizing the intestinal tracts of different host species harbor a gene coding for an unusual flagellin, FlaC, that is not involved in motility but is secreted and possesses a chimeric amino acid sequence composed of TLR5-activating and non-TLR5-activating flagellin sequences. Campylobacter jejuni FlaC activates cells to increase in cytokine expression in chicken and human cells, promotes cross-tolerance to TLR4 ligands, and alters chicken cecal microbiota. We propose that FlaC is a secreted effector flagellin that has specifically evolved to modulate the immune response in the intestinal tract in the presence of the resident microbiota and may contribute to bacterial persistence. The results also strengthen the role of the flagellar type III apparatus as a functional secretion system for bacterial effector proteins.

Long-range stabilization of anthrax protective antigen upon binding to CMG2

Protective antigen (PA) mediates entry of edema factor (EF) and lethal factor (LF) into the cytoplasmic space of the cells through the formation of a membrane-spanning pore. To do this, PA must initially bind to a host cellular receptor. Recent mass spectrometry analysis of PA using histidine hydrogen–deuterium exchange (His-HDX) has shown that binding of the von Willebrand factor A (vWA) domain of the receptor capillary morphogenesis protein-2 (CMG2) lowers the exchange rates of the imidazole C₂ hydrogen of several histidines, suggesting that receptor binding decreases the structural flexibility of PA. Here, using His-HDX and fluorescence as a function of denaturant, and protease susceptibility, we show that binding of the vWA domain of CMG2 largely increases the stability of PA and the effect reaches up to 70 Å from the receptor binding interface. We also show that the pK_a values and HDX rates of histidines located in separate domains change upon receptor binding. These results indicate that when one end of the protein is anchored, the structure of PA is tightened, noncovalent interactions are strengthened, and the global stability of the protein increases. These findings suggest that CMG2 may be used to stabilize PA in future anthrax vaccines.

Thioredoxin-dependent redox regulation of chloroplastic phosphoglycerate kinase from Chlamydomonas reinhardtii

In photosynthetic organisms, thioredoxin-dependent redox regulation is a well established mechanism involved in the control of a large number of cellular processes, including the Calvin-Benson cycle. Indeed, 4 of 11 enzymes of this cycle are activated in the light through dithiol/disulfide interchanges controlled by chloroplastic thioredoxin. Recently, several proteomics-based approaches suggested that not only four but all enzymes of the Calvin-Benson cycle may withstand redox regulation. Here, we characterized the redox features of the Calvin-Benson enzyme phosphoglycerate kinase (PGK1) from the eukaryotic green alga Chlamydomonas reinhardtii, and we show that C. reinhardtii PGK1 (CrPGK1) activity is inhibited by the formation of a single regulatory disulfide bond with a low midpoint redox potential (-335 mV at pH 7.9). CrPGK1 oxidation was found to affect the turnover number without altering the affinity for substrates, whereas the enzyme activation appeared to be specifically controlled by f-type thioredoxin. Using a combination of site-directed mutagenesis, thiol titration, mass spectrometry analyses, and three-dimensional modeling, the regulatory disulfide bond was shown to involve the not strictly conserved Cys(227) and Cys(361). Based on molecular mechanics calculation, the formation of the disulfide is proposed to impose structural constraints in the C-terminal domain of the enzyme that may lower its catalytic efficiency. It is therefore concluded that CrPGK1 might constitute an additional light-modulated Calvin-Benson cycle enzyme with a low activity in the dark and a TRX-dependent activation in the light. These results are also discussed from an evolutionary point of view.

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.

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 interplay between protein stability and dynamics in conformational diseases: the case of hPGK1 deficiency

Conformational diseases often show defective protein folding efficiency in vivo upon mutation, affecting protein properties such as thermodynamic stability and folding/unfolding/misfolding kinetics as well as the interactions of the protein with the protein homeostasis network. Human phosphoglycerate kinase 1 (hPGK1) deficiency is a rare inherited disease caused by mutations in hPGK1 that lead to loss-of-function. This disease offers an excellent opportunity to explore the complex relationships between protein stability and dynamics because of the different unfolding mechanisms displayed towards chemical and thermal denaturation. This work explores these relationships using two thermostable mutants (p.E252A and p.T378P) causing hPGK1 deficiency and WT hPGK1 using proteolysis and chemical denaturation. p.T378P is degraded ~30-fold faster at low protease concentrations (here, the proteolysis step is rate-limiting) and ~3-fold faster at high protease concentrations (where unfolding kinetics is rate-limiting) than WT and p.E252A, indicating that p.T378P is thermodynamically and kinetically destabilized. Urea denaturation studies support the decrease in thermodynamic stability and folding cooperativity for p.T378P, as well as changes in folding/unfolding kinetics. The present study reveals changes in the folding landscape of hPGK1 upon mutation that may affect protein folding efficiency and stability in vivo, also suggesting that native state stabilizers and protein homeostasis modulators may help to correct folding defects in hPGK1 deficiency. Moreover, detailed kinetic proteolysis studies are shown to be powerful and simple tools to provide deep insight into mutational effects on protein folding and stability in conformational diseases.

Folding the proteome

Protein folding is an essential prerequisite for protein function and hence cell function. Kinetic and thermodynamic studies of small proteins that refold reversibly were essential for developing our current understanding of the fundamentals of protein folding mechanisms. However, we still lack sufficient understanding to accurately predict protein structures from sequences, or the effects of disease-causing mutations. To date, model proteins selected for folding studies represent only a small fraction of the complexity of the proteome and are unlikely to exhibit the breadth of folding mechanisms used in vivo. We are in urgent need of new methods - both theoretical and experimental - that can quantify the folding behavior of a truly broad set of proteins under in vivo conditions. Such a shift in focus will provide a more comprehensive framework from which to understand the connections between protein folding, the molecular basis of disease, and cell function and evolution.

Structural and energetic basis of protein kinetic destabilization in human phosphoglycerate kinase 1 deficiency

Protein kinetic destabilization is a common feature of many human genetic diseases. Human phosphoglycerate kinase 1 (PGK1) deficiency is a rare genetic disease caused by mutations in the PGK1 protein, which often shows reduced kinetic stability. In this work, we have performed an in-depth characterization of the thermal stability of the wild type and four disease-causing mutants (I47N, L89P, E252A, and T378P) of human PGK1. PGK1 thermal denaturation is a process under kinetic control, and it is described well by a two-state irreversible denaturation model. Kinetic analysis of differential scanning calorimetry profiles shows that the disease-causing mutations decrease PGK1 kinetic stability from ~5-fold (E252A) to ~100000-fold (L89P) compared to that of wild-type PGK1, and in some cases, mutant enzymes are denatured on a time scale of a few minutes at physiological temperature. We show that changes in protein kinetic stability are associated with large differences in enthalpic and entropic contributions to denaturation free energy barriers. It is also shown that the denaturation transition state becomes more nativelike in terms of solvent exposure as the protein is destabilized by mutations (Hammond effect). Unfolding experiments with urea further suggest a scenario in which the thermodynamic stability of PGK1 at least partly determines its kinetic stability. ATP and ADP kinetically stabilize PGK1 enzymes, and kinetic stabilization is nucleotide- and mutant-selective. Overall, our data provide insight into the structural and energetic basis underlying the low kinetic stability displayed by some mutants causing human PGK1 deficiency, which may have important implications for the development of native state kinetic stabilizers for the treatment of this disease.

Expression, purification, crystallization and preliminary crystallographic analysis of the phosphoglycerate kinase from Acinetobacter baumannii

Acinetobacter baumannii is a common multidrug-resistant clinical pathogen that is often found in hospitals. The A. baumannii phosphoglycerate kinase (AbPGK) is involved in the key energy-producing pathway of glycolysis and presents a potential target for antibiotic development. AbPGK has been expressed and purified; it was crystallized using lithium sulfate as the precipitant. The AbPGK crystals belonged to space group P222(1). They diffracted to a resolution of 2.5 Å using synchrotron radiation at the Canadian Light Source.

Conformational properties of nine purified cystathionine β-synthase mutants

Protein misfolding due to missense mutations is a common pathogenic mechanism in cystathionine β-synthase (CBS) deficiency. In our previous studies, we successfully expressed, purified, and characterized nine CBS mutant enzymes containing the following patient mutations: P49L, P78R, A114V, R125Q, E176K, R266K, P422L, I435T, and S466L. These purified mutants exhibited full heme saturation, normal tetrameric assembly, and high catalytic activity. In this work, we used several spectroscopic and proteolytic techniques to provide a more thorough insight into the conformation of these mutant enzymes. Far-UV circular dichroism, fluorescence, and second-derivative UV spectroscopy revealed that the spatial arrangement of these CBS mutants is similar to that of the wild type, although the microenvironment of the chromophores may be slightly altered. Using proteolysis with thermolysin under native conditions, we found that the majority of the studied mutants is more susceptible to cleavage, suggesting their increased local flexibility or propensity for local unfolding. Interestingly, the presence of the CBS allosteric activator, S-adenosylmethionine (AdoMet), increased the rate of cleavage of the wild type and the AdoMet-responsive mutants, while the proteolytic rate of the AdoMet-unresponsive mutants was not significantly changed. Pulse proteolysis analysis suggested that the protein structure of the R125Q and E176K mutants is significantly less stable than that of the wild type and the other mutants. Taken together, the proteolytic data shows that the conformation of the pathogenic mutants is altered despite retained catalytic activity and normal tetrameric assembly. This study demonstrates that the proteolytic techniques are useful tools for the assessment of the biochemical penalty of missense mutations in CBS.

Cystathionine beta-synthase mutants exhibit changes in protein unfolding: conformational analysis of misfolded variants in crude cell extracts

Protein misfolding has been proposed to be a common pathogenic mechanism in many inborn errors of metabolism including cystathionine β-synthase (CBS) deficiency. In this work, we describe the structural properties of nine CBS mutants that represent a common molecular pathology in the CBS gene. Using thermolysin in two proteolytic techniques, we examined conformation of these mutants directly in crude cell extracts after expression in E. coli. Proteolysis with thermolysin under native conditions appeared to be a useful technique even for very unstable mutant proteins, whereas pulse proteolysis in a urea gradient had limited values for the study of the majority of CBS mutants due to their instability. Mutants in the active core had either slightly increased unfolding (p.A114V, p.E302K and p.G307S) or extensive unfolding with decreased stability (p.H65R, p.T191M, p.I278T and p.R369C). The extent of the unfolding inversely correlated with the previously determined degree of tetrameric assembly and with the catalytic activity. In contrast, mutants bearing aminoacid substitutions in the C-terminal regulatory domain (p.R439Q and p.D444N) had increased global stability with decreased flexibility. This study shows that proteolytic techniques can reveal conformational abnormalities even for CBS mutants that have activity and/or a degree of assembly similar to the wild-type enzyme. We present here a methodological strategy that may be used in cell lysates to evaluate properties of proteins that tend to misfold and aggregate and that may be important for conformational studies of disease-causing mutations in the field of inborn errors of metabolism.

Probing protein stability and proteolytic resistance by loop scanning: a comprehensive mutational analysis

Improvement in protein thermostability was often found to be associated with increase in its proteolytic resistance as revealed by comparative studies of homologous proteins from extremophiles or mutational studies. Structural elements of protein responsible for this association are not firmly established although loops are implicated indirectly due to their structural role in protein stability. To get a better insight, a detailed study of protein wide mutants and their influence on stability and proteolytic resistance would be helpful. To generate such a data set, a model protein, Bacillus subtilis lipase was subjected to loop scanning site-saturation mutagenesis on 86 positions spanning all loops including termini. Upon screening of ~16,000 clones, 17 single mutants with improved thermostability were identified with increment in apparent melting temperature (Tm(app) ) by 1-6°C resulting in an increase in free energy of unfolding (ΔG(unf) ) by 0.04-1.16 kcal/mol. Proteolytic resistance of all single mutants upon incubation with nonspecific protease, Subtilisin A, was determined. Upon comparison, post-proteolysis residual activities as well as kinetics of proteolysis of mutants showed excellent correlation with ΔG(unf) , (r > 0.9), suggesting that proteolysis was strongly correlated with the global stability of this protein. This significant correlation in this set, with least possible sequence changes (single aa substitution), while covering >60% of protein surface strongly argues for the covariance of these two variables. Compared to studies from extremophiles, with large sequence heterogeneity, the observed correlation in such a narrow sequence space (ΔΔG(unf) = 1.57 kcal⁻¹) justifies the robustness of this relation.

Unfolding simulations reveal the mechanism of extreme unfolding cooperativity in the kinetically stable alpha-lytic protease

Kinetically stable proteins, those whose stability is derived from their slow unfolding kinetics and not thermodynamics, are examples of evolution's best attempts at suppressing unfolding. Especially in highly proteolytic environments, both partially and fully unfolded proteins face potential inactivation through degradation and/or aggregation, hence, slowing unfolding can greatly extend a protein's functional lifetime. The prokaryotic serine protease α-lytic protease (αLP) has done just that, as its unfolding is both very slow (t_1/2 ∼1 year) and so cooperative that partial unfolding is negligible, providing a functional advantage over its thermodynamically stable homologs, such as trypsin. Previous studies have identified regions of the domain interface as critical to αLP unfolding, though a complete description of the unfolding pathway is missing. In order to identify the αLP unfolding pathway and the mechanism for its extreme cooperativity, we performed high temperature molecular dynamics unfolding simulations of both αLP and trypsin. The simulated αLP unfolding pathway produces a robust transition state ensemble consistent with prior biochemical experiments and clearly shows that unfolding proceeds through a preferential disruption of the domain interface. Through a novel method of calculating unfolding cooperativity, we show that αLP unfolds extremely cooperatively while trypsin unfolds gradually. Finally, by examining the behavior of both domain interfaces, we propose a model for the differential unfolding cooperativity of αLP and trypsin involving three key regions that differ between the kinetically stable and thermodynamically stable classes of serine proteases.

Proteins, synthesized as linear polymers of amino acids, fold up into compact native states, burying their hydrophobic amino acids into their interiors. Protein folding minimizes the non-specific interactions that unfolded protein chains can make, which include aggregation with other proteins and degradation by proteases. Unfortunately, even in the native state, proteins can partially unfold, opening up regions of their structure and making these adverse events possible. Some proteins, particularly those in harsh environments full of proteases, have evolved to virtually eliminate partial unfolding, significantly reducing their rate of degradation. This elimination of partial unfolding is termed “cooperative,” because unfolding is an all-or-none process. One class of proteins has diverged into two families, one bacterial and highly cooperative and the other animal and non-cooperative. We have used detailed simulations of unfolding for members of each family, α-lytic protease (bacterial) and trypsin (animal) to understand the unfolding pathways of each and the mechanism for the differential unfolding cooperativity. Our results explain prior biochemical experiments, reproduce the large difference in unfolding cooperativity between the families, and point to the interface between α-lytic protease's two domains as essential to establishing unfolding cooperativity. As seen in an unrelated protein family, generation of a cooperative domain interface may be a common evolutionary response for ensuring the highest protein stability.

Mapping transient partial unfolding by protein engineering and native-state proteolysis

Transient partial unfolding of proteins under native conditions may have significant consequences in the biochemical and biophysical properties of proteins. Native-state proteolysis offers a facile way to investigate the thermodynamic and kinetic accessibilities of partially unfolded forms (cleavable forms) under native conditions. However, determination of the structure of the cleavable form, which is populated only transiently, remains challenging. Although in some cases partially cleaved products from proteolysis provide information on the structure of this elusive form, proteolysis of many proteins does not accumulate detectable intermediates. Here, we describe a systematic approach to determining structures of cleavable forms by protein engineering and native-state proteolysis. By devising phi(c) analysis, which is analogous to conventional phi analysis, we have determined the structure of the cleavable form of Escherichia coli maltose-binding protein (MBP), which does not accumulate any partially cleaved products. We mutated 10 buried residues in MBP to alanine and determined phi(c) values from the effects of the mutations on global stability and proteolytic susceptibility. The result of this analysis suggests that two C-terminal helices in MBP are unfolded in their cleavable form. The effect of ligand binding on proteolytic susceptibility and C-terminal deletion mutations also confirms the proposed structure. Our approach and methodology are generally applicable not only in elucidating the mechanism of proteolysis but also in investigating other important processes involving partial unfolding under native conditions such as protein misfolding and aggregation.

Between-species variation in the kinetic stability of TIM proteins linked to solvation-barrier free energies

Theoretical, computational, and experimental studies have suggested the existence of solvation barriers in protein unfolding and denaturation processes. These barriers are related to the finite size of water molecules and can be envisioned as arising from the asynchrony between water penetration and breakup of internal interactions. Solvation barriers have been proposed to play roles in protein cooperativity and kinetic stability; therefore, they may be expected to be subject to natural selection. We study the thermal denaturation, in the presence and in the absence of chemical denaturants, of triosephosphate isomerases (TIMs) from three different species: Trypanosoma cruzi, Trypanosoma brucei, and Leishmania mexicana. In all cases, denaturation was irreversible and kinetically controlled. Surprisingly, however, we found large differences between the kinetic denaturation parameters, with T. cruzi TIM showing a much larger activation energy value (and, consequently, much lower room-temperature, extrapolated denaturation rates). This disparity cannot be accounted for by variations in the degree of exposure to solvent in transition states (as measured by kinetic urea m values) and is, therefore, to be attributed mainly to differences in solvation-barrier contributions. This was supported by structure-energetics analyses of the transition states and by application of a novel procedure to estimate from experimental data the solvation-barrier impact at the entropy and free-energy levels. These analyses were actually performed with an extended protein set (including six small proteins plus seven variants of lipase from Thermomyces lanuginosus and spanning a wide range of activation parameters), allowing us to delineate the general trends of the solvation-barrier contributions. Overall, this work supports that proteins sharing the same structure and function but belonging to different organisms may show widely different solvation barriers, possibly as a result of different levels of the selection pressure associated with cooperativity, kinetic stability, and related factors.

Conformational equilibria and rates of localized motion within hepatitis B virus capsids

Functional analysis of hepatitis B virus (HBV) core particles has associated a number of biological roles with the C terminus of the capsid protein. One set of functions require the C terminus to be on the exterior of the capsid, while others place this domain on the interior. According to the crystal structure of the capsid, this segment is strictly internal to the capsid shell and buried at a protein-protein interface. Using kinetic hydrolysis, a form of protease digestion assayed by SDS-PAGE and mass spectrometry, the structurally and biologically important C-terminal region of HBV capsid protein assembly domain (Cp149, residues 1-149) has been shown to be dynamic in both dimer and capsid forms. HBV is an enveloped virus with a T=4 icosahedral core that is composed of 120 copies of a homodimer capsid protein. Free dimer and assembled capsid forms of the protein are readily hydrolyzed by trypsin and thermolysin, around residues 127-128, indicating that this region is dynamic and exposed to the capsid surface. The measured conformational equilibria have an opposite temperature dependence between free dimer and assembled capsid. This work helps to explain the previously described allosteric regulation of assembly and functional properties of a buried domain. These observations make a critical connection between structure, dynamics, and function: made possible by the first quantitative measurements of conformational equilibria and rates of conversion between protein conformers for a megaDalton complex.

MeCP2-chromatin interactions include the formation of chromatosome-like structures and are altered in mutations causing Rett syndrome

hMeCP2 (human methylated DNA-binding protein 2), mutations of which cause most cases of Rett syndrome (RTT), is involved in the transmission of repressive epigenetic signals encoded by DNA methylation. The present work focuses on the modifications of chromatin architecture induced by MeCP2 and the effects of RTT-causing mutants. hMeCP2 binds to nucleosomes close to the linker DNA entry-exit site and protects approximately 11 bp of linker DNA from micrococcal nuclease. MeCP2 mutants differ in this property; the R106W mutant gives very little extra protection beyond the approximately 146-bp nucleosome core, whereas the large C-terminal truncation R294X reveals wild type behavior. Gel mobility assays show that linker DNA is essential for proper MeCP2 binding to nucleosomes, and electron microscopy visualization shows that the protein induces distinct conformational changes in the linker DNA. When bound to nucleosomes, MeCP2 is in close proximity to histone H3, which exits the nucleosome core close to the proposed MeCP2-binding site. These findings firmly establish nucleosomal linker DNA as a crucial binding partner of MeCP2 and show that different RTT-causing mutations of MeCP2 are correspondingly defective in different aspects of the interactions that alter chromatin architecture.