1. Introduction
Sickle cell disease (SCD) is an inherited hematological disorder caused by an A to T mutation in the adult β-globin gene leading to the production of adult hemoglobin S (HbS) which polymerizes under hypoxia conditions [
1,
2]. The polymerization of HbS is the primary event in the pathogenesis of SCD leading to chronic hemolysis, and oxidative stress which contribute to stress erythropoiesis [
3,
4,
5,
6,
7], systemic hypoxia [
8,
9,
10,
11], and abnormal metabolic programming [
12,
13,
14,
15,
16]. Although progress has been made recently, SCD treatment options remain limited. Two critical strategies to treat SCD include γ-globin gene reactivation and reversal of oxidative stress [
17]. For decades, hydroxyurea was the only FDA-approved treatment option for SCD [
18], until the recent approval of Endari (oral L-glutamine) [
19], crizanlizumab [
20] and voxoletor [
21]. Long-term experience with hydroxyurea has shown improvement in the clinical course, primarily through fetal hemoglobin (HbF) induction, but about 30% of SCD patients are not responsive. By contrast, Endari reverses oxidative stress by increasing the NADH and NAD redox potential and decreasing endothelial adhesion of sickle red blood cells [
22]. However, there are multiple sources of pro-oxidant events that contribute to the clinical phenotype of chronic systemic oxidative stress in SCD such as increased generation of free radicals [
23,
24], elevated mitochondrial respiratory chain activity [
25,
26], and blood cell auto-oxidation [
27]. Therefore additional therapeutic agents are under development for treatment of SCD [
28], mainly targeting HbF induction, anti-red blood cell adhesion, pain-alleviation, antioxidative stress, or anti-inflammation. Unfortunately, the complex pathophysiology of SCD makes it unlikely that a single agent will prevent all SCD complications [
29,
30].
The transcription factor Nrf2 (Nuclear factor (erythroid-derived 2)-like 2) is the master regulator of the cellular oxidative stress response, which regulates a variety of antioxidant genes [
31]. Nrf2 expression is predominantly regulated at the protein level through three distinct signaling pathways comprised of KEAP1 (Kelch-like ECH-associated protein 1), Hrd1, and β-TrCP (β-transducin repeat-containing protein), which promote Nrf2 proteasomal degradation by different mechanisms. KEAP1 is the main regulator for Nrf2 degradation through direct protein interactions. Under resting conditions, Nrf2 is sequestered in the cytoplasm by KEAP1 and directed for ubiquitination by the E3 ligase complex. Multiple small chemical compounds such as tert-Butylhydroquinone (tBHQ), diethylmaleate, and 2-cyano-3,12-dioxo-oleana-1,9(11)-dien-28-oic acid (CDDO-Im) modify KEAP1 and disrupt its interaction with Nrf2, leading to Nrf2 nuclear translocation and antioxidant gene activation [
32]. In addition, Hrd1- and β-TrCP mediate Nrf2 degradation, and pharmacologic inhibition of their expression by all-trans retinoic acid [
33], LS-102 [
34], SB216763 [
35] and 2,4-dihydropyrano[2,3-c]pyrazoles [
36] promote Nrf2 protein stability and subsequent activation of the antioxidation response. Of these three pathways, the KEAP1-Nrf2 interaction has been most extensively investigated for its role in erythropoiesis, iron hemostasis and heme metabolism through the expression of AMBP, ABCB6, FECH, HRG1, FTL, FTH1 [
37] and two critical regulators of globin gene expression BCL11A and KLF1 [
37,
38].
Notably, genetic evidence has shown that Nrf2 activation can specifically mitigate the severity of hemolytic anemia, and systemic and local inflammation in transgenic SCD mice via suppressing the proinflammatory response and reduction of reactive oxygen species (ROS) stress [
39,
40]. We demonstrated that Nrf2 ablation exacerbates the symptoms of SCD by decreasing expression of antioxidant factors and the globin genes [
41,
42]. Indeed, Nrf2 was found to mediate γ-globin gene transcription in human normal and sickle erythroid progenitors [
43] through binding the proximal promoter antioxidant response element (ARE) and locus control region hypersensitive sites to facilitate chromatin looping structure during erythropoiesis [
44]. We further investigated the important role of Nrf2 in mediating the ferroptosis stress response in SCD erythroblasts and transgenic animal models. Nrf2 ablation induced iron overload in SCD and metabolic reprogramming to induce histone hypermethylation in erythroblasts and suppressed antioxidant genes to increase lipid oxidation [
42]. These findings indicated that the Nrf2 signaling pathway contributes to the regulation of the globin genes and the antioxidant stress response during erythropoiesis.
To explore HbF induction and enhanced antioxidant capacity for SCD treatment through the Nrf2 signaling pathway, our group and others have conducted investigations to discover agents that mediated Nrf2 activation. However, the clinical trial of the Nrf2 activator Bardoxolone methyl (CDDO-Me) in patients with chronic kidney disease showed the development of adverse cardiovascular events [
45,
46]. Hitherto, the efficacy-toxicity profile generated by different Nrf2 activating drugs has unique pharmacokinetic, pharmacodynamic, toxicokinetic, and toxicodynamic profiles. Therefore, repurposing FDA-approved drugs as novel SCD treatment options has been our focus. Previously dimethyl fumarate (DMF), the methyl ester of fumaric acid that disrupts the interaction between Keap1 and Nrf2 was FDA-approved for the treatment of multiple sclerosis. We later demonstrated the ability of DMF to mediated Nrf2 release and nuclear-translocation to activate γ-globin and anti-oxidative genes [
32].
Other agents such as simvastatin and t-butylhydroquinone can induce γ-globin gene transcription through Nrf2 activation in tissue culture systems [
38,
47]. Moreover, simvastatin is an inhibitor for hydroxymethyl-glutaryl coenzyme A reductase, which activates the PI3K/Akt pathway to suppress GSK-3 activity and Nrf2 stabilization [
48,
49]. In addition to activating the Nrf2 signaling pathway, simvastatin mediates chromatin modifications via inhibition of histone deacetylase and DNA methyltransferase [
50]. Simvastatin was also found to reduce chromatin accessibility at transcriptional enhanced associate domain elements to repress endothelial-to-mesenchymal transition to confer protection against endothelial dysfunction [
51]. Moreover, simvastatin was demonstrated to inhibit HDAC1/2 and induce histone hyperacetylation to promote p21 expression and regulate tumor growth [
52,
53]. Notably, the efficacy of statins in the treatment of type 2 diabetes was attributed partially to DNA hypomethylation at the DHCR24, FAM50B, SC4MOL, AHRR, and ABCG1 gene loci [
54,
55]. The statins can also reverse subtelomeric methylation status in age-related diseases such as diabetes and inhibit Smad 6 and Smad 7 induction in autoimmune diseases [
56,
57]. However, whether simvastatin facilitates a disease modifying effect in SCD and the molecular mechanism involved has not been studied.
In this study, we aim to demonstrate the effect of simvastatin to mitigate HbS sickling and ROS stress levels in SCD erythroblasts in vitro and in SCD transgenic mice. Our mechanistic studies demonstrate the ability of simvastatin to upregulate Nrf2 signaling and silence enhancer of zeste homolog 2 (EZH2)-regulated histone methylation in mediating chromatin structure modifications to activate γ-globin and antioxidant gene transcription. Our findings support the development of the FDA-approved drug, simvastatin as a therapeutical option for treating SCD.
2. Materials and Methods
2.1. Chemicals and Reagents
Hydroxyurea (HU), dimethyl fumarate (DMF), dimethyl sulfoxide (DMSO), and H2O2 were obtained (Sigma-Aldrich, St. Louis, MO). Simvastatin (SIM) was purchased from Thermo Fisher Scientific (Hampton, NH).
2.2. Culture of Human SCD Patient Erythroblasts and H2O2 Exposure
Erythroid progenitors were generated from CD34+ stem cells isolated using microbeads selection of SCD patients peripheral blood mononuclear cells [
44]. The anonymous collection of discarded blood samples from SCD patients undergoing chronic transfusion was classified as exempt by the Institutional Review Board at Augusta University and did not require informed consent. SCD patient CD34+ cells were cultured in a modified Fibach two-phase liquid culture system to generate erythroblasts at different stages of maturation [
58]. Briefly, during phase 1, cells were grown in Iscove’s Modified Dulbecco medium with 15% fetal bovine serum, 15% human AB serum, 10 ng/mL interleukin-3, 50 ng/mL stem cell factor and 2 IU/mL erythropoietin (EPO; Sigma-Aldrich, St. Louis, MO). Phase 2 was initiated on day 7 with a similar medium without stem cell factor and interleukin 3. Chemical treatment of SCD erythroblasts was conducted on day 10 of culture and analyzed on day 12. Treatment with H
2O
2 at 40 µM for 8 h was completed where indicated, along with the different chemical treatments.
2.3. Animal Treatment
The preclinical Townes SCD mouse model (B6;129-
Hbatm1(HBA)TowHbbtm2(HBG1,HBB*)Tow/
Hbbtm3(HBG1,HBB)Tow/J) was characterized previously [
59]. The experiments using the SCD mouse model were performed after approval by the Augusta University Institutional Animal Care and Use Committee. Two to three-month-old SCD mice (both sexes) were treated by intraperitoneal injections of HU (50 mg/kg), DMF (50 mg/kg), SIM (7.5 mg/kg), or water daily, 5 days a week for 4 weeks. Peripheral blood was drawn at week 0, 2 and 4. At the end of the treatment, mice were euthanized, spleen and liver were harvested, snap frozen, and stored for iron, NADPH, NADP+, GSSG, and GSH levels or fixed with 1% formaldehyde for tissue pathology analysis.
2.4. Complete Blood Count and Differential
Peripheral blood was collected in BD Vacutainer EDTA tubes by tail bleeding and automated complete blood count and differentials were completed on the Micros 60 CS/CT machine (HORIBA Medical/ABX Diagnostics, Irvine, CA) according to the manufacturer′s protocol.
2.5. Isolation of Mouse Spleen CD71+ Erythroblasts
Single-cell suspensions of SCD mouse spleen were prepared and CD71+ erythroblasts isolated after labeling with biotinylated CD71 antibody followed by streptavidin MicroBeads separation with MS Columns (Miltenyi Biotec, Auburn, CA) as we previously reported [
42]. The single cells were processed for RNA, protein, and chromatin analysis.
2.6. Cell Proliferation and Viability
SCD erythroblasts at day 10 of culture were seeded at a density of 5 × 105 cells/ml in Phase 2 complete culture medium and treated with indicated concentrations of chemicals. On day 12 of the culture, cells were collected and evaluated for viability by trypan blue staining.
2.7. Flow Cytometry Analysis
Flow cytometry was performed to measure ROS using 2′-7′-dichlorodihydrofluorescein diacetate (DCFH-DA) (Invitrogen, Waltham, MA) staining. Cells were incubated with 10 μM DCFH-DA in complete culture medium for 30 min, washed twice with phosphate buffered saline, and ROS levels measured by flow cytometry analysis.
To measure the percentage of HbF expressing cells (%F-cells), erythroblasts were fixed with 1% formaldehyde, permeabilized with 0.4% Triton X100 and stained with FITC anti-HbF antibody (
Supplemental Table S1); isotype control IgG was used to detect non-specific staining. To measure the percentage of %F-cells in mouse peripheral blood, samples were fixed with 0.1% glutaraldehyde and permeabilized before staining with isotype control or FITC-HbF antibodies similarly.
2.8. RBC Sickling Analysis
In vitro RBC sickling analysis was performed as previously published [
43]. Briefly, Day 12 SCD erythroblasts with different chemical treatments were cultured under hypoxic conditions (1% O
2, 5% CO
2) for 14 h and fixed with 4% formaldehyde before moving to normal air. Using light microscopy, the number of sickle red blood cells was quantified by changes in cell morphology. Bright field images at 20× magnification were obtained on an EVOS Cell Imaging system (Thermo Fisher Scientifics, Hampton, NH). Areas were randomly selected where a single layer of erythrocytes were observed, and 500 cells/per treatment condition were counted in triplicate and the percentage of sickle cells calculated based on total red blood cells counted [
60].
2.9. Quantitative RT-PCR
Total RNA was extracted from cells using Trizol (Invitrogen, Waltham, MA), reverse transcribed, and analyzed by quantitative PCR (qPCR) using gene-specific primers (
Supplemental Table S2) with SYBR Green Supermix (Biorad). Relative levels of gene expression for both target and reference genes were calculated by the 2
-ΔΔCt method based on Ct values. Data were presented as the mean ± standard deviation (SD) of fold-change in expression. All mRNA levels were normalized to endogenous β-actin expression.
2.10. Western Blot Analysis
The whole cell extracts from cultured SCD erythroblasts or isolated SCD mouse spleen CD71+ cells were prepared in RIPA buffer with protease inhibitors (Sigma-Aldrich, St. Louis, MO). For histone extract preparation, cells were homogenized in Triton extraction buffer (PBS containing 0.5% Triton X-100 and 2 mM phenylmethylsulfonyl fluoride) for nuclei isolation, and nuclear proteins were prepared by 0.2 N HCl acid extraction followed by neutralization with 1/10 volume of 2 M NaOH. Protein concentrations of the whole-cell and histone extracts were determined using the Bradford assay.
Whole-cell protein extract (30 μg) or histone extract (1 μg) were denatured, subjected to SDS-PAGE, and transferred to nitrocellulose membranes. After blocking with 5% nonfat dry milk in Tris-buffered saline and 0.1% Tween (TBS-T), the membranes were incubated with the indicated antibodies overnight at 4oC. β-actin and histone H3 were used as loading controls. After 1 h incubation with horseradish peroxidase-conjugated anti-Mouse or anti-Rabbit secondary antibodies (Santa Cruz Biotechnology, Dallas, TX), the membrane was washed and visualized by chemiluminescence in a FUJIFILM LAS-3000 system or GE Amersham ImageQuant 800 and quantified with ImageJ software.
2.11. GSH and GSSG Assay
The levels of reduced (GSH) and oxidized (GSSG) glutathione were measured using a GSH/GSSG-Glo™ Assay Kit (Promega, Madison, WI). Briefly, cultured SCD erythroblasts were washed with PBS once and three times of freeze-thaw out in ice-cold PBS containing 2 mM EDTA. The extracts were centrifuged, and the supernatant was collected and used immediately for the assay according to the manufacturer’s instructions. Luminescence was read by a Victor3 Multilabel Counter (Perkin Elmer, Waltham, MA) and GSH and GSSG levels normalized to cell numbers.
2.12. NQO1 Activity Assay
The cultured SCD erythroblasts were washed twice with PBS and assessed for NQO1 activity in whole cell extracts following reduction of menadione with cofactor NADH and WST1. The samples were assayed at 440 nm absorbance using a NQO1 Activity Assay Kit (Abcam, Waltham, MA), according to the manufacturer’s instructions. The NQO1 activity was normalized to protein content.
2.13. Measurement of Intracellular Cysteine Levels
Day 12 cultured erythroblasts were harvested, washed with PBS buffer once, and resuspended in distilled water. The intracellular amino acids were extracted by boiling for 10 min. The supernatant was collected after centrifugation and protein content was determined. The levels of cysteine in the supernatant were measured by a Cysteine Assay Kit (Sigma-Aldrich, St. Louis, MO) along with the cysteine standards. Cysteine levels after different treatments were compared after normalization to cell numbers.
2.14. Cystine Uptake Assay
The cystine uptake assay was conducted with BioTracker Cystine-FITC Live Cell Dye (Sigma, St. Louis, MO). Briefly, SCD Day 10 erythroblasts were cultured in a complete IMDM medium with cytokines and treated with DMSO or simvastatin for 48 h. Cystine-FITC was added to cultured SCD erythroblasts at 200 µM for 45 min allowing uptake then the cells were collected and washed followed by quantification of cellular levels of Cystine-FITC by flow cytometry analysis.
2.15. EZH2 Gene Knockdown Analysis
On day 2 of culture, SCD progenitors were transduced with lentiviral EZH2 shRNA particles and selected with puromycin until analysis. Two different shRNAs targeting sequences including ugacuucugugagcucauugc and guuuguuggcggaagcgugua (
Supplemental Table S2) were constructed in pLKO.1 (Addgene, Watertown, MA, USA). A standard scramble sequence (shControl) targeting a nonspecific sequence ccuaagguuaagucgcccucg was used as a control. The sh
EZH2 knockdown efficacy was determined by Western blot analyses.
2.16. Chromatin Immunoprecipitation (ChIP)
ChIP assays were performed with indicated antibodies (
Supplemental Table S1) using gene-specific primers (
Supplemental Table S2) as previously reported with minor modifications [
44,
61]. Samples were immunoprecipitated with anti-NRF2, anti-TATA binding protein (TBP), anti-RNA polymerase II (Pol II), and H3K27Me3 antibodies (
Supplemental Table S1). ChIP DNA was quantified by qPCR with gene locus specific primers to determine the pull-down signals. Rabbit normal IgG was used as a non-specific antibody control.
2.17. NADP+/NADPH Levels
Mouse spleen CD71+ cells (2x10^6) were washed with ice-cold PBS and homogenized in NADP or NADPH extraction buffer before being assayed with EnzyChrom™ NADP/NADPH Assay Kit (BioAssay Systems, Hayward, CA) according to the manufacturer’s instructions. The extracts were heated at 60oC for 5 min and then 20 μl of assay buffer and 100 μl of opposite extraction buffer was added to neutralize the extracts. The levels of NADP+/NADPH in the supernatant were determined by adding enzymes provided in the kit. Absorbance at 570 nm was measured according to the manufacturer’s instructions using a Synergy H1 hybrid plate reader (Biotek, Winooski, VT) to calculate NADP+/NADPH concentrations.
2.18. Histology, Immunohistochemistry, and Image Analysis
Mouse spleen and liver tissues were fixed in 10% neutral buffered formalin, dehydrated, and embedded in paraffin. Tissue sections (5 μm thickness) were subjected to hematoxylin and eosin (H&E) staining and subsequent Prussian blue staining (IHCworld, Ellicott City, MD) to visualize iron deposition.
2.19. Iron Content
Approximate 20 mg of mouse spleen and liver tissues were homogenized in 1 ml of distilled water and supernatants were collected by centrifugation at 4OC. The iron content in the supernatant was measured by an Iron assay kit (Bioassay Systems, Hayward, CA) and iron levels were normalized to protein content.
2.20. Statistical Analysis
Data from at least three independent biological replicates were reported as the mean ± SD. Statistical differences were determined by an unpaired Student’s t-test or by two-way ANOVA with the corresponding two-tailed significance (p value) determined. Statistical analysis was performed using GraphPad Prism 9 software (GraphPad Software Inc., San Diego, CA) and differences were considered significant at p < 0.05.
4. Discussion
Sickle cell disease is produced by a single nucleotide mutation in the adult β-globin gene leading to HbS synthesis and polymerization under hypoxia conditions, resulting in intravascular hemolysis, cell adhesion, vascular occlusion, and ischemia–reperfusion injury from elevated oxidative stress. In addition to multiple approaches for γ-globin activation [
17,
28], other strategies to treat SCD include targeting the downstream long-term effects by reducing oxidative stress, cellular adhesion mediated by P-selectin interactions and increasing hemoglobin oxygen affinity to directly block HbS polymerization. These four strategies have the potential to inform drug combinations to best treat the clinical complications and severity of individual with SCD. In this study, we have chosen to explore the possibility of repurposing the novel FDA-approved agent SIM, which stabilized Nrf2 protein levels and γ-globin activation by creating open chromatin domains, and antioxidant and anti-inflammatory effects.
In the current study, we performed mechanistic studies which showed that SIM decreased global and gene loci specific H3K27me3 levels to affect heterochromatin for gene regulation in SCD erythroblasts (grown in culture) and SCD transgenic mice treatment in vivo. This effect of SIM facilitated an activation of the γ-globin and antioxidant genes at the mRNA and protein levels to reduce red blood cell sickling and ROS stress, demonstrating its dual effect in modifying phenotypic severity of SCD. Depending on the cellular environment, SIM can produce antioxidant effects by increasing the levels of reduced GSH such as in liver [
66], generate an pro-oxidative stress in malignant conditions [
67], or show no effect on oxidative stress regulation in healthy subjects [
68]. Here we demonstrated that SIM promotes cellular uptake of cystine, the precursor for GSH biosynthesis, through Slc7a11 upregulation. In addition, both the catalytic (GCLC) and mediator (GCLM) subunits of the glutamate–cysteine ligase were upregulated by SIM to accelerate GSH biosynthesis and scavenging of ROS stress. Together, our findings and previous observations suggest a cellular environment-dependent effect of SIM in regulating oxidative stress.
Previously, SIM was shown to inhibit histone deacetylase activity specifically targeting HDAC1/2 at a comparable capability to well-studied HDAC inhibitors such as Trichostatin A and suberoylanilide hydroxamic acid [
52,
69,
70]. Simvastatin also mediates induction of histone acetylation by activating the AMPK signaling pathway [
63,
71]. These HDAC inhibit functions of SIM affects the expression of genes whose products participate in cell proliferation. In addition, SIM modifies the expression of key epigenetic proteins such as DNMT1 and influences DNA methylation to regulate gene transcription [
72,
73]. Our results in SCD erythroblasts showed that SIM did not change histone acetylation levels, however, SIM reduced the suppressive histone code H3K27Me3, through EZH2 activation, but not through EZH1 or the histone demethylases KDM6a/b. Although we did not explore the mechanism of SIM-mediated reduction of EZH2, a previous study supports the ability of SIM to inhibit the function of METTL3 in downregulating EZH2 signals [
74]. Those previous reports suggest cellular environment-specific mechanisms of SIM regulating histone modifications [
75]. Of note, we previously showed that Nrf2 mediates metabolic reprogramming in erythroblasts to modify levels of the metabolite L-2-hydroxyglutarate, which is a competitor of α-ketoglutarate-dependent histone demethylases and thus regulating histone methylation [
42]. Since we demonstrated that SIM enhances Nrf2 levels, our results raise the question whether SIM could affect metabolic profiling, histone methylation and gene expression through the Nrf2 signaling pathway in SCD erythroblasts. Future work will focus on answering this question.
From ChIP assay analysis using primers spanning the gene-specific ARE motifs, our results support Nrf2 activation and binding to target ARE motifs to activate downstream target γ-globin and antioxidant genes through EZH2 suppression by SIM to produce an open chromatin structure and increase binding of the general transcription factors, TBP and RNA Pol II. Previously, Nrf2 was shown to produce ARE-independent effects in regulating the expression of inflammatory factors, through binding to the proximity of the interleukin-6 and interleukin-1β genes to inhibit RNA Pol II recruitment and gene transcription in macrophages [
76]. Though we are not able to exclude Nrf2 binding to regions devoid of ARE motifs to mediate γ-globin and antioxidant genes transcription, SIM reduced global and loci-specific H3K27Me3 levels, which led to chromatin structure changes at ARE motifs and other regions as supported by increased TBP and Pol II binding to promote gene transcription.
Related to the molecular function of EZH2, previous studies have shown the ability of EZH2 to act as a corepressor in cooperation with BCL11A, to achieve γ-globin gene silencing [
77]. Indeed, EZH2 is broadly distribute at high levels across the
HBB locus in non-erythroid cells compared to erythroid cells, suggesting its general suppressive effect on globin gene transcription [
78]. Moreover, the EZH2 inhibitor Tazemetostat can significantly increase fetal hemoglobin expression in healthy and β-thalassemic human primary erythroid cells [
79]. Therefore, those observations and our finding of SIM-mediated EZH2/H3K27Me3 suppression to activate γ-globin transcription could expand treatment options for SCD by inhibiting EZH2 expression.
In addition to activating HbF and antioxidant proteins, SIM suppressed pro-inflammatory signals in the SCD mouse model, which is in agreement with previous findings from human SCD clinical trials, where SIM administration reduced markers of inflammation and frequency of vaso-occlusive pain [
80]. In a single center pilot study of adolescents and adults with SCD, SIM treatment for up to 3 months led to a dramatic reduction in the rate of SCD related pain and oral analgesic use; there was also improvement in soluble biomarkers of inflammation and an acceptable safety profile [
81]. Inflammation and oxidative stress occur in chronic diseases and are closely correlated to each other, and mutually influenced. Therefore, the failure of antioxidant drug trials might result from failure to select appropriate agents that will target inflammation, and vice versa [
82], whereas targeting both inflammation and oxidative stress simultaneously might generate synergistic effects. Importantly, oxidative and inflammatory stresses result from HbS polymerization and chronic hemolysis in SCD, therefore SIM has the potential to produce a three-fold effect including HbF induction, and inflammation and oxidative stress reduction to ameliorate the clinical severity of SCD.
A recent meta-analysis showed that statins can significantly decrease serum ferritin levels, which might be beneficial for the prevention and progression of atherosclerotic cardiovascular disease [
83]. Indeed, iron overload causes inflammation and iron-dependent non-apoptotic cell death, i.e. ferroptosis. Individuals with SCD often suffer iron overload due to frequent blood transfusions, eventually requiring treatment with chelating agents to prevent liver and heart toxicity [
84]. Our in vivo studies in SCD mice showed that SIM reduced the iron content in spleen and liver tissues. However, a link between chronic SIM use and iron deficiency was detected in humans [
85]. Iron reduction by SIM is through the regulation of iron metabolism by HMOX1 expression [
86]. We also observed an increase in HMOX1 expression in SCD mice treated chronically with SIM. Although at the SIM dose of 7.5 mg/kg/day, we did not observe adverse effects in SCD mice treated up to 4 weeks, a lower dose of SIM might be more beneficial and should be evaluated.
Author Contributions
Conceptualization, C.X., B.S.P. and X.Z.; Data curation, C.X., C.P. and M.T.; Formal analysis, C.X., C.P.; H.S., A.H. and B.S.P.; Funding acquisition, B.S.P. and X.Z.; Investigation, C.X., C.P. and M.T.; Methodology, C.X., C.P. and M.T.; Project administration, X.Z.; Resources, A.H.; Supervision, X.Z.; Writing – original draft, C.X. and C.P.; Writing – review & editing, B.S.P. and X.Z. All authors have read and agreed to the published version of the manuscript.
Figure 1.
Simvastatin activates γ-globin expression to suppress sickling of SCD erythroblasts. (A) Two-phase culture of SCD patients’ peripheral blood CD34+ stem cells. (B) SCD erythroblasts at Day 10 of culture were treated with 0.1% DMSO, 50 µM hydroxyurea (HU), 100 µM dimethyl fumarate (DMF), 2.5-20 µM simvastatin (SIM), or untreated control (UT) for 48 h and protein expression of fetal hemoglobin (HbF) and adult sickle hemoglobin (HbS) was quantified. (C-D) SCD erythroblasts from the same culture were determined the percentage of HbF positive cells (F-cells %) by flow cytometry (C) and quantified (D). (E) Cell morphology was analyzed to measure the number of sickled erythroblasts under hypoxia for the different treatment conditions and the percentage of sickled cells quantified (F). Data represents the mean ± SD of three independent biological replicates. *, p < 0.05.
Figure 1.
Simvastatin activates γ-globin expression to suppress sickling of SCD erythroblasts. (A) Two-phase culture of SCD patients’ peripheral blood CD34+ stem cells. (B) SCD erythroblasts at Day 10 of culture were treated with 0.1% DMSO, 50 µM hydroxyurea (HU), 100 µM dimethyl fumarate (DMF), 2.5-20 µM simvastatin (SIM), or untreated control (UT) for 48 h and protein expression of fetal hemoglobin (HbF) and adult sickle hemoglobin (HbS) was quantified. (C-D) SCD erythroblasts from the same culture were determined the percentage of HbF positive cells (F-cells %) by flow cytometry (C) and quantified (D). (E) Cell morphology was analyzed to measure the number of sickled erythroblasts under hypoxia for the different treatment conditions and the percentage of sickled cells quantified (F). Data represents the mean ± SD of three independent biological replicates. *, p < 0.05.
Figure 2.
Simvastatin increases the expression of NRF2 and antioxidant proteins. Day 12 SCD erythroblasts after 48 h treatment with 0.1% DMSO, 50 µM HU, 100 µM DMF, 2.5-20 µM simvastatin (SIM), or untreated control (UT) were detected for the expression of antioxidant factors at the protein (A) and mRNA levels (B). Data represents the mean ± SD of three biological replicates. *, p < 0.05.
Figure 2.
Simvastatin increases the expression of NRF2 and antioxidant proteins. Day 12 SCD erythroblasts after 48 h treatment with 0.1% DMSO, 50 µM HU, 100 µM DMF, 2.5-20 µM simvastatin (SIM), or untreated control (UT) were detected for the expression of antioxidant factors at the protein (A) and mRNA levels (B). Data represents the mean ± SD of three biological replicates. *, p < 0.05.
Figure 3.
Simvastatin increased the antioxidant capacity of SCD erythroblasts. Incubation with H2O2 (40µM) for 8 h was completed with Day 12 SCD erythroblasts after 48 h treatment of 2.5-20 µM simvastatin (SIM), then cellular viability (A) and the cellular content of reduced (GSH), oxidative (GSSG) glutathione and their ratio (B-D) were determined. The cellular NQO1 activity (E) and cystine contents (F) in the same cells cultured as described in panel A were measured with assay kits. (G) The cystine uptake efficacy of SCD erythroblasts was determined by a BioTracker Cystine-FITC Live Cell Dye followed by flow cytometry analysis. Data represents the mean ± SD of three biological replicates. *, p < 0.05.
Figure 3.
Simvastatin increased the antioxidant capacity of SCD erythroblasts. Incubation with H2O2 (40µM) for 8 h was completed with Day 12 SCD erythroblasts after 48 h treatment of 2.5-20 µM simvastatin (SIM), then cellular viability (A) and the cellular content of reduced (GSH), oxidative (GSSG) glutathione and their ratio (B-D) were determined. The cellular NQO1 activity (E) and cystine contents (F) in the same cells cultured as described in panel A were measured with assay kits. (G) The cystine uptake efficacy of SCD erythroblasts was determined by a BioTracker Cystine-FITC Live Cell Dye followed by flow cytometry analysis. Data represents the mean ± SD of three biological replicates. *, p < 0.05.
Figure 4.
Simvastatin attenuates EZH2 expression and histone H3K27Me3 levels to modify chromatin structure and gene expression. (A) Shown are the effects of SIM on histone H3 methylation (H3K27Me3), acetylation (AcH3), and EZH2 protein expression. Total histone H3 and β-ACTIN were loading controls. (B) Quantitative RT-PCR analysis of the relative mRNA levels of H3K27Me3 modifiers EZH1/2 and KDM6a/b in SCD erythroblasts treated with 2.5-20 µM simvastatin (SIM) or untreated control (UT). Data represents mean ± SD of three biological replicates. *, p < 0.05.
Figure 4.
Simvastatin attenuates EZH2 expression and histone H3K27Me3 levels to modify chromatin structure and gene expression. (A) Shown are the effects of SIM on histone H3 methylation (H3K27Me3), acetylation (AcH3), and EZH2 protein expression. Total histone H3 and β-ACTIN were loading controls. (B) Quantitative RT-PCR analysis of the relative mRNA levels of H3K27Me3 modifiers EZH1/2 and KDM6a/b in SCD erythroblasts treated with 2.5-20 µM simvastatin (SIM) or untreated control (UT). Data represents mean ± SD of three biological replicates. *, p < 0.05.
Figure 5.
EZH2 regulates NRF2 expression to alter ARE motif chromatin structure in target genes. (A) The effect of shEZH2 mediated gene silencing (two different constructs) on protein levels of NRF2, HbF, antioxidant factors, and H3K27Me3 were investigated. Histone H3 and β-ACTIN were loading controls. (B-E) ChIP assay determined the association of Nrf2, H3K27Me3, TATA box binding protein (TBP) and RNA polymerase II (Pol II) to the HBG1 (B), HMOX1 (C), NQO1 (D) and SLC7A11 (E) gene loci for simvastatin (SIM, 5µM) treated or untreated (UT) SCD erythroblasts. Primers spanning the ARE motifs of each gene were used. Data represents the mean ± SD of three biological replicates. *, p < 0.05.
Figure 5.
EZH2 regulates NRF2 expression to alter ARE motif chromatin structure in target genes. (A) The effect of shEZH2 mediated gene silencing (two different constructs) on protein levels of NRF2, HbF, antioxidant factors, and H3K27Me3 were investigated. Histone H3 and β-ACTIN were loading controls. (B-E) ChIP assay determined the association of Nrf2, H3K27Me3, TATA box binding protein (TBP) and RNA polymerase II (Pol II) to the HBG1 (B), HMOX1 (C), NQO1 (D) and SLC7A11 (E) gene loci for simvastatin (SIM, 5µM) treated or untreated (UT) SCD erythroblasts. Primers spanning the ARE motifs of each gene were used. Data represents the mean ± SD of three biological replicates. *, p < 0.05.
Figure 6.
Simvastatin suppresses H3K27Me3 modifications as a mechanism of fetal γ-globin activation in preclinical SCD mice. (A) Peripheral blood samples of SCD mice after 4 weeks of simvastatin treatment were determined the HbF expressing cells (F-cells) by flow cytometry. The F-cell % and mean fluorescence intensity (MFI) were quantified. (B-C) The spleen CD71+ cells from SIM-treated SCD mice were determined the expression of H3K27Me3, HbF, and HbS at the protein (B) and mRNA levels (C). (D) ChIP assay determined the association of H3K27Me3, TBP, and Pol II to the HBG1 and HBA gene loci in spleen CD71+ cells from SIM treated SCD mice. One-way ANOVA with Bonferroni’s multiple comparison tests was used for statistical analysis (n=5). *, p < 0.05.
Figure 6.
Simvastatin suppresses H3K27Me3 modifications as a mechanism of fetal γ-globin activation in preclinical SCD mice. (A) Peripheral blood samples of SCD mice after 4 weeks of simvastatin treatment were determined the HbF expressing cells (F-cells) by flow cytometry. The F-cell % and mean fluorescence intensity (MFI) were quantified. (B-C) The spleen CD71+ cells from SIM-treated SCD mice were determined the expression of H3K27Me3, HbF, and HbS at the protein (B) and mRNA levels (C). (D) ChIP assay determined the association of H3K27Me3, TBP, and Pol II to the HBG1 and HBA gene loci in spleen CD71+ cells from SIM treated SCD mice. One-way ANOVA with Bonferroni’s multiple comparison tests was used for statistical analysis (n=5). *, p < 0.05.
Figure 7.
Simvastatin reduces organ damage and iron content in preclinical SCD mice. (A-B) Representative images of the spleen (A) and liver (B) of SCD mice treated with or without simvastatin (SIM, 7.5 mg/kg daily) for 4 weeks. Quantification of spleen and liver weights are shown in the graph. (C) Representative H&E and Prussian blue staining of SIM-treated SCD mouse spleens. (D) The iron content in spleen tissue was quantified. (E) Representative H&E and Prussian blue staining of SIM-treated SCD mouse livers. (F) The iron content in liver tissues was quantified. Scale bar, 50 μm. One-way ANOVA with Bonferroni’s multiple comparison tests was used for statistical analysis (n=5). *, p < 0.05.
Figure 7.
Simvastatin reduces organ damage and iron content in preclinical SCD mice. (A-B) Representative images of the spleen (A) and liver (B) of SCD mice treated with or without simvastatin (SIM, 7.5 mg/kg daily) for 4 weeks. Quantification of spleen and liver weights are shown in the graph. (C) Representative H&E and Prussian blue staining of SIM-treated SCD mouse spleens. (D) The iron content in spleen tissue was quantified. (E) Representative H&E and Prussian blue staining of SIM-treated SCD mouse livers. (F) The iron content in liver tissues was quantified. Scale bar, 50 μm. One-way ANOVA with Bonferroni’s multiple comparison tests was used for statistical analysis (n=5). *, p < 0.05.
Figure 8.
Simvastatin mitigates reactive oxygen species (ROS) and inflammatory stress in preclinical SCD mice. Spleen CD71+ cells of SCD mice were determined the cellular NADPH andNADP+ levels and their ratio (A), ROS levels by 2′-7′-dichlorodihydrofluorescein diacetate (DCFH-DA) staining (B), and relative expression of antioxidant genes (C). (D) ChIP assay determined the association of H3K27Me3, TBP, and Pol II to the antioxidant gene loci (Nqo1, Cat, Hmox1, Gclc, and Slc7a11) in spleen CD71+ cells of SIM-treated SCD mice. (E) Heatmap representation of proinflammatory factor gene transcripts in spleen CD71+ cells of SCD mice between SIM and vehicle treatment control.. One-way ANOVA with Bonferroni’s multiple comparison tests was used for statistical analysis (n=5 mice). *, p < 0.05.
Figure 8.
Simvastatin mitigates reactive oxygen species (ROS) and inflammatory stress in preclinical SCD mice. Spleen CD71+ cells of SCD mice were determined the cellular NADPH andNADP+ levels and their ratio (A), ROS levels by 2′-7′-dichlorodihydrofluorescein diacetate (DCFH-DA) staining (B), and relative expression of antioxidant genes (C). (D) ChIP assay determined the association of H3K27Me3, TBP, and Pol II to the antioxidant gene loci (Nqo1, Cat, Hmox1, Gclc, and Slc7a11) in spleen CD71+ cells of SIM-treated SCD mice. (E) Heatmap representation of proinflammatory factor gene transcripts in spleen CD71+ cells of SCD mice between SIM and vehicle treatment control.. One-way ANOVA with Bonferroni’s multiple comparison tests was used for statistical analysis (n=5 mice). *, p < 0.05.