Next Article in Journal
Biological Evaluation and Molecular Docking Studies of Novel 1,3,4-Oxadiazole Derivatives of 4,6-Dimethyl-2-sulfanylpyridine-3-carboxamide
Previous Article in Journal
Insights into the Structure and Protein Composition of Moorella thermoacetica Spores Formed at Different Temperatures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thymosin β4 Is an Endogenous Iron Chelator and Molecular Switcher of Ferroptosis

by
Joanna I. Lachowicz
1,
Giusi Pichiri
1,*,
Marco Piludu
2,*,
Sara Fais
3,
Germano Orrù
3,
Terenzio Congiu
1,
Monica Piras
1,
Gavino Faa
1,
Daniela Fanni
1,
Gabriele Dalla Torre
4,
Xabier Lopez
4,
Kousik Chandra
5,
Kacper Szczepski
5,
Lukasz Jaremko
5,
Mitra Ghosh
5,
Abdul-Hamid Emwas
6,
Massimo Castagnola
7,8,
Mariusz Jaremko
5,*,
Ewald Hannappel
9 and
Pierpaolo Coni
1
1
Department of Medical Sciences and Public Health, University of Cagliari, Cittadella Universitaria, 09042 Monserrato, Italy
2
Department of Biomedical Sciences, University of Cagliari, Cittadella Universitaria, 09042 Monserrato, Italy
3
Department of Surgical Science, OBL Oral Biotechnology Laboratory, University of Cagliari, 09124 Cagliari, Italy
4
Kimika Fakultatea, Euskal Herriko Unibertsitatea UPV/EHU, Donostia International Physics Center (DIPC), P.K. 1072 Donostia Euskadi, 20080 San Sebastian, Spain
5
Smart-Health Initiative (SHI) and Red Sea Research Center (RSRC), Division of Biological and Environmental Sciences and Engineering (BESE), King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
6
Core Labs, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
7
Institute of Chemistry of Molecular Recognition, National Research Council (Consiglio Nazionale delle Ricerche), 00185 Rome, Italy
8
Laboratory of Proteomics and Metabolomics, IRCCS, Santa Lucia Foundation, 00143 Rome, Italy
9
Institute of Biochemistry, Friedrich-Alexander-University Erlangen-Nuremberg, 91058 Erlangen, Germany
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(1), 551; https://doi.org/10.3390/ijms23010551
Submission received: 13 December 2021 / Revised: 1 January 2022 / Accepted: 3 January 2022 / Published: 4 January 2022
(This article belongs to the Section Biochemistry)

Abstract

:
Thymosin β4 (Tβ4) was extracted forty years agofrom calf thymus. Since then, it has been identified as a G-actin binding protein involved in blood clotting, tissue regeneration, angiogenesis, and anti-inflammatory processes. Tβ4 has also been implicated in tumor metastasis and neurodegeneration. However, the precise roles and mechanism(s) of action of Tβ4 in these processes remain largely unknown, with the binding of the G-actin protein being insufficient to explain these multi-actions. Here we identify for the first time the important role of Tβ4 mechanism in ferroptosis, an iron-dependent form of cell death, which leads to neurodegeneration and somehow protects cancer cells against cell death. Specifically, we demonstrate four iron2+ and iron3+ binding regions along the peptide and show that the presence of Tβ4 in cell growing medium inhibits erastin and glutamate-induced ferroptosis in the macrophage cell line. Moreover, Tβ4 increases the expression of oxidative stress-related genes, namely BAX, hem oxygenase-1, heat shock protein 70 and thioredoxin reductase 1, which are downregulated during ferroptosis. We state the hypothesis that Tβ4 is an endogenous iron chelator and take part in iron homeostasis in the ferroptosis process. We discuss the literature data of parallel involvement of Tβ4 and ferroptosis in different human pathologies, mainly cancer and neurodegeneration. Our findings confronted with literature data show that controlled Tβ4 release could command on/off switching of ferroptosis and may provide novel therapeutic opportunities in cancer and tissue degeneration pathologies.

Graphical Abstract

1. Introduction

Ferroptosis is an iron-dependent oxidative form of cell death that is morphologically, biochemically, and genetically distinct from other non-apoptotic cell death processes [1]. The latest definition of ferroptosis was presented in 2018 by the Nomenclature Committee on Cell Death (NCCD) [2] as “a form of regulated cell death (RCD) initiated by oxidative perturbations of the intracellular microenvironment that is under constitutive control by GPX4 and can be inhibited by iron chelators and lipophilic antioxidants. Ferroptosis initiates with lipid peroxidation due to ROS generation and iron availability”. Ferroptosis leads to morphological changes with a necrotic morphotype (e.g., predominant mitochondrial shrinking, an electron-dense ultrastructure, reduced/disappeared cristae, and ruptured outer mitochondrial membrane) [2]. Ferroptotic cell death is thought to require mutation of the oncogenes HRAS, KRAS, or BRAF, and is caused by lethal iron-dependent lipid peroxidation [3,4]. To provoke cell death, ferroptosis likely requires continuous iron-dependent ROS production over a prolonged period, which can be observed in neuronal, cardiac, and hepatic cells when exposed to high iron concentrations in neurodegeneration pathology [5], ischemic conditions [6], and hepatic chromatosis [7], respectively.
Iron, a highly redox-reactive element, is required in a variety of biological processes. Cancer cells develop different mechanisms both to uptake high quantities of iron and to protect themselves against iron toxicity [8]. Activation of the ferroptosis suppresses tumor cell growth and successively leads to cell death. The controlled on/off switching of ferroptosis may provide novel therapeutic opportunities in pathologies of cancer and tissue degeneration, respectively.
Thymosin β4 (Tβ4; Scheme 1) is a G-actin binding protein with different functions in the human body, including blood clotting [9], tissue regeneration [10], angiogenesis [11], prenatal and early childhood development [12], and tumor metastasis [13]. Tβ4 is also involved in anti-inflammatory and neurodegenerative processes [14]. At the subcellular level, Tβ4 is localized in the cytoplasm [15] and/or nucleus [16] and can be translocated between compartments after induction [16,17]. Thymosin β4-overexpression is correlated with different human cancer types, e.g., pancreatic [18] and colorectal cancer [19], and Tβ4 is considered a potential molecular target for anticancer therapy. Indeed, Tβ4 gene silencing decreases stemness and invasiveness in glioblastoma [20] and suppresses proliferation and invasion of non-small cell lung cancer cells [21].
The exact mechanisms of Tβ4 action are still unknown. Recently, we hypothesized that the regulatory function of Tβ4 in numerous physiological processes is correlated with its potential essential metal-binding ability [22]. Here, we discuss for the first time the formation of a metal complex with Tβ4 and its biological importance in ferroptosis. We observed that an increased concentration of Tβ4 can effectively inhibit the ferroptosis process, a mechanism that is used by cancer cells with enhanced Tβ4 expression. The inhibition of ferroptosis by augmenting extracellular concentrations of Tβ4 can stop cell degeneration, and such a strategy could be used in therapy.

2. Results

2.1. Thymosin β4 Metal Coordination

Thymosin β4 is a 43 amino acid peptide with an acetylated N-terminus (Scheme 1). It is an acidic peptide (pI 5.1) due to the presence of 11 amino acids with carboxylic groups in the side chain, which are the potential metal binding sites [23]. Moreover, nine lysine residues in the primary sequence have a high affinity to transition metal ions [23]. Thymosin β4 has no secondary structure in pure water, while in aqueous solution with the addition of TFE (or HFP), two alpha-helixes are formed in the 5–17 and 30–39 regions [24].
To our knowledge, no experimental or molecular dynamics study of the formation of essential metal complexes has been performed with Tβ4, despite there being numerous potential metal-binding sites. In order to investigate the metal-binding sites in the Tβ4 peptide, we chose nuclear magnetic resonance (NMR) techniques, which deliver valuable structural data of metal complexes [25].
The interaction of iron ions with Tβ4 in water was investigated in a residue-specific manner by monitoring of 2D 13C-1H HSQC (heteronuclear single quantum coherence) NMR (Figures S1–S3, S5 and S6), focused on the aliphatic region. The overlay of the HSQC data (Figures S2, S3, S5 and S6) showed that metal/peptide interactions are specific while formed in specific regions of the peptide, but there is no structural change upon interaction with iron2+ nor iron3+ ions.
The presence of metal ions in the peptide solution leads to the lowering of the peptide signals. The comparative analysis of the peptide signal pattern with and without metal ions permits the isolation of amino acids affected by the interactions with the metal ions (Figures S4 and S7). In the thymosin solution containing iron2+ and iron3+ ions, eighteen and fifteen residues, respectively, showed a major drop in intensity (Table 1). Since iron3+ is paramagnetic, we observed greater overall signal broadening of the residues, with the aforementioned residues showing major drops in intensity.
Circular dichroism (CD) structural analysis of free Tβ4 in water showed that the secondary structure is mainly unorganized, and only 20% could be associated with the helical form. The addition of iron ions to the peptide solution did not significantly change the secondary structure (Figure S9).
In general, NMR data precision is limited by the presence of paramagnetic species and does not permit detailed coordination analysis with iron ions. For this purpose, Al3+ ions, which form similar stoichiometry and geometry complexes as Fe3+, were used to model Fe3+/Tβ4 complexes.
The analysis of the overlay HSQC data (Figure 1) showed that metal/peptide interactions are specific while formed in specific regions of the peptide, but there were no significant changes in the secondary structure of thymosin upon interaction with Al3+ ions.
The chemical shift perturbation (CSP) was calculated considering the free peptide form and 1:5 (peptide:metal) as the maximum stoichiometry leading to the reduction of the peptide’s peak intensities and no further chemical shift changes. Only 48 isolated peaks were considered for calculating the CSP and intensity analysis. Analysis of the CSP plot (Figure 2) showed that Al3+ ions interact with N-terminal residues of Tβ4. The residues that most likely interact with the metal ion are listed in Table 1.
Figure 2 shows that residues Ile9 and Gln35 presented noticeable changes upon interaction with Al3+, however, the CSP value is too small to define the metal/peptide interaction as strong and with high affinity. Moreover, the changes in Ile9 and Gln35 signals could be associated with peptide bending upon metal coordination by adjacent residues (with carboxylic groups in amino acids side chains), rather than participation in the metal coordination core. The plot of signal intensities as a function of the Tβ4:metal ratio showed a consistent decrease for most of the residues in 1D (Figure S12) and in 2D (Figure S13). However, it is important to note that the intensity of 1D spectra is directly proportional only to the population, while the 2D spectra intensity is modulated by many different factors.
The diffusion NMR data at the different metallopeptide ratios showed the presence of only one species with a diffusion constant 1.25 × 10 10 m2/s, thus, the formation of the peptide aggregates upon metal binding can be excluded.
The NMR data analysis showed that metal ions interact with numerous amino acid residues, and four metal-binding sites could be distinguished within the Tβ4 structure (namely I–IV, Scheme 1, Table 1). The binding modes within four metal binding sites were modeled with the density functional theory (DFT) calculation.
Tβ4 contains 11 acidic residues: eight glutamate and three aspartate residues [27]. According to Pearson’s hard and soft acids and bases (HSAB) principle, hard metal ions (Fe3+ and Al3+) show a preference for binding hard bases, such as the negatively charged carboxylate side chains of glutamate and aspartate amino acids. Although there can be other donors to hard metal ions, such as backbone carbonyl and amide groups of proteins, by virtue of the partial covalent nature of the coordination bond, these donors are less competitive. Indeed, we recently demonstrated, from a thermodynamic perspective, that the preferred binding sites of Al3+ in proteins are the negatively charged side chains of amino acids rather than the backbone carbonyl and amide groups [28].
The chemical shift perturbation of the NMR data showed that the main changes in chemical shifts upon Al3+ addition arise in two different areas of the protein: within the N-terminal segment (metal-binding sites I and II) and the C-terminal segment (metal-binding sites III and IV).
The conformational analyses were prepared in order to build 3D models of the Al3+-Tβ4 complex and obtain some reliable initial configurations of this hard metal ion bound to Tβ4. According to the literature, Tβ4 in water solution exists mainly in a disordered state and it is unclear whether secondary structure elements such as extended helices in the N- and C-terminus of the protein arise upon interaction with other partners such as actin [22]. Thus, in the starting structure, the Tβ4 sequence was considered “unfolded”, without main secondary structure elements (helices and sheets).
Three major binding modes for the Al3+-Tβ4 system were obtained (Table 2): first, the metal interacting with both the N- and C-terminus of the protein (Tβ4N-C); second with the metal-bound to the N-terminal and in the middle of the protein (Tβ4N-mid); third with the metal-bound only in the N-terminus of the protein (whereas the rest of the sequence remains highly disordered; Tβ4N-N). Interestingly, all the three binding modes featured four negatively charged residues (side chains of either ASP or GLU) with a net charge of −4 surrounding the metal center, a situation that is coherent with the one proposed by Dudev and Lim [29], who stated that the optimal number of acetate-type donors in a metal-binding site is equal to q + 1, where q is the charge of the metal when first coordination shell residues are taken into account.
In Tβ4N-C binding mode, the metal is bound to GLU8 and GLU21 in a bidentate fashion, and to ASP5 and GLU35 by a single oxygen donor (Figure 3A). The six oxygen donors coordinate Al3+ through an octahedral geometry, with O-Al-O angles ≥150 degrees. During the simulation, after 3 μs, the carboxylate side chain of GLU35 coordinates the metal with both oxygen donors, whereas the side chain of GLU21 switches to a monodentate binding mode. The original state is then restored after 4 μs up to the end of the simulation. Overall, the TB4N-N protein-metal complex remains stable for the entire course of the simulation, and no other residues nor water molecules enter the first coordination shell of the metal.
The Tβ4N-mid binding mode is characterized by Al3+ ions coordinated to the side chains of ASP2 and ASP5 in a monodentate fashion and to the side chains of GLU10 and GLU24 in a bidentate fashion (Figure 3B). Moreover, in this case, there is no change in the coordination pattern of Tβ4N-mid during the entire simulation time of 5 μs.
Finally, the Tβ4N-N binding mode displayed an octahedral geometry, in which the first coordination shell of the metal is filled by one oxygen donor from the side chains of ASP5 and GLU8 and two oxygen donors from the side chains of ASP2 and GLU10 (Figure 3C). The metal center maintains this binding mode for the whole simulation.
Structural characterization revealed that Tβ4N-N displays the highest degree of disorder less compaction rate with the average root mean square deviation (RMSD) and radius of gyration (Rg) of 1.07 and 1.29, respectively, compared to Tβ4N-mid (0.89 and 1.11, respectively) and Tβ4N-C (0.92 and 1.15, respectively) as reported in Table S1. This is coherent with the fact that in the case of Tβ4N-mid and Tβ4N-C, the N- and C- termini are closer due to their respective metal binding modes that involve residues either in the medial or distal region of the protein, as also highlighted by their shorter end to end distance compared to TB4N-C (Table S1). The average amount of intramolecular hydrogen bonds are similar for Tβ4N-C and Tβ4N-N and slightly higher for Tβ4N-mid (Table S1); most populated H-bonds and salt bridges for each binding mode are reported in Table 2.
As can be deduced from Table S2, all three metal-Tβ4 complexes exist mainly as random coil structures; the prevalent structural motifs are bends and turns, consistent with the disordered nature of Tβ4. Interestingly, Tβ4N-C and Tβ4N-mid display a similar β-sheet content (16.8 and 17.6%, respectively), whereas Tβ4N-N shows a much lower content of this structural motif (2.5%). On the other hand, the latter binding mode displays a higher helicity propensity compared with the other two, in particular, it features 94.1% of the 310-helix motif. Such a situation could be due to the fact that the metal-binding activity of Tβ4N-N is confined to the N-terminal segment of the protein, while the remaining sequence is highly disordered and may increase the propensity to form helical patterns in its C-terminal regions.
All three binding modes are quite stable across the 5 μs time scale. We also observed that Tβ4 binding to aluminum could affect the overall secondary structure pattern of the thymosin. Since both glutamate and aspartate residues are the preferred acidic groups of hard metals, the three binding modes characterized herein could fit either Al3+ or Fe3+ metal ions.

2.2. Thymosin β4 and Ferroptosis

Macrophages, together with hepatocytes, play key roles in iron metabolism by mediating iron storage and recycling [30]. In hereditary hemochromatosis (HH), a representative disease that causes iron overload, hepatocytes are the first site of iron accumulation [31,32]. Recent studies [33,34] have shown that pathological iron overload in mouse macrophage cell line J774 induces ferroptosis and subsequent cell destruction. Therefore, we used J774 cells as a model to study the importance of macrophages in iron metabolism [35] and their consequent susceptibility to ferroptosis.
Ferroptosis was primarily characterized by condensed mitochondrial membrane densities and smaller mitochondria, as well as diminished or absent mitochondria crista and ruptured outer membrane [36]. Recent studies divided mitochondrial morphological alterations into three categories according to whether the mitochondria were fragmented or not, and whether they accumulated around the nucleus. Cells with a network of elongated mitochondria were assigned to category I, uniform-distributed fragmented mitochondria in cells were classified as category II, and cells with fragmented mitochondria mainly accumulating around the nucleus were classified as category III [36,37,38,39]. Mitochondrial morphological changes and metabolic regulation, especially the energetic, iron and aliphatic acid metabolism involved in ferroptosis, were extensively described in a recent review [40].
For a more detailed characterization of glutamate- and erastin-induced ferroptosis in J774 cells, including Tβ4 inhibition, we evaluated the morphological changes by TEM microscopy. The results obtained from this ultrastructural analysis are summarized in Figure 4 and Figure S19, showing cells from different experimental groups. In the control group (Figure 4A,B), the J774 cells appeared normal at the ultrastructural level. They were mainly characterized by a well-developed Golgi apparatus, rough endoplasmic reticulum, and the presence of numerous mitochondria. The nuclear compartment was separated from the cytoplasm by the double membrane of the nuclear envelope, and the cells appeared to be in interphase due to the presence of both heterochromatin and euchromatin. The presence of densely packed chromatin and the absence of the nuclear envelope represented the main features of mitotic cells, which were often observed in control samples (Figure S19).
Administration of Tβ4 (Figure 4C,D) did not appear to alter cell morphology: samples displayed structural integrity of all their cellular compartments and organelles, including mitochondria.
Following glutamate (Figure 4E,F) and erastin (Figure 4I,J) treatments, the cellular morphology of the J774 cells was dramatically altered. Specific ultrastructural alterations were mainly observed in the cytoplasmic compartment, where numerous mitochondria were found accumulating around the nucleus. These cellular organelles appeared drastically altered, showing mitochondrial membrane rupture and evident cristae disorganization. Moreover, the presence of lysosomal activity and disruption of the Golgi apparatus and rough endoplasmic reticulum indicated a poor cellular condition. However, none of these morphological alterations were related to an apoptotic form of cell death. Indeed, during our TEM investigation, no nuclear alterations were detected, and the nucleus appeared surrounded by a continuous nuclear envelope, displayed an evident nucleolus and well-organized chromatin, and no chromatin condensation or margination was detected in the J774 cells after either glutamate or erastin treatments. Moreover, both erastin and glutamate treatment seemed to induce cell detaching and a rounded cellular shape (Figure 4E). Their cell surface appeared continuous without any signs of plasma membrane rupture or disorganization, although we did observe thin protrusions. Most of these appeared to be continuous tunneling nanotubes connecting neighboring cells (Figure 4E, asterisk) previously described as structures involved in intercellular communication [41,42].
Administration of Tβ4 in glutamate (Figure 4G,H) and erastin (Figure 4K,L) treated samples seemed to gradually reverse the alterations described above. The cytoplasmic compartment showed reorganization of its cellular structures associated with an evident diminishing of lysosomal activity. Mitochondria showed partially reconstituted cristae and a continuous mitochondrial membrane. The rounded cell shape transitioned to a more flattened type.
Interestingly, in the glutamate and erastin-treated cells, we frequently detected the presence of tunneling nanotubes (TNT) connecting neighboring cells (Figure 4E, asterisk). TNT are considered dynamic structures whose formation and organization can be affected by different pathological conditions. It has been reported that oxidative stress may induce TNT formation in astrocytes and hippocampal neurons [43]. Inflammation and even viral infections have been observed to spread via tunneling nanotubes. They may also be involved in the cell-to-cell transfer of signal transduction molecules. Fas ligand (FasL) involved in programmed cell death has been demonstrated to stimulate TNT formation, promoting the propagation of cell death signals to neighboring cells [44]. Our results suggest a possible involvement of the ferroptotic process in the cell surface reorganization and TNT formation, which may facilitate diffusion of the ferroptosis process between connected cells. Considering that tunneling nanotubes are F-actin-containing protrusions of the cell surface, Tβ4 could modulate ferroptosis not only by its metal binding sites, as described above, but also through its G-actin binding properties. Tβ4, as a major G-actin sequestering protein, could prevent tunnelling nanotube formation in erastin and glutamate treated J774 cells inhibiting propagation of ferroptotic cell death signals to the connected cells.
These morphological results clearly show that glutamate and erastin induce ferroptosis in J774 cells, while Tβ4 administration after 4 h can inhibit the process and reverse the critical cellular condition.
In the next experiment, we examined the inhibitory activity of Tβ4 and compared it with the activity of Ferrostatin-1 (Fer-1), a known ferroptosis (erastin-induced) inhibitor previously described by Dixon et al. [3]. The activity of both inhibitors was measured by the cell viability assay (Figure 5) in the cell cultures incubated with erastin and increasing concentration of inhibitors. The dose-response effect of Tβ4 (Figure 5A) can be compared with the effect of Fer-1 in the same experimental conditions.
Fer-1 is a synthetic compound isolated by high-throughput screening of small molecule libraries as a potent inhibitor of ferroptosis [3]. Fer-1 prevents erastin-induced accumulation of cytosolic and lipid ROS [3]. Recent electrochemical studies showed that Fer-1 forms a complex with iron ions, supporting the theory that the anti-ferroptotic activity of Fer-1 is actually due to the scavenging of initiating alkoxyl radicals produced by ferrous iron from lipid hydroperoxides [45]. Indeed, Fer-1 derivatives substituting both amine nitrogen atoms, thus iron-binding sites, lack inhibitory activity [3]. The same inhibitory results (Figure 5), obtained for Fer-1 and Tβ4, suggest the same inhibitory mechanism through complex formation with iron ions. In addition, numerous studies showed that inhibitory activity of Fer-1 is correlated with the gene regulation of proteins involved in cellular stress processes [46]. For this reason, we examined the gene expression of selected proteins involved in stress-oxidative processes. In the next experiments, we focused on the inhibitory mechanisms of Tβ4 rather than Fer-1, while the activity of the latter is described elsewhere.
According to Dixon [1], ferroptosis could be defined as a sabotage mechanism where the cell inhibits its own normal processes. However, ferroptosis is related to oxidative stress, and consequently, the involvement of oxidative stress proteins is expected. Indeed, among thirty different genes studied in relation to ferroptosis [47], many are involved in oxidative stress. In our study, we chose four genes, namely BAX (Bcl-2 associated X-protein), HSP-70 (heat shock protein 07), which are involved in the oxidation process, but are poorly described in connection with ferroptosis; and HO-1 (hem oxygenase) and TXNRD-1 (Thioredoxin reductase 1) which were recently linked to ferroptosis [46].
Bax is a cytosol protein, which translocates into mitochondria when apoptosis occurs [48] (Scheme 2). Jungas et al. [49] studied glutamate toxicity and regulation in neuronal development. Their studies showed BAX-activation mechanism dependent on glutathione, with the downstream effect of glutathione depletion leading to the direct activation of BAX. Moreover, Bax and its transcript increase significantly in cerebral ischemia-injured neurons [50,51], suggesting that glutamate could be involved in the mechanism of BAX expression control. Of note, Tβ4 reduces the intracellular ROS levels induced by H2O2, decreases Bax/BCl2 ratio, and increases the expression of antiapoptotic proteins in cardiomyocytes [52] and cardiac fibroblasts [53].
Heat shock protein 70 (Hsp70) is a stress-responsive protein increasing under stress conditions, namely heat, hypoxia, and glucose deprivation. Increased Hsp70 concentration leads to higher stress tolerance and cytoprotection [54]. Hsp70–mediated cytoprotective effects were studied in vitro in numerous cell/tissue types under changing environmental conditions of stress. Noteworthily, Hsp70 has significant cytoprotective effects under ischemic conditions. For instance, induction of Hsp70 (e.g., by enforced overexpression with viral vectors or drugs) protects brain cells against ischemic injuries [55,56,57,58]. Oxidative stress and antioxidants can influence Hsp70 expression [59,60], while the reduction in Hsp70 expression leads to higher ROS generation and mitochondrial protein oxidation [61]. Guo et al. [54] showed that Hsp70 regulates cellular redox status by changing the activities of the GSH-related enzymes, glutathione peroxidase, and glutathione reductase, in response to hypoxic and ischemic stress. The switching of antioxidant enzyme activities could be a critical mechanism mediating the enhanced cytoprotection afforded by the enhanced expression of Hsp70. Despite numerous studies and different theories, the exact molecular mechanisms of Hsp70 action in ischemic conditions remains undefined
Heme oxygenase (HO)-1 metabolizes heme into biliverdin/bilirubin, carbon monoxide, and ferrous iron, and has been suggested to demonstrate cytoprotective effects under various stress-related conditions. HO-1 is commonly regarded as a survival molecule, exerting an important role in cancer progression, and its inhibition is considered beneficial in several cancers [62]. HO-1 is often upregulated in tumor tissues, and its expression is further increased in response to therapies [63]. The role of HO-1 in ferroptosis is extensively described in the recent review [62]. Depletion of cellular glutathione has been shown to increase HO-1 gene transcription in the mouse motor neuron-like hybrid cells, NSC34 cells [64]. The correlation between HO-1 expression and ferroptosis is unclear. In HT-1080 fibrosarcoma cells, Erastin induces a time- and dose-dependent increase of HO-1 expression [65]. However, HO-1 also functions as a negative regulator in erastin- and sorafenib-induced hepatocellular carcinoma, since knockdown of HO-1 expression enhanced cell growth inhibition by erastin and sorafenib [66].
The selenoprotein thioredoxin reductase 1 (TXNRD1 gene) is a cytoplasmic protein that decreases thioredoxin levels and protects against oxidative stress [67]. Recently, it was shown that TXNRD1 is a strong negative modulator of ferroptosis susceptibility in pancreatic cancer cells [68]. Moreover, a high-dose (25 mg/kg) of the anti-rheumatoid arthritis drug auranofin (AUR) upregulates hepcidin expression and induces ferroptosis, and causes lipid peroxidation through inhibition of thioredoxin reductase (TXNRD) activity, in C57BL/6J mice and a mouse model of hemochromatosis (Hfe−/− mice) [69].
Here, we evaluated the levels of BAX, HO-1, Hsp70, and TRNRD-1 mRNAs (Figure 6) under normal and ferroptosis conditions, which were induced by the extracellular overload of glutamate. In addition, we evaluated mRNA expression levels of these proteins in the presence of high extracellular concentrations of thymosin. Finally, we examined their mRNA levels in the cells treated with thymosin after 4 h of ferroptosis induction. We observed that glutamate overload led to the downregulation of all studied genes; conversely, the presence of Tβ4 upregulated mRNA expression, particularly of HO-1. Notably, the downregulation of glutamate could be inverted by the successive treatment with Tβ4.

3. Discussion

Thymosin β4 is highly conserved across the animal kingdom, from amphibians to mammals, and its amino acid sequence is well conserved [15]. Among mammalian species, the amino acid sequences of human, rat, murine, bovine, porcine, and ovine Tβ4 are identical [70]. Thymosin β4 is encoded by the TMSB4X gene. Mouse, rat, zebrafish, frog, and pig display differential mRNA expression of TMSB4X in organs, while in humans, it is widely distributed in tissues with high concentrations of different essential metal ions [71] (Scheme S1).
Iron is an essential metal ion with a fundamental role in diverse cellular processes, such as DNA synthesis, proliferation, cell cycle regulation, and the function of proteins containing iron–sulfur clusters [72]. Iron–sulfur cluster-containing proteins include enzymes that contribute to maintaining genomic stability, as well as respiratory function [72]. Indeed, rapidly reproducing cancer cells re-program iron metabolism in ways that result in net iron influx [8,72]. Iron is a highly reactive metal due to its redox potential, and when present in a free ion form, it contributes to the production of reactive oxygen species (ROS) in the Fenton reaction (Scheme 2). Oxidative and iron-dependent cell death with neurodegenerative consequences was previously described in epilepsy, stroke, and other trauma situations [73,74,75]. Importantly, increased levels of iron were found in the central nervous system of patients with neurodegeneration [5]. Surprisingly, elevated iron concentrations in cancer cells do not lead to cell disruption, but the exact mechanism of cell defense remains unknown. Cancer cells have developed many different iron-regulating mechanisms, which permit them to benefit from high iron concentrations without suffering from the harmful effects of free radicals. Many iron-regulating processes in tumor pathology have already been discovered [76], while others are unknown. Considering the iron-dependent etiology of ferroptosis, it is likely that some of the unknown iron-defense processes in tumor cells can inhibit ferroptosis-induced cell death.
The involvement of Tβ4 in different cancer types has been widely described in recent years (Table S4). It is noteworthy that among the 13 most common cancer types (https://www.healthline.com/health/most-common-cancers, accessed date: 1 February 2021), the nine tissues with the highest cancer incidence also have high thymosin mRNA expression (Scheme S1, Table S4), which is further enhanced in tumors (Table S4). Recently, it was shown that ferroptosis is fundamental in anticancer therapy, and efficient anticancer drugs can induce ferroptosis-cell death (Table S4).
Neurodegeneration is characterized by neuronal cell death. Detailed biochemical and morphological studies of neuronal tissues have linked neurodegeneration with increased iron concentrations in the tissue and successive production of free radicals [77,78]. Since ferroptosis was first described in 2012 [3], numerous research data have correlated neurodegenerative processes with iron-ROS induced cell death (Table S4) and glutamate-induced cell injury [79]. Up to now, Tβ4 has never been linked to ferroptosis, but different experimental studies have shown the beneficial effects of Tβ4 in neurodegeneration (Table S4).
Ferroptosis can be induced in vivo in rat hippocampal slices and primary oligodendrocyte models as well as in ischemia-reperfusion injury models, by high concentrations of extracellular iron or glutamate, or by depletion of extracellular cysteine. In vitro, ferroptosis can be induced by physiological conditions (e.g., high extracellular iron or glutamate) or small molecules (e.g., erastin), which blocks system Xc-mediated Cys2 import [1] (heterodimer composed of solute carrier family 7, member 11 (SLC7A11) and SLC3A2; Scheme 2), while glutamate receptors are not involved in ferroptosis [80]. Erastin and glutamate are termed type I ferroptosis inducers that inhibit cysteine uptake by the cysteine/glutamate antiporter (system xc), enabling glutathione synthesis (Scheme 2), and thereby lowering antioxidant defenses [3].
Glutamate-induced oxidative stress leads to lipid peroxidation-mediated ferroptosis (Scheme 2) and is a major contributor to neurodegenerative pathologies [81]. In recent studies, L-Glutamate (3 h, 5 mM) mimicked in vitro the consequences of stroke and neurodegenerative pathologies in cell cultures [82,83]. The addition of iron enhances glutamate-induced cell death [81], while inhibition of ferroptosis prevents glutamate-induced cell death in organotypic hippocampal slice cultures [3]. Importantly, glutamate activity can be inhibited by iron chelators (e.g., Desferal) [3].
Ferroptosis, in contrast to other RCDs processes can be inhibited, thus inhibiting molecules, namely iron chelators and lipophilic antioxidants, that could be used as a molecular switcher of the entire process. Our metal coordination studies showed that Tβ4 could bind both Fe2+ and Fe3+ metal ions in vitro in four different metal binging regions (Figure 3). Iron chelating ability of Tβ4 takes a part of iron cellular homeostasis, and increased concentration of Tβ4 leads to upregulation of HO-1 protein expression (Figure 6B), which metabolizes heme into ferrous iron.
Tβ4 as endogenous iron chelator influences anti-oxidative processes related to free iron ions and ROS production in Fenton reaction. Recently, Zhu et al. [84] showed that Tβ4 reduces oxidative stress and lipid peroxidation in the injured liver model in mice by inhibiting ferroptosis via up-regulation of GPX4 expression. In our experimental conditions, increased concentration of Tβ4 enhances the expression of oxidative stress-related genes: BAX, HSP70, and TRNRD-1, which are downregulated during ferroptosis. Noteworthily, the expression of these genes is even higher when Tβ4 concentration increases during the ferroptosis process (Figure 6A,D). Such observation suggests that Tβ4 expression and/or supplementation could be used to enhance anti-oxidative processes and further inhibit ferroptosis.
In our experimental conditions, erastin- and glutamate-induced ferroptosis leads to different ultrastructural alternations of cell morphology (Figure 4E,F,I,J). Some of them, e.g., those involving mitochondria, were previously described in ferroptosis processes, while the others—disruption of Golgi apparatus and rough endoplasmic reticulum, continuous tunneling nanotubes connecting neighboring cells—are less investigated in this context and need further studies. Importantly, Tβ4 administration during ferroptosis can gradually reverse the morphological alternations (Figure 4G,H,K,L).

4. Conclusions

Thymosin β4 is an iron chelator, which takes part in the iron homeostasis mechanism and iron-related oxidative stress. Administration of Tβ4 enhances the expression of oxidative stress genes (HO-1, BAX, HSP70, and TRNRD-1), which are suppressed by glutamate-induced ferroptosis. The characteristic morphological alternations of the ferroptosis process (e.g., predominant mitochondrial shrinking, an electron-dense ultrastructure, reduced/disappeared cristae, and ruptured outer mitochondrial membrane) can be reversed by the administration of Tβ4, which restores normal cytoplasmic compartment organization. In this context, Tβ4 can be considered a ferroptosis switcher. A better understanding of Tβ4 pathological functioning in cancer cells, and its beneficial effects in neurodegeneration, will enable us to switch different pathological states into physiological processes.

5. Experimental Procedures

5.1. Reagents

Human Thymosin β4 (HPLC > 98%) was purchased from Regene RX, Rockville, MA, USA. Erastin, L-glutamate, BODIPY-C11, formaldehyde, glutaraldehyde in 0.1 M sodium cacodylate buffer pH 7.4, osmium tetroxide, uranyl acetate, ethanol, and xylene were Sigma Aldrich products; purity 99%. All high-purity grade reagents were used without further purification.

5.2. EPR Experimental Details

All EPR spectra were recorded using a continuous wave Bruker ELEXSYS E500 electron paramagnetic resonance (EPR) spectrometer operating at the X band equipped with a superhigh Q (ER 4122 SHQ) resonator. To create comparable data, the spectra were measured using the same conditions with instrumental parameters as described previously [85,86]. Bruker Xenon software (Bruker BioSpin, Rheinstetten, Germany) was used for both data collection and spectral processing.

5.3. NMR

NMR samples were prepared by dissolving 1.0 mg of Tβ4 in 500 μL of D2O. The pD was adjusted by adding DCl or NaOD and calculated as pD = pH (pH meter reading) + 0.4 [87]. All NMR experiments were performed on 800 MHz spectrometers equipped with a triple resonance TCI probe and Z gradients coils. The temperature was kept at 25 °C and it was calibrated with an external standard mixture of ethylene glycol and DMSO-D6 by monitoring the chemical shift difference between the hydroxyl and methylene proton. The chemical shifts were calibrated with respect to standard DSS for proton, and 13C chemical shifts were calibrated indirectly. D2O was used for locking.
Apart from the one-dimensional spectrum with water suppression using excitation sculpting, two-dimensional homonuclear TOCSY, NOESY, hetero-nuclear HSQC and HSQC-TOCSY (with mixing 16 ms) experiments were recorded. For the heteronuclear experiments, 2048 and 360 data points were collected in proton and carbon dimensions, respectively, which corresponds to 125 ms and 17 ms acquisition time in the respective dimension. The data were apodized with sine function and zero filled with 4096 and 1024 points in proton and carbon dimensions, respectively. The recycle delay was kept for 1 and 1.5 s, respectively, for HSQC and TOCSY-HSQC. The total measurement time was 300 and 240 min for HSQC and TOCSY-HSQC, respectively. The sample was stable during measurement.
Five different samples were prepared with the addition of individually FeCl2 and FeCl3 with thymosin: metal ion ratios 1:0.05, 1:0.1, 1:0.2, 1:0.5, and 1:1; and five different samples were prepared with the addition of AlCl3 with Tβ4: metal ion ratios 1:0, 1:1, 1:2, 1:5 and 1:10. For the sample with some AlCl3, slightly visible precipitation was noticed.
The diffusion was measured using the pulsed field gradient (PFG) NMR on a Bruker Avance 800 MHz spectrometer equipped with a triple resonance TCI probe and Z gradients coils. The pulse sequence incorporates longitudinal eddy–current delay (LED) along with bipolar gradients. All the diffusion measurements were recorded in a pseudo-two-dimensional manner. In the process, we acquired 24 spectra over the gradient strengths ranging from 2% to 95% with a linear variation of the gradient strength. The diffusion time was kept 175 ms along with 32 scans and 1.5 s recycle delay. The decaying signal intensity was analyzed using the following equation,
I = I0 * exp{−Dγ2g2δ2(Δ − δ/3)}
I is the observed intensity, I0 is the reference intensity, D is the diffusion coefficient, γ is the gyromagnetic ratio of the observed nucleus, g is the gradient strength, δ is the length of the gradient, and ∆ is the diffusion time.
All NMR data were processed using Topspin 4.0.7 and analyzed with Sparky software. Nearly complete assignments of the 1H and 13C resonances for all the residues were carried out by application of a standard and well-established procedure using HSQC and TOCSY-HSQC [88]. We also used the PROSECCO (PROtein SEquences and Chemical shift Correlations) server, which predicts the chemical shift of both folded and unfolded proteins, which is sequence-based. Since Tβ4 is intrinsically disordered in D2O, we used the server accordingly. The earlier assignment of the apo form of Tβ4 was transferable.

5.4. DFT

System Setup and Conformational Analysis

Conformational analyses were carried out on the full-length N-acetylated thymosin β-4. The protein sequence was built with Maestro v. 12.5 as a random coil (no secondary structure elements) amino acid sequence, as this is the most likely Tβ4 structure in solution [22]. An acetyl group is added to the N-terminal serine residue, whereas the C-terminal serine residue is taken in its charged COO state, in agreement with the physiological form of Tβ4 [22]. Therefore, the initial system is composed of the random coil protein and one Fe3+ metal ion placed nearby the protein. The system was modeled with the OPLS-AA force field [89] as implemented in the MacroModel v12.4 package. Choice of the metal ion is dictated by the absence of Al3+ parameters in this force field. The MacroModel v12.4 [90] module was used along with the MonteCarlo Multiple Minimum (MCMM) algorithm [91]. Three different conformational analyses were performed; in the first one, the number of search iterations was set to 50,000 steps and the Polak–Ribiere conjugate gradient (PRCG) minimization algorithm was used [92] set to 10,000 minimization steps and the gradient of first derivatives convergence threshold set to 0.005. The second and third conformational analyses were run for a total of 10,000 search steps using truncated-Newton conjugate gradient minimization algorithm (TNCG) [93], with 20,000 minimization steps and the gradient of second derivatives convergence threshold set to 0.05. All torsional variables were turned on during the conformational search, with a maximum atom deviation cutoff of 0.8 Å to discriminate between redundant conformers. The energy window for saving structures was set to 20.0 kcal/mol.
Through visual inspection of all the conformers produced by the three conformational runs, we selected 11 structures whose metal-binding mode lies in one of the four metal anchoring areas predicted by NMR and CSP measurements (Figure 1). Three binding modes were identified: Tβ4N-N, in which the metal was bound to the N-terminal of the protein (ASP2, ASP5, GLU8 and GLU10 residues); Tβ4N-mid in which the metal was bound to the N-terminal and the middle of the protein (ASP2, ASP5, GLU10 and GLU24 residues); Tβ4N-C in which the metal was bound to the N-terminal and the C-terminal regions of the protein (ASP5, GLU8, GLU21 and GLU35 residues).
Five conformers for the Tβ4N-C and Tβ4N-mid binding modes were chosen, and one for the Tβ4N-N binding mode, for a total of eleven initial structures.

5.5. Molecular Dynamics Simulations

The eleven Fe3+-Tβ4 structures selected from the previous conformational analyses were immersed in a dodecahedron box with the solute kept to a distance of 20 Å from the box and solvated with ≃ 8000 TIP4P-ew [94] water molecules along with a salt concentration of 150 mM NaCl. The protein was modeled with the AMBER99SB*-ILDN-Q force field [95,96], suitable for intrinsically disordered proteins. Although conformational analyses have been performed with Fe3+ ion, MD simulations were performed using the Al3+ metal ion. This choice is due to the fact iron is a paramagnetic element, therefore, NMR and CSP data for Al(III) (Table 1) are more reliable. Parameters for the Al3+ ion were taken from Pengfei Li et al. [97] as they have been developed for the AMBER family of force fields and the TIP4P-ew water model.
Long-range electrostatic interactions were evaluated using a particle-mesh Ewald (PME) approach [98], and Lennard-Jones interactions were truncated at 12~Å with an atom-based force switching function which is effective at 10 Å. After 5000 steps of steepest descent minimization, the system was equilibrated in a multistep manner with decreasing force constants according to the scheme illustrated in Table S3. The production run was conducted in the NPT ensemble for 1 μs for each run. All atoms were constrained with the LINCS algorithm [99]; hydrogen atoms were treated as virtual interaction sites so that to allow a time step of 4.0 fs while maintaining energy conservation [100]. Simulations were performed at the temperature of 310 K and 1.0 atm pressure. The V-rescale thermostat (JCP 126, 014101) with a time constant of 0.1 ps was used to control the temperature. The isotropic Berendsen barostat was used for the equilibration stage, whereas the Parrinello-Rahman one [101] was used for production runs; the reference pressure was set to 1.0 atm, the time constant to 2.0 ps and the isothermal compressibility to 4.5 × 10−5 in order to maintain the pressure of the system. Dispersion correction was used to correct long range dispersion of energy and pressure.

5.6. Cell Cultures

Commercial cell line J774 (ICLC ATL98011) was obtained from the Istituto Nazionale per la Ricerca sul Cancro c/o CBA (ICLC, Genova). The culture medium used for this purpose was a mixture of MEM (EBSS), 10% fetal bovine serum (FBS), 100 units/mL penicillin, 100 mg/mL streptomycin, 2 mM L-Glutamine, 1% non-essential amino acids. To perform different experimental conditions, confluent cells were isolated using trypsin/EDTA and, for the experimental procedure, samples of 2–3 × 104 cells/cm2 J774 cells were plated on different glass coverslips at 37 °C, 5% CO2. The daily prepared water solution of Tβ4 was added to the cell growing medium to the final 10 μM concentration.

5.7. Cell Viability Measurements

The cell viability was evaluated using the LunaFL Cell Counter (www.logosbio.com, accessed date: 2 March 2021) as indicated in the manufacturer’s user manual. Cell viability under test conditions is reported as a percentage relative to the negative control treatment.
For each of three analyses (repeated experiments), three distinct biological replicas were made, and quantitative data were expressed as a mean. Cell viability quantitative data were expressed as mean + standard error. Statistical comparison of the cell viability quantitative data in different experimental conditions were performed, with 5.0 GraphPad Prism software, by the unpaired t-test. The data were also analysed with one-way ANOVA. p values < 0.05 was considered significant.

5.8. Gene Expression

The expression of the following genes at the mRNA level was evaluated by real-time PCR: BAX; HSP70; HO-1 and TXNRD-1. The cells were immediately frozen after treatment and stored at −80 °C until RNA extraction. The RNA was extracted with NorDiag ARROW RNA kit using the automated tool of the same company. Actin was used as a reference gene [102]. The Real-time reverse-transcriptase PCR was performed using the Roche Light Cycler system and SYBR Green I kit amplification kit (Roche Diagnostics). The following primers (b-actF 1⁄4 50-GCATGGGTCAGAAGG-30, b-actR 1⁄4 50-AGGCGTA-CAGGGATAG-30) were designed using the sequences of the human beta-actin mRNA (Gen-Bank accession no. NM_001101).
Real time was performed in 20 μL of final volume containing: 3 mM MgCl2, 0.25 mM of each primer, and 2 μL of RNA extract. Cycling was performed using the following amplification conditions: an initial reverse transcription at 55 °C for 10 min, denaturation at 95 °C for 30 s followed by 35 cycles at 95 °C for 10 s, 53 °C for 10 s and 72 °C for 8 s with subsequent melting analysis: heating to 95 °C for 20 s, cooling to 5 °C for 10 s and re-heating to 95 °C per second.
Fluorescence was detected at the end of the 81 °C segment in the PCR step (single mode) and at 45 °C segments in the melting step (continuous mode) in the F1 channel. The relative gene expression was analyzed by using the 2-DDCT method [103]. For each of three analyses (repeated experiments), three distinct biological replicas were made, and quantitative data were expressed as a mean. Folding change values in genes expression relative to the beta-actin have been represented as mean + standard error. Statistical comparison of the gene’s expression in different experimental conditions was performed with 5.0 GraphPad Prism software. The data were also analyzed with one-way ANOVA. p values < 0.05 was considered significant.

5.9. TEM Microscopy

Samples of J774 cells used in this study were fixed for 2 h in a solution of 1% (para)formaldehyde and 1.25% glutaraldehyde in 0.1 M sodium cacodylate buffer, pH 7.4. They were post-fixed in 1% osmium tetroxide for 1 h and stained in aqueous uranyl acetate 0.25% overnight at 4 °C. Cell cultures underwent a dehydration process in an ascending graded series of ethanol and xylene, followed by infiltration and embedding in Embed 812 resin. The specimens were transferred to flat embedding molds and polymerized in the oven at 60 °C for 24 h. 60–90 nm thin sections were cut with an LKB ultratome 8800 ultramicrotome, post-stained with uranyl acetate and bismuth subnitrate, observed and photographed using a transmission electron microscope (JEOL 1400 plus model, Tokyo, Japan) operating at 80 kV.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23010551/s1.

Author Contributions

Project idea and execution, wrote the manuscript, J.I.L.; research data of cell culturing and microscopy, G.P.; research data of TEM microscopy, M.P. (Marco Piludu), T.C. and M.P. (Monica Piras); research data of mRNA expression analysis, S.F. and G.O.; contributed to discussion and reviewed and edited the manuscript, M.P. (Marco Piludu), G.F. and D.F.; Molecular dynamics research data, G.D.T. and X.L.; research data of NMR spectroscopy, M.J., K.C., K.S. and L.J.; research data of EPR spectroscopy, M.G. and A.-H.E.; contributed to discussion and reviewed and edited the manuscript, E.H. and M.C.; research data of cell culturing and microscopy; contributed to discussion and reviewed and edited the manuscript, P.C. All authors have read and agreed to the published version of the manuscript.

Funding

Financial support from FIR 2019 and from Regione Autonoma della Sardegna (grant RASSR79857) is gratefully acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge the CeSAR (Centro Servizi d’Ateneo per la Ricerca) of the University of Cagliari, Italy for the transmission electron microscopic analysis performed with the JEOL 1400 plus TEM.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dixon, S.J. Ferroptosis: Bug or feature? Immunol. Rev. 2017, 277, 150–157. [Google Scholar] [CrossRef] [PubMed]
  2. Galluzzi, L.; Vitale, I.; Aaronson, S.A.; Abrams, J.M.; Adam, D.; Agostinis, P.; Alnemri, E.S.; Altucci, L.; Amelio, I.; Andrews, D.W.; et al. Molecular mechanisms of cell death: Recommendations of the Nomenclature Committee on Cell Death 2018. Cell Death Differ. 2018, 25, 486–541. [Google Scholar] [CrossRef]
  3. Dixon, S.J.; Lemberg, K.M.; Lamprecht, M.R.; Skouta, R.; Zaitsev, E.M.; Gleason, C.E.; Patel, D.N.; Bauer, A.J.; Cantley, A.M.; Yang, W.S.; et al. Ferroptosis: An Iron-Dependent Form of Nonapoptotic Cell Death. Cell 2012, 149, 1060–1072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Yagoda, N.; von Rechenberg, M.; Zaganjor, E.; Bauer, A.J.; Yang, W.S.; Fridman, D.J.; Wolpaw, A.J.; Smukste, I.; Peltier, J.M.; Boniface, J.J. RAS–RAF–MEK-dependent oxidative cell death involving voltage-dependent anion channels. Nature 2007, 447, 865–869. [Google Scholar] [CrossRef] [Green Version]
  5. Ndayisaba, A.; Kaindlstorfer, C.; Wenning, G.K. Iron in neurodegeneration–cause or consequence? Front. Neurosci. 2019, 13, 180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Carbonell, T.; Rama, R. Iron, oxidative stress and early neurological deterioration in ischemic stroke. Curr. Med. Chem. 2007, 14, 857–874. [Google Scholar] [CrossRef]
  7. Wang, H.; An, P.; Xie, E.; Wu, Q.; Fang, X.; Gao, H.; Zhang, Z.; Li, Y.; Wang, X.; Zhang, J.; et al. Characterization of ferroptosis in murine models of hemochromatosis. Hepatology 2017, 66, 449–465. [Google Scholar] [CrossRef]
  8. Bogdan, A.R.; Miyazawa, M.; Hashimoto, K.; Tsuji, Y. Regulators of Iron Homeostasis: New Players in Metabolism, Cell Death, and Disease. Trends Biochem. Sci. 2015, 41, 274–286. [Google Scholar] [CrossRef] [Green Version]
  9. Cassimeris, L.; Safer, D.; Nachmias, V.T.; Zigmond, S.H. Thymosin beta 4 sequesters the majority of G-actin in resting human polymorphonuclear leukocytes. J. Cell Biol. 1992, 119, 1261–1270. [Google Scholar] [CrossRef]
  10. Philp, D.; Kleinman, H.K. Animal studies with thymosin β4, a multifunctional tissue repair and regeneration peptide. Ann. N. Y. Acad. Sci. 2010, 1194, 81–86. [Google Scholar] [CrossRef] [Green Version]
  11. Smart, N.; Rossdeutsch, A.; Riley, P.R. Thymosin β4 and angiogenesis: Modes of action and therapeutic potential. Angiogenesis 2007, 10, 229–241. [Google Scholar] [CrossRef]
  12. Faa, G.; Piras, M.; Mancuso, L.; Coni, P.; Pichiri, G.; Orrù, G.; Fanni, D.; Gerosa, C.; Cao, G.; Taibi, R.; et al. Thymosin beta-4 prenatal administration improves fetal development and halts side effects due to preterm delivery. Curr. Med. Chem. 2021, 25, 431–437. [Google Scholar]
  13. Goldstein, A.L. Thymosin β4: A New Molecular Target for Antitumor Strategies; Oxford University Press: Oxford, UK, 2003. [Google Scholar]
  14. Pardon, M.-C. Anti-inflammatory potential of thymosin β4 in the central nervous system: Implications for progressive neu-rodegenerative diseases. Expert Opin. Biol. Ther. 2018, 18 (Suppl. 1), 165–169. [Google Scholar] [CrossRef]
  15. Erickson-Viitanen, S.; Ruggieri, S.; Natalini, P.; Horecker, B. Distribution of thymosin β4 in vertebrate classes. Arch. Biochem. Biophys. 1983, 221, 570–576. [Google Scholar] [CrossRef]
  16. Huff, T.; Rosorius, O.; Otto, A.M.; Müller, C.S.G.; Ballweber, E.; Hannappel, E.; Mannherz, H.G. Nuclear localisation of the G-actin sequestering peptide thymosin β4. J. Cell Sci. 2004, 117, 5333–5341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Piludu, M.; Piras, M.; Pichiri, G.; Coni, P.; Orrù, G.; Cabras, T.; Messana, I.; Faa, G.; Castagnola, M. Thymosin Beta 4 May Translocate from the Cytoplasm in to the Nucleus in HepG2 Cells following Serum Starvation. An Ultrastructural Study. PLoS ONE 2015, 10, e0119642. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, Y.; Feurino, L.W.; Zhai, Q.; Wang, H.; Fisher, W.E.; Chen, C.; Yao, Q.; Li, M. Thymosin beta 4 is overexpressed in human pancreatic cancer cells and stimulates proinflammatory cytokine secretion and JNK activation. Cancer Biol. Ther. 2008, 7, 419–423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Lee, S.Y.; Park, M.J.; Lee, H.K.; Son, H.J.; Kim, C.N.; Kim, J.H.; Kang, D.W. Increased expression of thymosin β4 is in-dependently correlated with hypoxia inducible factor-1α (HIF-1α) and worse clinical outcome in human colorectal cancer. J. Pathol. Transl. Med. 2017, 51, 9–16. [Google Scholar] [CrossRef]
  20. Wirsching, H.-G.; Krishnan, S.; Florea, A.-M.; Frei, K.; Krayenbühl, N.; Hasenbach, K.; Reifenberger, G.; Weller, M.; Tabatabai, G. Thymosin beta 4 gene silencing decreases stemness and invasiveness in glioblastoma. Brain 2014, 137, 433–448. [Google Scholar] [CrossRef] [Green Version]
  21. Huang, D.; Wang, S.; Wang, A.; Chen, X.; Zhang, H. Thymosin beta 4 silencing suppresses proliferation and invasion of non-small cell lung cancer cells by repressing Notch1 activation. Acta Biochim. Biophys. Sin. 2016, 48, 788–794. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lachowicz, J.I.; Jaremko, M.; Jaremko, L.; Pichiri, G.; Coni, P.; Piludu, M. Metal coordination of thymosin β4: Chemistry and possible implications. Coord. Chem. Rev. 2019, 396, 117–123. [Google Scholar] [CrossRef]
  23. Al-Harthi, S.; Lachowicz, J.I.; Nowakowski, M.E.; Jaremko, M.; Jaremko, Ł. Towards the functional high-resolution coor-dination chemistry of blood plasma human serum albumin. J. Inorg. Biochem. 2019, 198, 110716. [Google Scholar] [CrossRef]
  24. Zarbock, J.; Oschkinat, H.; Hannappel, E.; Kalbacher, H.; Voelter, W.; Holak, T.A. Solution Conformation of Thymosin beta4: A Nuclear Magnetic Resonance and Simulated Annealing Study. Biochemistry 1990, 29, 7814–\7821. [Google Scholar] [CrossRef]
  25. Emwas, A.-H.; Szczepski, K.; Poulson, B.G.; Chandra, K.; McKay, R.T.; Dhahri, M.; AlAhmari, F.; Jaremko, L.; Lachowicz, J.I.; Jaremko, M. NMR as a “Gold Standard” Method in Drug Design and Discovery. Molecules 2020, 25, 4597. [Google Scholar] [CrossRef]
  26. Williamson, M.P. Using chemical shift perturbation to characterise ligand binding. Prog. Nucl. Magn. Reson. Spectrosc. 2013, 73, 1–16. [Google Scholar] [CrossRef]
  27. Low, T.; Goldstein, A.L. Chemical characterization of thymosin beta 4. J. Biol. Chem. 1982, 257, 1000–1006. [Google Scholar] [CrossRef]
  28. Mujika, J.I.; Torre, G.D.; Formoso, E.; Grande-Aztatzi, R.; Grabowski, S.J.; Exley, C.; Lopez, X. Aluminum's preferential binding site in proteins: Sidechain of amino acids versus backbone interactions. J. Inorg. Biochem. 2018, 181, 111–116. [Google Scholar] [CrossRef] [PubMed]
  29. Dudev, T.; Lim, C. Metal Binding Affinity and Selectivity in Metalloproteins: Insights from Computational Studies. Annu. Rev. Biophys. 2008, 37, 97–116. [Google Scholar] [CrossRef] [PubMed]
  30. Ganz, T. Systemic iron homeostasis. Physiol. Rev. 2013, 93, 1721–1741. [Google Scholar] [CrossRef] [Green Version]
  31. Cazzola, M.; Della Porta, M.G.; Malcovati, L. Clinical Relevance of Anemia and Transfusion Iron Overload in Myelodysplastic Syndromes. Hematology 2008, 2008, 166–175. [Google Scholar] [CrossRef] [Green Version]
  32. Shenoy, N.; Vallumsetla, N.; Rachmilewitz, E.; Verma, A.; Ginzburg, Y. Impact of iron overload and potential benefit from iron chelation in low-risk myelodysplastic syndrome. Blood 2014, 124, 873–881. [Google Scholar] [CrossRef]
  33. Youssef, L.A.; Rebbaa, A.; Pampou, S.; Weisberg, S.P.; Stockwell, B.R.; Hod, E.A.; Spitalnik, S.L. Increased erythrophago-cytosis induces ferroptosis in red pulp macrophages in a mouse model of transfusion. Blood Am. J. Hematol. 2018, 131, 2581–2593. [Google Scholar]
  34. Imoto, S.; Kono, M.; Suzuki, T.; Shibuya, Y.; Sawamura, T.; Mizokoshi, Y.; Sawada, H.; Ohbuchi, A.; Saigo, K. Haemin-induced cell death in human monocytic cells is consistent with ferroptosis. Transfus. Apher. Sci. 2018, 57, 524–531. [Google Scholar] [CrossRef]
  35. Soares, M.P.; Hamza, I. Macrophages and Iron Metabolism. Immunity 2016, 44, 492–504. [Google Scholar] [CrossRef] [Green Version]
  36. Angeli, J.P.F.; Schneider, M.; Proneth, B.; Tyurina, Y.Y.; Tyurin, V.; Hammond, V.J.; Herbach, N.; Aichler, M.; Walch, A.; Eggenhofer, E.; et al. Inactivation of the ferroptosis regulator Gpx4 triggers acute renal failure in mice. Nat. Cell Biol. 2014, 16, 1180–1191. [Google Scholar] [CrossRef] [Green Version]
  37. Grohm, J.; Plesnila, N.; Culmsee, C. Bid mediates fission, membrane permeabilization and perinuclear accumulation of mitochondria as a prerequisite for oxidative neuronal cell death. Brain Behav. Immun. 2010, 24, 831–838. [Google Scholar] [CrossRef] [PubMed]
  38. Jelinek, A.; Heyder, L.; Daude, M.; Plessner, M.; Krippner, S.; Grosse, R.; Diederich, W.E.; Culmsee, C. Mitochondrial rescue prevents glutathione peroxidase-dependent ferroptosis. Free Radic. Biol. Med. 2018, 117, 45–57. [Google Scholar] [CrossRef] [PubMed]
  39. Neitemeier, S.; Jelinek, A.; Laino, V.; Hoffmann, L.; Eisenbach, I.; Eying, R.; Ganjam, G.K.; Dolga, A.; Oppermann, S.; Culmsee, C. BID links ferroptosis to mitochondrial cell death pathways. Redox Biol. 2017, 12, 558–570. [Google Scholar] [CrossRef] [PubMed]
  40. Wang, H.; Liu, C.; Zhao, Y.; Gao, G. Mitochondria regulation in ferroptosis. Eur. J. Cell Biol. 2019, 99, 151058. [Google Scholar] [CrossRef]
  41. Gerdes, H.-H.; Carvalho, R.N. Intercellular transfer mediated by tunneling nanotubes. Curr. Opin. Cell Biol. 2008, 20, 470–475. [Google Scholar] [CrossRef]
  42. Davis, D.M.; Sowinski, S. Membrane nanotubes: Dynamic long-distance connections between animal cells. Nat. Rev. Mol. Cell Biol. 2008, 9, 431–436. [Google Scholar] [CrossRef]
  43. Zhu, D.; Tan, K.S.; Zhang, X.; Sun, A.Y.; Sun, G.Y.; Lee, J.C.-M. Hydrogen peroxide alters membrane and cytoskeleton properties and increases intercellular connections in astrocytes. J. Cell Sci. 2005, 118, 3695–3703. [Google Scholar] [CrossRef] [Green Version]
  44. Luchetti, F.; Canonico, B.; Arcangeletti, M.; Guescini, M.; Cesarini, E.; Stocchi, V.; Degli Esposti, M.; Papa, S. Fas Signalling Promotes Intercellular Communication in T Cells. PLoS ONE 2012, 7, e35766. [Google Scholar] [CrossRef] [Green Version]
  45. Miotto, G.; Rossetto, M.; Di Paolo, M.L.; Orian, L.; Venerando, R.; Roveri, A.; Vučković, A.-M.; Travain, V.B.; Zaccarin, M.; Zennaro, L.; et al. Insight into the mechanism of ferroptosis inhibition by ferrostatin-1. Redox Biol. 2019, 28, 101328. [Google Scholar] [CrossRef]
  46. Kajarabille, N.; Latunde-Dada, G.O. Programmed Cell-Death by Ferroptosis: Antioxidants as Mitigators. Int. J. Mol. Sci. 2019, 20, 4968. [Google Scholar] [CrossRef] [Green Version]
  47. Stockwell, B.R.; Angeli, J.P.F.; Bayir, H.; Bush, A.; Conrad, M.; Dixon, S.J.; Fulda, S.; Gascón, S.; Hatzios, S.K.; Kagan, V.E.; et al. Ferroptosis: A Regulated Cell Death Nexus Linking Metabolism, Redox Biology, and Disease. Cell 2017, 171, 273–285. [Google Scholar] [CrossRef] [Green Version]
  48. Hou, Q.; Hsu, Y.-T. Bax translocates from cytosol to mitochondria in cardiac cells during apoptosis: Development of a GFP-Bax-stable H9c2 cell line for apoptosis analysis. Am. J. Physiol. Circ. Physiol. 2005, 289, H477–H487. [Google Scholar] [CrossRef] [PubMed]
  49. Jungas, T.; Motta, I.; Duffieux, F.; Fanen, P.; Stoven, V.; Ojcius, D.M. Glutathione Levels and BAX Activation during Apoptosis Due to Oxidative Stress in Cells Expressing Wild-type and Mutant Cystic Fibrosis Transmembrane Conductance Regulator. J. Biol. Chem. 2002, 277, 27912–27918. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Chen, J.; Zhu, R.L.; Nakayama, M.; Kawaguchi, K.; Jin, K.; Stetler, R.A.; Simon, R.P.; Graham, S.H. Expression of the Apoptosis-Effector Gene, Bax, Is Up-Regulated in Vulnerable Hippocampal CA1 Neurons Following Global Ischemia. J. Neurochem. 2002, 67, 64–71. [Google Scholar] [CrossRef] [PubMed]
  51. Krajewski, S.; Mai, J.K.; Krajewska, M.; Sikorska, M.; Mossakowski, M.J.; Reed, J.C. Upregulation of bax protein levels in neurons following cerebral ischemia. J. Neurosci. 1995, 15, 6364–6376. [Google Scholar] [CrossRef]
  52. Wei, C.; Kumar, S.; Kim, I.-K.; Gupta, S. Thymosin beta 4 protects cardiomyocytes from oxidative stress by targeting ant-ioxidative enzymes and anti-apoptotic genes. PLoS ONE 2012, 7, e42586. [Google Scholar] [CrossRef]
  53. Kumar, S.; Gupta, S. Thymosin beta 4 prevents oxidative stress by targeting antioxidant and anti-apoptotic genes in cardiac fibroblasts. PLoS ONE 2011, 6, e26912. [Google Scholar] [CrossRef]
  54. Guo, S.; Wharton, W.; Moseley, P.; Shi, H. Heat shock protein 70 regulates cellular redox status by modulating glutathi-one-related enzyme activities. Cell Stress Chaperones 2007, 12, 245. [Google Scholar] [CrossRef]
  55. Papadopoulos, M.C.; Sun, X.Y.; Cao, J.; Mivechi, N.F.; Giffard, R.G. Over-expression of HSP-70 protects astrocytes from combined oxygen-glucose deprivation. NeuroReport 1996, 7, 429–432. [Google Scholar] [CrossRef]
  56. Plumier, J.-C.L.; Krueger, A.M.; Currie, R.W.; Kontoyiannis, D.; Kollias, G.; Pagoulatos, G.N. Transgenic mice expressing the human inducible Hsp70 have hippocampal neurons resistant to ischemic injury. Cell Stress Chaperon 1997, 2, 162–167. [Google Scholar] [CrossRef] [Green Version]
  57. Aufricht, C.; Lu, E.; Thulin, G.; Kashgarian, M.; Siegel, N.J.; Van Why, S.K. ATP releases HSP-72 from protein aggregates after renal ischemia. Am. J. Physiol. Content 1998, 274, F268–F274. [Google Scholar] [CrossRef] [PubMed]
  58. Bidmon, B.; Endemann, M.; Müller, T.; Arbeiter, K.; Herkner, K.; Aufricht, C. Heat shock protein-70 repairs proximal tubule structure after renal ischemia. Kidney Int. 2000, 58, 2400–2407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Calabrese, V.; Copani, A.; Testa, D.; Ravagna, A.; Spadaro, F.; Tendi, E.; Nicoletti, V.; Stella, A.G. Nitric oxide synthase induction in astroglial cell cultures: Effect on heat shock protein 70 synthesis and oxidant/antioxidant balance. J. Neurosci. Res. 2000, 60, 613–622. [Google Scholar] [CrossRef]
  60. McLaughlin, B.; Hartnett, K.A.; Erhardt, J.A.; Legos, J.J.; White, R.F.; Barone, F.C.; Aizenman, E. Caspase 3 activation is essential for neuroprotection in preconditioning. Proc. Natl. Acad. Sci. USA 2003, 100, 715–720. [Google Scholar] [CrossRef] [Green Version]
  61. Yan, L.-J.; Christians, E.S.; Liu, L.; Xiao, X.; Sohal, R.S.; Benjamin, I.J. Mouse heat shock transcription factor 1 deficiency alters cardiac redox homeostasis and increases mitochondrial oxidative damage. EMBO J. 2002, 21, 5164–5172. [Google Scholar] [CrossRef] [Green Version]
  62. Chiang, S.-K.; Chen, S.-E.; Chang, L.-C. A dual role of heme oxygenase-1 in cancer cells. Int. J. Mol. Sci. 2018, 20, 39. [Google Scholar] [CrossRef] [Green Version]
  63. Jozkowicz, A.; Was, H.; Dulak, J. Heme oxygenase-1 in tumors: Is it a false friend? Antioxid. Redox Signal. 2007, 9, 2099–2118. [Google Scholar] [CrossRef] [Green Version]
  64. Hsu, Y.-Y.; Chen, C.-S.; Wu, S.-N.; Jong, Y.-J.; Lo, Y.-C. Berberine activates Nrf2 nuclear translocation and protects against oxidative damage via a phosphatidylinositol 3-kinase/Akt-dependent mechanism in NSC34 motor neuron-like cells. Eur. J. Pharm. Sci. 2012, 46, 415–425. [Google Scholar] [CrossRef] [PubMed]
  65. Kwon, M.-Y.; Park, E.; Lee, S.-J.; Chung, S.W. Heme oxygenase-1 accelerates erastin-induced ferroptotic cell death. Oncotarget 2015, 6, 24393–24403. [Google Scholar] [CrossRef] [Green Version]
  66. Sun, X.; Ou, Z.; Chen, R.; Niu, X.; Chen, D.; Kang, R.; Tang, D. Activation of the p62-Keap1-NRF2 pathway protects against ferroptosis in hepatocellular carcinoma cells. Hepatology 2015, 63, 173–184. [Google Scholar] [CrossRef] [PubMed]
  67. Gasdaska, P.Y.; Gasdaska, J.R.; Cochran, S.; Powis, G. Cloning and sequencing of a human thioredoxin reductase. FEBS Lett. 1995, 373, 5–9. [Google Scholar] [CrossRef] [Green Version]
  68. Cai, L.L.; Ruberto, R.A.; Ryan, M.J.; Eaton, J.K.; Schreiber, S.L.; Viswanathan, V.S. Modulation of ferroptosis sensitivity by TXNRD1 in pancreatic cancer cells. bioRxiv 2020. [Google Scholar] [CrossRef]
  69. Yang, L.; Wang, H.; Yang, X.; Wu, Q.; An, P.; Jin, X.; Liu, W.; Huang, X.; Li, Y.; Yan, S.; et al. Auranofin mitigates systemic iron overload and induces ferroptosis via distinct mechanisms. Signal Transduct. Target. Ther. 2020, 5, 138. [Google Scholar] [CrossRef] [PubMed]
  70. Low, T.L.; Mercer, R.C. Isolation and structural studies of porcine, ovine and murine thymosin β4 by high-performance liquid chromatography. J. Chromatogr. A 1984, 301, 221–239. [Google Scholar] [CrossRef]
  71. Snyder, W.; Cook, M.; Nasset, E.; Karhausen, L.; Howells, G.P.; Tipton, I. Report of the task group on reference man; Pergamon Oxford: Oxford, UK, 1975; Volume 23. [Google Scholar]
  72. Torti, S.V.; Torti, F.M. Iron and cancer: More ore to be mined. Nat. Cancer 2013, 13, 342–355. [Google Scholar] [CrossRef] [Green Version]
  73. Cheah, J.H.; Kim, S.F.; Hester, L.D.; Clancy, K.W.; Patterson, S.E., III; Papadopoulos, V.; Snyder, S.H. NMDA receptor-nitric oxide transmission mediates neuronal iron homeostasis via the GTPase Dexras1. Neuron 2006, 51, 431–440. [Google Scholar] [CrossRef] [Green Version]
  74. Choi, D.W. Glutamate neurotoxicity and diseases of the nervous system. Neuron 1988, 1, 623–634. [Google Scholar] [CrossRef]
  75. Murphy, T.H.; Miyamoto, M.; Sastre, A.; Schnaar, R.L.; Coyle, J.T. Glutamate toxicity in a neuronal cell line involves inhi-bition of cystine transport leading to oxidative stress. Neuron 1989, 2, 1547–1558. [Google Scholar] [CrossRef]
  76. Forciniti, S.; Greco, L.; Grizzi, F.; Malesci, A.; Laghi, L. Iron metabolism in cancer progression. Int. J. Mol. Sci. 2020, 21, 2257. [Google Scholar] [CrossRef] [Green Version]
  77. Hametner, S.; Wimmer, I.; Haider, L.; Pfeifenbring, S.; Brück, W.; Lassmann, H. Iron and neurodegeneration in the multiple sclerosis brain. Ann. Neurol. 2013, 74, 848–861. [Google Scholar] [CrossRef] [PubMed]
  78. Muhoberac, B.B.; Vidal, R. Abnormal iron homeostasis and neurodegeneration. Front. Aging Neurosci. 2013, 5, 32. [Google Scholar] [CrossRef] [Green Version]
  79. Yang, W.S.; Stockwell, B.R. Ferroptosis: Death by Lipid Peroxidation. Trends Cell Biol. 2015, 26, 165–176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Tan, S.; Schubert, D.; Maher, P. Oxytosis: A Novel Form of Programmed Cell Death. Curr. Top. Med. Chem. 2001, 1, 497–506. [Google Scholar] [CrossRef] [PubMed]
  81. Kang, Y.; Tiziani, S.; Park, G.; Kaul, M.; Paternostro, G. Cellular protection using Flt3 and PI3Kα inhibitors demonstrates multiple mechanisms of oxidative glutamate toxicity. Nat. Commun. 2014, 5, 3672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Morrison, B.; Pringle, A.; McManus, T.; Ellard, J.; Bradley, M.; Signorelli, F.; Iannotti, F.; E Sundstrom, L. L -Arginyl-3,4-Spermidine is neuroprotective in several in vitro models of neurodegeneration and in vivo ischaemia without suppressing synaptic transmission. Br. J. Pharmacol. 2002, 137, 1255–1268. [Google Scholar] [CrossRef]
  83. Sundstrom, L.; Pringle, A.; Morrison, B.; Bradley, M. Organotypic cultures as tools for functional screening in the CNS. Drug Discov. Today 2005, 10, 993–1000. [Google Scholar] [CrossRef]
  84. Zhu, Z.; Zhang, Y.; Huang, X.; Can, L.; Zhao, X.; Wang, Y.; Xue, J.; Cheng, M.; Zhu, L. Thymosin beta 4 alleviates non-alcoholic fatty liver by inhibiting ferroptosis via up-regulation of GPX4. Eur. J. Pharmacol. 2021, 908, 174351. [Google Scholar] [CrossRef]
  85. Alghrably, M.; Dudek, D.; Emwas, A.-H.; Jaremko, M.; Rowińska-Żyrek, M. Copper(II) and Amylin Analogues: A Complicated Relationship. Inorg. Chem. 2020, 59, 2527–2535. [Google Scholar] [CrossRef] [PubMed]
  86. Sharfalddin, A.A.; Emwas, A.H.; Jaremko, M.; Hussien, M.A. Transition metal complexes of 6-mercaptopurine: Characterization, Theoretical calculation, DNA-Binding, molecular docking, and anticancer activity. Appl. Organomet. Chem. 2021, 35, e6041. [Google Scholar] [CrossRef]
  87. Glasoe, P.K.; Long, F.A. USE OF GLASS ELECTRODES TO MEASURE ACIDITIES IN DEUTERIUM OXIDE1,2. J. Phys. Chem. 1960, 64, 188–190. [Google Scholar] [CrossRef]
  88. Willker, W.; Wollborn, U.; Leibfritz, D. Exact Measurement of 3JCH Coupling Constants Using Proton-Detected Editing and Selection Sequences. J. Magn. Reson. Ser. B 1993, 101, 83–86. [Google Scholar] [CrossRef]
  89. Jorgensen, W.L.; Maxwell, D.S.; Tirado-Rives, J. Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  90. MacroModel; Schrödinger LLC: New York, NY, USA, 2020.
  91. Chang, G.; Guida, W.C.; Still, W.C. An internal-coordinate Monte Carlo method for searching conformational space. J. Am. Chem. Soc. 1989, 111, 4379–4386. [Google Scholar] [CrossRef]
  92. Polak, E.; Ribiere, G. Note sur la convergence de méthodes de directions conjuguées. ESAIM Math. Model. Numer. Anal. 1969, 3, 35–43. [Google Scholar] [CrossRef]
  93. Ponder, J.W.; Richards, F.M. An efficient newton-like method for molecular mechanics energy minimization of large mole-cules. J. Comput. Chem. 1987, 8, 1016–1024. [Google Scholar] [CrossRef]
  94. Horn, H.W.; Swope, W.C.; Pitera, J.W.; Madura, J.; Dick, T.J.; Hura, G.L.; Head-Gordon, T. Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. J. Chem. Phys. 2004, 120, 9665–9678. [Google Scholar] [CrossRef]
  95. Best, R.B.; Hummer, G. Optimized Molecular Dynamics Force Fields Applied to the Helix−Coil Transition of Polypeptides. J. Phys. Chem. B 2009, 113, 9004–9015. [Google Scholar] [CrossRef] [Green Version]
  96. Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J.L.; Dror, R.O.; Shaw, D.E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins Struct. Funct. Bioinform. 2010, 78, 1950–1958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. Li, P.; Song, L.F.; Merz, K.M., Jr. Parameterization of highly charged metal ions using the 12-6-4 LJ-type nonbonded model in explicit water. J. Phys. Chem. B 2015, 119, 883–895. [Google Scholar] [CrossRef] [Green Version]
  98. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  99. Hess, B.; Bekker, H.; Berendsen, H.J.; Fraaije, J.G. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  100. Feenstra, K.A.; Hess, B.; Berendsen, H.J. Improving efficiency of large time-scale molecular dynamics simulations of hydrogen-rich systems. J. Comput. Chem. 1999, 20, 786–798. [Google Scholar] [CrossRef]
  101. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  102. Durzyńska, J.; Barton, E. IGF expression in HPV-related and HPV-unrelated human cancer cells. Oncol. Rep. 2014, 32, 893–900. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Scheme 1. Primary sequence and 3D structure of Tβ4 (PDB: 4PL7). Amino acids in parentheses were not defined by the X-ray structure and were added using Modeller v. 10.0 following the proper secondary structure elements (alpha helix). The metal-binding sites (I–IV) are signed in red. The scheme was created by Joanna Izabela Lachowicz, University of Cagliari. This image is original and designed specifically for the targeted publication. It should not be cropped, distorted, or in any way edited without the expressed consent of Research Publication Services.
Scheme 1. Primary sequence and 3D structure of Tβ4 (PDB: 4PL7). Amino acids in parentheses were not defined by the X-ray structure and were added using Modeller v. 10.0 following the proper secondary structure elements (alpha helix). The metal-binding sites (I–IV) are signed in red. The scheme was created by Joanna Izabela Lachowicz, University of Cagliari. This image is original and designed specifically for the targeted publication. It should not be cropped, distorted, or in any way edited without the expressed consent of Research Publication Services.
Ijms 23 00551 sch001
Figure 1. (a) The overlay of the two-dimensional 13C-1H HSQC spectra in free form (red) and with thymosin: Al3+ ratio 1:5 (green). The alpha and beta proton regions are highlighted in the red and blue dotted rectangles in (a) and are presented in (b,c), respectively. Noticeable chemical shift changes were observed for F12, I9 as indicated by red circles in (b). Otherwise, the chemical shifts remain the same, with only some decrease in the intensities, as highlighted by the red circles in (c).
Figure 1. (a) The overlay of the two-dimensional 13C-1H HSQC spectra in free form (red) and with thymosin: Al3+ ratio 1:5 (green). The alpha and beta proton regions are highlighted in the red and blue dotted rectangles in (a) and are presented in (b,c), respectively. Noticeable chemical shift changes were observed for F12, I9 as indicated by red circles in (b). Otherwise, the chemical shifts remain the same, with only some decrease in the intensities, as highlighted by the red circles in (c).
Ijms 23 00551 g001
Figure 2. The chemical shift perturbation (CSP) plot for the 48 distinct peaks from the overlay of the two-dimensional 13C-1H HSQC spectra in free form and with thymosin:Al3+ ratio 1:5.The formula used for CSP is Sqrt [1/2(δH2 + α ∗ δ13C2)] as suggested by Willamson [26]. The α was kept 0.3.
Figure 2. The chemical shift perturbation (CSP) plot for the 48 distinct peaks from the overlay of the two-dimensional 13C-1H HSQC spectra in free form and with thymosin:Al3+ ratio 1:5.The formula used for CSP is Sqrt [1/2(δH2 + α ∗ δ13C2)] as suggested by Willamson [26]. The α was kept 0.3.
Ijms 23 00551 g002
Figure 3. Representative structures of (A) Tβ4N-C; (B) Tβ4N-mid and (C) Tβ4N-N. Figure was created by Gabriele dalla Torre, Euskal Herriko Unibertsitatea UPV/EHU. This image is original and designed specifically for the targeted publication. It should not be cropped, distorted, or in any way edited without the expressed consent of Research Publication Services.
Figure 3. Representative structures of (A) Tβ4N-C; (B) Tβ4N-mid and (C) Tβ4N-N. Figure was created by Gabriele dalla Torre, Euskal Herriko Unibertsitatea UPV/EHU. This image is original and designed specifically for the targeted publication. It should not be cropped, distorted, or in any way edited without the expressed consent of Research Publication Services.
Ijms 23 00551 g003
Figure 4. Electron micrographs of J774cells. (A,B) The Control group displays a well-preserved ultrastructure. (C,D) Administration of Tβ4 (10 μM, 24 h) does not appear to alter cell morphology. Note in both control and Tβ4 treated samples the presence of several mitochondria showing a continuous mitochondria membrane and well-organized cristae. M = mitochondria, N = nucleus. (E,F) Glutamate treatment (5.5 mM, 24 h) induces numerous cellular alterations that affect especially the mitochondria (arrows). Initially, it appears that the membranes of the mitochondria open up, letting the contents pour out (arrowhead). Note the presence of tunneling nanotubes connecting adjacent cells (asterisk). (G,H) After Tβ4 administration (10 μM, 20 h), the cytoplasmic compartment displays more preserved cell structures. Mitochondria show a continuous mitochondrial membrane and partially reorganized cristae. N = nucleus. (I,J) Erastin treatment (0.5 μM; 24 h) induces numerous cellular alterations that affect especially the mitochondria (arrows). Mitochondria are frequently observed with ruptured mitochondrial membrane and altered cristae (arrowheads). (K,L) After Tβ4 administration (10 μM, 20 h), the cytoplasmic compartment displays more preserved cell structures associated with an evident decrease in lysosomal activity. N = nucleus.
Figure 4. Electron micrographs of J774cells. (A,B) The Control group displays a well-preserved ultrastructure. (C,D) Administration of Tβ4 (10 μM, 24 h) does not appear to alter cell morphology. Note in both control and Tβ4 treated samples the presence of several mitochondria showing a continuous mitochondria membrane and well-organized cristae. M = mitochondria, N = nucleus. (E,F) Glutamate treatment (5.5 mM, 24 h) induces numerous cellular alterations that affect especially the mitochondria (arrows). Initially, it appears that the membranes of the mitochondria open up, letting the contents pour out (arrowhead). Note the presence of tunneling nanotubes connecting adjacent cells (asterisk). (G,H) After Tβ4 administration (10 μM, 20 h), the cytoplasmic compartment displays more preserved cell structures. Mitochondria show a continuous mitochondrial membrane and partially reorganized cristae. N = nucleus. (I,J) Erastin treatment (0.5 μM; 24 h) induces numerous cellular alterations that affect especially the mitochondria (arrows). Mitochondria are frequently observed with ruptured mitochondrial membrane and altered cristae (arrowheads). (K,L) After Tβ4 administration (10 μM, 20 h), the cytoplasmic compartment displays more preserved cell structures associated with an evident decrease in lysosomal activity. N = nucleus.
Ijms 23 00551 g004
Figure 5. Dose-response relationship for inhibition of erastin-induced (10 μM, 24 h) death in J774 cells by Fer-1 (A) and Tβ4 (B). Cell viability quantitative data are shown as mean and standard error. Statistical comparison of the cell viability quantitative data in different experimental conditions were performed, with 5.0 GraphPad Prism software, by the unpaired t-test. The data were also analyzed with one-way ANOVA. p values < 0.05 was considered significant.
Figure 5. Dose-response relationship for inhibition of erastin-induced (10 μM, 24 h) death in J774 cells by Fer-1 (A) and Tβ4 (B). Cell viability quantitative data are shown as mean and standard error. Statistical comparison of the cell viability quantitative data in different experimental conditions were performed, with 5.0 GraphPad Prism software, by the unpaired t-test. The data were also analyzed with one-way ANOVA. p values < 0.05 was considered significant.
Ijms 23 00551 g005
Scheme 2. Tβ4 inhibitory activity of ferroptosis. An excess of L-Glu in the extracellular matrix leads to the inhibition of GSH synthesis and reduced anti-oxidative activity, which increases the concentration of free iron ions in the cell and enhances the production of reactive oxygen species (ROS). High ROS and free iron ion concentrations trigger ferroptosis. Ferroptosis leads to lipid membrane disruption by ROS, and silences the expression of HO-1, Hsp70, and BAX. Extracellular treatment with Tβ4 can inhibit ferroptosis by increasing the expression of HO-1, Hsp70, BAX, and TXNRD-1, and direct chelation of free iron ions. The scheme was created by Heno Hwang, scientific illustrator at King Abdullah University of Science and Technology (KAUST). This image is original and designed specifically for the targeted publication. It should not be cropped, distorted, or in any way edited without the expressed consent of Research Publication Services.
Scheme 2. Tβ4 inhibitory activity of ferroptosis. An excess of L-Glu in the extracellular matrix leads to the inhibition of GSH synthesis and reduced anti-oxidative activity, which increases the concentration of free iron ions in the cell and enhances the production of reactive oxygen species (ROS). High ROS and free iron ion concentrations trigger ferroptosis. Ferroptosis leads to lipid membrane disruption by ROS, and silences the expression of HO-1, Hsp70, and BAX. Extracellular treatment with Tβ4 can inhibit ferroptosis by increasing the expression of HO-1, Hsp70, BAX, and TXNRD-1, and direct chelation of free iron ions. The scheme was created by Heno Hwang, scientific illustrator at King Abdullah University of Science and Technology (KAUST). This image is original and designed specifically for the targeted publication. It should not be cropped, distorted, or in any way edited without the expressed consent of Research Publication Services.
Ijms 23 00551 sch002
Figure 6. mRNA expression levels (mean and standard error) of (A) BAX, (B) HSP70 (C) HO-1, and (D) TRNRD-1, which are related to stress conditions, in J774 cells. The amount of mRNA of each gene for each experiment was normalized with the respective value of the beta-actin housekeeping gene mRNA (X-fold of beta-actin expression). Statistical comparison of the gene’s expression in different experimental conditions was performed with 5.0 GraphPad Prism software. The data were also analyzed with one-way ANOVA. p values < 0.05 was considered significant.
Figure 6. mRNA expression levels (mean and standard error) of (A) BAX, (B) HSP70 (C) HO-1, and (D) TRNRD-1, which are related to stress conditions, in J774 cells. The amount of mRNA of each gene for each experiment was normalized with the respective value of the beta-actin housekeeping gene mRNA (X-fold of beta-actin expression). Statistical comparison of the gene’s expression in different experimental conditions was performed with 5.0 GraphPad Prism software. The data were also analyzed with one-way ANOVA. p values < 0.05 was considered significant.
Ijms 23 00551 g006
Table 1. Amino acid residues affected by metal/peptide interactions. Metal-binding regions are also signed in Scheme 1.
Table 1. Amino acid residues affected by metal/peptide interactions. Metal-binding regions are also signed in Scheme 1.
Binding RegionIIIIIIIV
Fe (II)S1K3P4-M6--K11F12D13S15L17T22P27L28P29S30T33I34--A40S43
Fe (III)S1K3P4-M6I9--F12D13S15-T22-L28--T33I34Q35E36A40-
Al (III)--P4D5M6I9E10-F12D13----L28---I34Q35E36--
Table 2. Percentage of hydrogen bond and salt bridge occupancies during the simulations. bb: backbone, sc: side chain.
Table 2. Percentage of hydrogen bond and salt bridge occupancies during the simulations. bb: backbone, sc: side chain.
TB4N-CTB4N-midTB4N-N
DonorAcceptorOcc. (%)DonorAcceptorOcc. (%)DonorAcceptorOcc. (%)
GLU35-NbbTHR33-Osc27.6LYS38-NbbGLU35-Obb26.0THR22-NbbLYS19-Obb23.8
LEU17-NbbLYS14-Obb25.8LEU28-NbbASN26-Osc28.2GLN23-NbbLYS18-Obb22.2
LEU28-NbbASN26-Osc36.6ILE34-NbbLYS31-Obb26.7GLU21-NbbLEU17-Obb23.4
THR33-NbbGLN36-Obb32.0LYS31-NbbLEU28-Obb21.1ILE9-NbbMET6-Obb53.0
LYS31-NbbLEU28-Obb22.5---GLU8-NbbASP5-Obb51.5
------LYS38-NbbGLU35-Obb20.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lachowicz, J.I.; Pichiri, G.; Piludu, M.; Fais, S.; Orrù, G.; Congiu, T.; Piras, M.; Faa, G.; Fanni, D.; Dalla Torre, G.; et al. Thymosin β4 Is an Endogenous Iron Chelator and Molecular Switcher of Ferroptosis. Int. J. Mol. Sci. 2022, 23, 551. https://doi.org/10.3390/ijms23010551

AMA Style

Lachowicz JI, Pichiri G, Piludu M, Fais S, Orrù G, Congiu T, Piras M, Faa G, Fanni D, Dalla Torre G, et al. Thymosin β4 Is an Endogenous Iron Chelator and Molecular Switcher of Ferroptosis. International Journal of Molecular Sciences. 2022; 23(1):551. https://doi.org/10.3390/ijms23010551

Chicago/Turabian Style

Lachowicz, Joanna I., Giusi Pichiri, Marco Piludu, Sara Fais, Germano Orrù, Terenzio Congiu, Monica Piras, Gavino Faa, Daniela Fanni, Gabriele Dalla Torre, and et al. 2022. "Thymosin β4 Is an Endogenous Iron Chelator and Molecular Switcher of Ferroptosis" International Journal of Molecular Sciences 23, no. 1: 551. https://doi.org/10.3390/ijms23010551

APA Style

Lachowicz, J. I., Pichiri, G., Piludu, M., Fais, S., Orrù, G., Congiu, T., Piras, M., Faa, G., Fanni, D., Dalla Torre, G., Lopez, X., Chandra, K., Szczepski, K., Jaremko, L., Ghosh, M., Emwas, A. -H., Castagnola, M., Jaremko, M., Hannappel, E., & Coni, P. (2022). Thymosin β4 Is an Endogenous Iron Chelator and Molecular Switcher of Ferroptosis. International Journal of Molecular Sciences, 23(1), 551. https://doi.org/10.3390/ijms23010551

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop