Key Points
β-thal erythroblasts display persistent compensatory expression of PKM2.
Mitapivat mitigates Cys-oxidative damage to β-thal erythroblast proteins, with normalization of Prdx2 subcellular compartmentalization.
Visual Abstract
β-thalassemia (β-thal) is a worldwide hereditary red cell disorder characterized by severe chronic anemia. Recently, the pyruvate kinase (PK) activator mitapivat has been shown to improve anemia and ineffective erythropoiesis in a mouse model of β-thal and in patients with non–transfusion-dependent thalassemia. Here, we showed that in vitro CD34+-derived erythroblasts from patients with β-thal (codb039) are characterized by persistent expression of 2 PK isoforms, PKR and PKM2, compared with healthy cells. Activation of PKR and PKM2 via mitapivat promoted significant metabolic reprogramming of β-thal erythroblasts, resulting in higher levels of high-energy phosphate compounds, including adenosine triphosphate (ATP) and triphosphate nucleoside pools. Proteomics analyses revealed an accumulation of PKR, suggesting a possible beneficial effect of mitapivat on the stability of PKs. Increased ATP availability was accompanied by a higher degree of protein phosphorylation, especially in proteins involved in cell cycle regulation at the transcriptional, translational, and posttranslational levels, supporting the effect of mitapivat on erythroid maturation. Upon treatment with mitapivat, β-thal erythroblasts showed decreased markers of oxidation, including cysteine oxidative posttranslational modifications, downregulation of heat shock protein 70 and peroxiredoxin-2 expression, and normalization of the redox-dependent subcellular distribution of the latter enzyme. Collectively, our data support a protective effect of mitapivat in β-thal erythropoiesis, an effect favored by its activation of persistently expressed PKR and PKM2. In addition to the anticipated benefits on energy metabolism, we report that mitapivat treatment mitigated the oxidative damage in β-thal erythropoiesis, ensuring improved β-thal erythroblast maturation and survival.
Introduction
β-thalassemia (β-thal) is a worldwide hereditary red cell disorder, characterized by chronic anemia, related to both reduced red cell life span and ineffective erythropoiesis.1,2 Chronic transfusion and intensive iron chelation are standard treatments for β-thal syndromes,1 but new therapeutic options are being developed, including gene therapy3 and drugs such as luspatercept.4 Recently, the pyruvate kinase (PK) activator mitapivat has been shown to improve anemia and ineffective erythropoiesis in a mouse model of β-thal and in patients with non–transfusion-dependent thalassemia. Here, we showed that in vitro CD34+-derived erythroblasts from patients with β-thal (codb039) are characterized by persistent expression of 2 PK isoforms, PKR and PKM2, compared with healthy cells. Activation of PKR and PKM2 via mitapivat promoted significant metabolic reprogramming of β-thal erythroblasts, resulting in higher levels of high-energy phosphate compounds, including adenosine triphosphate (ATP) and triphosphate nucleoside pools. Proteomics analyses revealed an accumulation of PKR, suggesting a possible beneficial effect of mitapivat on the stability of PKs. Increased ATP availability was accompanied by a higher degree of protein phosphorylation, especially in proteins involved in cell cycle regulation at the transcriptional, translational, and posttranslational levels, supporting the effect of mitapivat on erythroid maturation. Upon treatment with mitapivat, β-thal erythroblasts showed decreased markers of oxidation, including cysteine oxidative posttranslational modifications, downregulation of heat shock protein 70 (HSP70) and peroxiredoxin-2 (Prdx2) expression, and normalization of the redox-dependent subcellular distribution of the latter enzyme. Collectively, our data support a protective effect of mitapivat in β-thal erythropoiesis, an effect favored by its activation of persistently expressed PKR and PKM2. In addition to the anticipated benefits on energy metabolism, we report that mitapivat treatment mitigated oxidative damage in β-thal erythropoiesis, ensuring improved β-thal erythroblast maturation and survival.
Pathophysiologic studies have shown that erythroid oxidation plays a key role in the severity of the β-thal hematologic phenotype.5,6 Indeed, the expression and localization of typical molecular chaperones, such as HSP70, or atypical chaperones, such as Prdx2, play a crucial role in supporting pathologic erythropoiesis and in ensuring erythroid maturation.6-9 In the context of murine β-thal erythropoiesis, we recently demonstrated a cooperation between Prx2 and Nrf2, transcriptional factors involved in the expression of antioxidant and cytoprotective systems, such as heme-oxygenase-1 (HO-1).6-9 In addition, nuclear translocation of Prdx2 has been reported to protect against oxidation-induced DNA damage.10-13
In β-thal erythroid cells, a highly pro-oxidant environment results in ineffective erythropoiesis, a process in part explained by low ATP and altered activity of key enzymes in the glycolytic pathway.14-16 PK is crucial in cell ATP production, being the last enzyme in the payoff steps of the glycolytic pathway; as such, PK is important not only for the survival of mature red cells in circulation but also contributes to meet the energetic needs of a metabolically demanding process such as erythropoiesis.16,17 The PKR isoform is constitutively expressed in erythroid cells, whereas the PKM2 isoform can be induced by oxidative stress or hypoxia.18-21 The importance of PK in erythroid maturation is also supported by signs of ineffective erythropoiesis in patients with PK deficiency and in murine models of PK deficiency.18-22 Recently, the investigational PK activator mitapivat was shown to significantly improve anemia in a β-thal mouse model: specifically, PK activation was beneficial against both ineffective erythropoiesis and chronic hemolysis.17 Mitapivat significantly reduced oxidation during β-thal mouse erythropoiesis, with an improvement of mitochondrial function and progress in erythroid maturation.17 Indeed, PK activation in mitochondria-endowed erythroid precursors favors late glycolytic fluxes, which in turn fuels the generation of late glycolytic end products that can sustain the Krebs cycle. Notably, we found that both PKR and PKM2 were expressed at higher values in murine β-thal red cells fractioned according to red cell age vs wild-type controls, whereas PKR, but not PKM2, was observed in healthy erythrocytes.17 We suggested that this might be compensatory to overcome the severe oxidation of β-thal erythroid cells. As a proof of concept, we tested mitapivat in vitro in CD34+-derived erythroblasts from patients with β-thal (codb039). Mitapivat ameliorated the human erythroid cell maturation index and reduced the amount of Annexin-V–positive β-thal erythroblasts.17 These laboratory studies were also supported by corroborating evidence of improved anemia from a phase 2 study (NCT03692052) in patients with non–transfusion-dependent thalassemia.23,24
Here, we studied CD34+-derived erythroblasts from either healthy controls (HCs) or patients with β-thal (codb039). We showed a protective effect of mitapivat in β-thal erythropoiesis by targeting both PKR and PKM2 through mechanisms involving beneficial effects on both energy and redox metabolism during β-thal erythropoiesis.
Methods
Design of the study
Erythroid cultures obtained from the peripheral blood of 5 different healthy subjects and 5 patients with homozygous β-thal intermedia (β0cod39) were analyzed. The study was approved by the Ethical Committee of the Azienda Ospedaliera Integrata Verona, Italy (FGRF13IT), and written informed consent was obtained from patients with β-thal and HCs. Mitapivat (2 μM) was added to the cell cultures as previously described.17
In vitro human erythropoiesis
Light-density mononuclear cells were obtained by centrifugation on Lymphoprep density gradient (Nycomed Pharma, Oslo, Norway), as previously described.17 The CD34+ cells were positively selected by anti-CD34–tagged magnetic beads (Mini-MACS columns; Miltenyi Biotec, Auburn, CA) according to the manufacturer’s protocol, and erythroid precursors were cultured as previously described.17,25 Details are reported in the supplemental Methods.
Molecular and immunoblot analysis of human erythroid precursors
A total of 1.5 million erythroblasts were collected, solubilized, and analyzed as previously described,17 For quantitative real-time polymerase chain reaction, messenger RNA was isolated and reverse transcribed into high-purity complementary DNA using μMACS (Magnetic-Activated Cell Sorting) One-step cDNA Kit (Miltenyi Biotec) according to the manufacturer’s instructions. Details are reported in the supplemental Methods (supplemental Table 1).
Immunomicroscopy assay
Omics analysis
Biochemical studies on recombinant PKLR and PKM isoforms
Escherichia coli–expressed recombinant PKLR and PKM isoforms were purchased from R&D Systems Bio-Techne. Details of biochemical studies on PKLR and PKM2 are reported in the supplemental Methods.
Statistics, data processing, and visualization
Details are reported in the supplemental Methods.
Results
β-thal erythroblasts display persistent expression of PKM and abnormal metabolomic profile compared with healthy cells
We have recently characterized the maturation profile of healthy and β-thal erythroblasts exposed to mitapivat treatment.29 The ineffective erythropoiesis of β-thal can be recapitulated in vitro, starting from CD34+ cells collected from patients with β-thal.17,25 Compared with healthy cells (Figure 1A; supplemental Figure 1A), β-thal erythroblasts exhibited reduced cell growth and abnormal morphology after 14 days in cell culture (Figure 1A), which corresponds to the late phase of erythropoiesis. We then characterized the expression values of PK isoforms PKLR and PKM, the molecular targets of mitapivat. PKLR gene expression declined during erythroid maturation in both healthy and β-thal erythroblasts (Figure 1B). We found upregulation of PKM gene expression in the late phase of β-thal erythropoiesis (14 days of culture) compared with HCs (Figure 1B). These data were also confirmed at the protein expression level by immunoblot analysis (Figure 1C; supplemental Figure 1B).
In vitro β-thal erythroid precursors show reduced cell proliferation, persistence of PKM expression, and abnormal metabolomic and proteomic profiles compared with healthy cells. (A) May-Grunwald-Giemsa staining for erythroblast morphology (upper) of CD34+-derived erythroid precursors on day 14 of culture from HCs and patients with β-thal (cod β039) showing irregular nuclear shape and chromatin condensation. One representative image from 5 with similar results is shown. Original magnification, 100×. Quantification of abnormal erythroblasts is shown in supplemental Figure 1A. Cell growth (lower) of the erythroid precursors, as in upper panel. Data are mean ± standard error of the mean (SEM; n = 5). ∗P < .05 (compared with HC cells). (B) PKLR and PKM gene expression normalized on HBA1 expression (upper) in erythroid precursors on days 11 and 14 of culture from HCs and patients with β-thal. Data are mean ± SEM (n = 3-5). ∗P < .05 (compared with 11th day); °P < .05 (compared with HC cells). (C) Western blot analysis with specific antibodies against PKLR and PKM2 of erythroid precursors, as in panel A, on day 14 of culture. Catalase was used as loading control. One representative immunoblot of 3 others with similar results. Densitometric analysis is shown in supplemental Figure 1B. (D) Metabolomics profile of erythroid precursors, after 14 days in culture, from HCs and patients with β-thal (cod β039). Heat maps show the top 50 significant metabolites (see also supplemental Figure 8). (E) Heat map focusing on glutathione homeostasis, purine metabolites, glycolysis, Krebs cycle (tricarboxylic acid cycle [TCA]), and PPP. (F) A summary overview of the metabolic pathways in panel E; blue-white-red indicate lower to higher values of a given metabolite in β-thal erythroblasts compared with healthy cells.
In vitro β-thal erythroid precursors show reduced cell proliferation, persistence of PKM expression, and abnormal metabolomic and proteomic profiles compared with healthy cells. (A) May-Grunwald-Giemsa staining for erythroblast morphology (upper) of CD34+-derived erythroid precursors on day 14 of culture from HCs and patients with β-thal (cod β039) showing irregular nuclear shape and chromatin condensation. One representative image from 5 with similar results is shown. Original magnification, 100×. Quantification of abnormal erythroblasts is shown in supplemental Figure 1A. Cell growth (lower) of the erythroid precursors, as in upper panel. Data are mean ± standard error of the mean (SEM; n = 5). ∗P < .05 (compared with HC cells). (B) PKLR and PKM gene expression normalized on HBA1 expression (upper) in erythroid precursors on days 11 and 14 of culture from HCs and patients with β-thal. Data are mean ± SEM (n = 3-5). ∗P < .05 (compared with 11th day); °P < .05 (compared with HC cells). (C) Western blot analysis with specific antibodies against PKLR and PKM2 of erythroid precursors, as in panel A, on day 14 of culture. Catalase was used as loading control. One representative immunoblot of 3 others with similar results. Densitometric analysis is shown in supplemental Figure 1B. (D) Metabolomics profile of erythroid precursors, after 14 days in culture, from HCs and patients with β-thal (cod β039). Heat maps show the top 50 significant metabolites (see also supplemental Figure 8). (E) Heat map focusing on glutathione homeostasis, purine metabolites, glycolysis, Krebs cycle (tricarboxylic acid cycle [TCA]), and PPP. (F) A summary overview of the metabolic pathways in panel E; blue-white-red indicate lower to higher values of a given metabolite in β-thal erythroblasts compared with healthy cells.
Despite the limited blood volume that can be obtained from patients with β-thal, the recent advances in sensitivity of the latest generation Orbitrap and trapped ion mobility mass spectrometers allow for the analysis of samples with low cell counts. By leveraging metabolomics and proteomics approaches (Figures 1D and 2, respectively), we first characterized the baseline differences between β-thal and HC erythroblasts (supplemental Table 2). Significant differences between the 2 groups were noted. β-thal erythroblasts were characterized by significantly lower values of metabolites in the pyrimidine pools (Figure 1D); amino acid values (eg, L-alanine, L-asparagine, glycine, and L-glutamine) and 18-series fatty acids (linoleic and linolenic acid; 18:2; 18:3; Figure 1D). On the contrary, long-chain poly-unsaturated and highly unsaturated fatty acids were enriched in β-thal erythroblasts (especially arachidonate 20:4; eicosa- and docosa- pentaenoic and hexaenoic, 20:5, 20:6, 22:5; 22:6). Similar trends in enrichment of multiple acyl-carnitines (AcCa 18:0, 18:2, and 20:4) and lipid heads (glycerol 3-phosphate and ethanolamine phosphate) were observed. β-thal erythroblasts were also characterized by higher values of pentose phosphate pathway (PPP) metabolites (pentose phosphate isomers, reduced and oxidized nicontinamide adenine dinucleotide phosphate [NADPH], NAD phosphate, and sedoheptulose phosphate); metabolites involved in redox and glutathione homeostasis (ascorbate, glutamyl-alanine, and taurine); and purines and deamination products (adenosine and hypoxanthine; Figure 1D). Figure 1E provides an overview of glycolysis as a heat map, including the PPP, the Krebs cycle (unlike mature red blood cells, erythroblasts are endowed with functional mitochondria), and the adenylate pools. A clear effect on early glycolysis is seen, with accumulation of hexose phosphate compounds accompanied by an apparent rerouting (as inferred by steady-state values) toward the PPP. Elevated fluxes through the PPP have been previously reported in the context of β-thal due to elevated basal oxidant stress.14 Notably, β-thal erythroblasts had higher values of phosphoglycerate (2 or 3 isomers combined), suggestive of a potential bottleneck at the PK value. Consistently, downstream of pyruvate, lactate, citrate, fumarate, and malate (erythroblasts are endowed with mitochondria) were also lower in β-thal erythroblasts. Altered redox homeostasis would also contribute to the competition between reduced nicotinamide adenine dinucleotide (NADH)-dependent methemoglobin reductase and lactate dehydrogenase for NADH recycling. In part, this is justifying further decrease in pyruvate conversion to lactate, consistent with prior reports in other hematologic disorders such as glucose 6-phosphate dehydrogenase deficiency.27,30,31 Results are summarized in the map in Figure 1F.
β-thal erythroblasts are characterized by abnormal proteomic signature and increased overall cysteine oxidation, resulting in upregulation of classic chaperones, such as HSP90 and HSP70, and upregulation of antioxidants, such as peroxiredoxins. Proteomics profile of erythroid precursors from patients with β-thal vs HCs. Specifically, the heat maps show the top 500 (A) and 50 (B) significant proteins by t test. (C) The volcano plot indicates proteins whose values are increased (red) or decreased (blue) in β-erythroblasts compared with healthy counterparts. (D) Heat map showing the relative abundance of all peroxiredoxins (Prdx) in the 2 groups. (E) Western blot analysis with specific antibodies against HSP70 and Prdx-2 of erythroid precursors, as in panel A, on day 14 of culture. Catalase was used as loading control. One representative immunoblot of 3 others with similar results. Densitometric analysis is shown in supplemental Figure 2. (F) Heat map showing cysteine redox modifications (colors were coded as follows: oxidation, orange; dioxidation, red; cysteine to dehydroalanine, magenta) between the 2 groups.
β-thal erythroblasts are characterized by abnormal proteomic signature and increased overall cysteine oxidation, resulting in upregulation of classic chaperones, such as HSP90 and HSP70, and upregulation of antioxidants, such as peroxiredoxins. Proteomics profile of erythroid precursors from patients with β-thal vs HCs. Specifically, the heat maps show the top 500 (A) and 50 (B) significant proteins by t test. (C) The volcano plot indicates proteins whose values are increased (red) or decreased (blue) in β-erythroblasts compared with healthy counterparts. (D) Heat map showing the relative abundance of all peroxiredoxins (Prdx) in the 2 groups. (E) Western blot analysis with specific antibodies against HSP70 and Prdx-2 of erythroid precursors, as in panel A, on day 14 of culture. Catalase was used as loading control. One representative immunoblot of 3 others with similar results. Densitometric analysis is shown in supplemental Figure 2. (F) Heat map showing cysteine redox modifications (colors were coded as follows: oxidation, orange; dioxidation, red; cysteine to dehydroalanine, magenta) between the 2 groups.
Proteomics analyses confirmed a significant remodeling of the erythroblast proteome in β-thal (top 500 proteins by t test in Figure 2A; zoom in for the top 50 differences in Figure 2B). Results are also shown as a volcano plot in Figure 2C and are provided in tabulated form in supplemental Table 2. Of note, results indicated a significant decrease in arginase 2 (ARG2) and elevation in ARG1 in β-thal erythroblasts (Figure 2B-C). The expression values of the enzyme UDP-N-acetylglucosamine-peptide N-acetylglucosaminyl-transferase (OGT) were higher in β-thal erythroblasts, potentially explaining the lower values of its substrate UDP-N-acetylglucosamine in this group (ie, increased consumption as a function of altered glycosylation processes). Pathway analyses of the proteomic results overall indicate an upregulation of proteins involved in chromatin remodeling and cell cycle (CDK1, CEBPZ, and CENPE, etc) and downregulation of protein phosphorylation pathways in β-thal erythroblasts (Figure 2C; supplemental Figures 2 and 3). Higher values of multiple chaperones were observed in the β-thal group (eg, DNAJA1 in Figure 2C). Consistent with a higher value of oxidative stress in the β-thal erythroblast, all peroxiredoxins were found to be elevated in this group (Figure 2D). As an orthogonal validation of the proteomics data, western blot experiments confirmed higher values of HSP70 and Prxd2 in β-thal erythroblasts than in healthy cells (Figure 2E; supplemental Figure 4). Owing to the role of these enzymes in antioxidant defenses, we then focused on redox modifications to protein cysteine residues, as gleaned by proteomics approaches. Results identified a general increase in overall cysteine oxidation in β-thal erythroblasts compared with HCs, especially oxidation, deoxidation, and, to a lesser extent, irreversible β-elimination of the cysteine thiol group to generate dehydroalanine (Figure 2F). Several antioxidant enzymes (eg, catalase [CAT] and glutathione reductases [Gsr]), chaperones (HSP90AA1 and HSP90AAB1), and metabolic enzymes (aldolase [ALDOA], transketolase [TKT], and succinate dehydrogenase B [SDHB]) were found to be more oxidized at functional cysteine residues in the β-thal group (Figure 2F). For example, deoxidation of Cys101 of SDHB contributes to explaining lower fumarate-to-succinate ratios (Figure 1F) and inferred slower fluxes through mitochondrial electron transport chain complex II,32 as gleaned from steady-state data.33 Taken together, our data indicate that in β-thal erythropoiesis, the persistence of PKM2 expression might be compensatory to severe oxidation of β-thal erythroblasts to ensure cell survival and cell growth. However, it has been extensively shown that PKM2 has slower kinetics than the PKM1 or PKLR isoforms.34 Therefore, at the metabolic level, these molecular changes translate into the incapacity to maintain comparable pools of late payoff step values of glycolytic metabolites downstream of PK, suggesting that pharmacological interventions could be used to favor ATP synthesis in metabolically challenged β-thal erythroblasts.
Mitapivat metabolically reprograms human β-thal erythroblasts
Early studies have shown that mitapivat increased the activity of both PKR and PKM2 isoforms, resulting in increased ATP.35,36 Here, we confirmed the ability of mitapivat to activate both recombinant PKR and PKM2 (Figure 3A). The thermal stability profiles showed that PKR is more stable than PKM2. The binding of mitapivat significantly increases protein stability of both isoforms to approximately the same extent (Figure 3B-C). In addition, mitapivat improved resistance to thermal gradient–induced inactivation for both isoforms (Figure 3D-E), with PKM2 being slightly, even if not significant, more resistant to heat inactivation than PKR, consistent with previous reports.35,36 Of note, PKM2 recovered a markedly increased activity with respect to PKR (Figure 3E). This unique functional profile of mitapivat, targeting both PK isoforms, might be extremely useful in a model of pathologic erythropoiesis characterized by the persistence of Pkm2 expression, as in β-thal. Indeed, we have previously shown that mitapivat has limited effects on cell growth, but significantly increases β-thal erythroid maturation and reduces the amount of Annexin-V–positive β-thal cells.17 Here, we confirmed that mitapivat promoted a slight increase in cell growth in treated β-thal cells (Figure 3F). This was associated with an improvement of β-thal erythroblast morphology, as supported by the presence of cells with more regular nuclei and more condensate chromatin (Figure 3F; supplemental Figure 5A). No major changes in PKLR and PKM gene expression were observed in β-thal with or without mitapivat after 14 days in cell culture (supplemental Figure 5B). Proteomics analyses showed an increase in PKR values but not of PKM proteoforms37 upon mitapivat treatment (Figure 4A-B). Despite minor inconsistency with the western blot data due to the different technical approach, as shown in Figure 1A-B (ie, similar trends and different level significance), this result is consistent with the thermal protein stability data (Figure 3E). Thus, the accumulation of PKR in mitapivat-treated cells might represent a regulatory mechanism to reach high activity for both protein species. Mitapivat treatment had widespread significant effects on the erythroblast proteome (Figure 4C), promoting decreases in the values of proteins involved in inflammatory pathways (cytokine-mediated signaling) and JAK/STAT cascades (supplemental Figure 6), with concomitant increase in the values of proteins involved in mitophagy and autophagic responses (supplemental Figure 7). Functionally, mitapivat treatment significantly affected the metabolome of β-thal erythroblasts (Figure 4D), with significant upregulation of multiple high-energy phosphate nucleosides, such as ATP (Figure 4E), CTP, and GTP (Figure 4D). Because we did not observe significant changes in the values of Krebs cycle intermediates, these results are at least in part explained by the decreases in the early and increases in all the late intermediates of glycolysis, suggesting that mitapivat promotes fluxes toward the critical ATP-generating steps of this pathway (Figure 4F). Further studies with stable isotope-labeled substrates (eg, U-13C6-glucose) will be necessary to test this interpretation of the steady-state data. Specifically, mitapivat activation of PK was accompanied by increases in PGLY:DPG ratios, despite observed increases in the values of both metabolites. Significantly lower values of pentose phosphate compounds (isobaric isomers are unresolved by the high-throughput analytical workflow used here) were observed upon mitapivat treatment (Figure 4D), suggesting a decrease in PPP activation, consistent with decreases in oxidant stress in the mitapivat-treated group. Of note, β-thal erythroblasts are characterized by depletion of all pyrimidine pools, especially cytosine- and uridine-containing metabolites that are involved in sialic acid metabolism and membrane-protein glycosylation. This is relevant in light of the well-established differential glycosylation38 and sialic acid profiles of abundant erythrocyte membrane proteins (eg, GYPA) in β-thal.39 Treatment with mitapivat was found to replenish the triphosphate pools for these nucleosides as well.
Mitapivat activates both recombinant PKLR and PKM2, increases protein thermal stability, and ameliorates morphology of β-thal erythroblasts. (A) In vitro activity of recombinant PKLR and PKM2 enzymes incubated with increasing concentrations of mitapivat. Data are mean ± standard deviation (SD; n = 3). (B) Dichroic signal (at 222 nm) of thermal denaturation of PKLR and PKM2 recorded in the absence or presence of saturating concentration of mitapivat (10 μM). Data are mean ± SD (n = 3). (C) Comparison between the calculated melting temperatures and statistical analysis. Data are mean ± SD (n = 3). t test significance, ∗∗∗∗P ≤ 0.0001. (D) Residual activity of PKLR and PKM2 enzymes in the absence or presence of mitapivat (10 μM) after incubation at 53°C. Values are reported as percentage of the initial activity assayed at 25°C. (E) Residual activity of PKLR and PKM2 at 45°C. All error bars are SDs. (F) May-Grunwald-Giemsa staining for erythroblast morphology (upper) of CD34+-derived erythroid precursors, on day 14 of culture, from patients with β-thal (cod β039) with or without mitapivat, showing improvement of irregular nuclear shape and chromatin condensation. One representative image from 5 with similar results is shown. Original magnification, 100×. Quantification of abnormal erythroblasts is shown in supplemental Figure 3A. Cell growth (lower) of the erythroid precursors, as in panel A. Data are mean ± SEM (n= 5). ˆP < .05 (compared with vehicle-treated cells).
Mitapivat activates both recombinant PKLR and PKM2, increases protein thermal stability, and ameliorates morphology of β-thal erythroblasts. (A) In vitro activity of recombinant PKLR and PKM2 enzymes incubated with increasing concentrations of mitapivat. Data are mean ± standard deviation (SD; n = 3). (B) Dichroic signal (at 222 nm) of thermal denaturation of PKLR and PKM2 recorded in the absence or presence of saturating concentration of mitapivat (10 μM). Data are mean ± SD (n = 3). (C) Comparison between the calculated melting temperatures and statistical analysis. Data are mean ± SD (n = 3). t test significance, ∗∗∗∗P ≤ 0.0001. (D) Residual activity of PKLR and PKM2 enzymes in the absence or presence of mitapivat (10 μM) after incubation at 53°C. Values are reported as percentage of the initial activity assayed at 25°C. (E) Residual activity of PKLR and PKM2 at 45°C. All error bars are SDs. (F) May-Grunwald-Giemsa staining for erythroblast morphology (upper) of CD34+-derived erythroid precursors, on day 14 of culture, from patients with β-thal (cod β039) with or without mitapivat, showing improvement of irregular nuclear shape and chromatin condensation. One representative image from 5 with similar results is shown. Original magnification, 100×. Quantification of abnormal erythroblasts is shown in supplemental Figure 3A. Cell growth (lower) of the erythroid precursors, as in panel A. Data are mean ± SEM (n= 5). ˆP < .05 (compared with vehicle-treated cells).
In vitro, mitapivat ameliorates the metabolic and proteomic profiles of human β-thal erythroblasts, associated with decreased overall cysteine oxidation and downregulation of HSP70 expression. PK isoforms PKLR (A) and PKM (all proteoforms) (B) in β-thal erythroblasts as a function of mitapivat treatment. (C) Top 50 proteins (by t test) affected by mitapivat treatment in β-thal erythroblasts (see also supplemental Figure 9). (D) The heat map shows the metabolites that were significantly (t test) affected by mitapivat treatment of β-thal erythroblasts. (E) A detail of ATP, increasing significantly upon mitapivat treatment. (F) Overview of adenylate pools and glycolysis (group averages) in β-thal erythroblasts in the presence or absence of mitapivat (blue vs red indicates low vs high). (G) Elevated ATP values in the mitapivat-treated group corresponded to significant elevation in overall protein phosphorylation. (H) Heat map of the top 50 cysteine redox modifications in β-thal erythroblasts upon mitapivat treatment shows a significant decrease in overall oxidation (orange), deoxidation (red), and cysàdehydroalanine (magenta) in the treated group.
In vitro, mitapivat ameliorates the metabolic and proteomic profiles of human β-thal erythroblasts, associated with decreased overall cysteine oxidation and downregulation of HSP70 expression. PK isoforms PKLR (A) and PKM (all proteoforms) (B) in β-thal erythroblasts as a function of mitapivat treatment. (C) Top 50 proteins (by t test) affected by mitapivat treatment in β-thal erythroblasts (see also supplemental Figure 9). (D) The heat map shows the metabolites that were significantly (t test) affected by mitapivat treatment of β-thal erythroblasts. (E) A detail of ATP, increasing significantly upon mitapivat treatment. (F) Overview of adenylate pools and glycolysis (group averages) in β-thal erythroblasts in the presence or absence of mitapivat (blue vs red indicates low vs high). (G) Elevated ATP values in the mitapivat-treated group corresponded to significant elevation in overall protein phosphorylation. (H) Heat map of the top 50 cysteine redox modifications in β-thal erythroblasts upon mitapivat treatment shows a significant decrease in overall oxidation (orange), deoxidation (red), and cysàdehydroalanine (magenta) in the treated group.
Increased ATP availability can fuel kinase activity, which in turn could contribute to the widespread proteomics changes observed upon treatment with mitapivat. When focusing on phosphopeptide changes, results show a significant increase in total phosphorylated peptides at serine, threonine, and tyrosine residues in the mitapivat-treated group (Figure 4G). Phosphorylation targets include multiple histones (H1.5 at T11, a marker of mitosis),40,41 hematopoietic cell growth–regulating factors (hemogen [HEMGN] at S201 and MKI67 at S2105),42-44 transcriptional/translational regulators (eg, XRN2 and LARP1), and epigenetic regulators (including histone deacetylase 1 [HDAC1] and HDAC2 at functional S421 and S394, respectively, involved in histone complex formation32 and deacetylase function),45-47 heat shock proteins (HSP90AA1 at regulatory S263),48,49 and protein degradation (NEDD8, HUWE1, USP7, and PSMA3) and structural/functional components (glycophorin A [GYPA] at S130; ankyrin [ANK1] at S166; dematin [DMTN] at T27450; but not β spectrin [SPTB] at S2114 and T211051; Figure 4G). We show in supplemental Figures 8 and 9 the top 100 identified metabolites and proteins in all 3 different conditions (healthy cells and β-thal cells treated with either vehicle or mitapivat). Altogether, these results in part contribute to bridging the observed metabolic effects of the mitapivat treatment, consistent with its known molecular target PK, and the overall increase in maturation and longevity observed in treated β-thal erythroblasts.
Mitapivat-treated β-thal erythroblasts showed amelioration of cysteine redox proteome oxidation and downregulation of Prdx2 expression toward a normalization of Prdx2 cell distribution
To further expand on the potential benefits of mitapivat on β-thal erythroblasts with respect to the intracellular redox milieu, we performed redox proteomics analyses of β-thal cells as a function of mitapivat treatment (supplemental Table 2; Figure 4H). Results indicate an overall decrease in total cysteine oxidation, deoxidation, and conversion of cysteine to dehydroalanine (irreversible) in the treated group (heat map in Figure 4H). Specifically, decreases in cysteine oxidation were observed on multiple cysteine residues of several metabolic enzymes (eg, 6-phosphogluconate dehydrogenase [PGD] at C199, in proximity to the active site, and at C422) and structural proteins. The latter group included band 3 (SCL4A1 at C317, functional residue interacting with ankyrin),52 ANK1 (C864), and SPTA1 (at C965), suggestive of a beneficial effect of mitapivat on membrane protein interactions and stability.53 Prdx2, one of the more abundant cytoplasmic proteins of erythroid cells, plays an important role in red cells and during erythropoiesis; specifically, Prdx2 acts as an antioxidant, cytoprotective system, and atypical chaperone.5,6,54,55 As shown in Figure 5A, the improvement in the intracellular milieu of β-thal erythroblasts by mitapivat was also supported by the downregulation of PRDX2 and HMOX1 (heme-oxygenase-1) gene expression after 14 days in cell culture (Figure 5A; supplemental Figure 10A). Both systems were upregulated at both gene and protein levels in response to oxidative stress in vehicle-treated β-thal erythroblasts, in agreement with prior studies by our group on erythroid precursors (Figure 5A-B; supplemental Figure 10A-B).6,17,54 Previous reports in different cell types, such as cancer cell lines or neurons, have shown that differences in the compartmental distribution of Prdx2 between cytoplasm and nucleus affect the role of enzyme as a protective agent against oxidation-induced DNA damage.10-13 When we evaluated Prdx2 localization in β-thal erythroid precursors, we observed a differential localization as a function of mitapivat treatment. As shown in Figure 5C, Prdx2 distribution was essentially cytoplasmic in healthy cells, whereas Prdx2 had cytoplasmic and nuclear compartmentalization in β-thal cells treated with vehicle, compared with healthy cells. Mitapivat treatment reduced Prxdx2 expression (Figure 5B; supplemental Figure 10B) and normalized Prdx2 localization in β-thal cells compared with vehicle-treated β-thal erythroblasts (Figure 5C). In mitapivat-treated β-thal erythroblasts, the expression of HSP70 significantly decreased, reaching values similar to healthy cells, compared with vehicle-treated β-thal erythroid cells (supplemental Figure 10C).
In β-thal erythroblasts, mitapivat normalizes the stress-induced upregulation and nuclear recruitment of Prx-2. (A) PRDX2 gene expression normalized on HBA1 expression in CD34+-derived erythroid precursors, on days 7 (upper) and 14 (lower) of culture, from patients with β-thal (cod β039) in vitro treated with vehicle or mitapivat (2 μM). Data are mean ± SEM (n = 5); ˆP < .05 (compared with vehicle-treated cells). (B) Western blot analysis of Prdx2 in erythroid precursors, as in panel A. Catalase served as loading controls. One representative immunoblot of other 4 with similar results. (C) Prdx2 immunostaining of CD34+-derived erythroid precursors from HC and patients with β-thal (cod β039) in vitro treated with vehicle or mitapivat (2 μM). DAPI (4′,6-diamidino-2-phenylindole) was used to stain nuclei. Prx2 mean fluorescence in the nucleus and cytoplasm was measured using ImageJ. Four different squares and at least 150 cells were analyzed. Data are presented as median and minimum/maximum, with boxes indicating 25th to 75th percentiles; ∗P < .05 (compared with HC); ˆP < .05 (compared with vehicle-treated cells).
In β-thal erythroblasts, mitapivat normalizes the stress-induced upregulation and nuclear recruitment of Prx-2. (A) PRDX2 gene expression normalized on HBA1 expression in CD34+-derived erythroid precursors, on days 7 (upper) and 14 (lower) of culture, from patients with β-thal (cod β039) in vitro treated with vehicle or mitapivat (2 μM). Data are mean ± SEM (n = 5); ˆP < .05 (compared with vehicle-treated cells). (B) Western blot analysis of Prdx2 in erythroid precursors, as in panel A. Catalase served as loading controls. One representative immunoblot of other 4 with similar results. (C) Prdx2 immunostaining of CD34+-derived erythroid precursors from HC and patients with β-thal (cod β039) in vitro treated with vehicle or mitapivat (2 μM). DAPI (4′,6-diamidino-2-phenylindole) was used to stain nuclei. Prx2 mean fluorescence in the nucleus and cytoplasm was measured using ImageJ. Four different squares and at least 150 cells were analyzed. Data are presented as median and minimum/maximum, with boxes indicating 25th to 75th percentiles; ∗P < .05 (compared with HC); ˆP < .05 (compared with vehicle-treated cells).
Discussion
In the context of in vitro human β-thal erythropoiesis, we have previously reported that mitapivat improves the erythroid maturation index and reduces the amount of Annexin-V–positive cells.17 Here, we show that mitapivat, which targets both PKR and PKM2 isoforms, metabolically reprograms β-thal erythroblasts, favoring late glycolysis, which in turn boosts adenylate pools and promotes regulatory protein phosphorylation to enhance cell maturation and growth. Indeed, we found a persistent expression of PKM2 during β-thal erythropoiesis compared with healthy erythroblasts. Papoin et al recently reported a progressive downregulation of PKM2 expression from early to late stages of maturation in healthy human erythroid progenitor cells, suggesting PKR to be the dominant PK isoform in healthy erythroid progenitors entering the late stage of erythropoiesis.56 In human β-thal, we propose that the persistence of PKM2 expression mitigates the metabolic perturbations associated with severe oxidation in β-thal erythroblasts, mainly related to the intracellular pro-oxidant environment. This adaptation holds some relevant metabolic implications, owing to the redox sensitivity of PKM2,57 an important caveat in β-thal erythroblasts, in which basal values of oxidant stress are higher than in HCs. Depletion of glutathione pools and elevation in the expression values and cellular compartmentalization of all peroxiredoxins, especially Prdx2, support the severe and sustain pro-oxidant environment of β-thal erythroblasts. This agrees with our previous reports showing the important cytoprotective and adaptogenic role of Prdx2 against oxidation in β-thal ineffective erythropoiesis.6
Metabolically, mitapivat sustains an upregulation of fluxes of late glycolytic steps with increased ATP content, ensuring cell survival associated with the activation of protein networks involved in chromatin remodeling (supplemental Figure 2) and cell signaling toward cell maturation. Indeed, we observed increased phosphorylation state of multiple targets of cell cycle–regulating kinases, such as phosphorylation of H1.5 at T11, affecting its nuclear localization.40,58 We also observed increases in phosphorylation of HEMG, a regulator of stem cell proliferation via control of STAT1 phosphorylation status,46 and MKI67 (also known as KI-67 proliferation antigen),59 which is consistent with a role of mitapivat in promoting erythroblast maturation. Similarly, casein kinase II targets HDAC1 and HDAC2 are phosphorylated as a function of oxidant stress45 and proliferation state.60,61 Phosphorylation of structural proteins such as GYPA (at S130), ANK1 (at S166), or DMTN (at T274,62 a target of PKA activity) may be one of the metabolic benefits of mitapivat, which leads to improved erythroid maturation compared with vehicle-treated cells. On the contrary, 2 of the 4 common phosphorylation sites for SPTB (S2114 and T2110)51 were found to be less phosphorylated in the β-thal group, suggesting a finely tuned regulation of specific kinase pathways, such as PKA and casein kinase 2, rather than a general upregulation of all phosphorylation events. This is extremely interesting because previous reports have shown that improvements in intracellular cell signaling pathways intersecting EPO-STAT cascade support erythroid maturation.63-65
Another piece of evidence of the beneficial effect of mitapivat on cellular function and signaling is the modulation of autophagy-related proteins, which are involved in one of the key adaptative mechanisms to ensure cell survival. Indeed, previous reports linked impairment of autophagy to ineffective erythropoiesis.66-69 Mitapivat-treated human β-thal cells display activation of mitophagy/autophagy protein networks (supplemental Figure 7). This agrees with our previous report on β-thal mice treated with mitapivat showing improved mitophagy and increased mitochondrial biogenesis in sorted mouse β-thal erythroid precursors.17 The functional amelioration of the mitochondrial compartment contributes to the decreased oxidation observed in mitapivat-treated human β-thal cells compared with vehicle-treated cells. Indeed, mitapivat-treated β-thal cells display (1) lower values of oxidant stress, as gleaned by decreased rerouting of glucose oxidation fluxes to the PPP14; and (2) decreased overall cysteine oxidation. In agreement, we found downregulation of classical chaperone HSP70 and Prdx2 expression compared with vehicle-treated β-thal cells. In addition, mitapivat-treated β-thal cells display Prdx2 cellular distribution similar to that of healthy erythroblasts, characterized more by cytoplasmic than nuclear localization. Recent reports show that Prdx2 might translocate into the nucleus, bind to the promoter regions of ARE genes, and protect erythroblasts from DNA-damaging oxidation related to aging or radiation.5,70,71 The amelioration of intracellular environment of mitapivat-treated β-thal erythroblasts is also supported by the reduction in the expression of HO-1, which is a Nrf2-dependent cytoprotective system involved in heme catabolism.25,55 This is consistent with our previous observation of a functional interplay between Prdx2 and Nrf2, cooperating against oxidation during stress erythropoiesis.6 Of note, Nrf2 might also modulate the expression of HSP70.72 Thus, in β-thal erythroblasts, mitapivat metabolically reprograms β-thal cells and attenuates oxidation, preventing the induction of protective mechanisms such as Prdx2 and HSP70.
Although our data shed new light on the effects of mitapivat on human β-thal erythropoiesis, our study has 2 limitations that could be addressed in future studies. The first one is the cell model. Although it recapitulates in vivo human erythropoiesis and our proteomic/metabolomic profile of normal erythroblasts is similar to that reported in erythroid cells from bone marrow of healthy donors,73-75 future studies might be designed to profile metabolome from bone marrow cells of patients with β-thal. The second limitation concerns the observation on HSP70. Further studies will be conducted to evaluate HSP70 compartmentalization, HSP70 cochaperones, and their link to Nrf2 in mitapivat-treated β-thal erythroblasts.
In conclusion, although caveats are noted owing to the focus on ex vivo models rather than primary in vivo cells, our data show that both PKR and PKM2 isoforms are expressed in the late phase of β-thal erythropoiesis. The ability of mitapivat to stabilize and activate both isoforms might represent an unanticipated benefit of mitapivat treatment, favoring energy metabolism, mitigating oxidant stress in β-thal erythroblasts, and, ultimately, promoting cell maturation and growth.
Acknowledgments
A.D. was supported by funds by the National Heart, Lung, and Blood Institute (R01HL146442, R01HL149714, R01HL148151, and R01HL161004). This research was funded by Ministero dell'Università e della Ricerca PRIN2020-MIUR (2020Z22PM7) (L.D.F. and A.I.) and supported by Agios collaborative research grant (L.D.F.).
The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Authorship
Contribution: L.D.F., M.D., A.A., A.D., and A. Matte. generated data; L.D.F., A.D., and M.B. analyzed the data; L.D.F. and M.D. analyzed the omics; A.D. prepared the figures; C.B., L.D.F., and A.D. wrote the manuscript; A.S. and V.R. performed cell cultures; M.B. and G.B. performed in vitro PKLR/PKM2 experiments; A.I., J.C., A.O., L.R., A. Mattarei, and M.W.-R. critically revised the manuscript; A. Matte. and R.H.P. performed cytofluorimetric, immunomicroscopic, and immunoblot analyses; and R.R. performed molecular analyses, analyzed the data, and wrote the manuscript.
Conflict-of-interest disclosure: A.D. is a founder of Omix Technologies Inc, and is also a scientific advisory board member of Hemanext Inc and Macopharma Inc. L.R. and M.W.-R. are employees and shareholders of Agios Pharmaceuticals. L.D.F. received an Agios research grant. The remaining authors declare no competing financial interests.
Correspondence: Lucia De Franceschi, Department of Engineering for Innovation Medicine, University of Verona & AOUI Verona, P.Le L Scuro, 10, 37134 Verona, Italy; email: lucia.defranceschi@univr.it.
References
Author notes
All data are stored in the Nas-Synology-DS216se hard disk, located at the University of Verona, and will be available on request from the corresponding author, Lucia De Franceschi (lucia.defranceschi@univr.it).
The full-text version of this article contains a data supplement.